Journal of Applied Mathematics and Physics
Vol.07 No.07(2019), Article ID:93847,16 pages
10.4236/jamp.2019.77104
The Global Stability Analysis of a Mathematical Cellular Model of Hepatitis C Virus Infection with Non-Cytolytic Process
Alexis Nangue1, Cyprien Fokoue2, Raoue Poumeni1
1Department of Mathematics, Higher Teachers’ Training College, University of Maroua, Maroua, Cameroon
2Department of Mathematics and Computer Science, Faculty of Science, University of Maroua, Maroua, Cameroon
Copyright © 2019 by author(s) and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
Received: May 28, 2019; Accepted: July 21, 2019; Published: July 24, 2019
ABSTRACT
The aim of this work is to analyse the global dynamics of an extended mathematical model of Hepatitis C virus (HCV) infection in vivo with cellular proliferation, spontaneous cure and hepatocyte homeostasis. We firstly prove the existence of local and global solutions of the model and establish some properties of this solution as positivity and asymptotic behaviour. Secondly we show, by the construction of appropriate Lyapunov functions, that the uninfected equilibrium and the unique infected equilibrium of the mathematical model of HCV are globally asymptotically stable respectively when the threshold number and when .
Keywords:
HCV Model, Global Solutions, Non-Cytolytic Process, Invariant Set, Lyapunov Functions, Basic Reproduction Number, Equilibrium Points
1. Introduction
According to [1] [2] , approximately 200 million people worldwide are persistently infected with the hepatitis C virus (HCV) and are at risk of developing chronic liver disease, cirrhosis, and hepatocellular carcinoma. HCV infection therefore represents a significant global public health problem. HCV establishes chronic hepatitis in 60% - 80% of infected adults [3]. A vaccine against infection with HCV does not exist yet, and standard treatment with interferon-α and ribavirin has produced sustained virological response rates of approximately 50%, with no effective alternative treatment for nonresponders to this treatment protocol. A model of human immunodeficiency virus infection was adapted by Neumann et al. [4] to study the kinetics of chronic HCV infection during treatment and some mathematical analysis was done by [5]. Since then viral kinetics modeling has played an important role in the analysis of HCV RNA decay during antiviral therapy (see Perelson [6] , Perelson et al. [7] ). The original Neumann et al. model for HCV infection included three differential equations representing the populations of target cells, productively infected cells, and virus. In this paper we are going to study global dynamics of an HCV infection mathematical model with full logistic terms, antivirus treatments and homeostasis phenomenon. A similar work has been done by A. Nangue [8] concerning a mathematical intracellular HCV infection model with therapy.
1.1. The Compartmental Model
There are too many mathematical models of HCV dynamics amongst those, the original model or model of Newmann [4] and its extended models as those in [5] [9] for example. Each model can be represented by a compartmental scheme. A compartmental scheme is a scheme for estimating the variation in the number of individuals in each compartment over time. Figure 1 is the schematic representation of the extended model, which we will study, of HCV with cellular proliferation and spontaneous healing designed by T. C. Reluga et al. [10]. This model expands the viral dynamics of the original model of infection and the disappearance of HCV by incorporating the proliferation and death density dependence. In addition to cell proliferation, the number of uninfected hepatocytes may increase through immigration or differentiation of hepatocyte precursors that develop into hepatocytes at a constitutive rate of s or by spontaneous infected
Figure 1. Schematic representation of HCV infection models. T and I represent target and infected cells, respectively, and V represents free virus. The parameters shown in the figure are defined in the text. The original model of Neumann et al. [4] assumed that there is no proliferation of target and infected cells (i.e. ) and no spontaneous cure (i.e. ). The extended model of Dahari, Ribeiro, and Perelson [11] , which was used for predicting complex HCV kinetics under therapy, includes target and infected cell proliferation without cure ( , and ). A model including both proliferation and the spontaneous cure of infected cells (dashed line; ) was used to explain the kinetics of HCV in primary infection in chimpanzees [12].
hepatocyte healing by a non-cytolytic process at the rate q.
The model proposed by Dahari and coworkers [11] [13] expands on the standard HCV viral-dynamic model [4] of infection and clearance by incorporating density-dependent proliferation and death. Uninfected hepatocytes or noninfected hepatocytes, T, are infected at a rate per free virus per hepatocyte. Infected cells, I, produce free virus at rate p per cell but also die with rate . Free virus, V, is cleared at rate c by immune and other degradation processes. Besides infection processes, hepatocyte numbers are influenced by homeostatic processes. Uninfected hepatocytes die at rate . Both infected and uninfected hepatocytes proliferate logistically with maximum rates and , respectively, as long as the total number of hepatocytes is less than . Besides proliferation, uninfected hepatocytes may increase in number through immigration or differentiation of hepatocyte precursors that develop into hepatocytes at constitutive rate s, or by spontaneous cure of infected hepatocytes through a noncytolytic process at rate q. Treatment with antiviral drugs reduces the infection rate by a fraction and the viral production rate by a fraction . It should be noted that and are parameters which values are non-negative and less than one. A further comprehensive survey on the description of the model is given in [10] [14] [15]. Given the meanings of and , the term represents the mass action principle; is the rate of infection of healthy hepatocytes T by interaction with virus V.
1.2. The Mathematical Model
According to Reluga et al. [10] and more precisely according to the schematic representation of HCV infection model in Figure 1 we have the following dynamics:
· the variation of the healthy hepatocytes or uninfected hepatocytes, T, is expressed by the following equation:
(1)
· the variation of infected hepatocytes, I, is expressed by the following equation:
(2)
· the variation of free virions or virus, V, is expressed by the following equation:
(3)
Thus, the phenomenon described above is governed by the following mathematical model (4), which is a system of three differential autonomous equations:
(4)
To analyse the system (4) we need the following initial conditions:
(5)
For biological significance of the parameters, three assumptions are employed. (a) Due to the burden of supporting virus replication, infected cells may proliferate more slowly than uninfected cells, i.e. . (b) To have a physiologically realistic model, in an uninfected liver when is reached, liver size should no longer increase, i.e. . (c) Infected cells have a higher turnover rate than uninfected cells, i.e. . The interpretations and biologically plausible values of other parameters and a further comprehensive survey on the description of (4) is given in [10]. Besides HCV infection, the similar model of (4) is also used to describe the dynamics of HBV or HIV infection, in which the full logistic terms mean the proliferation of uninfected/infected hepatocytes [12] [16] [17] , or the mitotic transmission of uninfected/infected CD4+ T cells.
Our goal is therefore to analyze the stability of an extended model of HCV infection in a patient with cell proliferation and spontaneous healing given by (4) to reveal significant information on pathogenesis and dynamics of this virus. The paper is organized as follows: In Section 1, first focuses on some properties of the solutions of the model, then we calculate the basic reproduction ratio , which is an indispensable element in the study and analysis of the models. We theoretically analyze the local stability where we widely use the works of A. Nangue et al. [18] in Section 2. In Section 3 we theoretically analyze with some assumptions the global stability of the model by constructing appropriate Lyapounov’s functions.
2. Properties of Solutions to the Initial Value Problem (4), (5)
2.1. Positivity, Global Solutions and Asymptotic Behaviour
Theorem 1. Let . There exists and functions continuously differentiable such that is a solution of system (1) satisfying (4).
Theorem 2. Let be a solution of the system (1) over an interval such that et .
If are positive, then , and are also positive for all .
Proof. We are going to prove by contradiction. so suppose there is such that or or .
Let
Let also be the smallest of all t in the interval such that , , and for a certain i.
Then each of the equations of the system (4) can be written where is a non negative function and any function. As a consequence and , . A contradiction. ,
Theorem 3. [18] The solutions of the Cauchy problem (4), (5), with positive initial data, exist globally in time in the future that is on .
Theorem 4. For any positive solution of system (4), (5) we have:
where
Proof. Summing Equations (6) and (7), we get:
thus since .
Let , , , and let us solve the following equation
(6)
Coupled to Equation (6) the initial condition:
(7)
The solving of the problem (6), (7) gives for all ,
As for all , , it follows that:
Let
Therefore
Since T and I are positive and , so it follows that and .
Equation (3), according to Gromwall inequality, leads to:
where
This completes the proof of theorem 4. ,
2.2. Basic Reproduction Ratio , Invariant Set of the Model and Equilibria
Proposition 5. The uninfected equilibrium point of the system (4) is given by
where:
We use the method proposed in [19] [20] to compute the basic reproduction number .
Proposition 6. The expression of the basic reproduction number associated to the system (4) is given by:
(8)
where
Remark 1. denotes the overall effectiveness rate of the drug.
Remark 2. Henceforth, we will let and .
Theorem 7. Let and be a maximal solution of the Cauchy problem (1), (4) ( ). If and then the set:
where:
is a positively invariant set by system (4).
When it exists, the infected equilibrium point is given by: where , and are positive constants that we are going to determine.
Lemma 1. [18] exists if and only if
Lemma 2. [18] When it exists, is defined by:
where:
and
The combination of the lemma 1 and the lemma 2 leads to the following theorem:
Theorem 8. The model (4) admits a unique infected equilibrium if and only if , where
When the unique equilibrium is the uninfected equilibrium point or the infection-free steady state .
3. Local Stability Analyses
3.1. Case of the Uninfected Equilibrium Point or Infection-Free Steady State
Theorem 9. The infection-free steady state of model (4) is locally asymptotically stable if and unstable if .
Proof. See the appendice of [18]. ,
3.2. Case of Infected Equilibrium Point
We start this section by this lemma where the proof can be found in [18].
Lemma 3 The characteristic equation of the Jacobian matrix of the system (4) at is given by the following cubic equation:
where:
Proof. See [18]. ,
Now let:
According to lemma 3 combined with the Routh-Hurwitz criterion [21] , we have the following results where the proofs can be found in [18].
Theorem 10. For model (4), when is valid, the unique endemic equilibrium is locally asymptotically stable if and unstable if .
Especially, we have:
Corollary 1. The infected steady state during the therapy of the model (4) is locally asymptotically stable if and unstable if .
4. Global Stability Analyses
The global stability analysis of a dynamical system is usually a very complex problem. One of the most efficient methods to solve this problem is Lyapunov’s theory. To build the functions of Lyapunov we will follow the method proposed by A. Korobeinikov [22] [23] [24]. In the proofs of the results that follow, to simplify the writings, we can use differently or for the derivation with respect to time.
4.1. Case of Infection-Free Steady State
Theorem 11. The infection-free steady state of the model (4) is globally asymptotically stable if the basic reproduction number and unstable if .
Proof. Consider the Lyapunov function:
L is defined, continuous and positive definite for all , , . Also, the global minimum occurs at the infection free equilibrium . Further, function L, along the solutions of system (4), satisfies:
i.e.
yet
hence, Further collecting terms, we have:
Furthermore,
hence
Since and , we have and if and only if and simultaneously.
Therefore, the largest compact invariant subset of the set
is the singleton . By the Lasalle invariance principle [25] , the infection-free equilibrium is globally asymptotically stable if . We have seen previously that if , at least one of the eigenvalues of the Jacobian matrix evaluated at has a positive real part. Therefore, the infection-free equilibrium is unstable when . This completes the proof of the theorem.
Remark 3. The Lyapunov function defined in the proof of theorem 11 has been obtained following the general form giving by Korobonikov [22] [23] [24] for the dynamic virus fondamental model.
4.2. Case of Infected Equilibrium Point
We recall:
Remark 4. The infected equilibrium point satisfies:
(9)
(10)
(11)
Now we are stating and demonstrating one of the most important results of this work.
Theorem 12. Suppose that , and . Then the infected steady state during therapy of model (4) is globally asymptotically stable as soon as it exists.
Proof. Consider the Lyapunov function defined by:
Let us show that and if and only if , , simultaneously.
The time derivative of L along the trajectories of system (1) is:
Collecting terms, and canceling identical terms with opposite signs, yields:
(12)
Reporting equalities (9), (10) and (11) of the remark 4 into (12), we have:
Note that
and
Recall that:
furthermore,
hence,
By hypothesis, this leads to:
Yet
since the geometric mean is less than or equal to the arithmetic mean.
It should be noted that and holds if and only if take the steady states values . Therefore, By the Lasalle invariance principle [22] , the infected equilibrium point is globally asymptotically stable. This completes the proof of this theorem. ,
5. Concluding Remark
To understand the dynamics of HCV infection and its infectious processes, mathematical models are present as an important and unavoidable tool. Global stability analysis has been done, by the technique of Lyapunov, to the model of HCV infection with proliferation cell and spontaneous healing, for revealing significant information for making good decision for the fighting against hepatitis C. This work is a starting point to many interesting other future investigations. We plan to extend our study by focusing on more realistic models such as: 1) mathematical models with delay which involve delay ordinary differential equations. 2) mathematical models taking into account space which involve Partial differential equations. 3) mathematical models taking into account random phenomena which evolve stochastic differential equations. We also plan to focus on others methods of studying global stability like the geometric method that can provide results with fewer hypotheses on mathematical model (4).
Acknowledgements
We thank the Editor and the referee for their comments. We are grateful to Professor Alan Rendall for valuable and tremendous discussions about this paper. We wish to thank him for introducing us to Mathematical Biology and to its relationship with Mathematical Analysis. We also thank the Higher Teachers’ Training College of the University of Maroua were this paper were initiated.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Cite this paper
Nangue, A., Fokoue, C. and Poumeni, R. (2019) The Global Stability Analysis of a Mathematical Cellular Model of Hepatitis C Virus Infection with Non-Cytolytic Process. Journal of Applied Mathematics and Physics, 7, 1531-1546. https://doi.org/10.4236/jamp.2019.77104
References
- 1. NIH (2002) National Institutes of Health Consensus Development Conference. Management of Hepatitis C: 2002-June 10-12, 2002. Hepatology, 36, S3-S20.https://doi.org/10.1053/jhep.2002.37117
- 2. WHO (2017) Global Hepatitis Repport. http://www.who.int/mediacentre/factsheets/fs164/en/
- 3. Seeff, L.B. (2002) Natural History of Chronic Hepatitis C. Hepatology, 36, S35-S46.https://doi.org/10.1053/jhep.2002.36806
- 4. Neumann, A.U., Lam, N.P., Dahari, H., Gretch, D.R., Wiley, T.E., Layden, T.J. and Perelson, A.S. (1998) Hepatitis C Viral Dynamics in Vivo and the Antiviral Efficacy of Interferon-α Therapy. Science, 282, 103-107.https://doi.org/10.1126/science.282.5386.103
- 5. Chong, M.S.F., Crossley, L.S.M. and Madzvamuse, A. (2015) The Stability Analyses of the Mathematical Models of Hepatitis C Virus Infection. Modern Applied Science, 9, 250-271. http://sro.sussex.ac.uk/52312/ https://doi.org/10.5539/mas.v9n3p250
- 6. Perelson, A., Neumann, A., Markowitz, M., Leonard, J. and Ho, D. (1996) HIV-1 Dynamics in Vivo: Virion Clearance Rate, Infected Cell Life-Span, and Viral Generation Time. Science, 271, 1582-1586. https://doi.org/10.1126/science.271.5255.1582
- 7. Perelson, A.S. (2002) Modelling Viral and Immune System Dynamics. Nature Reviews Immunology, 2, 28-36. https://doi.org/10.1038/nri700
- 8. Nangue, A. (2019) Global Stabilty Analysis of a Cellular Model of Hepatitis C Virus Infection under Treatement. arXiv:1901.09984v2.
- 9. Banerjee, S., Keval, R., Gakkhar, R. (2011) Modeling the Dynamics of Hepatitis C Virus with Combined Antiviral Drug Therapy: Interferon and Ribavirin. Bulletin of Mathematical Biology, arXiv:1105.3669v1.
- 10. Reluga, T.C. Dahari, H. and Perelson, A.S. (2009) Analysis of Hepatitis C Virus Infection Models with Hepatocyte Homeostasis. SIAM Journal on Applied Mathematics, 69, 999-1023. http://www.siam.org/journals/ojsa.php https://doi.org/10.1137/080714579
- 11. Dahari, H., Ribeiro, R.M. and Perelson, A.S. (2007) Triphasic Decline of HCV RNA during Antiviral Therapy. Hepatology, 46, 16-21. https://doi.org/10.1002/hep.21657
- 12. Dahari, H., Shudo, H., Ribeiro, R.M. and Perelson, A.S. (2009) Modeling Complex Decay Profiles of Hepatitis B Virus during Antiviral Therapy. Hepatology, 49, 32-38.
- 13. Daharia, H., Loa, A., Ribeiroa, R.M. and Perelson, A.S. (2007) Modeling Hepatitis C Virus Dynamics: Liver Regeneration and Critical Drug Efficacy. Journal of Theoretical Biology, 247, 371-381. http://www.sciencedirect.com https://doi.org/10.1016/j.jtbi.2007.03.006
- 14. Dahari, H., Layden-Almer, J.E., Kallwitz, E., Ribeiro, R.M., Cotler, S.J., Layden, T.J. and Perelson, A.S. (2009) A Mathematical Model of Hepatitis C Virus Dynamics in Patients with High Baseline Viral Loads or Advanced Liver Disease. Gastroenterology, 136, 1402-1409. https://doi.org/10.1053/j.gastro.2008.12.060
- 15. Dahari, H., Major, M., Zhang, X., Mihalik, K., Rice, C.M., Perelson, A.S., Feinstone, S.M. and Neumann, A.U. (2005) Mathematical Modeling of Primary Hepatitis C Infection: Noncytolytic Clearance and Early Blockage of Virion Production. Gastroenterology, 128, 1056-1066. https://doi.org/10.1053/j.gastro.2005.01.049
- 16. Ciupe, S.M. Ribeiro, S.M. Nelson, P.W. and Perelson, A.S. (2007) Modeling the Mechanisms of Acute Hepatitis B Virus Infection. Journal of Theoretical Biology, 247, 23-35. https://doi.org/10.1016/j.jtbi.2007.02.017
- 17. Li, Q., Lu, F., Deng, G. and Wang, K. (2014) Modeling the Effects of Covalently Closed Circular DNA and Dendritic Cells in Chronic HBV Infection. Journal of Theoretical Biology, 357, 1-9. https://doi.org/10.1016/j.jtbi.2014.04.037
- 18. Nangue, A., Fokoue, C. and Poumeni, R. (2019) Mathematical Analysis of an Extended Cellular Model of the Hepatitis C Virus Infection with Non-Cytolytic Process. ArXiv: 1905.02561v1.
- 19. Diekmann, O., Heesterbeek, J.P.A. and Metz, J.A.J. (1990) On the Definition and the Computation of the Basic Reproduction Ratio R0 in Models for Infectious Diseases in Heterogeneous Populations. Journal of Mathematical Biology, 28, 365-382.
- 20. Van den Driessche, P. and Watmough, J. (2002) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48. https://doi.org/10.1016/S0025-5564(02)00108-6
- 21. Murray, J.M. (2001) Mathematical Biology I: An Introduction. Springer, Berlin, Heidelberg.
- 22. Korobeinikov, A. (2004) Global Properties of Basic Virus Dynamics Models. Bulletin of Mathematical Biology, 66, 879-883. https://doi.org/10.1016/j.bulm.2004.02.001
- 23. Korobeinikov, A. (2004) Lyapunov Functions and Global Properties for SEIR and SEIS Epidemic Models. Mathematical Medicine and Biology, 21, 75-83.https://doi.org/10.1093/imammb21.2.75
- 24. Korobeinikov, A. and Wake, G.C. (2002) Lyapunov Functions and Global Stability for SIR, SIRS, and SIS Epidemiological Models. Applied Mathematics Letters, 15, 955-960. https://doi.org/10.1016/S0893-9659(02)00069-1
- 25. Khalil, H. (2002) Nonlinear Systems. 3rd Edition, Prentice Hall, New York.