# Nonlinear and Robust Control Strategy Based on Chemotherapy to Minimize the HIV Concentration in Blood Plasma

Authors: Ricardo Aguilar-López; Rigel Valentín Gómez-Acata; Gerardo Lara-Cisneros; Ricardo Femat
Journal: Journal of Control Science and Engineering
 (2016)

---

## Abstract

A nonlinear PI-type control strategy is designed in order to minimize the HIV concentration in blood plasma, via medical drug injection, under the framework of bounded uncertain input disturbances. For control design it is considered a simplified mathematical model of the virus infection as a benchmark. The model is based on mass balances of healthy cells, infected cells, and the virus concentrations. The proposed controller contains a nonlinear feedback PI structure of bounded functions of the regulation error. The closed-loop stability of the system is analyzed via Lyapunov technique, in which robustness against system disturbances is demonstrated. Numerical experiments show a satisfactory performance of the proposed methodology as a HIV therapy, in which the virion particles and the infected CD4+T cells are minimized and, as an interesting result, the drug dosage can be suspended, thus avoiding drug resistance from the virus. Finally, the proposed controller is compared to a standard sliding-mode and hyperbolic tangent controllers showing better performance.


---

## Body

## 1. Introduction

The Human Immunodeficiency Virus (HIV) causes the named Acquired Immune Deficiency Syndrome (AIDS). This virus is distributed through the mucous membranes and is often transmitted through bodily fluids. It destroys the immune system, by infecting CD4+T cell and the macrophages. Consequently, a poor immune system makes it difficult for the organism to fight against opportunistic infections as pneumonia, meningitis, or tuberculosis causing death. A healthy person has a CD4+T count ranging from 500 to 1000 cells mm−3. A person with a count of 200 cells mm−3 or less of CD4+T is considered to have AIDS. There are rare patients who are infected with HIV, but they are able to suppress the virus without drugs therapy and hence do not develop AIDS; this status is known as long-term nonprogressor (LTNP) [1, 2].The HIV replication cycle can be summarized in six steps: binding and entry; no coating; reverse transcription; provirus integration; virus protein synthesis and assembly; and finally budding [3, 4]. A cure for the HIV infection still has not been found; however since the synthesis of zidovudine (AZT), several drugs have been synthetized belonging to 6 different inhibitor classes approved to fight against the virus and its complications, reducing the viral load. Unfortunately, the synthesis of new drugs is still needed due to the virus drug resistance ability [5]. The drug’s target is to prevent the HIV infection in the organism and delay the AIDS, by acting on enzymes and proteins that are vital for the virus to complete its replication cycle and as a result reaching a status in which the patient has a sufficient number of CD4+T cells to fight off HIV and other infections, keeping the immune system as healthy as possible. From experience, these drugs alone or in combination are not capable of dropping the viral load to be undetectable for sustained periods of time [5].One way to analyze the evolution of the HIV dynamic is the mathematic modeling considering the underlying infection mechanisms of the HIV/AIDS to understand and anticipate the spread and evaluate the potential effectiveness of different therapies to bring the infection under control [6]. In keeping with the above, several types of basic modeling approaches have been developed as new information about the HIV pathology and immune system, and therapeutic effects of drugs arise such as the deterministic models, the statistical models, and the stochastic models. In general, they can predict the early phases of disease but cannot predict the decrease in CD4+T cell counts and rise in peripheral virus seen in AIDS [7, 8].Additionally, the HIV mathematical modeling has allowed the analysis of the dynamic behavior of the virus infection, which has helped to propose regulation strategies to diminish indirectly the viral concentration [9–11]. An increasing number of works in the literature deal with the HIV virus regulation via control approaches and reduced HIV models [12–20]. Kim et al. [21] found through a bifurcation analysis the range of optimal constant drug dosing for the reverse transcription inhibitor and the protease inhibitor in order to achieve LTNP in HIV patients. Moreover, by controllability analysis and its singular values, they showed that a therapy with lower/moderate efficacy had a relatively higher efficiency in HIV treatment in terms of controlling the immune system. Barão and Lemos [22] studied the nonlinear control of HIV-1 infection for a three-state model, where the number of healthy and infected CD4+T cells and the number of virion particles are considered; the control law employed a reduced HIV model, comprising two nested loops, a linearized feedback and a LQ regulator. It showed that it is not possible to completely eliminate the infection, but it can drive the virion particles to a number below the threshold specification (50 virus copies/mm3 in plasma are undetectable in HIV patients). Also they showed that the controlled system is stable to moderate uncertainties in the HIV model parameters.In this work a novel nonlinear PI-type controller is proposed in order to provide systematic drug-scheduling which drives the HIV virus concentration to minimal values in simulation tests. Unlike the previous therapies based on control algorithms, the proposed control strategy has the advantage of interrupting the drug dosing permanently. The proposed control law contains a class of continuous bounded functions, in order to compensate the dynamic behavior of the virion concentration in plasma and possible unknown disturbances onto the system. Numerical experiments show the successful performance of the proposed methodology as an idealized HIV therapy.


