A Bayesian method for construction of Markov models to describe dynamics on various time-scales

Journal of Chemical Physics

Volumn 133, Issue 14, 2010, Pages

  Rains, Emily K ^a   Andersen, Hans C ^a  

^a STANFORD UNIVERSITY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

ALANINE DIPEPTIDES; ATOMISTIC SIMULATIONS; BAYESIAN ANALYSIS; BAYESIAN METHODS; BAYESIAN PROBABILITIES; BIOLOGICAL PROCESS; COMPUTATIONAL BARRIERS; CONFIGURATION SPACE; DIAGONAL ELEMENTS; MARKOV MODEL; MARKOV STATE MODELS; MESOSCOPICS; MOLECULAR DYNAMICS SIMULATIONS; OPTIMIZATION ALGORITHMS; RANDOM WALK; SIMULATION DATA; TEST SYSTEMS; TIME STEP; TIME-SCALES; TRANSITION PROBABILITIES; TRANSITION RATE MATRIX; TWO-DIMENSION;

AMINO ACIDS; BAYESIAN NETWORKS; DYNAMICS; MARKOV PROCESSES; MOLECULAR DYNAMICS; PROBABILITY DISTRIBUTIONS; SUPERCONDUCTING MATERIALS;

COMPUTER SIMULATION;

OLIGOPEPTIDE; PROTEIN;

ALGORITHM; ARTICLE; ARTIFACT; BAYES THEOREM; BIOLOGICAL MODEL; CHEMICAL STRUCTURE; CHEMISTRY; COMPUTER SIMULATION; MOLECULAR DYNAMICS; PROBABILITY; TIME;

ALGORITHMS; ARTIFACTS; BAYES THEOREM; COMPUTER SIMULATION; MARKOV CHAINS; MODELS, BIOLOGICAL; MODELS, MOLECULAR; MOLECULAR DYNAMICS SIMULATION; OLIGOPEPTIDES; PROBABILITY; PROTEINS; TIME FACTORS;

EID: 77958103362     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.3496438     Document Type: Article

References (48)

