Optimal control of a multiple delay model for hepatitis B virus infection of hepatocytes with DNA-containing capsids, immune response and drug therapy

Volume 35, Issue 2, pp 182--207 https://dx.doi.org/10.22436/jmcs.035.02.05

Publication Date: May 08, 2024 Submission Date: September 25, 2023 Revision Date: October 23, 2023 Accteptance Date: March 29, 2024

Download PDF

Download XML

• Export citations
  • CITATIONS & REFERENCES
  • EndNote (.ENW)
  • Mendeley, JabRef (.BIB)
  • Papers, Zotero, Reference Manager, RefWorks (.RIS)
  • ARTICLE CITATION
  • EndNote (.ENW)
  • Mendeley, JabRef (.BIB)
  • Papers, Zotero, Reference Manager, RefWorks (.RIS)
  • REFERENCES
  • EndNote (.ENW)
  • Mendeley, JabRef (.BIB)
  • Papers, Zotero, Reference Manager, RefWorks (.RIS)

• 33 Downloads
• 210 Views

Authors

P. Yosyingyong - Department of Mathematics, Faculty of Science, Naresuan University, Phitsanulok, 65000, Thailand. R. Viriyapong - Department of Mathematics, Faculty of Science, Naresuan University, Phitsanulok, 65000, Thailand. E. J. Moore - Department of Mathematics, Faculty of Applied Science, King Mongkut's University of Technology North Bangkok, Bangkok, 10800, Thailand.

Abstract

Since the number of patients with chronic hepatitis B virus (HBV) infection is still high globally, appropriate strategies for therapy control are required. In this study, we propose an HBV model consisting of six populations, namely uninfected hepatocytes, infected hepatocytes, HBV DNA-containing-capsids, free viruses, antibodies and cytotoxic T-lymphocytes. The model also includes an adaptive immune response, cure rate of infected hepatocytes, two time delays and two control therapies. The two time delays are delay of virus production after infection by HBV and delay of antigenic stimulation generating cytotoxic T-lymphocytes. The two optimal controls are a therapy for blocking new infection and a therapy for inhibiting viral production. The existence, uniqueness, non-negativity and boundedness of model solutions are first proved. The model is shown to have three equilibrium states, namely infection-free, immune-free and endemic. The basic reproduction number for the stability of the infection-free equilibrium is derived and a sensitivity analysis is carried out to determine the most important parameters to change to control the disease. Then optimal control and the Pontryagin maximum principle are used to maximize the concentrations of uninfected hepatocytes, antibodies and cytotoxic T-lymphocytes at minimum cost.

Share and Cite

Keywords

• Delay model
• drug therapy
• HBV-DNA containing capsids
• hepatitis B virus
• immune response
• cytotoxic T-lymphocyte (CTLs)
• Pontryagin maximum principle

MSC

•  00A71
•  34A34
•  49J15

References

• [1] R. Ahmed, D. Gray, Immunological memory and protective immunity: understanding their relation, Science, 272 (1996), 54–60
  • View Article
  • Google Scholar


• [2] A. Alberti, S. Diana, G. Sculard, A. Eddleston, R. Williams, Detection of a new antibody system reacting with dane particles in hepatitis B virus infection, Br. Med. J., 2 (1978), 1056–1058
  • View Article
  • Google Scholar


• [3] K. Allali, A. Meskaf, A. Tridane, Mathematical modeling of the adaptive immune responses in the early stage of the HBV infection, Int. J. Differ. Equ., 2018 (2018), 13 pages
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [4] H. Anton, C. Rorres, Elementary Linear Algebra, Wiley, New York (2019)
  • View Article
  • Google Scholar


• [5] A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model, J. Therm. Sci., 20 (2016), 763–769
  • View Article
  • Google Scholar


• [6] A. Atangana J. F. G´omez-Aguilar, A new derivative with normal distribution kernel; Theory, methods and applications, Phys. A, 476 (2017), 1–14
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [7] M. Bachraoui, K. Hattaf, N. Yousfi, Qualitative analysis of a fractional model for HBV infection with capsids and adaptive immunity, Int. J. Dyn. Syst. Differ. Equ., 12 (2022), 146–162
  • View Article
  • Google Scholar
  • MATH


