Universally valid reduction of multiscale stochastic biochemical systems using simple non-elementary propensities

Abstract

Biochemical systems consist of numerous elementary reactions governed by the law of mass action. However, experimentally characterizing all the elementary reactions is nearly impossible. Thus, over a century, their deterministic models that typically contain rapid reversible bindings have been simplified with non-elementary reaction functions (e.g., Michaelis-Menten and Morrison equations). Although the non-elementary reaction functions are derived by applying the quasi-steady-state approximation (QSSA) to deterministic systems, they have also been widely used to derive propensities for stochastic simulations due to computational efficiency and simplicity. However, the validity condition for this heuristic approach has not been identified even for the reversible binding between molecules, such as protein-DNA, enzyme-substrate, and receptor-ligand, which is the basis for living cells. Here, we find that the non-elementary propensities based on the deterministic total QSSA can accurately capture the stochastic dynamics of the reversible binding in general. However, serious errors occur when reactant molecules with similar levels tightly bind, unlike deterministic systems. In that case, the non-elementary propensities distort the stochastic dynamics of a bistable switch in the cell cycle and an oscillator in the circadian clock. Accordingly, we derive alternative non-elementary propensities with the stochastic low-state QSSA, developed in this study. This provides a universally valid framework for simplifying multiscale stochastic biochemical systems with rapid reversible bindings, critical for efficient stochastic simulations of cell signaling and gene regulation. To facilitate the framework, we provide a user-friendly open-source computational package, ASSISTER, that automatically performs the present framework.

Combined multiple transcriptional repression mechanisms generate ultrasensitivity and oscillations

“…Specifically, the concentration of transcriptional factors ( A T ) is 2×10−9–10−7 M as their number is 104–105 (i.e. 2×10−20–10−19 mol) and the typical mammalian cell volume is 10−11–10−12 l [32,38,39]. Thus, even the extremely high affinity protein whose dissociation constant is picomolar (i.e.…”

Systematic inference identifies a major source of heterogeneity in cell signaling dynamics: The rate-limiting step number

“…To analyze such complex dynamics with MBI, we extended the delayed birthdeath process by incorporating feedback inhibition with a fixed delay ( 2 ) using the sequestration-based repression function (Fig. 6A, ii) (49)(50)(51)(52). Then, we derived its time-varying exact low-order moments by using a transient Little's law (39) and the Chemical Fluctuation Theorem (see Supplementary Text F) (21).…”

Spatially coordinated collective phosphorylation filters spatiotemporal noises for precise circadian timekeeping

“…This allowed us to describe the effects of various noise sources on circadian rhythms. To use the Gillespie algorithm for a system containing non-elementary propensity functions (e.g., Hill functions) other than functions from mass-action kinetics, timescale separation is necessary [67][68][69] . Thus, we assumed that the synthesis and the degradation of PER are slower than other reactions, such as binding and unbinding of PER to activator proteins and each step of its phosphorylation and dephosphorylation.…”