Discrete Dynamics in Nature and Society

Volume 2015, Article ID 235420, 11 pages

http://dx.doi.org/10.1155/2015/235420

## Analysis of a Viral Infection Model with Delayed Nonlytic Immune Response

^{1}College of Science, Northwest A&F University, Yangling, Shaanxi 712100, China^{2}School of Economics and Management, Xidian University, Xi’an, Shaanxi 710071, China

Received 6 June 2014; Accepted 26 September 2014

Academic Editor: Kaifa Wang

Copyright © 2015 Mengye Chen 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 the dynamical behavior of a virus infection model with delayed nonlytic immune response. By analyzing corresponding characteristic equations, the local stabilities of two boundary equilibria are established. By using suitable Lyapunov functional and LaSalle’s invariance principle, we establish the global stability of the infection-free equilibrium. We find that the infection free equilibrium is globally asymptotically stable when , and the infected equilibrium without immunity is local asymptotically stable when . Under the condition we obtain the sufficient conditions to the local stability of the infected equilibrium with immunity . We show that the time delay can change the stability of and lead to the existence of Hopf bifurcations. The stabilities of bifurcating periodic solutions are studied and numerical simulations to our theorems are provided.

#### 1. Introduction

Mathematical models have been proven valuable in understanding the population dynamics of viral load in vivo. A proper model may play significant role in a better understanding of the disease and the various drug therapy strategies. Viral infection models have received great attention in recent years [1–7]. In most viral infections, cytotoxic T lymphocytes (CTLs) play a critical role in antiviral defense by providing a cell-mediated response to specific foreign antigens associated with cells. Therefore, the immune response after a viral infection is universal and necessary to eliminate or control the disease.

Recently, there have been a lot of papers on virus dynamics within host; some include the immune response directly [8–13]. During viral infections, the host immune system reacts with innate and antigen-specific immune response. Both types of response can be subdivided broadly into lytic and nonlytic components. Lytic components kill infected cells, whereas the nonlytic components inhibit viral replication through soluble mediators. As a part of the innate response, cytotoxic T lymphocytes (CTLs) kill infected cells, while antibodies neutralize free virus particles and thus inhibit the infection of susceptible cells. In addition, CD4+ and CD8+ T cells can secrete cytokines that inhibit viral replication [12]. In order to investigate the role of direct lytic and nonlytic inhibition of viral replication by immune cells in viral infections, a mathematical model was constructed to describe the basic dynamics of the interaction among susceptible host cells, a virus population, and immune response, which is described by the following differential equations [14, 15]: where , , and represent the densities of uninfected target cells, virus, and CTL cells at time , respectively. Uninfected cells are produced at rate , die at rate , and become infected by virus at rate without the immune response; to model nonlytic antiviral, viral replication is inhibited by the immune response at rate ; infected cells die at rate and are removed at rate by the CTL immune response. The virus-specific CTL cells proliferate at rate by contact with infected cells and die at rate . The parameter expresses the strength of the lytic component; the parameter expresses the efficacy of the nonlytic component.

Time delays cannot be ignored in models for immune response. Antigenic stimulation generating CTLs may need a period of time [16–18]; the CTL response at time may depend on the population of antigen at a previous time . Under the assumption of retarded immune response, Wang et al. [18] studied the effects of the time delay for immune response and assumed the time evolution of the population of CTL cells is governed by the delayed nonlinear differential equation . Li and Shu [19] and Xie et al. [20] investigated the effects of a time delay on a three-dimensional system with . However, both studies do not consider the nonlytic component. In this paper we propose a more general model, with the initial conditions where which is the Banach space of continuous functions mapping the interval into with the topology of uniform convergence. Clearly, the models studied in [18–20] correspond to the case in our general model (2). In this work we study (2) in the more general case .

This paper is organized as follows. In Section 2, we study the local asymptotic stability of the infection-free equilibrium and the immune-exhausted equilibrium of system (2), and the global asymptotic stability of the infection-free equilibrium also is investigated. In Section 3, we analyze the local stability of the positive equilibrium and the existence of Hopf bifurcations. In Section 4, the direction and stability of Hopf bifurcation are analyzed by the normal form theory and center manifold approach. In Section 5, we present numerical simulations to illustrate our results. Finally, in Section 6 we provide concluding remarks.

#### 2. The Stability Analysis of Equilibrium

Considering the existence of the three equilibria, then we have the following conclusions.

Proposition 1. *Let . Then, consider the following. *(1)*If , system (2) only has an infection-free equilibrium , where .*(2)*If , in addition to the infection-free equilibrium , system (2) has another unique equilibrium where and .*(3)*If , in addition to the infection-free equilibrium , system (2) has another infected equilibrium that corresponds to the survival of free virus and CTL, where
**By the similar method in [20], we have the following result.*

Theorem 2. *Under the above initial conditions (3), all solutions of system (2) are positive for and there exists , such that all the solutions satisfy , , and for all large .*

Theorem 3. *If , then the infection-free equilibrium is locally asymptotically stable, and is unstable if .*

