SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Singular Boundary Value Problems for Ordinary Differential Equations.

Open Access Research Article

A Viral Infection Model with a Nonlinear Infection Rate

Yumei Yu1, JuanJ Nieto2, Angela Torres3 and Kaifa Wang4*

Author Affiliations

1 School of Science, Dalian Jiaotong University, Dalian 116028, China

2 Departamento de Análisis Matemático, Facultad de Matemáticas, Universidad de Santiago de Compostela, 15782 Santioga de compostela, Spain

3 Departamento de Psiquiatría, Radiología y Salud Pública, Facultad de Medicina, Universidad de Santiago de Compostela, 15782 Santioga de compostela, Spain

4 Department of Computers Science, Third Military Medical University, Chongqing 400038, China

For all author emails, please log on.

Boundary Value Problems 2009, 2009:958016  doi:10.1155/2009/958016


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2009/1/958016


Received:28 February 2009
Revisions received:23 April 2009
Accepted:27 May 2009
Published:29 June 2009

© 2009 The Author(s)

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.

A viral infection model with a nonlinear infection rate is constructed based on empirical evidences. Qualitative analysis shows that there is a degenerate singular infection equilibrium. Furthermore, bifurcation of cusp-type with codimension two (i.e., Bogdanov-Takens bifurcation) is confirmed under appropriate conditions. As a result, the rich dynamical behaviors indicate that the model can display an Allee effect and fluctuation effect, which are important for making strategies for controlling the invasion of virus.

1. Introduction

Mathematical models can provide insights into the dynamics of viral load in vivo. A basic viral infection model [1] has been widely used for studying the dynamics of infectious agents such as hepatitis B virus (HBV), hepatitis C virus (HCV), and human immunodeficiency virus (HIV), which has the following forms:

(11)

where susceptible cells () are produced at a constant rate , die at a density-dependent rate , and become infected with a rate ; infected cells () are produced at rate and die at a density-dependent rate ; free virus particles () are released from infected cells at the rate and die at a rate . Recently, there have been many papers on virus dynamics within-host in different aspects based on the (1.1). For example, the influences of spatial structures on virus dynamics have been considered, and the existence of traveling waves is established via the geometric singular perturbation method [2]. For more literature, we list [3, 4] and references cited therein.

Usually, there is a plausible assumption that the amount of free virus is simply proportional to the number of infected cells because the dynamics of the virus is substantially faster than that of the infected cells, . Thus, the number of infected cells can also be considered as a measure of virus load (e.g., see [57]). As a result, the model (1.1) is reduced to

(12)

As for this model, it is easy to see that the basic reproduction number of virus is given by , which describes the average number of newly infected cells generated from one infected cell at the beginning of the infectious process. Furthermore, we know that the infection-free equilibrium is globally asymptotically stable if , and so is the infection equilibrium if .

Note that both infection terms in (1.1) and (1.2) are based on the mass-action principle (Perelson and Nelson [8]); that is, the infection rate per susceptible cell and per virus is a constant . However, infection experiments of Ebert et al. [9] and McLean and Bostock [10] suggest that the infection rate of microparasitic infections is an increasing function of the parasite dose and is usually sigmoidal in shape. Thus, as Regoes et al. [11], we take the nonlinear infection rate into account by relaxing the mass-action assumption that is made in (1.2) and obtain

(13)

where the infection rate per susceptible cell, , is a sigmoidal function of the virus (parasite) concentration because the number of infected cells can also be considered as a measure of virus load (e.g., see [57]), which is represented in the following form:

(14)

Here, denotes the infectious dose at which of the susceptible cells are infected, measures the slope of the sigmoidal curve at and approximates the average number of virus that enters a single host cell at the begin stage of invasion, measures the infection force of the virus, and measures the inhibition effect from the behavioral change of the susceptible cells when their number increases or from the production of immune response which depends on the infected cells.

In fact, many investigators have introduced different functional responses into related equations for epidemiological modeling, of which we list [1217] and references cited therein. However, a few studies have considered the influences of nonlinear infection rate on virus dynamics. When the parameter , [18, 19] considered a viral mathematical model with the nonlinear infection rate and time delay. Furthermore, some different types of nonlinear functional responses, in particular of the form or Holling-type functional response, were investigated in [2023].

Note that in (1.4). To simplify the study, we fix the slope in the present paper, and system (1.3) becomes

(15)

To be concise in notations, rescale (1.5) by . For simplicity, we still use variables instead of and obtain

(16)

