Discrete Dynamics in Nature and SocietyVolume 2011 (2011), Article ID 201274, 19 pagesdoi:10.1155/2011/201274

Research Article

Global Properties of Virus Dynamics Models with Multitarget Cells and Discrete-Time Delays

A. M. Elaiw1,2 and M. A. Alghamdi¹

¹Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia²Department of Mathematics, Faculty of Science, Al-Azhar University, Assiut, Egypt

Received 8 July 2011; Accepted 16 October 2011

Academic Editor: Yong Zhou

Copyright © 2011 A. M. Elaiw and M. A. Alghamdi. 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 propose a class of virus dynamics models with multitarget cells and multiple intracellular delays and study their global properties. The first model is a 5-dimensional system of nonlinear delay differential equations (DDEs) that describes the interaction of the virus with two classes of target cells. The second model is a ()-dimensional system of nonlinear DDEs that describes the dynamics of the virus, classes of uninfected target cells, and classes of infected target cells. The third model generalizes the second one by assuming that the incidence rate of infection is given by saturation functional response. Two types of discrete time delays are incorporated into these models to describe (i) the latent period between the time the target cell is contacted by the virus particle and the time the virus enters the cell, (ii) the latent period between the time the virus has penetrated into a cell and the time of the emission of infectious (mature) virus particles. Lyapunov functionals are constructed to establish the global asymptotic stability of the uninfected and infected steady states of these models. We have proven that if the basic reproduction number is less than unity, then the uninfected steady state is globally asymptotically stable, and if (or if the infected steady state exists), then the infected steady state is globally asymptotically stable.

1. Introduction

Nowadays, various types of viruses infect the human body and cause serious and dangerous diseases. Mathematical modeling and model analysis of virus dynamics have attracted the interests of mathematicians during the recent years, due to their importance in understanding the associated characteristics of the virus dynamics and guiding in developing efficient antiviral drug therapies. Several mathematical models have been proposed in the literature to describe the interaction of the virus with the target cells [1]. Some of these models are given by a system of nonlinear ordinary differential equations (ODEs). Others are given by a system of nonlinear delay differential equations (DDEs) to account the intracellular time delays. The basic virus dynamics model with intracellular discrete time delay has been proposed in [2] and given by where , and represent the populations of uninfected target cells, infected cells, and free virus particles at time , respectively. Here, represents the rate of which new target cells are generated from sources within the body, is the death rate constant, and is the infection rate constant. Equation (1.2) describes the population dynamics of the infected cells and shows that they die with rate constant . The virus particles are produced by the infected cells with rate constant , and are removed from the system with rate constant . The parameter accounts for the time between viral entry into the target cell and the production of new virus particles. The recruitment of virus-producing cells at time is given by the number of cells that were newly infected cells at time and are still alive at time . The probability of surviving the time period from to is , where is the constant death rate of infected cells but not yet virus-producing cells.

A great effort has been made in developing various mathematical models of viral infections with discrete or distributed delays and studying their basic and global properties, such as positive invariance properties, boundedness of the model solutions and stability analysis [3–19]. In [20–24], multiple inracellular delays have been incroporated into the virus dynamics model. Most of the existing models are based on the assumption that the virus attacks one class of target cells (e.g., CD4⁺ T cells in case of HIV or hepatic cells in case of HCV and HBV). Since the interactions of some types of viruses inside the human body is not very clear and complicated, therefore, the virus may attack more than one class of target cells. Hence, virus dynamics models describing the interaction of the virus with more than one class of target cells are needed. In case of HIV infection, Perelson et al. [25] observed that the HIV attack two classes of target cells, CD4⁺ T cells and macrophages. In [26, 27], an HIV model with two target cells has been proposed. In very recent works [28–30], we have proposed several HIV models with two target cells and investigated the global asymptotic stability of their steady states. In [31], we have proposed a class of virus dynamics models with multitarget cells. However, the intracellular time delay has been neglected in [26–31].

The purpose of this paper is to propose a class of virus dynamics models with multitarget cells and establish the global stability of their steady states. The first model considers the interaction of the virus with two classes of target cells. In the second model, we assume that the virus attacks classes of target cells. The third model generalizes the second one by assuming that the infection rate is given by saturation functional response. We incorporate two types of discrete time delays into these models describing (i) the time between the target cell is contacted by the virus particle and the contacting virus enters the cell, (ii) the time between the virus has penetrated into a cell and the emission of infectious (mature) virus particles. The global stability of these models is established using Lyapunov functionals, which are similar in nature to those used in [11, 20]. We prove that the global dynamics of these models are determined by the basic reproduction number . If , then the uninfected steady state is globally asymptotically stable (GAS). If (or if the infected steady state exists), then the infected steady state is GAS for all time delays.

