Computational and Mathematical Methods in Medicine
Volume 2019 (2019), Article ID 4059549, 12 pages
https://doi.org/10.1155/2019/4059549

Research Article

The Fractional Differential Model of HIV-1 Infection of CD4⁺ T-Cells with Description of the Effect of Antiviral Drug Treatment

Bijan Hasani Lichae,1,2 Jafar Biazar,³ and Zainab Ayati⁴

¹Department of Mathematics, Guilan Science and Research Branch, Islamic Azad University, Rasht, Iran
²Department of Mathematics, Rasht Branch, Islamic Azad University, Rasht, Iran
³Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, P.O. Box 41335-1914, Rasht, Guilan, Iran
⁴Department of Engineering Sciences, Faculty of Technology and Engineering East of Guilan, University of Guilan, P.C. 44891-63157, Rudsar-Vajargah, Iran

Correspondence should be addressed to Jafar Biazar; biazar@guilan.ac.ir

Received 13 June 2018; Revised 11 September 2018; Accepted 17 December 2018; Published 8 January 2019

Guest Editor: Katarzyna Kolasa

Copyright © 2019 Bijan Hasani Lichae 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, the fractional-order differential model of HIV-1 infection of CD4⁺ T-cells with the effect of drug therapy has been introduced. There are three components: uninfected CD4⁺ T-cells, , infected CD4⁺ T-cells, , and density of virions in plasma, . The aim is to gain numerical solution of this fractional-order HIV-1 model by Laplace Adomian decomposition method (LADM). The solution of the proposed model has been achieved in a series form. Moreover, to illustrate the ability and efficiency of the proposed approach, the solution will be compared with the solutions of some other numerical methods. The Caputo sense has been used for fractional derivatives.

1. Introduction

Human immunodeficiency virus (HIV) is a retrovirus that causes acquired immunodeficiency syndrome (AIDS) [1]. HIV infects, damages, and reduces CD4⁺ T-cells. Therefore, it causes to decrease the resistance of immune system [2]. The body becomes more gradually sensitive to infections and loses its safety. AIDS is one the most important and dangerous diseases in our time. According to UNAIDS 2017 annual report, “36.7 million people globally were living with HIV and 1.8 million people became newly infected with HIV and 1 million people died from AIDS-related illnesses in 2016.” In spite of the great progress in controlling the disease, no vaccine has been yet discovered for HIV. In the last two decades, a lot of efforts have been made to design and solve mathematical models that have essential rule in analyzing to control and prevent the spread of HIV-related diseases [3–13]. Usually almost all of these mathematical models explain the relation between HIV viruses and uninfected CD4⁺ cells and the effect of drug therapy to infected cells. Bonhoeffer et al. [4] presented a model for virus dynamics with two components and , where denotes the density of infected cells and shows the density of virus-producing cells.

The proposed model is as follows:

where is the rate of production of infected cells, is the natural death rate of infected cells, is the rate of virus-producing cells’ death, and is the rate of infection of uninfected cells. This model and many such models were inspired from Anderson’s model [14, 15]. Anderson’s model is one of the first and the most important models of infectious diseases. Tuckwell and Wan [16] introduced a modified model of equation (1) with three components: uninfected, infected CD4⁺ T-cells, and density of virions in plasma (, , and , respectively). The presented model with three equations is as follows:

with the initial conditions ,, and where , , , , , and are constant coefficients, is the rate of creation or production of CD4⁺ T-cells, is the natural death rate, is the rate of infected CD4⁺ cells from uninfected CD4⁺ cells, is the rate of virus-producing cells’ death, is the rate of creation of virions viruses by infected cells, and is the rate of virus particle death. For the sake of comparison and showing the ability of the proposed approach, we use the parameter values reported in references [6, 16]. The parameter values are as follows: (day/mm³), (day/mm³), (day/virion/mm³), (day/mm³), (virion/CLM/day), and (day). The rate of some coefficients will change if drug therapy is not 100% successful. When the drug treatment begins, infected cells which create virus components are affected. If the drug therapy is not effective, a part of infected cells will improve and remaining cells will begin to produce virus [8].