## 2. Mathematical Model for HIV Infection

The mathematical model to describe the dynamics of the HIV and CD4+T cells during the primary infection stage was taken from [23]. The model is based on the mass balance in the blood plasma taking into account the concentrations for uninfected CD4+T cells x
1, infected cells (
x
2
), and virion particles in plasma (
x
3
), all of them expressed in mg L−1 as follows:(1)
x
˙
1
=
k
1
-
k
2
x
1
-
k
3
x
1
x
3
,
x
˙
2
=
k
3
x
1
x
3
-
k
4
x
2
,
x
˙
3
=
k
5
x
2
-
k
2
x
3
+
u
twith x
=
x
1
x
2
x
3 as the state vector which belongs to a compact domain Ω
∈
R
+
3. The term k
3
x
1
x
3 corresponds to the rate of the infection of the CD4+T cells by the virus and k
i
x
2 is the rate of diminishment of infected cells. The parametric identification of this kind of systems is somewhat difficult, due to the quality of the available experimental data [24]; in particular, the measurements of the parameters and variables of the HIV medical tests depend, in some degree, on each commercial laboratory, such that each one has its own set of normal values and these values depend on the equipment or method used [25]. Generally the parameters for these models are reported in cell/microliter; however the drug administration is related with mass concentration units (mg L−1) [26]; therefore in order to be in agreement with the mass balances units, the following parameters are considered for this work: k
1
=
50.0 mg L−1 d−1, k
2
=
0.05 d−1, k
3
=
5
e
-
4 L mg−1 d−1, k
4
=
0.4 d−1, k
5
=
40.0 d−1, and k
6
=
9.0 d−1; all of these kinetic constants are nonnegative and k
1 represents the constant rate of CD4+T cell generation; k
2 is their death rate per capita. The term u
(
t
) is the exogenous control input given by the drug injection administration, under the assumption that the corresponding drug is able to eliminate the virion particles in plasma. From the above, the control input is imposed additively in the corresponding state equation.


## 3. Control Design

Let us consider a generalized space state representation of the system (1):(2)
x
˙
=
f
x
+
φ
t
+
g
x
u
t
,
(3)
y
=
h
x
=
C
x
,where x
∈
R
+
n is the corresponding state vector, y
∈
R
+
m is the measured system output vector, with m
≤
n,  f
:
R
n
→
R
n is a smooth vector field given by f
(
x
)
=
k
1
-
k
2
x
1
-
k
3
x
1
x
3
;
k
3
x
1
x
3
-
k
4
x
1
;
k
5
x
1
-
k
6
x
3
T, g
(
x
)
=
0
0
1
T, C
=
0
0
1
T is the measurement coefficient matrix, n
=
3 is the order of the system (2), m
=
1 is the dimension of the measured vector for the specific case of study, and φ
(
t
) is an additive bounded disturbance.Now, let us consider the following hypothesis:(H1)
The control inputu
t is bounded; that is, u
t
≤
u
-.
(H2)
φ
(
t
)
≤
ð, f
(
x
)
∈
C
∞.
(H3)
f
(
x
)
≤
F, g
(
x
)
≤
G, ∀
x
∈
R
+
3
,
F
,
G
,
ð
<
∞.Proposition 1.
Consider the dynamical system (1) on the compact domain Ω
∈
R
+
3 and suppose that H1–H3 are satisfied, then the following feedback control input achieves the regularization of (1) at the desired reference value x
s
p:(4)
u
t
=
-
g
1
e
t
1
+
e
t
-
g
2
sign
⁡
e
t
∫
0
t
e
x
p
-
e
σ
2
d
e
σ
,where e
(
t
)
=
C
x
-
C
x
s
p is the regulation error, x
s
p is the desired set point and g
1 and g
2 are the named control’s gains.The proposed controller is inspired on the framework of the sliding-mode theory, which provides some robustness properties against external disturbances [27]; however and despite its robust performance, the employment of discontinuous controllers leads to an important closed-loop performance degradation due to the chattering phenomena. From the above, some alternative designs as the high order sliding-mode and sigmoid controllers have been designed in order to provide enough smoothness to the control action [28, 29].Therefore, the proposed controller includes two bounded sigmoid functions as described by (4); the first term of the proposed methodology tries to compensate the nonlinear behavior of the system and the second one tries to compensate the additive disturbance. Note that the integral term is multiplied with the sign function in order to allow a sign change in this term and diminish the chattering from the standard sliding-mode controllers, smoothing the discontinuous nature of the sign function.


