Computational and Mathematical Methods in Medicine
Volume 2015 (2015), Article ID 758362, 8 pages
http://dx.doi.org/10.1155/2015/758362

Research Article

Global Dynamics of a Virus Dynamical Model with Cell-to-Cell Transmission and Cure Rate

Tongqian Zhang,1,2 Xinzhu Meng,1,2 and Tonghua Zhang³

¹College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, China
²State Key Laboratory of Mining Disaster Prevention and Control Co-Founded by Shandong Province and the Ministry of Science and Technology, Shandong University of Science and Technology, Qingdao 266590, China
³Department of Mathematics, Swinburne University of Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia

Received 16 May 2015; Revised 30 June 2015; Accepted 7 July 2015

Academic Editor: Konstantin Blyuss

Copyright © 2015 Tongqian Zhang 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.

The cure effect of a virus model with both cell-to-cell transmission and cell-to-virus transmission is studied. By the method of next generation matrix, the basic reproduction number is obtained. The locally asymptotic stability of the virus-free equilibrium and the endemic equilibrium is considered by investigating the characteristic equation of the model. The globally asymptotic stability of the virus-free equilibrium is proved by constructing suitable Lyapunov function, and the sufficient condition for the globally asymptotic stability of the endemic equilibrium is obtained by constructing suitable Lyapunov function and using LaSalle invariance principal.

1. Introduction and Model Formulation

Since tobacco mosaic virus, the first virus of the world discovered by Beijerinck in 1898 [1], more and more viruses have been discovered by biologists, biomedical scientists, and medical scientists and more than 5,000 viruses have been recorded in detail [2]. However, according to a recent study, there are at least 32,0000 viruses waiting to be discovered in the spread between mammalian species. Identifying diseases caused by these viruses, especially those that can infect people, perhaps can help us to prevent epidemic disease [3]. At the early stage of the study, it is generally accepted that because of the specificity of viruses, virus can only infect certain plant or animal species; however, more and more cases associated with emerging zoonoses have appeared, and with a deeper understanding of the virus, we found that most viruses can infect humans, such as Human Immunodeficiency Virus (HIV), Prions, Influenza Virus, Rabies Virus, Ebola Virus, and Middle East Respiratory Syndrome Coronavirus (MERSV) [3–6].

Generally, the basic process of viral infection and virus replication occurs in six main steps: attachment, penetration, uncoating, replication, assembly, and release [7]. After the whole replicative cycle, free viruses begin to diffuse and infect new host cell. Therefore, investigating the processes of viral growth and destruction of host cells so as to gain the insights into the evolutionary processes of virus and cell in body is very important. To this end, mathematical models and analysis are powerful tools.

Since mathematical models and method of mathematical analysis were used to study the dynamics of the virus, lots of models have been established to explain the evolution of the uninfected target cells, infected cells, and the free virus. In these models, early works belonged to Nowak et al. [8], Nowak and May [9], Perelson and Nelson [10], and Perelson et al. [11]. The general class of models that have been studied [8–12] have a form similar to

where , , and represent the concentrations of uninfected target cells, infected cells, and virus, respectively. For explanations of other parameters we refer to literature [12]. This model describes the processes of virus invading the target cells and the release of the virus due to the infected cell apoptosis. In the model, the authors use to represent the interaction between uninfected target and virus , which obey the principle of mass action. Based on model (1), more authors used nonlinear functions to describe the rate constant characterizing infection of cells, for example, in [13], in [14, 15], in [16], and in [17], and for details of more general nonlinear incidence rate functions please see [18–20]. Notice that there exists a potentially possible cure rate of the infected cells to the susceptible host cells in the infection process of some virus, such as Hepatitis B Virus (HBV) [21–25] and HIV [26–32]; recently, Hattaf et al. [19] adopted a general nonlinear incidence rate function with the form and introduced cure rate (denoted by ) into the following model:

In model (2), satisfies the following hypotheses:

, for all , for all , , and and , for all , , and

