Bistability in the chemical master equation for dual phosphorylation cycles

Abstract

Dual phospho/dephosphorylation cycles, as well as covalent enzymatic-catalyzed modifications of substrates are widely diffused within cellular systems and are crucial for the control of complex responses such as learning, memory and cellular fate determination. Despite the large body of deterministic studies and the increasing work aimed at elucidating the effect of noise in such systems, some aspects remain unclear. Here we study the stationary distribution provided by the two-dimensional Chemical Master Equation for a well-known model of a two step phospho/dephosphorylation cycle using the quasi-steady state approximation of enzymatic kinetics. Our aim is to analyze the role of fluctuations and the molecules distribution properties in the transition to a bistable regime. When detailed balance conditions are satisfied it is possible to compute equilibrium distributions in a closed and explicit form. When detailed balance is not satisfied, the stationary non-equilibrium state is strongly influenced by the chemical fluxes. In the last case, we show how the external field derived from the generation and recombination transition rates, can be decomposed by the Helmholtz theorem, into a conservative and a rotational (irreversible) part. Moreover, this decomposition, allows to compute the stationary distribution via a perturbative approach. For a finite number of molecules there exists diffusion dynamics in a macroscopic region of the state space where a relevant transition rate between the two critical points is observed. Further, the stationary distribution function can be approximated by the solution of a Fokker-Planck equation. We illustrate the theoretical results using several numerical simulations.

Irreversible thermodynamics of open chemical networks. I. Emergent cycles and broken conservation laws

“…Finally, and most importantly, the results presented in this paper are meant to be applicable to the stochastic thermodynamics associated with the chemical master equation, 10,[27][28][29] where one typically deals with force constraints (encoded in the rates) rather than with current constraints. Chemostats play a central role in the derivation of advanced results in the stochastic theory of chemical reactions, such as the Fluctuation Theorem for the currents and the Green-Kubo relations.…”

Modeling and inference of genetic interactions

“…The experimental data observation technology, the sampling rate of data extraction, and the purpose of modeling significantly affect the choice of the genetic regulatory network model. We have observed that there is an inherent stochasticity in the generation of RNA and protein molecules, and further developments in observation of singular cell RNA and protein expressions will provide an impetus to detailed stochastic modeling . Some instances of cellular behavior that are hard to capture using the commonly used ODE models are alteration of dynamic behavior due to intrinsic noise, bistable behavior in gene expression, effect of delay on oscillations, or noise‐induced signal amplification and noise‐driven divergence of cell fates .…”

Discrete flux and velocity fields of probability and their global maps in reaction systems

“…We then examined the steady state behavior of the system at two conditions of the copy numbers of species A and B: (a, b) = (10, 50) and (a, b) = (20,40). The probability landscape in − log(p(x, t)) for (a, b) = (10, 50) shown in Fig.…”