2. Virus Dynamics Model with Two Target Cells and Delays

In this section, we introduce a mathematical model of virus infection with two classes of target cells. This model can describe the HIV dynamics with two classes of target cells, CD4⁺ T cells and macrophages [26, 27]. This model can be considered as an extension of the models given in [11, 26, 27]: where and represent the populations of the two classes of uninfected target cells; and are the populations of the infected cells. The population of the target cells are described by (2.1) and (2.3), where and represent the rates of which new target cells are generated, and are the death rate constants, and and are the infection rate constants. Equations (2.2) and (2.4) describe the population dynamics of the two classes of infected cells and show that they die with rate constants and . The virus particles are produced by the two classes of infected cells with rate constants and and are cleared with rate constant . Here the parameter accounts for the time between the target cells of class are contacted by the virus particle and the contacting virus enters the cells. The recruitment of virus-producing cells at time is given by the number of cells that were newly infected cells at times and are still alive at time . Also, is assumed to be a constant death rate for infected target cells, but not yet virus-producing cells. Thus, the probability of surviving the time period from to is , . The time between the virus has penetrated into a target cell of class and the emission of infectious (matures) virus particles is represented by . The probability of survival of an immature virus is given by , where is constant.

2.1. Initial Conditions

The initial conditions for system (2.1)–(2.5) take the form where , the Banach space of continuous functions mapping the interval into , where .

By the fundamental theory of functional differential equations [32], system (2.1)–(2.5) has a unique solution() satisfying the initial conditions (2.6).

2.2. Nonnegativity and Boundedness of Solutions

In the following, we establish the nonnegativity and boundedness of solutions of (2.1)–(2.5) with initial conditions (2.6).

Proposition 2.1. Let be any solution of (2.1)–(2.5) satisfying the initial conditions (2.6), then , , , and are all nonnegative for and ultimately bounded.

Proof. From (2.1) and (2.3), we have which indicates that , for all . Now from (2.2), (2.4), and (2.5), we have confiming that , , for all . By a recursive argument, we obtain , , for all .To show the boundedness of the solutions, we let and , then where and . Hence, , and , where and . On the other hand, then , where . It follows that the solution is ultimately bounded.

2.3. Steady States

It will be explained in the following that the global behavior of model (2.1)–(2.5) crucially depends on the basic reproduction number given by where and . We observe that can be written as where are the basic reproduction numbers of each class of target cell dynamics separately (see [28]).

Following the same line as in [28], we can show that if , then system (2.1)–(2.5) has only one steady state which is called uninfected steady state, and if , then system (2.1)–(2.5) has two steady states and infected steady state . The coordinates of the infected steady state are given by where

2.4. Global Stability

In this section, we prove the global stability of the uninfected and infected steady states of system (2.1)–(2.5). The strategy of the proof is to use suitable Lyapunov functionals which are similar in nature to those used in [11, 20]. Next we will use the following notation: , for any . We also define a function as It is clear that for any and has the global minimum .

Theorem 2.2. (i) If , then is GAS for any . (ii) If , then is GAS for any .

Proof. (i) We consider a Lyapunov functional where where . We note that is defined, continuous, and positive definite for all . Also, the global minimum occurs at the uninfected steady state .The time derivatives of and are given by Similarly, is given by It follows that Since the arithmetical mean is greater than or equal to the geometrical mean, then the first two terms of (2.21) are less than or equal to zero. Therefore, if , then for all . By Theorem  5.3.1 in [32], the solutions of system (2.1)–(2.5) limit to , the largest invariant subset of . Clearly, it follows from (2.21) that if and only if , , . Noting that is invariant, for each element of we have , . From (2.5) we drive that Since and for all , then if and only if . Hence, if and only if , , . From LaSalle’s invariance principle, is GAS for any .(ii) Define a Lyapunov functional as Differentiating with respect to time yields Using the infected steady state conditions we obtain From (2.25), we see that the last term in (2.26) vanishes. Then, using the following equalities: we can rewrite (2.26) as Since the arithmetical mean is greater than or equal to the geometrical mean, then the first two terms of (2.28) are less than or equal to zero. It is easy to see that if , then . By [32,Theorem  5.3.1], the solutions of system (2.1)–(2.5) limit to , the largest invariant subset of . It can be seen that if and only if , and , that is, If , then from (2.29), we have and , and hence equal to zero at . LaSalle’s invariance principle implies global stability of .

3. Basic Virus Dynamics Model with Multitarget Cells and Delays

In this section, we propose a virus dynamics model which describes the interaction of the virus with classes of target cells. Two types of discrete-time delays are incorporated into the model. The model is a generalization of those of one class of target cells and two classes of target cells models presented, respectively, in [26, 33]. Moreover, it can be seen that when and , then the following model leads to the model presented in [11]. where and represent the populations of the uninfected target cells and infected cells of class , respectively, is the population of the virus particles. All the parameters of the model have the same biological meaning as given in the previous section.

