BioMed Research International
Volume 2017 (2017), Article ID 7805675, 19 pages
https://doi.org/10.1155/2017/7805675

Research Article

Analysis of a New Delayed HBV Model with Exposed State and Immune Response to Infected Cells and Viruses

Deshun Sun and Fei Liu

Control and Simulation Center, Harbin Institute of Technology, Harbin 15008, China

Correspondence should be addressed to Deshun Sun; sun_deshun@hit.edu.cn and Fei Liu; liufei@hit.edu.cn

Received 16 May 2017; Revised 2 September 2017; Accepted 14 September 2017; Published 16 November 2017

Academic Editor: Lucia Lopalco

Copyright © 2017 Deshun Sun and Fei Liu. 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 comprehensive delayed HBV model, which not only considers the immune response to both infected cells and viruses and a time delay for the immune system to clear viruses but also incorporates an exposed state and the proliferation of hepatocytes. We prove the positivity and boundedness of solutions and analyze the global stability of two boundary equilibria and then study the local asymptotic stability and Hopf bifurcation of the positive (infection) equilibrium and also the stability of the bifurcating periodic solutions. Moreover, we illustrate how the factors such as the time delay, the immune response to infected cells and viruses, and the proliferation of hepatocytes affect the dynamics of the model by numerical simulation.

1. Introduction

Hepatitis B virus (HBV) has become one of the serious infectious diseases threatening global human health, which can cause chronic liver infection and further result in liver inflammation, fibrosis, cirrhosis, or even cancer [1]. Each year more than 1 million people die of end-stage liver diseases like cancer due to the HBV infection [2].

Mathematical modeling and analysis of the dynamics of such infectious viruses as HBV play important roles in understanding the factors that govern the infectious disease progression and offering insights into developing treatment strategies and guiding antiviral drug therapies [3]. So far, there have been plenty of mathematical models proposed to describe and analyze virus infection, immune responses, and antiretroviral treatment [4–10].

Among these works, the development of virus models with immune responses is gaining much attention [3, 11, 12]. The immune system is essential in controlling the level of virus reproduction in terms of the strength of the Cytotoxic T Lymphocyte (CTL) response. A small change of the CTL response may have a large effect on virus production and infected cells load. As to this aspect, typical work can be summarized as follows. Chen et al. [12] indicated that the immunity system can not only clear free viruses but also kill infected cells. Elaiw and AlShamrani [3] proposed a four-dimensional model with humoral immunity response and general function and analyzed the global asymptotic stability of all equilibria based on the general function. However, both models in [3, 12] do not consider time delay. In order to characterize the time of a body’s immune response after the virus infection of target cells, time delay has been taken into account [13–16]. For example, Zhu et al. [16] proposed an HIV infection model with CTL response delay and analyzed the effect of time delay on the stability of equilibria. Besides, a latent period would be necessary to be incorporated into a virus model because when viruses infect a healthy organ like liver, it will not be pathogenetic at once, as it takes about six weeks to six months from the infection to the incidence [17–20]. For example, Medley et al. [17] proposed an HBV model with an exposed state, namely, infected but not yet infectious. Moreover, it was modeled in [21, 22] that the liver can regenerate cells and compensate the lost infected hepatocytes by the proliferation of hepatocytes.

In this paper, we will propose a more comprehensive model than those existing ones, which not only considers the immune response to both infected cells and viruses and a time delay for the immune system to clear viruses but also incorporates an exposed state and the proliferation of hepatocytes. We first discuss the existence of two boundary equilibria and one positive (infection) equilibrium. We then analyze the global stability of the two boundary equilibria, the local asymptotic stability and Hopf bifurcation of the positive equilibrium and also the stability of the bifurcating periodic solutions. Moreover, we perform numerical simulations to illustrate some of the theoretical results we obtain and also illustrate how the factors such as the immune response to infected cells and viruses and the proliferation of hepatocytes affect the dynamics of the model under time delay.

