Computational and Theoretical Analysis of Human Diseases Associated with Infectious Pathogens
View this Special IssueResearch Article  Open Access
Abdon Atangana, Emile Franc Doungmo Goufo, "Computational Analysis of the Model Describing HIV Infection of CD4^{+}T Cells", BioMed Research International, vol. 2014, Article ID 618404, 7 pages, 2014. https://doi.org/10.1155/2014/618404
Computational Analysis of the Model Describing HIV Infection of CD4^{+}T Cells
Abstract
An analysis of the model underpinning the description of the spread of HIV infection of CD4^{+}T cells is examined in detail in this work. Investigations of the disease free and endemic equilibrium are done using the method of Jacobian matrix. An iteration technique, namely, the homotopy decomposition method (HDM), is implemented to give an approximate solution of nonlinear ordinary differential equation systems. The technique is described and illustrated with numerical examples. The approximated solution obtained via HDM is compared with those obtained via other methods to prove the trustworthiness of HDM. Moreover, the lessening and simplicity in calculations furnish HDM with a broader applicability.
1. Introduction
The cause of the progressive depletion of CD4^{+}T cells in HIVinfected people is one of the most fundamental and controversial issues in AIDS research. HIV infects and kills CD4^{+}T cells. The infection results in high Tcell activation and turnover. An immediately intuitive assumption is that HIVmediated destruction of CD4^{+} cells directly reduces the number of these cells and that the high turnover rates of T cells and the slow progression to AIDS reflect a long but eventually lost struggle of the immune system to replace killed cells in its effort to maintain Tcell homeostasis [1ā4]. However, HIV mainly infects activated CD4^{+} cells, and activated cells normally follow different dynamics than cells that belong to resting populations whose numbers are controlled by homeostatic mechanisms.
In this study, we consider that the HIV infection model of CD4^{+}T cells is examined in [5]. This model is given by the components of the basic threecomponent model which are the concentration of susceptible CD4^{+}T cells, CD4^{+}T cells infected by the HIV viruses, and free HIV virus particles in the blood. CD4^{+}T cells are also called leukocytes or T helper cells. These with order cells in human immunity systems fight against diseases. HIV use cells in order to propagate. In a healthy person, the number of CD4^{+}T cells is 800/1200āmm^{3}. This model is characterized by a system of the nonlinear differential equations: subject to the initial conditions Here, is any positive constant, , , and show the concentration of susceptible CD4^{+}T cells, CD4^{+}T cells infected by the HIV viruses, and free HIV virus particles in the blood, respectively, , , and stand for natural turnover rates of uninfected Tells, infected T cells, and virus particles, respectively, describes the logistic growth of the healthy CD4^{+}T cells, and proliferation of infected CD4^{+}T cells is neglected. For is the infection rate, the term describes the incidence of HIV infection of healthy CD4^{+}T cells. Each infected CD4^{+}T cell is assumed to produce l virus particles during its lifetime, including any of its daughter cells. The body is believed to produce CD4^{+}T cells from precursors in the bone marrow and thymus at a constant rate . T cells multiply through mitosis with a rate when T cells are stimulated by antigen or mitogen. denotes the maximum CD4^{+}T cell concentration in the body [6ā9].
As Andrianov and Manevitch wrote in the foreword of a popular science book Asymtotology: Ideas, Methods, and Applications, asymptotic methods belong to the, perhaps, most romantic area of modern mathematics [10]. Though computer science is growing very fast and numerical simulation is applied everywhere, nonnumerical issues will still play a large role [11ā14]. There exist some alternative analytical asymptotic approaches, such as the nonperturbative method, modified LindstedtPoincarĆ© method [13], variational iteration method [14], Adomian decomposition method [15], homotopy perturbation method [16, 17], and bookkeeping artificial parameter perturbation method [18].
The purpose of this paper is to derive analytical solutions of model for HIV infection of CD4^{+}T cells (1) via the relatively new analytical solution homotopy decomposition method (HDM). The HDM was recently used to solve one and twodimensional fractional heatlike problem, Burgers equation, the KleinGordon equation, a coupled Burgers equation, 2D and 3D Poisson equations and biharmonic equations, the groundwater flow equation [19], the HamiltonJacobiBellman equation, high evenorder differential equation, parabolic equations and coupled Van der Polās nonlinear partial differential equation, and coupled attractor onedimensional KellerSegel equations. This method displays some advantages over the existing methods.
2. Stability Analysis
The aim of this section is to present a stability analysis of the system equation (1) that will be achieved via the use the eigenvalues obtained via the Jacobian matrix. We will start by providing the equilibrium point and also the desease free equilibrium. To obtain the equilibrium points, we assume that the system does not depend on the parameter ; this will further imply that After the first manipulations, we obtain that From the above equation, we can obtain However, replacing (5) into the first equation of system (3) and rearranging it, we obtain the following: Therefore, the equilibrium points are given as The above equilibrium points are valid if the following conditions are satisfied: The disease free equilibrium is obtained by solving the following equation: which has the following solution: Therefore, the diseases free equilibrium is given as
With the above information on hand, we will now find the eigenvalues associated with this problem, which will allow us to give stability of the system. The Jacobian matrix associated with this problem is given by Now, using the equilibrium free point, we obtain the following: In order to find the eigenvalues associate, we will solve the following equation: Using the standard method, we obtain the following solution: With the Eigen values in hand, we conclude that there is stability if and only if the following conditions are observed:
3. Some Useful Information regarding the Methodology of (HDM)
To illustrate the basic idea of this method, we consider a general nonlinear nonhomogeneous differential equation with initial conditions of the following form (see also in [16, 19]): Subjected to the initial condition, is the order of the derivative.
is a known function, is the general nonlinear differential operator, represents a linear differential operator, and is the order of the derivative. The methodās first step here is to apply the inverse operator on both sides of (17) to obtain The multiintegral in (17) can be transformed to so that (19) can be reformulated as Using the homotopy scheme, the solution of the above integral equation is given in series form as follows: and the nonlinear term can be decomposed as where is an embedding parameter. is the Heās polynomials that can be generated by The homotopy decomposition method is obtained by the graceful coupling of decomposition method with Heās polynomials and is given by with Comparing the terms of the same powers of , we obtain solutions of various orders. The initial guess of the approximation is ; this is actually the Taylor series of the exact solution of order . Note that this initial guess insures the uniqueness of the series decompositions.
4. Applications
In this section, we applied the HDM to solve the system of the nonlinear differential equations (1). To be consistent in comparison with the existing methods, we chose the following initial condition and parameters as in [20, 21]: Following carefully the steps of the HDM, we obtain the following equations: Comparing the terms of the same power of , we obtain the following integral equations that are much easier to solve. Note that, with the homotopy perturbation method (HPM), one will obtain a set of ordinary differential equations, after comparing the terms of the same power of , which something is very heavy to compute in the case of high order ODE. Consider We will give the general algorithm in order to accommodate scholars using computers program; the algorithm will then be used to derive special solution of the system of equation for a given set of theoretical parameters.
Algorithm 1. (i)input: , and as initial guest;(ii)number terms in the rough calculation;(iii)output: the approximate solutions.
Step 1. Put , , , .
Step 2. From to , do Steps 3 and Step 4.
Step 3. Compute
Step 4. Compute , , and āāāā.
Stop.
We will now make use of the above algorithm to derive the special solution. We therefore obtain the following series solutions: Using the package Mathematica, in the same manner, one can obtain the rest of the components. But, in this case, 9 terms were computed and the asymptotic solution is given by
4.1. Numerical Applications
To test the effectiveness and the accuracy of the HDM for solving this type of problem, we compare the approximated solutions obtained via other methods and HDM, and the results are shown in Tables 1, 2, and 3.



