Stability and Hopf Bifurcation in a Delayed HIV Infection Model with General Incidence Rate and Immune Impairment.

We investigate the dynamical behavior of a delayed HIV infection model with general incidence rate and immune impairment. We derive two threshold parameters, the basic reproduction number R 0 and the immune response reproduction number R 1. By using Lyapunov functional and LaSalle invariance principle, we prove the global stability of the infection-free equilibrium and the infected equilibrium without immunity. Furthermore, the existence of Hopf bifurcations at the infected equilibrium with CTL response is also studied. By theoretical analysis and numerical simulations, the effect of the immune impairment rate on the stability of the infected equilibrium with CTL response has been studied.

1. Introduction

In recent years, mathematical models have been proved to be valuable in understanding the dynamics of viral infection (see, e.g., [1-8]). In most virus infections, cytotoxic T lymphocyte (CTL) cells play a significant role in antiviral defense by attacking virus-infected cells. In order to study the role of the population dynamics of the viral infection with CTL response, Nowak and Bangham et al. proposed a basic viral infection model describing the interactions between a replicating virus population and a specific antiviral CTL response, which takes into account four populations: uninfected cells, actively infected cells, free virus, and CTL cells (see, e.g., [1–4, 9, 10]). Now, the population dynamics of viral infection with CTL response has been paid much attention and many properties have been investigated (see, e.g., [11-16]).

Furthermore, the state of latent infection cannot be ignored in many biological models. The infected cells are separated into two distinct compartments, latently infected and actively infected. These latently infected cells do not produce virus and can evade from viral cytopathic effects and host immune mechanisms (see, e.g., [17-20]). Recently, the following model with latent infection and CTL response has been proposed (see, e.g., [11]):where x(t), u(t), y(t), v(t), and z(t) represent the numbers of uninfected cells, latently infected cells, actively infected cells, free virus, and CTLs at time t, respectively. Uninfected cells are produced at the rate λ, die at the rate μ 1, and become infected at the rate β. The constant σ is the rate of latently infected cells translating to actively infected cells and μ 3 is the death rate of actively infected cells. The constant μ 2 represents the death rate of latently infected cells. The constant p is the rate of CTL-mediated lysis and q is the rate of CTL proliferation. The constant k is the rate of production of virus by infected cells and μ 4 is the clearance rate of free virus. The removal rate of CTLs is μ 5.

However, in plenty of previous papers, many models are constructed under the assumption that the presence of antigen can stimulate immunity and ignore the immune impairment (see, e.g., [8, 11, 16, 17]). In fact, some pathogens can also suppress immune response or even destroy immunity especially when the load of pathogens is too high such as HIV, HBV (see, e.g., [15, 21–25]). Regoes et al. consider an ordinary differential equation (ODE) model with an immune impairment term myz (see, e.g., [12, 26, 27]), where m denotes the immune impairment rate. Time delay should be considered in models for CTL response. It is shown that time delay plays an important role to the dynamic properties in models for CTL response (see, e.g., [1, 5, 6, 8, 15]). In fact, antigenic stimulation generating CTLs may need a period of time t; that is, the CTL response at time t may depend on the numbers of CTLs and infected cells at time t − τ, for a time lag τ > 0 (see, e.g., [1, 5, 13]).

Motivated by the above works, in this paper, we will study a delay differential equation (DDE) model of HIV infection with immune impairment and delayed CTL response. Furthermore, we know that the actual incidence rate is probably not linear over the entire range of x and v. Based on the works mentioned above (see, e.g., [21, 28–31]), we propose the following system with general incidence function:where the state variables x(t), u(t), y(t), v(t), and z(t) and the parameters λ, σ, p, k, q, μ 1, μ 2, μ 3, μ 4, and μ 5 have the same biological meaning as in system (1). m is the immune impairment rate. Suppose all the parameters are nonnegative. We assume the incidence rate is the general incidence function f(x, v)v, where f ∈ C 1([0, +∞]×[0, +∞], R) satisfies the following hypotheses:

f(x, v)v ≥ 0, for all x ≥ 0 and v ≥ 0; f(x, v) = 0 if and only if x = 0;

∂f(x, v)/∂x > 0, for all x ≥ 0 and v ≥ 0;

∂f(x, v)/∂v ≤ 0, for all x ≥ 0 and v ≥ 0;

