- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Abstract and Applied Analysis

Volume 2014 (2014), Article ID 187897, 9 pages

http://dx.doi.org/10.1155/2014/187897

## Global Stability for a Viral Infection Model with Saturated Incidence Rate

^{1}School of Mathematics and Information Sciences, Guangzhou University, Guangzhou 510006, China^{2}Key Laboratory of Mathematics and Interdisciplinary Science of Guangdong, Higher Education Institutes, Guangzhou University, Guangzhou 510006, China

Received 13 January 2014; Accepted 18 February 2014; Published 6 April 2014

Academic Editor: Chuangxia Huang

Copyright © 2014 Huaqin Peng and Zhiming Guo. 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

A viral infection model with saturated incidence rate and viral infection with delay is derived and analyzed; the incidence rate is assumed to be a specific nonlinear form . The existence and uniqueness of equilibrium are proved. The basic reproductive number is given. The model is divided into two cases: with or without delay. In each case, by constructing Lyapunov functionals, necessary and sufficient conditions are given to ensure the global stability of the models.

#### 1. Introduction

In recent years, study of infectious disease model has been a hot issue; the main cause of infectious disease is the virus invasion. As we know, viral cytopathicity within target cells is very common. A number of mathematical models have been used to study virus dynamics. In 1996, Nowak et al. [1] designed a simple but natural mathematical model based on ordinary differential equation. The model is as follows: where denotes the number of uninfected cells, the numbers of infected cells, and the numbers of free viral particles at time , respectively. In model (1), uninfected target cells are assumed to be produced at a constant rate and died at rate . Infection of target cells by in-host free viruses is assumed to occur at a bilinear rate ; infected cells are lost at a rate . Free viruses are produced by infected cells at a rate , in which is the average number of viral particles produced over the lifetime of a single infected cell. Free viral particles die at a rate . For model (1), Korobeinikov [2] established the condition of global stability in 2004. Some other viral dynamical models were proposed by later researchers; see for example [3–8].

In [7], Wodarz and Levy pointed out that the term in model (1) should consist of two parts: one is the natural death of infected cells, the other is viral cytopathicity. In 2012, Li et al. [4] assumed that infected cells burst and then release viral particles (i.e., viral cytopathicity occurs) after uninfected cells were infected by a constant period of time ; that is, the time period of viral cytopathicity within target cells is . They incorporated the delay of viral cytopathicity within target cells and built a new model: By constructing Lyapunov functionals, necessary and sufficient conditions were obtained ensuring the global stability of the model.

In models (1) and (2), the researcher studied the viral dynamics with bilinear incidence rate . As we know, as the viral particles diffuse in the body, the person often takes some actions when gets large. In order to describe the inhibitory effect from the uninfected cells when the number of viral cytopathicity is large enough, following the idea of [9], we propose an incidence rate , where measures the infection force of the viral, and reflects the level of inhibitory action.

Similar to the discussions in [4], we assume that the viral cytopathicity has time delay. When the delay of viral cytopathicity within target cells is and the natural death rate of per target cell is , the number of infected cells at time can be represented by where is the probability that target cells survive from time to time , and is the number of target cells being infected at time and still surviving at time .

Differentiating of (3), we get where is the transfer rate of the infected cells being used to produce free viruses at time ; the recruitment rate of free virus at time is , in which is the average number of viral particles produced by an infected target cell when viral cytopathicity occurs, which implies that the recruitment of virus at time depends on the number of target cells that were newly infected at time and still alive at time . Therefore following the model (2), we obtain a basic viral dynamical model of delay differential equations:

Since the variable does not appear in the first and the third equations of (5), we only focus on the following equations: which has the same dynamics with system (5).

Let , by (6), we have where all the parameters are assumed to be positive.

The rest of this paper is organized as follows. In the next section we will derive the infection-free equilibrium and the infection equilibrium. In Section 3, we carry out a qualitative analysis of the model, and stability conditions for the infection-free equilibrium and the infection equilibrium are derived, respectively. A brief conclusion will be given in Section 4.

#### 2. Positive Solutions and Equilibria

Due to the biological meaning of the components , we consider system (7) with the following initial conditions: Equation (8) is a system of retarded differential equations in . is a Banach space of continuous mappings from into with norm for . We denote As usual, for any continuous function and any given , is defined as , , for any .

Theorem 1. *All the solutions of (8) under the initial conditions are positive on .*

*Proof. *Assume that there is a such that ; then by and the continuity of , there is such that for . Then we have . However, by the first equation of (8); this is a contradiction. Therefore for all .

From the second equation of (8)
Multiplying in both sides of the above equation and integrating it from to , we have
Let
Since , , and for , then for .

Further, when , we have
By the fact that and for , then
Also and for ; then
Consequently, for , which implies that holds true for .

We assume that for a positive integer , for . When , we have
Then similar discussions show that for . Hence, for all .

Theorem 2. *All the solutions of (8) under the initial conditions are ultimately bounded.*

*Proof. *For any solution of (8), define a function . Then the derivative of is
where .

Integrating both sides of inequality above from to , we have
It means that for any as long as . Also,
From the first equation of (8), we have
Similar discussion shows that
Then for any as long as . Moreover,
Thus, the region
is an invariant set and an attractor of system (8) with initial condition .

In what follows, we study the existence of equilibria.

We consider algebraic equations It is easy to see that system (24) always has an infection-free equilibrium . To find the other equilibrium, we assume . By the second equation of (24), we have

Then

We put into the first equation of (24); then Thus the positive root exists if and only if . For system (8), define the basic reproduction number [10] as follows

It is easy to see that(i)If , then system (8) has a unique equilibrium , which corresponds to the case that viruses die out, and it is called infection-free equilibrium.(ii)If , then system (8) has two equilibria, one is the infection-free equilibrium and the other is a positive equilibrium .

#### 3. Stability of the Equilibrium

In this section, we consider the stability of the equilibrium. There are two cases, and .

##### 3.1. Local Stability of Equilibria

First we consider the case of . In this case system (8) is reduced to a system of ordinary differential equations. In order to examine local stability of an equilibrium, we should compute the eigenvalues of the linearized operator for system (8) at the equilibrium.

By a direct computation, the Jacobian matrix is as follows: Consider infection-free equilibrium . The characteristic equation is obtained by the standard method as follows.

It is obvious that , and are the characteristic roots of the characteristic equation. Therefore, we have the following theorem.

Theorem 3. *(i) If , then infection-free equilibrium is locally asymptotically stable.**(ii) If , then infection-free equilibrium is unstable.**(iii) If , then infection-free equilibrium is degenerated.*

Now, local stability of the infection equilibrium is considered. As we know, infection equilibrium , exists if and only if .

Theorem 4. *If , then the infection equilibrium is locally asymptotically stable.*

*Proof. *Set . Then system (8) at the equilibrium has Jacobian matrix
A direct computation shows that the characteristic equation is
By Hurwitz criterion, all of the eigenvalues of characteristic equation have negative real parts if and only if
Indeed,
This implies that all the eigenvalues of characteristic equation have negative real parts. Then the infection equilibrium is locally asymptotically stable. This completes the proof of theorem.

Now we consider the case . By linearizing system (8) at the infection-free equilibrium , we obtain the characteristic equation as follows:

It is easy to see that ; hence we only need to discuss the roots of the following equation:

Theorem 5. *When , then*(i)*If , then the infection-free equilibrium is locally asymptotically stable.*(ii)*If , then the infection-free equilibrium is unstable.*(iii)*If , then the infection-free equilibrium is degenerated.*

*Proof. *(i) By implicit function theorem for complex variables, we know that the roots of (35) are continuous on the parameter .

If , then is not a root of (35) for all . Note that all complex roots of (35) must come in conjugate pairs and the root of (35) is negative for . Thus, all roots of (35) have negative real parts for small ; that is, . Suppose that there exists a positive number such that (35) has a pair of purely imaginary roots ; here is a positive number. We have
Then
Summing up the square of both equations in (37) we obtain
When , then . It is a contradiction with which leads to the nonexistence of . This contradiction proves the result.

(ii) When , and , then
Therefore equation must have a positive real root for all .

(iii) If , it is easy to know that is a root of (35) for all , which leads to conclusion. (iii) This completes the proof of theorem.

Now we consider the local stability of the infection equilibrium . As we know, the infection equilibrium exists if and only if . By computation, the associated transcendental characteristic equation of (8) at becomes where

Theorem 6. *When , if , then the infection equilibrium is locally asymptotically stable.*

*Proof. *By implicit function theorem for complex variables, we know that the root of (40) is continuous on the parameter . If , then all roots of (40) have negative real parts as and (40) has no zero root for all . Thus, all roots of (40) have negative real parts for very small ; that is, . Assume that there exists a positive such that (40) has a pair of purely imaginary roots , . Then must satisfy
Separating the real and imaginary parts, we have
which implies that
Direct computation shows that
Let
Note that ,

By Hurwitz criterion, (44) has no positive roots, which implies the nonexistence of . Thus all roots of (40) have negative real parts for .

##### 3.2. Global Stability of Equilibria

