We consider stochastic reaction networks modeled by continuous-time Markov chains. Such reaction networks often contain many reactions, potentially occurring at different time scales, and have unknown parameters (kinetic rates, total amounts). This makes their analysis complex. We examine stochastic reaction networks with non-interacting species that often appear in examples of interest (e.g. in the two-substrate Michaelis Menten mechanism). Non-interacting species typically appear as intermediate (or transient) chemical complexes that are depleted at a fast rate. We embed the Markov process of the reaction network into a one-parameter family under a two time-scale approach, such that molecules of non-interacting species are degraded fast. We derive simplified reaction networks where the non-interacting species are eliminated and that approximate the scaled Markov process in the limit as the parameter becomes small. Then, we derive sufficient conditions for such reductions based on the reaction network structure for both homogeneous and time-varying stochastic settings, and study examples and properties of the reduction.
Citation: Linard Hoessly, Carsten Wiuf. Fast reactions with non-interacting species in stochastic reaction networks[J]. Mathematical Biosciences and Engineering, 2022, 19(3): 2720-2749. doi: 10.3934/mbe.2022124
We consider stochastic reaction networks modeled by continuous-time Markov chains. Such reaction networks often contain many reactions, potentially occurring at different time scales, and have unknown parameters (kinetic rates, total amounts). This makes their analysis complex. We examine stochastic reaction networks with non-interacting species that often appear in examples of interest (e.g. in the two-substrate Michaelis Menten mechanism). Non-interacting species typically appear as intermediate (or transient) chemical complexes that are depleted at a fast rate. We embed the Markov process of the reaction network into a one-parameter family under a two time-scale approach, such that molecules of non-interacting species are degraded fast. We derive simplified reaction networks where the non-interacting species are eliminated and that approximate the scaled Markov process in the limit as the parameter becomes small. Then, we derive sufficient conditions for such reductions based on the reaction network structure for both homogeneous and time-varying stochastic settings, and study examples and properties of the reduction.
| [1] | J. D. Murray, Mathematical Biology I. An Introduction, volume 17 of Interdisciplinary Applied Mathematics, Springer, New York, 3 edition, (2002). doi: 10.1007/b9886 |
| [2] | E. Weinan, Principles of Multiscale Modeling, Cambridge University Press, (2011). |
| [3] |
L. Segal, M. Slemrod, The quasi-steady-state assumption: A case study in perturbation, SIAM Rev., 31 (1989), 446–477. doi: 10.1137/1031091 doi: 10.1137/1031091
|
| [4] |
N. Fenichel, Geometric singular perturbation theory for ordinary differential equations, J. Diff. Eqns., 31 (1979), 53–98. doi: 10.1016/0022-0396(79)90152-9 doi: 10.1016/0022-0396(79)90152-9
|
| [5] | A. N. Tikhonov, Systems of differential equations containing a small parameter multiplying the derivative (in Russian), Math. Sb., 31 (1952), 575–586. |
| [6] |
A. Goeke, S. Walcher, E. Zerz, Classical quasi-steady state reduction – a mathematical characterization, Phys. D Nonlinear Phenom., 345 (2017), 11–26. doi: 10.1016/j.physd.2016.12.002 doi: 10.1016/j.physd.2016.12.002
|
| [7] | F. G. Heineken, H. M. Tsuchiya, R. Aris, On the mathematical status of the pseudo-steady state hypothesis of biochemical kinetics, Math. Biosci.., 1 (1967), 95–113. |
| [8] |
H.-W. Kang, W. R. KhudaBukhsh, H. Koeppl, G. A. Rempała, Quasi-steady-state approximations derived from the stochastic model of enzyme kinetics, Bull. Math. Biol., 81 (2019), 1303–1336. doi: 10.1007/s11538-019-00574-4 doi: 10.1007/s11538-019-00574-4
|
| [9] |
E. Feliu, C. Wiuf, Variable elimination in chemical reaction networks with mass-action kinetics, SIAM J. Appl. Math., 72 (2012), 959–981. doi: 10.1137/110847305 doi: 10.1137/110847305
|
| [10] | E. Feliu, S. Walcher, C. Wiuf, Quasi-steady state and singular perturbation reduction for reaction networks with non-interacting species, SIAM J. Appl. Dyn. Syst., In press. |
| [11] |
A. Gorban, Model reduction in chemical dynamics: Slow invariant manifolds, singular perturbations, thermodynamic estimates, and analysis of reaction graph, Curr. Opin. Chem. Eng., 21 (2018), 48–59. doi: 10.1016/j.coche.2018.02.009 doi: 10.1016/j.coche.2018.02.009
|
| [12] |
X. Kan, Chang Hyeong Lee, H. G. Othmer, A multi-time-scale analysis of chemical reaction networks: Ii. stochastic systems, J. Math. Biol., 73 (2016), 1081–1129. doi: 10.1007/s00285-016-0980-x doi: 10.1007/s00285-016-0980-x
|
| [13] |
D. Schnoerr, G. Sanguinetti, R. Grima, Approximation and inference methods for stochastic biochemical kinetics—a tutorial review, J. Phys. A Math., 50 (2017), 093001. doi: 10.1088/1751-8121/aa54d9 doi: 10.1088/1751-8121/aa54d9
|
| [14] |
J. A. M. Janssen, The elimination of fast variables in complex chemical reactions. ii. mesoscopic level (reducible case), J. Stat. Phys., 57 (1989), 171–185. doi: 10.1007/BF01023639 doi: 10.1007/BF01023639
|
| [15] |
J. A. M. Janssen, The elimination of fast variables in complex chemical reactions. iii. mesoscopic level (irreducible case), J. Stat. Phys., 57 (1989), 187–198. doi: 10.1007/BF01023640 doi: 10.1007/BF01023640
|
| [16] | T. G. Kurtz, Approximation of Population Processes, Society for Industrial and Applied Mathematics, (1981). doi: 10.1137/1.9781611970333 |
| [17] |
K. Ball, T. G. Kurtz, L. Popovic, G. Rempala, Asymptotic analysis of multiscale approximations to reaction networks, Ann. Appl. Probab., 16 (2006), 1925–1961. doi: 10.1214/105051606000000420 doi: 10.1214/105051606000000420
|
| [18] |
H.-W. Kang, T. G. Kurtz, Separation of time-scales and model reduction for stochastic reaction networks, Ann. Appl. Probab., 23 (2013), 529–583. doi: 10.1214/12-AAP841 doi: 10.1214/12-AAP841
|
| [19] |
P. Pfaffelhuber, L. Popovic, Scaling limits of spatial compartment models for chemical reaction networks, Ann. Appl. Probab., 25 (2015), 3162–3208. doi: 10.1214/14-AAP1070 doi: 10.1214/14-AAP1070
|
| [20] |
D. Cappelletti, C. Wiuf, Elimination of intermediate species in multiscale stochastic reaction networks, Ann. Appl. Probab., 26 (2016), 2915–2958. doi: 10.1214/15-AAP1166 doi: 10.1214/15-AAP1166
|
| [21] |
M. Sáez, C. Wiuf, E. Feliu, Graphical reduction of reaction networks by linear elimination of species, J. Math. Biol., 74 (2017), 195–237. doi: 10.1007/s00285-016-1028-y doi: 10.1007/s00285-016-1028-y
|
| [22] | G.G. Yin, Q. Zhang, Continuous-Time Markov Chains and Applications: A Two-Time-Scale Approach, Stochastic Modelling and Applied Probability, Springer New York, (2012). doi: 10.1007/978-1-4614-4346-9 |
| [23] | D. Freedman, Approximating Countable Markov Chains, Springer New York, (2012). |
| [24] |
X. Chen, C. Jia, Limit theorems for generalized density-dependent Markov chains and bursty stochastic gene regulatory networks, J. Math. Biol., 80 (2020), 959–994. doi: 10.1007/s00285-019-01445-1 doi: 10.1007/s00285-019-01445-1
|
| [25] |
S. Be'er, M. Assaf, Rare events in stochastic populations under bursty reproduction, J. Stat. Mech–Theory E., (2016), 113501. doi: 10.1088/1742-5468/2016/11/113501 doi: 10.1088/1742-5468/2016/11/113501
|
| [26] | B. Ingalls, Mathematical Modeling in Systems Biology, Cambridge, Massachusetts: MIT Press, (2013). |
| [27] | D. F. Anderson, T. G. Kurtz, Stochastic Analysis of Biochemical Systems, Springer Publishing Company, Incorporated, (2015). doi: 10.1007/978-3-319-16895-1 |
| [28] | J. R. Norris, Markov Chains, Cambridge University Press, Cambridge, (1997). doi: 10.1017/CBO9780511810633 |
| [29] | L. Hoessly, C. Wiuf, P. Xia, On the sum of chemical reactions, (2021). arXiv: 2105.04353. |
| [30] |
D. F. Anderson, A modified next reaction method for simulating chemical systems with time dependent propensities and delays, J. Chem. Phys., 127 (2007), 214107. doi: 10.1063/1.2799998 doi: 10.1063/1.2799998
|
| [31] | A. Cornish-Bowden, Fundamentals of Enzyme Kinetics, Wiley, (2013). |
| [32] |
C. Jia, Reduction of markov chains with two-time-scale state transitions, Stochastics, 88 (2016), 73–105. doi: 10.1080/17442508.2015.1036433 doi: 10.1080/17442508.2015.1036433
|
| [33] | G. Yin, Q. Zhang, G. Badowski, Asymptotic properties of a singularly perturbed markov chain with inclusion of transient states, Ann. Appl. Probab., 10 (2000), 549–572. |
| [34] |
A. Jakubowski, A non-Skorohod topology on the Skorohod space, Electron. J. Probab., 2 (1997), 1–21. doi: 10.1214/EJP.v2-18 doi: 10.1214/EJP.v2-18
|
| [35] | C. Xu, M. C. Hansen, C. Wiuf, Dynamics of continuous time markov chains with applications to stochastic reaction networks, (2019). arXiv: 1909.12825. |
| [36] |
T. Kurtz, The relationship between stochastic and deterministic models for chemical reactions, J. Chem. Phys.., 57 (1972), 2976–2978. doi: 10.1063/1.1678692 doi: 10.1063/1.1678692
|
| [37] | G. Pavliotis, A. Stuart, Multiscale Methods: Averaging and Homogenization, volume 53. Springer, 01 (2008). doi: 10.1007/978-0-387-73829-1 |
| [38] | J. G. Kemeny, J. L. Snell, Finite Markov Chains: With a new appendix "Generalization of a Fundamental Matrix". Undergraduate Texts in Mathematics, Springer New York, (1983). |
| [39] |
C. Meyer, Stochastic complementation, uncoupling markov chains, and the theory of nearly reducible systems, SIAM Rev.., 31 (1995), 09. doi: 10.1137/1031050 doi: 10.1137/1031050
|
| [40] |
R. J. Plemmons, M-matrix characterizations.i—nonsingular m-matrices, Linear Algebra Its Appl., 18 (1977), 175–188. doi: 10.1016/0024-3795(77)90073-8 doi: 10.1016/0024-3795(77)90073-8
|
Linard Hoessly, Carsten Wiuf. Fast reactions with non-interacting species in stochastic reaction networks[J]. Mathematical Biosciences and Engineering, 2022, 19(3): 2720-2749. doi: 10.3934/mbe.2022124
DownLoad: