# 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] Ritu Agarwal, Kritika, Sunil Dutt Purohit, Devendra Kumar. Mathematical modelling of cytosolic calcium concentration distribution using non-local fractional operator. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021017 [2] Xiaohong Li, Mingxin Sun, Zhaohua Gong, Enmin Feng. Multistage optimal control for microbial fed-batch fermentation process. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021040 [3] Andrey Kovtanyuk, Alexander Chebotarev, Nikolai Botkin, Varvara Turova, Irina Sidorenko, Renée Lampe. Modeling the pressure distribution in a spatially averaged cerebral capillary network. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021016 [4] Meiqiao Ai, Zhimin Zhang, Wenguang Yu. First passage problems of refracted jump diffusion processes and their applications in valuing equity-linked death benefits. Journal of Industrial & Management Optimization, 2021  doi: 10.3934/jimo.2021039 [5] Alberto Bressan, Ke Han, Franco Rampazzo. On the control of non holonomic systems by active constraints. Discrete & Continuous Dynamical Systems - A, 2013, 33 (8) : 3329-3353. doi: 10.3934/dcds.2013.33.3329 [6] Jianping Gao, Shangjiang Guo, Wenxian Shen. Persistence and time periodic positive solutions of doubly nonlocal Fisher-KPP equations in time periodic and space heterogeneous media. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2645-2676. doi: 10.3934/dcdsb.2020199 [7] Cécile Carrère, Grégoire Nadin. Influence of mutations in phenotypically-structured populations in time periodic environment. Discrete & Continuous Dynamical Systems - B, 2020, 25 (9) : 3609-3630. doi: 10.3934/dcdsb.2020075 [8] Paula A. González-Parra, Sunmi Lee, Leticia Velázquez, Carlos Castillo-Chavez. A note on the use of optimal control on a discrete time model of influenza dynamics. Mathematical Biosciences & Engineering, 2011, 8 (1) : 183-197. doi: 10.3934/mbe.2011.8.183 [9] Guillermo Reyes, Juan-Luis Vázquez. Long time behavior for the inhomogeneous PME in a medium with slowly decaying density. Communications on Pure & Applied Analysis, 2009, 8 (2) : 493-508. doi: 10.3934/cpaa.2009.8.493 [10] Xinyuan Liao, Caidi Zhao, Shengfan Zhou. Compact uniform attractors for dissipative non-autonomous lattice dynamical systems. Communications on Pure & Applied Analysis, 2007, 6 (4) : 1087-1111. doi: 10.3934/cpaa.2007.6.1087 [11] Wei-Jian Bo, Guo Lin, Shigui Ruan. Traveling wave solutions for time periodic reaction-diffusion systems. Discrete & Continuous Dynamical Systems - A, 2018, 38 (9) : 4329-4351. doi: 10.3934/dcds.2018189 [12] Kin Ming Hui, Soojung Kim. Asymptotic large time behavior of singular solutions of the fast diffusion equation. Discrete & Continuous Dynamical Systems - A, 2017, 37 (11) : 5943-5977. doi: 10.3934/dcds.2017258 [13] Pascal Noble, Sebastien Travadel. Non-persistence of roll-waves under viscous perturbations. Discrete & Continuous Dynamical Systems - B, 2001, 1 (1) : 61-70. doi: 10.3934/dcdsb.2001.1.61 [14] Linlin Li, Bedreddine Ainseba. Large-time behavior of matured population in an age-structured model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2561-2580. doi: 10.3934/dcdsb.2020195 [15] Vieri Benci, Marco Cococcioni. The algorithmic numbers in non-archimedean numerical computing environments. Discrete & Continuous Dynamical Systems - S, 2021, 14 (5) : 1673-1692. doi: 10.3934/dcdss.2020449 [16] Liqin Qian, Xiwang Cao. Character sums over a non-chain ring and their applications. Advances in Mathematics of Communications, 2021  doi: 10.3934/amc.2020134 [17] Zhiming Guo, Zhi-Chun Yang, Xingfu Zou. Existence and uniqueness of positive solution to a non-local differential equation with homogeneous Dirichlet boundary condition---A non-monotone case. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1825-1838. doi: 10.3934/cpaa.2012.11.1825 [18] Fumihiko Nakamura. Asymptotic behavior of non-expanding piecewise linear maps in the presence of random noise. Discrete & Continuous Dynamical Systems - B, 2018, 23 (6) : 2457-2473. doi: 10.3934/dcdsb.2018055 [19] Tomáš Roubíček. An energy-conserving time-discretisation scheme for poroelastic media with phase-field fracture emitting waves and heat. Discrete & Continuous Dynamical Systems - S, 2017, 10 (4) : 867-893. doi: 10.3934/dcdss.2017044 [20] Xiaomao Deng, Xiao-Chuan Cai, Jun Zou. A parallel space-time domain decomposition method for unsteady source inversion problems. Inverse Problems & Imaging, 2015, 9 (4) : 1069-1091. doi: 10.3934/ipi.2015.9.1069

2019 Impact Factor: 1.27