In the section, we study the global stability of equilibria; we first consider the infection-free equilibrium .

Theorem 7. *When ,*(i)*If , then infection-free equilibrium is globally asymptotically stable.*(ii)*If , then infection equilibrium is globally asymptotically stable.*

*Proof. *(i) Define a Lyapunov function as what follows
It means that is negative semidefinite as . Moreover, the last equality of the above equation shows that the largest invariant set of system (8) on the region is the singleton . Therefore, the infection-free equilibrium is global asymptotically stability.

(ii) We rewrite the system (8)
Choose a Dulac function
We have
Thus system (49) does not have nontrivial periodic orbits in . The conclusion follows.

Theorem 8. *When , if , then infection-free equilibrium is globally asymptotically stable.*

*Proof. *Define a functional as follows:
For any , , then
Then,
When , we have . It is easy to know, when the largest invariant set of system (8) on the region is the singleton . By Lassalle invariant principle for autonomous retarded differential equations [11], infection-free equilibrium is globally asymptotically stable. This completes the proof.

Theorem 9. *When , if , then the infection equilibrium is globally asymptotically stable.*

*Proof. *Let and ; system (8) becomes

Denote . Evaluating both sides of (55) at , we obtain
Define a Lyapunov functional as follows:
where
Thus, with equality if and only if .

For any , , then

We calculate derivatives of , , and with respect to the system (55):
We obtain
By adding and subtracting the quantity , we have
Then,
It is easy to know that , and if and only if . It follows that , and if and only if . By classical stability theory for functional differential equations, is globally asymptotically stable. This completes the proof.

#### 4. Conclusion

The viral infection model addressed in this paper has saturated incidence rate and viral infection with delay. The basic reproductive number is given. When , for the model with or without delay time, the infection-free equilibrium is globally asymptotically stable, which implies that the viral infection goes extinct eventually; when , the infection equilibrium is globally asymptotically stable, which implies that the viral infection persists in the body of the host.

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work was partially supported by National Natural Science Foundation of China (no. 11371107 and no. 11031002), Research Fund for the Doctoral Program of Higher Education of China (no. 20124410110001), and Program for Changjiang Scholars and Innovative Research Team in University (IRT1226).

#### References

- M. A. Nowak, S. Bonhoeffer, A. M. Hill, R. Boehme, H. C. Thomas, and H. Mcdade, “Viral dynamics in hepatitis B virus infection,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 93, no. 9, pp. 4398–4402, 1996. View at Publisher · View at Google Scholar · View at Scopus - A. Korobeinikov, “Global properties of basic virus dynamics models,”
*Bulletin of Mathematical Biology*, vol. 66, no. 4, pp. 879–883, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - A. M. Elaiw, “Global properties of a class of HIV models,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 4, pp. 2253–2263, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li, X. Song, and F. Gao, “Global stability of a viral infection model with two delays and two types of target cells,”
*The Journal of Applied Analysis and Computation*, vol. 2, no. 3, pp. 281–292, 2012. View at Google Scholar · View at MathSciNet - R. Qesmi, J. Wu, J. Wu, and J. M. Heffernan, “Influence of backward bifurcation in a model of hepatitis B and C viruses,”
*Mathematical Biosciences*, vol. 224, no. 2, pp. 118–125, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Sansonno, A. R. Iacobelli, V. Cornacchiulo, et al., “Detection of hepatitis C virus (HIV) proteins by immunouorescence and HCV RNA genomic sequences by nonisotopic in situ hybridization in bone marrow and peripheral blood mononnuclear cells of chronically HCV infected,”
*Clinical & Experimental Immunology*, vol. 103, pp. 414–421, 1996. View at Google Scholar - D. Wodarz and D. N. Levy, “Human immunodeficiency virus evolution towards reduced replicative fitness in vivo and the development of AIDS,”
*Proceedings of the Royal Society B: Biological Sciences*, vol. 274, no. 1624, pp. 2481–2490, 2007. View at Publisher · View at Google Scholar · View at Scopus - X. Zhou, X. Song, and X. Shi, “A differential equation model of HIV infection of
*CD*4^{+}T-cells with cure rate,”*Journal of Mathematical Analysis and Applications*, vol. 342, no. 2, pp. 1342–1355, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Capasso and G. Serio, “A generalization of the Kermack-McKendrick deterministic epidemic model,”
*Mathematical Biosciences*, vol. 42, no. 1-2, pp. 43–61, 1978. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,”
*Mathematical Biosciences*, vol. 180, pp. 29–48, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Hale,
*Theory of Functional Differential Equations*, Springer, New York, NY, USA, 2nd edition, 1977. View at MathSciNet