Abstract and Applied Analysis

Volume 2013 (2013), Article ID 163484, 12 pages

http://dx.doi.org/10.1155/2013/163484

Research Article

Global Stability of HIV-1 Infection Model with Two Time Delays

Hui Miao, Xamxinur Abdurahman, and Ahmadjan Muhammadhaji

College of Mathematics and System Science, Xinjiang University, Xinjiang, Urumqi 830046, China

Received 8 October 2013; Accepted 19 November 2013

Academic Editor: Ming Mei

Copyright © 2013 Hui Miao 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.

Abstract.

We investigate global dynamics for a system of delay differential equations which describes a virus-immune interaction in vivo. The model has two time delays describing time needed for infection of cell and CTLs generation. Our model admits three possible equilibria: infection-free equilibrium, CTL-absent infection equilibrium, and CTL-present infection equilibrium. The effect of time delay on stability of the equilibria of the CTL immune response model has been studied.

1. Introduction

As well known, in recent years the population dynamics of infectious diseases have been extensively studied [1–20]. Particularly, the HIV (human immunodeficiency virus) has been extensively studied in [1–5, 7–11, 14, 15, 20] and became a global problem. The HIV infection is characterized by three different phases, namely, the primary infection, clinically asymptomatic stage (chronic infection), and acquired immunodeficiency syndrome (AIDS) or drug therapy. During primary infection, viral load in the peripheral blood experiences a substantial increase to the peak level, followed by decline to the steady state, which is referred to as the viral set point. Extremely high viral load during primary infection leads to the activation of T cells, which are recognized as cytotoxic T cells (CTL) capable of suppressing viral replication. Viral decline from the peak is due to the control by these immune cells and/or limited target cell availability. The viral set point has been shown to be predictive for the pace of disease development [8]. Clinical research combined with mathematical modeling has enhanced progress in the understanding of HIV-1 infection [4]. This is because mathematical models can offer a way to study the dynamics of viral load in vivo and can be very useful in understanding the interaction between virus and host cell.

On the other hand, in the real situation, there may be a lag between the time target cells are contacted by the virus particles and the time the contacted cells become actively affected meaning that the contacting virions enter cells. This can be explained by the initial (or eclipse) phase of the virus life cycle, which include all stages from viral attachment until the time that the host cell contains the infectious viral particles in its cytoplasm [3]. Research [8] has shown that models of HIV-1 infection that include intracellular delays are more accurate representations of the biology and change the estimated values of kinetic parameters when compared to models without delays. Therefore, we should introduce time delays into model foundation, which will have more resemblance to the real ecosystem. In the last decade, the HIV-infection models with time delay have been studied by many authors, and time delays of one type or another have been incorporated into biological models by many authors (e.g., [1–5, 8–11, 14, 15, 20] and the references cited therein). Here we include an intracellular delay as well as immune delay. There are some models which include an intracellular delay [1–5]; some authors believe that time delays can not be ignored in models for immune response [6, 7].

The salient features of the mechanism of the immune response during viral infection are as follows. First, the free virus enters its target, a susceptible cell. Inside this cell it replicates itself. And this susceptible cell becomes an infected cell. Then the infected cell dies and releases new viruses; these viruses begin to infect other susceptible cells. During the process of viral infection, the host is induced which is initially rapid and nonspecific (natural killer cells, macrophage cell, etc.) and then delayed and specific (cytotoxic T lymphocyte cells, antibody cell). But in most virus infections cytotoxic T lymphocyte (CTL) cells which attack infected cells, and antibody cell which attack viruses, play a critical part in antiviral defense. In order to investigate the role of the population dynamics of viral infection with CTL response, Pawelek et al. [8] constructed a mathematical model describing the basic dynamics of the interaction between the uninfected target cells, productively infected cells, free virus, and the CTL response cells, which is described by the following differential equation:

where the uninfected target cells are denoted by  , productively infected cells are denoted by  , free virus is denoted by  , and the CTL response cells are denoted by  . The parameter represents the rate at which new target cells are created, is the death rate of uninfected target cells, is the infection rate of uninfected cells by virus, is the death rate of productively infected cells, represents the killing rate of infected cells by CTL response cells, is the rate of the virus particles produced by infected cells, is the viral clearance rate constant, is the rate at which the CTL response is produced, and is the death rate of the CTL response.

In order to incorporate the intracellular phase of the virus life cycle, we assume that virus production occurs after the virus entry by the constant delay . The recruitment of virus-producing cells at time is given by the number of the uninfected T cells that were newly infected at time and are still alive at time [9, 10]. If we assume constant death rates for infected T cells but not yet producing virus particles, the probability of surviving the time period from to is . And the immune response plays an important role in eliminating or controlling the disease after human body is infected by virus. Antigenic stimulation generating CTLs may need a period of time ; that is, the CTL response at time may depend on the population of antigen at a period time [11].

