Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Comparison of the kinetics of different Markov models for ligand binding under varying conditions

Journal Article · · Journal of Chemical Physics
DOI:https://doi.org/10.1063/1.4908531· OSTI ID:22416202
 [1];  [2]
  1. Max Planck Institute for Developmental Biology, Tübingen (Germany)
  2. Felix Bernstein Institute for Mathematical Statistics in the Biosciences, University of Göttingen, Göttingen (Germany)

We recently derived a Markov model for macromolecular ligand binding dynamics from few physical assumptions and showed that its stationary distribution is the grand canonical ensemble [J. W. R. Martini, M. Habeck, and M. Schlather, J. Math. Chem. 52, 665 (2014)]. The transition probabilities of the proposed Markov process define a particular Glauber dynamics and have some similarity to the Metropolis-Hastings algorithm. Here, we illustrate that this model is the stochastic analog of (pseudo) rate equations and the corresponding system of differential equations. Moreover, it can be viewed as a limiting case of general stochastic simulations of chemical kinetics. Thus, the model links stochastic and deterministic approaches as well as kinetics and equilibrium described by the grand canonical ensemble. We demonstrate that the family of transition matrices of our model, parameterized by temperature and ligand activity, generates ligand binding kinetics that respond to changes in these parameters in a qualitatively similar way as experimentally observed kinetics. In contrast, neither the Metropolis-Hastings algorithm nor the Glauber heat bath reflects changes in the external conditions correctly. Both converge rapidly to the stationary distribution, which is advantageous when the major interest is in the equilibrium state, but fail to describe the kinetics of ligand binding realistically. To simulate cellular processes that involve the reversible stochastic binding of multiple factors, our pseudo rate equation model should therefore be preferred to the Metropolis-Hastings algorithm and the Glauber heat bath, if the stationary distribution is not of only interest.

OSTI ID:
22416202
Journal Information:
Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 9 Vol. 142; ISSN JCPSA6; ISSN 0021-9606
Country of Publication:
United States
Language:
English

Similar Records

A parallel evolutionary multiple-try metropolis Markov chain Monte Carlo algorithm for sampling spatial partitions
Journal Article · Mon Jan 11 23:00:00 EST 2021 · Statistics and Computing · OSTI ID:1777811

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo
Journal Article · Mon Jun 01 00:00:00 EDT 2020 · International Journal for Uncertainty Quantification · OSTI ID:1670752

Accelerating Markov Chain Monte Carlo sampling with diffusion models
Journal Article · Tue Dec 12 23:00:00 EST 2023 · Computer Physics Communications · OSTI ID:2281802