Exploring RNA structure and dynamics through enhanced sampling simulations.

RNA function is intimately related to its structural dynamics. Molecular dynamics simulations are useful for exploring biomolecular flexibility but are severely limited by the accessible timescale. Enhanced sampling methods allow this timescale to be effectively extended in order to probe biologically relevant conformational changes and chemical reactions. Here, we review the role of enhanced sampling techniques in the study of RNA systems. We discuss the challenges and promises associated with the application of these methods to force-field validation, exploration of conformational landscapes and ion/ligand-RNA interactions, as well as catalytic pathways. Important technical aspects of these methods, such as the choice of the biased collective variables and the analysis of multi-replica simulations, are examined in detail. Finally, a perspective on the role of these methods in the characterization of RNA dynamics is provided.

[1]  Darrin M. York,et al.  Ribozyme Catalysis with a Twist: Active State of the Twister Ribozyme in Solution Predicted from Molecular Simulation. , 2016, Journal of the American Chemical Society.

[2]  Stefano Piana,et al.  Demonstrating an Order-of-Magnitude Sampling Enhancement in Molecular Dynamics Simulations of Complex Protein Systems. , 2016, Journal of chemical theory and computation.

[3]  Sharon Hammes-Schiffer,et al.  Role of the Active Site Guanine in the glmS Ribozyme Self-Cleavage Mechanism: Quantum Mechanical/Molecular Mechanical Free Energy Simulations , 2014, Journal of the American Chemical Society.

[4]  Jihua Wang,et al.  Ligand Selectivity Mechanism and Conformational Changes in Guanine Riboswitch by Molecular Dynamics Simulations and Free Energy Calculations , 2017, J. Chem. Inf. Model..

[5]  Hamed S. Hayatshahi,et al.  Computational Assessment of Potassium and Magnesium Ion Binding to a Buried Pocket in GTPase-Associating Center RNA , 2016, The journal of physical chemistry. B.

[6]  D. Marx,et al.  Pressure modulates the self-cleavage step of the hairpin ribozyme , 2017, Nature Communications.

[7]  Yuan-yan Wu,et al.  Multivalent ion-mediated nucleic acid helix-helix interactions: RNA versus DNA , 2015, Nucleic acids research.

[8]  Paolo Carloni,et al.  Molecular view of ligands specificity for CAG repeats in anti-Huntington therapy. , 2015, Journal of chemical theory and computation.

[9]  Sharon Hammes-Schiffer,et al.  The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme. , 2017, Nature chemical biology.

[10]  Manho Lim,et al.  Predicting RNA Structures via a Simple van der Waals Correction to an All-Atom Force Field. , 2017, Journal of chemical theory and computation.

[11]  Eric Westhof.  Twenty years of RNA crystallography. , 2015, RNA.

[12]  Alexander D. MacKerell,et al.  Characterization of Mg2+ Distributions around RNA in Solution , 2016, ACS omega.

[13]  Peter C Anderson,et al.  Molecular simulations and Markov state modeling reveal the structural diversity and dynamics of a theophylline-binding RNA aptamer in its unbound state , 2017, PloS one.

[14]  Harald Schwalbe,et al.  Mapping the landscape of RNA dynamics with NMR spectroscopy. , 2011, Accounts of chemical research.

[15]  L. Kay,et al.  Simultaneous NMR characterisation of multiple minima in the free energy landscape of an RNA UUCG tetraloop. , 2017, Physical chemistry chemical physics : PCCP.

[16]  Angel E García,et al.  Free-energy landscape of a hyperstable RNA tetraloop , 2016, Proceedings of the National Academy of Sciences.

[17]  Giovanni Bussi,et al.  Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies. , 2016, Journal of chemical theory and computation.

[18]  Giovanni Bussi,et al.  Kissing loop interaction in adenine riboswitch: insights from umbrella sampling simulations , 2015, BMC Bioinformatics.

[19]  A. Serganov,et al.  A Decade of Riboswitches , 2013, Cell.

[20]  Giovanni Bussi,et al.  Enhanced Conformational Sampling Using Replica Exchange with Collective-Variable Tempering , 2015, Journal of chemical theory and computation.

[21]  Carlo Camilloni,et al.  Structure of a low-population binding intermediate in protein-RNA recognition , 2016, Proceedings of the National Academy of Sciences.