∂(f(x, v)v)/∂v > 0, for all x > 0 and v ≥ 0.

Clearly, the hypotheses can be satisfied by different types of the incidence rate including the mass action, the Holling type II function, the saturation incidence, Beddington-DeAngelis incidence function, Crowley-Martin incidence function, and the more generalized incidence functions (see, e.g., [4, 6, 17, 32, 33]). Further, in order to study the global stability of the equilibria of system (2) by the method of Lyapunov functionals, we assume the following hypotheses hold (see, e.g., [28]):

x − x 0 − ∫ (f(x 0, 0)/f(s, 0))ds → +∞, as x → +∞ or x → 0+;

x − x 1 − ∫ (f(x 1, v 1)/f(s, v 1))ds → +∞, as x → +∞ or x → 0+;

x − x − ∫ (f(x , v )/f(s, v ))ds → +∞, as x → +∞ or x → 0+.

The main purpose of this paper is to carry out a complete theoretical analysis on the global stability of the equilibria of system (2). The organization of this paper is as follows. In Section 2, we consider the nonnegativity and boundedness of the solutions and the existence of the equilibria of system (2). In Section 3, we consider the global stability of the infection-free equilibrium E 0 and the infected equilibrium without immunity E 1 by constructing suitable Lyapunov functionals and using LaSalle invariance principle. In Section 4, we discuss the local stability of the infected equilibrium with CTL response E and the existence of Hopf bifurcations. Finally, in Section 5, the brief conclusions are given and some numerical simulations are carried out to illustrate the main results.

2. Basic Results

2.1. The Nonnegativity and Boundedness of the Solutions

According to biological meanings, the initial condition of system (2) is given as follows:where θ ∈ [−τ, 0] and (φ 1, φ 2, φ 3, φ 4, φ 5) ∈ C = C([−τ, 0], R + 5) and C is the Banach space of the continuous functions mapping the interval [−τ, 0] into R + 5, R + 5 = {(x 1, x 2, x 3, x 4, x 5)∣x ≥ 0, i = 1,2, 3,4, 5}.

Under the initial condition (3), it easily shows that the solution of system (2) is unique and nonnegative for all t ≥ 0 and ultimately bounded. It has the following result.

Proposition 1 .

Under the initial condition (3), the solution of system (2) is unique and nonnegative for all t ≥ 0 and also ultimately bounded, when (H1)–(H7) are satisfied.

Proof

The uniqueness and nonnegativity of the solution (x(t), u(t), y(t), v(t), z(t)) can be easily proved by using the theorems in [34, 35].

Next, for t ≥ 0, define By the nonnegativity of the solutions, it follows that, for t ≥ 0, where γ = min{μ 1, μ 2, μ 3/2, μ 4, μ 5}. Thus, it has that limsup L(t) ≤ λ/γ, from which it has that the solution (x(t), u(t), y(t), v(t), z(t)) is ultimately bounded.

2.2. The Existence of the Equilibria

Next, we consider the existence of the equilibria. The equilibrium of system (2) satisfies

If u = 0, y = 0, v = 0, and z = 0, system (2) has only one equilibrium, that is, the infection-free equilibrium E 0 = (x 0, 0,0, 0,0), where x 0 = λ/μ 1.

If u ≠ 0, y ≠ 0, v ≠ 0, and z = 0, we have Since v > 0, we have that x < λ/μ 1. Hence, we only need to consider the case of x < λ/μ 1.

Consider the following function defined on the interval (0, λ/μ 1) by Under hypotheses (H2) and (H3), we have We know that the function F(x) is strictly monotonically increasing with respect to x. Denote the basic reproduction number R 0 of system (2) by Clearly, we have It has that there exists a unique x 1 ∈ (0, λ/μ 1) such that F(x 1) = 0, if R 0 > 1. Then we can compute u 1, y 1 and v 1 by (8). Hence, we get the unique infected equilibrium without immunity E 1 = (x 1, u 1, y 1, v 1, 0).

If z ≠ 0 and q > m, we get the following equations:

Since z > 0, we have , where Hence, the existence of the equilibrium requires and (13) has a solution on the interval .

Denote Hence, if , it has . Denote Under hypothesis (H2), we know that the function G(x) is strictly monotonically increasing with respect to x. Clearly, we have where Hence, we have that there exists such that G(x ) = 0, if and R 1 > 1. Then we can compute u , y , v , and z by (14) and (15).

