Research Article  Open Access
Yongmei Su, Sinuo Liu, Shurui Song, Xiaoke Li, Yongan Ye, "Stability Analysis and Clinic Phenomenon Simulation of a FractionalOrder HBV Infection Model", Complexity, vol. 2020, Article ID 8864403, 11 pages, 2020. https://doi.org/10.1155/2020/8864403
Stability Analysis and Clinic Phenomenon Simulation of a FractionalOrder HBV Infection Model
Abstract
In this paper, a fractionalorder HBV model was set up based on standard mass action incidences and quasisteady assumption. The basic reproductive number and the cytotoxic T lymphocytes’ immuneresponse reproductive number were derived. There were three equilibrium points of the model, and stable analysis of each equilibrium point was given with corresponding hypothesis about or . Some numerical simulations were also given based on HBeAg clinical data, and the simulation showed that there existed positive logarithmic correlation between the number of infected cells and HBeAg, which was consistent with the clinical facts. The simulation also showed that the clinical individual differences should be reflected by the fractionalorder model.
1. Introduction
Viral infection is a major global problem, and mathematical models are an important tool for the study of biological phenomenon [1–5] and viral infectious disease [6–9] because they can help us to understand the dynamics of some infectious diseases and some chronic viral infections.
Mathematical models were also used to interpret experimental and clinical results in the fields of (anti) HIV, HBV, and HCV infections [10–12]. These models were all set up with ordinary differential equation. In recent years, fractional differential equation models were often used in biology because the researchers found that the biological cell membranes have electron conductivity, which could be classified as a fractionalorder model [13, 14]. In addition, some biological models established by fractional differential equations have proved to be more advantageous than integers [14]. In particular, the biggest difference between the fractionalorder model and the integerorder model is that the fractionalorder model has the memory, while the characteristic of the immune response contains the memory [14].
So, when we discuss virus immune models, fractional mathematical models have become important tools. Arafa et al. [14] proposed a fractionalorder HIV infection model, Wang et al. [15] proposed a fractionalorder HIV infection model, considering the logistic growth of the healthy CD4 cells, and Yan and Kou [16] had further proposed the following HIV model:
Withinhost HBV models of fractional order were also discussed [17, 18], and the model in paper [17] was as follows:where represent the concentration of uninfected hepatocytes cells, infected cells, and viruses, the death rate of them is , the uninfected hepatocytes are supposed to be produced at rate , is the infection rate at which the uninfected cell becomes infected, and the infected hepatocytes are cured by noncytolytic processes at a constant rate per cell. Based on model (2), Cardoso et al. [18] also disused a fractional model of hepatitis B with drug therapy by representing in model (2) with and , in which and represent the drug efficacy.
It should be pointed out that the bilinear incidences were used to describe the infection between uninfected cells and virus in [14–18], but Min et al. changed it to standard incidence function to describe the HBV infection model, which seemed more reasonable because it is independent of the number of total cells of liver [19]. On the other hand, Bartholdy et al. and Wodarz et al. [20] found that the turnover of free virus was much faster than that of infected cells, which means quasisteady assumption could be used, that is, the amount of free virus is simply proportional to the number of infected cells, so the number of infected cells can also be considered as a measure of virus load. Based on quasisteady assumption, Guo and Cai [21] discussed a HBV infection model by instead of , but the model is integral order. On the other hand, model (2) does not include the immune cell, while the cytotoxic T lymphocyte (CTL) immune response after viral infection is universal and necessary to eliminate or control the disease.
Though fractional differential equations have proved to be a good choice to describe biological phenomenon, most discussions only focus on mathematic analysis and numerical simulation, and so far, almost no papers have used fractionalorder models to explain clinic phenomenon about HBV. In this paper, based on the above discussion and quasisteady assumption, a HBV model was set up as follows:in which , the meanings of x and y are the same as those in model (2), z represents the number of CTL, where the immune response is assumed to get stronger at a rate and decays exponentially at a rate , which is proportional to their current concentration, and the parameter expresses the efficacy of the nonlytic component.
This paper is organized as follows. In Section 2, some definitions and lemmas of fractionalorder differential equation are cited. In Section 3, we mainly discuss the existence and uniqueness of positive solutions. The stability analysis is given in Section 4. The simulation is given in Section 5. This paper ends with a conclusion in Section 6 and discussion in Section 7.
2. Primary Concept and Lemma
It is known that fractional derivative has a variety of definitions [17, 22]. In this paper, we used the Caputo fractional derivatives which are defined as follows.
Definition 1. (see [14]). The Caputo fractional derivatives of order are defined aswhere
Definition 2. (see [20]). The discriminant of a polynomialis defined by , where is the derivative of . is the determinant of the corresponding Sylvester matrix. The Sylvester matrix is formed by filling the matrix beginning with the upper left corner with the coefficients of and then shifting down one row and one column to the right and filling in the coefficients starting there until they hit the right side. The process is then repeated for the coefficients of .
The following lemmas were useful to our arguments.
Lemma 1 (see [23]). For a fractionalorder system:with , the equilibrium point of the system is locally asymptotically stable if all eigenvalues of Jacobian matrix evaluated at the equilibrium point satisfy
Lemma 2 (see [24]). For the polynomial equation,
The conditions which make all the roots of (8) satisfy (7) are displayed as follows:(i)For n = 1, the condition is (ii)For n = 2, the conditions are either Routh–Hurwitz conditions or(iii)For n = 3,(a)If the discriminant of is positive, then Routh–Hurwitz conditions are the necessary and sufficient conditions, that is, (b)If then (7) holds when (c)If and then (7) holds when (d)If and , then (7) holds for all
Lemma 3 (see [25]). Assume that the vector function f: satisfies the following conditions:(1)Function is Lebesgue measurable with respect to on (2)Function is continuous with respect to on (3) is continuous with respect to on (4)Here and are two positive constants; then, initial value problems , , have a unique solution.
Lemma 4 (see [26]). Let and assume that , ; then, .
3. The Existence and Uniqueness of Positive Solutions
For the proof of the existence and uniqueness about the positive solution, we firstly prove that there exists a positively invariant region for system (3).
Letand we havein which , so we have
Let ; it is easy to see that is a positively invariant region for model (3).
Theorem 1. For any initial condition in D, system (3) has a unique solution , and the solution will remain nonnegative for all .
Proof 1. Firstly, we prove the existence and uniqueness of the solution. We denoteObviously satisfies conditions (1)–(3) of Lemma 3, and we only prove that system (3) satisfies the last condition (4) of Lemma 3. LetHence,in which . By Lemma 3, system (3) has a unique solution.
We now prove the solution is nonnegative for all . From model (3), we knowBy Lemma 4, the solution is nonnegative. In summary, the system has a unique nonnegative solution.
4. Stable Analysis
In this section, we will discuss the stability of model (3). This system always has an infectionfree equilibrium , where .
The basic reproduction number is
When system (3) will have an immuneabsence equilibrium wherewhich means the uninfected cells and infected cells coexist but the immune response is not activated yet, that is, . Further, we will give the immuneresponse reproductive number when , that is, , which means immune response is activated. So, when there is another immuneresponse equilibrium , wherein which . Note that , , and must be positive, and we can prove holds only when . Hence, exists if and only if .
Now, we introduce the main theorem.
Theorem 2. For system (3), if the equilibrium is globally asymptotically stable. If , the equilibrium is unstable.
Proof 2. Jacobian matrix J at isThe characteristic equation for the infectionfree equilibrium is given as follows:We can see that the characteristic roots , and when , , all the characteristic roots satisfied , which shows that is locally asymptotically stable by Lemma 1. If is unstable.
Letand we haveSince , we have . Let ; obviously,. Let be the largest positively invariant subset of ; obviously, E is nonempty since . Let (x(t), y(t), z(t)) be the solution of system (3) with initial value ; then, By the first and the third equations of (3), we haveand when, , , by the invariance of E, , . By the Lyapunov–LaSalle theorem [26], when , all solutions in the set approach . Noting that is locally asymptotically stable, the infectionfree equilibrium is globally asymptotically stable.
Theorem 3. For system (3), when ,(1)The equilibrium is locally asymptotically stable if (2)If , the equilibrium is unstable
Proof 3. The characteristic equation for the equilibrium is given asthat is,We can see , and when , is negative and .
Next, we consider the following equation:We letSince , we can see ; hence, (27) has two negative real roots, and we denote them by and ( and ), so when , the equilibrium is locally asymptotically stable. When , is positive, and , is unstable.
Finally, we discuss the local stability of the immuneresponse equilibrium state , and the Jacobian matrix at is given byThe corresponding characteristic equation iswhereIt is easy to verify thatBased on Definition 2, we obtainBy Lemma 2, we have the following theorem.
Theorem 4. For system (3), when ,(1)If the discriminant of is positive, the equilibrium is locally asymptotically stable for (2)If then the equilibrium is locally asymptotically stable for
5. Numerical Simulation
5.1. Simulation of the ImmuneResponse Equilibrium
We first simulate the stability of , and the locally stable condition of is given by Theorem 4. The parameters we used are listed in Table 1 which can make , hold. We choose , and , which satisfy and the simulation is shown in Figure 1. The dynamic routes of uninfected and infected cells are shown in Figures 1(a) and 1(b), and the CTL cells are shown in Figure 1(c). The simulation shows that the immuneresponse equilibrium is locally asymptotically stable for which is consistent with the conclusion of Theorem 4.