Therefore, based on the discussion above, the model can be written in the following form:

The organization of this paper is as follows. In the next section we deal with some basic properties such as positivity and boundedness of the solutions and existence of equilibria of system (2). In Section 3, we prove the local stability of three possible equilibria. Further, by using the well-known Lyapunov-Lasalle invariance principle, we prove the global asymptotic stability of the infection-free equilibrium, CTL-absent infection equilibrium, and a special case of CTL-present equilibrium. By bifurcation theory, we also can prove that there is a stability switch for another special case for CTL-present equilibrium. In Section 4, one example is given to illustrate that our main results are applicable. In the final section, we offer a brief conclusion.

2. Basic Results

Based on the biological meaning, we will consider the system (2) with the following initial conditions:

Let be the Banach space of continuous functions mapping the internal into equipped with the sup-norm, where . Based on the existence and uniqueness theory of solution for functional differential equations [12, 13], it is easy to show that there is a unique solution to system (2) with initial condition (3).

First we will discuss the positivity and boundedness of the solution.

Theorem 1. Let be the solution of (2) satisfying conditions (3); then , and are positive and ultimately bounded.

Proof. From (2), we have

This means that all the solutions of system (2) with initial condition (3) are positive. To prove the boundedness of the solution, denote

Calculating the derivative of along the solution of system (2) and by positivity of the solutions, we have

where . This implies that is bounded for large . So , and are ultimately bounded.

Next we will discuss the equilibria of system (2). The following is the basic reproductive ratio of system (2) which describes the average number of newly infected cells generated from one infected cell at the beginning of the infectious process:

By direct calculation we have that system (2) has three equilibria. Infection-free equilibrium . If , there is a CTL-absent infection equilibrium

If , then there is a CTL-present infection equilibrium

3. Stability Analysis of Delay Model (2)

In this section, we will analyse locally and globally asymptotic stability of the three equilibria.

3.1. Stability of Infection-Free Equilibrium

Theorem 2. For the infection-free equilibrium of system (2);

(i)if , then is locally asymptotically stable; (ii)if , then is unstable.

Proof. First, we will prove the local stability of . At the infection-free equilibrium , the characteristic equation for the corresponding linearized system of (2) becomes

Two of the roots of the characteristic equation (10) are , and . The remaining two roots are obtained by considering the following equation:

If has a nonnegative real part, then the modulus of the left-hand side of (11) satisfies

while the modulus of the right-hand side (11) satisfies

This leads to a contradiction. Thus, when , all the eigenvalues have negative real parts, and hence the infection-free steady state is locally asymptotically stable.

When , let

Then , . By the continuity of , there exists at least one positive root of . Thus, the infection-free equilibrium is unstable if .

Theorem 3. If , then infection-free steady state is globally asymptotically stable.

Proof. Define a Lyapunov functional

where . Calculating the time derivative of along the solution of system (2), we obtain

Note that if and only if . By the second equation of (2), we also have . Therefore, the maximal compact invariant set in is the singleton . By the Lasalle Invariance principle [12], the infection-free steady state is globally attracting. Therefore, is globally asymptotically stable.

3.2. Stability of the CTL-Absent Infection Equilibrium

Theorem 4. For the CTL-absent infection equilibrium of system (2),

(i)if , then is locally asymptotically stable; (ii)if , then is unstable.

Proof. At the CTL-absent infection equilibrium , the characteristic equation for the corresponding linearized system of (2) is

where

First we consider

When , (19) becomes

As , , , by Routh-Hurwitz criterion, we only need to show the following statement:

Thus all the roots of (19) have negative real parts when .

Now we consider the distribution of roots of (19) while . If is a solution of (19), separating real and imaginary parts, it follows that

Squaring and adding the two equations of (22) yield

Letting , (23) becomes

where

By the Routh-Hurwitz criterion, (24) has no positive roots. This shows that (19) can not have a purely imaginary root.

Next, we analyze the transcendental equation

For , if , we have . This shows that the root of (26) is negative for .

Now we only need to consider (26) in the case . By letting be a purely imaginary root of (26) for some , we have

which implies that . Note that ; then , which is a contradiction.

Therefore, we conclude that the characteristic equation (17) does not have any root with nonnegative real part. By the general theory of delay differential equations from Kuang [12], we see that if , the equilibrium is locally asymptotically stable.

Theorem 5. If , then the CTL-absent infection equilibrium of system (2) is globally asymptotically stable.

Proof. Denote , . Define a Lyapunov functional