Numerical comparison shows that the approximated solutions obtained via LADMPadĆ© [20] are in good agreement with the results obtained via HDM. More precisely, the approximated solutions are exactly the same. However, in the technique used in [20], one needs first to apply the Laplace transform on the system, following by the ADM, and finally take the inverse Laplace transform to obtain the approximated solutions, which is much time consuming and sometimes can lead to a very difficult situation, for example, if the inverse Laplace transform cannot be obtained.
5. Conclusions
In this paper, homotopy decomposition method has been developed for finding approximate solutions of HIV infection model of CD4^{+}T which is a class of nonlinear ordinary differential equation systems. We have demonstrated the accuracy and efficiency of the present technique with an example and comparison of the approximate solution obtained via the technique with those obtained with other methods. Comparing the methodology HDM to homotopy perturbation method (HPM), Adomian decomposition method (ADM), variational iteration method (VIM), and homotopy analysis method (HAM) has advantages. Disparate the ADM, the HDM is free from the need of Adomian polynomials. In this method, we do not need the Lagrange multiplier, correction functional, stationary conditions, or calculation of heavy integrals; the solution obtained is noise free, which eliminate the complications that exist in the VIM [23]. In contract to HPM, we do not need to continuously deform a difficult problem to another that is easier to solve. In contrast to LADM, the technique does not need any change of space, which sometime can be very difficult situation to handle, for example, in case the inverse Laplace transform cannot be found. We can easily conclude that the homotopy decomposition method is an efficient tool to solve approximate solution of nonlinear system partial differential equations.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 D. D. Ho, A. U. Neumann, A. S. Perelson, W. Chen, J. M. Leonard, and M. Markowitz, āRapid turnover of plasma virions and CD4 lymphocytes in HIV1 infection,ā Nature, vol. 373, no. 6510, pp. 123ā126, 1995. View at: Publisher Site  Google Scholar
 H. Mohri, S. Bonhoeffer, S. Monard, A. S. Perelson, and D. D. Ho, āRapid turnover of T lymphocytes in SIVinfected rhesus macaques,ā Science, vol. 279, no. 5354, pp. 1223ā1227, 1998. View at: Publisher Site  Google Scholar
 X. Wei, S. K. Ghosh, M. E. Taylor et al., āViral dynamics in human immunodeficiency virus type 1 infection,ā Nature, vol. 373, no. 6510, pp. 117ā122, 1995. View at: Publisher Site  Google Scholar
 A. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, and D. D. Ho, āHIV1 dynamics in vivo: virion clearance rate, infected cell lifespan, and viral generation time,ā Science, vol. 271, no. 5255, pp. 1582ā1586, 1996. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 A. S. Perelson and P. W. Nelson, āMathematical analysis of HIV1 dynamics in vivo,ā SIAM Review, vol. 41, no. 1, pp. 3ā44, 1999. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 B. Asquith and C. R. M. Bangham, āThe dynamics of Tcell fratricide: application of a robust approach to mathematical modelling in immunology,ā Journal of Theoretical Biology, vol. 222, no. 1, pp. 53ā69, 2003. View at: Publisher Site  Google Scholar  MathSciNet
 M. A. Nowak and R. M. May, āMathematical biology of HIV infections: antigenic variation and diversity threshold,ā Mathematical Biosciences, vol. 106, no. 1, pp. 1ā21, 1991. View at: Publisher Site  Google Scholar
 I. V. Andrianov and L. I. Manevitch, Asymptotology: Ideas, Methods, and Applications, Kluwer Academic Publishers, 2003.
 I. Andrianov and J. Awrejcewicz, āConstruction of periodic solutions to partial differential equations with nonlinear boundary conditions,ā International Journal of Nonlinear Sciences and Numerical Simulation, vol. 1, no. 4, pp. 327ā332, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 C. M. Bender, K. A. Milton, S. S. Pinsky, and J. Simmons, āA new perturbative approach to nonlinear problems,ā Journal of Mathematical Physics, vol. 30, no. 7, pp. 1447ā1455, 1989. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B. Delamotte, āNonperturbative method for solving differential equations and finding limit cycles,ā Physical Review Letters, vol. 70, no. 22, pp. 3361ā3364, 1993. View at: Publisher Site  Google Scholar  MathSciNet
 G. Wu and D. Baleanu, āVariational iteration method for fractional calculus—a universal approach by Laplace transform,ā Advances in Difference Equations, vol. 2013, article 18, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 G. Adomian, āA review of the decomposition method in applied mathematics,ā Journal of Mathematical Analysis and Applications, vol. 135, no. 2, pp. 501ā544, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Atangana and A. Secer, āThe timefractional coupledKortewegdeVries equations,ā Abstract and Applied Analysis, vol. 2013, Article ID 947986, 8 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 G. C. Wu and D. Baleanu, āVariational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers,ā Applied Mathematical Modelling, vol. 37, no. 9, pp. 6183ā6190, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 G. C. Wu, āNew trends in the variational iteration method,ā Communications in Fractional Calculus, vol. 2, pp. 59ā75, 2011. View at: Google Scholar
 A. Atangana and J. F. Botha, āAnalytical solution of the groundwater flow equation obtained via homotopy decomposition method,ā Journal of Earth Science & Climatic Change, vol. 3, p. 115, 2012. View at: Publisher Site  Google Scholar
 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. 56, pp. 597ā603, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 A. Atangana and E. Alabaraoye, āSolving a system of fractional partial differential equations arising in the model of HIV infection of CD4^{+} cells and attractor onedimensional KellerSegel equations,ā Advance in Difference Equation, vol. 2013, article 94, 2013. View at: Publisher Site  Google Scholar
 Ş. 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. View at: Publisher Site  Google Scholar  MathSciNet
 A. Atangana, āNew class of boundary value problems,ā Applied Mathematics & Information Sciences, vol. 1, no. 2, pp. 67ā76, 2012. View at: Google Scholar
Copyright
Copyright © 2014 Abdon Atangana and Emile Franc Doungmo Goufo. 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.