Mathematical Analysis of a Transformed ODE from a PDE Multiscale Model of Hepatitis C Virus Infection
Tóm tắt
Mathematical modeling has revealed the quantitative dynamics during the process of viral infection and evolved into an important tool in modern virology. Coupled with analyses of clinical and experimental data, the widely used basic model of viral dynamics described by ordinary differential equations (ODEs) has been well parameterized. In recent years, age-structured models, called “multiscale model,” formulated by partial differential equations (PDEs) have also been developed and become useful tools for more detailed data analysis. However, in general, PDE models are considerably more difficult to subject to mathematical and numerical analyses. In our recently reported study, we successfully derived a mathematically identical ODE model from a PDE model, which helps to overcome the limitations of the PDE model with regard to clinical data analysis. Here, we derive the basic reproduction number from the identical ODE model and provide insight into the global stability of all possible steady states of the ODE model.
Tài liệu tham khảo
Althaus CL, De Vos AS, De Boer RJ (2009) Reassessing the human immunodeficiency virus type 1 life cycle through age-structured modeling: life span of infected cells, viral generation time, and basic reproductive number, R0. J Virol 83:7659–7667. https://doi.org/10.1128/jvi.01799-08
Dixit NM, Layden-Almer JE, Layden TJ, Perelson AS (2004) Modelling how ribavirin improves interferon response rates in hepatitis C virus infection. Nature 432:922–924. https://doi.org/10.1038/nature03153
Frohlich F, Kaltenbacher B, Theis FJ, Hasenauer J (2017) Scalable parameter estimation for genome-scale biochemical reaction networks. PLoS Comput Biol 13:e1005331. https://doi.org/10.1371/journal.pcbi.1005331
Guedj J, Neumann AU (2010) Understanding hepatitis C viral dynamics with direct-acting antiviral agents due to the interplay between intracellular replication and cellular infection dynamics. J Theor Biol 267:330–340. https://doi.org/10.1016/j.jtbi.2010.08.036
Guedj J et al. (2013) Modeling shows that the NS5A inhibitor daclatasvir has two modes of action and yields a shorter estimate of the hepatitis C virus half-life. In: Proceedings of the National Academy of Sciences of the United States of America, vol 110, pp 3991–3996. https://doi.org/10.1073/pnas.1203110110
Ho DD, Neumann AU, Perelson AS, Chen W, Leonard JM, Markowitz M (1995) Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature 373:123–126. https://doi.org/10.1038/373123a0
Huang G, Liu X, Takeuchi Y (2012) Lyapunov functions and global stability for age-structured HIV infection model. SIAM J Appl Math 72:25–38
Inaba H (2012) On a new perspective of the basic reproduction number in heterogeneous environments. J Math Biol 65:309–348. https://doi.org/10.1007/s00285-011-0463-z
Iwami S, Takeuchi Y, Liu X (2007) Avian-human influenza epidemic model. Math Biosci 207:1–25. https://doi.org/10.1016/j.mbs.2006.08.001
Iwami S, Takeuchi Y, Liu X, Nakaoka S (2009) A geographical spread of vaccine-resistance in avian influenza epidemics. J Theor Biol 259:219–228. https://doi.org/10.1016/j.jtbi.2009.03.040
Iwami S, Sato K, De Boer RJ, Aihara K, Miura T, Koyanagi Y (2012) Identifying viral parameters from in vitro cell cultures. Front Microbiol 3:319. https://doi.org/10.3389/fmicb.2012.00319
Iwasa Y, Andreasen V, Levin S (1987) Aggregation in model ecosystems. I. Perfect aggregation. Ecol Modell 37:287–302
Iwasa Y, Levin SA, Andreasen V (1989) Aggregation in model ecosystems II. Approximate aggregation mathematical medicine and biology. J IMA 6:1–23
Kitagawa K, Nakaoka S, Asai Y, Watashi K, Iwami S (2018) A PDE multiscale model of hepatitis C virus infection can be transformed to a system of ODEs. J Theor Biol 448:80–85. https://doi.org/10.1016/j.jtbi.2018.04.006
Magal P, McCluskey C, Webb G (2010) Lyapunov functional and global asymptotic stability for an infection-age model. Appl Anal 89:1109–1140
Martyushev A, Nakaoka S, Sato K, Noda T, Iwami S (2016) Modelling Ebola virus dynamics: implications for therapy. Antivir Res 135:62–73. https://doi.org/10.1016/j.antiviral.2016.10.004
McCluskey CC (2012) Global stability for an SEI epidemiological model with continuous age-structure in the exposed and infectious classes. Math Biosci Eng MBE 9:819–841. https://doi.org/10.3934/mbe.2012.9.819
Melnik AV, Korobeinikov A (2013) Lyapunov functions and global stability for SIR and SEIR models with age-dependent susceptibility. Math Biosci Eng MBE 10:369–378. https://doi.org/10.3934/mbe.2013.10.369
Nelson PW, Gilchrist MA, Coombs D, Hyman JM, Perelson AS (2004) An age-structured model of HIV infection that allows for variations in the production rate of viral particles and the death rate of productively infected cells. Math Biosci Eng MBE 1:267–288
Neumann AU, Lam NP, Dahari H, Gretch DR, Wiley TE, Layden TJ, Perelson AS (1998) Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferon-alpha therapy. Science 282:103–107
Nowak MA, May RM (2000) Virus dynamics. Oxford University Press, Oxford
Peifer M, Timmer J (2007) Parameter estimation in ordinary differential equations for biochemical processes using the method of multiple shooting. IET Syst Biol 1:78–88
Perelson AS (2002) Modelling viral and immune system dynamics. Nat Rev Immunol 2:28–36. https://doi.org/10.1038/nri700
Perelson AS, Nelson PW (1999) Mathematical analysis of HIV-1 dynamics in vivo. SIAM Rev 41:3–44
Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD (1996) HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science 271:1582–1586
Perelson AS et al (1997) Decay characteristics of HIV-1-infected compartments during combination therapy. Nature 387:188–191. https://doi.org/10.1038/387188a0
Pinilla LT, Holder BP, Abed Y, Boivin G, Beauchemin CA (2012) 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:10651–10660. https://doi.org/10.1128/jvi.07244-11
Quintela BM, Conway JM, Hyman JM, Guedj J, Dos Santos RW, Lobosco M, Perelson AS (2018) A new age-structured multiscale model of the Hepatitis C virus life-cycle during infection and therapy with direct-acting antiviral agents. Front Microbiol 9:601. https://doi.org/10.3389/fmicb.2018.00601
Ramratnam B et al (1999) Rapid production and clearance of HIV-1 and hepatitis C virus assessed by large volume plasma apheresis. Lancet 354:1782–1785. https://doi.org/10.1016/s0140-6736(99)02035-8
Reinharz V, Dahari H, Barash D (2018) Numerical schemes for solving and optimizing multiscale models with age of hepatitis C virus dynamics. Math Biosci 300:1–13. https://doi.org/10.1016/j.mbs.2018.03.011
Rong L, Perelson AS (2013) Mathematical analysis of multiscale models for hepatitis C virus dynamics under therapy with direct-acting antiviral agents. Math Biosci 245:22–30. https://doi.org/10.1016/j.mbs.2013.04.012
Rong L, Feng Z, Perelson AS (2007) Mathematical analysis of age-structured HIV-1 dynamics with combination antiretroviral therapy. SIAM J Appl Math 67:731–756
Rong L, Guedj J, Dahari H, Coffield DJ Jr, Levi M, Smith P, Perelson AS (2013) Analysis of hepatitis C virus decline during treatment with the protease inhibitor danoprevir using a multiscale model. PLoS Comput Biol 9:e1002959. https://doi.org/10.1371/journal.pcbi.1002959
Wang J, Zhang R, Kuniya T (2015) Global dynamics for a class of age-infection HIV models with nonlinear infection rate. J Math Anal Appl 432:289–313
Wei X et al (1995) Viral dynamics in human immunodeficiency virus type 1 infection. Nature 373:117–122. https://doi.org/10.1038/373117a0