A Bayesian method for construction of Markov models to describe dynamics on various time scales [electronic resource]
- Emily Kathleen Rains.
- 2010, c2011.
- Physical description
- 1 online resource.
Also available at
At the library
All items must be viewed on site
Request items at least 2 days before you visit to allow retrieval from off-site storage. You can request at most 5 items per day.
|3781 2010 R||In-library use|
- The dynamics of many biological processes of interest, such as the folding of a protein, are slow and complicated enough that a single molecular dynamics simulation trajectory of the entire process is difficult to obtain in any reasonable amount of time. Moreover, one such simulation may not be sufficient to develop an understanding of the mechanism of the process, and multiple simulations may be necessary. One approach to circumvent this computational barrier is the use of Markov state models. These models are useful because they can be constructed using data from a large number of shorter simulations instead of a single long simulation. This thesis presents a new Bayesian method for the construction of Markov models from simulation data. A Markov model is specified by (t, P, T), where t is the mesoscopic time step, P is a partition of configuration space into mesostates, and T is an N x N transition rate matrix for transitions between the mesostates in one mesoscopic time step, where N is the number of mesostates in P. The method presented here is different from previous Bayesian methods in several ways. 1. The method uses Bayesian analysis to determine the partition as well as the transition probabilities. 2. The method allows the construction of a Markov model for any chosen mesoscopic time scale t. 3. It constructs Markov models for which the diagonal elements of T are all equal to or greater than 0.5. Such a model will be called a 'consistent mesoscopic Markov model' (or CMMM). Such models have important advantages for providing an understanding of the dynamics on a mesoscopic time scale. The Bayesian method uses simulation data to find a posterior probability distribution for (P, T) for any chosen t. This distribution can be regarded as the Bayesian probability that the kinetics observed in the atomistic simulation data on the mesoscopic time scale t was generated by the CMMM specified by (P, T). An optimization algorithm is used to find the most probable CMMM for the chosen mesoscopic time step. We applied this method of Markov model construction to several toy systems (random walks in one and two dimensions) as well as the dynamics of alanine dipeptide in water and of trpzip2 in water. The resulting Markov state models were indeed successful in capturing the dynamics of our test systems on a variety of mesoscopic time scales.
- Publication date
- Copyright date
- Submitted to the Department of Chemistry.
- Ph.D. Stanford University 2011
Browse related items
Start at call number: