Research Article  Open Access
Stability of an HIV/AIDS Treatment Model with Different Stages
Abstract
An HIV/AIDS treatment model with different stages is proposed in this paper. The stage of the HIV infection is divided into two stages, that is, HIVpositive in the asymptomatic stage of HIV infection and HIVpositive individuals in the preAIDS stage. The fact that some individuals with HIVpositive individuals after the treatment can be transformed into the compartment of HIVpositive individuals in the asymptomatic stage of HIV infection, the compartment of HIVpositive individuals in the preAIDS stage, or the compartment of individuals with fullblown AIDS is also considered. Mathematical analyses establish the idea that the global dynamics of the HIV/AIDS model are determined by the basic reproduction number . The diseasefree equilibrium is globally asymptotically stable if . The endemic equilibrium is globally asymptotically stable if for a special case. Numerical simulations are also conducted to support the analytic results.
1. Introduction
Human immunodeficiency virus (HIV) is one of the major life threatening viruses that are spreading worldwide. AIDS is caused by the human immunodeficiency virus (HIV), which has developed into a global pandemic since the first patient was identified in 1981, making it one of the most destructive epidemics in history [1, 2]. The HIV prevalence among pregnant women attending antenatal clinics rose from 0.8% in 1990 to 30.2% in 2005 and declined slightly to 29.1% in 2006 [3, 4].
Because of the fact that there is no vaccine, there are many obstacles in the AIDS treatment. Recently, the most prevalent treatment strategy for HIV infected patients is highly active antiretroviral therapies (HAART), which can prolong the lifespans and improve their life quality of infected individuals [5]. Current antiretroviral HIV/AIDS therapy involves the simultaneous administration of two or more antiviral drugs, typically chosen from two major classes: reverse transcriptase inhibitors (RTIs) and protease inhibitors (PIs) [6]. With the help of the antiretroviral therapy, patients can prolong the lifespan during which they are free of HIVrelated symptoms and therefore improve their life quality. At least 3.0 millions of people have been saved in the United States, which highlights the significant advances made in HIV disease treatment [7]. That is the reason why the early treatment of AIDS is very meaningful and of great urgency.
It is well known that mathematical modeling is a very important approach to understand the dynamics of any epidemic and to further develop various control and prevention policies [8, 9]. For the sake of the importance in biological meanings, dynamical properties of HIV/AIDS models and general theory on such dynamical systems have been studied by many authors in recent years [10–13] and references cited therein. They propose mathematical models that described the interactions of the HIV virus, CD4 cells, and CTLs within host, and these mathematical models can provide some insights into the dynamics of HIV viral load in vivo and may play a significant role in the development of a better understanding of HIV/AIDs and drug therapies.
Huo and Chen [14] introduce an HIV/AIDS epidemic model with antiretroviral treatment and establish the idea that the global dynamics of the HIV/AIDS model are determined by the basic reproduction number . They show that the early treatment of AIDS is necessary and meaningful. But they do not consider the compartments of HIVpositive individuals in the asymptomatic stage of HIV infection and the compartment of HIVpositive individuals in the preAIDS stage not receiving antiretroviral treatment. Yusuf and Benyah [15] present a deterministic model for controlling the spread of the disease using change in sexual habits and antiretroviral (ARV) therapy as control measures and formulate a fixed time optimal control problem subject to the model dynamics with the goal of finding the optimal combination of the two control measures that will minimize the cost of the control efforts as well as the incidence of the disease. But they neglect the fact that some individuals with HIVpositive after the treatment can be transformed into the compartment of HIVpositive individuals in the asymptomatic stage of HIV infection, the compartment of HIVpositive individuals in the preAIDS stage, or the compartment of individuals with fullblown AIDS.
Motivated by the above, in this paper, we will divide the stage of the HIV infection into two stages, that is, HIVpositive individuals in the asymptomatic stage of HIV infection and HIVpositive individuals in the preAIDS stage. We also consider the fact that some individuals with HIVpositive after the treatment can be transformed into the compartment of HIVpositive individuals in the asymptomatic stage of HIV infection, the compartment of HIVpositive individuals in the preAIDS stage, or the compartment of individuals with fullblown AIDS.
The organization of this paper is as follows. In the next section, we propose an HIV/AIDS model. In Section 3, the existence and the global stability of equilibria are investigated. Some numerical simulations are given in Section 4. Some discussions are given in last section.
2. The Model
2.1. System Description
The total population is divided into six compartments: , , , , , and . represents the number of susceptible individuals; represents the number of HIVpositive individuals in the asymptomatic stage of HIV infection; represents the number of HIVpositive individuals in the preAIDS stage not receiving treatment; represents the number of individuals with fullblown AIDS not receiving treatment; represents the number of individuals who are receiving ARV treatment; represents the number of individuals who change their sexual habits and maintain the habits for the rest of their lives. The total number of population at time is given by The model structure is shown in Figure 1. The transfer diagram leads to the following system of ordinary differential equations: For simplicity, denote , , , , , and , and then system 2 becomes where is recruitment rate, is the natural death rate, is the contact rate between the susceptible and the infected individuals. is the contact rate between the susceptible and the preAIDS stages. is the rate at which individuals leave the infection class and become individuals in the preAIDS stage. is the rate at which individuals leave the preAIDS class and become individuals with fullblown AIDS. is the rate at which an individual in the asymptomatic stage of HIV infection receives treatment. is the rate at which an individual in the preAIDS stage receives treatment. and are the diseaseinduced death rates for individuals in compartments and . is the rate at which susceptible individuals change their sexual habits per unit time. is the rate at which a treated individual leaves compartment . is the rate at which treatment of an individual with HIV fails. is the rate at which an individual in the preAIDS stage is treated successfully. is the rate at which tratment of an individual in the preAIDS stage fails. is the rate at which treatment of an individual in the preAIDS stage fails. All the parameters can be found in Table 1.

