Biochemical systems have numerous practical applications, in particular to the study of critical intracellular processes. Frequently, biochemical kinetic models depict cellular processes as systems of chemical reactions. Many biological processes in a cell are inherently stochastic, due to the existence of some low molecular amounts. These stochastic fluctuations may have a great effect on the biochemical system’s behaviour. In such cases, stochastic models are necessary to accurately describe the system’s dynamics. Biochemical systems at the cellular level may entail many species or reactions and their mathematical models may be non-linear and with multiple scales in time. In this work, we provide a numerical technique for simplifying stochastic discrete models of well-stirred biochemical systems, which ensures that the main properties of the original system are preserved. The proposed technique employs sensitivity analysis and requires solving an optimization problem. The numerical tests on several models of practical interest show that our model reduction strategy performs very well.
McAdams, H. and Arkin, A. (1999) It’s a Noisy Business, Genetic Regulation at the Nanomolar Scale. Trends in Genetics, 15, 65-69. https://doi.org/10.1016/S0168-9525(98)01659-X
Raser, J.M. and O’Shea, E.K. (2004) Control of Stochasticity in Eukaryotic Gene Expression. Science, 304, 1811-1814. https://doi.org/10.1126/science.1098641
Gillespie, D.T. (1992) A Rigorous Derivation of the Chemical Master Equation. Physica A: Statistical Mechanics and its Applications, 188, 402-425. https://doi.org/10.1016/0378-4371(92)90283-V
Gillespie, D.T. (1976) A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions. Journal of Computational Physics, 22, 403-434. https://doi.org/10.1016/0021-9991(76)90041-3
Gillespie, D.T. (1977) Exact Stochastic Simulation of Coupled Chemical Reactions. The Journal of Physical Chemistry, 81, 2340-2361. https://doi.org/10.1021/j100540a008
Huang, H., Fairweather, M., Griffiths, J.F., Tomlin, A.S. and Brad, R.B. (2005) A Systematic Lumping Approach for the Reduction of Comprehensive Kinetic Models. Proceedings of the Combustion Institute, 30, 1309-1316. https://doi.org/10.1016/j.proci.2004.08.001
Li, G. and Rabitz, H. (1989) A General Analysis of Exact Lumping in Chemical Kinetics. Chemical Engineering Science, 44, 1413-1430. https://doi.org/10.1016/0009-2509(89)85014-6
Petzold, L.R. and Zhu, W.J. (1999) Model Reduction for Chemical Kinetics: An Optimization Approach. AIChE Journal, 45, 869-886. https://doi.org/10.1002/aic.690450418
Kumar, A., Christofides, P.D. and Daoutidis, P. (1998) Singular Perturbation Modeling of Non-linear Processes with Non-Explicit Time-Scale Separation. Chemical Engineering Science, 53, 1491-1504. https://doi.org/10.1016/S0009-2509(98)00006-2
Vora, N. and Daoutidis, P. (2001) Nonlinear Model Reduction of Chemical Reaction Systems. AIChE Journal, 47, 2320-2332. https://doi.org/10.1002/aic.690471016
Varma, A., Morbidelli, M. and Wu, H. (1999) Parametric Sensitivity in Chemical Systems. Cambridge University Press, Cambridge. https://doi.org/10.1017/CBO9780511721779
Cao, Y., Gillespie, D.T. and Petzold, L. (2005) The Slow-Scale Stochastic Simulation Algorithm. The Journal of Chemical Physics, 122, 014116-014134. https://doi.org/10.1063/1.1824902
Rao, C.V. and Arkin, A.P. (2003) Stochastic Chemical Kinetics and the Quasi-steady-state Assumption: Application to the Gillespie Algorithm. The Journal of Chemical Physics, 118, 4999-5010. https://doi.org/10.1063/1.1545446
Samant, A. and Vlachos, D. (2005) Overcoming Stiffness in Stochastic Simulation Stemming from Partial Equilibrium: A Multiscale Monte-Carlo Algorithm. Journal of Chemical Physics, 123, 144114-144122. https://doi.org/10.1063/1.2046628
Sotiropoulos, V., Contou-Carrere, M.N., Daoutidis, P. and Kaznessis, Y.N. (2009) Model Reduction of Multiscale Chemical Langevin Equations: A Numerical Case Study, IEEE/ACM Transactions on Computational Biology and Bioinformatics, 6, 470-482. https://doi.org/10.1109/TCBB.2009.23
Ilie, S. and Gholami, S. (2013) Simplifying Stochastic Mathematical Models of Biochemical Systems. Applied Mathematics, 4, 248-256. https://doi.org/10.4236/am.2013.41A038
Anderson, D.F. (2012) An Efficient Finite Difference Method for Parameter Sensitivities of Continuous Time Markov Chains. SIAM Journal on Numerical Analysis, 50, 2237-2258. https://doi.org/10.1137/110849079
Rathinam, M., Sheppard, P.W. and Khammash, M. (2010) Efficient Computation of Parameter Sensitivities of Discrete Stochastic Chemical Reaction Networks. The Journal of Chemical Physics, 132, Article No. 034103. https://doi.org/10.1063/1.3280166
Kurtz, T.G. (1982) Representation and Approximation of Counting Processes. In: Fleming, W.H. and Gorostiza, L.G., Eds., Advances in Filtering and Optimal Stochastic Control, Vol. 42, Springer, Berlin, 177-191. https://doi.org/10.1007/BFb0004537
Morshed, M., Ingalls, B. and Ilie, S. (2017) An Efficient Finite-Difference Strategy for Sensitivity Analysis of Stochastic Models of Biochemical Systems. BioSystems, 151, 43-52. https://doi.org/10.1016/j.biosystems.2016.11.006
Morshed, M., Ingalls, B. and Ilie, S. (2018) An Effective Implicit Finite-difference Method for Sensitivity Analysis of Stiff Stochastic Discrete Biochemical Systems. IET Systems Biology, 12, 123-130. https://doi.org/10.1049/iet-syb.2017.0048
Beal, L.D.R., Hill, D., Martin, R.A. and Hedengren, J.D. (2018) GEKKO Optimization Suite. Processes, 6, Article No. 106. https://doi.org/10.3390/pr6080106
Janke, T. (2011) On Reduced Models for Chemical Master Equation. Multiscale Modeling & Simulation, 9, 1646-1676. https://doi.org/10.1137/110821500
Pettigrew, M.F. and Resat, H. (2007) Multinomial Tau-leaping Method for Stochastic Kinetic Simulations. The Journal of Chemical Physics, 126, Article ID: 084101. https://doi.org/10.1063/1.2432326
Sayyidmousavi, A. and Ilie, S. (2017) An Efficient Hybrid Method for Stochastic Reaction-Diffusion Biochemical Systems with Delay. AIP Advances, 7, Article ID: 125305. https://doi.org/10.1063/1.5001760
Kahramanogullari, O., Fantaccini, G., Lecca, P., Morpurgo, D. and Priami, C. (2012) Algorithmic Modeling Quantifies the Complementary Contribution of Metabolic Inhibitions to Gemcitabine Efficacy. PLoS ONE, 7, e50176. https://doi.org/10.1371/journal.pone.0050176
Thanh, V.H. (2019) A Critical Comparison of Rejection-Based Algorithms for Simulation of Large Biochemical Reaction Networks. Bulletin of Mathematical Biology, 81, 3053-3073. https://doi.org/10.1007/s11538-018-0462-y