Denote the immune response reproduction number of system (2) as R 1. Therefore, we have that there exists a unique infected equilibrium with CTL response E = (x , u , y , v , z ), if and R 1 > 1. This proves the following theorem.

Theorem 2 .

Suppose that hypotheses (H1)–(H4) are satisfied; the following conclusions hold.

System (2) always has an infection-free equilibrium E 0.

System (2) has an infected equilibrium without immunity E 1 if R 0 > 1.

System (2) has an infected equilibrium with immunity E if and R 1 > 1.

From hypotheses (H1)–(H3), it is clear that R 1 < R 0. In order to study the global stability of the infected equilibrium E 1 in the next section, we give the following remark.

Remark 3 .

Suppose that is satisfied; then the following results hold:

If R 1 > 1, then (q − m)y 1/μ 5 > 1.

If R 1 ≤ 1, then (q − m)y 1/μ 5 ≤ 1.

Let us give the proof of Remark 3. Firstly, for Case (i), since R 1 > 1, then Since the function F(x) is strictly monotonically increasing with respect to x and F(x 1) = 0, we have . Therefore Then

Secondly, for Case (ii), since R 1 ≤ 1, then We have . Therefore Then

3. The Global Stability of the Equilibria

In this section, we study the global stability of the equilibria of system (2). Firstly, we analyze the global stability of the infection-free equilibrium E 0.

Theorem 4 .

Suppose that hypotheses (H1)–(H7) are satisfied. If R 0 ≤ 1, then the infection-free equilibrium E 0 is globally asymptotically stable for any time delay τ ≥ 0. If R 0 > 1, then the infection-free equilibrium E 0 is unstable for any time delay τ ≥ 0.

Let (x(t), u(t), y(t), v(t), z(t)) be a positive solution of system (2) with the initial condition (3) for t ≥ 0. Motivated by the works in [14, 28, 31, 36, 37], we consider the following Lyapunov functional:where λ = μ 1 x 0. By (H1)–(H5), it is obvious that V 1 is positive definite with respect to E 0. For t ≥ 0, the time derivative of V 1 along the solutions of system (2) is

Since hypotheses (H1)–(H3) and R 0 ≤ 1, we haveTherefore, if R 0 ≤ 1. Then it follows from stability theorems in [34, 35] that the infection-free equilibrium E 0 is stable for any time delay τ ≥ 0 if R 0 ≤ 1.

Furthermore, note that, for each t ≥ 0, implies that x(t) = x 0, z(t) = 0. Let M be the largest invariant set in the set We have from the first four equations of system (2) and the invariance of M that M = {E 0}. Since any solution of system (2) is bounded, it follows from LaSalle invariance principle (see, e.g., [34, 35]) that the infection-free equilibrium E 0 is also globally attractive for any time delay τ ≥ 0 if R 0 ≤ 1.

The characteristic equation of system (2) at the infection-free equilibrium E 0 isClearly, if R 0 > 1, (31) has at least a positive real root. Thus, the infection-free equilibrium E 0 is unstable.

Next we study the global stability of the infected equilibrium without immunity E 1.

Theorem 5 .

Suppose that hypotheses (H1)–(H7) and are satisfied. If R 0 > 1 ≥ R 1, then the infected equilibrium without immunity E 1 is globally asymptotically stable for any time delay τ ≥ 0. If R 1 > 1, then the infected equilibrium without immunity E 1 is unstable for any time delay τ ≥ 0.

Let (x(t), u(t), y(t), v(t), z(t)) be a positive solution of system (2) with the initial condition (3) for t ≥ 0. Consider the following Lyapunov functional:

Let ψ(x) = x − x 1 − ∫ (f(x 1, v 1)/f(s, v 1))ds. Then, ψ(x) has the global minimum at x = x 1 and ψ(x 1) = 0. Furthermore, ψ(x) > 0 for x > 0. Hence, V 2 is positive definite with respect to E 1. For t ≥ 0, the time derivative of V 2 along the solutions of system (2) is

Note that λ = f(x 1, v 1)v 1 + μ 1 x 1, f(x 1, v 1)v 1 = (σ + μ 2)u 1, and μ 3 y 1 = σu 1; we have

Since the arithmetic mean is greater than or equal to the geometric mean, it has From hypotheses (H3)-(H4), we have