### 3.1. Sketch of Proof of Proposition1

Consider the Lyapunov candidate function:(5)
℧
=
1
2
e
2
.The time derivative along the trajectories of (5) is(6)
℧
˙
=
e
e
˙
.Here is considered a regulation case, where the set point is constant; that is, x
˙
s
p
=
0, in consequence e
˙
=
C
x
˙ and φ
(
t
).Therefore,(7)
℧
˙
=
e
C
f
x
+
φ
+
g
x
u
t
.Substituting the proposed feedback control input in (7),(8)
℧
˙
=
e
C
f
x
+
φ
+
g
x
-
g
1
e
1
+
e
-
g
2
sign
⁡
e
⁡
∫
0
t
e
x
p
-
e
2
d
e
σ
.Now, considering that(9)
e
1
+
e
≤
1
,
s
i
g
n
e
∫
0
t
e
x
p
-
e
2
d
e
σ
≤
1
.The above, because the function e
/
1
+
e is a class of sigmoid functions, s
i
g
n
e is the standard sign function and ∫
0
t
e
x
p
-
e
2
d
e
(
σ
) is an erf-type (sigmoid) function.Now, considering the hypotheses (H1) to (H3),(10)
℧
˙
≤
e
C
F
+
ð
+
G
-
g
1
e
1
+
e
-
g
2
s
i
g
n
e
∫
0
t
e
x
p
-
e
2
d
e
σ
.Finally,(11)
℧
˙
≤
C
F
-
g
1
G
+
ð
-
g
2
G
e
.If(12)
g
1
G
>
F
,
g
1
>
G
-
1
F
,
g
2
G
>
ð
,
g
2
>
G
-
1
ð
,therefore,(13)
℧
˙
≤
0
.Therefore, it can be concluded that the closed-loop system (2) is stable at e
(
t
)
=
0, such that x
=
x
s
p.In particular, the integral term is considered to compensate the effect of the additive bounded disturbanceφ
(
t
), as is showed in the stability analysis under the Lyapunov framework (see (11)).


## 3.1. Sketch of Proof of Proposition1

Consider the Lyapunov candidate function:(5)
℧
=
1
2
e
2
.The time derivative along the trajectories of (5) is(6)
℧
˙
=
e
e
˙
.Here is considered a regulation case, where the set point is constant; that is, x
˙
s
p
=
0, in consequence e
˙
=
C
x
˙ and φ
(
t
).Therefore,(7)
℧
˙
=
e
C
f
x
+
φ
+
g
x
u
t
.Substituting the proposed feedback control input in (7),(8)
℧
˙
=
e
C
f
x
+
φ
+
g
x
-
g
1
e
1
+
e
-
g
2
sign
⁡
e
⁡
∫
0
t
e
x
p
-
e
2
d
e
σ
.Now, considering that(9)
e
1
+
e
≤
1
,
s
i
g
n
e
∫
0
t
e
x
p
-
e
2
d
e
σ
≤
1
.The above, because the function e
/
1
+
e is a class of sigmoid functions, s
i
g
n
e is the standard sign function and ∫
0
t
e
x
p
-
e
2
d
e
(
σ
) is an erf-type (sigmoid) function.Now, considering the hypotheses (H1) to (H3),(10)
℧
˙
≤
e
C
F
+
ð
+
G
-
g
1
e
1
+
e
-
g
2
s
i
g
n
e
∫
0
t
e
x
p
-
e
2
d
e
σ
.Finally,(11)
℧
˙
≤
C
F
-
g
1
G
+
ð
-
g
2
G
e
.If(12)
g
1
G
>
F
,
g
1
>
G
-
1
F
,
g
2
G
>
ð
,
g
2
>
G
-
1
ð
,therefore,(13)
℧
˙
≤
0
.Therefore, it can be concluded that the closed-loop system (2) is stable at e
(
t
)
=
0, such that x
=
x
s
p.In particular, the integral term is considered to compensate the effect of the additive bounded disturbanceφ
(
t
), as is showed in the stability analysis under the Lyapunov framework (see (11)).