The initial conditions for system (3.1) take the form where and is the Banach space of continuous functions mapping the interval into .

Similar to the previous section, the nonnegativity and the boundedness of the solutions of system (3.1) can be shown.

3.1. Steady States

It is clear that system (3.1) has an uninfected steady state , where , , , and . The system can also have a positive infected steady state . The coordinates of the infected steady state, if they exist, satisfy the equalities: The basic reproduction number of system (3.1) is given by where is the basic reproduction number for the dynamics of the interaction of the virus only with the target cells of class .

3.2. Global Stability

In the following theorem, the global stability of the uninfected and infected steady states of system (3.1) will be established.

Theorem 3.1. (i) If , then is GAS for any , . (ii) If exists, then it is GAS for any , .

Proof. (i) Define a Lyapunov functional as follows: where . The time derivative of along the solution of system (3.1) satisfies Since the arithmetical mean is greater than or equal to the geometrical mean, then the first term of (3.8) is less than or equal to zero. Therefore, if , then for all . Similar to the previous section, one can show that the maximal compact invariant set in is the singleton when . The global stability of follows from LaSalle’s invariance principle.To prove (ii), we consider the Lyapunov functional Differentiating with respect to time yields Using the infected steady state conditions (3.3)–(3.5), and the following equality: we obtain From (3.5), we can see that the last term in (3.12) vanishes. Then, by using the following equality: we can rewrite (3.12) as It is easy to see that if , , then for all (the arithmetical mean is greater than or equal to the geometrical mean and ). Clearly, the singleton is the only invariant set in . LaSalle's invariance principle implies global stability of .

4. Virus Dynamics Model with Saturation Infection Rate

In this section, we proposed a virus dynamics model which describes the interaction of the virus with classes of target cells taking into account the saturation infection rate and multiple intracellular delays: where , are positive constants. The variables and parameters of the model have the same definitions as given in Section 2. We mention that if and , then model (4.1) leads to the model presented in [9], and if and , , then model (4.1) leads to the model presented in [34].

4.1. Steady States

It is clear that system (4.1) has an uninfected steady state , where , , and . The system can also have a positive infected steady state . The coordinates of the infected steady state, if they exist, satisfy the equalities: The basic reproduction number for system (4.1) is the same as given by (3.6).

4.2. Global Stability

In this section, we study the global stability of the uninfected and infected steady states of system (4.1).

Theorem 4.1. (i) If , then is GAS for any , . (ii) If exists, then it is GAS for any , .

Proof. (i) Define a Lyapunov functional as follows: The time derivative of along the trajectories of (4.1) satisfies It is clear that if , then for all , where equality occurs at . The global stability of follows from LaSalle’s invariance principle.To prove (ii), we consider the Lyapunov functional: Differentiating with respect to time yields Using the infected steady state conditions (4.2), and the following equality: we obtain Then using the following equalities: we obtain It is easy to see that if , then for all . Clearly, the singleton is the only invariant set in . LaSalle’s invariance principle implies global stability of .

5. Conclusion

In this paper, we have studied the global properties of a class of virus dynamics models with multitarget cells and multiple delays. First, we have introduced a model with two classes of target cells (CD4⁺ T and macrophages in case of HIV). Then, we have proposed a model describing the interaction of the virus with classes of target cells. A model with multitarget cells taking into account the saturation infection rate is also studied. Two types of discrete time delays have been incorporated into these models to take into account (i) the latent period between the time the target cell is contacted by the virus particle and the time virus enters the cell, (ii) the latent period between the time the virus has penetrated into a cell and the time of the emission of infectious (mature) virus particles. The global stability of the uninfected and infected steady states has been established by using suitable Lyapunov functionals and LaSalle invariant principle. We have proven that, if the basic reproduction number is less than unity, then the uninfected steady state is GAS and if (or the infected steady state exists) then the infected steady state is GAS for all time delays.

References

M. A. Nowak and R. M. May, Virus Dynamics: Mathematical Principles of Immunology and Virolog, Oxford University Press, Oxford, UK, 2000.

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.

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.

N. M. Dixit and A. S. Perelson, “Complex patterns of viral load decay under antiretroviral therapy: influence of pharmacokinetics and intracellular delay,” Journal of Theoretical Biology, vol. 226, no. 1, pp. 95–109, 2004.

J. E. Mittler, B. Sulzer, A. U. Neumann, and A. S. Perelson, “Influence of delayed viral production on viral dynamics in HIV-1 infected patients,” Mathematical Biosciences, vol. 152, no. 2, pp. 143–163, 1998.

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.

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.

