A Method of Directly Defining the inverse Mapping for a HIV infection of CD4+ T-cells model


 In 2015, Shijun Liao introduced a new method of directly defining the inverse mapping (MDDiM) to approximate analytically a nonlinear differential equation. This method, based on the Homotopy Analysis Method (HAM) was proposed to reduce the time it takes in solving a nonlinear equation. Very recently, Dewasurendra, Baxter and Vajravelu (Applied Mathematics and Computation 339 (2018) 758–767) extended the method to a system of two nonlinear differential equations. In this paper, we extend it further to obtain the solution to a system of three nonlinear differential equations describing the HIV infection of CD4+ T-cells. In addition, we analyzed the advantages of MDDiM over HAM, in obtaining the numerical results. From these results, we noticed that the infected CD4+ T-cell density increases with the number of virions N; but decreases with the blanket death rate μI
 .


Introduction
Mathematical models become an important tool in analyzing the dynamics of human immunodeficiency virus (HIV), hepatitis B virus (HBV), and hepatitis C virus (HCV) infection [1]- [2]. The HIV mainly targets a host's CD4+ T-cells, which are the most abundant white blood cells of the immune system. Further, HIV wreaks the most havoc on the CD4+ T-cells by causing their destruction and decline, and decreasing the body ability to fight infection.
Chronic HIV infection causes gradual depletion of the CD4+ T-cell pool. Hence, progressively compromises the host's immune response to opportunistic infections, leading to acquired immunodeficiency syndrome (AIDS). For this reason, the count of CD4+ T-cells is a primary indicator used to measure progression of HIV infection.
In 1989, Perelson et al. [3] developed a simple model for the infection between the human immune system and HIV. Perelson et al. [4] then extended further the model, and observed that the model exhibits many of the symptoms of AIDS seen clinically: The long latency period, low levels of free virus in the body, and the depletion of CD4+ T-cells. They defined the model by considering four separate variables: Uninfected cells, latently infected cells, actively infected cells, and the concentration of free virus particles. Then, described the dynamics of these populations by a system of four coupled nonlinear differential equations. This model was then simplified by Culshaw and Ruan [5] by assuming that all of the infected cells are capable of producing the virus. In other words, they combined the latently infected cells and the actively infected ones.
This simplification resulted in a new system of three differential equations, which we will approximate using the expanded form of MDDiM. The simplified model is as follows: In this model, the variables T(t), I(t), and V(t) are the uninfected CD4+ T-cells, the infected CD4+ T-cells, and the concentration free HIV particles at time t respectively. Here, s is the source of the CD4+ T-cells from precursors, µ T is the natural death rate of CD4+ T-cells, r is their growth rate ( r > µ T in general), and T max is their carry capacity. The parameter k 1 represents the rate of infection of T-cells with free virus, and is the source term for infected cells. k 1 is the rate at which infected cells become actively infected, µ 1 is a blanket death term for infected cells, µ b is the lytic death rate for infected cells, µ V is the loss rate of virus, N is the number of viral particles are released by each lysing cell. Table 1 summarizes parameters and variables.
This particular model and its variations have been solved multiple times by different methods including the HAM by M. Goreishi et al. [6], Homotopy Perturbation Method [7] and the Laplace Adomian decomposition method (LADM) by Ongun [8]. In the present paper we solve the model by MDDiM. This is the first time we used MDDiM to solve a system of three nonlinear coupled ordinary differential equations arise in the field of mathematical biology, using a single inverse linear map to solve the three deformation equations (2.21)-(2.23).

HAM and MDDiM approach
In this section, we first discuss the space that solution and base functions come from, and then will introduce the Method of Directly Defining inverse Mapping Method (MDDiM). MDDiM is an extension to the Optimal Homotopy Analysis Method (OHAM). In OHAM we solve an infinite number of linear ordinary differential equations to obtain series solution but in MDDiM with the help of directly defined inverse map J we solve systems of linear equations. (2.1) Next we define a space for the initial guesses by taking all possible linear combinations of the set S * = {1,t}.
That is, Here we use the first two base functions of the set S ∞ to define S * because each equation in our coupled nonlinear system (1.1) has two initial conditions. Next, we defineV Similarly, let S R = {ψ 0 (t), ψ 1 (t), ...} be an infinite set of base functions that are linearly independent, and define by taking all possible linear combinations of S R .
Now we define three nonlinear operators N 1 , N 2 , N 3 : V → U as (2.10) Here L 1 , L 2 , L 3 are linear operators, T 0 (t), I 0 (t),V 0 (t) ∈ V * are initial guesses, c 0 is the convergence control parameter and q ∈ [0, 1] is the homotopy parameter. The convergence control parameter will be used to optimize the function approximation in the next section.
Assuming that the solutions to the zeroth-order deformation equations are analytic at q = 0, we write Maclaurin series for Φ 1 (t; q), Φ 2 (t; q) and Φ 3 (t; q) with respect to q as where, Here D k is called k th -order homotopy-derivative operator (see [9] for details).
The benefit of OHAM is that we have a considerable freedom to choose linear operators L 1 , L 2 and L 3 so that the corresponding form of equations (2.16)-(2.18) can be determined. Then the number of terms desired in T k , I k and V k of the series solutions can be determine iteratively. OHAM has been successfully used in science and engineering applications (see [10]- [14]), and the choice of axiliary linear operators and convergent control parameters were studied in [11]. The only drawback of OHAM is it takes lot of CPU time to calculate inverse linear operator. To overcome this obstacle, Liao [15] and Dewasurendra et al. [16] introduced the Method of Directly Defining inverse Mapping. Now applying the inverse linear map to OHAM deformation equations (2.16)-(2.18). We obtain the following deformation equations in the frame of MDDiM subject to the initial conditions Here we use only one inverse linear operator to all three deformation equations which leads to less complicated solutions. However, a different inverse linear map could be used to obtain series solutions in different structures. Here we define the inverse linear map J : U → V by where A is a parameter which will be used to optimize the square residual error functions.

Results and Error Analysis
The approximate series solution for the coupled nonlinear system (1.    and define a total error function by taking affine combination of three square residual functions  in table 3 -table 4. From these tables, it is clear that the MDDiM gives an analytical solution with high accuracy only with few iterations when compare to OHAM solution.

Conclusions
In the present study, MDDiM has been developed and used to solve the model of HIV infection of CD4+ T-cells. Approximate analytical solutions for the uninfected CD4+T-cell population size, infected CD4+ T-cell density and initial density of HIV RNA were found. Our analytical solutions are in good agreement with the results of Ghoreishi et al. [4], and with the numerical results, we obtained by Runge-Kutta method (see . From these results, we noticed that the infected CD4+ T-cell density increases with the number of virions N; but decreases with the blanket death rate µ I (see Figures 6,9 and 12).
Since, inverse mapping is directly defined, approximate series solutions are obtained with less CPU time compare with OHAM solution. Also, it is investigated that selected inverse map leads converge series solutions with total error 10 −10 ( see the Table 2 and Figure 4).
In this study, we used a single inverse linear operator for all three deformation equations (2.21)-(2.23), but it is open to use two, or three inverse maps instead. Further, convergence of the series solutions may depend on the choice of inverse linear map and also with the defined solution space. So, it is worth to investigate properties of inverse linear maps in the frame of MDDiM.
To the best of our knowledge, this is the first time someone has used this method to solve HIV infection model. This novel method is more general and can be used to analyze non-linear systems of differential equations arising in science and engineering problems.