Note Remark 3, we have y 1 ≤ μ 5/(q − m). Therefore, if R 1 ≤ 1. Then it follows from stability theorems in [34, 35] that the infected equilibrium without immunity E 1 is stable for any time delay τ ≥ 0 if R 1 ≤ 1.

Furthermore, note that, for each t ≥ 0, implies that x(t) = x 1, u(t) = u 1, y(t) = y 1, and v(t) = v 1. Let M be the largest invariant set in the set We have from system (2) and the invariance of M that M = {E 1}. Since any solution of system (2) is bounded, it follows from LaSalle invariance principle (see, e.g., [34, 35]) that the infected equilibrium without immunity E 1 is also globally attractive for any time delay τ ≥ 0 if R 1 ≤ 1.

The characteristic equation of system (2) at E 1 takes the formwhere ψ 0(s) is a polynomial with respect to s. Let Thus we have lim ψ 1(s) > 0 and ψ 1(0) = μ 5 − (q − m)y 1. From Remark 3, we have that (q − m)y 1/μ 5 > 1 if R 1 > 1. Thus, ψ 1(0) < 0 if R 1 > 1. Hence, if R 1 > 1, then ψ 1(s) = 0 has at least a positive real root; that is, (38) has at least a positive real root. Therefore, the infected equilibrium without immunity E 1 is unstable.

4. The Local Stability of the Infected Equilibrium and Hopf Bifurcation

The characteristic equation of system (2) at the infected equilibrium with CTL response E is given bywhereWhen τ = 0, (40) becomeswhere Denote Since μ 4(σ + μ 2)(μ 3 + pz ) = σkf(x , v ) and (H4), we have σkH = μ 4 DE, μ 4 DE + Bkσ > 0, and ADEμ 4 − G > 0.

Thus, where

Assume further that

E ≥ max{μ 4, D}; that is, μ 3 + pz ≥ μ 4 and μ 3 + pz ≥ σ + μ 2.

We have Therefore, Δ4 > 0, Δ5 > 0. By Routh-Hurwitz criterion, all the roots of (42) have negative real parts. Hence we have the following result.

Proposition 6 .

When τ = 0, if , R 1 > 1, and (H8) hold, then the infected equilibrium with CTL response E is locally asymptotically stable.

In fact, when τ = 0, we can show that if and R 1 > 1 hold, the infected equilibrium with CTL response E is globally asymptotically stable by constructing suitable Lyapunov function.

Proposition 7 .

Suppose that hypotheses (H1)–(H7) and are satisfied. If R 1 > 1, then the infected equilibrium with CTL response E is globally asymptotically stable when τ = 0.

By the following Lyapunov function, V 3 is positive definite with respect to E . For t ≥ 0, the time derivative of V 3 along the solutions of system (2) is Note that λ = f(x , v )v + μ 1 x , f(x , v )v = (σ + μ 2)u , and μ 3 y = σu − py z ; we have Since the arithmetic mean is greater than or equal to the geometric mean, it has From hypotheses (H3)-(H4), we have Therefore, if R 1 > 1. Then it follows from stability theorems in [34, 35] that the infected equilibrium CTL response E is stable for τ = 0 if R 1 > 1. Similarly, by LaSalle invariance principle, we can show that the infected equilibrium CTL response E is also globally attractive for τ = 0 if R 1 > 1.

Next, we consider the case when τ > 0. Since α 5 > 0, s = 0 is not a root of (40). We suppose (40) has a purely imaginary root s = iω  (ω > 0) for some τ > 0. Substituting s = iω into (40) and separating the real and imaginary parts, we haveSquaring and adding the two equations of (53), it follows thatwhere Letting ν = ω 2, (54) can be written asThen we haveDenote

By a similar argument as that in [38], we have the following results.

Lemma 8 .

For the polynomial equation (56), the following results hold.

Equation (56) has at least one positive root, if one of the following conditions (a)–(d) holds:

C 5 < 0.

C 5 ≥ 0, q 1 = 0, Θ0 ≥ 0, and p 1 < 0 or r 1 ≤ 0 and there exists ν ∈ {ν 1, ν 2, ν 3, ν 4} such that ν > 0 and h(ν ) ≤ 0, where ν = y − (1/5)C 1  (i = 1,2, 3,4), and