• [8] G. Birkhoff, G.-C. Rota, Ordinary differential equations, John Wiley & Sons, New York (1989)
  • View Article
  • Google Scholar


• [9] V. Bruss, Envelopment of the hepatitis B virus nucleocapsid, Virus Res., 106 (2004), 199–209
  • View Article
  • Google Scholar


• [10] M. Caputo, Linear model of dissipation whose Q is almost frequency independent. II, Geophys. J. Int., 13 (1967), 529–539
  • View Article
  • Google Scholar


• [11] M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1 (2015), 73–85
  • View Article
  • Google Scholar


• [12] L. Cesari, Existence theorems for weak and usual optimal solutions in Lagrange problems with unilateral constraints. I, Trans. Amer. Math. Soc., 124 (1966), 369–412
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [13] F. F. Chenar, Y. Kyrychko, K. Blyuss, Mathematical model of immune response to hepatitis B, J. Theoret. Biol., 447 (2018), 98–110
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [14] F. V. Chisari, Cytotoxic T cells and viral hepatitis, J. Clin. Invest., 99 (1997), 1472–1477
  • View Article
  • Google Scholar


• [15] S. M. Ciupe, R. M. Ribeiro, A. S. Perelson, Antibody responses during hepatitis B viral infection, PLoS Comput. Biol., 10 (2014), 1–16
  • View Article
  • Google Scholar


• [16] J. Danane, K. Allali, Mathematical analysis and treatment for a delayed hepatitis B viral infection model with the adaptive immune response and DNA–containing capsids, High-throughput, 7 (2018), 35–50
  • View Article
  • Google Scholar


• [17] J. Danane, A. Meskaf, K. Allali, Optimal control of a delayed hepatitis B viral infection model with HBV DNA-containing capsids and CTL immune response, Optimal Control Appl. Methods, 39 (2018), 1262–1272
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [18] R. D. Driver, Ordinary and delay differential equations, Springer-Verlag, New York-Heidelberg (1977)
  • View Article
  • Google Scholar


• [19] S. Eikenberry, S. Hews, J. D. Nagy, Y. Kuang, The dynamics of a delay model of hepatitis B virus infection with logistic hepatocyte growth, Math. Biosci. Eng., 6 (2009), 283–299
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [20] E. M. Farah, S. Amine, K. Allali, Dynamics of a time-delayed two-strain epidemic model with general incidence rates, Chaos Solitons Fractals, 153 (2021), 15 pages
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [21] W. H. Fleming, R. W. Rishel, Deterministic and Stochastic Optimal Control, Springer-Verlag, New York (2012)
  • View Article
  • Google Scholar


• [22] J. E. Forde, S. M. Ciupe, A. Cintron-Arias, S. Lenhart, Optimal control of drug therapy in a hepatitis B model, Appl. Sci., 6 (2016), 219–236
  • View Article
  • Google Scholar


• [23] D. Ganem, A. M. Prince, Hepatitis B virus infection-natural history and clinical consequences, N. Engl. J. Med., 350 (2004), 1118–1129
  • View Article
  • Google Scholar


• [24] L. G¨ollmann, H. Maurer, Theory and applications of optimal control problems with multiple time-delays, J. Ind. Manag. Optim., 10 (2014), 413–441
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [25] S. A. Gourley, Y. Kuang, J. D. Nagy, Dynamics of a delay differential equation model of hepatitis B virus infection, J. Biol. Dyn., 2 (2008), 140–153
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [26] S. J. Hadziyannis, N. C. Tassopoulos, E. J. Heathcote, T.-T. Chang, G. Kitis, M. Rizzetto, P. Marcellin, S. G. Lim, Z. Goodman, M. S. Wulfsohn, Shelly Xiong, J. Fry, C. L. Brosgart, Adefovir dipivoxil for the treatment of hepatitis B e antigen-negative chronic hepatitis B, N. Engl. J. Med., 348 (2003), 800–807
  • View Article
  • Google Scholar


• [27] S. Hagiwara, N. Nishida, M. Kudo, Antiviral therapy for chronic hepatitis B: Combination of nucleoside analogs and interferon, World J. Hepatol., 7 (2015), 2427–2431
  • View Article
  • Google Scholar


• [28] S. Harroudi, A. Meskaf, K. Allali, Modelling the adaptive immune response in HBV infection model with HBV DNAcontaining capsids, Differ. Equ. Dyn. Syst., 31 (2020), 371–393
  • View Article
  • Google Scholar


• [29] K. Hattaf, Spatiotemporal dynamics of a generalized viral infection model with distributed delays and ctl immune response, Computation, 7 (2019), 1–16
  • View Article
  • Google Scholar


• [30] K. Hattaf, On the stability and numerical scheme of fractional differential equations with application to biology, Computation, 10 (2022), 1–12
  • View Article
  • Google Scholar


• [31] K. Hattaf, N. Yousfi, A generalized HBV model with diffusion and two delays, Comput. Math. Appl., 69 (2015), 31–40
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [32] S. Hews, S. Eikenberry, J. D. Nagy, Y. Kuang, Rich dynamics of a hepatitis B viral infection model with logistic hepatocyte growth, J. Math. Biol., 60 (2010), 573–590
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [33] R. Hilfer, Applications of Fractional Calculus in Physics, World Scientific Publishing Co., River Edge, NJ (2000)
  • View Article
  • MathSciNet
  • Google Scholar


• [34] F. Ibrahim, K. Hattaf, F. A. Rihan, S. Turek, Numerical method based on extended one-step schemes for optimal control problem with time-lags, Int. J. Dyn. Control, 5 (2017), 1172–1186
  • View Article
  • MathSciNet
  • Google Scholar


• [35] U. N. Katugampola, A new approach to generalized fractional derivatives, Bull. Math. Anal. Appl., 6 (2014), 1–15
  • View Article
  • Google Scholar


• [36] K. Hattaf, A new class of generalized fractal and fractal-fractional derivatives with non-singular kernels, Fractal Fract., 7 (2023), 1–16
  • View Article
  • Google Scholar


• [37] G. L. Kharatishvili, Maximum principle in the theory of optimum time-delay processes, Dokl. Akad. Nauk SSSR, 136 (1961), 39–42
  • View Article
  • Google Scholar


• [38] M. S. Khatun, M. H. A. Biswas, Optimal control strategies for preventing hepatitis B infection and reducing chronic liver cirrhosis incidence, Infect. Dis. Model., 5 (2020), 91–110
  • View Article
  • Google Scholar


• [39] A. A. Kilbas, Hadamard-type fractional calculus, J. Korean Math. Soc., 38 (2001), 1191–1204
  • View Article
  • Google Scholar
  • MATH


• [40] A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, New York (2006)
  • View Article
  • Google Scholar


• [41] T. Kisela, Fractional differential equations and their applications, Faculty of Mechanical Engineering, Institute of Mathematics, Brno (2018)
  • View Article
  • Google Scholar


• [42] S. Lenhart, J. T. Workman, Optimal Control Applied to Biological Models, Chapman & Hall/CRC, Boca Raton, FL (2007)
  • View Article
  • Google Scholar
  • MATH


• [43] S. R. Lewin, R. M. Ribeiro, T. Walters, G. K. Lau, S. Bowden, S. Locarnini, A. S. Perelson, Analysis of hepatitis B viral load decline under potent therapy: complex decay profiles observed, Hepatology, 34 (2001), 1012–1020
  • View Article
  • Google Scholar


• [44] D. G. Luenberger, Introduction to Dynamic Systems: Theory, Models & Applications, John Wiley & Sons, New York (1979)
  • View Article
  • Google Scholar


