Research Article  Open Access
Nigar Ali, Gul Zaman, Abdullah, Aisha M. Alqahtani, Ali Saleh Alshomrani, "The Effects of Time Lag and Cure Rate on the Global Dynamics of HIV1 Model", BioMed Research International, vol. 2017, Article ID 8094947, 11 pages, 2017. https://doi.org/10.1155/2017/8094947
The Effects of Time Lag and Cure Rate on the Global Dynamics of HIV1 Model
Abstract
In this research article, a new mathematical model of delayed differential equations is developed which discusses the interaction among CD T cells, human immunodeficiency virus (HIV), and recombinant virus with cure rate. The model has two distributed intracellular delays. These delays denote the time needed for the infection of a cell. The dynamics of the model are completely described by the basic reproduction numbers represented by , , and . It is shown that if , then the infectionfree equilibrium is locally as well as globally stable. Similarly, it is proved that the recombinant absent equilibrium is locally as well as globally asymptotically stable if . Finally, numerical simulations are presented to illustrate our theoretical results. Our obtained results show that intracellular delay and cure rate have a positive role in the reduction of infected cells and the increasing of uninfected cells due to which the infection is reduced.
1. Introduction
Human immunodeficiency virus (HIV) is a virus that causes the condition of acquired immunodeficiency syndrome (AIDS). The virus attacks a particular type of immune system cell in the body, known as CD4 helper lymphocyte cells. HIV destroys these cells, making it harder for body to fight off other infections. Without treatment HIV1 infection passes through three different phases for HIV1 infection without treatment. The first one is the primary infection, the second is chronic infection, and the third is acquired immunodeficiency syndrome (AIDS). Several scientists and researchers are working globally to investigate an effective way to cure AIDS but they failed to completely eliminate immunodeficiency virus from the human body.
In the recent past, mathematical modeling was often used to study in vivo infection dynamics of many viruses such as HIVI, HBV, and HCV. Researchers have gained much knowledge from these models about the mechanism of the interactions of different components such as infected cells and immune system within a host and have thereby enhanced the progress in understanding the HIV1 infection. Such understanding in turn may offer guidance for developing new drugs and for designing optimal combination of existing therapies. The basic and simple model of HIV1 infection consisting of three populations, uninfected cells, infected cells, and viral particles, is governed by the following threedimensional model of nonlinear ordinary differential equations (ODEs) [1–3]:The different densities of uninfected cells and infected cells and the density of virus have been denoted by , , and , respectively. is the rate at which new susceptible cell is generated. is natural death rate of uninfected cells and is the rate of infection. is the death rate of infected cells which produce new virus particles at a rate . Rong et al. [4] extended the basic model of HIV1 infection to fourdimensional ordinary differential equation model, where latent period for the infected cells is included and a portion of these cells is reverted to the uninfected class. This improvement makes the use of nonlinear stability methods nontrivial. The authors only established the local asymptotic stability of the equilibria. However, they left the global stability of the model as an open problem. Buonomo and VargasdeLeón [5] resolved elegantly the issue left in [4] and obtained the conditions of global stability of the equilibrium states by using two distinct techniques: Lyapunov direct method and Li and Muldowney’s geometric approach. In fact, to build up a more beneficial understanding of a virus dynamics in vivo, so many authors have been devoted to studying the mechanism of infected cells reverting to the uninfected state by loss of all cccDNA from their nucleus, (see Tian and Liu [6]). Conclusive evidences that infected CD4 T cells could be cured by chemotherapy can be found in [7, 8], which is one of our motivations in the modeling for viral dynamics.
In different control measures, recombinant virus is one which is used for controlling the infection of HIV1 [9–13]. The decline of HIV1 load about fold has been proved using recombinant in vitro studies. But the efficacy of this control strategy for decreasing the viral load in AIDS patients is unknown. Genetic engineering offers an alternative approach, featuring modification of a viral genome to produce recombinant capable of controlling infections by other viruses [11]. This method has been used to modify rhabdoviruses, including the rabies and the vesicular stomatitis viruses (VSV), making them capable of infecting and killing cells previously attacked by HIV1. The engineered virus codifies the preceptor pair CD4 and CXCR4 of the host cell membrane and bind to the protein complex gp of HIV1 expressed on the surface of infected cells [9]. A basic estimation using a currently engineered virus indicated an HIV1 load reduction of percent and a recovery of host cells to percent of their normal level. Greater success ( percent HIV reduction, percent host cells recovery) is expected as more competent engineered viruses are designed. These results suggest that therapy using viruses could be an alternative to extend the survival of AIDS patients. The purpose of introducing this virus is to fight with HIV to control this infection. In [11], a new virus was introduced into model (1) and the model modified to the following form:Here the new variables and stand for recombinant virus and doubleinfected cells, respectively. The rate of production of doubleinfected cells is . The removal rate of recombinant is denoted by . is the death rate of doubleinfected cells which release recombinant at a rate . The authors of the above model analyzed the structure of equilibrium solutions and presented some simulations. Further, Jiang et al. [12] completely analyzed this model. Yu and Zou [13] modified model (2) by incorporating a control parameter to measure the injection rate of the recombinant for controlling/eliminating the HIV virus. Tian et al. [14] modified this model further by introducing the time lag into model (2) because there is time lag in infection process. They extended model (2) by introducing time delay and studied the effect of delay in controlling this infection.
In this paper, we consider that the contact process between the uninfected and virusproducing cells is not instantaneous. Thus, we include a delay, similar to the disease transmission term, in the rate of contact term. Further we also incorporate recovery rate of unproductively infected cells to uninfected cells. The recovery of these cells to uninfected cells is due to loss of all DNA from their nucleus by using drugs therapy [15, 16]. Our proposed model is extended to the following model after incorporating the abovementioned terms:where is the rate of reversion of infected cells. denotes time lag in contact and infection process. We present the dynamical behavior of the proposed model and show how delays and cure rate influence stability. We prove the wellposedness of the proposed model and study the effect of delay and cure rate in controlling HIV1. We find the basic reproduction numbers. It is shown that infectionfree equilibrium is locally as well as globally asymptotically stable. It is also shown that (recombinant absent equilibrium) is locally as well as globally asymptotically stable.
We have divided this paper into the following sections. The wellposedness and positivity of the solution are discussed in the next section. In Section 3, local and global stabilities of infectionfree equilibrium are discussed. The stability of recombinant absent equilibrium is presented in Section 4. Numerical simulation is discussed in Section 5. Finally, we have given conclusion in Section 6.
2. Positivity and WellPosedness of the Solution
This section discusses the positivity and wellposedness of system (3).
Theorem 1. All the solutions of system (3) are nonnegative provided the initial conditions are nonnegative and bounded.
Proof. Consider to be the Banach space of continuous mapping. These are the mappings from to equipped with the supnorm. For system (4), consider the initial conditions , satisfyingThere exists unique solution of system (3) under the given initial conditions (4). By using constant of variation formula, we get the following solution of system (3): which show the positivity of the solution. For boundedness of the solution , we considerThe derivative of (6) yields Here . This means that is bounded, so , and are bounded.
System (4) has the following three possible biologically meaningful equilibria [17, 18]: diseasefree equilibrium , recombinant absent equilibrium , and recombinant present equilibrium are given by Each equilibrium point can be interpreted as follows. is an infectionfree equilibrium corresponding to maximal levels of healthy CD T cells. The second equilibrium corresponds to positive levels of healthy CD T cells, infected cells, and virus, but no recombinant virus. The third equilibrium corresponds to positive levels of healthy CD T cells, infected cells, virus, and recombinant virus.
The basic reproduction number (see [19]) is obtained from the proposed model as follows: For , is the only equilibrium which is biologically meaningful. If , there is another equilibrium point . But exists if and only if , where Suppose that , and if and only if .
3. Stability of the DiseaseFree Equilibrium
The dynamical behavior of system (4) at is discussed in this section.
Theorem 2. For , the diseasefree equilibrium is locally asymptotically stable while, for , becomes unstable and the recombinant absent equilibrium occurs.
Proof. After liberalization around system (4) becomesThe characteristic equation corresponding to the Jacobian matrix of the linearized system (11) is given bywhere stands for eigenvalue. The first factor of the above equation has three negative roots and the nature of the roots of the second factor is discussed in the following:The modulus of the left hand side of (13) satisfies provided that has nonnegative real part. The modulus of the right hand side of (13) gives But this is contradiction. Thus, when , then all the eigenvalues have negative real part. Thus the infectionfree state is locally asymptotically stable. For , we have Now and . There exists at least one positive root of . Therefore, the infectionfree equilibrium is unstable if (see [20]).
Theorem 3. The diseasefree equilibrium is globally asymptotically stable when ,
Proof. Considerwhere stands for Lyapunov function. The derivative of (17) and the use of system (4) yield After further simplification, the above equation becomesThus, when , then (19) implies that and equality holds if and only if , , , , . Thus, by using LaSalle’s invariance principle (see [21]), we conclude that is globally asymptotically stable when .
4. Stability of Recombinant Absent Equilibrium
This section is devoted to the analysis of .
Theorem 4. If , then the recombinant present equilibrium is locally asymptotically stable while becomes unstable for .
Proof. The linearized form of model (4) at becomesLet be the characteristic equation of the Jacobian matrix of system (20), where Now can be simplified as which indicates that has two roots with negative real part if and only if (i.e., ) or one positive and one negative root if (i.e., ). Therefore, if , then the single infection equilibrium is unstable. After some simplification can be written aswhere It is clear that is not a root of (24) if . When , (24) becomesUsing RouthHurwitz criterion (see [22]), we can prove that Similarly, Thus, any root of (24) has negative real part when . Now we consider the distribution of the roots when . Let be the pure imaginary root of (24). Then, we obtain The modulus of the above equation gives the following result:Since we see that all the coefficients of the above equation are positive which implies that the function is monotonically increasing for with . Therefore, (30) has no positive roots if . Equation (24) has all the roots with negative real part if and .
Theorem 5. For , the recombinant present equilibrium is globally asymptotically stable.
Proof. Let us construct the Lyapunov functionalThe derivative of (32) yieldsModel (1) at single infection equilibrium becomes If is very large, that is, when the time delay in the contact of uninfected targeted cells and pathogen virus is large and the latent period is very large, then the rate of infection will be very small and contrarily if is very small, then the infection will spread more rapidly. Therefore, we suppose that delay is very large, and taking limit we getUsing the above identities and assumption (33) in (35), we getThe following inequalities hold [14]: Therefore, by using the above inequalities, (36) implies that when , and the equality holds when , , , , Therefore, by LaSalle’s invariance principle [21], we conclude that is globally asymptotically stable.
5. Numerical Simulation
In this section, we present the numerical simulations by using MATLAB to illustrate our theoretical results. The drugs therapy can control the HIV1. Using drugs therapy, the infected cells revert to the uninfected cells. For numerical simulation, we consider the values of the parameters presented in Table 1 [11].

