A Fractional Order Model for Viral Infection with Cure of Infected Cells and Humoral Immunity

International Journal of Differential Equations, Dec 2018

In this paper, we study the dynamics of a viral infection model formulated by five fractional differential equations (FDEs) to describe the interactions between host cells, virus, and humoral immunity presented by antibodies. The infection transmission process is modeled by Hattaf-Yousfi functional response which covers several forms of incidence rate existing in the literature. We first show that the model is mathematically and biologically well-posed. By constructing suitable Lyapunov functionals, the global stability of equilibria is established and characterized by two threshold parameters. Finally, some numerical simulations are presented to illustrate our theoretical analysis.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://downloads.hindawi.com/journals/ijde/2018/1019242.pdf

A Fractional Order Model for Viral Infection with Cure of Infected Cells and Humoral Immunity

A Fractional Order Model for Viral Infection with Cure of Infected Cells and Humoral Immunity Adnane Boukhouima,1 Khalid Hattaf,1,2 and Noura Yousfi1 1Laboratory of Analysis, Modeling and Simulation (LAMS), Faculty of Sciences Ben M’sik, Hassan II University, P.O Box 7955 Sidi Othman, Casablanca, Morocco 2Centre Régional des Métiers de l’Education et de la Formation (CRMEF), 20340 Derb Ghalef, Casablanca, Morocco Correspondence should be addressed to Khalid Hattaf; rf.oohay@fattah.k Received 9 September 2018; Accepted 7 November 2018; Published 2 December 2018 Guest Editor: Nurcan B. Savaşaneril Copyright © 2018 Adnane Boukhouima et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract In this paper, we study the dynamics of a viral infection model formulated by five fractional differential equations (FDEs) to describe the interactions between host cells, virus, and humoral immunity presented by antibodies. The infection transmission process is modeled by Hattaf-Yousfi functional response which covers several forms of incidence rate existing in the literature. We first show that the model is mathematically and biologically well-posed. By constructing suitable Lyapunov functionals, the global stability of equilibria is established and characterized by two threshold parameters. Finally, some numerical simulations are presented to illustrate our theoretical analysis. 1. Introduction The immune response plays an important role to control the dynamics of viral infections such as human immunodeficiency virus (HIV), hepatitis B virus (HBV), hepatitis C virus (HCV), and human T-cell leukemia virus (HTLV). Therefore, many mathematical models have been developed to incorporate the role of immune response in viral infections. Some of these models considered the cellular immune response mediated by cytotoxic T lymphocytes (CTL) cells that attack and kill the infected cells [1–5] and the others considered the humoral immune response based on the antibodies which are produced by the B-cells and are programmed to neutralize the viruses [6–11]. However, all these models have been formulated by using ordinary differential equations (ODEs) in which the memory effect is neglected while the immune response involves memory [12, 13]. Fractional derivative is a generalization of integer derivative and it is a suitable tool to model real phenomena with memory which exists in most biological systems [14–16]. The fractional derivative is a nonlocal operator in contrast to integer derivative. This means that if we want to compute the fractional derivative at some point , it is necessary to take into account the entire history from the starting point up to the point . For these reasons, modeling some real process by using fractional derivative has drawn attention of several authors in various fields [17–22]. In biology, it has been shown that the fractional derivative is useful to analyse the rheological proprieties of cells [23]. Furthermore, it has been deduced that the membranes of cells of biological organism have fractional order electrical conductance [24]. Recently, much works have been done on modeling the dynamics of viral infections with FDEs [25–31]. These works ignored the impact of the immune response and the majority of them deal only with the local stability. In some viral infections, the humoral immune response is more effective than cellular immune response [32]. For this reason, we improve the above ODE and FDE models by proposing a new fractional order model that describes the interactions between susceptible host cells, viral particles, and the humoral immune response mediated by the antibodies; that is,where , , , , and are the concentrations of susceptible host cells, latently infected cells (infected cells which are not yet able to produce virions), productive infected cells, free virus particles, and antibodies at time , respectively. Susceptible host cells are assumed to be produced at a constant rate , die at the rate , and become infected by virus at the rate . Latently infected cells die at the rate and return to the uninfected state by loss of all covalently closed circular DNA (cccDNA) from their nucleus at the rate . Productive infected cells are produced from latently infected cells at the rate and die at the rate . Free virus particles are produced from productive infected cells at the rate , cleared at the rate , and are neutralized by antibodies at the rate . Antibodies are activated against virus at the rate and die at the rate . In system (1), represents the Caputo fractional derivative of order defined for an arbitrary function by with [33]. Further, the infection transmission process in (1) is modeled by Hattaf-Yousfi functional response [34] which was recently used in [35, 36] and has the form , where are the saturation factors measuring the psychological or inhibitory effect and is the infection rate. In addition, this functional response generalizes many common types existing in the literature such as the specific functional response proposed by Hattaf et al. in [37] and used in [2, 31] when ; the Crowley-Martin functional response introduced in [38] and used in [39] when and ; and the Beddington-DeAngelis functional response proposed in [40, 41] and used in [3, 4, 10] when and . Also, the Hattaf-Yousfi functional response is reduced to the saturated incidence rate used in [9] when and and the standard incidence function used in [27] when and , and it was simplified to the bilinear incidence rate used in [5, 6] when and . On the other hand, system (1) becomes a model with ODEs when , which improves and generalizes the ODE model with bilinear incidence rate [42], the ODE model with saturated incidence rate [43], and the ODE model with specific functional response [44]. The rest of the paper is organized as follows. The next section deals with some basic proprieties of the solutions and the existence of equilibria. The global stability of equilibria is established in Section 3. To verify our theoretical results, we provide some numerical simulations in Section 4, and we conclude in Section 5. 2. Basic Properties and Equilibria In this section, we will show that our model is well-posed and we discuss the existence of equilibria. Since system (1) describes the evolution of cells, then we need to prove that the cell numbers should remain nonnegative and bounded. For biological considerations, we assume that the initial conditions of (1) satisfyThen we have the following result. Theorem 1. Assume that the initial conditions satisfy (3). Then there exists a unique solution of system (1) defined on . Moreover, this solution remains nonnegative and bounded for all . Proof. First, system (1) can be written as follows:whereIt is important to note that when , (4) becomes a system with ODEs. In this case, we refer the reader to [45] for the existence of solutions and to the works [46–50] for the stability of equilibria. In the case of FDEs, we will use Lemma 2.4 in [31] to prove the existence and uniqueness of solutions. Hence, we put We discuss four cases:(i)If , can be formulated as follows:  where  Hence,(ii)If , we can write in the form  where  Moreover, we get (iii)If , we have  where  Further, we obtain (iv)If , we have  where  Then Hence, the conditions of Lemma 2.4 in [31] are verified. Then system (1) has a unique solution on . Now, we show the nonnegativity of solutions. By (1), we have As in [31, Theorem 2.7], we deduce that the solution of (1) is nonnegative. Finally, we prove the boundedness of solutions. We define the function Then, we have where . Thus, we obtain Since , we get This completes the proof. Now, we discuss the existence of equilibria. It is clear that system (1) has always an infection-free equilibrium . Then the basic reproduction number of (1) is as follows: To find the other equilibria, we solve the following system:From (29), we get or . Then we discuss two cases. If , by (25)-(28), we have , , , andSince , , and , then . Consequently, there is no equilibrium when . We define the function on by We have , , and . Hence if , (30) has a unique root . As a result, when there exists an equilibrium satisfying , , , and . If , then . By (25)-(27), we obtain , , , andSince , , and , we have . Hence, there is no equilibrium if . We define the function on by We have , , and . Let us introduce the reproduction number for humoral immunity as follows: which denotes the average life expectancy of antibodies and is the number of free viruses at . For the biological significance, represents the average number of the antibodies activated by virus. If , we have and Therefore, there is no equilibrium when . If , then and In this case, (32) has one root . Consequently, when , there exists an equilibrium satisfying , , , , and . When , We summarize the above discussions in the following theorem. Theorem 2.  (i)If , then system (1) has one infection-free equilibrium of the form , where .(ii)If , then system (1) has an infection equilibrium without humoral immunity of the form , where , , , and .(iii)If , then system (1) has an infection equilibrium with humoral immunity of the form , where , , , , and . 3. Global Stability of Equilibria In this section, we focus on the global stability of equilibria. Theorem 3. If , then the infection-free equilibrium is globally asymptotically stable and it becomes unstable if . Proof. The proof of the first part of this theorem is based on the construction of a suitable Lyapunov functional that satisfies the conditions given in [51, Lemma 4.6]. Hence, we define a Lyapunov functional as follows: where for . It is not hard to show that the functional is nonnegative. In fact, the function has a global minimum at . Consequently, for all . Calculating the fractional derivative of along solutions of system (1) and using the results in [52], we get Using , we obtain Hence if , then . In addition, the equality holds if and only if , and . If , then . If , from (1), we get which implies that . Consequently, the largest invariant set of is the singleton . Therefore, by the LaSalle’s invariance principle [51], is globally asymptotically stable. The proof of the instability of is based on the computation of the Jacobean matrix of system (1) and the results presented in [53–55]. The Jacobean matrix of (1) at any equilibrium is given byWe recall that is locally asymptotically stable if the all eigenvalues of (40) satisfy the following condition [53–55]: From (40), the characteristic equation at is given as follows:where Obviously, (42) has the roots and . If , we have and . Then, there exists satisfying . In addition, we have . Consequently, when , is unstable. Theorem 4.  (i)The infection equilibrium without humoral immunity is globally asymptotically stable if , , and(ii)When , is unstable. Proof. Define a Lyapunov functional as follows:Calculating the fractional derivative of , we get Using , , , , and , we obtain Hence,Using the arithmetic-geometric inequality, we haveSince , we have if . It is easy to see that this condition is equivalent to (44). Furthermore, if and only if , and . We discuss two cases: If , then . If , from (1), we get , and then . Hence, the largest invariant set of is the singleton . By the LaSalle’s invariance principle, is globally asymptotically stable. At , the characteristic equation of (40) is given as follows:whereWe can easily see that (50) has the root . Then, when , we have . In this case, is unstable. Theorem 5. The infection equilibrium with humoral immunity is globally asymptotically stable if and Proof. Consider the following Lyapunov functional:Computing the fractional derivative of and using , , , , and , we get From (49), we have when . This condition is equivalent to (52). In addition, if , and . Further, ; then . Consequently, the largest invariant set of is the singleton . By the LaSalle’s invariance principle, is globally asymptotically stable. It is important to note that when is sufficiently small or is sufficiently large, the two conditions (44) and (52) are satisfied. Then, we have the following corollary. Corollary 6. Assume that . When is sufficiently small or is sufficiently large, then we have the following:(i)The infection equilibrium without humoral immunity is globally asymptotically stable if .(ii)The infection equilibrium with humoral immunity is globally asymptotically stable if . 4. Numerical Simulations In this section, we validate our theoretical results to HIV infection. Firstly, we take the parameter values as shown in Table 1. Table 1: Parameter values of system (1). By calculation, we have . Then system (1) has an infection-free equilibrium . By Theorem 3, the solution of (1) converges to (see Figure 1). Consequently, the virus is cleared and the infection dies out. Figure 1: Stability of the infection-free equilibrium . Now, we choose and we keep the other parameter values. Hence, we obtain , , and Consequently, condition (44) is satisfied. Therefore, the infection equilibrium without humoral immunity is globally asymptotically stable. Figure 2 demonstrates this result. In this case, the infection becomes chronic. Figure 2: Stability of the infection equilibrium without humoral immunity . Next, we take and do not change the other parameter values. In this case, we have , , and . Hence, condition (52) is satisfied. Consequently, system (1) has an infection equilibrium with humoral immunity which is globally asymptotically stable. Figure 3 illustrates this result. We can observe that the activation of the humoral immune response increases the healthy cells and decreases the productive infected cells and viral load to a lower levels but it is not able to eradicate the infection. Figure 3: Stability of the infection equilibrium with humoral immunity . 5. Conclusion In the present paper, we have studied the dynamics of a viral infection model by taking into account the memory effect represented by the Caputo fractional derivative and the humoral immunity. We have proved that the solutions of the model are nonnegative and bounded which assure the well-posedness. We have shown that the proposed model has three infection equilibriums, namely, the infection-free equilibrium , the infection equilibrium without humoral immunity , and the infection equilibrium with humoral immunity . By constructing suitable Lyapunov functionals, the global stability of these equilibria is fully determined by two threshold parameters and . More precisely, when , is globally asymptotically stable, whereas if , it becomes unstable and another equilibrium point appears, that is, , which is globally asymptotically stable whenever and condition (44) is satisfied. In the case that , becomes unstable and there exists another equilibrium point which is globally asymptotically stable when condition (52) is satisfied. In addition, we remarked that when is sufficiently small or is sufficiently large, conditions (44) and (52) are verified, and then the global stability of and is characterized only by and . From our theoretical and numerical results, we deduce that the order of the fractional derivative has no effect on the dynamics of the model. However, when the value of decreases (long memory), the solutions of our model converge rapidly to the steady states (see Figures 1–3). This behavior can be explained by the memory term included in the fractional derivative which represents the time needed for the interaction between cells and viral particles and the time needed for the activation of humoral immune response. In fact, the knowledge about the infection and the activation of the humoral immune response in an early stage can help us to control the infection. Data Availability The data used to support the findings of this study are available from the corresponding author upon request. Conflicts of Interest The authors declare that they have no conflicts of interest. References K. Hattaf, N. Yousfi, and A. Tridane, “Global stability analysis of a generalized virus dynamics model with the immune response,” Canadian Applied Mathematics Quarterly, vol. 20, no. 4, pp. 499–518, 2012. View at Google Scholar · View at MathSciNetM. Maziane, K. Hattaf, and N. Yousfi, “Global stability for a class of HIV infection models with cure of infected cells in eclipse stage and CTL immune response,” International Journal of Dynamics and Control, 2016. View at Publisher · View at Google ScholarX. Wang, Y. Tao, and X. Song, “Global stability of a virus dynamics model with Beddington-DeAngelis incidence rate and CTL immune response,” Nonlinear Dynamics, vol. 66, no. 4, pp. 825–830, 2011. View at Publisher · View at Google Scholar · View at MathSciNetC. Lv, L. Huang, and Z. Yuan, “Global stability for an HIV-1 infection model with Beddington-DeAngelis incidence rate and CTL immune response,” Communications in Nonlinear Science and Numerical Simulation, vol. 19, no. 1, pp. 121–127, 2014. View at Publisher · View at Google Scholar · View at MathSciNetM. A. Nowak and C. R. M. Bangham, “Population dynamics of immune responses to persistent viruses,” Science, vol. 272, no. 5258, pp. 74–79, 1996. View at Publisher · View at Google Scholar · View at ScopusA. Murase, T. Sasaki, and T. Kajiwara, “Stability analysis of pathogen-immune interaction dynamics,” Journal of Mathematical Biology, vol. 51, no. 3, pp. 247–267, 2005. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusM. A. Obaid and A. M. Elaiw, “Stability of Virus Infection Models with Antibodies and Chronically Infected Cells,” Abstract and Applied Analysis, vol. 2014, Article ID 650371, 12 pages, 2014. View at Publisher · View at Google Scholar · View at MathSciNetM. A. Obaid, “Dynamical behaviors of a nonlinear virus infection model with latently infected cells and immune response,” Journal of Computational Analysis and Applications, vol. 21, no. 1, pp. 182–193, 2016. View at Google Scholar · View at MathSciNet · View at ScopusH. F. Huo, Y. L. Tang, and L. X. Feng, “A Virus Dynamics Model with Saturation Infection and Humoral Immunity,” Int. Journal of Math. Analysis, vol. 6, no. 40, 2012. View at Google ScholarA. M. Elaiw, “Global stability analysis of humoral immunity virus dynamics model including latently infected cells,” Journal of Biological Dynamics, vol. 9, no. 1, pp. 215–228, 2015. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusA. M. Elaiw and N. H. AlShamrani, “Global properties of nonlinear humoral immunity viral infection models,” International Journal of Biomathematics, vol. 8, no. 5, 1550058, 53 pages, 2015. View at Publisher · View at Google Scholar · View at MathSciNetJ. X. Velasco-Herna'ndez, J. A. Garci'a, and D. E. Kirschner, “Remarks on modeling host-viral dynamics and treatment, Mathematical Approaches for Emerging and Reemerging Infectious Diseases,” An Introduction to Models, Methods, and Theory, vol. 125, pp. 287–308, 2002. View at Google ScholarA. S. Perelson, “Modelling viral and immune system dynamics,” Nature Reviews Immunology, vol. 2, no. 1, pp. 28–36, 2002. View at Publisher · View at Google Scholar · View at ScopusR. L. Magin, “Fractional calculus models of complex dynamics in biological tissues,” Computers & Mathematics with Applications, vol. 59, no. 5, pp. 1586–1593, 2010. View at Publisher · View at Google Scholar · View at ScopusA. A. Stanislavsky, “Memory effects and macroscopic manifestation of randomness,” Physical Review E: Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, vol. 61, no. 5, pp. 4752–4759, 2000. View at Google Scholar · View at ScopusM. Saeedian, M. Khalighi, N. Azimi-Tafreshi, G. R. Jafari, and M. Ausloos, “Memory effects on epidemic evolution: The susceptible-infected-recovered epidemic model,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 95, no. 2, Article ID 022409, 2017. View at Publisher · View at Google Scholar · View at ScopusY. A. Rossikhin and M. V. Shitikova, “Applications of fractional calculus to dynamic problems of linear and nonlinear hereditary mechanics of solids,” Applied Mechanics Reviews, vol. 50, no. 1, pp. 15–67, 1997. View at Publisher · View at Google Scholar · View at ScopusR. J. Marks and M. W. Hall, “Differintegral Interpolation from a Bandlimited Signal's Samples,” IEEE Transactions on Signal Processing, vol. 29, no. 4, pp. 872–877, 1981. View at Publisher · View at Google Scholar · View at ScopusG. L. Jia and Y. X. Ming, “Study on the viscoelasticity of cancellous bone based on higher-order fractional models,” in Proceedings of the 2nd International Conference on Bioinformatics and Biomedical Engineering (ICBBE '08), pp. 1733–1736, Shanghai, China, May 2008. View at Publisher · View at Google Scholar · View at ScopusR. Magin, “Fractional calculus in bioengineering,” Cretical reviews in biomedical engineering, vol. 32, pp. 13–77, 2004. View at Google ScholarE. Scalas, R. Gorenflo, and F. Mainardi, “Fractional calculus and continuous-time finance,” Physica A: Statistical Mechanics and its Applications, vol. 284, no. 1–4, pp. 376–384, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusR. Capponetto, G. Dongola, L. Fortuna, and I. Petras, “Fractional order systems: Modelling and control applications,” World Scientific Series in Nonlinear Science, Series A, vol. 72, 2010. View at Google ScholarV. D. Djordjević, J. Jarić, B. Fabry, J. J. Fredberg, and D. Stamenović, “Fractional derivatives embody essential features of cell rheological behavior,” Annals of Biomedical Engineering, vol. 31, no. 6, pp. 692–699, 2003. View at Publisher · View at Google Scholar · View at ScopusK. S. Cole, “Electric conductance of biological systems,” Cold Spring Harbor Symposium on Quantitative Biology, vol. 1, pp. 107–116, 1933. View at Publisher · View at Google ScholarF. A. Rihan, M. Sheek-Hussein, A. Tridane, and R. Yafia, “Dynamics of hepatitis C virus infection: mathematical modeling and parameter estimation,” Mathematical Modelling of Natural Phenomena, vol. 12, no. 5, pp. 33–47, 2017. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusN. Khodabakhshi, S. M. Vaezpour, and D. Baleanu, “On Dynamics Of a Fractional-order Model Of HCV Infection,” Journal of Mathematical Analysis, vol. 8, no. 1, pp. 16–27, 2017. View at Google Scholar · View at MathSciNetX. Zhou and Q. Sun, “Stability analysis of a fractional-order HBV infection model,” International Journal of Advances in Applied Mathematics and Mechanics, vol. 2, no. 2, pp. 1–6, 2014. View at Google Scholar · View at MathSciNetS. S. M. and Y. A. M., “On a fractional-order 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 · View at Google Scholar · View at MathSciNetY. Ding and H. Ye, “A fractional-order differential equation model of HIV infection of CD4+ T-cells,” Mathematical and Computer Modelling, vol. 50, no. 3-4, pp. 386–392, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at ScopusC. M. Pinto and A. R. Carvalho, “The role of synaptic transmission in a HIV model with memory,” Applied Mathematics and Computation, vol. 292, pp. 76–95, 2017. View at Publisher · View at Google Scholar · View at MathSciNetA. Boukhouima, K. Hattaf, and N. Yousfi, “Dynamics of a Fractional Order HIV Infection Model with Specific Functional Response and Cure Rate,” International Journal of Differential Equations, vol. 2017, Article ID 8372140, 8 pages, 2017. View at Publisher · View at Google Scholar · View at MathSciNetJ. A. Deans and S. Cohen, “Immunology of malaria,” Annual Review of Microbiology, vol. 37, pp. 25–49, 1983. View at Publisher · View at Google Scholar · View at ScopusI. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and Some of their Applications, vol. 198, Academic press, 1998. K. Hattaf and N. Yousfi, “A class of delayed viral infection models with general incidence rate and adaptive immune response,” International Journal of Dynamics and Control, vol. 4, no. 3, pp. 254–265, 2016. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusD. Riad, K. Hattaf, and N. Yousfi, “Dynamics of Capital-labour Model with Hattaf-Yousfi Functional Response,” British Journal of Mathematics & Computer Science, vol. 18, no. 5, pp. 1–7, 2016. View at Publisher · View at Google ScholarM. Mahrouf, K. Hattaf, and N. Yousfi, “Dynamics of a stochastic viral infection model with immune response,” Mathematical Modelling of Natural Phenomena, vol. 12, no. 5, pp. 15–32, 2017. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusK. Hattaf, N. Yousfi, and A. Tridane, “Stability analysis of a virus dynamics model with general incidence rate and two delays,” Applied Mathematics and Computation, vol. 221, pp. 514–521, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusP. Crowley and E. Martin, “Functional responses and interference within and between year classes of a dragonfly population,” Journal of the North American Benthological Society, vol. 8, pp. 211–221, 1989. View at Publisher · View at Google ScholarS. Xu, “Global stability of the virus dynamics model with Crowley-Martin functional response,” Electronic Journal of Qualitative Theory of Differential Equations, No. 9, 10 pages, 2012. View at Publisher · View at Google Scholar · View at MathSciNetJ. R. Beddington, “Mutual interference between parasites or predat ors and its effect on searching efficiency,” Journal of Animal Ecology, vol. 44, pp. 331–340, 1975. View at Publisher · View at Google ScholarD. L. DeAngelis, R. A. Goldstein, and R. V. ONeill, “A model for trophic interaction,” Ecology, pp. 881–892, 1975. View at Google ScholarL. Rong, M. A. Gilchrist, Z. Feng, and A. S. Perelson, “Modeling within-host HIV-1 dynamics and the evolution of drug resistance: trade-offs between viral enzyme function and drug susceptibility,” Journal of Theoretical Biology, vol. 247, no. 4, pp. 804–818, 2007. View at Publisher · View at Google Scholar · View at MathSciNetZ. Hu, W. Pang, F. Liao, and W. Ma, “Analysis of a cd4+ t cell viral infection model with a class of saturated infection rate,” Discrete and Continuous Dynamical Systems - Series B, vol. 19, no. 3, pp. 735–745, 2014. View at Publisher · View at Google Scholar · View at ScopusM. Maziane, E. M. Lotfi, K. Hattaf, and N. Yousfi, “Dynamics of a Class of HIV Infection Models with Cure of Infected Cells in Eclipse Stage,” Acta Biotheoretica, vol. 63, no. 4, pp. 363–380, 2015. View at Publisher · View at Google Scholar · View at ScopusJ. K. Hale and S. M. Verduyn Lunel, Introduction to Functional-Differential Equations, Springer, Berlin, Germany, 1993. View at Publisher · View at Google Scholar · View at MathSciNetJ. K. Hale, “Sufficient conditions for stability and instability of autonomous functional-differential equations,” Journal of Differential Equations, vol. 1, pp. 452–482, 1965. View at Publisher · View at Google Scholar · View at MathSciNet · View at ScopusT. A. Burton, Stability and Periodic sSolutions of Ordinary and Functional Differential Equations, Academic Press, Orlando, Fla, USA, 1985. View at MathSciNetB. S. Ogundare, “Stability and boundedness properties of solutions to certain fifth order nonlinear differential equations,” Matematicki Vesnik, vol. 61, no. 4, pp. 257–268, 2009. View at Google Scholar · View at MathSciNetC. Tunç, “A study of the stability and boundedness of the solutions of nonlinear differential equations of the fifth order,” Indian Journal of Pure and Applied Mathematics, vol. 33, no. 4, pp. 519–529, 2002. View at Google Scholar · View at MathSciNetC. Tunç, “New results on the stability and boundedness of nonlinear differential equations of fifth order with multiple deviating arguments,” Bulletin of the Malaysian Mathematical Sciences Society, vol. 36, no. 3, pp. 671–682, 2013. View at Google Scholar · View at MathSciNetJ. 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 · View at Google Scholar · View at MathSciNetC. Vargas-De-Leon, “Volterra-type Lyapunov functions for fractional-order epidemic systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 24, no. 1–3, pp. 75–85, 2015. View at Publisher · View at Google Scholar · View at MathSciNetD. Matignon, “Stability results for fractional differential equations with applications to control processing,” in Computational Eng. in Sys. Appl, vol. 2, p. 963, Lille, France, 1996. View at Google ScholarE. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “On some Routh-Hurwitz 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 · View at Google Scholar · View at Zentralblatt MATH · View at ScopusE. Ahmed, A. M. A. El-Sayed, and H. A. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,” Journal of Mathematical Analysis and Applications, vol. 325, no. 1, pp. 542–553, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus


This is a preview of a remote PDF: http://downloads.hindawi.com/journals/ijde/2018/1019242.pdf

Adnane Boukhouima, Khalid Hattaf, Noura Yousfi. A Fractional Order Model for Viral Infection with Cure of Infected Cells and Humoral Immunity, International Journal of Differential Equations, 2018, DOI: 10.1155/2018/1019242