Dynamics of an adenine-adenine RNA conformational switch from discrete path sampling.

The study of "rare event" dynamics can be challenging despite continuing advances in computer hardware. A wide variety of methods based on the master equation approach have been developed to tackle such problems, where the focus is on Markovian dynamics between appropriately defined states. In this contribution, we employ the discrete path sampling approach to characterize pathways and rates for an adenine-adenine RNA conformational switch. The underlying free energy landscape supports competing structures separated by relatively high barriers, with the two principal funnels leading to the major and minor conformations identified by NMR experiments. The interconversion time scale is predicted to be a few hundred seconds, consistent with the experimental lower bound estimates. We find that conformational switching occurs via stacked intermediates, through a sliding mechanism, in agreement with a previous simulation study. By retaining full dimensionality and avoiding low-dimensional projections, the mechanism can be described at an atomistic level of detail.

[1]  David J Wales,et al.  Energy Landscapes of Mini-Dumbbell DNA Octanucleotides. , 2018, Journal of chemical theory and computation.

[2]  Darrin M York,et al.  Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin , 2015, RNA.

[3]  David J Wales,et al.  Decoding the energy landscape: extracting structure, dynamics and thermodynamics , 2012, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[4]  M. Karplus,et al.  Hidden complexity of free energy surfaces for peptide (protein) folding. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[5]  D. Wales,et al.  Probing helical transitions in a DNA duplex. , 2017, Physical chemistry chemical physics : PCCP.

[6]  Robert M. Dirks,et al.  RNA force field with accuracy comparable to state-of-the-art protein force fields , 2018, Proceedings of the National Academy of Sciences.

[7]  D. Wales,et al.  A doubly nudged elastic band method for finding transition states. , 2004, The Journal of chemical physics.

[8]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[9]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[10]  Eric Vanden-Eijnden,et al.  Markovian milestoning with Voronoi tessellations. , 2009, The Journal of chemical physics.

[11]  Daniel M Zuckerman,et al.  The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures. , 2008, The Journal of chemical physics.

[12]  P. Beroza,et al.  Application of a pairwise generalized Born model to proteins and nucleic acids: inclusion of salt effects , 1999 .

[13]  David J. Wales,et al.  Symmetrization of the AMBER and CHARMM force fields , 2010, J. Comput. Chem..

[14]  Vijay S. Pande,et al.  Screen Savers of the World Unite! , 2000, Science.

[15]  Hashim M. Al-Hashimi,et al.  Functional complexity and regulation through RNA dynamics , 2012, Nature.

[16]  R. Stephen Berry,et al.  Statistical interpretation of topographies and dynamics of multidimensional potentials , 1995 .

[17]  Angel E García,et al.  High-resolution reversible folding of hyperstable RNA tetraloops using molecular dynamics simulations , 2013, Proceedings of the National Academy of Sciences.

[18]  Andrés Marzal,et al.  Computing the K Shortest Paths: A New Algorithm and an Experimental Comparison , 1999, WAE.

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

[20]  Scott D Kennedy,et al.  An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA. , 2006, Biochemistry.

[21]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[22]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[23]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.

[24]  J. Doye,et al.  THE DOUBLE-FUNNEL ENERGY LANDSCAPE OF THE 38-ATOM LENNARD-JONES CLUSTER , 1998, cond-mat/9808265.

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

[26]  Giovanni Bussi,et al.  Predicting the Kinetics of RNA Oligonucleotides Using Markov State Models. , 2016, Journal of chemical theory and computation.

[27]  D. Wales,et al.  Energy landscapes, folding mechanisms, and kinetics of RNA tetraloop hairpins. , 2014, Journal of the American Chemical Society.

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

[29]  Richard A. Cunha,et al.  RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview , 2018, Chemical reviews.

[30]  Graeme Henkelman,et al.  Unification of algorithms for minimum mode optimization. , 2014, The Journal of chemical physics.