Mathematical modeling of many problems in biology and other branches of sciences appears as differential equations in fractional order. Because the fractional-order differential equations save memory on themselves and are related to fractals [8, 17–19], we prefer to use the fractional-order form of the model (2) as follows:

with the same initial conditions, where , . There are some numerical approaches for solving these types of mathematical models. Some of these methods are as follows: homotopy analysis, optimal homotopy asymptotic, homotopy perturbation, Adomian decomposition, and variational iteration [20–30]. In [8], system of fractional equation (3) has been solved by homotopy analysis method (HAM) and generalized Euler method (GEM). In [31], equation (3) has also been solved by homotopy perturbation method (HPM). Adomian [32], introduced a decomposition method (ADM) which is a powerful method to get analytic approximate solution of differential equations. Using Laplace transform method with couple of ADM (LADM) to solve systems of differential equations leads to an effective method that finds many applications in applied mathematics. In this paper, we will solve equation (3) by LADM and will compare the results with the results achieved by generalized Euler, homotopy analysis, homotopy perturbation, and Runge–Kutta methods. The structure of the paper is as follows: in Section 2, we will present a brief review of fractional calculus. In Sections 3 and 4, we will solve the fractional-order HIV-1 model by LADM. In Section 5, the convergence of the method will be discussed. In the last section, we present the conclusion.

2. Fractional Calculus

The purpose of this section is to recall a few preliminaries about what appears in this research.

Definition 2.1. The Riemann–Liouville fractional integral of order for a function is defined as

where (see [33]).

Definition 2.2. The Caputo fractional derivative for a function on the closed interval is defined as

where is the integer part of . Another presentation of the Caputo fractional derivative can be shown as follows (see [33]):

Lemma 2.1. If and , then the following result holds for fractional calculus:

Proof. See [33, 34].

Definition 2.3. The Laplace transform of Caputo fractional derivative is defined as follows:

3. Solution of Model (3)

In this section, LADM has been implemented to solve system of fractional equation (3) with the initial conditions.

We apply Laplace transform on both sides of each equation of equation (3):

which implies that

Substitution of the initial conditions in equation (10) and applying inverse Laplace transform results in

To apply ADM, let us consider , , and as the following series:

To decompose the nonlinear term , let us follow an alternate algorithm [35] to get,

where is as the following equation:

substituting equations (12)–(14) into (11) reads

We take inverse Laplace transform on both sides of each equation of equation (15):

We have calculated four terms of the infinite series of , , and as an approximate solution. To get any desired accuracy, one is able to proceed the process and obtain more terms. Finally, the solution of mathematical model can be obtained as follows:

4. Numerical Simulation

In this section, constants and initial values are substituted in equation (16) to obtain an approximate solution.

Substituting the following values:  = 0.272 (day/mm³),  = 0.00136 (day/mm³),  = 0.00027 (day/virion/mm³),  = 0.33 (day/mm³),  = 50 (virion/CLM/day), and  = 2.0 (day) and the initial conditions , and in equation (16), we get,

Three terms approximations can be written as the following form:

Let us take , , and equal to , so the solution of fractional-order of model (3) is obtained as follows:

For , the solution of equation (3) will be as follows:

In Tables 1–3, one can compare the approximate solution of fractional-order of model (3) with the results of GEM, HAM, RK4 in [8], and HPM in [31] using traditional order . The results of LADM are more accurate than the results obtained by other methods.

Table 1: Numerical results of (uninfected CD4⁺ T-cells).

t
LADM
GEM
HPM
HAM
RK4

0
100
100
100
100
100

0.2
100.023
100.023
100.023
100.023
100.023

0.4
100.046
100.047
100.047
100.047
100.047

0.6
100.070
100.071
100.071
100.071
100.071

0.8
100.092
100.097
100.097
100.096
100.097

1.0
100.112
100.122
100.123
100.122
100.122

