Search

Analysis of the role of electrostatic interactions in the amide hydrolysis reaction shows the assistance of an ancillary water molecule in formamide hydrolysis. A neutral reaction and an acid-catalyzed reaction are studied to obtain a mechanism. Assisted reaction involves a hydrogen bonded reaction intermediate between water and formamide and is due to a cooperative effect in the water dimer. H-acceptor in the water dimer favors the pyramidalization of the nitrogen atom to activate carbon nucleophilic attack.

Flavonoids are secondary metabolites ubiquitously found in plants. Their antioxidant properties make them highly interesting natural compounds for use in pharmacology. Therefore, unravelling the mechanisms of flavonoid biosynthesis is an important challenge. Among all the enzymes involved in this biosynthetic pathway, dihydroflavonol-4-reductase (DFR) plays a key role in the production of anthocyanins and proanthocyanidins. Here, we provide new information on the mechanism of action of this enzyme by using QM/MM-MD simulations applied to both dihydroquercetin (DHQ) and dihydrokaempferol (DHK) substrates. The consideration of these very similar compounds shed light on the major role played by the enzyme on the stabilization of the transition state but also on the activation of the substrate before the reaction through near-attack conformer effects., (© 2021 Wiley-VCH GmbH.)

Sugar replacement is still an active issue in the food industry. The use of structure-taste relationships remains one of the most rational strategy to expand the chemical space associated to sweet taste. A new machine learning model has been setup based on an update of the SweetenersDB and on open-source molecular features. It has been implemented on a freely accessible webserver. Cellular functional assays show that the sweet taste receptor is activated in vitro by a new scaffold of natural compounds identified by the in silico protocol. The newly identified sweetener belongs to the lignan chemical family and opens a new chemical space to explore., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Positive allosteric modulators (PAMs) of metabotropic glutamate receptor type 5 (mGluR 5 ), a prototypical class C G protein-coupled receptor (GPCR), have shown therapeutic potential for various neurological disorders. Understanding the allosteric activation mechanism is essential for the rational design of mGluR 5 PAMs. We studied the actions of positive and negative allosteric modulators within the transmembrane domain of mGluR 5 , using enhance-sampling all-atom molecular dynamics simulations. We found dual binding modes of the PAM, associated with distinct shapes of the allosteric pocket. The negative allosteric modulators, in contrast, showed only one binding mode. The simulations revealed the mechanism by which the PAM activated the receptor, in the absence of the orthosteric agonist (the so-called allosteric agonism). The mechanism relied on dynamic communications between amino-acid motifs that are highly conserved across class C GPCRs. The findings may guide structure-based design and virtual screening of allosteric modulators for mGluR 5 as well as for other class C GPCRs.

Mammalian sensory systems detect sweet taste through the activation of a single heteromeric T1R2/T1R3 receptor belonging to class C G-protein-coupled receptors. Allosteric ligands are known to interact within the transmembrane domain, yet a complete view of receptor activation remains elusive. By combining site-directed mutagenesis with computational modeling, we investigate the structure and dynamics of the allosteric binding pocket of the T1R3 sweet-taste receptor in its apo form, and in the presence of an allosteric ligand, cyclamate. A novel positively charged residue at the extracellular loop 2 is shown to interact with the ligand. Molecular dynamics simulations capture significant differences in the behavior of a network of conserved residues with and without cyclamate, although they do not directly interact with the allosteric ligand. Structural models show that they adopt alternate conformations, associated with a conformational change in the transmembrane region. Site-directed mutagenesis confirms that these residues are unequivocally involved in the receptor function and the allosteric signaling mechanism of the sweet-taste receptor. Similar to a large portion of the transmembrane domain, they are highly conserved among mammals, suggesting an activation mechanism that is evolutionarily conserved. This work provides a structural basis for describing the dynamics of the receptor, and for the rational design of new sweet-taste modulators., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Lavender essential oils (EOs) of higher quality are produced by a few Lavandula angustifolia cultivars and mainly used in the perfume industry. Undesirable compounds such as camphor and borneol are also synthesized by lavender leading to a depreciated EO. Here, we report the cloning of bornyl diphosphate synthase of lavender (LaBPPS), an enzyme that catalyzes the production of bornyl diphosphate (BPP) and then by-products such as borneol or camphor, from an EST library. Compared to the BPPS of Salvia officinalis, the functional characterization of LaBPPS showed several differences in amino acid sequence, and the distribution of catalyzed products. Molecular modeling of the enzyme's active site suggests that the carbocation intermediates are more stable in LaBPPS than in SoBPPS leading probably to a lower efficiency of LaBPPS to convert GPP into BPP. Quantitative RT-PCR performed from leaves and flowers at different development stages of L. angustifolia samples show a clear correlation between transcript level of LaBPPS and accumulation of borneol/camphor, suggesting that LaBPPS is mainly responsible of in vivo biosynthesis of borneol/camphor in fine lavender. A phylogenetic analysis of terpene synthases (TPS) pointed out the basal position of LaBPPS in the TPSb clade, suggesting that LaBPPS could be an ancestor of others lavender TPSb. Finally, borneol could be one of the first monoterpenes to be synthesized in the Lavandula subgenus. Knowledge gained from these experiments will facilitate future studies to improve the lavender oils through metabolic engineering or plant breeding. Accession numbers: LaBPPS: KM015221., (Copyright © 2017. Published by Elsevier Ltd.)

