Estimating first‐passage time distributions from weighted ensemble simulations and non‐Markovian analyses. (9th September 2015)

• Record Type: Journal Article
• Authors: Suárez, Ernesto; Pratt, Adam J.; Chong, Lillian T.; Zuckerman, Daniel M.
• Is Part Of: Protein science. Volume 25:Number 1(2016:Jan.)
• Page Start: 67
• Page End: 78
• Publication Date: 2015-09-09
• Subjects: weighted ensemble -- rare event -- first‐passage time -- non‐Markovian; Proteins -- Periodicals; 572.6
• DOI: 10.1002/pro.2738
• Languages: English
• ISSNs: 0961-8368

Abstract

Abstract: First‐passage times (FPTs) are widely used to characterize stochastic processes such as chemical reactions, protein folding, diffusion processes or triggering a stock option. In previous work (Suarez et al ., JCTC 2014;10:2658‐2667), we demonstrated a non‐Markovian analysis approach that, with a sufficient subset of history information, yields unbiased mean first‐passage times from weighted‐ensemble (WE) simulations. The estimation of the distribution of the first‐passage times is, however, a more ambitious goal since it cannot be obtained by direct observation in WE trajectories. Likewise, a large number of events would be required to make a good estimation of the distribution from a regular "brute force" simulation. Here, we show how the previously developed non‐Markovian analysis can generate approximate, but highly accurate, FPT distributions from WE data. The analysis can also be applied to any other unbiased trajectories, such as from standard molecular dynamics simulations. The present study employs a range of systems with independent verification of the distributions to demonstrate the success and limitations of the approach. By comparison to a standard Markov analysis, the non‐Markovian approach is less sensitive to the user‐defined discretization of configuration space.