Reduced Stochastic Models for Complex Molecular Systems

We present a new numerical method for the identification of the most important metastable states of a system with complicated dynamical behavior from time series information. The approach is based on the representation of the effective dynamics of the full system by a Markov jump process between metastable states, and the dynamics within each of these metastable states by rather simple stochastic differential equations (SDEs). Its algorithmic realization exploits the concept of hidden Markov models (HMMs) with output behavior given by SDEs. A first complete algorithm including an explicit Euler–Maruyama-based likelihood estimator has already been presented in Horenko et al. (MMS, 2006a). Herein, we present a semi-implicit exponential estimator that, in contrast to the Euler–Maruyama-based estimator, also allows for reliable parameter optimization for time series where the time steps between single observations are large. The performance of the resulting method is demonstrated for some generic examples, in detail compared to the Euler–Maruyama-based estimator, and finally applied to time series originating from a 100 ns B-DNA molecular dynamics simulation.

[1]  P. Deuflharda,et al.  Identification of almost invariant aggregates in reversible nearly uncoupled Markov chains , 2000 .

[2]  H Frauenfelder,et al.  Spectroscopic evidence for conformational relaxation in myoglobin. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Van Nostrand,et al.  Error Bounds for Convolutional Codes and an Asymptotically Optimum Decoding Algorithm , 1967 .

[4]  L. Baum,et al.  An inequality and associated maximization technique in statistical estimation of probabilistic functions of a Markov process , 1972 .

[5]  P. Deuflhard,et al.  Robust Perron cluster analysis in conformation dynamics , 2005 .

[6]  O. Junge,et al.  On the Approximation of Complicated Dynamical Behavior , 1999 .

[7]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[8]  C. Lubich,et al.  A Poisson Integrator for Gaussian Wavepacket Dynamics , 2006 .

[9]  Jeff A. Bilmes,et al.  A gentle tutorial of the em algorithm and its application to parameter estimation for Gaussian mixture and hidden Markov models , 1998 .

[10]  Alexander Fischer,et al.  Hybrid Monte Carlo with adaptive temperature in mixed-canonical ensemble: Efficient conformational analysis of RNA , 1998, J. Comput. Chem..

[11]  C. Schütte,et al.  From Simulation Data to Conformational Ensembles , 1999 .

[12]  Alexander Fischer,et al.  Identification of biomolecular conformations from incomplete torsion angle observations by hidden markov models , 2007, J. Comput. Chem..

[13]  Alexander Fischer,et al.  Automated Model Reduction for Complex Systems Exhibiting Metastability , 2006, Multiscale Model. Simul..

[14]  Hybrid Monte Carlo with adaptive temperature in mixed‐canonical ensemble: Efficient conformational analysis of RNA , 1998 .

[15]  A. Bandrivskyy,et al.  Model reconstruction of nonlinear dynamical systems driven by noise , 2003 .

[16]  Wilhelm Huisinga,et al.  From simulation data to conformational ensembles: Structure and dynamics‐based methods , 1999 .

[17]  M. Karplus,et al.  Multiple conformational states of proteins: a molecular dynamics analysis of myoglobin. , 1987, Science.

[18]  P. J. Steinbach,et al.  Conformational relaxation in heme proteins: Ligand rebinding above the glass transition , 1990 .

[19]  Michael Dellnitz,et al.  Congestion and Almost Invariant Sets in Dynamical Systems , 2001, SNSC.

[20]  P. Nguyen,et al.  Energy landscape of a small peptide revealed by dihedral angle principal component analysis , 2004, Proteins.