Based on the most exhaustive database of sweeteners with known sweetness values, a new quantitative structure-activity relationship model for sweetness prediction has been set up. Analysis of the physico-chemical properties of sweeteners in the database indicates that the structure of most potent sweeteners combines a hydrophobic scaffold functionalized by a limited number of hydrogen bond sites (less than 4 hydrogen bond donors and 10 acceptors), with a moderate molecular weight ranging from 350 to 450g·mol -1 . Prediction of sweetness, bitterness and toxicity properties of the largest database of natural compounds have been performed. In silico screening reveals that the majority of the predicted natural intense sweeteners comprise saponin or stevioside scaffolds. The model highlights that their sweetness potency is comparable to known natural sweeteners. The identified compounds provide a rational basis to initiate the design and chemosensory analysis of new low-calorie sweeteners., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

All sweet-tasting compounds are detected by a single G-protein coupled receptor (GPCR), the heterodimer T1R2-T1R3, for which no experimental structure is available. The sweet taste receptor is a class C GPCR, and the recently published crystallographic structures of metabotropic glutamate receptor (mGluR) 1 and 5 provide a significant step forward for understanding structure-function relationships within this family. In this article, we recapitulate more than 600 single point site-directed mutations and available structural data to obtain a critical alignment of the sweet taste receptor sequences with respect to other class C GPCRs. Using this alignment, a homology 3D-model of the human sweet taste receptor is built and analyzed to dissect out the role of key residues involved in ligand binding and those responsible for receptor activation. Proteins 2017; 85:332-341. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

In the course of metabolite formation, some multienzymatic edifices, the so-called metabolon, are formed and lead to a more efficient production of these natural compounds. One of the major features of these enzyme complexes is the facilitation of direct transfer of the metabolite between enzyme active sites by substrate channelling. Biophysical insights into substrate channelling remain scarce because the transient nature of these macromolecular complexes prevents the observation of high resolution structures. Here, using molecular modelling, we describe the substrate channelling of a flavonoid compound between DFR (dihydroflavonol-4-reductase) and LAR (leucoanthocyanidin reductase). The simulation presents crucial details concerning the kinetic, thermodynamic, and structural aspects of this diffusion. The formation of the DFR-LAR complex leads to the opening of the DFR active site giving rise to a facilitated diffusion, in about 1 μs, of the DFR product towards LAR cavity. The theoretically observed substrate channelling is supported experimentally by the fact that this metabolite, i.e. the product of the DFR enzyme, is not stable in the media. Moreover, along this path, the influence of the solvent is crucial. The metabolite remains close to the surface of the complex avoiding full solvation. In addition, when the dynamic behaviour of the system leads to a loss of interaction between the metabolite and the enzymes, water molecules through bridging H-bonds prevent the former from escaping to the bulk.

