Abstract
Virus infection dynamics is governed by the processes on multiple scales: on the whole organism level, tissue level, and intracellular level. In this paper, we develop a multi-scale multi-compartment model of HIV infection in a simplified setting and the computational methods for numerical realization of the model. The multiscale model describes the processes from various scales and of different nature (cell motility, virus diffusion, intracellular virus replication). Intracellular replication model is based on a Markov chain with time-inhomogeneous propensities that depend on the extracellular level of virions. Reaction diffusion equations used to model free virion diffusion in the lymphoid tissue have moving sources, which are determined by the positions of the infected cells (immune cell motility model) and the rate of virion secretion from them (intracellular model). Immune cell motility model parameterizes the intercellular interaction forces, friction and the stochastic force of active cell motility. Together, this allows for a proper description of the intracellular stochasticity that propagates across multiple scales. A hybrid discrete-continuous stochastic-deterministic algorithm for simulation of the multiscale model based on the uniformization Monte Carlo method is implemented.
Funding statement: The research was supported by the Russian Science Foundation grant number 18-11-00171.
Acknowledgment
The author is grateful to Gennady A. Bocharov for fruitful discussions and valuable comments.
References
[1] D. F. Anderson, A modified next reaction method for simulating chemical systems with time dependent propensities and delays. J. Chem. Phys 127 (2007), 21410710.1063/1.2799998Search in Google Scholar PubMed
[2] A. Auger, P. Chatelain, and P. Koumoutsakos, R-leaping: Accelerating the stochastic simulation algorithm by reaction leaps. J. Chem. Phys 125 (2006), 08410310.1063/1.2218339Search in Google Scholar PubMed
[3] T. A. Averina, Statistical simulation methods for a nonhomogeneous Poisson ensemble. Numer. Anal. Appl 2 (2009), No. 4, 289–30110.1134/S1995423909040016Search in Google Scholar
[4] T. A. Averina and G. A. Mikhailov, Algorithms for exact and approximate statistical simulation of Poisson ensembles. Comp. Math. Math. Phys 50 (2010), No. 6, 951–96210.1134/S0965542510060035Search in Google Scholar
[5] C. H. L. Beentjes and R. E. Baker, Uniformization techniques for stochastic simulation of chemical reaction networks. J. Chem. Phys 150 (2019), 15410710.1063/1.5081043Search in Google Scholar PubMed
[6] H. Belanger, D. Mancusi, and A. Zoia, Review of Monte Carlo methods for particle transport in continuously-varying media. Europ. Phys. J. Plus 135 (2020), 87710.1140/epjp/s13360-020-00731-ySearch in Google Scholar
[7] J. L. Blue, I. Beichl, and F. Sullivan, Faster Monte Carlo simulations. Phys. Review E 51 (1995), R867–R86810.1103/PhysRevE.51.R867Search in Google Scholar PubMed
[8] G. Bocharov, A. Danilov, Y. Vassilevski, G. Marchuk, V. Chereshnev, and B. Ludewig, Reaction–diffusion modelling of interferon distribution in secondary lymphoid organs. Math. Modelling Natur. Phenom 6 (2011), 13–2610.1051/mmnp/20116702Search in Google Scholar
[9] G. A. Bocharov, D. S. Grebennikov, and R. S. Savinkov, Mathematical immunology: from phenomenological to multiphysics modelling. Russ. J. Numer. Anal. Math. Modelling 35 (2020), No. 4, 203–21310.1515/rnam-2020-0017Search in Google Scholar
[10] A. Bortz, M. Kalos, and J. Lebowitz, A new algorithm for Monte Carlo simulation of Ising spin systems. J. Comput. Phys 17 (1975), 10–1810.1016/0021-9991(75)90060-1Search in Google Scholar
[11] A. Bouchnita, G. Bocharov, A. Meyerhans, and V. Volpert, Towards a multiscale model of acute HIV infection. Computation 5 (2017), No. 1, 610.3390/computation5010006Search in Google Scholar
[12] X. Cai and Z. Xu, K-leap method for accelerating stochastic simulation of coupled chemical reactions. J. Chem. Phys 126 (2007), 07410210.1063/1.2436869Search in Google Scholar PubMed
[13] Y. Cao, H. Li, and L. Petzold, Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. J. Chem. Phys 121 (2004), 4059–406710.1063/1.1778376Search in Google Scholar PubMed
[14] W. A. Coleman, Mathematical verification of a certain Monte Carlo sampling technique and applications of the technique to radiation transport problems. Nuclear Sci. Engrg 32 (1968), 76–8110.13182/NSE68-1Search in Google Scholar
[15] M. A. Gibson and J. Bruck, Efficient exact stochastic simulation of chemical systems with many species and many channels. J. Phys. Chem. A 104 (2000), 1876–188910.1021/jp993732qSearch in Google Scholar
[16] D. T. Gillespie, A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. J. Comput. Phys 22 (1976), 403–43410.1016/0021-9991(76)90041-3Search in Google Scholar
[17] D. T. Gillespie, The chemical Langevin equation. J. Chem. Phys 113 (2000), 297–30610.1063/1.481811Search in Google Scholar
[18] D. T. Gillespie, Approximate accelerated stochastic simulation of chemically reacting systems. J. Chem. Phys 115 (2001), 1716–173310.1063/1.1378322Search in Google Scholar
[19] D. T. Gillespie, Stochastic simulation of chemical kinetics. Annual Review of Physical Chemistry 58 (2007), 35–5510.1146/annurev.physchem.58.032806.104637Search in Google Scholar PubMed
[20] D. T. Gillespie and L. R. Petzold, Improved leap-size selection for accelerated stochastic simulation. J. Chem. Phys 119 (2003), 8229–823410.1063/1.1613254Search in Google Scholar
[21] D. Grebennikov, A. Bouchnita, V. Volpert, N. Bessonov, A. Meyerhans, and G. Bocharov, Spatial limphocyte dynamics in limph nodes predicts the cytotoxic T cell frequency needed for HIV infection control. Frontiers in Immunology 10 (2019), 0121310.3389/fimmu.2019.01213Search in Google Scholar PubMed PubMed Central
[22] D. S. Grebennikov and G. A. Bocharov, Spatially resolved modelling of immune responses following a multiscale approach: from computational implementation to quantitative predictions. Russ. J. Numer. Anal. Math. Modelling 34 (2019), No. 5, 253–26010.1515/rnam-2019-0021Search in Google Scholar
[23] T. G. Kurtz, The relationship between stochastic and deterministic models for chemical reactions. J. Chem. Phys 57 (1972), 2976–297810.1063/1.1678692Search in Google Scholar
[24] P. A. W. Lewis and G. S. Shedler, Simulation of nonhomogeneous poisson processes by thinning. Naval Research Logistics Quarterly 26 (1979), 403–41310.1002/nav.3800260304Search in Google Scholar
[25] J. M. McCollum, G. D. Peterson, C. D. Cox, M. L. Simpson, and N. F. Samatova, The sorting direct method for stochastic simulation of biochemical systems with varying reaction execution behavior. Comput. Biology and Chemistry 30 (2006), 39–4910.1016/j.compbiolchem.2005.10.007Search in Google Scholar PubMed
[26] H. Resat, H. S. Wiley, and D. A. Dixon, Probability-weighted dynamic Monte Carlo method for reaction kinetics simulations. J. Phys. Chem. B 105 (2001), 11026–1103410.1021/jp011404wSearch in Google Scholar
[27] R. Savinkov, D. Grebennikov, D. Puchkova, V. Chereshnev, I. Sazonov, and G. Bocharov, Graph theory for modelling and analysis of the human lymphatic system. Mathematics 8 (2020), 223610.3390/math8122236Search in Google Scholar
[28] I. Sazonov, D. Grebennikov, M. Kelbert, and G. Bocharov, Modelling stochastic and deterministic behaviours in virus infection dynamics. Math. Modelling Natur. Phenom 12 (2017), 63–7710.1051/mmnp/201712505Search in Google Scholar
[29] I. Sazonov, D. Grebennikov, M. Kelbert, A. Meyerhans, and G. Bocharov, Viral infection dynamics model based on a Markov process with time delay between cell infection and progeny production. Mathematics 8 (2020), 120710.3390/math8081207Search in Google Scholar
[30] I. Sazonov, D. Grebennikov, A. Meyerhans, and G. Bocharov, Markov chain-based stochastic modelling of HIV-1 life cycle in a CD4 T cell. Mathematics 9 (2021), 202510.3390/math9172025Search in Google Scholar
[31] I. Sazonov, D. Grebennikov, R. Savinkov, A. Soboleva, K. Pavlishin, A. Meyerhans, and G. Bocharov, Stochastic modelling of HIV-1 replication in a CD4 T cell with an IFN response. Viruses 15 (2023), 29610.3390/v15020296Search in Google Scholar PubMed PubMed Central
[32] T. P. Schulze, Efficient kinetic Monte Carlo simulation. J. Comput. Phys 227 (2008), 2455–246210.1016/j.jcp.2007.10.021Search in Google Scholar
[33] S. A. Serebrinsky, Physical time scale in kinetic Monte Carlo simulations of continuous-time Markov chains. Phys. Review E 83 (2011), 03770110.1103/PhysRevE.83.037701Search in Google Scholar PubMed
[34] O. Shcherbatova, D. Grebennikov, I. Sazonov, A. Meyerhans, and G. Bocharov, Modelling of the HIV-1 life cycle in productively infected cells to predict novel therapeutic targets. Pathogens 9 (2020), 25510.3390/pathogens9040255Search in Google Scholar PubMed PubMed Central
[35] A. Slepoy, A. P. Thompson, and S. J. Plimpton, A constant-time kinetic Monte Carlo algorithm for simulation of large biochemical reaction networks. J. Chem. Phys 128 (2008), 20510110.1063/1.2919546Search in Google Scholar PubMed
[36] V. H. Thanh and C. Priami, Simulation of biochemical reactions with time-dependent rates by the rejection-based algorithm. J. Chem. Phys 143 (2015), 05410410.1063/1.4927916Search in Google Scholar PubMed
[37] V. H. Thanh, C. Priami, and R. Zunino, Efficient rejection-based simulation of biochemical reactions with stochastic noise and delays. J. Chem. Phys 141 (2014), 13411610.1063/1.4896985Search in Google Scholar PubMed
[38] V. H. Thanh, R. Zunino, and C. Priami, On the rejection-based algorithm for simulation and analysis of large-scale reaction networks. J. Chem. Phys 142 (2015), 24410610.1063/1.4922923Search in Google Scholar PubMed
[39] V. H. Thanh, C. Priami, and R. Zunino, Accelerating rejection-based simulation of biochemical reactions with bounded acceptance probability. J. Chem. Phys 144 (2016), 22410810.1063/1.4953559Search in Google Scholar PubMed
[40] A. P. A. Van Moorsel and W. H. Sanders, Adaptive uniformization. Communications in Statistics. Stochastic Models 10 (1994), 619–64710.1080/15326349408807313Search in Google Scholar
[41] N. van Dijk, S. van Brummelen, and R. Boucherie, Uniformization: Basics, extensions and applications. Performance Evaluation 118 (2018), 8–3210.1016/j.peva.2017.09.008Search in Google Scholar
[42] C. L. Vestergaard and M. Génois. Temporal Gillespie algorithm: Fast simulation of contagion processes on time-varying networks. PLOS Computational Biology 11 (2015), e100457910.1371/journal.pcbi.1004579Search in Google Scholar PubMed PubMed Central
[43] M. Voliotis, P. Thomas, R. Grima, and C. G. Bowsher, Stochastic simulation of biomolecular networks in dynamic environments. PLOS Computational Biology 12 (2016), e100492310.1371/journal.pcbi.1004923Search in Google Scholar PubMed PubMed Central
© 2023 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Contents
- Effect of electron temperature on formation of travelling waves in plasma: Kinetic and hydrodynamic models
- Computational methods for multiscale modelling of virus infection dynamics
- Dependence of optimal disturbances on periodic solution phases for time-delay systems
- Sketching for a low-rank nonnegative matrix approximation: Numerical study
Articles in the same Issue
- Contents
- Effect of electron temperature on formation of travelling waves in plasma: Kinetic and hydrodynamic models
- Computational methods for multiscale modelling of virus infection dynamics
- Dependence of optimal disturbances on periodic solution phases for time-delay systems
- Sketching for a low-rank nonnegative matrix approximation: Numerical study