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 cusptype with codimension two (i.e., BogdanovTakens 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:
where susceptible cells () are produced at a constant rate , die at a densitydependent rate , and become infected with a rate ; infected cells () are produced at rate and die at a densitydependent 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 withinhost 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 [5–7]). As a result, the model (1.1) is reduced to
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 infectionfree 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 massaction 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 massaction assumption that is made in (1.2) and obtain
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 [5–7]), which is represented in the following form:
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 [12–17] 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 Hollingtype functional response, were investigated in [20–23].
Note that in (1.4). To simplify the study, we fix the slope in the present paper, and system (1.3) becomes
To be concise in notations, rescale (1.5) by . For simplicity, we still use variables instead of and obtain
where . Note that is the average life time of susceptible cells and is the average lifetime 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 Alleetype dynamics which result from the appearance of bistable states or saddlenode state in (1.6). The bifurcation analysis indicates that (1.6) undergoes a BogdanovTakens bifurcation at the degenerate singular infection equilibrium which includes a saddlenode 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
The objective of this section is to perform a qualitative analysis of system (1.6) and derive the Alleetype 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 infectionfree 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
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 infectionfree equilibrium always exists and is a stable hyperbolic node because the corresponding characteristic equation is .
In order to find the positive (infection) equilibria, set
then we have the equation
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,
Thus, the surface
is a SaddleNode 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
After some calculations, we have
Since in this case, is valid. Thus, and the equilibrium is a saddle.
The Jacobian matrix at is
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
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
and its sign is determined by
Note that
which means that is a monotone decreasing function of variable .
Clearly,
Note that implies that
Squaring (2.15) we find that
Thus,
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 righthand sides of (1.6) by and , respectively. We have
which is negative if . Hence , we can obtain the following result.
Theorem 2.2.
There is no limit cycle in (1.6) if
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
Indeed, since
we have (2.19) that is equivalent to
that is,
Thus,
On the other hand, squaring (2.23) we find that
which is equivalent to
The combination of (2.24) and (2.26) yields (2.20).
Furthermore,
is equivalent to , both
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
The determinant of is
and the trace of is
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 BogdanovTakens 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
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
In order to translate the positive equilibrium to origin, we set and obtain
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 BogdanovTakens singularity.
Proof.
Under assumptions (H1) and (H2), it is clear that the linearized matrix of (3.3)
has two zero eigenvalues. Let . Since the parameters satisfy the assumptions (H1) and (H2), after some algebraic calculations, (3.3) is transformed into
where , , are smooth functions in variables at least of the third order. Using an affine translation to (3.5), we obtain
where , , are smooth functions in variables at least of order three. To obtain the canonical normal forms, we perform the transformation of variables by
Then, (3.6) becomes
where , , are smooth functions in at least of the third order.
Obviously,
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:
where and are positive constants while (H1) and (H2) are satisfied. That is to say,
and are in the small neighborhood of ; and are in the small neighborhood of , where
Clearly, if is the degenerate equilibrium of (1.6). Substituting into (3.10) and using Taylor expansion, we obtain
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
where
, , 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
where , , are smooth functions of and at least of the third order in variables ,
Let . Then (3.16) becomes
where
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
where is smooth function in variables at least of the third order and
Now, introducing a new time variable to (3.20), which satisfies , and still writing as , we have
where is smooth function of , and at least of three order in variables . Setting to (3.22), we obtain
where is smooth function of , and at least of order three in variables and
If and , it is easy to obtain the following results:
By setting and , and rewriting as , we obtain
where
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 saddlenode bifurcation curve ;
(ii) there is a Hopf bifurcation curve ;
(iii) there is a homoclinicloop 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.
Figure 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 semiorbits 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 socalled 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.
Figure 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 saddlenode bifurcation, Hopf bifurcation, homoclinic bifurcation, and bifurcation of cusptype with codimension two (i.e., BogdanovTakens 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 [29–31]). 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

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

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

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

Srivastava, PKr, Chandra, P: Modeling the dynamics of HIV and T cells during primary infection.

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

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

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

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

Ebert, D, ZschokkeRohringer, CD, Carius, HJ: Dose effects and densitydependent regulation of two microparasites of Daphnia magna. Oecologia. 122(2), 200–209 (2000). Publisher Full Text

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

Regoes, RR, Ebert, D, Bonhoeffer, S: Dosedependent 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

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

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

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

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

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

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

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

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

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

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

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

Yang, J: Dynamics behaviors of a discrete ratiodependent predatorprey system with Holling type III functional response and feedback controls. Discrete Dynamics in Nature and Society. 2008, (2008)

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

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

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

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

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

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

Deng, GH, Wang, ZL, Wang, YM, Wang, KF, 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)

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