Odorant receptors (ORs) are the largest subfamily within class A G protein-coupled receptors (GPCRs). No experimental structural data of any OR is available to date and atomic-level insights are likely to be obtained by means of molecular modeling. In this article, we critically align sequences of ORs with those GPCRs for which a structure is available. Here, an alignment consistent with available site-directed mutagenesis data on various ORs is proposed. Using this alignment, the choice of the template is deemed rather minor for identifying residues that constitute the wall of the binding cavity or those involved in G protein recognition., (© 2015 The Protein Society.)

A joint approach combining free-energy calculations and calcium-imaging assays on the broadly tuned human 1G1 olfactory receptor is reported. The free energy of binding of ten odorants was computed by means of molecular-dynamics simulations. This state function allows separating the experimentally determined eight agonists from the two non-agonists. This study constitutes a proof-of-principle for the computational deorphanization of olfactory receptors., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Hydrogenases are promising candidates for the catalytic production of green energy by means of biological ways. The major impediment to such a production is rooted in their inhibition under aerobic conditions. In this work, we model dioxygen migration rates in mutants of a hydrogenase of Desulfovibrio fructusovorans. The approach relies on the calculation of the whole potential of mean force for O2 migration within the wild-type as well as in V74M, V74F, and V74Q mutant channels. The three free-energy barriers along the entire migration pathway are converted into chemical rates through modeling based on Transition State Theory. The use of such a model recovers the trend of O2 migration rates among the series.

We have investigated O₂ and H₂ transport across a NiFe hydrogenase at the atomic scale by means of computational methods. The Wild Type protein has been compared with the V74Q mutant. Two distinct methodologies have been applied to study the gas access to the active site. Temperature locally enhanced sampling simulations have emphasized the importance of protein dynamics on gas diffusion. The O₂ diffusion free energy profiles, obtained by umbrella sampling, are in agreement with the known kinetic data and show that in the V74Q mutant, the inhibition process is lowered from both a kinetic and a thermodynamic point of view., (Copyright © 2011 Wiley Periodicals, Inc.)

Quercetinase enzymatic activity consists in the addition of dioxygen onto flavonoids, some natural polyphenol compounds, leading to the production of both molecular carbon monoxide and to the structurally related depside compound. Experimental studies have reported degradation rates of various flavonoids by such enzymes that can not be directly correlated neither to the number nor to the place of the hydroxyl groups. In order to decipher the role of these functions, we have theoretically characterised the stationary points of various flavonoids oxygenolysis mechanisms by density functional quantum methods. Thus in the present study are reported the main energetic, structural and electronic features that drive this degradation. Together with previous analysis from MD simulations taking into account the dynamic behaviour of the substrate embedded in the enzyme cavity, the present results show that the role of the enzyme, in terms of structural and electronic effects, can not be neglected. Thus, we propose here that deformations of the substrate induced by the enzyme could originate the differences in the degradation rates experimentally observed.

In this article we report calculations dedicated to estimate the selectivity of the Bombyx mori pheromone binding protein towards the two closely related pheromonal components Bombykol and Bombykal. The selectivity is quantified by the binding free-energy difference, obtained either by the thermodynamic integration or by the MM-GBSA approach. In the latter, the selectivity is decomposed on a per-residue basis, which identifies the residues considered crucial for the selectivity of the protein for Bombykol over Bombykal. A discussion on the role of Bombyx mori pheromone binding protein is provided on the basis of these results.

Inhibition of lipoxygenases (LOXs) by flavonoid compounds is now well documented, but the description of the associated mechanism remains controversial due to a lack of information at the molecular level. For instance, X-ray determination of quercetin/LOX-3 system has led to a structure where the enzyme was cocrystallized with a degradation product of the substrate, which rendered the interpretation of the reported interactions between this flavonoid compound and the enzyme difficult. Molecular modeling simulations can in principle allow obtaining precious insights that could fill this lack of structural information. Thus, in this study, we have investigated various binding modes of quercetin to LOX-3 enzyme in order to understand the first step of the inhibition process, that is the association of the two entities. Molecular dynamics simulations and free energy calculations suggest that quercetin binds the metal center via its 3-hydroxychromone function. Moreover, enzyme/substrate interactions within the cavity impose steric hindrances to quercetin that may activate a direct dioxygen addition on the substrate., ((c) 2008 Wiley-Liss, Inc.)