where . Note that is the average life time of susceptible cells and is the average life-time of infected cells. Thus, is always valid by means of biological detection. If , the virus does not kill infected cells. Therefore, the virus is non cytopathic in vivo. However, when , which means that the virus kills infected cells before its average life time, the virus is cytopathic in vivo.

The main purpose of this paper is to study the effect of the nonlinear infection rate on the dynamics of (1.6). We will perform a qualitative analysis and derive the Allee-type dynamics which result from the appearance of bistable states or saddle-node state in (1.6). The bifurcation analysis indicates that (1.6) undergoes a Bogdanov-Takens bifurcation at the degenerate singular infection equilibrium which includes a saddle-node bifurcation, a Hopf bifurcation, and a homoclinic bifurcation. Thus, the nonlinear infection rate can induce the complex dynamic behaviors in the viral infection model.

The organization of the paper is as follows. In Section 2, the qualitative analysis of system (1.6) is performed, and the stability of the equilibria is obtained. The results indicate that (1.6) can display an Allee effect. Section 3 gives the bifurcation analysis, which indicates that the dynamics of (1.6) is more complex than that of (1.1) and (1.2). Finally, a brief discussion on the direct biological implications of the results is given in Section 4.

2. Qualitative Analysis

Since we are interested in virus pathogenesis and not initial processes of infection, we assume that the initial data for the system (1.6) are such that

(21)

The objective of this section is to perform a qualitative analysis of system (1.6) and derive the Allee-type dynamics. Clearly, the solutions of system (1.6) with positive initial values are positive and bounded. Let , and note that (1.6) has one and only one infection-free equilibrium . Then by using the formula of a basic reproduction number for the compartmental models in van den Driessche and Watmough [24], we know that the basic reproduction number of virus of (1.6) is

(22)

which describes the average number of newly infected cells generated from one infected cell at the beginning of the infectious process as zero. Although it is zero, we will show that the virus can still persist in host.

We start by studying the equilibria of (1.6). Obviously, the infection-free equilibrium always exists and is a stable hyperbolic node because the corresponding characteristic equation is .

In order to find the positive (infection) equilibria, set

(23)

then we have the equation

(24)

Based on (2.4), we can obtain that

(i)there is no infection equilibria if ;

(ii)there is a unique infection equilibrium if ;

(iii)there are two infection equilibria and if .

Here,

(25)

Thus, the surface

(26)

is a Saddle-Node bifurcation surface, that is, on one side of the surface system (1.6) has not any positive equilibria; on the surface system (1.6) has only one positive equilibrium; on the other side of the surface system (1.6) has two positive equilibria. The detailed results will follow.

Next, we determine the stability of and . The Jacobian matrix at is

(27)

After some calculations, we have

(28)

Since in this case, is valid. Thus, and the equilibrium is a saddle.

The Jacobian matrix at is

(29)

By a similar argument as above, we can obtain that . Thus, the equilibrium is a node, or a focus, or a center.

For the sake of simplicity, we denote

(210)

We have the following results on the stability of .

Theorem 2.1.

Suppose that equilibrium exists; that is, . Then is always stable if . When , we have

(i) is stable if ;

(ii) is unstable if ;

(iii) is a linear center if .

Proof.

After some calculations, the matrix trace of is

(211)

and its sign is determined by

(212)

Note that

(213)

which means that is a monotone decreasing function of variable .

Clearly,

(214)

Note that implies that

(215)

Squaring (2.15) we find that

(216)

Thus,

(217)

This means that . Thus, under the condition of and the sign of , is always valid if . When , if , if , and if .

For (1.6), its asymptotic behavior is determined by the stability of if it does not have a limit cycle. Next, we begin to consider the nonexistence of limit cycle in (1.6).

Note that is a saddle and is a node, a focus, or a center. A limit cycle of (1.6) must include and does not include . Since the flow of (1.6) moves toward down on the line where and and moves towards up on the line where and , it is easy to see that any potential limit cycle of (1.6) must lie in the region where . Take a Dulac function , and denote the right-hand sides of (1.6) by and , respectively. We have

(218)

which is negative if . Hence , we can obtain the following result.

Theorem 2.2.

There is no limit cycle in (1.6) if

(219)

Note that as long as it exists. Thus, inequality (2.19) is always valid if . When , using the expression of in (2.5), we have that inequality (2.19) that is equivalent to

(220)

Indeed, since

(221)

we have (2.19) that is equivalent to

(222)

that is,

(223)

Thus,

(224)

On the other hand, squaring (2.23) we find that

(225)

which is equivalent to

(226)

The combination of (2.24) and (2.26) yields (2.20).

Furthermore,

(227)

is equivalent to , both

(228)

are equivalent to . Consequently, we have the following.

