Gaussian Markov transition models of molecular kinetics.
The slow processes of molecular dynamics (MD) simulations--governed by dominant eigenvalues and eigenfunctions of MD propagators--contain essential information on structures of and transition rates between long-lived conformations. Existing approaches to this problem, including Markov state models and the variational approach, represent the dominant eigenfunctions as linear combinations of a set of basis functions. However the choice of the basis functions and their systematic statistical estimation are unsolved problems. Here, we propose a new class of kinetic models called Markov transition models (MTMs) that approximate the transition density of the MD propagator by a mixture of probability densities. Specifically, we use Gaussian MTMs where a Gaussian mixture model is used to approximate the symmetrized transition density. This approach allows for a direct computation of spectral components. In contrast with the other Galerkin-type approximations, our approach can automatically adjust the involved Gaussian basis functions and handle the statistical uncertainties in a Bayesian framework. We demonstrate by some simulation examples the effectiveness and accuracy of the proposed approach.
Department for Mathematics and Computer Science, FU Berlin, Arnimallee 6, 14195 Berlin, Germany.
Hao
H
Wu
Department for Mathematics and Computer Science, FU Berlin, Arnimallee 6, 14195 Berlin, Germany.
Frank
F
Noé
D000409:Alanine / Q000737:chemistry; D000465:Algorithms; D004058:Diffusion; D007700:Kinetics; D016013:Likelihood Functions; D008390:Markov Chains; D056004:Molecular Dynamics Simulation; D010455:Peptides / Q000737:chemistry
D010455:Peptides; D000409:Alanine