Free Energy Landscape of GAGA and UUCG RNA Tetraloops

Abstract

We report the folding thermodynamics of ccUUCGgg and ccGAGAgg RNA tetraloops using atomistic molecular dynamics simulations. We obtain a previously unreported estimation of the folding free energy using parallel tempering in combination with well-tempered metadynamics. A key ingredient is the use of a recently developed metric distance, eRMSD, as a biased collective variable. We find that the native fold of both tetraloops is not the global free energy minimum using the Amber\c{hi}OL3 force field. The estimated folding free energies are 30.2kJ/mol for UUCG and 7.5 kJ/mol for GAGA, in striking disagreement with experimental data. We evaluate the viability of all possible one-dimensional backbone force field corrections. We find that disfavoring the gauche+ region of α and ζ angles consistently improves the existing force field. The level of accuracy achieved with these corrections, however, cannot be considered sufficient by judging on the basis of available thermodynamic data and solution experiments.

arXiv commentJournal of Physical Chemistry Letters (2016)

References in corpus (9)

• Canonical sampling through velocity-rescaling
• Escaping free-energy minima
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• Metadynamics with adaptive Gaussians
• The Role of Nucleobase Interactions in RNA Structure and Dynamics
• Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies
• Empirical corrections to the Amber RNA force field with Target Metadynamics
• RNA Folding Pathways in Stop Motion

Cited by in corpus (12)

• Using the Maximum Entropy Principle to Combine Simulations and Solution Experiments
• Combining simulations and solution experiments as a paradigm for RNA force field refinement
• Exploring RNA structure and dynamics through enhanced sampling simulations
• Analyzing and biasing simulations with PLUMED
• A nucleobase-centered coarse-grained representation for structure prediction of RNA motifs
• Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles
• Predicting the kinetics of RNA oligonucleotides using Markov state models
• The mechanism of RNA base fraying: molecular dynamics simulations analyzed with core-set Markov state models
• Conformational Entropy as Collective Variable for Proteins
• Effects and limitations of a nucleobase-driven backmapping procedure for nucleic acids using steered Molecular Dynamics
• Boosting ensemble refinement with transferable force field corrections: synergistic optimization for molecular simulations
• Automatic learning of hydrogen-bond fixes in an AMBER RNA force field