Recently, Tian and Liu [20] improved model (2) by proposing a more general nonlinear incidence rate function with the form and investigated the following model:

In model (3), satisfies the following hypotheses:

, for all and , and , for all and , for all , , and , for all , , and and , for all , , and

However, many researches show that direct cell-to-cell spread can happen in some enveloped viruses (e.g., Human Immunodeficiency Virus type-1 (HIV-1) [33–37], Human T-Lymphotropic Virus Type-1 (HTLV-1) [38–41], Herpes Simplex Virus (HSV) [42], and Measles [43–45]). Cell-to-cell spread not only facilitates rapid viral dissemination, but may also promote immune evasion and influence disease [46]. Moreover, a recent study has shown that cell-to-cell spread of HIV-1 can reduce the sensitivity to the antiretroviral drugs by multiple infections of target cells and, as a result, the efficacy of antiretroviral therapy is reduced [47].

Motivated by the works [18–20, 48], we propose a virus dynamical model with both cell-to-virus infection and cell-to-cell transmission and cure rate as follows:

where , , and denote the number of host cells, infected cells, and free virus, respectively. And , , and are the death rates of them, respectively. Free virus is produced by infected cells at a rate . represents the regeneration rate of host cells. is the cure rate. represents the total infection rate of host cells, which is divided into two parts and The former represents the part where infected cells infect host cells by direct contact, and the latter means that host cells are infected by the free virus. For more detail, please see [48]. In the present model, we can see , for all and , and for all , , and , which do not satisfy conditions in model (2) and conditions and in model (3). For biological considerations, we will study system (4) in the closed set

The main goal of the present paper is to investigate the globally asymptotic stability of the equilibria of (4). This work is structured as follows. In Section 1, we give the motivation and study the background of the model. In Section 2, the existence of virus-free equilibrium and the endemic equilibrium is shown based on the basic reproduction number. And the local stability of the two equilibria is discussed in Section 3. We focus on the globally asymptotic stability of the two equilibria in Section 4. Finally, a brief conclusion and discussion are given in Section 5.

2. Basic Reproduction Number and Equilibria

The basic reproduction number [49, 50] of model (4) is given as

Based on the basic reproduction number , we have Theorem 1.

Theorem 1. Model (4) always has a virus-free equilibrium , where If , model (4) has a unique endemic equilibrium , where

3. Local Stability of the Two Equilibria

In this section, we shall show the local stability of equilibria and

Theorem 2. For model (4), we have the following conclusion:

(i) is locally stable if and unstable if (ii) is locally stable if

Proof. We firstly prove (i). Notice the Jacobian of model (4) evaluated is given by

Obviously, has an eigenvalue , and the other two eigenvalues and satisfy

Then, when , , and , all the eigenvalues of have negative real parts and is locally asymptotically stable. And when and , has a positive eigenvalue and is unstable.
Next, we prove (ii). The Jacobian of model (4) evaluated is

from which we have the characteristic equation

where

Obviously, And noticing that , we have

where and are used. Then, by the Routh-Hurwitz Criterion [51], we know that all the roots of (10) always have negative real parts. Thus, the epidemic equilibrium is locally asymptotically stable for

4. Global Stability of the Two Equilibria

In this section, we study the global behaviors of model (4) by constructing Lyapunov functions. Firstly, we show the global stability of

Theorem 3. If , the virus-free equilibrium is globally asymptotically stable.

Proof. Define a Lyapunov function on as follows:

here, is a constant to be determined. It follows from (4) and (13) that

Since , we have ; then, we can choose such that Hence, we have that Then, is globally asymptotically stable.

Next, we study the global stability of the endemic equilibrium

Theorem 4. If , the epidemic equilibrium is globally asymptotically stable, where

Proof. If , we define a Lyapunov function as follows:

It follows from (4) and (15) that

Since the arithmetic mean is greater than or equal to the geometric mean, it follows that

The above equality holds only for , , and Clearly, if and , then Note that can be formulated as

Since if and only if , , and , by LaSalle invariance principle [52], the equilibrium is globally asymptotically stable.