Figures 1–3 are the oscillations of uninfected cells, infected cells, doubleinfected cells, pathogen virus, and recombinant virus. Figure 1 shows the dynamical behavior of HIV1 infection for the delay term and for different recovery rates and represents that as the value of the recovery rate increases the density of of uninfected cells increases and the concentration of infected cells decreases. Figure 2 shows that by varying time delay and keeping the values of constant, the amplitude of oscillation increases and the rate towards stability decreases. Figure 3 shows that if we further reduce the delay time , then amplitudes of oscillations increases. More importantly, it is noted that the amplitudes of the oscillations in Figure 3 are almost double of that in Figure 1 though their frequencies are almost not changed. These figures show that introducing even very small time delay in the model can produce significant quantitative changes in solutions, which cannot be observed from the model without delay. Also, as the value of recovery rate increases the infected cells revert to the healthy cells more rapidly and converge to stable equilibrium. We can see that the infection would always keep stability when the cure rate γ is larger. Therefore, we can also claim that the cure rate γ is a very important parameter and by improving the cure rate, we will control the disease. Moreover, the significant qualitative changes due to existence of delay can be observed. These results also suggest that the delay is very important fact which should not be missed.
6. Conclusion
In this paper, a delayed HIV1 model with drugs therapy is presented. The improved model with delay has three equilibrium solutions , , and . It has been shown that is locally as well as globally asymptotically stable for , which loses its stability at . Then, bifurcates into . Next, it is also proved that is also locally and globally asymptotically stable for . Delay, as the bifurcation parameter, plays a very important role in determining the dynamic behavior of the system. Delay may change the dynamical behavior quantitatively, even in the normal range of values. This indeed suggests that delay is a very important fact which should not be missed in HIV1 modeling. The drugs therapy also has an important effect on model (3). As the value of recovery rate increased the infected cells revert to the uninfected cells resulting in decrease in infected cells and increase in healthy cells. And this infection can easily be controlled if we improve the cure rate.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 A. S. Perelson, D. E. Kirschner, and R. D. Boer, “Dynamics of HIV infection of CD4^{+} T cells,” Mathematical Biosciences, vol. 114, no. 1, pp. 81–125, 1993. View at: Publisher Site  Google Scholar
 A. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, and D. D. Ho, “HIV1 dynamics in vivo: virion clearance rate, infected cell lifespan, and viral generation time,” Science, vol. 271, no. 5255, pp. 1582–1586, 1996. View at: Publisher Site  Google Scholar
 A. S. Perelson and P. W. Nelson, “Mathematical analysis of HIV1 dynamics in vivo,” Society for Industrial and Applied Mathematics, vol. 41, no. 1, pp. 3–44, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 L. Rong, M. A. Gilchrist, Z. Feng, and A. S. Perelson, “Modeling withinhost HIV1 dynamics and the evolution of drug resistance: tradeoffs between viral enzyme function and drug susceptibility,” Journal of Theoretical Biology, vol. 247, no. 4, pp. 804–818, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 B. Buonomo and C. VargasdeLeón, “Global stability for an HIV1 infection model including an eclipse stage of infected cells,” Journal of Mathematical Analysis and Applications, vol. 385, no. 2, pp. 709–720, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Tian and X. Liu, “Global dynamics of a virus dynamical model with general incidence rate and cure rate,” Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal, vol. 16, pp. 17–26, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 J. A. Zack, A. M. Haislip, P. Krogstad, and I. S. Y. Chen, “Incompletely reversetranscribed human immunodeficiency virus type 1 genomes in quiescent cells can function as intermediates in the retroviral life cycle,” Journal of Virology, vol. 66, no. 3, pp. 1717–1725, 1992. View at: Google Scholar
 J. A. Zack, S. J. Arrigo, S. R. Weitsman, A. S. Go, A. Haislip, and I. S. Y. Chen, “HIV1 entry into quiescent primary lymphocytes: Molecular analysis reveals a labile, latent viral structure,” Cell, vol. 61, no. 2, pp. 213–222, 1990. View at: Publisher Site  Google Scholar
 M. J. Schnell, J. E. Johnson, L. Buonocore, and J. K. Rose, “Construction of a novel virus that targets HIV1infected cells and controls HIV1 infection,” Cell, vol. 90, no. 5, pp. 849–857, 1997. View at: Publisher Site  Google Scholar
 C. P. Nolan, “Harnessing viral devices as pharmaceuticals: Fighting HIV1's fire with fire,” Cell, vol. 90, no. 5, pp. 821–824, 1997. View at: Publisher Site  Google Scholar
 T. Revilla and G. GarcíaRamos, “Fighting a virus with a virus: a dynamical model for HIV1 therapy,” Mathematical Biosciences, vol. 185, no. 2, pp. 191–203, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 X. Jiang, P. Yu, Z. Yuan, and X. Zou, “Dynamics of an HIV1 therapy model of fighting a virus with another virus,” Journal of Biological Dynamics, vol. 3, no. 4, pp. 387–409, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 P. Yu and X. Zou, “Bifurcation analysis on an HIV1 model with constant injection of recombinant,” International Journal of Bifurcation and Chaos, vol. 22, no. 3, Article ID 1250062, 21 pages, 2012. View at: Publisher Site  Google Scholar
 Y. Tian, Y. Bai, and P. Yu, “Impact of delay on HIV1 dynamics of fighting a virus with another virus,” Mathematical Biosciences and Engineering, vol. 11, no. 5, pp. 1181–1198, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 K. Hattaf and N. Yousfi, “Two optimal treatments of HIV infection model,” World Journal of Modelling and Simulation, vol. 8, no. 1, pp. 27–36, 2012. View at: Google Scholar
 B. El Boukari and Y. Noura, “A delay differential equation model of HIV infection, with therapy and CTL response,” Bulletin of Mathematical Sciences and Applications, vol. 9, pp. 53–68, 2014. View at: Publisher Site  Google Scholar
 N. Ali, G. Zaman, and O. Algahtani, “Stability analysis of HIV1 model with multiple delays,” Advances in Difference Equations, vol. 88, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 R. V. Culshaw, S. Ruan, and G. Webb, “A mathematical model of celltocell spread of HIV1 that includes a time delay,” Journal of Mathematical Biology, vol. 46, no. 5, pp. 425–444, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 A. V. M. Herz, S. Bonhoeffer, R. M. Anderson, R. M. May, and M. A. Nowak, “Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 14, pp. 7247–7251, 1996. View at: Publisher Site  Google Scholar
 J. E. Mittler, M. Markowitz, D. D. Ho, and A. S. Perelson, “Improved estimates for HIV1 clearance rate and intracellular delay,” AIDS, vol. 13, no. 11, pp. 1415–1417, 1999. View at: Publisher Site  Google Scholar
 J. P. LaSalle, The Stability of Dynamical Systems, SIAM, Philadelphia, Pa, USA, 1976. View at: MathSciNet
 F. R. Gantmacher, The Theory of Matrices, vol. 1, Chelsea Publishing, New York, NY, USA, 1960. View at: MathSciNet
Copyright
Copyright © 2017 Nigar Ali et al. 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.