Table 2: Numerical results of (infected CD4⁺ T-cells).

t
LADM
GEM
HPM
HAM
RK4

0
0
0
0
0
0

0.2
0.00436
0.00434
0.00434
0.00434
0.00434

0.4
0.00750
0.00715
0.00721
0.00714
0.00715

0.6
0.01074
0.00908
0.00934
0.00909
0.00908

0.8
0.01336
0.01049
0.01117
0.01063
0.01049

1.0
0.01866
0.01161
0.01276
0.01194
0.01161

Table 3: Numerical results of (density of virions in plasma).

t
LADM
GEM
HPM
HAM
RK4

0
1
1
1
1
1

0.2
0.69581
0.69030
0.69071
0.69059
0.69070

0.4
0.53851
0.51152
0.51208
0.51237
0.51190

0.6
0.46097
0.41069
0.41394
0.40994
0.41103

0.8
0.39607
0.35656
0.37749
0.35148
0.35684

1.0
0.27671
0.33053
0.42419
0.32869
0.33073

Figures 1–3 show the results for different values of , and the results can be compared.

Figure 1: Dynamics of uninfected CD4⁺ T-cells for various values of .

Figure 2: Dynamics of infected CD4⁺ T-cells for various values of .

Figure 3: Dynamics of density of virions in plasma for various values of .

5. Convergence Analysis of the Method

In this section, convergence of the proposed method, using the idea presented in [36], is studied.

Consider the following functional equation:

where is a functional operator and can be decomposed as and is a known function.

is a Caputo fractional derivative operator, is a linear operator, and is a nonlinear analytic operator, respectively. So equation (22) can be written as follows:

The goal is to find a function satisfying equation (22). By applying the Laplace transform on both sides of equation (23) reads

By using definition (2.3), equation (23) can be written as follows:

by considering and using inverse of Laplace transform on both sides of equation (25) results in

which implies that

By implementing ADM and assuming the solution as an infinite series say, , and writing the nonlinear term based on Adomian polynomials such as

where

Equation (27) can be written as the following form:

So, we have

From which we can define

Theorem 5.1. LADM for equation (21) with the solution (19) is equivalent to

By using the following iterative scheme:

where

Proof. For , from equation (33):

Then, by assumption of equation (33), we have

for :

We know that , so we obtain

By strong induction, let us have

and prove the following for ,

Then, from equation (33), we derive

This entails the statement is true and the theorem is proved.

Theorem 5.2. Let be a Banach space.

(i) resulted from equation (31), convergence to , if (ii) satisfies in

Proof.

, we have

This means that , therefore, is a Cauchy sequence in the Banach space of and is convergent. So, .
From equation (33), we derive

From equation (35), we have

So,

Lemma 5.1. Equation (43) is equivalent to equation (22).

Proof. By using Laplace transform on both sides of equation (43) reads to

so we can write

In virtue of definition 2.3 and linearity of the Laplace transform, equation (50) can be written as the follows:

By applying the inverse of Laplace transform on both sides of equation (51), we derive

Considering , one gets equation (22). So, the solution of equation (43) is the same as the solution of equation (22).

6. Conclusion

In this paper, a fractional-order model of HIV-1 with three components has been introduced. When , then ; therefore, the fractional-order of presented model reduces to traditional model. By applying Laplace transform and Adomian decomposition method, or LADM for short, which is a strong approach to compute numerical solution of fractional differential equations, we gain an approximate solution of the proposed model. The accuracy of the proposed approach has made it a reliable method. We have calculated four terms of the infinite series of , , and as an approximate solution. The result of LADM has been compared with the results of some other methods such as GEM, HAM, RK4 [8], and HPM [31]. The results are presented in Tables 1–3. Figures 1–3, show that the uninfected CD4⁺ T-cells, , infected CD4⁺ T-cells, , and density of virions in plasma, , depend on the various values of , the various values of the parameters, and the time fractional derivative. A comparison of the approximate solutions shows that LADM can work more accurate than other methods. Convergence of the proposed method is studied. Because of the fact that obtaining the exact solution for system of fractional equation is difficult or impossible, we would like to suggest such an easy and reliable approach for further research, in the future.