Corollary 2.3.

There is no limit cycle in (1.6) if either of the following conditions hold:

(i) and ;

(ii) and .

When , system (1.6) has a unique infection equilibrium . The Jacobian matrix at is

(229)

The determinant of is

(230)

and the trace of is

(231)

Thus, is a degenerate singular point. Since its singularity, complex dynamic behaviors may occur, which will be studied in the next section.

3. Bifurcation Analysis

In this section, the Bogdanov-Takens bifurcation (for short, BT bifurcation) of system (1.6) is studied when there is a unique degenerate infection equilibrium .

For simplicity of computation, we introduce the new time by , rewrite as , and obtain

(31)

Note that (3.1) and (1.6) are -equivalent; both systems have the same dynamics (only the time changes).

As the above mentioned, assume that

(H1)

Then (3.1) admits a unique positive equilibrium , where

(32)

In order to translate the positive equilibrium to origin, we set and obtain

(33)

Since we are interested in codimension bifurcation, we assume further that

(H2)

Then, after some transformations, we have the following result.

Theorem 3.1.

The equilibrium of (1.6) is a cusp of codimension if (H1) and (H2) hold; that is, it is a Bogdanov-Takens singularity.

Proof.

Under assumptions (H1) and (H2), it is clear that the linearized matrix of (3.3)

(34)

has two zero eigenvalues. Let . Since the parameters satisfy the assumptions (H1) and (H2), after some algebraic calculations, (3.3) is transformed into

(35)

where , , are smooth functions in variables at least of the third order. Using an affine translation to (3.5), we obtain

(36)

where , , are smooth functions in variables at least of order three. To obtain the canonical normal forms, we perform the transformation of variables by

(37)

Then, (3.6) becomes

(38)

where , , are smooth functions in at least of the third order.

Obviously,

(39)

This implies that the origin of (3.3), that is, of (1.6), is a cusp of codimension by in [25, Theorem , Section ].

In the following we will investigate the approximating BT bifurcation curves. The parameters and are chosen as bifurcation parameters. Consider the following perturbed system:

(310)

where and are positive constants while (H1) and (H2) are satisfied. That is to say,

(311)

and are in the small neighborhood of ; and are in the small neighborhood of , where

(312)

Clearly, if is the degenerate equilibrium of (1.6). Substituting into (3.10) and using Taylor expansion, we obtain

(313)

where , , , are smooth functions of and at least of order three in variables . Making the change of variables to (3.13) and noting the conditions in (3.11) and expressions in (3.12), we have

(314)

where

(315)

, , are smooth functions in variables at least of the third order, and the coefficients depend smoothly on and .

Let . Using (3.11) and (3.12), after some algebraic calculations, we obtain

(316)

where , , are smooth functions of and at least of the third order in variables ,

(317)

Let . Then (3.16) becomes

(318)

where

(319)

is smooth function in variables at least of order three, and all the coefficients depend smoothly on and .

By setting to (3.18), we obtain

(320)

where is smooth function in variables at least of the third order and

(321)

Now, introducing a new time variable to (3.20), which satisfies , and still writing as , we have

(322)

where is smooth function of , and at least of three order in variables . Setting to (3.22), we obtain

(323)

where is smooth function of , and at least of order three in variables and

(324)

If and , it is easy to obtain the following results:

(325)

By setting and , and rewriting as , we obtain

(326)

where

(327)

and is smooth function of , and at least of order three in variables .