D. Li and W. Ma, “Asymptotic properties of a HIV-1 infection model with time delay,” Journal of Mathematical Analysis and Applications, vol. 335, no. 1, pp. 683–691, 2007.

R. Xu, “Global stability of an HIV-1 infection model with saturation infection and intracellular delay,” Journal of Mathematical Analysis and Applications, vol. 375, no. 1, pp. 75–81, 2011.

Y. Nakata, “Global dynamics of a viral infection model with a latent period and Beddington-DeAngelis response,” Nonlinear Analysis. Theory, Methods & Applications, vol. 74, no. 9, pp. 2929–2940, 2011.

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.

G. Huang, W. Ma, and Y. Takeuchi, “Global analysis for delay virus dynamics model with Beddington-DeAngelis functional response,” Applied Mathematics Letters, vol. 24, no. 7, pp. 1199–1203, 2011.

Y. Wang, Y. Zhou, J. Wu, and J. Heffernan, “Oscillatory viral dynamics in a delayed HIV pathogenesis model,” Mathematical Biosciences, vol. 219, no. 2, pp. 104–112, 2009.

X. Wang, Y. Tao, and X. Song, “A delayed HIV-1 infection model with Beddington-DeAngelis functional response,” Nonlinear Dynamics, vol. 62, no. 1-2, pp. 67–72, 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.

X. Song, S. Wang, and J. Dong, “Stability properties and Hopf bifurcation of a delayed viral infection model with lytic immune response,” Journal of Mathematical Analysis and Applications, vol. 373, no. 2, pp. 345–355, 2011.

X. Shi, X. Zhou, and X. Song, “Dynamical behavior of a delay virus dynamics model with CTL immune response,” Nonlinear Analysis. Real World Applications, vol. 11, no. 3, pp. 1795–1809, 2010.

Z. Hu, X. Liu, H. Wang, and W. Ma, “Analysis of the dynamics of a delayed HIV pathogenesis model,” Journal of Computational and Applied Mathematics, vol. 234, no. 2, pp. 461–476, 2010.

Q. Xie, D. Huang, S. Zhang, and J. Cao, “Analysis of a viral infection model with delayed immune response,” Applied Mathematical Modelling, vol. 34, no. 9, pp. 2388–2395, 2010.

G. Huang, Y. Takeuchi, and W. Ma, “Lyapunov functionals for delay differential equations model of viral infections,” SIAM Journal on Applied Mathematics, vol. 70, no. 7, pp. 2693–2708, 2010.

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.

R. Ouifki and G. Witten, “Stability analysis of a model for HIV infection with RTI and three intracellular delays,” BioSystems, vol. 95, no. 1, pp. 1–6, 2009.

Y. Nakata, “Global dynamics of a cell mediated immunity in viral infection models with distributed delays,” Journal of Mathematical Analysis and Applications, vol. 375, no. 1, pp. 14–27, 2011.

R. Xu, “Global dynamics of an HIV-1 infection model with distributed intracellular delays,” Computers & Mathematics with Applications, vol. 61, no. 9, pp. 2799–2805, 2011.

A. S. Perelson, P. Essunger, Y. Cao et al., “Decay characteristics of HIV-1-infected compartments during combination therapy,” Nature, vol. 387, no. 6629, pp. 188–191, 1997.

D. S. Callaway and A. S. Perelson, “HIV-1 infection and low steady state viral loads,” Bulletin of Mathematical Biology, vol. 64, no. 1, pp. 29–64, 2002.

A. S. Perelson and P. W. Nelson, “Mathematical analysis of HIV-1 dynamics in vivo,” SIAM Review, vol. 41, no. 1, pp. 3–44, 1999.

A. M. Elaiw, “Global properties of a class of HIV models,” Nonlinear Analysis. Real World Applications, vol. 11, no. 4, pp. 2253–2263, 2010.

A. M. Elaiw, I. A. Hassanien, and S. A. Azoz, “Global properties of a class of HIV infection models with Beddington-DeAngelis functional response,” Submitted to Mathematical Methods in the Applied Sciences

A. M. Elaiw and X. Xia, “HIV dynamics: analysis and robust multirate MPC-based treatment schedules,” Journal of Mathematical Analysis and Applications, vol. 359, no. 1, pp. 285–301, 2009.

A. M. Elaiw, “Global properties of a class of virus infection models with multitarget cells,” Nonlinear Dynamics, In Press

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

M. A. Nowak, R. Anderson, M. Boerlijst, S. Bonhoeffer, R. May, and A. McMichael, “HIV-1 evolution and disease progression,” Science, vol. 274, pp. 1008–1010, 1996.

A. M. Elaiw, I. A. Hassanien, and S. A. Azoz, “Global stability of HIV infection models with intracellular delays,” Journal of the Korean Mathematical Society