Constructing multi-resolution Markov State Models (MSMs) to elucidate RNA hairpin folding mechanisms.
Simulating biologically relevant timescales at atomic resolution is a challenging task since typical atomistic simulations are at least two orders of magnitude shorter. Markov State Models (MSMs) provide one means of overcoming this gap without sacrificing atomic resolution by extracting long time dynamics from short simulations. MSMs coarse grain space by dividing conformational space into long-lived, or metastable, states. This is equivalent to coarse graining time by integrating out fast motions within metastable states. By varying the degree of coarse graining one can vary the resolution of an MSM; therefore, MSMs are inherently multi-resolution. Here we introduce a new algorithm Super-level-set Hierarchical Clustering (SHC), to our knowledge, the first algorithm focused on constructing MSMs at multiple resolutions. The key insight of this algorithm is to generate a set of super levels covering different density regions of phase space, then cluster each super level separately, and finally recombine this information into a single MSM. SHC is able to produce MSMs at different resolutions using different super density level sets. To demonstrate the power of this algorithm we apply it to a small RNA hairpin, generating MSMs at four different resolutions. We validate these MSMs by showing that they are able to reproduce the original simulation data. Furthermore, long time folding dynamics are extracted from these models. The results show that there are no metastable on-pathway intermediate states. Instead, the folded state serves as a hub directly connected to multiple unfolded/misfolded states which are separated from each other by large free energy barriers.
Huang
Xuhui
X
Department of Chemistry, The Hong Kong University of Science & Technology, Kowloon, Hong Kong, China. xuhuihuang@gmail.com
Yao
Yuan
Y
Bowman
Gregory R
GR
Sun
Jian
J
Guibas
Leonidas J
LJ
Carlsson
Gunnar
G
Pande
Vijay S
VS
D000465:Algorithms; D019295:Computational Biology; D008390:Markov Chains; D008958:Models, Molecular; D056004:Molecular Dynamics Simulation; D009690:Nucleic Acid Conformation; D012313:RNA / Q000737:chemistry*; D059370:RNA Folding
D012313:RNA