The paper is structured as follows. In Section 2, a delayed mathematical model is proposed, and the positivity and boundedness of solutions, existence of two boundary equilibria, and one positive equilibrium are discussed, followed by the global stability analysis of these two boundary equilibria and the local asymptotic stability and Hopf bifurcation of the positive equilibrium in Section 3. The stability of the bifurcating periodic solutions is studied in Section 4. In Section 5, some numerical simulations and discussions are given. Finally, a conclusion is given in Section 6.

2. Mathematical Model

Wang et al. [23] proposed a virus infection model of four-dimensional equations with delayed humoral immune response, which, however, does not involve an exposed state and consider the proliferation of hepatocytes. Although it considers the immune response to viruses, it does not involve the immune response to infected cells.

Based on this model, we propose a new and comprehensive HBV model, which not only considers the immune response to both infected cells and viruses and a time delay for the immune system to clear viruses but also incorporates an exposed state and the proliferation of hepatocytes. To better understand our model, we illustrate its mechanism in Figure 1.

Figure 1: The mechanism of our model.

The model is then given as follows:

where ,  ,  ,  ,  , and denote the number of uninfected cells, exposed cells, infected cells, free viruses, CTLs, and alanine aminotransferases (ALT), respectively. The parameter represents the natural production rate of uninfected cells. is a new term which is introduced to represent the proliferation of hepatocytes, where is the proliferation rate. Parameters , (and the following) ,  ,  , and represent the natural death rate of uninfected cells, exposed cells, infected cells, free viruses, CTLs, and ALT, respectively. represents the infection rate from uninfected cells to exposed cells and the transfer rate from exposed cells to infected cells. The production rate of free viruses from infected cells is denoted by , and the production rate of CTLs by . represents the production rate of ALT from the extrahepatic tissue and the production rate of ALT when the infected hepatocytes are killed by CTL. The immunity-induced clearance for infected cells is modeled by a term , where represents the clearance rate of infected cells. Similarly, the immunity-induced clearance for free viruses is modeled by , where represents the clearance rate of free viruses. is time delay. All the parameters in this paper are positive and . For convenience, we define new parameter .

2.1. Positivity and Boundedness of Solutions

In this subsection, we prove the positivity and the boundedness of solutions of system (1).

We denote ,  ,  ,  ,  ,  ,  . From the first equation of system (1), we have ; therefore, for if . Next, we consider the second, third, and fourth equation in system (1) as a nonautonomous system for :

Based on Theorem in [24], we have ,  ,   if ,  ,  .

, we have ,   if ,  .

, because , so we have ,   if .

Hence, the nonnegative is proved. In what follows, we will study the boundedness of solutions. We define as a linear combination of ,  ,  ,  ,  :

Therefore, we obtain , namely, . So we have Because the boundedness of ,  .

The boundedness is proved.

2.2. Equilibrium

In this subsection, we study the equilibria of system (1). The method to obtain equilibria is setting and computes the following:

The system (1) has two boundary equilibria (an infection-free equilibrium in which ,  ,   and an equilibrium without immune response in which ,  ,  ,  ,  ,  ) and a positive (infection) equilibrium in which ,  ,  ,  ,  ,  .

The infection-free equilibrium is , where , and , and the basic reproductive number is obtained by the following method.

Based on integral operator spectral radius, the basic reproductive number is , where

Hence, we have the basic reproductive number being .

The equilibrium without immune response is , where ,  ,  ,  , and . Similarly, we have the basic reproductive number is at .

The infected positive equilibrium is , where

3. Analysis

3.1. Global Stability Analysis of the Two Boundary Equilibria

In this section, we will employ the direct Lyapunov method and LaSalle’s invariance principle to establish the global asymptotic stability of the two boundary equilibria.

Theorem 1. The infection-free equilibrium is globally asymptotically stable if and only if .
See Appendix A for proof.