• [45] K. Manna, Dynamics of a delayed diffusive HBV infection model with capsids and CTL immune response, Int. J. Appl. Comput. Math., 4 (2018), 16 pages
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [46] K. Manna, S. P. Chakrabarty, Chronic hepatitis B infection and HBV DNA-containing capsids: Modeling and analysis, Commun. Nonlinear Sci. Numer. Simul., 22 (2015), 383–395
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [47] K. Manna, S. P. Chakrabarty, Global stability of one and two discrete delay models for chronic hepatitis B infection with HBV DNA-containing capsids, J. Comput. Appl. Math., 36 (2017), 525–536
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [48] K. Manna, K. Hattaf, Spatiotemporal dynamics of a generalized HBV infection model with capsids and adaptive immunity, Int. J. Appl. Comput. Math., 5 (2019), 29 pages
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [49] K. Manna, K. Hattaf, A generalized distributed delay model for hepatitis b virus infection with two modes of transmission and adaptive immunity. A mathematical study, Math. Methods Appl. Sci., 45 (2022), 11614–11634
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [50] A. Meskaf, Optimal control of a delayed hepatitis B viral infection infection model with DNA-containing capsids, the adaptive immune response and cure rate, Int. J. Open Probl. Comput. Sci. Math., 12 (2019), 18–33
  • View Article
  • Google Scholar


• [51] A. Meskaf, K. Allali, Y. Tabit, Optimal control of a delayed hepatitis B viral infection model with cytotoxic T-lymphocyte and antibody responses, Int. J. Dyn. Control, 5 (2017), 893–902
  • View Article
  • MathSciNet
  • Google Scholar


• [52] F. N. Ngoteya Y. Nkansah-Gyekye, Sensitivity analysis of parameters in a competition model, Appl. Comput. Math., 4 (2015), 363–368
  • View Article
  • Google Scholar


• [53] I. Petr´aˇ, Fractional-order Nonlinear Systems, Springer, Heidelberg (2011)
  • View Article
  • Google Scholar


• [54] I. Podlubny, Fractional differential equations, Academic Press, San Diego (1999)
  • View Article
  • Google Scholar


• [55] L. S. Pontryagin, The Mathematical Theory of Optimal Processes, Routledge, London (2018)
  • View Article
  • Google Scholar


• [56] Md. Samsuzzoha, M. Singh, D. Lucy, Uncertainty and sensitivity analysis of the basic reproduction number of a vaccinated epidemic model of influenza, Appl. Math. Model., 37 (2013), 903–915
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [57] D. Sun F. Liu, Analysis of a new delayed HBV model with exposed state and immune response to infected cells and viruses, BioMed Res. Int., 2017 (2017), 19 pages
  • View Article
  • Google Scholar


• [58] A. Tan, S. Koh, A. Bertoletti, Immune response in hepatitis B virus infection, Cold Spring Harb. Perspect. Med., (2015), 1–18
  • View Article
  • Google Scholar


• [59] P. van den Driessche, J.Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [60] WHO, Hepatitis B, , (2019),
  • View Article


• [61] J. Xu, G. Huang, Global stability and bifurcation annalysis of a virus infection model with nonlinear incidence and multiple delays, Fractal Fract., 7 (2023), 1–35
  • View Article
  • Google Scholar


• [62] P. Yosyingyong, R. Viriyapong, Global stability and optimal control for a hepatitis B virus infection model with immune response and drug therapy, J. Appl. Math. Comput., 60 (2019), 537–565
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [63] P. Yosyingyong, R. Viriyapong, Global dynamics of multiple delays within-host model for a hepatitis b virus infection of hepatocytes with immune response and drug therapy, Math. Biosci. Eng., 20 (2023), 7349–7386
  • View Article
  • MathSciNet
  • Google Scholar


• [64] N. Yousfi, K. Hattaf, A. Tridane, Modeling the adaptive immune response in HBV infection, J. Math. Biol., 63 (2011), 933–957
  • View Article
  • MathSciNet
  • Google Scholar
  • MATH


• [65] M.-F. Yuen, D.-S. Chen, G. M. Dusheiko, H. L. A. Janssen, D. T. Lau, S. A. Locarnini, M. G. Peters, C.-L. Lai, Hepatitis B virus infection, Nat. Rev. Dis. Primers, 4 (2018), 1–20
  • View Article
  • Google Scholar


• [66] X. Zhou, X. Song, X. Shi, A differential equation model of hiv infection of CD4+ T-cells with cure rate, J. Math. Anal. Appl., 342 (2008), 1342–1355
  • View Article
  • Google Scholar
  • MATH


• [67] K. Zhuang, Dynamical analysis of a delayed hepatitis b virus model with immune responses, WSEAS Trans. Math., 15 (2016), 325–334
  • View Article
  • Google Scholar