C 5 ≥ 0, q 1 ≠ 0, s > p 1, Θ2 ≥ 0, or Θ3 ≥ 0 and there exists ν ∈ {ν 1 , ν 2 , ν 3 , ν 4 } such that ν > 0 and h(ν ) ≤ 0, where ν = y − (1/5)C 1  (i = 1,2, 3,4), and

, and , where .

If the conditions (a)–(d) of (i) are all not satisfied, then (56) has no positive real root.

Suppose that h(ν) = 0 has positive real roots. Without loss of generality, we may assume that (56) has positive real roots, denoted, respectively, as . Then, (54) has positive real roots . From (40), we get Therefore, letwhere . Then ±iω are a pair of purely imaginary roots of (54) with τ = τ (.

DefineLet s(τ) = ξ(τ) + iω(τ) be a root of (40) satisfying ξ(τ () = 0 and ω(τ () = ω . Differentiating the two sides of (40) with respect to τ and noticing that s is a function of τ, it follows that Thus, we get From (40), we attain ThenTherefore, it follows that Since ν > 0, we can know that Re[ds (τ)/dτ| = τ (] and h′(ν ) have the same sign.

From the above analysis, we have the following results.

Theorem 9 .

Let τ (, τ 0, and ω 0 be defined by (62) and (63). If and R 1 > 1 are satisfied, then the following results hold:

If the conditions (a)–(d) of Lemma 8 are all not satisfied, then the infected equilibrium with CTL response E is locally asymptotically stable for all time delay τ > 0.

If one of the conditions (a)–(d) of Lemma 8 is satisfied, then the infected equilibrium with CTL response E is locally asymptotically stable for τ ∈ [0, τ 0) and unstable for τ > τ 0.

If all the conditions as stated in (ii) hold and h′(ν ) ≠ 0, then system (2) undergoes a Hopf bifurcation at E when τ = τ (  (j = 0,1, 2,…).

5. Conclusion and Numerical Simulations

In this paper, we proposed a class of delayed HIV infection model (2) with general incidence rate and immune impairment. This general incidence only satisfies some general hypotheses and includes many types of special incidence functions as special cases. First, we discussed the nonnegativity and boundedness of the solutions and the existence of equilibria of system (2). Then, by constructing suitable Lyapunov functionals and using Lyapunov-LaSalle invariance principle and Hopf bifurcation theorem, we proved the following results.

If R 0 ≤ 1, the infection-free equilibrium E 0 is globally asymptotically stable for any time delay τ ≥ 0; that is, any solution (x(t), u(t), y(t), v(t), z(t)) → E 0 = (x 0, 0,0, 0,0). In biology, this means that the virus can be finally cleared from the body and the disease dies out. At the same time, as the time t increases, the numbers of latently infected cells, actively infected cells, and CTLs trends to zero and the number of uninfected cells trends to a constant x 0.

If R 0 > 1 ≥ R 1 and , the infected equilibrium without immunity E 1 is globally asymptotically stable for any time delay τ ≥ 0; that is, any solution (x(t), u(t), y(t), v(t), z(t)) → E 1 = (x 1, u 1, y 1, v 1, 0). In biology, this indicates that the HIV infection will finally become chronic with no persistent CTL response.

If R 1 > 1 and , there exists a unique infected equilibrium with CTL response E . The result of Theorem 9 implies that the time delay τ can destabilize the stability of the infected equilibrium with CTL response E and leads to the occurrence of Hopf bifurcations.

If the time delay τ ∈ [0, τ 0), the infected equilibrium with CTL response E is locally asymptotically stable. In biology, this implies that the HIV infection may become chronic and the CTL immune response may be persistent. When the time delay τ passes through the critical value τ 0, the infected equilibrium with CTL response E will become unstable and a Hopf bifurcation occurs under some conditions. In biology, this suggests that as the time delay τ increases, the numbers of the uninfected cells, latently infected cells, actively infected cells, free virus, and CTLs will first attend constant values and then become oscillated.

We now give numerical simulations to illustrate the main results in Sections 3 and 4.

Let us choose f(x, v)v = βxv. Then we have that R 0 = kλβσ/(μ 1 μ 3 μ 4(σ + μ 2)) and R 1 = R 0 − kβμ 5/(μ 1 μ 4(q − m)). Based on the numerical simulations in [12, 19, 21, 22, 27], let us take the following data: Direct calculations show that R 0 = 0.8354 < 1 and R 1 = 0.3146 < 1; system (2) has the infection-free equilibrium E 0 = (13500,0, 0,0, 0). By Theorem 4, the infection-free equilibrium E 0 is globally asymptotically stable for any time delay τ ≥ 0. Figure 1 gives the phase trajectories of system (2) with suitable initial condition.

Figure 1

Phase trajectories of system (2) with R 0 ≤ 1.

Next, let us choose the following data: Direct computations show that R 0 = 1.6544 > 1, , and R 1 = 0.8211 < 1; system (2) has the infected equilibrium without immunity E 1 = (8160,1047.0588,2.6176,13.0882,0). Therefore, by Theorem 5, the infected equilibrium without immunity E 1 is globally asymptotically stable for any time delay τ ≥ 0. Figure 2 gives the phase trajectories of system (2) with suitable initial condition.

Figure 2

Phase trajectories of system (2) with R 0 > 1 ≥ R 1.

Furthermore, let us choose the following data: Then we have that R 0 = 3.2452 > 1, , and R 1 = 2.2837 > 1 and (56) has no positive root. System (2) has the infected equilibrium with CTL response E = (6882.3529,1272.6244,3.84615,19.2308,13.0882). From Theorem 9(i), the infected equilibrium with CTL response E is locally asymptotically stable for any time delay τ > 0. Figure 3 gives the phase trajectories of system (2) with suitable initial condition.

Figure 3

Phase trajectories of system (2) with R 1 > 1, τ = 28. The initial condition is (7970, 1000, 2.6, 0.2, 17).

Finally, let us choose the following data: Then we have that R 0 = 3.2452 > 1, , and R 1 = 2.4119 > 1, (56) has two positive roots, and h′(ν ) ≠ 0. By simple computations, we have ω 0 ≈ 0.0394 and τ 0 ≈ 27.2546. From Theorem 9(ii), the infected equilibrium with CTL response E = (7363.6364,1180.0699,3.3333,16.6667,15.4021) is asymptotically stable if 0 < τ < τ 0 and unstable if τ > τ 0. Figure 4 gives the phase trajectories of system (2) with τ < τ 0 and suitable initial condition. Figure 5 gives the phase trajectories of system (2) with τ > τ 0 and suitable initial condition and shows the occurrence of the Hopf bifurcations.

Figure 4

Phase trajectories of system (2) with R 1 > 1, τ = 25 < τ 0. The initial condition is (7970, 1000, 2.6, 0.2, 17).

Figure 5

Phase trajectories of system (2) with R 1 > 1, τ = 28 > τ 0. The initial condition is (7970, 1000, 2.6, 0.2, 17).

Since v = kμ 5/(μ 4(q − m)) and z = μ 1 μ 3 μ 4(q − m)(R 1 − 1)/(pμ 1 μ 4(q − m) + pβkμ 5), it is easy to see that the number of free viruses is increased and the number of CTLs is decreased with respect to the immune impairment rate m. For example, if we choose m = 0.015, then v 1 = 25 and z 1 = 8.8462. If we choose m = 0.016, then v 2 = 27.7778 and z 2 = 7.1691. If we choose m = 0.017, then v 3 = 31.25 and z 3 = 5.3283. Figure 6 shows that, for any initial conditions, when m = 0.015,0.016,0.017, the numbers of free viruses trend to v 1 , v 2 , and v 3 , respectively. Figure 7 shows that, for any initial conditions, when m = 0.015,0.016,0.017, the numbers of CTLs trend to z 1 , z 2 , and z 3 , respectively. In Figures 6 and 7, all the data are chosen as in Figure 4 except the immune impairment rate m.

Figure 6

The curves of the free virus of system (2) with τ = 20, m = 0.015, 0.016, and 0.017. The initial conditions are chosen as (7000,1000,8, 35,10), (7000,1000,8, 20,10), and (7000,1000,8, 8,10), respectively.

Figure 7

The curves of the CTLs of system (2) with τ = 20, m = 0.015, 0.016, and 0.017. The initial conditions are chosen as (7000,1000,8, 2,8), (7000,1000,8, 2,13), and (7000,1000,8, 2,18), respectively.

As immune impairment rate m increases, the CTL response gradually becomes weak and the individuals eventually develop AIDS. Thus, in order to control the HIV infection, we should decrease the value of m. Numerical simulations show the similar known results (see, e.g., [22]).