Theorem 2. The equilibrium without immune response is globally asymptotically stable if and only if .
See Appendix B for proof.

3.2. Local Asymptotic Stability and Hopf Bifurcation of the Positive Equilibrium

In this section, we will discuss the local asymptotic stability and Hopf bifurcation of the positive equilibrium .

The characteristic equation of system (1) at is as follows:

Define

Then the characteristic equation above becomes

When , (9) further becomes

where

Using the Routh-Hurwitz criterion [23], we obtain the following lemma.

Lemma 3. If (10) satisfies ,  ,   and , then the positive equilibrium is locally asymptotically stable when .

Proof. By the Routh-Hurwitz criterion, if the four conditions are satisfied, then all roots of (10) have negative real parts. Therefore, the positive equilibrium is locally asymptotically stable when . For more details, we refer the readers to [25, 26].
From Lemma 3, we know that all roots of lie to the left of the imaginary axis when . However, with increasing from zero, some of its roots may cross the imaginary axis to the right. In this case, there are some roots having positive real parts, and therefore the equilibrium becomes unstable. Next, we will discuss the stability of system (1) at when .
We first divide (9) into two parts and obtain

Suppose (9) has a purely imaginary root . Substituting into (12) yields

By separating the real part and imaginary part, the following real part is obtained:

where

Let

Therefore, if (9) has a purely imaginary root , it is equivalent to the fact that has a positive real root .

Theorem 4. If has no positive real roots, then the positive equilibrium is locally asymptotically stable for any .

Proof. If has no positive real roots, then obviously (9) has no positive real roots. Therefore, the positive equilibrium is locally asymptotically stable for any .
Substituting into (22), we obtain the real part,

and imaginary part,

Assuming that has positive real roots, denoted by . As , we then have

Let

where and .
Therefore, the characteristic equation has a pair of purely imaginary roots . For every integer and , define as the root of (9) near , satisfying and . Then the following theorem is obtained.

Theorem 5. If has some positive real roots, then is locally asymptotically stable for , where

Proof. For ,   has no positive real roots when , which means that all the roots of (9) have strictly negative real parts when . Therefore, is locally asymptotically stable for .

Theorem 6. If is a simple root of , then there is a Hopf bifurcation for the system as increases past .
See Appendix C for proof.

4. Stability of the Bifurcating Periodic Solutions

In this section, we will continue to derive the explicit formulas for determining the stability, direction, and other properties of the Hopf bifurcation at a critical value by means of the normal form and the center manifold theory [27].

First, we make the following hypotheses.

(1) Equation (9) has a pair of purely imaginary roots at , where .

(2) The remaining roots of (22) have strictly negative real parts.

(3) is a simple root of .

We use to represent a new bifurcation parameter. Let , and , where . Therefore, system (1) can be written as the following functional differential equation:

where

By the Riesz representation theorem [28], there exists a matrix-valued function such that

where .

For , we further define

Then system (22) can be written as

For , define

and an inner product of

where and . Then and are adjoint operators.

Let and be the eigenvectors of and corresponding to the eigenvalues and , respectively. We choose and as

so that is satisfied. We give the detailed computation of (29) in Appendix D.

In the following, we will compute the coefficients, ,  ,  , and , using the method given in [27]. The detailed computation of ,  ,  , and is presented in Appendix E.

Then the following values can be computed:

The signs of ,   determine the direction of the Hopf bifurcation and the stability of bifurcating periodic solutions, respectively [27]. From (C.11) of Appendix C, we obtain .

Let . We obtain the following theorem.

Theorem 7. Assume the hypotheses (1), (2), and (3) at the beginning of Section 4 hold.
(1) If , then the bifurcating periodic solutions exist for in a -neighborhood.
(2) If , the bifurcating periodic solutions are orbitally asymptotically stable as .

