Impact of 2′‐hydroxyl sampling on the conformational properties of RNA: Update of the CHARMM all‐atom additive force field for RNA
Top Cited Papers
- 5 April 2011
- journal article
- research article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 32 (9), 1929-1943
- https://doi.org/10.1002/jcc.21777
Abstract
Here, we present an update of the CHARMM27 all‐atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson‐Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2′‐hydroxyl moiety to almost exclusively sample the O3′ orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3′ orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2′‐hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3′ orientations of the 2′‐hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J‐coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2′‐hydroxyl moiety and support a model whereby the 2′‐hydroxyl can enhance the probability of conformational transitions in RNA. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011Keywords
This publication has 95 references indexed in Scilit:
- Toward a Consensus View of Duplex RNA FlexibilityBiophysical Journal, 2010
- Free Energy Profile of RNA Hairpins: A Molecular Dynamics Simulation StudyBiophysical Journal, 2010
- Crystallographic Studies of Chemically Modified Nucleic Acids: A Backward GlanceChemistry & Biodiversity, 2010
- CHARMM: The biomolecular simulation programJournal of Computational Chemistry, 2009
- Nonlocal Helix Formation Is Key to Understanding S-Adenosylmethionine-1 Riboswitch FunctionBiophysical Journal, 2009
- Additive empirical force field for hexopyranose monosaccharidesJournal of Computational Chemistry, 2008
- Molecular dynamics simulations of RNA: An in silico single molecule approachPeptide Science, 2006
- Scalable molecular dynamics with NAMDJournal of Computational Chemistry, 2005
- Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulationsJournal of Computational Chemistry, 2004
- Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparisonProteins: Structure, Function, and Bioinformatics, 1988