(a)
(b)
(c)
5.2. Simulation of Correlation between HBeAg and Infected Cells
A positive correlation was found between logHBV DNA and logHBeAg in [27], that is, the higher the HBeAg, the higher the HBV DNA load. Since by quasisteady assumption, the amount of free virus is simply proportional to the number of infected cells, the number of infected cell should also be positive logarithmic correlative with HBeAg. Since we could not get the specific number of infected cells from clinic, HBeAg can be quantified. In the following part, we will testify the positive correlation between the infected cell and the HBeAg based on the clinical data of HBeAg by numerical simulation. We use four CHB patients’ HBeAg clinical data of 108 weeks with the treatment of entecavir from Dongzhimen Hospital; two patients were from successfultreatment group, and the other two patients were from unsuccessfultreatment group; the data are shown in Tables 2 and 3.


By Theorem 2, the infectionfree equilibrium is asymptotically stable if . When simulating the successfultreatment group, the parameters we choose should make . Considering the biological significance of the parameters, we choose the parameters in Table 4 as follows.

A human liver contains approximately hepatocytes [28], and a patient has a total of about 3000 ml plasma. Usually, clinical testing quantification is based on one millilitre. Consequently, we can assume that . Since the halflife of a hepatocyte is about half a year, we choose similar to [28]. The halflife of CTL is about 77 days [29], so we choose . Other parameters were chosen by simulation. The parameters in Table 4 can make hold. We give three simulation results with , , and ; the simulations are shown in Figure 2. The dynamic routs of uninfected cells and CTL cells are shown in Figures 2(a) and 2(b), and the HBeAg and infected cells are shown in Figure 2(c). The simulation shows that the infectionfree equilibrium is asymptotically stable with different which is consistent with the conclusion of Theorem 2, but when , the dynamic route is obviously different. From the simulation in Figure 2(c), we can also observe the positive correlation between loginfected cells and logHBeAg for different , which is consistent with clinical fact. On the other hand, from the clinical data of HBeAg, we can see that the change rate was different between different patient, though they have same change trend, which shows that individual differences do exist in clinical therapy even with the same medicine and the same dose; the clinical phenomenon may be explained by our model to some extent. From the simulation, we can see that the larger the , the slower the decrease rate of infected cells and CTL in the former stage, but the faster the decrease rate in the later stage; the individual difference may be reflected by different , and the fractionalorder model should be a good tool to describe HBV infection.
(a)
(b)
(c)
When we simulate the unsuccessfultreatment group in Table 3, the patients’ ALT level was lower than normal standard, so we can think that the patient has not yet developed an immune response to the hepatitis B virus. When , there exists the equilibrium without immune activities, so we choose the immune parameter b larger than that in Table 4 and the immune parameter c smaller than that in Table 4. The death rate of infected cells should decrease without immune activities, so we choose parameter a which is smaller than that in Figure 2. Since is asymptotically stable, if , the parameter we choose in Table 5 can make and hold. We also used , , and to simulate, and the simulation is shown in Figure 3. The dynamic routes of uninfected cells and CTL cells are shown in Figures 3(a) and 3(b), and the HBeAg and infected cells are shown in Figure 3(c). From the simulation, the dynamic routes of uninfected cells, infected cells, and CTL cells show that the immuneabsence equilibrium is asymptotically stable with different which is consistent with the conclusion of Theorem 3. From the simulation, there still exists the same change trend but different change rate with different , and the positive correlation between loginfected cells and logHBeAgis reflected roughly. The positive correlation is not so good for patient 2’s HBeAg data, but the trend is consistent between the infected cells and HBeAg. We can see there also exists some difference with different which may reflect the individual difference.

