The most common viral dynamics models for analyzing viral infections assume a well-mixed spatial distribution between viral particles and uninfected target cells. However, throughout an infection, the spatial distribution of virus and cells changes. Initially, virus and infected cells are localized so that a target cell in an area with lower virus presence will be less likely to be infected than a cell close to a location of viral production. A density-dependent infection rate has the potential to improve models that treat the cellular infection probability as a constant. Building on previous work that used epidemiological models, which introduced Saturated Incidence, Beddington-DeAngelis, and Crowley Martin models to simulate spatial heterogeneity in population-level models, we implemented these density-dependent models in within-host viral kinetics models to understand how a density dependent infection rate might impact the predicted severity of an influenza infection. The parameter values that govern the strength of density dependence were varied to understand the implications of density dependent infection rate on disease severity and the potential impacts on patient treatment. For low density dependence, a steeper increase in the virus and a higher viral peak was predicted. We find that the initial localization of infected cells slows the progression of the infection, thus suggesting that accounting for density dependence when analyzing influenza infection severity can result in an altered expectation for viral progression.
Citation: Hope Sage, Hana M. Dobrovolny. Exploring the effect of density dependent infection rate in a within host viral infection model[J]. Mathematical Biosciences and Engineering, 2026, 23(6): 1687-1707. doi: 10.3934/mbe.2026061
The most common viral dynamics models for analyzing viral infections assume a well-mixed spatial distribution between viral particles and uninfected target cells. However, throughout an infection, the spatial distribution of virus and cells changes. Initially, virus and infected cells are localized so that a target cell in an area with lower virus presence will be less likely to be infected than a cell close to a location of viral production. A density-dependent infection rate has the potential to improve models that treat the cellular infection probability as a constant. Building on previous work that used epidemiological models, which introduced Saturated Incidence, Beddington-DeAngelis, and Crowley Martin models to simulate spatial heterogeneity in population-level models, we implemented these density-dependent models in within-host viral kinetics models to understand how a density dependent infection rate might impact the predicted severity of an influenza infection. The parameter values that govern the strength of density dependence were varied to understand the implications of density dependent infection rate on disease severity and the potential impacts on patient treatment. For low density dependence, a steeper increase in the virus and a higher viral peak was predicted. We find that the initial localization of infected cells slows the progression of the infection, thus suggesting that accounting for density dependence when analyzing influenza infection severity can result in an altered expectation for viral progression.
| [1] |
M. A. Rolfes, I. M. Foppa, S. Garg, B. Flannery, L. Brammer, J. A. Singleton, et al., Annual estimates of the burden of seasonal influenza in the united states: A tool for strengthening influenza surveillance and preparedness, Influenza Other Respir. Viruses, 12 (2018), 132–137. https://doi.org/10.1111/irv.12486 doi: 10.1111/irv.12486
|
| [2] | Centers for Disease Control and Prevention, Seasonal Flu Microsite, 2022. Available from: https://www.cdc.gov/flu/season/2022-2023.html. |
| [3] |
C. de Courville, S. M. Cadarette, E. Wissinger, F. P. Alvarez, The economic burden of influenza among adults aged 18 to 64: A systematic literature review, Influenza Other Respir. Viruses, 16 (2022), 376–385. https://doi.org/10.1111/irv.12963 doi: 10.1111/irv.12963
|
| [4] |
P. Baccam, C. Beauchemin, C. A. Macken, F. G. Hayden, A. S. Perelson, Kinetics of influenza a virus infection in humans, J. Virol., 80 (2006), 7590–7599. https://doi.org/10.1128/JVI.01623-05 doi: 10.1128/JVI.01623-05
|
| [5] |
G. González-Parra, H. M. Dobrovolny, Assessing uncertainty in A2 respiratory syncytial virus viral dynamics, Comput. Math. Meth. Med., 2015 (2015), 567589. https://doi.org/10.1155/2015/567589 doi: 10.1155/2015/567589
|
| [6] |
C. P. de Mello, X. Tao, T. H. Kim, M. Vicchiarelli, J. B. Bulitta, A. Kaushik, et al., Clinical regimens of favipiravir inhibit zika virus replication in the hollow-fiber infection model, Antimicrob. Agents Chemother., 62 (2018), e00967–18. https://doi.org/10.1128/AAC.00967-18 doi: 10.1128/AAC.00967-18
|
| [7] |
E. A. Hernandez-Vargas, J. X. Velasco-Hernandez, In-host modelling of COVID-19 kinetics in humans, Ann. Rev. Contr., 50 (2020), 448–456. https://doi.org/10.1101/2020.03.26.20044487 doi: 10.1101/2020.03.26.20044487
|
| [8] |
V. K. Nguyen, S. C. Binder, A. Boianelli, M. Meyer-Hermann, E. A. Hernandez-Vargas, Ebola virus infection modeling and identifiability problems, Front. Microbiol., 6 (2015), 257. https://doi.org/10.3389/fmicb.2015.00257 doi: 10.3389/fmicb.2015.00257
|
| [9] |
C. Zitzmann, C. Daechert, B. Schmid, H. Van der Schaar, M. van Hemert, A. S. Perelson, et al., Mathematical modeling of plus-strand RNA virus replication to identify broad-spectrum antiviral treatment strategies, PLoS Comput. Biol., 19 (2023), e1010423. https://doi.org/10.1371/journal.pcbi.1010423 doi: 10.1371/journal.pcbi.1010423
|
| [10] |
K. Melville, T. Rodriguez, H. M. Dobrovolny, Investigating different mechanisms of action in combination therapy for influenza, Front. Pharmacol., 9 (2018), 1207. https://doi.org/10.3389/fphar.2018.01207 doi: 10.3389/fphar.2018.01207
|
| [11] |
C. A. Beauchemin, J. J. McSharry, G. L. Drusano, J. T. Nguyen, G. T. Went, R. M. Ribeiro, et al., Modeling amantadine treatment of influenza A virus in vitro, J. Theor. Biol., 254 (2008), 439–451. https://doi.org/10.1016/j.jtbi.2008.05.031 doi: 10.1016/j.jtbi.2008.05.031
|
| [12] |
C. A. A. Beauchemin, A. Handel, A review of mathematical models of influenza a infections within a host or cell culture: Lessons learned and challenges ahead, BMC Public Health, 11 (2011), S7. https://doi.org/10.1186/1471-2458-11-S1-S7 doi: 10.1186/1471-2458-11-S1-S7
|
| [13] |
H. M. Dobrovolny, M. B. Reddy, M. A. Kamal, C. R. Rayner, C. A. Beauchemin, Assessing mathematical models of influenza infections using features of the immune response, PLoS One, 8 (2013), e57088. https://doi.org/10.1371/journal.pone.0057088 doi: 10.1371/journal.pone.0057088
|
| [14] |
C. Zitzmann, L. Kaderali, Mathematical analysis of viral replication dynamics and antiviral treatment strategies: From basic models to age-based multi-scale modeling, Front. Microbiol., 9 (2018), 1546. https://doi.org/doi:10.3389/fmicb.2018.01546 doi: 10.3389/fmicb.2018.01546
|
| [15] |
E. A. Hernandez-Vargas, E. Wilk, L. Canini, F. R. Toapanta, S. C. Binder, A. Uvarovskii, et al., Effects of aging on influenza virus infection dynamics, J. Virol., 88 (2014), 4123–4131. https://doi.org/10.1128/JVI.03644-13 doi: 10.1128/JVI.03644-13
|
| [16] |
L. Canini, A. S. Perelson, Viral kinetic modeling: State of the art, J. Theor. Biol., 242 (2006), 464–477. https://doi.org/10.1016/j.jtbi.2006.03.014 doi: 10.1016/j.jtbi.2006.03.014
|
| [17] |
B. Fain, H. M. Dobrovolny, Initial inoculum and the severity of COVID-19: A mathematical modeling study of the dose-response of SARS-CoV-2 infections, Epidemiologia, 1 (2020), 5–15. https://doi.org/10.3390/epidemiologia1010003 doi: 10.3390/epidemiologia1010003
|
| [18] |
D. S. Quan, Y. Weiming, Mathematical modeling of interaction between innate and adaptive immune responses in COVID-19 and implications for viral pathogenesis, J. Med. Virol., 92 (2020), 1615–1628. https://doi.org/10.1002/jmv.25866 doi: 10.1002/jmv.25866
|
| [19] |
B. P. Holder, C. A. Beauchemin, Exploring the effect of biological delays in kinetic models of influenza within a host or cell culture, BMC Public Health, 11 (2011), S10. https://doi.org/10.1186/1471-2458-11-S1-S10 doi: 10.1186/1471-2458-11-S1-S10
|
| [20] |
Y. Kakizoe, S. Nakaoka, C. A. Beauchemin, S. Morita, H. Mori, T. Igarashi, et al., A method to determine the duration of the eclipse phase for in vitro infection with a highly pathogenic SHIV strain, Sci. Rep., 5 (2015), 10371. https://doi.org/10.1038/srep10371 doi: 10.1038/srep10371
|
| [21] |
B. Fain, H. M. Dobrovolny, GPU acceleration and data fitting: Agent-based models of viral infections can now be parameterized in hours, J. Comput. Sci., 61 (2022), 101662. https://doi.org/10.1016/j.jocs.2022.101662 doi: 10.1016/j.jocs.2022.101662
|
| [22] |
T. J. Sego, E. D. Mochan, G. B. Ermentrout, J. A. Glazier, A multiscale multicellular spatiotemporal model of local influenza infection and immune response, J. Theor. Biol., 532 (2022), 110918. https://doi.org/10.1016/j.jtbi.2021.110918 doi: 10.1016/j.jtbi.2021.110918
|
| [23] |
B. P. Holder, L. E. Liao, P. Simon, G. Boivin, C. A. A. Beauchemin, Design considerations in building in silico equivalents of common experimental influenza virus assays and the benefits of such an approach, Autoimmunity, 44 (2011), 282–293. https://doi.org/10.3109/08916934.2011.523267 doi: 10.3109/08916934.2011.523267
|
| [24] |
M. Keeling, The effects of local spatial structure on epidemiological invasions, Proc. R. Soc. B Biol. Sci., 266 (1999), 859–867. https://doi.org/10.1098/rspb.1999.0716 doi: 10.1098/rspb.1999.0716
|
| [25] |
S. H. Paull, S. Song, K. M. McClure, L. C. Sackett, A. M. Kilpatrick, P. T. J. Johnson, From superspreaders to disease hotspots: linking transmission across hosts and space, Front. Ecol. Environ., 10 (2012), 75–82. https://doi.org/10.1890/110111 doi: 10.1890/110111
|
| [26] |
E. S. Nightingale, O. J. Brady, L. Yakob, C. C. W. Grp, The importance of saturating density dependence for population-level predictions of SARS-CoV-2 resurgence compared with density-independent or linearly density-dependent models, England, 23 March to 31 July 2020, Eurosurveillance, 26 (2021), 2001809. https://doi.org/10.2807/1560-7917.ES.2021.26.49.2001809 doi: 10.2807/1560-7917.ES.2021.26.49.2001809
|
| [27] |
M. Hochberg, Nonlinear transmission rates and the dynamics of infectious disease, J. Theor. Biol., 153 (1991), 301–321. https://doi.org/10.1016/S0022-5193(05)80572-7 doi: 10.1016/S0022-5193(05)80572-7
|
| [28] |
H. Laarabi, E. H. Labriji, M. Rachik, A. Kaddar, Optimal control of an epidemic model with a saturated incidence rate, Nonlinear Anal. Model. Control, 17 (2012), 448–459. https://doi.org/10.15388/NA.17.4.14050 doi: 10.15388/NA.17.4.14050
|
| [29] |
A. Kumar, Dynamic behavior of an SIR epidemic model along with time delay; Crowley-Martin type incidence rate and Holling type Ⅱ treatment rate, Int. J. Nonlinear Sci. Numer. Simul., 20 (2019), 757–771. https://doi.org/10.1515/ijnsns-2018-0208 doi: 10.1515/ijnsns-2018-0208
|
| [30] |
X. Wang, Y. Tao, X. Song, A delayed HIV-1 infection model with Beddington-DeAngelis functional response, Nonlinear Dyn., 62 (2010), 67–72. https://doi.org/10.1007/s11071-010-9699-1 doi: 10.1007/s11071-010-9699-1
|
| [31] |
K. Hattaf, N. Yousfi, A. Tridane, Mathematical analysis of a virus dynamics model with general incidence rate and cure rate, Nonlinear Anal. Real World Appl., 13 (2012), 1866–1872. https://doi.org/10.1016/j.nonrwa.2011.12.015 doi: 10.1016/j.nonrwa.2011.12.015
|
| [32] |
M. K. Singh, A. Anjali, B. K. Singh, C. Cattani, Impact of general incidence function on three-strain SEIAR model, Math. Biosci. Eng., 20 (2023), 19710–19731. https://doi.org/10.3934/mbe.2023873 doi: 10.3934/mbe.2023873
|
| [33] |
H. Hethcote, P. Vandendriessche, Some epidemiologic models with nonlinear incidence, J. Math. Biol., 29 (1991), 271–287. https://doi.org/10.1007/BF00160539 doi: 10.1007/BF00160539
|
| [34] |
G. Huang, Y. Takeuchi, W. Ma, D. Wei, Global stability for delay SIR and SEIR epidemic models with nonlinear incidence rate, Bull. Math. Biol., 72 (2010), 1192–1207. https://doi.org/10.1007/s11538-009-9487-6 doi: 10.1007/s11538-009-9487-6
|
| [35] |
A. P. Smith, D. J. Moquin, V. Bernhauerova, A. M. Smith, Influenza virus infection model with density dependence supports biphasic viral decay, Front. Microbiol., 9 (2018), 1554. https://doi.org/10.3389/fmicb.2018.01554 doi: 10.3389/fmicb.2018.01554
|
| [36] |
A. Gerg, H. M. Dobrovolny, Quantifying impact of HIV receptor surface density reveals differences in fusion dynamics of HIV strains, Viruses, 17 (2025), 583. https://doi.org/10.3390/v17040583 doi: 10.3390/v17040583
|
| [37] |
A. Farrell, T. Phan, C. B. Brooke, K. Koelle, R. Ke, Semi-infectious particles contribute substantially to influenza virus within-host dynamics when infection is dominated by spatial structure, Virus Evol., 9 (2023), vead020. https://doi.org/10.1093/ve/vead020 doi: 10.1093/ve/vead020
|
| [38] |
A. M. Elaiw, S. A. Azoz, Global properties of a class of HIV infection models with beddington-deangelis functional response, Math. Method Appl. Sci., 36 (2013), 383–394. https://doi.org/10.1002/mma.2596 doi: 10.1002/mma.2596
|
| [39] |
G. Huang, W. Ma, Y. Takeuchi, Global analysis for delay virus dynamics model with beddington–deangelis functional response, Appl. Math. Lett., 24 (2011), 1199–1203. https://doi.org/https://doi.org/10.1016/j.aml.2011.02.007 doi: 10.1016/j.aml.2011.02.007
|
| [40] |
Z. T. H. Miao, X. Abdurahman, Stability and hopf bifurcation for a five-dimensional virus infection model with beddington–deangelis incidence and three delays, J. Biol. Dyn., 12 (2018), 146–170. https://doi.org/10.1080/17513758.2017.1408861 doi: 10.1080/17513758.2017.1408861
|
| [41] |
G. Huang, W. Ma, Y. Takeuchi, Global properties for virus dynamics model with beddington–deangelis functional response, Appl. Math. Lett., 22 (2009), 1690–1693. https://doi.org/https://doi.org/10.1016/j.aml.2009.06.004 doi: 10.1016/j.aml.2009.06.004
|
| [42] |
H. Hu, K. Nigmatulina, P. Eckhoff, The scaling of contact rates with population density for the infectious disease models, Math. Biosci., 244 (2013), 125–134. https://doi.org/10.1016/j.mbs.2013.04.013 doi: 10.1016/j.mbs.2013.04.013
|
| [43] |
J. Tripathi, S. Bugalia, V. Tiwari, Y. Kang, A predator–prey model with crowley–martin functional response: a nonautonomous study, Nat. Resour. Model., 33 (2020), e12287. https://doi.org/10.1111/nrm.12287 doi: 10.1111/nrm.12287
|
| [44] |
N. Kumar, Dynamic behavior of an SIR epidemic model along with time delay; crowley-martin type incidence rate and holling type Ⅱ treatment rate, Int. J. Nonlinear Sci. Num., 20 (2019), 757–771. https://doi.org/10.1515/ijnsns-2018-0208 doi: 10.1515/ijnsns-2018-0208
|
| [45] |
O. Diekmann, J. A. P. Heesterbeek, M. G. Roberts, The construction of next-generation matrices for compartmental epidemic models, J. R. Soc. Interface, 7 (2010), 873–885. https://doi.org/10.1098/rsif.2009.0386 doi: 10.1098/rsif.2009.0386
|
| [46] |
H. M. Dobrovolny, R. Gieschke, B. E. Davies, N. L. Jumbe, C. A. A. Beauchemin, Neuraminidase inhibitors for treatment of human and avian strain influenza: A comparative study, J. Theor. Biol., 269 (2011), 234–244. https://doi.org/10.1016/j.jtbi.2010.10.017 doi: 10.1016/j.jtbi.2010.10.017
|
| [47] |
G. González-Parra, F. De Ridder, D. Huntjens, D. Roymans, G. Ispas, H. M. Dobrovolny, A comparison of RSV and influenza in vitro kinetic parameters reveals differences in infecting time, PLoS One, 13 (2018), e0192645. https://doi.org/10.1371/journal.pone.0192645 doi: 10.1371/journal.pone.0192645
|
| [48] |
L. Pinky, H. M. Dobrovolny, The impact of cell regeneration on the dynamics of viral coinfection, Chaos, 27 (2017), 063109. https://doi.org/10.1063/1.4985276 doi: 10.1063/1.4985276
|
| [49] |
A. M. Smith, F. R. Adler, A. S. Perelson, An accurate two-phase approximate solution to an acute viral infection model, J. Math. Biol., 60 (2010), 711–726. https://doi.org/10.1007/s00285-009-0281-8 doi: 10.1007/s00285-009-0281-8
|
| [50] |
N. F. Beggs, H. M. Dobrovolny, Determining drug efficacy parameters for mathematical models of influenza, J. Biol. Dyn., 9 (2015), 332–346. https://doi.org/10.1080/17513758.2015.1052764 doi: 10.1080/17513758.2015.1052764
|
| [51] |
J. Palmer, H. M. Dobrovolny, C. A. Beauchemin, The in vivo efficacy of neuraminidase inhibitors cannot be determined from the decay rates of influenza viral titers observed in treated patients, Sci. Rep., 7 (2017), 40210. https://doi.org/10.1038/srep40210 doi: 10.1038/srep40210
|
| [52] |
M. E. Gallagher, C. B. Brooke, R. Ke, K. Koelle, Causes and consequences of spatial within-host viral spread, Viruses, 10 (2018), 627. https://doi.org/10.3390/v10110627 doi: 10.3390/v10110627
|
| [53] |
C. Beauchemin, Probing the effects of the well-mixed assumption on viral infection dynamics, J. Theor. Biol., 242 (2006), 464–477. https://doi.org/10.1016/j.jtbi.2006.03.014 doi: 10.1016/j.jtbi.2006.03.014
|
| [54] |
D. Wodarz, C. N. Chan, B. Trinite, N. L. Komarova, D. N. Levy, On the laws of virus spread through cell populations, J. Virol., 88 (2014), 13240–13248. https://doi.org/10.1128/JVI.02096-14. doi: 10.1128/JVI.02096-14
|
| [55] |
X. Tong, J. Chen, H. Miao, T. Li, L. Zhang, Development of an agent-based model (ABM) to simulate the immune system and integration of a regression method to estimate the key ABM parameters by fitting the experimental data, PLoS One, 10 (2015), e0141295. https://doi.org/10.1371/journal.pone.0141295 doi: 10.1371/journal.pone.0141295
|
| [56] |
J. Whitman, A. Dhanji, F. Hayot, S. C. Sealfon, C. Jayaprakash, Spatio-temporal dynamics of host-virus competition: A model study of influenza A, J. Theor. Biol., 484 (2020), 110026. https://doi.org/10.1016/j.jtbi.2019.110026 doi: 10.1016/j.jtbi.2019.110026
|
| [57] |
A. Goyal, J. M. Murray, Modelling the impact of cell-to-cell transmission in hepatitis B virus, PLoS One, 11 (2016), e0161978. https://doi.org/10.1371/journal.pone.0161978 doi: 10.1371/journal.pone.0161978
|
| [58] |
S. Wasik, P. Jackowiak, M. Figlerowicz, J. Blazewicz, Multi-agent model of hepatitis C virus infection, Artif. Intell. Med., 60 (2014), 123–131. https://doi.org/10.1016/j.artmed.2013.11.001 doi: 10.1016/j.artmed.2013.11.001
|
| [59] |
D. LaBarre, R. Lowy, Improvements in methods for calculating virus titer estimates from $TCID_50$ and plaque assays, J. Virol. Meth., 96 (2001), 107–126. https://doi.org/10.1016/S0166-0934(01)00316-0 doi: 10.1016/S0166-0934(01)00316-0
|
| [60] |
L. T. Pinilla, B. P. Holder, Y. Abed, G. Boivin, C. A. A. Beauchemin, The H275Y neuraminidase mutation of the pandemic A/H1N1 influenza virus lengthens the eclipse phase and reduces viral output of infected cells, potentially compromising fitness in ferrets, J. Virol., 86 (2012), 10651–10660. https://doi.org/10.1128/JVI.07244-11 doi: 10.1128/JVI.07244-11
|