*Proof. *The characteristic equation about is given by
that is,
It is clear that (6) has negative roots and . If , then . This shows that all roots of (6) have the negative real part; infection-free equilibrium is locally asymptotically stable. If , then . Equation (6) has a positive root; infection-free equilibrium is unstable.

Theorem 4. *If , then immune-exhausted equilibrium is locally asymptotically stable, and is unstable if .*

*Proof. *The characteristic equation about is given by
that is,
It is easy to see that
If , then and are both negative real roots, and is determined by the following equation:
If and , then . Therefore, if and , then is locally asymptotically stable. If , let () be a solution of (10). Substituting into (10) and separating the real and imaginary parts, we have
Squaring and adding the two equations of (11) together, we obtain , which has solution
Since , then we have and ; thus we get , which contradicts with our original assumption . Therefore, (10) must have no pure imaginary root. Hence, in the case of , is locally asymptotically stable for all . If , then , and . By the intermediate value theorem of continuous functions, we see that equation has at least one positive real root; thus the immune-exhausted equilibrium is unstable.

Theorem 5. *If , the infection-free equilibrium is globally asymptotically stable.*

*Proof. *Define the Lyapunov functional
we calculate the derivative of along positive solutions of system (2) then substitute for to obtain
Clearly, in the case of , we have since , , and are positive, and if and only if . Thus is globally asymptotically stable by Lyapunov-LaSalle invariance principle.

#### 3. The Stability of Infected Equilibrium and Hopf Bifurcation

Let , apply to the characteristic equation of system (2) at the positive equilibrium , and we have that is, where If , (16) becomes where Then, By the Routh-Hurwitz criterion, we know that all roots of (18) have negative real parts. From the above analysis, the following theorem holds.

Theorem 6. *If , the equilibrium is locally asymptotically stable for .*

From Theorem 6, when , all roots of (16) lie to the left of imaginary axis. But, as is increased from zero, some of the roots may cross the imaginary axis to the right. Then the equilibrium becomes unstable. Now we suppose (16) has a pure imaginary root (). Obviously, () is a root of (16) if and only if satisfies Separating the real and imaginary parts, we have From (21), we also obtain So we have where Denote Therefore, if (16) has a pure imaginary root , then equation has a positive real root .

Suppose that (27) has , , positive real roots, which are , , respectively. From (22), we have Let where and . Then is a pair of pure imaginary roots of (16).

Theorem 7 (see [21]). *Suppose the characteristic equation is the form
**
where and are continuously differentiable with respect to . One of the roots to (31) is , where is continuously differentiable with respect to , satisfying and for a positive real number . Denote ; then we have .*

Let be the root of (16) near satisfying , . Then, by Theorem 7, we have the following theorem.

Theorem 8. * and have the same sign.*

We choose , where is defined by (30). Summarizing the above analysis and combining [22], the following theorem holds.

Theorem 9. * For system (2), its equilibrium is asymptotically stable for if (27) has some positive real roots;** if is a simple positive real root of (27), then system (2) undergoes a Hopf bifurcation at the equilibrium when .*

#### 4. Direction and Stability of Hopf Bifurcation

In this section, we will study the direction of the Hopf bifurcations and stability of bifurcating periodic solutions by applying the normal theory and the center manifold theorem introduced in [22]. We always assume that system (2) undergoes Hopf bifurcation at the positive equilibrium for , and then are pure imaginary roots of the characteristic equation at the positive equilibrium .

Let , , and , ; then is the Hopf bifurcation value. The Taylor expansion of (2) at the equilibrium point is given by We denote the above system by where , and the maps , are given by By the Riesz representation theorem, there exists a function of bounded variation for such that In fact, we can choose where denote the Dirac delta function. For , let then system (33) is equivalent to where for .

For , the adjoint operator of is defined by and define a bilinear inner product where . By the discussion in Section 3 we know that are eigenvalues of . Hence, they are also eigenvalues of . We define and to be the eigenvectors of and corresponding to the eigenvalues and , respectively.

Assume that is the eigenvector of corresponding to ; then . It follows from the definition of that we have Then we obtain Similarly, let be the eigenvector of corresponding to ; by the definition of we get In order to assure , we need to determine the value of . From (40) we have Hence Using the method in [22] we obtain the following coefficients: where Then we get which determine the quantities of bifurcating periodic solutions reduced on the central manifold at the critical value . More precisely, determines the direction of Hopf bifurcation and determines the stability of bifurcating periodic solution. Summarizing the above discussion, we have the following main result.

Theorem 10. *If (), then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for (). Moreover, determines the stability of the bifurcating periodic solutions. The bifurcating periodic solutions are stable (unstable) if ().*

#### 5. Numerical Simulation

In order to demonstrate our results and find complex dynamic behavior of system (2), we provide numerical calculations for different birth rates of susceptible cells and time delays . The parameter values are chosen from literatures [14, 15, 20]. If we choose the following data set: , , , , , , , and , then we get , . Hence, by Theorem 5, is globally asymptotically stable. Figure 1 illustrates this fact.