[31]  A. Garcia,et al.  Simulation of the pressure and temperature folding/unfolding equilibrium of a small RNA hairpin. , 2008, Journal of the American Chemical Society.

[32]  David J Wales,et al.  Calculating rate constants and committor probabilities for transition networks by graph transformation. , 2009, The Journal of chemical physics.

[33]  David J. Wales,et al.  Exploring biomolecular energy landscapes. , 2017, Chemical communications.

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

[35]  David J. Wales,et al.  The Energy Landscape as a Computational Tool , 2008 .

[36]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[37]  David J. Wales,et al.  Some further applications of discrete path sampling to cluster isomerization , 2004 .

[38]  Peter G. Bolhuis,et al.  A novel path sampling method for the calculation of rate constants , 2003 .

[39]  Sergei V. Krivov,et al.  Free energy disconnectivity graphs: Application to peptide models , 2002 .

[40]  Frank Noé,et al.  Markov models of molecular kinetics: generation and validation. , 2011, The Journal of chemical physics.

[41]  Mark A. Miller,et al.  Archetypal energy landscapes , 1998, Nature.

[42]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[43]  R. Elber,et al.  Computing time scales from reaction coordinates by milestoning. , 2004, The Journal of chemical physics.

[44]  D. Wales,et al.  Energy Landscape and Pathways for Transitions between Watson-Crick and Hoogsteen Base Pairing in DNA. , 2018, The journal of physical chemistry letters.

[45]  D. Frenkel,et al.  Simulating rare events in equilibrium or nonequilibrium stochastic systems. , 2005, The Journal of chemical physics.

[46]  M. Hoare,et al.  Statistical mechanics and morphology of very small atomic clusters , 1976 .

[47]  J. Šponer,et al.  Refinement of the Cornell et al. Nucleic Acids Force Field Based on Reference Quantum Chemical Calculations of Glycosidic Torsion Profiles , 2011, Journal of chemical theory and computation.

[48]  E. Westhof,et al.  Geometric nomenclature and classification of RNA base pairs. , 2001, RNA.

[49]  David J Wales,et al.  Communication: Analysing kinetic transition networks for rare events. , 2014, The Journal of chemical physics.

[50]  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.

[51]  Lindsey J. Munro,et al.  DEFECT MIGRATION IN CRYSTALLINE SILICON , 1999 .

[52]  R. MacKay,et al.  Aggregation of Markov flows I: theory , 2018, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

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

[54]  D. Thirumalai,et al.  Kinetics of Folding of Proteins and RNA , 1996 .

[55]  David J Wales,et al.  Transforming the Energy Landscape of a Coiled-Coil Peptide via Point Mutations. , 2017, Journal of chemical theory and computation.

[56]  F. Noé,et al.  Transition networks for modeling the kinetics of conformational change in macromolecules. , 2008, Current opinion in structural biology.

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

[58]  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.

[59]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[60]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[61]  M. Parrinello,et al.  From metadynamics to dynamics. , 2013, Physical review letters.

[62]  Michal Otyepka,et al.  Computer Folding of RNA Tetraloops? Are We There Yet? , 2013, Journal of chemical theory and computation.

[63]  Ron Elber,et al.  Extending molecular dynamics time scales with milestoning: example of complex kinetics in a solvated peptide. , 2007, The Journal of chemical physics.

[64]  David J Wales,et al.  Folding pathways and rates for the three-stranded beta-sheet peptide Beta3s using discrete path sampling. , 2008, The journal of physical chemistry. B.

[65]  R. S. Berry,et al.  USING CLUSTERS TO RELATE TOPOGRAPHY AND DYNAMICS OF MULTIDIMENSIONAL POTENTIALS , 1996 .

[66]  Ron Elber,et al.  Kinetics of helix unfolding: molecular dynamics simulations with milestoning. , 2009, The journal of physical chemistry. A.

[67]  D. Wales Discrete path sampling , 2002 .

[68]  V. Pande,et al.  Markov State Models: From an Art to a Science. , 2018, Journal of the American Chemical Society.

