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.

Molecular Synthesis for Probability Theory and Stochastic Process

Synthesis of Probability Theory Based on Molecular Computation

Synthesizing Markov Chain With Reversible Unimolecular Reactions

Molecular computing for Markov chains

Markovian Approaches to Modeling Intracellular Reaction Processes with Molecular Memory

DNA-based Chemical Reaction Networks.

DNA‐Based Chemical Reaction Networks

Stochastic Modelling of Biochemical Systems of Multi-Step Reactions Using A Simplified Two-Variable Model

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

A Moment-Convergence Method for Stochastic Analysis of Biochemical Reaction Networks

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Synthesis of Stochastic Flow Networks

Stochastic reaction networks within interacting compartments

Stochastic Chemical Reaction Networks with Discontinuous Limits and AIMD processes

Robust control of biochemical reaction networks via stochastic morphing

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

A Stochastic Framework for the Design of Transient and Steady State Behavior of Biochemical Reaction Networks

Throwing and Manipulating and Cheating with a DNA Nano-Dice

Synthetic likelihood method for reaction network inference

Effective computational methods for hybrid stochastic gene networks