5. Conclusion and Discussion

In this paper, we considered the cure effect of a virus model with both cell-to-cell transmission and cell-to-virus transmission. By the method of next generation matrix, the basic reproduction number is obtained. Firstly the locally asymptotic stability of the virus-free equilibrium and the endemic equilibrium is considered. Then, the globally asymptotic stability of the virus-free equilibrium is proved by constructing suitable Lyapunov function, and the sufficient condition for the globally asymptotic stability of the endemic equilibrium is obtained by constructing suitable Lyapunov function and using LaSalle invariance principal. By analyzing the condition for the globally asymptotic stability of the endemic equilibrium, we have that if , from Theorem 4, the conditions can ensure the global stability of the equilibrium , While if , by the numerical simulations (see Figures 1 and 2), we find that in Theorem 4 is not necessary and can be dropped.

(a) Time series of

(b) Time series of

(c) Time series of

Figure 1: Illustration of numerical solution of system (4) with , , , , , , , and , and , , and By calculation, one gets that and ; it is easy to verify ; then, the equilibrium is globally asymptotically stable.

(a) Time series of

(b) Time series of

(c) Time series of

Figure 2: Illustration of numerical solution of system (4) with , , , , , , , and , and , , and By calculation, one gets that and ; it is easy to verify , while the equilibrium is also globally asymptotically stable.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Tongqian Zhang and Xinzhu Meng are supported by the National Natural Science Foundation of China (no. 11371230), Shandong Provincial Natural Science Foundation, China (no. ZR2012AM012, ZR2015AQ001), a Project for Higher Educational Science and Technology Program of Shandong Province of China (no. J13LI05), Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources, Shandong Province of China, and SDUST Research Fund (2014TDJH102).

References

M. W. Beijerinck, “Über ein Contagium vivum fluidum als Ursache der Fleckenkrankheit der Tabaksblätter,” in Verhandelingen der Koninklyke akademie van Wettenschappen te Amsterdam, vol. 65, pp. 1–22, 1898, (German), Translated into English: J. Johnson, Ed., Phytopathological Classics, no. 7, pp. 33–52, American Phytopathological Society, St. Paul, Minn, USA, 1942.

M. Breitbart and F. Rohwer, “Here a virus, there a virus, everywhere the same virus?” Trends in Microbiology, vol. 13, no. 6, pp. 278–284, 2005.

S. J. Anthony, J. H. Epstein, K. A. Murray et al., “A strategy to estimate unknown viral diversity in mammals,” mBio, vol. 4, no. 5, Article ID e00598-13, 2013.

World Health Organization, Zoonoses and the Human-Animal-Ecosystems Interface, WHO, 2014.

L. H. Taylor, S. M. Latham, and M. E. J. Woolhouse, “Risk factors for human disease emergence,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 356, no. 1411, pp. 983–989, 2001.

P. A. Marx, C. Apetrei, and E. Drucker, “AIDS as a zoonosis? Confusion over the origin of the virus and the origin of the epidemics,” Journal of Medical Primatology, vol. 33, no. 5-6, pp. 220–226, 2004.

L. Collier, A. Balows, and M. Sussman, Topley and Wilson's Microbiology and Microbial Infections, vol. 1 of Virology, Arnold, 9th edition, 1998, edited by: B. Mahy, L. Collier.

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.

M. A. Nowak and R. M. May, Viral Dynamics, Oxford University Press, Oxford, UK, 2000.

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. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, and D. D. Ho, “HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time,” Science, vol. 271, no. 5255, pp. 1582–1586, 1996.

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.

L. Min, Y. Su, and Y. Kuang, “Mathematical analysis of a basic virus infection model with application to HBV infection,” Rocky Mountain Journal of Mathematics, vol. 38, no. 5, pp. 1573–1585, 2008.

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.

X. Song and A. U. Neumann, “Global stability and periodic solution of the viral dynamics,” Journal of Mathematical Analysis and Applications, vol. 329, no. 1, pp. 281–297, 2007.

