scholarly journals Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations and NMR Experiments

2020 ◽  
Author(s):  
Zhengyue Zhang ◽  
Jennifer Vögele ◽  
Klaudia Mráziková ◽  
Holger Kruse ◽  
Xiaohui Cang ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Rna Structure ◽  
Rna Stability ◽  
Field Performance ◽  
Md Simulations ◽  
Dispersion Component ◽  
Rna Backbone ◽  
Dynamics Simulations ◽  
The Individual

AbstractPhosphorothioates (PTs) are important chemical modifications of the RNA backbone where a single non-bridging oxygen of the phosphate is replaced with a sulphur atom. PT can stabilize RNAs by protecting them from hydrolysis and is commonly used as tool to explore their function. It is, however, unclear what basic physical effects PT has on RNA stability and electronic structure. Here, we present Molecular Dynamics (MD) simulations, quantum mechanical (QM) calculations, and NMR spectroscopy measurements, exploring the effects of PT modifications in the structural context of the Neomycin-sensing riboswitch (NSR). The NSR is the smallest biologically functional riboswitch with a well-defined structure stabilized by a U-turn motif. Three of the signature interactions of the U-turn; an H-bond, an anion-π interaction and a potassium binding site; are formed by RNA phosphates, making the NSR an ideal model for studying how PT affects RNA structure and dynamics. By comparing with high-level QM calculations, we reveal the distinct physical properties of the individual interactions facilitated by the PT. The sulphur substitution, besides weakening the direct H-bond interaction, reduces the directionality of H-bonding while increasing its dispersion and induction components. It also reduces the induction and increases dispersion component of the anion-π stacking. The sulphur force-field parameters commonly employed in the literature do not reflect these distinctions, leading to unsatisfactory description of PT in simulations of the NSR. We show that it is not possible to accurately describe the PT interactions using one universal set of van der Waals sulphur parameters and provide suggestions for improving the force-field performance.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

MOLECULAR DYNAMICS SIMULATIONS OF HELIX BUNDLE PROTEINS USING UNRES FORCE FIELD AND ALL-ATOM FORCE FIELD

2012 ◽  
Vol 11 (06) ◽  
pp. 1201-1215 ◽  
Author(s):  
KAIFU GAO ◽  
MINGHUI YANG
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Protein A ◽  
Md Simulations ◽  
Time Frame ◽  
Coarse Grained ◽  
Staphylococcal Protein A ◽  
Villin Headpiece ◽  
Helix Bundle ◽  
Dynamics Simulations

We have investigated the folding of two helix-bundle proteins, 36-residue Villin headpiece and 56-residue E-domain of Staphylococcal protein A, by combining molecular dynamics (MD) simulations with Coarse-Grained United-Residue (UNRES) Force Field and all-atom force field. Starting from extended structures, each of the proteins was folded to a stable structure within a short time frame using the UNRES model. However, the secondary structures of helices were not well formed. Further refinement using MD simulations with the all-atom force field was able to fold the protein structure into the native-like state with the smallest main-chain root-mean-square deviation of around 3 Å. Detailed analysis of the folding trajectories was presented and the performance of GPU-based MD simulations was also discussed.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2019 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

scholarly journals Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

Predicting NMR Relaxation of Proteins from Molecular Dynamics Simulations with Accurate Methyl Rotation Barriers

2020 ◽  
Author(s):  
Falk Hoffmann ◽  
Frans Mulder ◽  
Lars V. Schäfer
Keyword(s):  
Molecular Dynamics ◽  
Force Field ◽  
Force Fields ◽  
Nmr Relaxation ◽  
Md Simulations ◽  
Quantitative Agreement ◽  
Relaxation Rates ◽  
Protein Force Fields ◽  
Dynamics Simulations ◽  
Methyl Rotation

The internal dynamics of proteins occurring on time scales from picoseconds to nanoseconds can be sensitively probed by nuclear magnetic resonance (NMR) spin relaxation experiments, as well as by molecular dynamics (MD) simulations. This complementarity offers unique opportunities, provided that the two methods are compared at a suitable level. Recently, several groups have used MD simulations to compute the spectral density of backbone and side-chain molecular motions, and to predict NMR relaxation rates from these. Unfortunately, in the case of methyl groups in protein side-chains, inaccurate energy barriers to methyl rotation were responsible for a systematic discrepancy in the computed relaxation rates, as demonstrated for the AMBER ff99SB*-ILDN force field (and related parameter sets), impairing quantitative agreement between simulations and experiments. However, correspondence could be regained by emending the MD force field with accurate coupled cluster quantum chemical calculations. Spurred by this positive result, we tested whether this approach could be generally applicable, in spite of the fact that different MD force fields employ different water models. Improved methyl group rotation barriers for the CHARMM36 and AMBER ff15ipq protein force fields were derived, such that the NMR relaxation data obtained from the MD simulations now also display very good agreement with experiment. Results herein showcase the performance of present-day MD force fields, and manifest their refined ability to accurately describe internal protein dynamics.

Assessing the Antimalarial Potentials of Phytochemicals: Virtual Screening, Molecular Dynamics and In-Vitro Investigations

2019 ◽  
Vol 16 (3) ◽  
pp. 291-300
Author(s):  
Saumya K. Patel ◽  
Mohd Athar ◽  
Prakash C. Jha ◽  
Vijay M. Khedkar ◽  
Yogesh Jasrai ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Integrity ◽  
Md Simulations ◽  
Tylophora Indica ◽  
Multidrug Resistance Protein 1 ◽  
Receptor Complexes ◽  
Resistance Protein ◽  
Dynamics Simulations ◽  
Stable Trajectory

Background: Combined in-silico and in-vitro approaches were adopted to investigate the antiplasmodial activity of Catharanthus roseus and Tylophora indica plant extracts as well as their isolated components (vinblastine, vincristine and tylophorine). </P><P> Methods: We employed molecular docking to prioritize phytochemicals from a library of 26 compounds against Plasmodium falciparum multidrug-resistance protein 1 (PfMDR1). Furthermore, Molecular Dynamics (MD) simulations were performed for a duration of 10 ns to estimate the dynamical structural integrity of ligand-receptor complexes. </P><P> Results: The retrieved bioactive compounds viz. tylophorine, vinblastin and vincristine were found to exhibit significant interacting behaviour; as validated by in-vitro studies on chloroquine sensitive (3D7) as well as chloroquine resistant (RKL9) strain. Moreover, they also displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations. </P><P> Conclusion: We anticipate that the retrieved phytochemicals can serve as the potential hits and presented findings would be helpful for the designing of malarial therapeutics.

scholarly journals Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau
Keyword(s):  
Carbon Dioxide ◽  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Hydrodynamic Radius ◽  
Md Simulations ◽  
Molecular Models ◽  
Experimental Values ◽  
Carbon Dioxide Co2 ◽  
Dynamics Simulations ◽  
Apparent Disagreement

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.

Sign in / Sign up

Export Citation Format

Share Document