Proof. (1) If , where , and the hypotheses (1), (2), and (3) hold, then, from Theorem 6, we draw the conclusion that the existence and stability of the bifurcating periodic solutions are only determined by .
(2) If , namely, , then there exist stable periodic solutions for in a -neighborhood. So the bifurcating periodic solutions are orbitally asymptotically stable as .

5. Simulation and Discussions

In this section, we will numerically illustrate the theoretical results obtained above and also discuss how the factors such as the immune response to infected cells and viruses and the proliferation of hepatocytes affect the dynamics of the model under time delay.

For the following simulations, we choose the parameter values for system (1) as follows:

We set the initial values to ,  ,  ,  ,  , and , where .

5.1. Hopf Bifurcation and the Stability of Periodic Solutions

With the parameter values given in (31), we have the positive equilibrium and the critical time value .

When , we obtain stable bifurcating periodic solutions. For example, when , the simulation result is shown Figures 2 and 3. Figure 2 indicates that a stable limit cycle is obtained as expected and Figure 3 indicates the state dynamics of uninfected cells, exposed cells, infected cells, free viruses, CTLs, and ALT, which are periodically oscillating.

Figure 2: The positive equilibrium bifurcates into a periodic solution at (a limit cycle).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: The positive equilibrium bifurcates into a periodic solution at .

When , the bifurcating periodic solutions are unstable. For example, when , the simulation result is shown in Figures 4 and 5. From Figures 4 and 5, we know that the positive equilibrium is asymptotically stable and the system will converge to .

Figure 4: The positive equilibrium at .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: The positive equilibrium remains stale at .

With the increasing of time delay , the radius of limit cycle will increase. The simulation result is shown Figure 6.

Figure 6: The different periodic solutions at and (three limit cycles).

5.2. The Immune Response to Infected Cells

Here we will investigate the effect of the immune response to infected cells on the model dynamics under time delay. When we change to by fixing other values given in (31), we obtain a simulation result at , illustrated in Figure 7. Comparing Figure 2 when and Figure 7 when , we can see that, with the decrease of , the stable periodic solution becomes unstable, that is, asymptotically stable.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 7: The positive equilibrium remains stable at when (blue line) and exists at stable periodic solutions when (red line).

5.3. The Immune Response to Viruses

We continue to investigate the effect of the immune response to viruses on the model dynamics under time delay. When we change to by fixing other values given in (31), we obtain a simulation result at , illustrated in Figure 8. Similarly, comparing Figures 2 and 8, we can see that, with the decrease of , the stable periodic solution also becomes unstable, that is, asymptotically stable. We further can see that the effect of the immune response to infected cells on the model dynamics is similar to that of the immune response to viruses.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 8: The positive equilibrium remains stable at when (blue line) and exists at stable periodic solutions when (red line).

5.4. Proliferation of Hepatocytes

Then we investigate the effect of proliferation of hepatocytes on the model dynamics under time delay. For this, we still keep and change the value of parameter . When we change to by fixing other values given in (31), we obtain a simulation result at , illustrated in Figure 9. Figure 9 shows that when , the bifurcating periodic solution is stable, compared with Figure 2 when . We also try other values of parameter and obtaining similar results. Thus, we can see parameter has a small effect on the model dynamics, which reflect in periodicity and the positive equilibrium .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 9: The positive equilibrium remains stable at when (blue line) and exists at stable periodic solutions when (red line).

6. Conclusions

In this paper, we consider a comprehensive delayed HBV model. Different from other existing models, our model not only considers the immune response to both infected cells and viruses and a time delay for the immune system to clear viruses but also incorporates an exposed state and the proliferation of hepatocytes.

We then prove the positivity and boundedness of solutions and analyze the global stability of two boundary equilibria and investigate the local asymptotic stability and Hopf bifurcation of the positive (infection) equilibrium and also the stability of the bifurcating periodic solutions. We also numerically illustrate the Hopf bifurcation and the stability of the bifurcating periodic solutions.