## 4. Results and Discussion

Numerical experiments were done in order to show the open-loop and closed-loop behavior of the model for HIV infection. From the above, a PC computer with Intel Core i7 processor was used and the ordinary differential equations systems solver libraries, in particular, 23s ODE solver from MATLAB™, were employed. The system was simulated considering the following initial conditions: x
1,0
=
350 mg/L, x
2,0
=
125 mg/L, and x
3,0
=
75 mg/L; the disturbance considered in the virion particles in plasma is given by φ
(
t
)
=
5
s
i
n
(
10
t
). Firstly, Figure 1 shows the behavior of uncontrolled system, where an important increase for the HIV concentration compared to the initial condition can be observed, while the final mass concentration of the uninfected cells diminished to 180.5 mg/L and the final infected cell concentration is 102 mg/L; these values provide an important evidence that the infection has not yet reached the AIDS stage. Figure 2 shows the corresponding 3D phase portrait, where the evolution of the trajectory described above is observed. The chemotherapy treatment via drug injection started at t
=
75 days; the treatment was conducted by the proposed controller; the parameters of the proposed controller are g
1
=
1100 d−1 and g
2
=
10 d−2. In Figure 3 an immediate response of the HIV concentration reaching the required set point (x
s
p
=
0 mg/L) in 15 days from the open-loop HIV concentration of 450 mg/L is observed. Besides, Figure 4 is related with the immediate response of the increase of the uninfected cell concentration (from 180.5 mg/L to 977 mg/L) and the decrease of the infected cells (from 102.5 mg/L to 0 mg/L); under this frame the corresponding closed-loop trajectory is presented in a 3D phase portrait (Figure 5). An important goal of the proposed methodology was the named control’s effort; when the chemotherapy was started, the proposed methodology injected 1100 ng/L of the hypothetic drug; on the other hand, when the HIV concentration falls, the proposed methodology decreased the concentration of the injected drug to reach a null concentration. Therefore it is not necessary to maintain the drug application in order to keep the HIV concentration in a small enough value; this can be observed in Figure 6. For comparison purposes, a standard sliding-mode and hyperbolic tangent feedback [30–32] are implemented; they are compared with the proposed methodology via performance index proposed in [33]. The performance index is related with the integral evolution of the average of the square of the regulation error; the methodology with the lower index is considered to be of the best performance. As observed in Figure 7 the methodology proposed in this work presents the better performance.Figure 1
Open-loop behavior of the HIV dynamics.Figure 2
Phase portrait of the open-loop system.Figure 3
Closed-loop dynamics of the HIV concentration.Figure 4
Closed-loop behavior of the uncontrolled concentrations.Figure 5
Closed-loop phase portrait of HIV dynamics.Figure 6
Control’s effort for the proposed controller.Figure 7
Performance index for the considered controllers.


## 5. Conclusions

In this work an alternative therapy based on a nonlinear PI feedback control law to minimize the HIV concentration in blood plasma was presented. The proposed methodology contains a class of bounded functions in order to provide a satisfactory performance in the regulation of the number of virion particles and infected CD4+T cells using a single input model with unknown disturbances. This approach shows by simulation an effective diminishment of the HIV concentration under the drug dosing scheme, which drives the cells near to health, improving the quality of life and longevity. Also this approach can avoid the drug resistance due to its minimal drug dosing values. Moreover, a better performance of the current methodology was achieved in contrast to other control laws.



---

*Source: 6284574-2016-05-19.xml*