Error analysis and efficient sampling in Markovian state models for molecular dynamics

被引:138
作者
Singhal, N
Pande, VS [1 ]
机构
[1] Stanford Univ, Dept Comp Sci, Stanford, CA 94305 USA
[2] Stanford Univ, Dept Chem, Stanford, CA 94305 USA
基金
美国国家科学基金会; 美国国家卫生研究院;
关键词
D O I
10.1063/1.2116947
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
In previous work, we described a Markovian state model (MSM) for analyzing molecular-dynamics trajectories, which involved grouping conformations into states and estimating the transition probabilities between states. In this paper, we analyze the errors in this model caused by finite sampling. We give different methods with various approximations to determine the precision of the reported mean first passage times. These approximations are validated on an 87 state toy Markovian system. In addition, we propose an efficient and practical sampling algorithm that uses these error calculations to build a MSM that has the same precision in mean first passage time values but requires an order of magnitude fewer samples. We also show how these methods can be scaled to large systems using sparse matrix methods.
引用
收藏
页数:13
相关论文
共 21 条
[1]   PARTITIONING, TEARING AND MODIFICATION OF SPARSE LINEAR-SYSTEMS [J].
BUNCH, JR ;
ROSE, DJ .
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS, 1974, 48 (02) :574-593
[2]   Tryptophan zippers:: Stable, monomeric β-hairpins [J].
Cochran, AG ;
Skelton, NJ ;
Starovasnik, MA .
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA, 2001, 98 (10) :5578-5583
[3]   PROTEIN UNFOLDING PATHWAYS EXPLORED THROUGH MOLECULAR-DYNAMICS SIMULATIONS [J].
DAGGETT, V ;
LEVITT, M .
JOURNAL OF MOLECULAR BIOLOGY, 1993, 232 (02) :600-619
[4]  
Devroye L., 1986, NONUNIFORM RANDOM VA
[5]  
Duff IS, 1986, DIRECT METHODS SPARS
[6]   Foldamer dynamics expressed via Markov state models. I. Explicit solvent molecular-dynamics simulations in acetonitrile, chloroform, methanol, and water [J].
Elmer, SP ;
Park, S ;
Pande, VS .
JOURNAL OF CHEMICAL PHYSICS, 2005, 123 (11)
[7]   Foldamer simulations: Novel computational methods and applications to poly-phenylacetylene oligomers [J].
Elmer, SP ;
Pande, VS .
JOURNAL OF CHEMICAL PHYSICS, 2004, 121 (24) :12760-12771
[8]  
GILL PE, 1974, MATH COMPUT, V28, P505, DOI 10.1090/S0025-5718-1974-0343558-6
[9]  
Golub G. H., 1996, MATRIX COMPUTATIONS
[10]  
HUNTER EJ, 1978, MILITARY FAMILIES AD