Data Availability

Data used to support this study are available at DOI: https://doi.org/10.1038/sj.icb.7100056. These prior studies (and datasets) are cited at relevant places within the text as references [6, 8, 16].

Disclosure

This research did not receive specific funding, but has resulted from a Ph.D. program, by Mr. Bijan Hasani Lichaee, under supervision of Prof. Jafar Biazar and cosupervision of Dr. Zainab Ayati.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

R. Weiss, “How does HIV cause AIDS?” Science, vol. 260, no. 5112, pp. 1273–1279, 1993.

D. C. Douek, M. Roederer, and R. A. Koup, “Emerging concepts in the immunopathogenesis of AIDS,” Annual Review of Medicine, vol. 60, no. 1, pp. 471–484, 2009.

A. V. Herz, S. Bonhoeffer, R. M. Anderson, R. M. May, and M. A. Nowak, “Viral dynamics in vivo: limitations on estimates of intracellular delay and virus decay,” Proceedings of the National Academy of Sciences, vol. 93, no. 14, pp. 7247–7251, 1996.

S. Bonhoeffer, J. M. Coffin, and M. A. Nowak, “Human immunodeficiency virus drug therapy and virus load,” Journal of Virology, vol. 71, no. 4, pp. 3275–3278, 1997.

R. V. Culshaw and S. Ruan, “A delay-differential equation model of HIV infection of CD4⁺ T-cells,” Mathematical biosciences, vol. 165, no. 1, pp. 27–39, 2000.

H. C. Tuckwell and F. Y. M. Wan, “On the behavior of solutions in viral dynamical models,” BioSystems, vol. 73, no. 3, pp. 157–161, 2004.

C. Lv and Z. Yuan, “Stability analysis of delay differential equation models of HIV-1 therapy for fighting a virus with another virus,” Journal of Mathematical Analysis and Applications, vol. 352, no. 2, pp. 672–683, 2009.

A. A. M. Arafa, S. Z. Rida, and M. Khalil, “The effect of anti-viral drug treatment of human immunodeficiency virus type 1 (HIV-1) described by a fractional order model,” Applied Mathematical Modelling, vol. 37, no. 4, pp. 2189–2196, 2013.

A. A. M. Arafa, S. Z. Rida, and M. Khalil, “Fractional modeling dynamics of HIV and CD4+ T-cells during primary infection,” Nonlinear biomedical physics, vol. 6, no. 1, p. 1, 2012.

Y. Yan and C. Kou, “Stability analysis for a fractional differential model of HIV infection of CD4⁺ T-cells with time delay,” Mathematics and Computers in Simulation, vol. 82, no. 9, pp. 1572–1585, 2012.

L. Wang and M. Y. Li, “Mathematical analysis of the global dynamics of a model for HIV infection of CD4⁺ T cells,” Mathematical Biosciences, vol. 200, no. 1, pp. 44–57, 2006.

A. S. Perelson, D. E. Kirschner, and R. De Boer, “Dynamics of HIV infection of CD4⁺ T cells,” Mathematical Biosciences, vol. 114, no. 1, pp. 81–125, 1993.

J. Biazar, “Solution of the epidemic model by Adomian decomposition method,” Applied Mathematics and Computation, vol. 173, no. 2, pp. 1101–1106, 2006.

R. M. Anderson and R. M. May, “Population biology of infectious diseases: part I,” Nature, vol. 280, no. 5721, pp. 361–367, 1979.

R. M. May and R. M. Anderson, “Population biology of infectious diseases: part II,” Nature, vol. 280, no. 5722, pp. 455–461, 1979.

H. C. Tuckwell and F. Y. M. Wan, “Nature of equilibria and effects of drug treatments in some simple viral population dynamical models,” IMA Journal of Mathematical Control and Information, vol. 17, no. 4, pp. 311–327, 2000.