By calculating the derivative of along the solution of system (2), we obtain that

Remember that

The equation above can be rewritten as

From (31), it follows that for all . if and only if . Then the globally asymptotic attractivity of follows from Lyapunov-Lasalle invariance principle [12]. Therefore, is globally asymptotically stable.

3.3. Stability of CTL-Present Infection Equilibrium

On the stability analysis of CTL-present infection equilibrium , we only consider the two special cases, that is, and .

Firstly, we consider the local and global stability of for the first case, and we have the following results.

Theorem 6. If , and , then the CTL-present infection equilibrium is locally asymptotically stable.

Proof. At equilibrium , the characteristic equation for the corresponding linearized system of (2) is

where . Further let ; then (32) becomes

By the continuous dependence of roots of the characteristic equation on , we know that the curve of the roots must cross the imaginary axis as decreases sufficiently close to 1. That is, the characteristic equation (33) has a pure imaginary root if and only if the following statement is true:

We claim that the following inequality holds:

In fact, we have

It follows from , , and the inequality (35) that the modulus of the left-hand side of (34) is greater than the modulus of the right-hand side. This leads to a contradiction. Therefore, we conclude that (33) does not have any root with nonnegative real part. Thus, the CTL-present infection equilibrium of system (2) is locally asymptotically stable when in the case of and .

Theorem 7. If , , and , then the is globally asymptotically stable.

Proof. Define a Lyapunov functional

where , , , and satisfy the following equations:

Calculating the derivative of along the solution of system (2) and using the similar method with the proof of Theorem 5, we have

Since , it follows that for all and if and only if . Then the global attractivity of follows from Lyapunov-Lasalle invariance principle [12]. Therefore the CTL-present infection equilibrium of system (2) is globally asymptotically stable.

The following discussions focus on the stability of the equilibrium in the second case. Let in (32); it follows that that

where

Let

Then we can rewrite (40) as follows:

Assume that and are analytic functions in the right half-plane , where .

Lemma 8. Consider (43); and satisfy the following conditions:

(i); (ii); (iii) for any ; (iv) for real , has at most a finite number of real zeros; (v)each positive root of is continuous and differentiable at whenever it exists.

Proof. (i). This means is not a characteristic root of (43).

(ii) Consider

(iii) Since

we have for any .

(iv) Since

It is obvious that property (iv) is satisfied.

(v) The last assertion is valid because is a quadratic polynomial in and the fact that are all continuous functions of . This completes the proof.

Let and denote , , , and ; then (46) becomes

On the distribution of positive roots of (46), we have the following results.

Lemma 9. If , then (46) has at least one positive root.

Proof. Since and , . Hence, there exists at least an such that .

In the following, we will analyze the case . Let ; then (47) can be rewritten as

where , , . If , then it is easy to obtain the four roots of (48) as follows

where . Thus are the roots of (47). Then we have the following result.

Lemma 10. Suppose that and .

(i)If , then (47) has no positive real roots. (ii)If , and , then (47) has no positive real roots.

Based on Theorem  4.1 in [12] and Lemmas 8–10, we have the following results.

Theorem 11. (i) If , ;

(a)(b), , and ,

the equation has no positive roots. Then system (2) is stable for all .

(ii) If , the equation has at least one positive root and each positive root is simple; then stability switches may occur as increases. That is, there exists a positive number such that (2) is unstable for all . As varies from 0 to , at most a finite number of stability switches may occur.

4. Numerical Simulations

In the previous sections, we studied dynamical behaviors of the system (2) and obtained some important results. In this section, we perform a numerical analysis of the model based on Theorem 11.

Example 12. By corresponding to system (2), we consider the following system:

where . By direct calculation we can get CTL-present infection equilibrium for system (50).

From the parameters given in system (50) and the values of , we can see that there are two critical values of the delay , denoted by and , . Simple numerical simulations show that the CTL-present infection equilibrium of system (50) is globally asymptotically stable for (see Figure 1). In this case, we take . The above CTL-present infection equilibrium of system (50) is unstable for (see Figure 2). In this case, we take . The above CTL-present infection equilibrium of system (50) is globally asymptotically stable for (see Figure 3). In this case, we take . From above analysis we can conclude that Hopf bifurcation occurs when .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 1: The time histories and the phase trajectories of system (50) before Hopf bifurcation occurs for .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 2: The time histories and the phase trajectories of the system (50) when Hopf bifurcation occurs for .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 3: The time histories and the phase trajectories of the system (2) after Hopf bifurcation occurs for .

5. Conclusions