Moreover, we numerically illustrate how the factors such as the time delay, the immune response to infected cells and viruses, and the proliferation of hepatocytes affect the dynamics of the model, which shows that the former two factors have a big effect on the model dynamics, while the latter one does not have a big effect.

Appendix

A. The Proof of Theorem 1

We construct a Lyapunov functional as follows:

Calculating the derivative of along with the trajectories of system (1), we obtain

Therefore, when ,   holds true. Furthermore, if and only if , and , the Lyapunov functional satisfies . According to LaSalle’s invariance principle [29], the infection-free equilibrium is globally asymptotically stable when . Apparently, when ,   is unstable.

B. The Proof of Theorem 2

Similarly, we define the following Lyapunov functional:

The derivative of along with the trajectories of system (1) can be calculated as

Therefore, ensures that holds true. Furthermore, if and only if ,  ,  ,  ,  , and , the Lyapunov functional satisfies . Using the LaSalle’s invariance principle, we can see that the equilibrium without immune response is globally asymptotically stable.

C. The Proof of Theorem 6

The characteristic equation (9) can be written into the following form:

where

and and are continuously differentiable to .

Suppose that one of the roots of (C.1) is , satisfying and for a positive real number .

Let

Calculating the derivative of to , we have

Then we have

Because , we have

Calculating the derivative of both sides of (C.1) to , we have

Then we have

Since , we obtain

Thus, we have

When , obviously, we have

As is a simple root of , we know . From (C.11), we further know . If ; we obtain that the roots of (9) have positive real part when , which contrasted with Theorem 5. Hence, we can see that . When , except for the pair of purely imaginary roots, the remaining roots of have strictly negative real parts, so the system has Hopf bifurcation.

D. The Detailed Computation of (29)

In what follows, we will compute (29).

Substituting into (26), we have

Therefore, is the characteristic equation of .

As , we have

where

So we have

The eigenvectors of are

Similarly,

The eigenvectors of are

Since

we have that

E. The Detailed Computation of , , , and

In the following, we will compute coefficients, ,  ,  , and , using the method given in [27].

Let

where

Then we have

We further obtain

Using the method [27], we obtain the following coefficients:

where

where

So far, we have obtained ,  ,  , and .

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by National Natural Science Foundation of China (61273226).

References

F. M. Lu and H. Zhuang, “Management of hepatitis B in China,” Chinese Medical Journal, vol. 122, pp. 3-4, 2009.

World Health Organization, “Hepatitis B,” 2008, http://www.who.int/mediacentre/factsheets/fs204/en/.

A. M. Elaiw and N. H. AlShamrani, “Global stability of humoral immunity virus dynamics models with nonlinear infection rate and removal,” Nonlinear Analysis: Real World Applications, vol. 26, pp. 161–190, 2015.

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.

K. Manna and S. P. Chakrabarty, “Chronic hepatitis B infection and HBV DNA-containing capsids: modeling and analysis,” Communications in Nonlinear Science and Numerical Simulation, vol. 22, no. 1-3, pp. 383–395, 2015.

Y. Su, D. Sun, and L. Zhao, “Global analysis of a humoral and cellular immunity virus dynamics model with the Beddington-DeAngelis incidence rate,” Mathematical Methods in the Applied Sciences, vol. 38, no. 14, pp. 2984–2993, 2015.

S. Wang, X. Feng, and Y. He, “Global asymptotical properties for a diffused HBV infection model with CTL immune response and nonlinear incidence,” Acta Mathematica Scientia B, vol. 31, no. 5, pp. 1959–1967, 2011.

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.

C. Monica and M. Pitchaimani, “Analysis of stability and Hopf bifurcation for HIV-1 dynamics with PI and three intracellular delays,” Nonlinear Analysis: Real World Applications, vol. 27, pp. 55–69, 2016.

