On Discretization Method for Optimization Control Model for The Treatment of Pathogenic Induced HIV Infection

Of utmost concern to the scientist, is the lack clear indication for proper documentation of cogent cures for the dreaded human immune deficiency virus (HIV), which usually metamorphoses into acquired immune deficiency syndrome (AIDS). Therefore, continuous research work predominantly on the suppression and prevention of the disease remains an area of paramount interest. Furthermore, the multiplicities of diverse new cases of HIV/AIDS have in no small measure put researchers under an intense pressure. Instances abound where patients are infected with more than a single virus (i.e. HIV with tuberculosis; HIV with hepatitis C; HIV with parasitoid-pathogen, etc.).


Background
Of utmost concern to the scientist, is the lack clear indication for proper documentation of cogent cures for the dreaded human immune deficiency virus (HIV), which usually metamorphoses into acquired immune deficiency syndrome (AIDS). Therefore, continuous research work predominantly on the suppression and prevention of the disease remains an area of paramount interest. Furthermore, the multiplicities of diverse new cases of HIV/AIDS have in no small measure put researchers under an intense pressure. Instances abound where patients are infected with more than a single virus (i.e. HIV with tuberculosis; HIV with hepatitis C; HIV with parasitoid-pathogen, etc.).
Formulation and analysis of the mentioned cases involves a number of parameters. Yet, not much or very little attention has been focus on the estimation of these parameters [1]. Exceptional works on this direction include only on virus clearance rate and the death rate of infected CD4+ T cells [2,3]; the analysis based on the quasi-steady state of the asymptomatic period before it is disturbed by chemotherapy [4]. In [1,5,6], attempts were made on the estimation of parameters of single virus infection-HIV. The model [2,7], discussed the infectious disease at the cellular level to describe the immune system and the host-pathogen interaction.
In this present paper, we consider the multiplicity situation with particular attention on HIV-infection induced by parasitoidpathogen. Here, we formulate a mathematical model which studies the parametric estimation of dual infection (viral load and pathogen), as against several works on single infection [1,2,[5][6][7]. The population density, with which the model is formulated, consists of 5 different subgroups, giving rise to basic 5-dimensional (5D) ordinary differential equations, presented as an optimization control problem. The model is subjected to reverse transcriptase inhibitor (RTI), as a therapeutic factor with the aim of minimizing the systemic cost of treatment.
The objective of optimizing the model is the affordability of predicting infection early stage-called 'set point", a vital breakthrough in treatment decision making [1,5]. Therefore,

Current Trends in Clinical & Medical Imaging
the present work does not only account for the estimation of all parameters in the dynamics of the dual infection, but also account for viral load replication and pathogen resistivity.
The model explores a well-known numerical method-the discretization method called A.V. Kamyad method in solving the formulated 5D-optimal control model with numerical illustrations using Runge-Kutter of order 4 in a Mathcad environment. Basic 3D models have been formulated [3,6].
Other notable models involving discretization method in the optimization control of parameters for single infection includes [7][8][9]. The deployment of highly anti-retroviral therapy (HAART) regimen on treatment and suppression of viral replication and immune system recovery can be found in [6,10,11]. The models [12,13], treated optimal control of HIV infection using Boubakar polynomials expansion scheme and Fuzzy dynamical system respectively. The paper [14], had discussed the impact of numerical stability of the treatment of vertical transmitted HIV/AIDS infections; while Global Convergence and Impact of Multistage and Padé Techniques for Iterative Methods in Nonlinear HIV/AIDS Preventive Chain Model can be found in [15]. The existence and uniqueness of optimal control pair for treatments of an HIV immunology model was studied [16].
The entire model is composed of the introductory aspect as section 1. Definition of the model problem statement and mathematical formulation are absorbed in section 2. Section 3, is devoted to the method and parameter estimation of the model, which explored the AVK discretization methodfor the analysis of the optimal control strategy.Numerical simulations of model illustrations and discussion were designated for section 4. Finally, in section 5, we summarized the work, following the observatory conclusion and remarks base on the formulation and extensive analysis of the model. It is hope that from our present 5D-model, more insight into the optimality of parameter estimation for dual infectious diseases using AVK method of discretizationwill be uncover.

Problem Statement and Model Formulation
In formulating this model, we shall bring to fore, relating ideas of ordinary differential equation in mathematical modeling, structured as problem statement of optimization control, which are solvable using well known method from numerical methods. We construct our model from a considered population density consisting of five different subpopulations giving rise to a set of five ordinary differential equations. Physiologically, let u T denote the uninfected CD4+ T cells, v V -HIV virus (viral load), respectively. The biological descriptionsof the parameters with which these variables interact are defined by [17] and its corresponding ODE's are derived as: Biologically, equation (2.1) assumed a relatively steady viral level during the asymptotic stage of HIV and pathogen infection known as "set-point". At this initial set-point, the body develops an immune system called, the innate immune system, which act as antibodies against HIV-infection and pathogen barrier preventing mechanism. However, the replication of viral load and the rapid adaptivity of pathogen make it impossible for easy detection and neutralization by this innate immune system and gradually leading to full blown AIDS [1,18,19]. Furthermore, it can be shown mathematically from equation (2.1), that the amount of and in the set-point are given by the equilibrium of virus and parasitoid-pathogen depicted i.e.
Thus, model (2.1) adequately reflected the disease progression from initial infection to an asymptotic stage [6].

Methods and Parameter Estimation
In this section, we outline our choice method-the discretization method, which is aimed at estimating all the parameters of HIV and pathogen as involved in our basic model (2.1). This is followed by the optimization control of the uninfected CD4 + T cells and minimizing the cost of treatment factor.

Discretization of the Model
Discretization method affords us the opportunity to transform equation as in model (2.1) into solvable discrete form. Clearly, with discretization method, all the twelve parameters in model (2.1) can be estimated (or measured). As a limitation, we shall deliberately omit the variables ( ) v I and ( ) p I in the initial process due to higher cost quantification of these infected cells [1]. The emphatic nature of this limitation is also stress by [5], on the ground that, the number of infected CD4 + T cells is found to be too small compared to the number of healthy CD4 + T cells.
Then, the higher order differential equations of the output can be derived as: , which also holds for 0 δ ≠ and 2 e τ < [5]. Then, the following inverse map can be derived as: Now, observing equations (2.3)-(2.5), we see that each of them contains four parameters, therefore, it is necessary to further generate at least three other equations on the bases of each of them. This is achieved by differentiating equations (2.3)-(2.5) more times in order to produce the upper derivatives of , 1 2 , y y . 3 y Then, it becomes obvious that at least, four measurement of 1 y , CD4 + T cell count, five measurement of y 2 , viral load and six measurement of y 3, pathogens are needed for a complete determination of equation model (2.1) parameters [5].
As a lid way, we assume that the following measurements are available: Then, by discretization of equations (2.3)-(2.5) and substitution of approximate values of the first derivative of y 1 , the first and second derivatives of y 2 ; and the first to third derivatives of y 3 for CD4 + T cells, viral load and pathogen respectively, we can establish the measurement of the variables at different time intervals i.e.
So, we could as well obtain similar matrices for equations (2.9)-(2.10). Thus, the variables i l , i = 1, 2… 12 and all the parameters from equation (2.6) of the basic model, can be calculated. A case study of some of the parameter estimation as in basic model (2.1) can be found in [5]. Now, applying the bases of equation (2.7), we establish (Table 1) and 0( ) 0.2 p P = is determined using well-known numerical method like RK4 in a Mathcad environment.The graphs simulation of the susceptible CD4 + T cells, virus -infected CD4 + T cells, viral lead, pathogenic-infected CD4 + T cells and the parasitoid-pathogen, will be presented in section four.

Optimal Control of HIV and Pathogen Infection
The increasing therapeutic drugs for the treatment of AIDS which also serve as treatment of most pathogens by the activation of required immune system (adaptive immune system) can be conveniently grouped into three major groups-Reverse transcriptase inhibitor (RTI), protease and integrase enzyme inhibitor. In this present model, we conduct our studies using the (RTI). This choice for RTI is not unconnected with the dual characteristics of RTI on both viral load and the activation Current Trends in Clinical & Medical Imaging of adaptive immune system against parasitoid-pathogen. RTI is responsible for the prevention of uninfected lymphocyte cells from infection by viral load and as well elimination of infected pathogen cell [16]. Corroborating the above assertion is the commensurate reaction of these two infectors with their host as reference in Table 1: Therefore, the utmost concern of the present work is to control model (2.1), in order to prevent any transmission of uninfected cells to virus and pathogenic infected cells. This control function is given as ( ) w t , where 0 ( ) 1 w t ≤ ≤ .The most efficiency ofdrug in this case, is when , implying that CD4 + T cells are not infected by viral load and pathogen anymore. At the tail side, if , ( ) 0 w t ≡ then drug is said to be inefficacious and therefore, disease progression does not change. On the basis of this constrains, we derive from model (2.1), the control system as follows: We now define the objective functional of equation (2.12), as: using [19], with the weight factor 10 ψ = . Here, we aimed at maximizing the objective functional (2.13) subject to the control system (2.12); that is, maximizing the total count of CD4 + T cells and minimizing the cost of treatment via RTI drugs. We resolve this optimal control problem by numerical methods [7], using a special discretization method called A.V. Kamyad [19,21]. Applying the AVK method, we let, Subject to: , for i = 0,1,……,r-1 and .
Step 5. Taking  We then apply this technique to calculate the control problem of equation (2.12), with the objective functional (2.13). The parameters used in the basic control model (2.12) are in agreement with that in [5]. So, we can assume that the treatment starts when CD4 + T cells was at minimum count as pathogen is on the decrease. Results are then confirmed by extrapolation from our treatment intervals. Therefore, from Figure 1, for , is said to be minimum count of CD4 + T cells. This shows that treatment interval is [3,30] years. Furthermore, we finally partition [3,30] years into r parts with length g. Thus, the discretization of equation (2.18) becomes: We present numerical illustration and results of this optimization problem which is obtained using RK4 in a Mathcad environment in the coming section.

Numerical Results and Discussion
Using an optimal weight factor 10 ψ = , numerical solutions for the uninfected CD4 + T cells, viral load, and pathogen concentration, as well as the optimal treatment strategy have been generated. From Figure 2 below, the solution for [3,30] t ∈ , having initial variables  Table 2 below:is used for the determination of treatment "set point" and the concentration of various population density understudies. It is seen from the graphical representation that at (3) 0.8 T = , is the minimum count of CD4 + T cells. So, the treatment interval is [3,30] years.
Noted that this interval includes: .   Table 3(a) below: represents the calculation for the maximization of the uninfected CD4+ T cells while minimizing the systemic cost of therapeutic treatment (RTI). We observed that from Figure 2 The concentration of the uninfected CD4+ T cells is seen increasing rapidly, which affirmed the overall healthy population density. From Table 3(b) below, as illustrated by Figure 2(b) above, we see that as optimal therapeutic treatment is consistently given at a steady systemic cost, the uninfected CD4+ T cells increases in a logistic way (manner), explaining why in Figure 2(a), the viral load and pathogen decreases in concentration in the blood plasma. Figure 2(b) shows that active therapeutic treatment terminates after 18 years, which again ascertained the necessity for our initial pre-condition for a definite time interval for a realistic therapeutic investigation. Mathematically, from the state system (2.12)-(2.13), the choice for the time interval [3,30] t ∈ conveniently accommodated the clinical latency stage, which has a maximum period of 10 years [22]. Moreover, the weight factor ψ , determined the level at which therapeutic treatment is administered within the period of time. The implication is that if systemic cost of treatment increases, the patient receives the maximal therapy for a finite period of ψ factor.

Conclusion
In this paper, numerical optimization of nonlinear 5-Dimensional mathematical HIV dynamic model for the parametric estimation of pathogenic induced HIV infections Current Trends in Clinical & Medical Imaging was formulated. The dynamics of the infection was presented as optimization control problem, solvable by a special discretization method called AVK method. Numerical computations showed that, with defined control function, a realistic optimal control minimal therapeutic cost can be achieved within measured time interval of 18 years, while assiduously maximizing healthy CD4 + T cells concentration as a measure of the parameter estimation. The study therefore, suggests the use of the model for the optimization control problems of other related infectious diseases.