G. Huang, W. Ma, and Y. Takeuchi, “Global properties for virus dynamics model with Beddington-DeAngelis functional response,” Applied Mathematics Letters, vol. 22, no. 11, pp. 1690–1693, 2009.

X. Zhou and J. Cui, “Global stability of the viral dynamics with Crowley-Martin functional response,” Bulletin of the Korean Mathematical Society, vol. 48, no. 3, pp. 555–574, 2011.

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.

K. Hattaf, N. Yousfi, and A. Tridane, “Mathematical analysis of a virus dynamics model with general incidence rate and cure rate,” Nonlinear Analysis: Real World Applications, vol. 13, no. 4, pp. 1866–1872, 2012.

Y. Tian and X. Liu, “Global dynamics of a virus dynamical model with general incidence rate and cure rate,” Nonlinear Analysis: Real World Applications, vol. 16, no. 1, pp. 17–26, 2014.

L. G. Guidotti, R. Rochford, J. Chung, M. Shapiro, R. Purcell, and F. V. Chisari, “Viral clearance without destruction of infected cells during acute HBV infection,” Science, vol. 284, no. 5415, pp. 825–829, 1999.

K. Hattaf and N. Yousfi, “Hepatitis B virus infection model with logistic hepatocyte growth and cure rate,” Applied Mathematical Sciences, vol. 5, no. 45–48, pp. 2327–2335, 2011.

K. Wang, A. Fan, and A. Torres, “Global properties of an improved hepatitis B virus model,” Nonlinear Analysis: Real World Applications, vol. 11, no. 4, pp. 3131–3138, 2010.

H. Dahari, E. Shudo, R. M. Ribeiro, and A. S. Perelson, “Modeling complex decay profiles of hepatitis B virus during antiviral therapy,” Hepatology, vol. 49, no. 1, pp. 32–38, 2009.

Q. Huan, P. Ning, and W. Ding, “Global stability for a dynamic model of hepatitis B with antivirus treatment,” Journal of Applied Analysis and Computation, vol. 3, no. 1, pp. 37–50, 2013.

X. Jiang, X. Zhou, X. Shi, and X. Song, “Analysis of stability and Hopf bifurcation for a delay-differential equation model of HIV infection of CD4⁺ T-cells,” Chaos, Solitons and Fractals, vol. 38, no. 2, pp. 447–460, 2008.

X. Zhou, X. Song, and X. Shi, “A differential equation model of HIV infection of CD4⁺T-cells with cure rate,” Journal of Mathematical Analysis and Applications, vol. 342, no. 2, pp. 1342–1355, 2008.

P. K. Srivastava and P. Chandra, “Modeling the dynamics of HIV and CD4⁺ T cells during primary infection,” Nonlinear Analysis: Real World Applications, vol. 11, no. 2, pp. 612–618, 2010.

K. Hattaf and N. Yousfi, “Dynamics of HIV infection model with therapy and cure rate,” International Journal of Tomography & Statistics, vol. 16, no. 11, pp. 74–80, 2011.

X. Liu, H. Wang, Z. Hu, and W. Ma, “Global stability of an HIV pathogenesis model with cure rate,” Nonlinear Analysis: Real World Applications, vol. 12, no. 6, pp. 2947–2961, 2011.

Y. Xiao, H. Miao, S. Tang, and H. Wu, “Modeling antiretroviral drug responses for HIV-1 infected patients using differential equation models,” Advanced Drug Delivery Reviews, vol. 65, no. 7, pp. 940–953, 2013.

J. Li, X. Song, and F. Gao, “Global stability of a virus infection model with two delays and two types of target cells,” Journal of Applied Analysis and Computation, vol. 2, no. 3, pp. 281–292, 2012.

D. McDonald, L. Wu, S. M. Bohks, V. N. KewalRamani, D. Unutmaz, and T. J. Hope, “Recruitment of HIV and its receptors to dendritic cell-T cell junctions,” Science, vol. 300, no. 5623, pp. 1295–1297, 2003.

