# American Institute of Mathematical Sciences

## Analysis of non-Markovian effects in generalized birth-death models

 Key Laboratory of Computational Mathematics, Guangdong Province and School of, Mathematics, Sun Yat-sen University, Guangzhou, 510275, China

* Corresponding author: Tianshou Zhou

Received  October 2019 Revised  June 2020 Published  August 2020

Birth-death processes are a fundamental reaction module for which we can find its prototypes in many scientific fields. For such a kind of module, if all the reaction events are Markovian, the reaction kinetics is simple. However, experimentally observable quantities are in general consequences of a series of reactions, implying that the synthesis of a macromolecule in general involve multiple middle reaction steps with some reactions that would not be specified by experiments. This multistep process can create molecular memory between reaction events, leading to non-Markovian behavior. Based on the theoretical framework established in a recent paper published in [39], we find that the effect of non-Markovianity is equivalent to the introduction of a feedback, non-Markovianity always amplifies the mean level of the product if the death reaction is non-Markovian but always reduces the mean level if the birth reaction is non-Markovian, and in contrast to Markovianity, non-Markovianity can reduce or amplify the product noise, depending on the details of waiting-time distributions characterizing reaction events. Examples analysis indicates that non-Markovianity, whose effects were neglected in previous studies, can significantly impact gene expression.

Citation: Zhenquan Zhang, Meiling Chen, Jiajun Zhang, Tianshou Zhou. Analysis of non-Markovian effects in generalized birth-death models. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2020254
##### References:

show all references

##### References:
Schematic diagram for birth-death reaction process (A), where $\psi_{1} \left(t\right)$ and $\psi_{2} \left(t;n\right)$ are interreaction waiting-time distributions, which may be exponential (B) or non-exponential (C). Here $n$ represents the number of molecules of reactive species $X$
Difference between dynamic probability distributions obtained in the original non-Markovian model (curves with empty circles) and in the constructed Markovian model (colored curves), both being used in the modeling of a gene transcription process. Parameter values are set as $k_{1} = 3,\lambda_{1} = 90,\lambda_{2} = 1$, and the initial conditions are set as the same in two cases. In the diagram, '$t = \inf$' means that the distributions after $t>10$ are approximately the same, implying that the stationary distribution exists
Characteristics of stationary protein distribution in a generalized model of constitutive gene expression (A), where solid lines represent theoretical predictions and empty circles represent numerical results obtained by the Gillespie algorithm [11]. Parameter values are set as: (B) $k_{d} = 1,\lambda_{d} = 1$ and (C) $k_{b} = 1,\lambda_{b} = 10$. In (B) and (C), we use 1000 realizations to obtain numerical results and error bars to indicate the error ranges of the numerical results
Characteristics of stationary protein distribution in a generalized model of bursty gene expression (A), where solid lines represent theoretical predictions and empty circles represent numerical results obtained by the Gillespie algorithm [11]. Parameter values are set as: (B) $\mu = 10,\delta = 1$, $b = 2,k_{on} = 1,\lambda_{on} = 0.5,k_{off} = 1,\lambda_{off} = 0.1$ for Markovianity, whereas $k_{on} = 4,\lambda_{on} = 4$ and the other parameter values are kept unchanged for non-Markovianity; (C)$\mu = 10,\delta = 1,b = 2$, $k_{on} = 1,\lambda_{on} = 0.5$, $k_{off} = 1,\lambda_{off} = 0.1$ for Markovianity, whereas $k_{off} = 4,\lambda_{off} = 0.4$ and the other parameter values are kept unchanged for non-Markovianity. In (B) and (C), we use 1000 realizations to obtain numerical results and error bars to indicate the error ranges of the numerical results
 [1] Sören Bartels, Jakob Keck. Adaptive time stepping in elastoplasticity. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 71-88. doi: 10.3934/dcdss.2020323 [2] Emre Esentürk, Juan Velazquez. Large time behavior of exchange-driven growth. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 747-775. doi: 10.3934/dcds.2020299 [3] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020449 [4] Héctor Barge. Čech cohomology, homoclinic trajectories and robustness of non-saddle sets. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020381 [5] Ying Lin, Qi Ye. Support vector machine classifiers by non-Euclidean margins. Mathematical Foundations of Computing, 2020, 3 (4) : 279-300. doi: 10.3934/mfc.2020018 [6] Sergey Rashkovskiy. Hamilton-Jacobi theory for Hamiltonian and non-Hamiltonian systems. Journal of Geometric Mechanics, 2020, 12 (4) : 563-583. doi: 10.3934/jgm.2020024 [7] Noufel Frikha, Valentin Konakov, Stéphane Menozzi. Well-posedness of some non-linear stable driven SDEs. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 849-898. doi: 10.3934/dcds.2020302 [8] Lars Grüne, Matthias A. Müller, Christopher M. Kellett, Steven R. Weller. Strict dissipativity for discrete time discounted optimal control problems. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020046 [9] Serena Dipierro, Benedetta Pellacci, Enrico Valdinoci, Gianmaria Verzini. Time-fractional equations with reaction terms: Fundamental solutions and asymptotics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 257-275. doi: 10.3934/dcds.2020137 [10] Guido Cavallaro, Roberto Garra, Carlo Marchioro. Long time localization of modified surface quasi-geostrophic equations. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020336 [11] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [12] Awais Younus, Zoubia Dastgeer, Nudrat Ishaq, Abdul Ghaffar, Kottakkaran Sooppy Nisar, Devendra Kumar. On the observability of conformable linear time-invariant control systems. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020444 [13] Veena Goswami, Gopinath Panda. Optimal customer behavior in observable and unobservable discrete-time queues. Journal of Industrial & Management Optimization, 2021, 17 (1) : 299-316. doi: 10.3934/jimo.2019112 [14] Chongyang Liu, Meijia Han, Zhaohua Gong, Kok Lay Teo. Robust parameter estimation for constrained time-delay systems with inexact measurements. Journal of Industrial & Management Optimization, 2021, 17 (1) : 317-337. doi: 10.3934/jimo.2019113 [15] Yangrong Li, Shuang Yang, Qiangheng Zhang. Odd random attractors for stochastic non-autonomous Kuramoto-Sivashinsky equations without dissipation. Electronic Research Archive, 2020, 28 (4) : 1529-1544. doi: 10.3934/era.2020080 [16] Pengyu Chen. Non-autonomous stochastic evolution equations with nonlinear noise and nonlocal conditions governed by noncompact evolution families. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020383 [17] Lin Shi, Xuemin Wang, Dingshi Li. Limiting behavior of non-autonomous stochastic reaction-diffusion equations with colored noise on unbounded thin domains. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5367-5386. doi: 10.3934/cpaa.2020242 [18] Dan Zhu, Rosemary A. Renaut, Hongwei Li, Tianyou Liu. Fast non-convex low-rank matrix decomposition for separation of potential field data using minimal memory. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020076 [19] Bernard Bonnard, Jérémy Rouot. Geometric optimal techniques to control the muscular force response to functional electrical stimulation using a non-isometric force-fatigue model. Journal of Geometric Mechanics, 2020  doi: 10.3934/jgm.2020032 [20] Hoang The Tuan. On the asymptotic behavior of solutions to time-fractional elliptic equations driven by a multiplicative white noise. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020318

2019 Impact Factor: 1.27