(a)
(b)
(c)
6. Conclusions
In this paper, based on the fact that immune response has memory, we discussed a fractionalorder HBV model with standard mass action incidences, and we obtained the basic reproductive numbers and the cytotoxic T lymphocytes’ immuneresponse reproductive number . When , we have proved that is globally asymptotically stable which meant the infected person can eventually recover automatically even when infected with a large number of HBV. When and , was locally asymptotically stable which meant the one infected by HBV with no immune response would be infected persistently. When , we also gave the local stable condition of , which meant the infected person has immune response to HBV, but the persistent infection still existed. Furthermore, we gave some simulations with different to test our theoretical results and some clinical phenomena. The simulation showed that our model can simulate the positive correlation between loginfected cells and logHBeAg to some extent. On the other hand, even with the same initial conditions and the same parameters, we also found that there existed some difference in the dynamic routes with different order , which may reflect the individual difference. From the analysis and simulation, we can see the fractionalorder model maybe more reasonable to describe HBV immune course.
7. Discussion
In our paper, we found that the fractionalorder model can reflect the characteristics of immune memory and clinical individual differences, but it should be pointed out that there are only three variables in our model; the model only described HBV infection roughly. Furthermore, we could not obtain the number of uninfected cells, infected cells, and CTL; when we performed the simulation, we only chose the parameters from biological meaning, not from real clinical data, so there are still many works to be done in the future.
Data Availability
The HBeAg data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This study was jointly supported by 2015 National Traditional Medicine Clinical Research Base Business Construction Special Topics (JDZX2015299) and Fundamental Research Funds for the Central Universities (FRFBR16019A).
References
 Y. Xiao, X. Xu, and S. Tang, “Sliding mode control of outbreaks of emerging infectious diseases,” Bulletin of Mathematical Biology, vol. 74, no. 10, pp. 2403–2422, 2012. View at: Publisher Site  Google Scholar
 Y. Zheng, L. Q. Min, Y. Ji, Y. Su, and Y. Kuang, “Global stability of endemic equilibrium point of basic virus infection model with application to HBV infection,” Journal of Systems Science and Complexity, vol. 23, no. 6, pp. 1221–1230, 2009. View at: Publisher Site  Google Scholar
 G. Huang, W. Ma, and Y. Takeuchi, “Global properties for virus dynamics model with BeddingtonDeAngelis functional response,” Applied Mathematics Letters, vol. 22, no. 11, pp. 1690–1693, 2009. View at: Publisher Site  Google Scholar
 G. Huang, H. Yokoi, Y. Takeuchi, T. Kajiwara, and T. Sasaki, “Impact of intracellular delay, immune activation delay and nonlinear incidence on viral dynamics,” Japan Journal of Industrial and Applied Mathematics, vol. 28, no. 3, pp. 383–411, 2011. View at: Publisher Site  Google Scholar
 T. Q. Zhang, J. L. Wang, Y. Li, Z. Jiang, and X. F. Han, “Dynamics analysis of a delayed virus model with two different transmission methods and treatments,” Advances in Difference Equations, vol. 2020, no. 1, pp. 1–17, 2020. View at: Publisher Site  Google Scholar
 W. Wang, W. B. Ma, and Z. S. Feng, “Dynamics of reaction–diffusion equations for modelling CD4^{+} T cells decline with general infection mechanism and distinct dispersal rate,” Nonlinear Analysis: Real World Applications, vol. 51, no. 2, Article ID 102976, 2020. View at: Publisher Site  Google Scholar
 C. Yan and W. Wang, “Modeling HIV dynamics under combination therapy with inducers and antibodies,” Bulletin of Mathematical Biology, vol. 81, no. 7, pp. 2625–2648, 2019. View at: Publisher Site  Google Scholar
 R. Q. Shi, T. Lu, and C. H. Wang, “Dynamic analysis of a fractional order delayed model for hepatitis B virus with CTL immune response,” Virus Research, vol. 277, no. 2, Article ID 197841, 2020. View at: Publisher Site  Google Scholar
 S. He, S. Y. Tang, S. Tang, and L. Rong, “A discrete stochastic model of the COVID19 outbreak: forecast and control,” Mathematical Biosciences and Engineering, vol. 17, no. 4, pp. 2792–2804, 2020. View at: Publisher Site  Google Scholar
 A. S. Perelson, “Modelling viral and immune system dynamics,” Nature Reviews Immunology, vol. 2, no. 1, pp. 28–36, 2002. View at: Publisher Site  Google Scholar
 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, vol. 93, no. 9, pp. 4398–4402, 1996. View at: Publisher Site  Google Scholar
 G. K. Lau, M. Tsiang, J. Hou et al., “Combination therapy with lamivudine and famciclovir for chronic hepatitis Binfected Chinese patients: a viral dynamics study,” Hepatology, vol. 32, no. 2, pp. 394–399, 2000. View at: Publisher Site  Google Scholar
 L. Zhang, G. Huang, A. Liu, and R. Fan, “Stability analysis for a fractional HIV infection model with nonlinear incidence,” Discrete Dynamics in Nature and Society, vol. 2015, Article ID 563127, 11 pages, 2015. View at: Publisher Site  Google Scholar
 A. A. M. Arafa, S. Z. Rida, and M. Khalil, “A fractionalorder model of HIV infection with drug therapy effect,” Journal of the Egyptian Mathematical Society, vol. 22, no. 3, pp. 538–543, 2014. View at: Publisher Site  Google Scholar
 Y. Wang, L. Liu, X. Zhang, and Y. Wu, “Positive solutions of an abstract fractional semipositone differential system model for bioprocesses of HIV infection,” Applied Mathematics and Computation, vol. 258, pp. 312–324, 2015. View at: Publisher Site  Google Scholar
 Y. Yan and C. Kou, “Stability analysis for a fractional differential model of HIV infection of CD4^{+} Tcells with time delay,” Mathematics and Computers in Simulation, vol. 82, no. 9, pp. 1572–1585, 2012. View at: Publisher Site  Google Scholar
 S. M. Salman and A. M. Yousef, “On a fractionalorder model for HBV infection with cure of infected cells,” Journal of the Egyptian Mathematical Society, vol. 25, no. 4, pp. 445–451, 2017. View at: Publisher Site  Google Scholar
 L. C. Cardoso, F. L. P. Dos Santos, and R. F. Camargo, “Analysis of fractionalorder models for hepatitis B,” Computational and Applied Mathematics, vol. 37, no. 4, pp. 4570–4586, 2018. View at: Publisher Site  Google Scholar
 Y. Su, L. Zhao, and L. Min, “Analysis and simulation of an adefovir antihepatitis B virus infection therapy immune model with alanine aminotransferase,” IET Systems Biology, vol. 7, no. 5, pp. 205–213, 2013. View at: Publisher Site  Google Scholar
 K. Wang, W. Wang, H. Pang, and X. Liu, “Complex dynamic behavior in a viral model with delayed immune response,” Physica D: Nonlinear Phenomena, vol. 226, no. 2, pp. 197–208, 2007. View at: Publisher Site  Google Scholar
 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 & Engineering, vol. 8, no. 3, pp. 689–694, 2012. View at: Publisher Site  Google Scholar
 D. Tavares, R. Almeida, and D. F. M. Torres, “Caputo derivatives of fractional variable order: numerical approximations,” Communications in Nonlinear Science and Numerical Simulation, vol. 35, pp. 69–87, 2016. View at: Publisher Site  Google Scholar
 Y. Emvudu, D. Bongor, and R. Koïna, “Mathematical analysis of HIV/AIDS stochastic dynamic models,” Applied Mathematical Modelling, vol. 40, no. 2122, pp. 9131–9151, 2016. View at: Publisher Site  Google Scholar
 E. Ahmed, A. M. A. ElSayed, and H. A. A. ElSaka, “On some RouthHurwitz conditions for fractional order differential equations and their applications in Lorenz, Rössler, Chua and Chen systems,” Physics Letters A, vol. 358, no. 1, pp. 1–4, 2006. View at: Publisher Site  Google Scholar
 J. Huo, H. Zhao, and L. Zhu, “The effect of vaccines on backward bifurcation in a fractional order HIV model,” Nonlinear Analysis: Real World Applications, vol. 26, pp. 289–305, 2015. View at: Publisher Site  Google Scholar
 M. Çiçek, C. Yakar, and B. Oğur, “Stability, boundedness, and Lagrange stability of fractional differential equations with initial time difference,” The Scientific World Journal, vol. 2014, Article ID 939027, 2014. View at: Publisher Site  Google Scholar
 W. J. Yin, F. Wang, L. P. Shen et al., “The positive correlation of log hepatitis B virus (HBV) deoxyribonucleic acid load and log HBeAg and log antibody to HBeAg,” Chinese Journal of Vaccines and Immunization, vol. 18, no. 4, pp. 316–319, 2012. View at: Google Scholar
 L. Q. Min, Y. M. 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. View at: Publisher Site  Google Scholar
 M. Hellerstein, M. B. Hanley, D. Cesar et al., “Directly measured kinetics of circulating T lymphocytes in normal and HIV1infected humans,” Nature Medicine, vol. 5, no. 1, pp. 83–89, 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Yongmei Su 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.