We present a comparison of various computational approaches aiming at predicting the binding free energy in ligand-protein systems where the ligand is located within a highly hydrophobic cavity. The relative binding free energy between similar ligands is obtained by means of the thermodynamic integration (TI) method and compared to experimental data obtained through isothermal titration calorimetry measurements. The absolute free energy of binding prediction was obtained on a similar system (a pyrazine derivative bound to a lipocalin) by TI, potential of mean force (PMF) and also by means of the MMPBSA protocols. Although the TI protocol performs poorly either with an explicit or an implicit solvation scheme, the PMF calculation using an implicit solvation scheme leads to encouraging results, with a prediction of the binding affinity being 2 kcal mol(-1) lower than the experimental value. The use of an implicit solvation scheme appears to be well suited for the study of such hydrophobic systems, due to the lack of water molecules within the binding site.

Molecular dynamics simulations, using the AMBER force field, were performed to study Quercetin 2,3-Dioxygenase enzyme (Quercetinase or 2,3QD). We have analyzed the structural modifications of the active site and of the linker region between the native enzyme and the enzyme-substrate complex. New structural informations, such as an allosteric effect in the presence of the substrate, as well as description of the enzyme-substrate interactions and values of binding free energies were brought out. All these results confirm the idea that the linker encloses the substrate in the active site and also enlighten the recognition role of the substrate B-ring by the enzyme. Moreover, a specific interaction scheme has been proposed to explain the relative degradation rate of various flavonoid compounds under the oxygenolysis reaction catalyzed by the Quercetin 2,3-Dioxygenase enzyme., (2007 Wiley-Liss, Inc.)

The cerium(IV)-mediated oxidation of 3-hydroxy-4'-methylflavone (1) proceeds by H-atom abstraction forming the flavonoxy radical (7), and the subsequent combination of its resonance forms leads to the 3-hydroxy-4'-methylflavone dehydro dimer (9). The above system serves as direct evidence for the intermediacy of the flavonoxy radical, its spin delocalization, and also indirect evidence for valence tautomerism as a key step on the substrate activation both in the quercetinase and its biomimic model system.