[22]  Jia Sheng,et al.  Thermodynamic insights into 2-thiouridine-enhanced RNA hybridization , 2015, Nucleic acids research.

[23]  The structure and folding of kink turns in RNA , 2012, Wiley interdisciplinary reviews. RNA.

[24]  K. Asai,et al.  Predicting RNA Duplex Dimerization Free-Energy Changes upon Mutations Using Molecular Dynamics Simulations. , 2015, The journal of physical chemistry letters.

[25]  Giovanni Bussi,et al.  Unraveling Mg2+–RNA binding with atomistic molecular dynamics , 2016, RNA.

[26]  Wei Wang,et al.  Free Energy Landscape and Multiple Folding Pathways of an H-Type RNA Pseudoknot , 2015, PloS one.

[27]  Michal Otyepka,et al.  Insights into Stability and Folding of GNRA and UNCG Tetraloops Revealed by Microsecond Molecular Dynamics and Well-Tempered Metadynamics. , 2015, Journal of chemical theory and computation.

[28]  David H Mathews,et al.  Physics‐based all‐atom modeling of RNA energetics and structure , 2017, Wiley interdisciplinary reviews. RNA.

[29]  Niel M. Henriksen,et al.  Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields , 2015, RNA.

[30]  D. Lilley.  How RNA acts as a nuclease: some mechanistic comparisons in the nucleolytic ribozymes. , 2017, Biochemical Society transactions.

[31]  L. Zhihong,et al.  Folding of SAM-II riboswitch explored by replica-exchange molecular dynamics simulation. , 2015, Journal of theoretical biology.

[32]  K. Schulten,et al.  Molecular Mechanism of Processive 3' to 5' RNA Translocation in the Active Subunit of the RNA Exosome Complex. , 2016, Journal of the American Chemical Society.

[33]  C. Roland,et al.  Structure and Dynamics of DNA and RNA Double Helices of CAG and GAC Trinucleotide Repeats. , 2017, Biophysical journal.

[34]  Zhaoxi Sun,et al.  Protonation-dependent base flipping in the catalytic triad of a small RNA , 2017 .

[35]  Giovanni Bussi,et al.  Empirical Corrections to the Amber RNA Force Field with Target Metadynamics , 2016, Journal of chemical theory and computation.

[36]  David H. Mathews,et al.  Revised RNA Dihedral Parameters for the Amber Force Field Improve RNA Molecular Dynamics , 2017, Journal of chemical theory and computation.

[37]  Ming Huang,et al.  Nucleic acid reactivity: Challenges for next‐generation semiempirical quantum models , 2015, J. Comput. Chem..

[38]  Michele Parrinello,et al.  Enhancing Important Fluctuations: Rare Events and Metadynamics from a Conceptual Viewpoint. , 2016, Annual review of physical chemistry.

[39]  Sharon Hammes-Schiffer,et al.  Inverse Thio Effects in the Hepatitis Delta Virus Ribozyme Reveal that the Reaction Pathway Is Controlled by Metal Ion Charge Density , 2015, Biochemistry.

[40]  Darrin M. York,et al.  Comparison of structural, thermodynamic, kinetic and mass transport properties of Mg2+ ion models commonly used in biomolecular simulations , 2015, J. Comput. Chem..

[41]  Kyle E. Watters,et al.  Using in-cell SHAPE-Seq and simulations to probe structure–function design principles of RNA transcriptional regulators , 2016, RNA.

[42]  Giovanni Bussi,et al.  Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides. , 2017, RNA.

[43]  Tai-Sung Lee,et al.  Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM , 2015, RNA.

[44]  Michal Otyepka,et al.  How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists. , 2013, Methods.

[45]  Thomas E Cheatham,et al.  Improved Force Field Parameters Lead to a Better Description of RNA Structure. , 2015, Journal of chemical theory and computation.

[46]  T. Bandyopadhyay,et al.  Water isotope effect on the thermostability of a polio viral RNA hairpin: A metadynamics study. , 2017, The Journal of chemical physics.

[47]  Alessandra Magistrato,et al.  Who Activates the Nucleophile in Ribozyme Catalysis? An Answer from the Splicing Mechanism of Group II Introns. , 2016, Journal of the American Chemical Society.

