Adaptive Markov State Model estimation using short reseeding trajectories

Abstract

In the last decade, advances in molecular dynamics (MD) and Markov State Model (MSM) methodologies have made possible accurate and efficient estimation of kinetic rates and reactive pathways for complex biomolecular dynamics occurring on slow timescales. A promising approach to enhanced sampling of MSMs is to use so-called "adaptive" methods, in which new MD trajectories are "seeded" preferentially from previously identified states. Here, we investigate the performance of various MSM estimators applied to reseeding trajectory data, for both a simple 1D free energy landscape, and for mini-protein folding MSMs of WW domain and NTL9(1-39). Our results reveal the practical challenges of reseeding simulations, and suggest a simple way to reweight seeding trajectory data to better estimate both thermodynamic and kinetic quantities.

References in corpus (7)

• Identification of slow molecular order parameters for Markov model construction
• Variational cross-validation of slow dynamical modes in molecular kinetics
• Statistically optimal analysis of state-discretized trajectory data from multiple thermodynamic states
• Markov State Models from short non-Equilibrium Simulations - Analysis and Correction of Estimation Bias
• Efficient estimation of rare-event kinetics
• The "weighted ensemble" path sampling method is statistically exact for a broad class of stochastic processes and binning procedures
• A maximum-caliber approach to predicting perturbed folding kinetics due to mutations

Cited by in corpus (3)

• Coarse-Grained Modelling Out of Equilibrium
• Iterative trajectory reweighting for estimation of equilibrium and non-equilibrium observables
• Unbiased estimation of equilibrium, rates, and committors from Markov state model analysis