Predicting the kinetics of RNA oligonucleotides using Markov state models
arXiv:1612.07470 · doi:10.1021/acs.jctc.6b00982
Abstract
Nowadays different experimental techniques, such as single molecule or relaxation experiments, can provide dynamic properties of biomolecular systems, but the amount of detail obtainable with these methods is often limited in terms of time or spatial resolution. Here we use state-of-the-art computational techniques, namely atomistic molecular dynamics and Markov state models, to provide insight into the rapid dynamics of short RNA oligonucleotides, in order to elucidate the kinetics of stacking interactions. Analysis of multiple microsecond-long simulations indicates that the main relaxation modes of such molecules can consist of transitions between alternative folded states, rather than between random coils and native structures. After properly removing structures that are artificially stabilized by known inaccuracies of the current RNA AMBER force field, the kinetic properties predicted are consistent with the timescales of previously reported relaxation experiments.
arXiv commentAccepted for publication on J. Chem. Theory Comput
References in corpus (11)
- Canonical sampling through velocity-rescaling
- Identification of slow molecular order parameters for Markov model construction
- The Role of Nucleobase Interactions in RNA Structure and Dynamics
- Projected and Hidden Markov Models for calculating kinetics and metastable states of complex molecules
- Combining simulations and solution experiments as a paradigm for RNA force field refinement
- Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering
- Free Energy Landscape of GAGA and UUCG RNA Tetraloops
- Empirical corrections to the Amber RNA force field with Target Metadynamics
- Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments
- RNA Folding Pathways in Stop Motion
- RNA unwinding from reweighted pulling simulations