N. Heymans and J.-C. Bauwens, “Fractal rheological models and fractional differential equations for viscoelastic behavior,” Rheologica Acta, vol. 33, no. 3, pp. 210–219, 1994.

M. Giona and H. E. Roman, “Fractional diffusion equation on fractals: one-dimensional case and asymptotic behaviour,” Journal of Physics A: Mathematical and General, vol. 25, no. 8, pp. 2093–2105, 1999.

F. B. Tatom, “The relationship between fractional calculus and fractals,” Fractals, vol. 3, no. 1, pp. 217–229, 1995.

S. J. Liao, “The proposed homotopy analysis technique for the solution of nonlinear problems,” Shanghai Jiao Tong University, Shanghai, China, 1992, Ph.D. thesis.

V. S. Ertürk, Z. M. Odibat, and S. Momani, “An approximate solution of a fractional order differential equation model of human T-cell lymphotropic virus I (HTLV-I) infection of CD4⁺ T-cells,” Computers and Mathematics with Applications, vol. 62, no. 3, pp. 996–1002, 2011.

Y. 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.

F. Haq, K. Shah, G. U. Rahman et al., “Numerical solution of fractional order smoking model via laplace Adomian decomposition method,” Alexandria Engineering Journal, vol. 57, no. 2, pp. 1061–1069, 2017.

M. Merdan, A. Gökdoğan, and A. Yildirim, “On the numerical solution of the model for HIV infection of CD4⁺ T cells,” Computers and Mathematics with Applications, vol. 62, no. 1, pp. 118–123, 2011.

Y. Khan, H. Vazquez-Leal, and Q. Wu, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications, vol. 23, no. 3-4, pp. 677–682, 2013.

Ş. Yüzbaşı, “A numerical approach to solve the model for HIV infection of CD4⁺ T cells,” Applied Mathematical Modelling, vol. 36, no. 12, pp. 5876–5890, 2012.

M. Ghoreishi, A. M. Ismail, and A. Alomari, “Application of the homotopy analysis method for solving a model for HIV infection of CD4⁺ T-cells,” Mathematical and Computer Modelling, vol. 54, no. 11-12, pp. 3007–3015, 2011.

A. Atangana and E. F. D. Goufo, “On the mathematical analysis of Ebola hemorrhagic fever: deathly infection disease in West African countries,” BioMed Research International, vol. 2014, Article ID 261383, 7 pages, 2014.

M. Y. Ongun, “The Laplace adomian decomposition method for solving a model for HIV infection of CD4⁺ T cells,” Mathematical and Computer Modelling, vol. 53, no. 5-6, pp. 597–603, 2011.

A. A. M. Arafa, S. Rida, and M. Khalil, “A fractional-order model of HIV infection with drug therapy effect,” Journal of the Egyptian Mathematical Society, vol. 22, no. 3, pp. 538–543, 2014.

M. Merdan and T. Khan, “Homotopy perturbation method for solving viral dynamical model,” CÜ Fen-Edebiyat Fakültesi, Fen Bilimleri Dergisi, vol. 31, pp. 65–77, 2010.

G. Adomian, Nonlinear Stochastic Systems Theory and Applications to Physics, vol. 46, Springer Science & Business Media, Berlin, Germany, 1988.

K. Diethelm, The Analysis of Fractional Differential Equations: An Application-Oriented Exposition using Differential Operators of Caputo Type, Springer, Berlin, Germany, 2010.

A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Preface, Elsevier, New York City, NY, USA, 2006.

J. Biazar, E. Babolian, G. Kember, A. Nouri, and R. Islam, “An alternate algorithm for computing Adomian polynomials in special cases,” Applied Mathematics and Computation, vol. 138, no. 2-3, pp. 523–529, 2003.

Z. Ayati and J. Biazar, “On the convergence of Homotopy perturbation method,” Journal of the Egyptian Mathematical Society, vol. 23, no. 2, pp. 424–428, 2015.