The mechanism of RNA base fraying: molecular dynamics simulations analyzed with core-set Markov state models
arXiv:1811.12144 · doi:10.1063/1.5083227
Abstract
The process of RNA base fraying (i.e. the transient opening of the termini of a helix) is involved in many aspects of RNA dynamics. We here use molecular dynamics simulations and Markov state models to characterize the kinetics of RNA fraying and its sequence and direction dependence. In particular, we first introduce a method for determining biomolecular dynamics employing core-set Markov state models constructed using an advanced clustering technique. The method is validated on previously reported simulations. We then use the method to analyze extensive trajectories for four different RNA model duplexes. Results obtained using D. E. Shaw research and AMBER force fields are compared and discussed in detail, and show a non-trivial interplay between the stability of intermediate states and the overall fraying kinetics.
arXiv commentSupporting information included in ancillary files
References in corpus (12)
- Canonical sampling through velocity-rescaling
- Identification of slow molecular order parameters for Markov model construction
- Estimating the intrinsic dimension of datasets by a minimal neighborhood information
- The Role of Nucleobase Interactions in RNA Structure and Dynamics
- Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies
- Free Energy Landscape of GAGA and UUCG RNA Tetraloops
- Unravelling Mg-RNA binding with atomistic molecular dynamics
- Exploring RNA structure and dynamics through enhanced sampling simulations
- 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
- Predicting the kinetics of RNA oligonucleotides using Markov state models