Odorant binding proteins (OBP's) are small hydrophilic proteins, belonging to the lipocalin family dedicated to bind and transport small hydrophobic ligands. Despite many works, the mechanism of ligand binding, together with the functional role of these proteins remains a topic of debate and little is known at the atomic level. The present work reports a computational study of odorants capture and release by an OBP, using both constrained and unconstrained simulations, giving a glimpse on the molecular mechanism of chemoreception. The residues at the origin of the regulation of the protein door opening are identified and a tyrosine amino-acid together with other nearby residues appear to play a crucial role in allowing this event to occur. The simulations reveal that this tyrosine and the protein's L5 loop are implicated in the ligand contact with the protein and act as an anchoring point for the ligand. The protein structural features required for the ligand entry are highly conserved among many transport proteins, suggesting that this mechanism could somewhat be extended to some members of the larger family of lipocalin., ((c) 2007 Wiley-Liss, Inc.)

Quercetin, one of the most representative flavonoid compounds, is involved in antiradical, antioxidant, and prooxidant biological processes. Despite a constant increase of knowledge on both positive and negative activities of quercetin, it is unclear which activated form (quinone, semiquinone, or deprotonated) actually plays a role in each of these processes. Structural, electronic, and energetic characteristics of quercetin, as well as the influence of a copper ion on all of these parameters, are studied by means of quantum chemical electronic structure calculations. Introduction of thermodynamic cycles together with the role of coreactive compounds, such as reactive oxygen species, gives a glimpse of the most probable reaction schemes. Such a theoretical approach provides another hint to clarify which reaction is likely to occur within the broad range of quercetin biological activities.

Molecular dynamics simulations performed on quercetin 2,3-dioxygenase have shown the existence of a channel linking the bulk solvent and the cavity of the enzyme. Although much is known about the the oxygenolysis reaction catalyzed by this enzyme, the way dioxygen enters the active site has not been firmly established. The size, orientation and hydrophobic character of this channel suggests that it could provide an entrance for molecular dioxygen into the cavity. Free energy calculations show that such a process is likely to occur., ((c) 2006 Wiley-Liss, Inc.)

Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.

Flavonoids are naturally occurring phenol derivatives present in substantial amounts in a large variety of plants, fruits and vegetables daily eaten by humans. Most of these compounds exhibit several interesting biological activities, such as antiradical and antioxidant actions. Indeed, by complexation with specific enzymes, flavonoids are notably liable to metabolize molecular dioxygen. On the basis of experimental results describing oxygenolysis of the flavonoid quercetin, activated by the enzyme quercetin 2,3-dioxygenase (2,3-QD),ur attention has focused on the role of metal center in the activation of the substrate quercetin. Thus, in the present study, by means of DFT calculations at the B3LYP/ 6-31(+)G* level on model molecular systems, we describe different mechanisms for dioxygen metabolization by quercetin. Stationary points are described, and energetic and structural analyses along the reaction paths are reported. Our calculations show that the copper cation must act as an oxidant towards the substrate and that the reaction proceeds through a 1,3-cycloaddition.

The dynamic behavior of the HCV IRES IIId domain is analyzed by means of a 2.6-ns molecular dynamics simulation, starting from an NMR structure. The simulation is carried out in explicit water with Na+ counterions, and particle-mesh Ewald summation is used for the electrostatic interactions. In this work, we analyze selected patterns of the helix that are crucial for IRES activity and that could be considered as targets for the intervention of inhibitors, such as the hexanucleotide terminal loop (more particularly its three consecutive guanines) and the loop-E motif. The simulation has allowed us to analyze the dynamics of the loop substructure and has revealed a behavior among the guanine bases that might explain the different role of the third guanine of the GGG triplet upon molecular recognition. The accessibility of the loop-E motif and the loop major and minor groove is also examined, as well as the effect of Na+ or Mg2+ counterion within the simulation. The electrostatic analysis reveals several ion pockets, not discussed in the experimental structure. The positions of these ions are useful for locating specific electrostatic recognition sites for potential inhibitor binding.

The intramolecular transfer of the allyl group of functionalized allyl aryl ethers to an aldehyde group in the presence of Ni0 complexes was studied from chemical, electrochemical and theoretical points of view. The chemical reaction involves the addition of Ni0 to the allyl ether followed by stoichiometric allylation. The electrochemical process is catalytic in nickel and involves the reduction of intermediate eta3-allylnickel(II) complexes.

A spiro attack on a peroxo group is calculated to be the preferred reaction pathway for olefin epoxidation with the catalytic system CH 3 ReO 3 /H 2 O 2 (see picture). This finding is supported by density functional calculations on more than ten transition states for the most probable mechanisms. Hydration has significant effects on various reaction species: it stabilizes the intermediates and destabilizes, with one exception, the transition states., (© 1998 WILEY-VCH Verlag GmbH, Weinheim, Fed. Rep. of Germany.)

Neprilysin (neutral endopeptidase-24.11, EC 3.4.24.11) is a mammalian zinc-endopeptidase involved in the degradation of biologically active peptides. Although no atomic structure is available for this enzyme, site-directed mutagenesis studies have shown that its active site resembles closely that of the bacterial zinc-endopeptidase, thermolysin (EC 3.4.24.27). One active site residue of thermolysin, Arg-203, is involved in inhibitor binding by forming hydrogen bonds with the carbonyl group of a residue in the P1 position and also participates in a hydrogen bond network involving Asp-170. Sequence alignment data shows that Arg-717 of neprilysin could play a similar role to Arg-203 of thermolysin. This was investigated by site-directed mutagenesis with Arg-203 of thermolysin and Arg-717 of neprilysin being replaced by methionine residues. This led, in both cases, to decreases in kcat/Km values, of 122-fold for neprilysin and 2300-fold for thermolysin, essentially due to changes in kcat. The Ki values of several inhibitors were also increased for the mutated enzymes. In addition, the replacement of Asp-170 of thermolysin by Ala residue resulted in a decrease in kcat/Km of 220-fold. The results, coupled with a molecular modeling study, suggest that Arg-717 of neprilysin corresponds to Arg-203 of thermolysin and that in both enzymes a hydrogen bond network exists, involving His-142, Asp-170, and Arg-203 in thermolysin and His-583, Asp-650, and Arg-717 in neprilysin, which is crucial for hydrolytic activity.