Refine
Language
- English (21)
Keywords
- molecular dynamics (3)
- Markov models (2)
- Bayesian error analysis (1)
- Brownian dynamics (1)
- FCS (1)
- FRET (1)
- MSM (1)
- Markov (state) models (1)
- Molecular kinetics (1)
- Phosphate Binding Protein (1)
Application Area
- A (21)
Large-scale stochastic models are relevant in many different fields such as com- putational biology, finance, social sciences, communication and traffic networks. In order to both efficiently simulate and analyze such models and to understand the essential properties of the sys- tem, it is desirable to have model reduction techniques that much reduce the dimensionality of the model while at the same time preserving the system’s essential dynamical properties. In this paper, a general model reduction technique for the class of discrete space and time Hidden Markov Models is presented, thereby also including the more special class discrete Markov Chains. The method is illustrated on some model applications.
Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe
long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great
amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into
microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data
to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current
perspective of Markov models and explain in short their theoretical basis and assumptions. We
describe transition path theory which allows the entire ensemble of protein folding pathways to be
investigated and that combines naturally with Markov models. Experimental observations can be
naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally
observable timescales can be equipped with an understanding of the structural rearrangement
processes that take place at these timescales. The concepts of this paper are illustrated by a
simple kinetic model of protein folding.
Discrete-state Markov (or master equation) models provide a useful simplified representation for
characterizing the long-time statistical evolution of biomolecules in a manner that allows direct
comparison with experiments as well as the elucidation of mechanistic pathways for an inherently
stochastic process. A vital part of meaningful comparison with experiment is the characterization of
the statistical uncertainty in the predicted experimental measurement, which may take the form of
an equilibrium measurement of some spectroscopic signal, the time-evolution of this signal following
a perturbation, or the observation of some statistic (such as the correlation function) of the equilib-
rium dynamics of a single molecule. Without meaningful error bars (which arise due to the finite
quantity of data used to construct the model), there is no way to determine whether the deviations
between model and experiment are statistically meaningful. Previous work has demonstrated that
a Bayesian method that enforces microscopic reversibility can be used to characterize the correlated
uncertainties in state-to-state transition probabilities (and functions thereof) for a model inferred from
molecular simulation data. Here, we extend this approach to include the uncertainty in observables
that are functions of molecular conformation (such as surrogate spectroscopic signals) characteriz-
ing each state, permitting the full statistical uncertainty in computed spectroscopic experiments to be
assessed. We test the approach in a simple model system to demonstrate that the computed uncer-
tainties provide a useful indictor of statistical variation, and then apply it to the computation of the
fluorescence autocorrelation function measured for a dye-labeled peptide previously studied by both
experiment and simulation.
Optimal Identification of Semi-Rigid Domains in Macromolecules from Molecular Dynamics Simulation
(2012)
Biological function relies on the fact that biomolecules can switch between different conformations and aggregation states.
Such transitions involve a rearrangement of parts of the biomolecules involved that act as dynamic domains. The reliable
identification of such domains is thus a key problem in biophysics. In this work we present a method to identify semi-rigid
domains based on dynamical data that can be obtained from molecular dynamics simulations or experiments. To this end
the average inter-atomic distance-deviations are computed. The resulting matrix is then clustered by a constrained
quadratic optimization problem. The reliability and performance of the method are demonstrated for two artificial peptides.
Furthermore we correlate the mechanical properties with biological malfunction in three variants of amyloidogenic
transthyretin protein, where the method reveals that a pathological mutation destabilizes the natural dimer structure of the
protein. Finally the method is used to identify functional domains of the GroEL-GroES chaperone, thus illustrating the
efficiency of the method for large biomolecular machines.
In this paper, we present a Gaussian Markov random field (GMRF) model for the transition
matrices (TMs) of Markov chains (MCs) by assuming the existence of a neighborhood relationship
between states, and develop the maximum a posteriori (MAP) estimators under different obser-
vation conditions. Unlike earlier work on TM estimation, our method can make full use of the
similarity between different states to improve the estimated accuracy, and the estimator can be
performed very efficiently by solving a convex programming problem. In addition, we discuss the
parameter choice of the proposed model, and introduce a Monte Carlo cross validation (MCCV)
method. The numerical simulations of a diffusion process are employed to show the effectiveness
of the proposed models and algorithms.
Hybrid systems are often used to describe many complex dynamic phenomena by combining multiple modes of
dynamics into whole systems. In this paper, we present a flat Dirichlet process switching (FDPS) model that defines
a prior on mode switching dynamics of hybrid systems. Compared with the classical Markovian jump system (MJS)
models, the FDPS model is nonparametric and can be applied to the hybrid systems with an unbounded number of
potential modes. On the other hand, the probability structure of the new model is simpler and more flexible than the
recently proposed hierarchical Dirichlet process (HDP) based MJS. Furthermore, we develop a Markov chain Monte
Carlo (MCMC) method for estimating the states of hybrid systems with FDPS prior. And the numerical simulations
of a hybrid system in different conditions are employed to show the effectiveness of the proposed approach.
Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility.
Here, we present a model of the actin filament (F-actin) that incorporates the global structure
of the recently published model by Oda et al. but also conserves internal stereochemistry. A
comparison is made using molecular dynamics simulation of the model with other recent F-
actin models. A number of structural determents such as the protomer propeller angle, the
number of hydrogen bonds and the structural variation among the protomers are analyzed.
The MD comparison is found to reflect the evolution in quality of actin models over the last
six years. In addition, simulations of the model are carried out in states with both ADP or
ATP bound and local hydrogen-bonding differences characterized. The results point to the
significance of a direct interaction of Gln137 with ATP for activation of ATPase activity after
the G-to-F-actin transition.
Single-molecule force spectroscopy has proven to be a powerful tool for studying the kinetic be- havior of biomolecules. Through application of an external force, conformational states with small or transient populations can be stabilized, allowing them to be characterized and the statistics of in- dividual trajectories studied to provide insight into biomolecular folding and function. Because the observed quantity (force or extension) is not necessarily an ideal reaction coordinate, individual ob- servations cannot be uniquely associated with kinetically distinct conformations. While maximum- likelihood schemes such as hidden Markov models have solved this problem for other classes of single-molecule experiments by using temporal information to aid in the inference of a sequence of distinct conformational states, these methods do not give a clear picture of how precisely the model parameters are determined by the data due to instrument noise and finite-sample statistics, both sig- nificant problems in force spectroscopy. We solve this problem through a Bayesian extension that allows the experimental uncertainties to be directly quantified, and build in detailed balance to fur- ther reduce uncertainty through physical constraints. We illustrate the utility of this approach in characterizing the three-state kinetic behavior of an RNA hairpin in a stationary optical trap.
While seemingly straightforward in principle, the reliable estimation of rate constants is seldom easy in practice. Numerous issues, such as the complication of poor reaction coordinates, cause obvious approaches to yield unreliable estimates. When a reliable order parameter is available, the reactive flux theory of Chandler allows the rate constant to be extracted from the plateau region of an appropriate reactive flux correlation function. However, when applied to real data from single- molecule experiments or molecular dynamics simulations, the reactive flux correlation function requires the numerical differentiation of a noisy empirical correlation function, which can result in an unacceptably poor estimate of the rate and pathological dependence on the sampling interval. We present a modified version of this theory which does not require numerical derivatives, allowing rate constants to be robustly estimated from the time-correlation function directly. We illustrate the approach using single-molecule passive force spectroscopy measurements of an RNA hairpin.
Folding and conformational changes of macromolecules often require the generation of large amounts of simulation data that are difficult to ana- lyze. Markov state models (MSMs) address this challenge by providing a systematic way to decompose the state space of the molecular system into substates and to estimate a transition matrix containing the transi- tion probabilities between these substates. This transition matrix can be analyzed to reveal the metastable, i.e. long-living, states of the system, its slowest relaxation timescales and transition pathways and rates e.g. from unfolded to folded, or from dissociated to bound states. To reduce the technical burden of constructing such MSMs we provide the software framework EMMA (available at https://simtk.org/home/emma) to con- struct, validate and analyse such Markov State Models.