Understanding the molecular mechanism of umami recognition by T1R1-T1R3 using molecular dynamics simulations
Ontology highlight
ABSTRACT:
SUBMITTER: Hai Liu
PROVIDER: S-BSST248 | biostudies-other |
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (Egap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design.
| S-EPMC6138447 | biostudies-literature
Project description:Our skin constitutes an effective permeability barrier that protects the body from exogenous substances but concomitantly severely limits the number of pharmaceutical drugs that can be delivered transdermally. In topical formulation design, chemical permeation enhancers (PEs) are used to increase drug skin permeability. In vitro skin permeability experiments can measure net effects of PEs on transdermal drug transport, but they cannot explain the molecular mechanisms of interactions between drugs, permeation enhancers, and skin structure, which limits the possibility to rationally design better new drug formulations. Here we investigate the effect of the PEs water, lauric acid, geraniol, stearic acid, thymol, ethanol, oleic acid, and eucalyptol on the transdermal transport of metronidazole, caffeine, and naproxen. We use atomistic molecular dynamics (MD) simulations in combination with developed molecular models to calculate the free energy difference between 11 PE-containing formulations and the skin's barrier structure. We then utilize the results to calculate the final concentration of PEs in skin. We obtain an RMSE of 0.58 log units for calculated partition coefficients from water into the barrier structure. We then use the modified PE-containing barrier structure to calculate the PEs' permeability enhancement ratios (ERs) on transdermal metronidazole, caffeine, and naproxen transport and compare with the results obtained from in vitro experiments. We show that MD simulations are able to reproduce rankings based on ERs. However, strict quantitative correlation with experimental data needs further refinement, which is complicated by significant deviations between different measurements. Finally, we propose a model for how to use calculations of the potential of mean force of drugs across the skin's barrier structure in a topical formulation design.
| S-EPMC10428223 | biostudies-literature
Project description:BackgroundLysine Specific Demethylase (LSD1 or KDM1A) in complex with its co-repressor protein CoREST catalyzes the demethylation of the H3 histone N-terminal tail and is currently one of the most promising epigenetic targets for drug discovery against cancer and neurodegenerative diseases. Models of non-covalent binding, such as lock and key, induced-fit, and conformational selection could help explaining the molecular mechanism of LSD1/CoREST-H3-histone association, thus guiding drug discovery and design efforts. Here, we quantify the extent to which LSD1/CoREST substrate binding is consistent with these hypothetical models using LSD1/CoREST conformational ensembles obtained through extensive explicit solvent molecular dynamics (MD) simulations.ResultsWe find that an induced-fit model is the most representative of LSD1/CoREST-H3-histone non-covalent binding and accounts for the local conformational changes occurring in the H3-histone binding site. We also show that conformational selection - despite in principle not ruled out by this finding - is minimal, and only relevant when global properties are considered, e.g. the nanoscale motion of the LSD1/CoREST clamp.ConclusionThe induced-fit mechanism revealed by our MD simulation study will aid the inclusion of protein dynamics for the discovery and design of LSD1 inhibitors targeting the H3-histone binding region. On a general basis, our study indicates the importance of using multiple metrics or selection schemes when testing alternative hypothetical mechanistic models of non-covalent binding.
| S-EPMC4175114 | biostudies-literature
Project description:Although proteins have represented the molecular target of choice in the development of new drug candidates, the pharmaceutical importance of ribonucleic acids has gradually been growing. The increasing availability of structural information has brought to light the existence of peculiar three-dimensional RNA arrangements, which can, contrary to initial expectations, be recognized and selectively modulated through small chemical entities or peptides. The application of classical computational methodologies, such as molecular docking, for the rational development of RNA-binding candidates is, however, complicated by the peculiarities characterizing these macromolecules, such as the marked conformational flexibility, the singular charges distribution, and the relevant role of solvent molecules. In this work, we have thus validated and extended the applicability domain of SuMD, an all-atoms molecular dynamics protocol that allows to accelerate the sampling of molecular recognition events on a nanosecond timescale, to ribonucleotide targets of pharmaceutical interest. In particular, we have proven the methodological ability by reproducing the binding mode of viral or prokaryotic ribonucleic complexes, as well as that of artificially engineered aptamers, with an impressive degree of accuracy.
| S-EPMC7057144 | biostudies-literature
Project description:Glycolipids like phosphatidylinositol hexamannosides (PIM6) and lipoglycans, such as lipomannan (LM) and lipoarabinomannan (LAM), play crucial roles in virulence, survival, and antibiotic resistance of various mycobacterial species. Phosphatidyl-myo-inositol mannosyltransferase A (PimA) catalyzes the transfer of the mannose moiety (M) from GDP-mannose (GDPM) to phosphatidyl-myo-inositol (PI) to synthesize GDP and phosphatidyl-myo-inositol monomannoside (PIM). This PIM is mannosylated, acylated, and further modified to give rise to the higher PIMs, LM, and LAM. It is yet to be known how PI, PIM, PI-GDPM, and PIM-GDP interact with PimA. Here, we report the docked structures of PI and PIM to understand how the substrates and the products interact with PimA. Using molecular dynamics (MD) simulations for 300 ns, we have investigated how various ligand-bound conformations change the dynamics of PimA. Our studies demonstrated the "open to closed" motions of PimA. We observed that PimA is least dynamic when bound to both GDPM and PI. MD simulations indicated that the loop residues 59-70 and the α-helical residues 73-86 of PimA play important roles while interacting with both PI and PIM. MD analyses also suggested that the residues Y9, P59, R68, L69, N97, R196, R201, K202, and R228 of PimA play significant roles in the mannose transfer reaction. Overall, docking studies and MD simulations provide crucial insights to design future therapeutic drugs against mycobacterial PimA.
| S-EPMC9202250 | biostudies-literature
Project description:We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.
| S-EPMC1301720 | biostudies-other
Project description:Because of high stability and slow unfolding rates of G-quadruplexes (G4), cells have evolved specialized helicases that disrupt these non-canonical DNA and RNA structures in an ATP-dependent manner. One example is DHX36, a DEAH-box helicase, which participates in gene expression and replication by recognizing and unwinding parallel G4s. Here, we studied the molecular basis for the high affinity and specificity of DHX36 for parallel-type G4s using all-atom molecular dynamics simulations. By computing binding free energies, we found that the two main G4-interacting subdomains of DHX36, DSM and OB, separately exhibit high G4 affinity but they act cooperatively to recognize two distinctive features of parallel G4s: the exposed planar face of a guanine tetrad and the unique backbone conformation of a continuous guanine tract, respectively. Our results also show that DSM-mediated interactions are the main contributor to the binding free energy and rely on making extensive van der Waals contacts between the GXXXG motifs and hydrophobic residues of DSM and a flat guanine plane. Accordingly, the sterically more accessible 5'-G-tetrad allows for more favorable van der Waals and hydrophobic interactions which leads to the preferential binding of DSM to the 5'-side. In contrast to DSM, OB binds to G4 mostly through polar interactions by flexibly adapting to the 5'-terminal guanine tract to form a number of strong hydrogen bonds with the backbone phosphate groups. We also identified a third DHX36/G4 interaction site formed by the flexible loop missing in the crystal structure.
| S-EPMC8114077 | biostudies-literature
Project description:p53 is a representative transcription factor that activates multiple target genes. To realize stimulus-dependent specificities, p53 has to recognize targets with structural variety, of which molecular mechanisms are largely unknown. Here, we conducted a series of long-time scale (totally more than 100-ms) coarse-grained molecular dynamics simulations, uncovering structure and dynamics of full-length p53 tetramer that recognizes its response element (RE). We obtained structures of a full-length p53 tetramer that binds to the RE, which is strikingly different from the structure of p53 at search. These structures are not only consistent with previous low-resolution or partial structural information, but also give access to previously unreachable detail, such as the preferential distribution of intrinsically disordered regions, the contacts between core domains, the DNA bending, and the connectivity of linker regions. We also explored how the RE variation affects the structure of the p53-RE complex. Further analysis of simulation trajectories revealed how p53 finds out the RE and how post-translational modifications affect the search mechanism.
| S-EPMC4656996 | biostudies-literature
Project description:Limit dextrinase (LD) is the only endogenous starch-debranching enzyme in barley (Hordeum vulgare, Hv), which is the key factor affecting the production of a high degree of fermentation. Free LD will lose its activity in the mashing process at high temperature in beer production. However, there remains a lack of understanding on the factor affecting the themostability of HvLD at the atomic level. In this work, the molecular dynamics simulations were carried out for HvLD to explore the key factors affecting the thermal stability of LD. The higher value of root mean square deviation (RMSD), radius of gyration (R g), and surface accessibility (SASA) suggests the instability of HvLD at high temperatures. Intra-protein hydrogen bonds and hydrogen bonds between protein and water decrease at high temperature. Long-lived hydrogen bonds, salt bridges, and hydrophobic contacts are lost at high temperature. The salt bridge interaction analysis suggests that these salt bridges are important for the thermostability of HvLD, including E568-R875, D317-R378, D803-R884, D457-R214, D468-R395, D456-R452, D399-R471, and D541-R542. Root mean square fluctuation (RMSF) analysis identified the thermal-sensitive regions of HvLD, which will facilitate enzyme engineering of HvLD for enhanced themostability.
| S-EPMC7241666 | biostudies-literature
Project description: Not available
| S-EPMC7571292 | biostudies-literature