2.2. Basic Properties
2.2.1. Invariant Region
It is necessary to prove that all solutions of system 2 with positive initial data will remain positive for all times . This will be established by the following lemma.
Lemma 1. The feasible region defined by with initial condition , , , , , and , is positive invariant for system 2.
Proof. Adding the equations of system 2, we obtain It follows that where represents the initial values of the total population. Thus . It implies that the region is a positively invariant set for system 2.
2.2.2. Positivity of Solutions
Lemma 2. If , , , , , and , the solutions , , , , , and of system 3 are positive for all .
Proof. Under the given initial conditions, we need to prove that the solutions of system 3 are positive [16]; if not, we assume that there is a contradiction: there exists a first time such that
there exists a such that
there exists a such that
there exists a such that
there exists a such that
and there exists a such that
In the first case, we have
which is a contradiction, so that we can derive that , . In the second case, we have
which is a contradiction, so that we can derive that , . Similarly, it can be shown that , , , and for all .
Thus, the solutions , , , , , and of system 3 remain positive for all .
3. Analysis of the Model
There exist one diseasefree equilibrium and one endemic equilibrium for system 3.
3.1. DiseaseFree Equilibrium and the Reproduction Number
The model has a diseasefree equilibrium given by In the following, the basic reproduction number of system 3 will be obtained by the next generation matrix method formulated in [17, 18].
Letting , then system 3 can be written as where The Jacobian matrices of and at the free equilibrium are, respectively, where The model reproduction number, denoted by , is thus given by
3.2. Global Stability of
Theorem 3. The diseasefree equilibrium is globally asymptotically stable for and unstable otherwise.
Proof. According to [19, 20], we introduce the Lyapunov candidate function Differentiating with respect to time yields Substituting the model equations into 22, we get With and , we have In order to get the value of and , we can use From 25, we can calculate Taking and into 24, we can get Since we can obtain when , While , only . Substituting into system 3, we can have , , and as . According to LaSalle’s invariance principle [21], the diseasefree equilibrium is globally asymptotically stable.
3.3. Endemic Equilibrium
3.3.1. Existence of the Endemic Equilibrium
If , system 3 has a unique endemic equilibrium , where
3.3.2. Global Stability of the Endemic Equilibrium
Theorem 4. If , , the endemic equilibrium is globally asymptotically stable.
Proof. To study the globally asymptotic stability of the endemic equilibrium, motivated by [22, 23], we use a Lyapunov function as follows: Applying the replacements the derivative of is given by Since satisfies the equation of the model with system 3, We choose the coefficients of , , , , and which are equal to zero; that is, From 36, we have From 37, the derivative of is given by According to [24], since the arithmetical mean is greater than or equal to the geometrical mean, then for and if and only if ; for , , and if and only if , and if and only if . Futhermore, if and only if . The maximum invariant set of system 3 on the set is the singleton . Thus, for system 3, the endemic equilibrium is globally asymptotically stable if by LaSalle invariance principle.
Remark 5. Biologically speaking, implies that the HIVpositive individuals in the asymptomatic stage of HIV infection do not receive treatment. Theorem 4 tells us that when HIVpositive individuals in the asymptomatic stage of HIV infection do not receive treatment, the endemic equilibrium is globally asymptotically stable. This means that AIDS is persistent. When , the same result is validated by numerical simulations (see Figure 4).
4. Numerical Simulation
In this section, we present some numerical simulations of system 3 to illustrate our results. To illustrate the analytic results obtained above, we give some simulations using the parameter values in Table 1.
First, when , we choose , , numerical simulation gives , and the diseasefree equilibrium is globally asymptotically stable (see Figure 2). Figure 2 indicates the effect of treatment on stability of the equilibria of the HIV/AIDS infection model. At first, The number of patients who do not receive treatment is increasing constantly. After treatment, The number of patients is decreasing and tends to be zero at the end, which indicates that the spread of infection can be slowed down.
Second, we choose , , numerical simulation gives , and the diseasefree equilibrium is globally asymptotically stable (Figure 3). If we give a proper treatment, the number of individuals with AIDS/HIV will be controlled to some extent. Hence, the disease will be eradicated in the future.
Third, we choose , , numerical simulation gives , and the endemic equilibrium is globally asymptotically stable (Figure 4). Figure 4 indicates that when the basic reproduction is greater than one, as time goes on, the number of people finally tends to be constant, which means that the endemic equilibrium is globally asymptotically stable.
Finally, we show the relation between , , and . When the value of and is increasing at the time, the basic reproduction is reducing sharply. What is more, Figure 5 confirms that, with the increase in the rate of treatment, the number of individuals with HIVpositive and fullblown AIDS but not receiving treatment and who are receiving ARV treatment tends to zero, which indicates that treatment is important and can prolong their lives.
5. Discussion
A mathematical model is proposed and analyzed to study the spread of HIV/AIDS with treatment. In this paper we have analyzed a stage structured model for HIV and the effect of treatment has also been studied. In the presented model, we get two equilibria: the diseasefree equilibrium and the endemic equilibrium. We further consider global asymptotic stability of the diseasefree equilibrium by using the wellknown LyapunovLaSalle invariance principal. It is found that the diseasefree equilibrium is globally asymptotically stable when the basic reproduction number is less than one. When the basic reproduction number is greater than one and , the endemic equilibrium is globally asymptotically stable. Our results show that early treatment for individuals in asymptomatic stage of HIV infection or the preAIDS stage is very important.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was partially supported by the NNSF of China (11461041), the NSF of Gansu Province of China (2013GS09485, 1107RJZA088), the NSF for Distinguished Young Scholars of Gansu Province of China (1111RJDA003), the Special Fund for the Basic Requirements in the Research of University of Gansu Province of China, and the Development Program for HongLiu Distinguished Young Scholars in Lanzhou University of Technology.
References
 World Health Organization, HIV/AIDS, http://www.who.int/topics/hiv_aids/en/.
 N. Kaur, M. Ghosh, and S. S. Bhatia, “Modeling the spread of HIV in a stage structured population: effect of awareness,” International Journal of Biomathematics, vol. 5, no. 5, Article ID 1250040, 2012. View at: Publisher Site  Google Scholar
 HIV and AIDS in South Africa, http://www.avert.org/aidssouthafrica.htm.
 H. F. Huo and L. X. Feng, “Global stability for an HIV/AIDS epidemic model with different latent stages and treatment,” Applied Mathematical Modelling, vol. 37, no. 3, pp. 1480–1489, 2012. View at: Google Scholar
 A. S. Perelson and P. W. Nelson, “Mathematical analysis of HIV1 dynamics in vivo,” SIAM Review, vol. 41, no. 1, pp. 3–44, 1999. View at: Publisher Site  Google Scholar
 J. Lou, Y. Lou, and J. Wu, “Threshold virus dynamics with impulsive antiretroviral drug effects,” Journal of Mathematical Biology, vol. 65, no. 4, pp. 623–652, 2012. View at: Publisher Site  Google Scholar
 R. P. Walensky, A. D. Paltiel, E. Losina et al., “The survival benefits of AIDS treatment in the United States,” The Journal of Infectious Diseases, vol. 194, no. 1, pp. 11–19, 2006. View at: Publisher Site  Google Scholar
 L. Cai, X. Li, M. Ghosh, and B. Guo, “Stability analysis of an HIV/AIDS epidemic model with treatment,” Journal of Computational and Applied Mathematics, vol. 229, no. 1, pp. 313–323, 2009. View at: Publisher Site  Google Scholar
 X. Lin, H. W. Hethcote, and P. van den Driessche, “An epidemiological model for HIV/AIDS with proportional recruitment,” Mathematical Biosciences, vol. 118, no. 2, pp. 181–195, 1993. View at: Publisher Site  Google Scholar
 D. Li and W. Ma, “Asymptotic properties of a HIV1 infection model with time delay,” Journal of Mathematical Analysis and Applications, vol. 335, no. 1, pp. 683–691, 2007. View at: Publisher Site  Google Scholar
 X. Tian and R. Xu, “Global stability and Hopf bifurcation of an HIV1 infection model with saturation incidence and delayed CTL immune response,” Applied Mathematics and Computation, vol. 237, pp. 146–154, 2014. View at: Publisher Site  Google Scholar
 S. Bonhoeffer, J. M. Coffin, and M. A. Nowak, “Human immunodeficiency virus drug therapy and virus load,” Journal of Virology, vol. 71, no. 4, pp. 3275–3278, 1997. View at: Google Scholar
 Z. Yuan and X. Zou, “Global threshold dynamics in an hiv virus model with nonlinear infection rate and distributed invasion and production delays,” Mathematical Biosciences and Engineering, vol. 10, no. 2, pp. 483–498, 2013. View at: Publisher Site  Google Scholar
 H. F. Huo and R. Chen, “Modelling the effect of treatment on an HIV/AIDS epidemic model,” submitted. View at: Google Scholar
 T. T. Yusuf and F. Benyah, “Optimal strategy for controlling the spread of HIV/AIDS disease: a case study of South Africa,” Journal of Biological Dynamics, vol. 6, no. 2, pp. 475–494, 2012. View at: Publisher Site  Google Scholar
 H.F. Huo, S.J. Dang, and Y.N. Li, “Stability of a twostrain tuberculosis model with general contact rate,” Abstract and Applied Analysis, vol. 2010, Article ID 293747, 31 pages, 2010. View at: Publisher Site  Google Scholar
 P. van den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at: Publisher Site  Google Scholar
 M. BaniYaghoub, R. Gautam, Z. Shuai, P. van den Driessche, and R. Ivanek, “Reproduction numbers for infections with freeliving pathogens growing in the environment,” Journal of Biological Dynamics, vol. 6, no. 2, pp. 923–940, 2012. View at: Publisher Site  Google Scholar
 G. H. Li and Z. Jin, “Global stability of a SEIR epidemic model with infectious force in latent, infected and immune period,” Chaos, Solitons and Fractals, vol. 25, no. 5, pp. 1177–1184, 2005. View at: Publisher Site  Google Scholar
 Y. Cheng, J. Wang, and X. Yang, “On the global stability of a generalized cholera epidemiological model,” Journal of Biological Dynamics, vol. 6, no. 2, pp. 1088–1104, 2012. View at: Publisher Site  Google Scholar
 J. P. LaSalle, The Stability of Dynamical Systems, vol. 25 of Regional Conference Series in Applied Mathmatics, SIAM, 1976.
 R. Xu, “Global stability of an HIV1 infection model with saturation infection and intracellular delay,” Journal of Mathematical Analysis and Applications, vol. 375, no. 1, pp. 75–81, 2011. View at: Publisher Site  Google Scholar
 J. Q. Li, Y. L. Yang, and Y. C. Zhou, “Global stability of an epidemic model with latent stage and vaccination,” Nonlinear Analysis: Real World Applications, vol. 12, no. 4, pp. 2163–2173, 2011. View at: Publisher Site  Google Scholar
 J. Q. Li, Z. A. Ma, and F. Q. Zhang, “Stability analysis for an epidemic model with stage structure,” Nonlinear Analysis: Real World Applications, vol. 9, no. 4, pp. 1672–1679, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 HaiFeng Huo and Rui Chen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.