Abstract:Significance Modeling intracellular processes has long relied on the markovian assumption. However, as soon as a reactant interacts with its environment, molecular memory definitely exists and its effects cannot be neglected. Since the Markov theory cannot translate directly to modeling and analysis of nonmarkovian processes, this leads to many significant challenges. We develop a formulation, namely the stationary generalized chemical-master equation, to model intracellular processes with molecular memory. This formulation converts a nonmarkovian question to a markovian one while keeping the stationary probabilistic behavior unchanged. Both a stationary generalized Fokker–Planck equation and a generalized linear noise approximation are further developed for the fast evaluation of fluctuations. These formulations can have broad applications and may help us discover new biological knowledge. Many cellular processes are governed by stochastic reaction events. These events do not necessarily occur in single steps of individual molecules, and, conversely, each birth or death of a macromolecule (e.g., protein) could involve several small reaction steps, creating a memory between individual events and thus leading to nonmarkovian reaction kinetics. Characterizing this kinetics is challenging. Here, we develop a systematic approach for a general reaction network with arbitrary intrinsic waiting-time distributions, which includes the stationary generalized chemical-master equation (sgCME), the stationary generalized Fokker–Planck equation, and the generalized linear-noise approximation. The first formulation converts a nonmarkovian issue into a markovian one by introducing effective transition rates (that explicitly decode the effect of molecular memory) for the reactions in an equivalent reaction network with the same substrates but without molecular memory. Nonmarkovian features of the reaction kinetics can be revealed by solving the sgCME. The latter 2 formulations can be used in the fast evaluation of fluctuations. These formulations can have broad applications, and, in particular, they may help us discover new biological knowledge underlying memory effects. When they are applied to generalized stochastic models of gene-expression regulation, we find that molecular memory is in effect equivalent to a feedback and can induce bimodality, fine-tune the expression noise, and induce switch.

Computation of stationary distributions in stochastic models of cellular processes with molecular memory

Markovian Approaches to Modeling Intracellular Reaction Processes with Molecular Memory

Analysis of singularly perturbed stochastic chemical reaction networks motivated by applications to epigenetic cell memory

Computationally Efficient Modelling of Stochastic Spatio-Temporal Dynamics in Biomolecular Networks

Analytical Results for Non-Markovian Models of Bursty Gene Expression

Effective computational methods for hybrid stochastic gene networks

Analytical Results for a Generalized Model of Bursty Gene Expression with Molecular Memory

Inference of Spatiotemporal Effects on Cellular State Transitions from Time-Lapse Microscopy

Accurate estimates of dynamical statistics using memory

A multiscale compartment-based model of stochastic gene regulatory networks using hitting-time analysis

Stochastic reaction networks within interacting compartments

Optimal Switching Paths of a Non-Markovian Model of Gene Bursty Expression

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Modelling non-Markovian dynamics in biochemical reactions

The Roles of Space and Stochasticity in Computational Simulations of Cellular Biochemistry: Quantitative Analysis and Qualitative Insights

Molecular computing for Markov chains

Cellular automaton models for time-correlated random walks: derivation and analysis

Exact Results for Queuing Models of Stochastic Transcription with Memory and Crosstalk

Hybrid Master Equation for Jump-Diffusion Approximation of Biomolecular Reaction Networks