In this paper, we have discussed HIV infection model with intracellular delay and CTL-response delay. We assume that the production of CTLs depends on the infected cells and CTL cells for some important biological meanings. Dynamical analysis of system (2) shows that intracellular delay and immune delay play different roles in the stability of the equilibrium. The results show that when , the infection-free equilibrium is globally asymptotically stable, which means that the viruses are cleared and immune is not active. When , the CTL-absent infection equilibrium exists and is globally asymptotically stable, which means that the CTL immune response would not be activated and viral infection becomes vanished. When and , the CTL-present infection equilibrium is globally asymptotically stable. Actually, the intracellular delay does not affect the stability of the system. When and , system (2) may undergo a stability switch.

Disclosure

The authors declare that they have no financial or personal relationships with other people or organizations that can inappropriately influence their work; there is no professional or other personal interest of any nature or kind in any product, service and/or company that could be construed as influencing the position presented in, or the paper.

Acknowledgment

This work was supported by the National Natural Science Foundation of China (Grants nos. 11261056, 11261058, and 11271312).

References

R. V. Culshaw and S. Ruan, “A delay-differential equation model of HIV infection of CD4⁺ T-cells,” Mathematical Biosciences, vol. 165, no. 1, pp. 27–39, 2000.

R. V. Culshaw, S. Ruan, and G. Webb, “A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay,” Journal of Mathematical Biology, vol. 46, no. 5, pp. 425–444, 2003.

P. W. Nelson and A. S. Perelson, “Mathematical analysis of delay differential equation models of HIV-1 infection,” Mathematical Biosciences, vol. 179, no. 1, pp. 73–94, 2002.

P. W. Nelson, J. D. Murray, and A. S. Perelson, “A model of HIV-1 pathogenesis that includes an intracellular delay,” Mathematical Biosciences, vol. 163, no. 2, pp. 201–215, 2000.

H. Zhu and X. Zou, “Impact of delays in cell infection and virus production on HIV-1 dynamics,” Mathematical Medicine and Biology, vol. 25, no. 2, pp. 99–112, 2008.

A. A. Canabarro, I. M. Gléria, and M. L. Lyra, “Periodic solutions and chaos in a non-linear model for the delayed cellular immune response,” Physica A, vol. 342, no. 1-2, pp. 234–241, 2004.

K. Wang, W. Wang, H. Pang, and X. Liu, “Complex dynamic behavior in a viral model with delayed immune response,” Physica D, vol. 226, no. 2, pp. 197–208, 2007.

K. A. Pawelek, S. Liu, F. Pahlevani, and L. Rong, “A model of HIV-1 infection with two time delays: mathematical analysis and comparison with patient data,” Mathematical Biosciences, vol. 235, no. 1, pp. 98–109, 2012.

M. Y. Li and H. Shu, “Global dynamics of an in-host viral model with intracellular delay,” Bulletin of Mathematical Biology, vol. 72, no. 6, pp. 1492–1505, 2010.

H. Zhu and X. Zou, “Dynamics of a HIV-1 infection model with cell-mediated immune response and intracellular delay,” Discrete and Continuous Dynamical Systems B, vol. 12, no. 2, pp. 511–524, 2009.

C. A. Janeway, P. Travers, M. Walport, and M. J. Schlomchik, Immunobiology, Garland Science Publishing, 6th edition, 2005.

Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Academic Press, San Diego, Calif, USA, 1993.

J. K. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 1997.

X. Wang, A. Elaiw, and X. Song, “Global properties of a delayed HIV infection model with CTL immune response,” Applied Mathematics and Computation, vol. 218, no. 18, pp. 9405–9414, 2012.

H. Zhu, Y. Luo, and M. Chen, “Stability and Hopf bifurcation of a HIV infection model with CTL-response delay,” Computers and Mathematics with Applications, vol. 62, no. 8, pp. 3091–3102, 2011.

E. Beretta and Y. Kuang, “Geometric stability switch criteria in delay differential systems with delay dependent parameters,” SIAM Journal on Mathematical Analysis, vol. 33, no. 5, pp. 1144–1165, 2002.

M. A. Nowak and C. R. M. Bangham, “Population dynamics of immune responses to persistent viruses,” Science, vol. 272, no. 5258, pp. 74–79, 1996.

J. K. Hale and S. Lunel, Introduction to Functional Differential Equations, Springer, New York, NY, USA, 1993.

T. Zhang, J. Liu, and Z. Teng, “Stability of Hopf bifurcation of a delayed SIRS epidemic model with stage structure,” Nonlinear Analysis: Real World Applications, vol. 11, no. 1, pp. 293–306, 2010.

X. Song, X. Zhou, and X. Zhao, “Properties of stability and Hopf bifurcation for a HIV infection model with time delay,” Applied Mathematical Modelling, vol. 34, no. 6, pp. 1511–1523, 2010.