Exploring Dihydroflavonol-4-Reductase Reactivity and Selectivity by QM/MM-MD Simulations.

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.

Novel scaffold of natural compound eliciting sweet taste revealed by machine learning.

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.

Allosteric Modulation Mechanism of the mGluR5 Transmembrane Domain.

Positive allosteric modulators (PAMs) of metabotropic glutamate receptor type 5 (mGluR5), 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 mGluR5 PAMs. We studied the actions of positive and negative allosteric modulators within the transmembrane domain of mGluR5, 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 mGluR5 as well as for other class C GPCRs.

Conserved Residues Control the T1R3-Specific Allosteric Signaling Pathway of the Mammalian Sweet-Taste Receptor.

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.

Update of the ATTRACT force field for the prediction of protein-protein binding affinity.

Determining the protein-protein interactions is still a major challenge for molecular biology. Docking protocols has come of age in predicting the structure of macromolecular complexes. However, they still lack accuracy to estimate the binding affinities, the thermodynamic quantity that drives the formation of a complex. Here, an updated version of the protein-protein ATTRACT force field aiming at predicting experimental binding affinities is reported. It has been designed on a dataset of 218 protein-protein complexes. The correlation between the experimental and predicted affinities reaches 0.6, outperforming most of the available protocols. Focusing on a subset of rigid and flexible complexes, the performance raises to 0.76 and 0.69, respectively. © 2017 Wiley Periodicals, Inc.

Bornyl-diphosphate synthase from Lavandula angustifolia: A major monoterpene synthase involved in essential oil quality.

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.

Sweetness prediction of natural compounds.

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.

The anatomy of mammalian sweet taste receptors.

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.

Fine-tuning of microsolvation and hydrogen bond interaction regulates substrate channelling in the course of flavonoid biosynthesis.

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.

G protein-coupled odorant receptors: From sequence to structure.

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.

Discrimination between olfactory receptor agonists and non-agonists.

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.

O2migration rates in [NiFe] hydrogenases. A joint approach combining free-energy calculations and kinetic modeling.

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.

Kinetics and thermodynamics of gas diffusion in a NiFe hydrogenase.

We have investigated O2 and H2 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 O2 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.

Theoretical investigations of the role played by quercetinase enzymes upon the flavonoids oxygenolysis mechanism.

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.

Deciphering the selectivity of Bombyx mori pheromone binding protein for bombykol over bombykal: a theoretical approach.

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.

Molecular simulations enlighten the binding mode of quercetin to lipoxygenase-3.

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.

Cerium(IV)-mediated oxidation of flavonol with relevance to flavonol 2,4-dioxygenase. Direct evidence for spin delocalization in the flavonoxy radical.

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.

Molecular simulations bring new insights into flavonoid/quercetinase interaction modes.

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.