J. Li, K. Wang, and Y. Yang, “Dynamical behaviors of an HBV infection model with logistic hepatocyte growth,” Mathematical and Computer Modelling, vol. 54, no. 1-2, pp. 704–711, 2011.

M. Prakash and P. Balasubramaniam, “Bifurcation analysis of macrophages infection model with delayed immune response,” Communications in Nonlinear Science and Numerical Simulation, vol. 35, pp. 1–16, 2016.

X. Chen, L. Q. Min, and Q. L. Sun, “Dynamics analysis and numerical simulation of an amended HBV infection model,” Journal of Biomathematics, vol. 28, no. 2, pp. 278–284, 2013.

C. Sun, Z. Cao, and Y. Lin, “Analysis of stability and Hopf bifurcation for a viral infectious model with delay,” Chaos, Solitons & Fractals, vol. 33, no. 1, pp. 234–245, 2007.

M. Y. Li and H. Shu, “Global dynamics of a mathematical model for HTLV-I infection of CD4⁺ T cells with delayed CTL response,” Nonlinear Analysis: Real World Applications, vol. 13, no. 3, pp. 1080–1092, 2012.

X. Tian and R. Xu, “Global stability and Hopf bifurcation of an HIV-1 infection model with saturation incidence and delayed CTL immune response,” Applied Mathematics and Computation, vol. 237, pp. 146–154, 2014.

H. Zhu, Y. Luo, and M. Chen, “Stability and Hopf bifurcation of a HIV infection model with CTL-response delay,” Computers & Mathematics with Applications, vol. 62, no. 8, pp. 3091–3102, 2011.

T. Zhang and Z. Teng, “On a nonautonomous SEIRS model in epidemiology,” Bulletin of Mathematical Biology, vol. 69, no. 8, pp. 2537–2559, 2007.

K. W. Blayneh, A. B. Gumel, S. Lenhart, and T. Clayton, “Backward bifurcation and optimal control in transmission dynamics of West Nile virus,” Bulletin of Mathematical Biology, vol. 72, no. 4, pp. 1006–1028, 2010.

B.-Z. Guo and L.-M. Cai, “A note for the global stability of a delay differential equation of hepatitis B virus infection,” Mathematical Biosciences and Engineering, vol. 8, no. 3, pp. 689–694, 2011.

S. Zhang and Y. Zhou, “Dynamic analysis of a hepatitis B model with three-age-classes,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 7, pp. 2466–2478, 2014.

H. Dahari, A. Lo, R. M. Ribeiro, and A. S. Perelson, “Modeling hepatitis C virus dynamics: liver regeneration and critical drug efficacy,” Journal of Theoretical Biology, vol. 247, no. 2, pp. 371–381, 2007.

G. Pachpute and S. P. Chakrabarty, “Dynamics of hepatitis C under optimal therapy and sampling based analysis,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 8, pp. 2202–2212, 2013.

T. Wang, Z. Hu, and F. Liao, “Stability and Hopf bifurcation for a virus infection model with delayed humoral immunity response,” Journal of Mathematical Analysis and Applications, vol. 411, no. 1, pp. 63–74, 2014.

H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperatiive Systems, vol. 41 of Mathematical Surveys and Monographs, American Mathematical Society, Providence, RI, USA, 1995.

M. Iannelli, Mathematical Theory of Age-Structured Population Dynamics, Giardini Editori E Stampatori, Pisa, Italy, 1995.

G. F. Webb, Theory of Nonlinear Age-Dependent Population Dynamics, vol. 89, Marcel Dekker, New York, NY, USA, 1985.

B. D. Hassard, N. D. Kazarinoff, and Y.-H. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, 1981.

R. K. Goodrich, “A Riesz representation theorem,” Proceedings of the American Mathematical Society, vol. 24, pp. 629–636, 1970.

A. M. Lyapunov, “The general problem of the stability of motion,” International Journal of Control, vol. 55, no. 3, pp. 521–790, 1992.