Q. Sattentau, “Avoiding the void: cell-to-cell spread of human viruses,” Nature Reviews Microbiology, vol. 6, no. 11, pp. 815–826, 2008.

B. S. Stein, S. D. Gowda, J. D. Lifson, R. C. Penhallow, K. G. Bensch, and E. G. Engleman, “pH-independent HIV entry into CD4-positive T cells via virus envelope fusion to the plasma membrane,” Cell, vol. 49, no. 5, pp. 659–668, 1987.

D. M. Phillips, “The role of cell-to-cell transmission in HIV infection,” AIDS, vol. 8, no. 6, pp. 719–731, 1994.

H. Sato, J. Orensteint, D. Dimitrov, and M. Martin, “Cell-to-cell spread of HIV-1 occurs within minutes and may not involve the participation of virus particles,” Virology, vol. 186, no. 2, pp. 712–724, 1992.

T. Igakura, J. C. Stinchcombe, P. K. C. Goon et al., “Spread of HTLV-I between lymphocytes by virus-induced polarization of the cytoskeleton,” Science, vol. 299, no. 5613, pp. 1713–1716, 2003.

A.-M. Pais-Correia, M. Sachse, S. Guadagnini et al., “Biofilm-like extracellular viral assemblies mediate HTLV-1 cell-to-cell transmission at virological synapses,” Nature Medicine, vol. 16, no. 1, pp. 83–89, 2010.

C. R. M. Bangham, “The immune control and cell-to-cell spread of human T-lymphotropic virus type 1,” Journal of General Virology, vol. 84, no. 12, pp. 3177–3189, 2003.

Y. Satou, J.-I. Yasunaga, M. Yoshida, and M. Matsuoka, “HTLV-I basic leucine zipper factor gene mRNA supports proliferation of adult T cell leukemia cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 3, pp. 720–725, 2006.

K. S. Dingwell, C. R. Brunetti, R. L. Hendricks et al., “Herpes simplex virus glycoproteins E and I facilitate cell-to-cell spread in vivo and across junctions of cultured cells,” Journal of Virology, vol. 68, no. 2, pp. 834–845, 1994.

W. P. Duprex, S. McQuaid, L. Hangartner, M. A. Billeter, and B. K. Rima, “Observation of measles virus cell-to-cell spread in astrocytoma cells by using a green fluorescent protein-expressing recombinant virus,” Journal of Virology, vol. 73, no. 11, pp. 9568–9575, 1999.

S. Krantic, C. Gimenez, and C. Rabourdin-Combe, “Cell-to-cell contact via measles virus haemagglutinin-CD46 interaction triggers CD46 downregulation,” Journal of General Virology, vol. 76, no. 11, pp. 2793–2800, 1995.

T. F. Wild, E. Malvoisin, and R. Buckland, “Measles virus: both the haemagglutinin and fusion glycoproteins are required for fusion,” Journal of General Virology, vol. 72, no. 2, pp. 439–442, 1991.

C. Jolly, “Cell-to-cell transmission of retroviruses: innate immunity and interferon-induced restriction factors,” Virology, vol. 411, no. 2, pp. 251–259, 2011.

A. Sigal, J. T. Kim, A. B. Balazs et al., “Cell-to-cell spread of HIV permits ongoing replication despite antiretroviral therapy,” Nature, vol. 477, no. 7362, pp. 95–99, 2011.

X. Lai and X. Zou, “Modeling HIV-1 virus dynamics with both virus-to-cell infection and cell-to-cell transmission,” SIAM Journal on Applied Mathematics, vol. 74, no. 3, pp. 898–917, 2014.

O. Diekmann, J. A. P. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.

P. Van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 1, pp. 29–48, 2002.

E. J. Routh, W. K. Clifford, C. Sturm, and M. Bocher, Stability of Motion, Taylor & Francis, London, UK, 1975.

J. P. La Salle, The Stability of Dynamica Systems, SIAM Publications, Philadelphia, Pa, USA, 1976.