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.

Synthesizing Markov Chain With Reversible Unimolecular Reactions

Molecular computing for Markov chains

Molecular Synthesis for Probability Theory and Stochastic Process

Implementation Of Mealy Machine With Molecular Reactions

Markovian Approaches to Modeling Intracellular Reaction Processes with Molecular Memory

Synthesis of Probability Theory Based on Molecular Computation

ReacNetGen: an Automatic Reaction Network Generator for Reactive Molecular Dynamic Simulations

Markov Chain Monte Carlo in the Analysis of Single-Molecule Experimental Data

Universality of a Surface Chemical Reaction Network Using Only Bi-Molecular Reactions

ReacNetGenerator: an Automatic Reaction Network Generator for Reactive Molecular Dynamic Simulations

Hierarchical Markov State Model Building to Describe Molecular Processes

Programming and training rate-independent chemical reaction networks

Arithmetic Computations Based on Chemical Reaction Networks

Markov field models: Scaling molecular kinetics approaches to large molecular machines

A Cellular Automaton Approach for Efficient Computing on Surface Chemical Reaction Networks

Reviving the Two-state Markov Chain Approach (Technical Report)

DNA-based Chemical Reaction Networks.

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

Markovian Dynamics on Complex Reaction Networks

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

DNA‐Based Chemical Reaction Networks