[69]  Ron Elber,et al.  On the assumptions underlying milestoning. , 2008, The Journal of chemical physics.

[70]  David J Wales,et al.  Thermodynamics and kinetics of aggregation for the GNNQQNY peptide. , 2007, Journal of the American Chemical Society.

[71]  Berry,et al.  Topography and Dynamics of Multidimensional Interatomic Potential Surfaces. , 1995, Physical review letters.

[72]  Harry A. Stern,et al.  Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs , 2011, Journal of chemical theory and computation.

[73]  D. Thirumalai,et al.  RNA and protein folding: common themes and variations. , 2005, Biochemistry.

[74]  Ron Elber,et al.  Calculating Iso-Committor Surfaces as Optimal Reaction Coordinates with Milestoning , 2017, Entropy.

[75]  David J Wales,et al.  Improving Computational Predictions of Single-Stranded RNA Tetramers with Revised α/γ Torsional Parameters for the Amber Force Field. , 2017, The journal of physical chemistry. B.

[76]  D Thirumalai,et al.  How do metal ions direct ribozyme folding? , 2015, bioRxiv.

[77]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[78]  Peter Salamon,et al.  Observation time scale, free-energy landscapes, and molecular symmetry , 2013, Proceedings of the National Academy of Sciences.

[79]  A. J. Ballard,et al.  Intrinsically Disordered Energy Landscapes , 2015, Scientific Reports.

[80]  K. Laidler,et al.  Symmetries of activated complexes , 1968 .

[81]  Karl Heinz Hoffmann,et al.  Atomic clusters and nanoscale particles: From coarse-grained dynamics to optimized annealing schedules , 1998 .

[82]  Diego Prada-Gracia,et al.  Exploring the Free Energy Landscape: From Dynamics to Networks and Back , 2009, PLoS Comput. Biol..

[83]  M. Karplus,et al.  The topology of multidimensional potential energy surfaces: Theory and application to peptide structure and kinetics , 1997 .

[84]  Vijay S Pande,et al.  Using path sampling to build better Markovian state models: predicting the folding rate and mechanism of a tryptophan zipper beta hairpin. , 2004, The Journal of chemical physics.

[85]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[86]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[87]  Changbong Hyeon,et al.  RNA under tension: Folding Landscapes, Kinetic partitioning Mechanism, and Molecular Tensegrity. , 2012, The journal of physical chemistry letters.

[88]  David J Wales,et al.  Energy landscapes: some new horizons. , 2010, Current opinion in structural biology.

[89]  David J Wales,et al.  Refined kinetic transition networks for the GB1 hairpin peptide. , 2009, Physical chemistry chemical physics : PCCP.

[90]  D. L. Freeman,et al.  Phase changes in 38-atom Lennard-Jones clusters. I. A parallel tempering study in the canonical ensemble , 2000, physics/0003068.

[91]  J. Šponer,et al.  Refinement of the AMBER Force Field for Nucleic Acids: Improving the Description of α/γ Conformers , 2007 .

[92]  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.

[93]  David J Wales,et al.  Finding pathways between distant local minima. , 2005, The Journal of chemical physics.

[94]  Scott D Kennedy,et al.  Molecular Mechanics Investigation of an Adenine-Adenine Non-Canonical Pair Conformational Change. , 2011, Journal of chemical theory and computation.

[95]  C. Dellago,et al.  Transition path sampling and the calculation of rate constants , 1998 .

[96]  Gerhard Hummer,et al.  Peptide folding kinetics from replica exchange molecular dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[97]  C. Brooks,et al.  Hierarchy of RNA functional dynamics. , 2014, Annual review of biochemistry.

[98]  David Chandler,et al.  Transition path sampling: throwing ropes over rough mountain passes, in the dark. , 2002, Annual review of physical chemistry.

[99]  Samuela Pasquali,et al.  Multifunctional energy landscape for a DNA G-quadruplex: An evolved molecular switch. , 2017, The Journal of chemical physics.