[48]  Giovanni Bussi,et al.  Free Energy Landscape of GAGA and UUCG RNA Tetraloops. , 2016, The journal of physical chemistry letters.

[49]  Tsuyoshi Murata,et al.  {m , 1934, ACML.

[50]  S. Ranganathan,et al.  Advances in RNA molecular dynamics: a simulator's guide to RNA force fields , 2017, Wiley interdisciplinary reviews. RNA.

[51]  J. L. Suter,et al.  Structure, dynamics, and function of the hammerhead ribozyme in bulk water and at a clay mineral surface from replica exchange molecular dynamics. , 2015, Langmuir : the ACS journal of surfaces and colloids.

[52]  Darrin M York,et al.  Divalent Metal Ion Activation of a Guanine General Base in the Hammerhead Ribozyme: Insights from Molecular Simulations. , 2017, Biochemistry.

[53]  Giovanni Bussi,et al.  Combining Simulations and Solution Experiments as a Paradigm for RNA Force Field Refinement. , 2016, Journal of chemical theory and computation.

[54]  K. Morris,et al.  The rise of regulatory RNA , 2014, Nature Reviews Genetics.

[55]  Sharon Hammes-Schiffer,et al.  Assessing the Potential Effects of Active Site Mg2+ Ions in the glmS Ribozyme–Cofactor Complex , 2016, The journal of physical chemistry letters.

[56]  Angel E García,et al.  Equilibrium Denaturation and Preferential Interactions of an RNA Tetraloop with Urea. , 2017, The journal of physical chemistry. B.

[57]  Rafael C. Bernardi,et al.  Enhanced sampling techniques in molecular dynamics simulations of biological systems. , 2015, Biochimica et biophysica acta.

[58]  Shantenu Jha,et al.  Characterization of the three-dimensional free energy manifold for the uracil ribonucleoside from asynchronous replica exchange simulations. , 2015, Journal of chemical theory and computation.

[59]  T. Cheatham,et al.  Stem-Loop V of Varkud Satellite RNA Exhibits Characteristics of the Mg2+ Bound Structure in the Presence of Monovalent Ions , 2015, The journal of physical chemistry. B.

[60]  Asaminew H. Aytenfisu,et al.  Structural analysis of a class III preQ1 riboswitch reveals an aptamer distant from a ribosome-binding site regulated by fast dynamics , 2015, Proceedings of the National Academy of Sciences.

[61]  J Andrew McCammon,et al.  CRISPR-Cas9 conformational activation as elucidated from enhanced molecular simulations , 2017, Proceedings of the National Academy of Sciences.

[62]  Hamed S. Hayatshahi,et al.  Investigating the ion dependence of the first unfolding step of GTPase-Associating Center ribosomal RNA , 2018, Journal of biomolecular structure & dynamics.

[63]  Gianvito Grasso,et al.  Free energy landscape of siRNA-polycation complexation: Elucidating the effect of molecular geometry, polymer flexibility, and charge neutralization , 2017, PloS one.

[64]  Ilyas Yildirim,et al.  Crystallographic and Computational Analyses of AUUCU Repeating RNA That Causes Spinocerebellar Ataxia Type 10 (SCA10). , 2015, Biochemistry.

[65]  Giovanni Bussi,et al.  RNA folding pathways in stop motion , 2016, Nucleic acids research.

[66]  R. Corradini,et al.  Focus on PNA Flexibility and RNA Binding using Molecular Dynamics and Metadynamics , 2017, Scientific Reports.

[67]  George C. Schatz,et al.  Computational Investigation of RNA CUG Repeats Responsible for Myotonic Dystrophy 1 , 2015, Journal of chemical theory and computation.

[68]  Michal Otyepka,et al.  How to understand atomistic molecular dynamics simulations of RNA and protein–RNA complexes? , 2017, Wiley interdisciplinary reviews. RNA.

[69]  H. Al‐Hashimi,et al.  Visualizing Transient Low-Populated Structures of RNA , 2012, Nature.

[70]  Adam Roth,et al.  Ribozyme speed limits. , 2003, RNA.

[71]  Darrin M York,et al.  Force Field for Mg(2+), Mn(2+), Zn(2+), and Cd(2+) Ions That Have Balanced Interactions with Nucleic Acids. , 2015, The journal of physical chemistry. B.