RNA Folding Pathways in Stop Motion

Abstract

We introduce a method for predicting RNA folding pathways, with an application to the most important RNA tetraloops. The method is based on the idea that ensembles of three-dimensional fragments extracted from high-resolution crystal structures are heterogeneous enough to describe metastable as well as intermediate states. These ensembles are first validated by performing a quantitative comparison against available solution NMR data of a set of RNA tetranucleotides. Notably, the agreement is better with respect to the one obtained by comparing NMR with extensive all-atom molecular dynamics simulations. We then propose a procedure based on diffusion maps and Markov models that makes it possible to obtain reaction pathways and their relative probabilities from fragment ensembles. This approach is applied to study the helix-to-loop folding pathway of all the tetraloops from the GNRA and UNCG families. The results give detailed insights into the folding mechanism that are compatible with available experimental data and clarify the role of intermediate states observed in previous simulation studies. The method is computationally inexpensive and can be used to study arbitrary conformational transitions.

Accepted for publication on Nucleic Acids Research

Cited by in corpus (10)

• Folding of guanine quadruplex molecules -- funnel-like mechanism or kinetic partitioning? An overview from MD simulation studies
• Combining simulations and solution experiments as a paradigm for RNA force field refinement
• Free Energy Landscape of GAGA and UUCG RNA Tetraloops
• Unravelling Mg${}^{2+}$-RNA binding with atomistic molecular dynamics
• Exploring RNA structure and dynamics through enhanced sampling simulations
• Empirical corrections to the Amber RNA force field with Target Metadynamics
• 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
• Understanding In-line Probing Experiments by Modeling Cleavage of Non-reactive RNA Nucleotides
• Automatic learning of hydrogen-bond fixes in an AMBER RNA force field