1. J. D. Chodera, N. Singhal, V. S. Pande, K. A. Dill, and W. C. Swope, J. Chem. Phys. JCPSA6 0021-9606 126, 155101 (2007). 10.1063/1.2714538 (Pubitemid 46644960)
2. N. Singhal, C. D. Snow, and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 121, 415 (2004). 10.1063/1.1738647
3. S. Muff and A. Catflisch, Proteins 18PLAF 0887-3585 70, 1185 (2008). 10.1002/prot.21565 (Pubitemid 351304079)
4. A. M. A. West, R. Elber, and D. Shalloway, J. Phys. Chem. JPCHAX 0022-3654 127, 145104 (2007).
5. A. F. Voter, Phys. Rev. B PLRBAQ 0556-2805 57, R13985 (1998). 10.1103/PhysRevB.57.R13985
6. A. F. Voter, J. Chem. Phys. JCPSA6 0021-9606 106, 4665 (1996).
7. G. Hummer and I. G. Kevrekidis, J. Chem. Phys. JCPSA6 0021-9606 118, 10762 (2003). 10.1063/1.1574777
8. M. R. Sørensen and A. F. Voter, J. Chem. Phys. JCPSA6 0021-9606 112, 9599 (2000).
9. A. Laio and M. Parrinello, Proc. Natl. Acad. Sci. U.S.A. PNASA6 0027-8424 99, 12562 (2002). 10.1073/pnas.202427399
10. H. Grubmüller, Phys. Rev. E PLEEE8 1063-651X 52, 2893 (1995). 10.1103/PhysRevE.52.2893
11. J. Doye and D. Wales, Z. Phys. D: At., Mol. Clusters ZDACE2 0178-7683 40, 194 (1997).
12. D. J. Wales, Mol. Phys. MOPHAM 0026-8976 100, 3285 (2002). 10.1080/00268970210162691 (Pubitemid 35446132)
13. D. Moroni, T. S. van Erp, and P. G. Bolhuis, Physica A PHYADX 0378-4371 340, 395 (2002).
14. G. Henkelman and H. Jónsson, J. Chem. Phys. JCPSA6 0021-9606 113, 9978 (2000). 10.1063/1.1323224 (Pubitemid 32076892)
15. G. Barkema and N. Mousseau, Phys. Rev. Lett. PRLTAO 0031-9007 77, 4358 (1996). 10.1103/PhysRevLett.77.4358
16. J. D. Chodera, W. C. Swope, J. W. Pitera, and K. A. Dill, Multiscale Model. Simul. MMSUBT 1540-3459 5, 1214 (2006). 10.1137/06065146X (Pubitemid 350285620)
17. S. P. Elmer, S. Park, and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 123, 114902 (2005). 10.1063/1.2001648 (Pubitemid 41377093)
18. H. Grubmüller and P. Tavan, J. Chem. Phys. JCPSA6 0021-9606 101, 5047 (1994). 10.1063/1.467427
19. F. Noe and S. Fisher, Curr. Opin. Struct. Biol. COSBEF 0959-440X 18, 154 (2008). 10.1016/j.sbi.2008.01.008
20. G. Jayachandran, V. Vishal, and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 124, 164902 (2006). 10.1063/1.2186317
21. M. Shirts and V. S. Pande, Science SCIEAS 0036-8075 290, 1903 (2000). 10.1126/science.290.5498.1903
22. F. Allen, G. Almasi, W. Andreoni, D. Beece, B. J. Berne, A. Bright, J. Brunheroto, C. Cascaval, J. Castanos, P. Coteus, P. Crumley, A. Curioni, M. Denneau, W. Donath, M. Eleftheriou, B. Fitch, B. Fleischer, C. J. Georgiou, R. Germain, M. Giampapa, D. Gresh, M. Gupta, R. Haring, H. Ho, P. Hochschild, S. Hummel, T. Jonas, D. Lieber, G. Martyna, K. Maturu, J. Moreira, D. Newns, M. Newton, R. Philhower, T. Picunko, J. Pitera, M. Pitman, R. Rand, A. Royyuru, V. Salapura, A. Sanomiya, R. Shah, Y. Sham, S. Singh, M. Snir, F. Suits, R. Swetz, W. C. Swope, N. Vishnumurthy, T. J. C. Ward, H. Warren, and R. Zhou, IBM Syst. J. IBMSA7 0018-8670 40, 310 (2001).
23. P. Deuflhard, W. Huisinga, A. Fischer, and C. Schütte, Linear Algebr. Appl. LAAPAW 0024-3795 315, 39 (2000). 10.1016/S0024-3795(00)00095-1
24. S. Sriraman, I. G. Kevrekidis, and G. Hummer, J. Phys. Chem. B JPCBFK 1089-5647 109, 6479 (2005). 10.1021/jp046448u (Pubitemid 40538414)
25. V. Schultheis, T. Hirschberger, H. Carstens, and P. Tavan, J. Chem. Theory Comput. JCTCCE 1549-9618 1, 515 (2005). 10.1021/ct050020x
26. S. Kumar, D. Bouzida, R. H. Swendsen, P. A. Kollman, and J. H. Rosenberg, J. Comput. Chem. JCCHDD 0192-8651 13, 1011 (1992). 10.1002/jcc.540130812
27. W. C. Swope, J. W. Pitera, and F. Suits, J. Phys. Chem. B JPCBFK 1089-5647 108, 6571 (2004). 10.1021/jp037421y
28. O. M. Becker, Proteins 18PLAF 0887-3585 27, 213 (1997). 10.1002/(SICI)1097-0134(1997 02)27:2<213::AID-PROT8>3.0.CO;2-G
29. M. Andrec, A. K. Felts, E. Gallicchio, and R. M. Levy, Proc. Natl. Acad. Sci. U.S.A. PNASA6 0027-8424 102, 6801 (2005). 10.1073/pnas.0408970102 (Pubitemid 40675396)
30. B. L. de Groot, X. Daura, A. E. Mark, and H. Grubmuller, J. Mol. Biol. JMOBAK 0022-2836 309, 299 (2001). 10.1006/jmbi.2001.4655 (Pubitemid 32553525)
31. M. Weber and T. Galliat, Konrad-Zuse-Zentrum für Informationstechnik in Berlin (ZIB) Technical Report No. 02-12, 2002.
32. P. Deuflhard and M. Weber, Linear Algebr. Appl. LAAPAW 0024-3795 398, 161 (2005). 10.1016/j.laa.2004.10.026 (Pubitemid 40220042)
33. M. Weber, Konrad-Zuse-Zentrum für Informationstechnik in Berlin (ZIB) Technical Report No. 03-04, 2003.
34. F. Nó, I. Horenko, C. Schütte, and J. C. Smith, J. Chem. Phys. JCPSA6 0021-9606 126, 155102 (2007). 10.1063/1.2714539 (Pubitemid 46644961)
35. S. Kube and P. Deuflhard, Konrad-Zuse-Zentrum für Informationstechnik in Berlin (ZIB) Technical Report No. 03-19, 2006.
36. S. Bacallado, J. D. Chodera, and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 131, 045106 (2009). 10.1063/1.3192309
37. W. C. Swope, J. W. Pitera, F. Suits, M. Pitman, M. Eleftheriou, B. G. Fitch, R. S. Germain, A. Rayshubski, T. J. C. Ward, Y. Zhestkov, and R. Zhou, J. Phys. Chem. B JPCBFK 1089-5647 108, 6582 (2004). 10.1021/jp037422q
38. N. Singhal and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 123, 204909 (2005). 10.1063/1.2116947 (Pubitemid 41721321)
39. N. S. Hinrichs and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 126, 244101 (2007). 10.1063/1.2740261 (Pubitemid 47141397)
40. N. -V. Buchete and G. Hummer, J. Phys. Chem. B JPCBFK 1089-5647 112, 6057 (2008). 10.1021/jp0761665
41. P. Diaconis and S. Rolles, Ann. Stat. ASTSC7 0090-5364 34, 1270 (2006). 10.1214/009053606000000290 (Pubitemid 44231163)
42. S. Park and V. S. Pande, J. Chem. Phys. JCPSA6 0021-9606 124, 054118 (2006).
43. In the case of dynamical models that are defined as random walks on a microscopic time-scale, a mesoscopic Markov model should not be confused with the random walk microscopic dynamics. Even for such a dynamical system, we can ask whether a mesoscopic Markov model exists for any particular value of τmes. Moreover, Markov behavior on the mesoscopic time-scale for mesostates is not a trivial consequence of random walk microscopic dynamics.
44. We do not assume that the Markov model is reversible (i.e., that it satisfies the detailed balance condition with regard to some equilibrium distribution). We have not been able to construct a suitable conjugate prior probability for reversible Markov chains that facilitates the Bayesian calculation as much as does the Dirichlet prior mentioned just below. In practice, for the test cases in the next section, the transition probability matrix for the "best CMMM," i.e., the one that maximizes the posterior probability (see below) satisfies the detailed balance condition except for statistical error related to the fact that a finite number of transitions occur in the original simulation data.
45. See supplementary material at http://dx.doi.org/10.1063/1.3496438 E-JCPSA6-133-015038 for a more detailed description of the partitioning algorithm.
46. 2) 1/2, where Δφ and Δψ are the differences between the two angles in the two configurations.
47. P. G. Bolhuis, C. Dellago, and D. Chandler, Proc. Natl. Acad. Sci. U.S.A. PNASA6 0027-8424 97, 5877 (2000). 10.1073/pnas.100127697 (Pubitemid 30367492)
48. P. Smith, J. Chem. Phys. JCPSA6 0021-9606 111, 5568 (1999). 10.1063/1.479860