Hierarchical Nyström methods for constructing Markov state models for conformational dynamics

Yuan Yao, Raymond Z. Cui, Gregory R. Bowman, Daniel Adriano Silva, Jian Sun, Xuhui Huang

Research output: Contribution to journalArticlepeer-review

37 Scopus citations

Abstract

Markov state models (MSMs) have become a popular approach for investigating the conformational dynamics of proteins and other biomolecules. MSMs are typically built from numerous molecular dynamics simulations by dividing the sampled configurations into a large number of microstates based on geometric criteria. The resulting microstate model can then be coarse-grained into a more understandable macrostate model by lumping together rapidly mixing microstates into larger, metastable aggregates. However, finite sampling often results in the creation of many poorly sampled microstates. During coarse-graining, these states are mistakenly identified as being kinetically important because transitions to/from them appear to be slow. In this paper, we propose a formalism based on an algebraic principle for matrix approximation, i.e., the Nyström method, to deal with such poorly sampled microstates. Our scheme builds a hierarchy of microstates from high to low populations and progressively applies spectral clustering on sets of microstates within each level of the hierarchy. It helps spectral clustering identify metastable aggregates with highly populated microstates rather than being distracted by lowly populated states. We demonstrate the ability of this algorithm to discover the major metastable states on two model systems, the alanine dipeptide and trpzip2 peptide.

Original languageEnglish
Article number174106
JournalJournal of Chemical Physics
Volume138
Issue number17
DOIs
StatePublished - May 7 2013

Fingerprint

Dive into the research topics of 'Hierarchical Nyström methods for constructing Markov state models for conformational dynamics'. Together they form a unique fingerprint.

Cite this