Using stochastic models calibrated from nanosecond nonequilibrium simulations to approximate mesoscale information.
J Chem Phys
; 130(14): 144908, 2009 Apr 14.
Article
en En
| MEDLINE
| ID: mdl-19368472
ABSTRACT
We demonstrate how the surrogate process approximation (SPA) method can be used to compute both the potential of mean force along a reaction coordinate and the associated diffusion coefficient using a relatively small number (10-20) of bidirectional nonequilibrium trajectories coming from a complex system. Our method provides confidence bands which take the variability of the initial configuration of the high-dimensional system, continuous nature of the work paths, and thermal fluctuations into account. Maximum-likelihood-type methods are used to estimate a stochastic differential equation (SDE) approximating the dynamics. For each observed time series, we estimate a new SDE resulting in a collection of SPA models. The physical significance of the collection of SPA models is discussed and methods for exploiting information in the population of estimated SPA models are demonstrated and suggested. Molecular dynamics simulations of potassium ion dynamics inside a gramicidin A channel are used to demonstrate the methodology, although SPA-type modeling has also proven useful in analyzing single-molecule experimental time series [J. Phys. Chem. B 113, 118 (2009)].
Texto completo:
1
Bases de datos:
MEDLINE
Asunto principal:
Potasio
/
Biología Computacional
/
Difusión
/
Gramicidina
Idioma:
En
Revista:
J Chem Phys
Año:
2009
Tipo del documento:
Article
País de afiliación:
Estados Unidos