By the theorem of Bogdanov in [26, 27] and the result of Perko in [25], we obtain the following local representations of bifurcation curves in a small neighborhood of the origin (i.e., of (1.6).

Theorem 3.2.

Let the assumptions (H1) and (H2) hold. Then (1.6) admits the following bifurcation behaviors:

(i) there is a saddle-node bifurcation curve ;

(ii) there is a Hopf bifurcation curve ;

(iii) there is a homoclinic-loop bifurcation curve HL.

Concretely, as the statement in [28, Chapter ], when , the orbital topical structure of the system (3.26) at origin (corresponding system (1.6) at ) is shown in Figure 1.

thumbnailFigure 1. The bifurcation set and the corresponding phase portraits of system (3.26) at origin.

4. Discussion

Note that most infection experiments suggest that the infection rate of microparasitic infections is an increasing function of the parasite dose, usually sigmoidal in shape. In this paper, we study a viral infection model with a type of nonlinear infection rate, which was introduced by Regoes et al. [11].

Qualitative analysis (Theorem 2.1) implies that infection equilibrium is always stable if the virus is noncytopathic, , or cytopathic in vivo but its cytopathic effect is less than or equal to an appropriate value, . When the cytopathic effect of virus is greater than the threshold value, , the stability of the infection equilibrium depends on the value of parameter , which is proportional to the birth rate of susceptible cells and is in inverse proportion to the infectious dose . The infection equilibrium is stable if and becomes unstable if . When gets to the critical value, , the infection equilibrium is a linear center, so the oscillation behaviors may occur.

If our model (1.6) does not have a limit cycle (see Theorem 2.2 and Corollary 2.3), its asymptotic behavior is determined by the stability of . When is stable, there is a region outside which positive semiorbits tend to as tends to infinity and inside which positive semi-orbits tend to as tends to infinity; that is, the virus will persist if the initial position lies in the region and disappear if the initial position lies outside this region. Thus, besides the value of parameters, the initial concentration of the virus can also affect the result of invasion. An invasion threshold may exist in these cases, which is typical for the so-called Allee effect that occurs when the abundance or frequency of a species is positively correlated with its growth rate (see [11]). Consequently, the unrescaled model (1.5) can display an Allee effect (see Figure 2), which is an infrequent phenomenon in current viral infection models though it is reasonable and important in viral infection process.

thumbnailFigure 2. Illustrations of the Allee effect for (1.5). Here, . is stable, is a saddle point, is stable. Note that SM is the stable manifolds of (solid line), UM is the unstable manifolds of (dash line), and the phase portrait of (1.6) is divided into two domains of extinction and persistence of the virus by SM.

Furthermore, when infection equilibrium becomes a degenerate singular point, we have shown that the dynamics of this model are very rich inside this region (see Theorems 3.1 and 3.2 and Figure 1). Static and dynamical bifurcations, including saddle-node bifurcation, Hopf bifurcation, homoclinic bifurcation, and bifurcation of cusp-type with codimension two (i.e., Bogdanov-Takens bifurcation), have been exhibited. Thus, besides the Allee effect, our model (1.6) shows that the viral oscillation behaviors can occur in the host based on the appropriate conditions, which was observed in chronic HBV or HCV carriers (see [2931]). These results inform that the viral infection is very complex in the development of a better understanding of diseases. According to the analysis, we find that the cytopathic effect of virus and the birth rate of susceptible cells are both significant to induce the complex and interesting phenomena, which is helpful in the development of various drug therapy strategies against viral infection.

Acknowledgments

This work is supported by the National Natural Science Fund of China (nos. 30770555 and 10571143), the Natural Science Foundation Project of CQ CSTC (2007BB5012), and the Science Fund of Third Military Medical University (06XG001).

References

  1. Nowak, MA, May, RM: Virus Dynamics,p. xii+237. Oxford University Press, Oxford, UK (2000)

  2. Wang, K, Wang, W: Propagation of HBV with spatial dependence. Mathematical Biosciences. 210(1), 78–95 (2007). PubMed Abstract | Publisher Full Text OpenURL

  3. Campos, D, Méndez, V, Fedotov, S: The effects of distributed life cycles on the dynamics of viral infections. Journal of Theoretical Biology. 254(2), 430–438 (2008). PubMed Abstract | Publisher Full Text OpenURL

  4. Srivastava, PKr, Chandra, P: Modeling the dynamics of HIV and T cells during primary infection. Nonlinear Analysis: Real World Applications. In press

  5. Bartholdy, C, Christensen, JP, Wodarz, D, Thomsen, AR: Persistent virus infection despite chronic cytotoxic T-lymphocyte activation in gamma interferon-deficient mice infected with lymphocytic choriomeningitis virus. Journal of Virology. 74(22), 10304–10311 (2000). PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Bonhoeffer, S, Coffin, JM, Nowak, MA: Human immunodeficiency virus drug therapy and virus load. Journal of Virology. 71(4), 3275–3278 (1997). PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Wodarz, D, Christensen, JP, Thomsen, AR: The importance of lytic and nonlytic immune responses in viral infections. Trends in Immunology. 23(4), 194–200 (2002). PubMed Abstract | Publisher Full Text OpenURL

  8. Perelson, AS, Nelson, PW: Mathematical analysis of HIV-1 dynamics in vivo. SIAM Review. 41(1), 3–44 (1999). Publisher Full Text OpenURL

  9. Ebert, D, Zschokke-Rohringer, CD, Carius, HJ: Dose effects and density-dependent regulation of two microparasites of Daphnia magna. Oecologia. 122(2), 200–209 (2000). Publisher Full Text OpenURL

  10. McLean, AR, Bostock, CJ: Scrapie infections initiated at varying doses: an analysis of 117 titration experiments. Philosophical Transactions of the Royal Society B. 355(1400), 1043–1050 (2000). Publisher Full Text OpenURL

  11. Regoes, RR, Ebert, D, Bonhoeffer, S: Dose-dependent infection rates of parasites produce the Allee effect in epidemiology. Proceedings of the Royal Society B. 269(1488), 271–279 (2002). PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Gao, S, Chen, L, Nieto, JJ, Torres, A: Analysis of a delayed epidemic model with pulse vaccination and saturation incidence. Vaccine. 24(35-36), 6037–6045 (2006). PubMed Abstract | Publisher Full Text OpenURL

  13. Ruan, S, Wang, W: Dynamical behavior of an epidemic model with a nonlinear incidence rate. Journal of Differential Equations. 188(1), 135–163 (2003). Publisher Full Text OpenURL

  14. Ruan, S, Xiao, D: Global analysis in a predator-prey system with nonmonotonic functional response. SIAM Journal on Applied Mathematics. 61(4), 1445–1472 (2001). Publisher Full Text OpenURL

  15. Sharomi, O, Gumel, AB: Re-infection-induced backward bifurcation in the transmission dynamics of Chlamydia trachomatis. Journal of Mathematical Analysis and Applications. 356(1), 96–118 (2009). Publisher Full Text OpenURL

  16. Wang, W: Epidemic models with nonlinear infection forces. Mathematical Biosciences and Engineering. 3(1), 267–279 (2006). PubMed Abstract OpenURL

  17. Zhang, H, Chen, L, Nieto, JJ: A delayed epidemic model with stage-structure and pulses for pest management strategy. Nonlinear Analysis: Real World Applications. 9(4), 1714–1726 (2008). Publisher Full Text OpenURL

  18. Li, D, Ma, W: Asymptotic properties of a HIV-1 infection model with time delay. Journal of Mathematical Analysis and Applications. 335(1), 683–691 (2007). Publisher Full Text OpenURL

  19. Song, X, Neumann, AU: Global stability and periodic solution of the viral dynamics. Journal of Mathematical Analysis and Applications. 329(1), 281–297 (2007). Publisher Full Text OpenURL

  20. Cai, L, Wu, J: Analysis of an HIV/AIDS treatment model with a nonlinear incidence. Chaos, Solitons & Fractals. 41(1), 175–182 (2009). PubMed Abstract | Publisher Full Text OpenURL

  21. Wang, W, Shen, J, Nieto, JJ: Permanence and periodic solution of predator-prey system with Holling type functional response and impulses. Discrete Dynamics in Nature and Society. 2007, (2007)

  22. Wang, X, Song, X: Global stability and periodic solution of a model for HIV infection of T cells. Applied Mathematics and Computation. 189(2), 1331–1340 (2007). Publisher Full Text OpenURL

  23. Yang, J: Dynamics behaviors of a discrete ratio-dependent predator-prey system with Holling type III functional response and feedback controls. Discrete Dynamics in Nature and Society. 2008, (2008)

  24. van den Driessche, P, Watmough, J: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical Biosciences. 180, 29–48 (2002). PubMed Abstract | Publisher Full Text OpenURL

  25. Perko, L: Differential Equations and Dynamical Systems, Texts in Applied Mathematics,p. xiv+519. Springer, New York, NY, USA (1996)

  26. Bogdanov, R: Bifurcations of a limit cycle for a family of vector fields on the plan. Selecta Mathematica Sovietica. 1, 373–388 (1981)

  27. Bogdanov, R: Versal deformations of a singular point on the plan in the case of zero eigenvalues. Selecta Mathematica Sovietica. 1, 389–421 (1981)

  28. Zhang, Z, Li, C, Zheng, Z, Li, W: The Base of Bifurcation Theory about Vector Fields, Higher Education Press, Beijing, China (1997)

  29. Chun, YK, Kim, JY, Woo, HJ: No significant correlation exists between core promoter mutations, viral replication, and liver damage in chronic hepatitis B infection. Hepatology. 32(5), 1154–1162 (2000). PubMed Abstract | Publisher Full Text OpenURL

  30. Deng, G-H, Wang, Z-L, Wang, Y-M, Wang, K-F, Fan, Y: Dynamic determination and analysis of serum virus load in patients with chronic HBV infection. World Chinese Journal of Digestology. 12(4), 862–865 (2004)

  31. Pontisso, P, Bellati, G, Brunetto, M: Hepatitis C virus RNA profiles in chronically infected individuals: do they relate to disease activity?. Hepatology. 29(2), 585–589 (1999). PubMed Abstract | Publisher Full Text OpenURL