1. Introduction
The genetic code of the HIV virus consists of a ribonucleic acid, the RNA. HIV belongs to the retroviruses family, characterized by the presence of an enzyme, DNA-polymerase RNA-independent, capable of transcribing the genetic code RNA into DNA. This ability allows the virus to integrate its genome into the one of the cells it infects, so that the integrated virus would not be defeated nor by the immune response nor by drugs.
HIV has a marked tendency to mutate: mutations are errors that individual viral particles make in replicative cycles. Each error leads to the appearance of a virus, which is more or less analogous to the original one. Mutations are mostly “disadvantageous” for the viral species, as mutated viruses tend to disappear. Nevertheless some mutations are “beneficial” and allow mutated viruses to acquire drug resistance and immune response. One of the most effective mechanisms used by HIV to evade the immune response and treatment is, in fact, its aptitude to change.
The main HIV target cells are T lymphocytes of type CD4, fundamental in the adaptive response against a variety of pathogens and oncogenes.
The overall function of the immune system is to prevent or limit infections. An immune response is generally divided into innate and adaptive immunity. Innate immunity occurs immediately, when circulating innate cells recognize a problem. Adaptive immunity occurs later, as it relies on the coordination and expansion of specific adaptive immune cells.
In particular, CD4+ lymphocytes are very important cells in the immune system, in fact they recognize the various uninvited guests organism (viruses, bacteria, protozoa, fungi, worms and cancer cells) through biochemical messages, and they activate the areas of the immune system most suitable to counter their presence. A large number of CD4+ paralyzes the immune system, exposing the body to any risk of infection and cancer [
1].
Actually, HIV can be suppressed by a combination of antiretroviral therapy (ART) consisting of three or more antiretroviral (ARV) drugs. In principle, ART cannot eradicate HIV infection from the “sanctuaries” (mainly lymph-nodes and lymphatic system); it controls viral replication within a body reducing viral burden, thus in turn it allows a functional improvement of an individual’s immune system that basically regains the capacity to fight off infections [
2]. It has been proved that people with HIV, subject to ART therapy, can have a healthy and productive life, in fact after the advent of the ART therapy the mortality curve of HIV-infected patients started to decline.
On the other hand, the lack of a definitive eradication of viral reservoirs determines two important consequences: the anti-HIV treatment should be life-long; the longer the treatment, the more likely the development of drug-resistance. On this basis, an intense research activity has been developed in last decades, usefully enlarging the anti-HIV molecules repertoire. Actually there are more than 20 approved antiretroviral drugs, divided into 6 different pharmacological classes. Anti-retroviral drugs are broadly classified through the phase of the retrovirus life-cycle that they inhibit, as follows
Nucleoside/Nucleotide Reverse Transcriptase Inhibitors (NRTI);
Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTI);
Protease inhibitors (PI);
Fusion Inhibitors (FI);
Coreceptor Antagonists (CA);
Integrase Inhibitors (II).
Typical combinations involve the conjunctive use of either two NRTIs and one NNRTI, or one PI (with or without Ritonavir) in combination with two NRTIs.
HAART (Highly active antiretroviral therapy) is an abbreviation for all protocols involving combinations of drugs, which are active against different molecular targets in the life cycle of HIV. These medications are administered in the form of the high concentration cocktails. This approach was born in 1995–1996 with the introduction of the second class of antiretroviral drugs: Protease inhibitor (PI), administered in combination with drugs of the first generation, or Reverse Transcriptase Inhibitors (NRTI and NNRTI: see zidovudine or AZT azidothymidine).
Among the benefits of treatment there is a decreased risk of progression to AIDS and a decreased risk of death. The adverse effects and the complexity of treatment regimens (due to the high number of pills and to dosing frequency) may reduce patient’s compliance [
3].
A correct administration of antiretroviral drugs should be daily scheduled and uninterrupted; it requires care and precision, even in relation to meals and to subministration of other drugs. The omission of a few doses leads to a reduction of drug concentration in blood, therefore the residual level could become lower than the one necessary to inhibit the virus. Not only does this allow the resumption of viral replication, but it also facilitates the inexorable emergence of resistant virus. Resistant mutants tend to persist, making the drug ineffective, even if the intake of drugs is carried out according to the rhythms and the correct doses. A scarce adherence to the prescriptions and its consequent emergence of resistant mutants may cause failure of the treatment.
The therapy must prevent viral replication in the body, rather than the complete eradication of the infection, that remains chronic, so as to reduce the damage caused to the immune system and to allow greater survival and quality of life.
The study of HIV biological evolution and of its relationship with the immune system is used to determine a therapy policy that can defeat the viruses. Despite it seems counterintuitive, the aim is to determine a proper dosage of drugs, which defends against viruses only moderately and which may benefit both the host and the virus, i.e., without creating strong environmental pressures. A stable long-term coexistence can be reached, hopefully avoiding drug-resistance development.
Mathematical modeling of HIV infection has proven to be instrumental for the modern understanding basis of the AIDS pathogenesis [
4]. There exist several attempts to formalize the evolution of HIV and the use of drugs to limit its diffusion. Most of the related literature is linked to the study of the HIV dynamics in the body and only few of them apply the Optimal control theory considering also the possibility to control the drugs action. An exhaustive survey that collects all these works and classifies them according to the most relevant criteria con be found in [
5].
In this work we present an application of the differential Game Theory to a medical-therapeutic context for the HIV treatment A first attempt has been performed by Wu and Zhang in [
6]. Here we present a differential game which considers the classes of antiretroviral drugs currently most used and different immune cells dynamics, with the aim of representing as much as possible the real setting of this problem. We consider two players: The HIV virus and the immune system, supported by antiretroviral therapy. We look for an optimal therapeutic treatment in terms of a Nash Equilibrium, with the aim of finding a proper dosage of drugs which defends against viruses, such that a stable long-term life expectation may be obtained.
This paper is structured as follows. In
Section 2 we introduce the involved variables and parameters and we formalize the problem in terms of a differential game. In
Section 3 we characterizes the Open Loop Nash Equilibrium. In
Section 4 we present some numerical simulations to show two different situations that may occur and how the algorithm performs.
Section 5. concludes with some comments. All analytical computations are reported in
Appendix A within the Theorems’ proofs. In
Appendix B are listed the parameter values used in the numerical simulations.
2. The Model
We present a generalized model within the Game Theory approach, in order to determine the optimal antiretroviral treatment against the HIV infection. We take into accounts 4 types of cells: HIV viruses, CD4+ T Helper cells (adaptive immunity), macrophage and monocytes cells (innate immunity) and immune precursor/effector Cytotoxic T lymphocyte cells.
After the primary HIV infection has established, an acute HIV syndrome appears with a wide dissemination of viruses and with seeding of lymphoid organs. This phase can last between three and nine weeks. From the ninth week the clinical latency phase may start, here the AIDS symptoms are not macroscopically manifested, nevertheless the number of T lymphocytes begins to decrease and the viral load increases. This stage can persist many years. Finally, at a last stage, there is the occurrence of the AIDS symptoms and of other opportunistic infections. In addition, this stage can last for years and ends with the death of the host [
7]. As there exist so many strains of viruses, in order to take into account all of them, without making the model intractable, we distinguish the HIV viruses between “sensitive”, (
) and “resistant” (
) to the therapy.
In the model we adopt the idea of Herod
et al. [
8] and Nowak [
9], assuming that the immune response (
) to a viral infection creates some subpopulations of lymphocytes (
) that are specific only for sensitive viruses, some others (
) that are specific only for resistant viruses, and finally the subpopulation
that is effective against both sensitive and resistant viruses. Furthermore, the mutation of the initial viral infection may cause the death of the entire population of lymphocytes.
Given the fact that HIV actually cannot be eradicated, our aim is to formulate a differential game between HIV and the Immune System in which the existence of an equilibrium would lead to an optimal drug therapy, which enables a stable long-term coexistence between virus and host, in other terms, a longer life expectancy.
In the following tables we present, in alphabetical order for the symbols, the variables and the parameters of the model. In particular:
Table 1 lists the state functions;
Table 2 lists the control functions and finally
Table 3 and
Table 4 list the parameters.
Table 1.
State functions (Unit of measure: ).
Table 1.
State functions (Unit of measure: ).
Symbol | Description |
---|
| Immune Precursor Cytotoxic T lymphocytes |
| Immune Effector Cytotoxic T lymphocytes |
| Uninfected macrophages |
| Macrophage infected by viruses j |
| Uninfected CD4+ T Helper cells |
| Uninfected CD4+ T Helper cells strain specific for |
| Uninfected CD4+ T Helper cells strain specific for |
| Uninfected CD4+ T Helper cells unspecific for and |
| Latently infected CD4+ cells infected by virus j |
| Actively infected CD4+ cells infected by virus j |
| HIV drugs sensitive viruses |
| HIV drugs resistant viruses |
Table 2.
Control functions (bounded in ).
Table 2.
Control functions (bounded in ).
Symbol | Description |
---|
| Dosage of coreceptor antagonists |
| Dosage of fusion inhibitors |
| Immune boosting |
| Dosage of integrase inhibitors |
| Dosage of protease inhibitors |
| Dosage of reverse transcriptase inhibitors |
| Mutation rate from resistant virus () to other resistant virus () |
| Mutation rate from resistant virus () to sensitive virus () |
| Mutation rate from sensitive virus () to resistant virus () |
| Mutation rate from sensitive virus () to other sensitive virus () |
Table 3.
Parameters (Unit of measure: ).
Table 3.
Parameters (Unit of measure: ).
Symbol | Description |
---|
c | Cytotoxic T lymphocyte (CTL) activation rate |
| Death rate of cells infected by viruses due to immune response |
k | Rate at which T cells convert to specific immune reaction cells ( or ) |
| Rate at which T cells convert to unspecific immune reaction cells |
| Rate at which latently infected cells convert to actively infected cells |
| Natural death rate of type i cells |
| Growth rate of sensitive viruses () |
| Growth rate of resistant viruses () |
q | Growth rate of due to infected cells and |
r | Growth rate for CD4+ T cell population |
Table 4.
Parameters.
Symbol | Description | U.m./Value |
---|
| Weight on the benefit i and cost of therapy i | |
| Weight of viral mutation from strain i to strain j | |
| Half saturation constant | |
g | Input rate of external viral source | |
| Initial value of cells | |
| Initial value of cells | |
| Rate at which viruses type j infect T cells | |
| Rate at which viruses type j infect M cells | |
| Rate at which infected macrophages infect T cells | |
| Initial value of cells | |
| Initial value of actively cells infected by viruses | |
| Scaling parameter for type i cells |
| Specific immune response rate against viruses | |
| Percentage of i cells that recognize the virus j | |
| Average number of virions j infecting a cell | |
| Source/production of type j cells | |
| Unspecific immune response rate against viruses | |
| Initial value of cells | |
| Initial value of latently cells infected by viruses | |
| Initial value of actively cells infected by viruses | |
| Maximum CD4+ T Helper cell population level | |
| Initial value of viruses | |
Let us collect all the state variables in the following array:
The last four controls are related to the HIV player, nevertheless, in line with most notation of the cited literature, we can limit to the two controls and only, observing that and
Let’s collect the immune system-therapy controls in the array:
and the viruses controls in the array:
In
Table 3 we show the parameters with the same unit of measure (
), while all the other parameters are in
Table 4.
HAART therapy can prolong the patient’s life, however, it has many side effects, so that a long-term administration may become difficult and it is necessary to determine the minimal dose of drugs that prevents the viral replication.
The immune system has to maximize the number of uninfected macrophages, uninfected T-cells, and immuno precursors/effectors, minimizing, at the same time, the side effects of drugs. Such an object is represented by the following payoff:
where
and
.
On the other hand, HIV tends to maximize the number of sensitive and resistant viruses, the number of latently and actively infected T-cells, minimizing, at the same time, the mutation costs. Its payoff fuction is:
where
and
. Observe that drug toxicity costs and HIV mutation costs are assumed quadratic in the controls to represent their increasing increments of scale.
The T-cells undergo to different evolutions depending on whether they have been infected or not by the virus. In what follows we distinguish between uninfected and infected cells. In particular, we expose the dynamics of the following five classes of cells:
Uninfected CD4+ T-helper cells;
Uninfected macrophages and CD8+ CTL precursors and effectors;
Latently and actively infected CD4+ T cells by sensitive and resistant viruses;
Infected macrophages by sensitive and resistant viruses;
Sensitive and resistant viruses;
Writing first the terms which contribute positively to their growth, and after the negative terms. We also stress the features that characterize each group of cells.
The evolution of every type of cell is subject to a natural death rate that can be represented in its dynamics by a (negative) decaying term. Namely, let be a particular type of cell, such a decaying term is and appears at the end of each considered dynamics.
In the evolution of the
uninfected T-Helper cells we consider features taken from different models by Kirschner
et al. [
10], Caetano
et al. [
11] and Joshi [
12]. The dynamics are the following
One common characteristic among the above cells is constituted by the presence of the immune busting effect that causes an increase of the immune barriers, see [
12]. It is proportional to the number of particular uninfected
T-Helper cells and it depends on the control
Note the positive term
that appears at the end of the first line in the dynamics of each type of cell
Moreover, some cells are infected by sensitive viruses: Their number is proportional to the number of sensitive viruses () and to the dimension of the community of cells itself. Coefficient represents the infection rate, the factors and describe the actions of fusion inhibitors and coreceptor antagonists against the infection of the cells. These drugs interfere with binding, fusion and entry of HIV to the host cell. Finally is the portion of generic CD4+ T Helper cells that does not recognize and does not obstruct sensitive viruses. Similarly, cells are infected also by resistant viruses, and their number is proportional to the number of resistant viruses () and to the dimension of the community of cells itself. What makes the difference is that fusion inhibitors and coreceptor antagonists can not counteract cells infection by resistant viruses (), therefore the drugs action is null. The factor is the portion of generic CD4+ T Helper cells that does not recognize and does not obstruct resistant viruses. These effects are represented by the second lines of the dynamics.
Observe that the negative terms
which appear in the third line of Equation (
3) represent the number of
cells converted to
,
or
cells respectively. Each one of such terms obviously appears with a positive sign in Equations (4)–(6) for uninfected
and
cells.
The term gives the portion of cells contaminated by infected macrophages, and the sensitivity/resistant differentiation in it constitutes a novelty of our model.
Once explained the common features of the
T-Helper cells, let’s explicate the terms which characterise the dynamics of the naive T lymphocytes (
) in Equation (
3).
The term:
gives the proliferation of uninfected CD4+
cells. It includes both an external (not plasma) contribution of cells from sources, such as the thymus and lymph nodes, and an internal (plasma) contribution from CD4+
cells differentiation.
The term:
represents the production of
-cells due to cloning in the presence of an antigen, taking into account the maximum number of lymphocytes,
.
The dynamics for the
uninfected macrophages, CTL precursors and effectors are the following:
According with [
10], macrophages dynamics Equation (
7) consider a constant proliferation source (
) and sensitive and resistant infections elements
Here reverse transcriptase inhibitors (RTI) and protease inhibitors (PI) counteract the action of sensitive viruses. In particular the first ones inhibit reverse transcription, and the second ones block the viral protease enzyme necessary to produce mature virions upon budding from the host membrane. Each one of such terms obviously appears with a positive sign in Equations (
26) and (27).
About immune precursors/effectors (
) dynamics Equations (8) and (9) we extend the model presented by Wodarz and Nowak in [
13] by differentiating infected T-cells with actively infected
and
cells (
).
If a
cells is infected it becomes either
latently or actively infected. The latently infected cells (
with
) can be activated and become actively infected (
with
) their activation rate is denoted by
. The actively infected cells are short living and will normally be killed upon activation with a high death rate
with
. Latently and actively infected T-cells dynamics are the following:
The positive terms
in Equation (
10) represent the number of healthy
-cells infected by sensitive viruses and healthy
-cells contaminated by actively infected macrophages (
), respectively. Observe that in the first one
can be counteracted by fusion and co-receptor inhibitors, in the second one drugs cannot obstruct
-action.
Negative terms is the number of infected cells that convert from latently to actively infected. This process is counteract by reverse transcriptase and integrase inhibitors which inhibit reverse transcription and integration of viral DNA respectively. This element appears with a positive sign in Equation (11) for actively -cells infected by sensitive viruses.
The negative contribute is a novelty of our model, it represents the number of actively infected naive T cells () killed by CTL effectors (), at a constant rate . We assume this number to be proportional to the number of immune effectors ().
Similar considerations can be observed for Equations (12) and (13) with the variation that there isn’t any drug that can counteract resistant viruses. Moreover, the action of CTL effectors occurs at a constant rate .
The dynamics for
infected cells are the following:
The dynamics for
infected cells are the following:
The dynamics for
infected cells are the following:
As said for naive T lymphocytes (
) the same holds for Equations (
14)–(17) of
cells, for Equations (
18)–(21) of
cells and for Equations (
22)–(25) of
cells.
The dynamics for
infected macrophages are the following:
We assumed that there is no latently infected macrophage population since the virus seems to replicate once inside them. We also assume that macrophages produce virus at a slow constant rate, sparing the host cell, so there is only natural death, not death by bursting like that for infected T cells [
14]. These equations present two infections rate, the former (
) related to sensitive viruses (
), the latter related to resistant viruses (
).
Finally, we assume that reverse transcriptase and protease inhibitors can counteract macrophages infected by sensitive viruses [
15,
16], while there are not drugs to counteract resistant viruses.
The dynamics for
sensitive () and resistant () viruses are the following:
The term:
in the sensitive viruses Equation (
28) is a source of virus that accounts for viral contributions to the plasma from both external compartments, such as the lymph system, as well as virus produced by infected cells in the plasma [
17]. This source is counteracted by the protease inhibitors therapy
.
The term:
represents the growth of sensitive viruses at a constant rate
. This growth is counteracted by protease inhibitors, and it is proportional to the number of actively infected cells contaminated by sensitive viruses and to the mutation rate from sensitive to other sensitive viruses (
). Finally:
represents the growth (at a constant rate
) of resistant viruses that mutate into sensitive viruses (at a rate
). Also this growth is counteracted by protease inhibitors, and it is proportional to the number of actively infected cells, contaminated by resistant viruses. Parameters
are assumed positive and constant.
The negative contributions:
indicate that viruses which infect lymphocytes are not free to circulate in the blood, and so they can not infect other cells at the same time. The constant
indicates the average number of sensitive virions infecting a cell: in our numerical simulations we set this number equal to 1, but the model permits also other values for this parameter, that is a cell could be infected by several viruses.
The last term:
represents the number of sensitive virions blocked by specific immune response (represented by
cells) and by unspecific immune response (represented by
cells).
For what concerns resistant viruses dynamics Equation (
28), considerations are the same just seen for sensitive viruses, with the only difference that protease inhibitors can not counteract
.
The
boundary conditions include the initial states:
and the final values:
Each therapy control, that represents the various drugs dosage, is assumed to vary within the range
the same holds for the HIV-controls, which represent the mutation rates of viruses, so that:
and recalling the control
definition, we have:
Summarizing, the payoffs Equation (
1) and Equation (
2) together with the associated state functions Equations (
3)–(29) and with the constraints Equation (
30)–Equation (
33) constitute a differential game characterized by:
8 controls (
Table 2): 6 controls for the therapy and 2 controls for the virus;