Next Article in Journal
From Cayley-Dickson Algebras to Combinatorial Grassmannians
Next Article in Special Issue
Modeling ITNs Usage: Optimal Promotion Programs Versus Pure Voluntary Adoptions
Previous Article in Journal
Construction of Periodic Wavelet Frames Generated by the Walsh Polynomials
Previous Article in Special Issue
The San Francisco MSM Epidemic: A Retrospective Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HIV vs. the Immune System: A Differential Game

1
Department of Mathematics, University of Padova, Via Trieste 63, Padova 35121, Italy
2
IOV-IRCCS, Via Gattamelata 64, Padova 35128, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2015, 3(4), 1139-1170; https://doi.org/10.3390/math3041139
Submission received: 17 July 2015 / Revised: 6 November 2015 / Accepted: 19 November 2015 / Published: 3 December 2015
(This article belongs to the Special Issue Optimal Control and Management of Infectious Diseases)

Abstract

:
A differential game is formulated in order to model the interaction between the immune system and the HIV virus. One player is represented by the immune system of a patient subject to a therapeutic treatment and the other player is the HIV virus. The aim of our study is to determine the optimal therapy that allows to prevent viral replication inside the body, so as to reduce the damage caused to the immune system, and allow greater survival and quality of life. We propose a model that considers all the most common classes of antiretroviral drugs taking into account different immune cells dynamics. We validate the model with numerical simulations, and determine optimal structured treatment interruption (STI) schedules for medications.

Graphical Abstract

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”, ( V s ) and “resistant” ( V r ) to the therapy.
In the model we adopt the idea of Herod et al. [8] and Nowak [9], assuming that the immune response ( T h ) to a viral infection creates some subpopulations of lymphocytes ( T s ) that are specific only for sensitive viruses, some others ( T r ) that are specific only for resistant viruses, and finally the subpopulation T g 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: c e l l s / m m 3 ).
Table 1. State functions (Unit of measure: c e l l s / m m 3 ).
SymbolDescription
I p Immune Precursor Cytotoxic T lymphocytes
I e Immune Effector Cytotoxic T lymphocytes
M h Uninfected macrophages
M h , j Macrophage infected by viruses j
T h Uninfected CD4+ T Helper cells
T s Uninfected CD4+ T Helper cells strain specific for V s
T r Uninfected CD4+ T Helper cells strain specific for V r
T g Uninfected CD4+ T Helper cells unspecific for V s and V r
T i , j l Latently infected CD4+ T i cells infected by virus j
T i , j a Actively infected CD4+ T i cells infected by virus j
V s HIV drugs sensitive viruses
V r HIV drugs resistant viruses
Table 2. Control functions (bounded in [ 0 , 1 ] ).
Table 2. Control functions (bounded in [ 0 , 1 ] ).
SymbolDescription
u c o Dosage of coreceptor antagonists
u f i Dosage of fusion inhibitors
u i b Immune boosting
u i i Dosage of integrase inhibitors
u p i Dosage of protease inhibitors
u r t Dosage of reverse transcriptase inhibitors
u r , r v Mutation rate from resistant virus ( V r ) to other resistant virus ( V r )
u r , s v Mutation rate from resistant virus ( V r ) to sensitive virus ( V s )
u s , r v Mutation rate from sensitive virus ( V s ) to resistant virus ( V r )
u s , s v Mutation rate from sensitive virus ( V s ) to other sensitive virus ( V s )
Table 3. Parameters (Unit of measure: 1 / d a y ).
Table 3. Parameters (Unit of measure: 1 / d a y ).
SymbolDescription
cCytotoxic T lymphocyte (CTL) activation rate
i j Death rate of cells infected by viruses V j due to immune response
kRate at which T cells convert to specific immune reaction cells ( T s or T r )
k Rate at which T cells convert to unspecific immune reaction cells
k 2 Rate at which latently infected cells convert to actively infected cells
μ i Natural death rate of type i cells
π s Growth rate of sensitive viruses ( V s )
π r Growth rate of resistant viruses ( V r )
qGrowth rate of I e due to infected cells and I p
rGrowth rate for CD4+ T cell population
Table 4. Parameters.
Table 4. Parameters.
SymbolDescriptionU.m./Value
A i Weight on the benefit i and cost of therapy i [ 0 , + ]
B i , j Weight of viral mutation from strain i to strain j [ 0 , + ]
b 2 Half saturation constant m m 3
gInput rate of external viral source c e l l s m m 3 × d a y
I 0 e Initial value of I e cells c e l l s m m 3
I 0 p Initial value of I p cells c e l l s m m 3
k 1 , j Rate at which viruses type j infect T cells m m 3 v i r i o n s × d a y
k 1 , j M Rate at which viruses type j infect M cells m m 3 v i r i o n s × d a y
k 3 Rate at which infected macrophages infect T cells m m 3 c e l l s × d a y
M h 0 Initial value of M h cells c e l l s m m 3
M h , j a , 0 Initial value of actively M h cells infected by viruses V j c e l l s m m 3
θ i Scaling parameter for type i cells
p j Specific immune response rate against viruses V j c e l l s m m 3 × d a y
Q i , j Percentage of i cells that recognize the virus j [ 0 , 1 ]
ρ j Average number of virions j infecting a cell v i r i o n s c e l l
s 1 , j Source/production of type j cells c e l l s m m 3 × d a y
s j Unspecific immune response rate against viruses V j c e l l s m m 3 × d a y
T i 0 Initial value of T i cells c e l l s m m 3
T i , j l , 0 Initial value of latently T i cells infected by viruses V j c e l l s m m 3
T i , j a , 0 Initial value of actively T i cells infected by viruses V j c e l l s m m 3
T m a x Maximum CD4+ T Helper cell population level c e l l s m m 3
V j 0 Initial value of viruses V j c e l l s m m 3
Let us collect all the state variables in the following array:
x ( t ) = ( T h ( t ) , T s ( t ) , T r ( t ) , T g ( t ) , M h ( t ) , I p ( t ) , I e ( t ) , T h , s l ( t ) , T h , s a ( t ) , T h , r l ( t ) , T h , r a ( t ) , T s , s l ( t ) , T s , s a ( t ) , T s , r l ( t ) , T s , r a ( t ) , T r , s l ( t ) , T r , s a ( t ) , T r , r l ( t ) , T r , r a ( t ) , T g , s l ( t ) , T g , s a ( t ) , T g , r l ( t ) , T g , r a ( t ) , M h , s a ( t ) , M h , r a ( t ) , V s ( t ) , V r ( t ) ) .
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 u s , s v and u r , r v only, observing that u s , s v = 1 - u s , r v and u r , r v = 1 - u r , s v .
Let’s collect the immune system-therapy controls in the array:
u 1 ( t ) = u c o ( t ) , u f i ( t ) , u i b ( t ) , u i i ( t ) , u p i ( t ) , u r t ( t )
and the viruses controls in the array:
u 2 ( t ) = u s , s v ( t ) , u r , r v ( t )
In Table 3 we show the parameters with the same unit of measure ( 1 / d a y ), 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:
J T h e r u 1 ( t ) = t 0 t f M h ( t ) + I p ( t ) + I e ( t ) + i Γ T i ( t ) - 1 2 k Ω A k u k 2 ( t ) d t
where Γ = { h , s , r , g } and Ω { c o , f i , i b , i i , p i , r t } .
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:
J H I V u 2 ( t ) = t 0 t f [ V s ( t ) + V r ( t ) + i Γ j Θ T i , j l ( t ) + T i , j a ( t ) - 1 2 B s , r 1 - u s , s v ( t ) 2 + B r , s 1 - u r , r v ( t ) 2 ] d t
where Γ = { h , s , r , g } and Θ = { s , r } . 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 x i x be a particular type of cell, such a decaying term is - μ x i x i and appears at the end of each considered dynamics.
In the evolution of the uninfected C D 4 + 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
T ˙ h ( t ) = s 1 , T h θ T h θ T h + V s ( t ) + V r ( t ) + r T h ( t ) 1 - T i ( t ) + T i , j l ( t ) + T i , j a ( t ) T m a x + u i b ( t ) T h ( t ) - k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T h , V s V s ( t ) T h ( t ) - k 1 , r 1 - Q T h , V r V r ( t ) T h ( t ) - k V s ( t ) T h ( t ) - k V r ( t ) T h ( t ) - k V s ( t ) + V r ( t ) T h ( t ) - k 3 M h , s a ( t ) + M h , r a ( t ) T h ( t ) - μ T h T h ( t )
T ˙ s ( t ) = k V s ( t ) T h ( t ) + u i b ( t ) T s ( t ) - k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T s , V s V s ( t ) T s ( t ) - k 1 , r 1 - Q T s , V r V r ( t ) T s ( t ) - k 3 M h , s a ( t ) + M h , r a ( t ) T s ( t ) - μ T s T s ( t )
T ˙ r ( t ) = k V r ( t ) T h ( t ) + u i b ( t ) T r ( t ) - k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T r , V s V s ( t ) T r ( t ) - k 1 , r 1 - Q T r , V r V r ( t ) T r ( t ) - k 3 M h , s a ( t ) + M h , r a ( t ) T r ( t ) - μ T r T r ( t )
T ˙ g ( t ) = k V s + V r ( t ) T h ( t ) + u i b ( t ) T g ( t ) - k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T g , V s V s ( t ) T g ( t ) - k 1 , r 1 - Q T g , V r V r ( t ) T g ( t ) - k 3 M h , s a ( t ) + M h , r a ( t ) T g ( t ) - μ T g T g ( t )
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 C D 4 + T-Helper cells and it depends on the control u i b . Note the positive term u i b ( t ) T i , i Γ = { h , s , r , g } that appears at the end of the first line in the dynamics of each type of cell T i .
Moreover, some cells are infected by sensitive viruses: Their number is proportional to the number of sensitive viruses ( V s ) and to the dimension of the community of cells itself. Coefficient k 1 , s represents the infection rate, the factors 1 - u f i ( t ) and 1 - u c o ( t ) 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 ( 1 - Q T i , V s ) 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 ( V r ) and to the dimension of the community of cells itself. What makes the difference is that fusion inhibitors and coreceptor antagonists can not counteract T i cells infection by resistant viruses ( V r ), therefore the drugs action is null. The factor ( 1 - Q T i , V r ) 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 - k V s ( t ) T h ( t ) - k V r ( t ) T h ( t ) - k V s ( t ) + V r ( t ) T h ( t ) which appear in the third line of Equation (3) represent the number of T h cells converted to T s , T r or T g cells respectively. Each one of such terms obviously appears with a positive sign in Equations (4)–(6) for uninfected T s , T r and T g cells.
The term - k 3 M h , s a ( t ) + M h , r a ( t ) gives the portion of T i 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 C D 4 + T-Helper cells, let’s explicate the terms which characterise the dynamics of the naive T lymphocytes ( T h ) in Equation (3).
The term:
s 1 , T h θ T h θ T h + V s ( t ) + V r ( t )
gives the proliferation of uninfected CD4+ T h 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+ T h cells differentiation.
The term:
r T h ( t ) 1 - T i ( t ) + T i , j l ( t ) + T i , j a ( t ) T m a x
represents the production of T h -cells due to cloning in the presence of an antigen, taking into account the maximum number of lymphocytes, T m a x .
The dynamics for the uninfected macrophages, C D 8 + CTL precursors and effectors are the following:
M ˙ h ( t ) = s 1 , M h - k 1 , s M 1 - u r t ( t ) 1 - u p i ( t ) V s ( t ) M h ( t ) - k 1 , r M V r ( t ) M h ( t ) - μ M h M h ( t )
I ˙ p ( t ) = c I p ( t ) T i ( t ) T i , j a ( t ) - c q I p ( t ) T i , j a ( t ) - μ I p I p ( t )
I ˙ e ( t ) = c q I p ( t ) T i , j a ( t ) - μ I e I e ( t )
According with [10], macrophages dynamics Equation (7) consider a constant proliferation source ( s 1 , M h ) and sensitive and resistant infections elements - k 1 , s M 1 - u r t ( t ) 1 - u p i ( t ) V s ( t ) M h ( t ) - k 1 , r M V r ( t ) M h ( t ) . 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 ( I p , I e ) dynamics Equations (8) and (9) we extend the model presented by Wodarz and Nowak in [13] by differentiating infected T-cells with actively infected T h , T s , T r and T g cells ( T i , j a ).
If a T h cells is infected it becomes either latently or actively infected. The latently infected cells ( T h , j l with j { s , r } ) can be activated and become actively infected ( T h , j a with j { s , r } ) their activation rate is denoted by k 2 . The actively infected cells are short living and will normally be killed upon activation with a high death rate μ T h , j a with j { s , r } . Latently and actively infected T-cells dynamics are the following:
T ˙ h , s l ( t ) = k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T h , V s V s ( t ) T h ( t ) + k 3 M h , s a ( t ) T h ( t ) - k 2 1 - u r t ( t ) 1 - u i i ( t ) T h , s l ( t ) - μ T h , s l T h , s l ( t )
T ˙ h , s a ( t ) = k 2 1 - u r t ( t ) 1 - u i i ( t ) T h , s l ( t ) - i s I e ( t ) T h , s a ( t ) - μ T h , s a T h , s a ( t )
 
T ˙ h , r l ( t ) = k 1 , r 1 - Q T h , V r V r ( t ) T h ( t ) + k 3 M h , r a ( t ) T h ( t ) - k 2 T h , r l ( t ) - μ T h , r l T h , r l ( t )
T ˙ h , r a ( t ) = k 2 T h , r l ( t ) - i r I e ( t ) T h , r a ( t ) - μ T h , r a T h , r a ( t )
The positive terms + k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T h , V s V s ( t ) T h ( t ) + k 3 M h , s a ( t ) T h ( t ) in Equation (10) represent the number of healthy T h -cells infected by sensitive viruses and healthy T h -cells contaminated by actively infected macrophages ( M h , s a ), respectively. Observe that in the first one V s can be counteracted by fusion and co-receptor inhibitors, in the second one drugs cannot obstruct M h , s a ( t ) -action.
Negative terms - k 2 1 - u r t ( t ) 1 - u i i ( t ) T h , s l ( t ) 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 T h -cells infected by sensitive viruses.
The negative contribute - i s I e ( t ) T h , s a ( t ) is a novelty of our model, it represents the number of actively infected naive T cells ( T h , s a ) killed by CTL effectors ( I e ), at a constant rate i s . We assume this number to be proportional to the number of immune effectors ( I e ).
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 i r .
The dynamics for infected T s cells are the following:
T ˙ s , s l ( t ) = k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T s , V s V s ( t ) T s ( t ) + k 3 M h , s a ( t ) T s ( t ) - k 2 1 - u r t ( t ) 1 - u i i ( t ) T s , s l ( t ) - μ T s , s l T s , s l ( t )
T ˙ s , s a ( t ) = k 2 1 - u r t ( t ) 1 - u i i ( t ) T s , s l ( t ) - i s I e ( t ) T s , s a ( t ) - μ T s , s a T s , s a ( t )
T ˙ s , r l ( t ) = k 1 , r 1 - Q T s , V r V r ( t ) T s ( t ) + k 3 M h , r a ( t ) T s ( t ) - k 2 T s , r l ( t ) - μ T s , r l T s , r l ( t )
T ˙ s , r a ( t ) = k 2 T s , r l ( t ) - i r I e ( t ) T s , r a ( t ) - μ T s , r a T s , r a ( t )
The dynamics for infected T r cells are the following:
T ˙ r , s l ( t ) = k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T r , V s V s ( t ) T r ( t ) + k 3 M h , s a ( t ) T r ( t ) - k 2 1 - u r t ( t ) 1 - u i i ( t ) T r , s l ( t ) - μ T r , s l T r , s l ( t )
T ˙ r , s a ( t ) = k 2 1 - u r t ( t ) 1 - u i i ( t ) T r , s l ( t ) - i s I e ( t ) T r , s a ( t ) - μ T r , s a T r , s a ( t )
T ˙ r , r l ( t ) = k 1 , r 1 - Q T r , V r V r ( t ) T r ( t ) + k 3 M h , r a ( t ) T r ( t ) - k 2 T r , r l ( t ) - μ T r , r l T r , r l ( t )
T ˙ r , r a ( t ) = k 2 T r , r l ( t ) - i r I e ( t ) T r , r a ( t ) - μ T r , r a T r , r a ( t )
The dynamics for infected T g cells are the following:
T ˙ g , s l ( t ) = k 1 , s 1 - u f i ( t ) 1 - u c o ( t ) 1 - Q T g , V s V s ( t ) T g ( t ) + k 3 M h , s a ( t ) T g ( t ) - k 2 1 - u r t ( t ) 1 - u i i ( t ) T g , s l ( t ) - μ T g , s l T g , s l ( t )
T ˙ g , s a ( t ) = k 2 1 - u r t ( t ) 1 - u i i ( t ) T g , s l ( t ) - i s I e ( t ) T g , s a ( t ) - μ T g , s a T g , s a ( t )
T ˙ g , r l ( t ) = k 1 , r 1 - Q T g , V r V r ( t ) T g ( t ) + k 3 M h , r a ( t ) T g ( t ) - k 2 T g , r l ( t ) - μ T g , r l T g , r l ( t )
T ˙ g , r a ( t ) = k 2 T g , r l ( t ) - i r I e ( t ) T g , r a ( t ) - μ T g , r a T g , r a ( t )
As said for naive T lymphocytes ( T h ) the same holds for Equations (14)–(17) of T s cells, for Equations (18)–(21) of T r cells and for Equations (22)–(25) of T g cells.
The dynamics for infected macrophages are the following:
M ˙ h , s a ( t ) = k 1 , s M 1 - u r t ( t ) 1 - u p i ( t ) V s ( t ) M h ( t ) - μ M h , s a M h , s a
M ˙ h , r a ( t ) = k 1 , r M V r ( t ) M h ( t ) - μ M h , r a M h , r a
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 ( k 1 , s M ) related to sensitive viruses ( V s ), the latter related to resistant viruses ( V r ).
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 ( V s ) and resistant ( V r ) viruses are the following:
V ˙ s ( t ) = g 1 - u p i ( t ) V s ( t ) b 2 + V s ( t ) + π s 1 - u p i ( t ) u s , s v ( t ) M h , s a ( t ) + T i , s a ( t ) + π r 1 - u p i ( t ) 1 - u r , r v ( t ) M h , r a ( t ) + T i , r a ( t ) - k 1 , s ρ s 1 - u f i ( t ) 1 - u c o ( t ) T i ( t ) 1 - Q T i , V s V s ( t ) - k 1 , s M ρ s 1 - u r t ( t ) 1 - u p i ( t ) M h ( t ) V s ( t ) - s s T g ( t ) + p s T s ( t ) V s ( t ) μ V s V s ( t ) - μ V s V s ( t )
V ˙ r ( t ) = g V r ( t ) b 2 + V r ( t ) + π s 1 - u s , s v ( t ) M h , s a ( t ) + T i , s a ( t ) + π r u r , r v ( t ) M h , r a ( t ) + T i , r a ( t ) - k 1 , r ρ r T i ( t ) 1 - Q T i , V r V r ( t ) - k 1 , r M ρ r M h V r ( t ) - s r T g ( t ) + p r T r ( t ) V r ( t ) - μ V r V r ( t )
The term:
g 1 - u p i ( t ) V s ( t ) b 2 + V s ( t )
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 u p i ( t ) .
The term:
π s 1 - u p i ( t ) u s , s v ( t ) M h , s a ( t ) + T i , s a ( t )
represents the growth of sensitive viruses at a constant rate π s . 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 ( u s , s v ). Finally:
π r 1 - u p i ( t ) 1 - u r , r v ( t ) M h , r a ( t ) + T i , r a ( t )
represents the growth (at a constant rate π r ) of resistant viruses that mutate into sensitive viruses (at a rate 1 - u r , r v = u r , s v ). Also this growth is counteracted by protease inhibitors, and it is proportional to the number of actively infected cells, contaminated by resistant viruses. Parameters b 2 , g , π s , π r are assumed positive and constant.
The negative contributions:
- k 1 , s ρ s 1 - u f i ( t ) 1 - u c o ( t ) T i ( t ) 1 - Q T i , V s V s ( t ) - k 1 , s M ρ s 1 - u r t ( t ) 1 - u p i ( t ) M h ( t ) V s ( t )
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 ρ s 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:
- s s T g ( t ) + p s T s ( t ) V s ( t )
represents the number of sensitive virions blocked by specific immune response (represented by T s cells) and by unspecific immune response (represented by T g 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 V r .
The boundary conditions include the initial states:
T h ( t 0 ) = T h 0 T s ( t 0 ) = T s 0 T r ( t 0 ) = T r 0
T g ( t 0 ) = T g 0 M h ( t 0 ) = M h 0 I p ( t 0 ) = I 0 p
I e ( t 0 ) = I 0 e T h , s l ( t 0 ) = T h , s l , 0 T h , s a ( t 0 ) = T h , s a , 0
T h , r l ( t 0 ) = T h , r l , 0 T h , r a ( t 0 ) = T h , r a , 0 T s , s l ( t 0 ) = T s , s l , 0
T s , s a ( t 0 ) = T s , a a , 0 T s , r l ( t 0 ) = T s , r l , 0 T s , r a ( t 0 ) = T s , r a , 0
T r , s l ( t 0 ) = T r , s l , 0 T r , s a ( t 0 ) = T r , s a , 0 T r , r l ( t 0 ) = T r , r l , 0
T r , r a ( t 0 ) = T r , r a , 0 T g , s l ( t 0 ) = T g , s l , 0 T g , s a ( t 0 ) = T g , s a , 0
T g , r l ( t 0 ) = T g , r l , 0 T g , r a ( t 0 ) = T g , r a , 0 M h , s a ( t 0 ) = M h , s a , 0
M h , r a ( t 0 ) = M h , r a , 0 V s ( t 0 ) = V s 0 V r ( t 0 ) = V r 0
and the final values:
T i ( t f ) , M h ( t f ) , I p ( t f ) , I e ( t f ) , T i , j l ( t f ) , T i , j a ( t f ) , V j ( t f ) R +
Each therapy control, that represents the various drugs dosage, is assumed to vary within the range [ 0 , 1 ] , the same holds for the HIV-controls, which represent the mutation rates of viruses, so that:
u 1 ( t ) [ 0 , 1 ] 6 , u 2 ( t ) [ 0 , 1 ] 2
and recalling the control u 2 ( t ) definition, we have:
u s , s v ( t ) + u s , r v ( t ) = 1 , u r , r v ( t ) + u r , s v ( t ) = 1
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:
  • 27 state variables (Table 1);
  • 8 controls (Table 2): 6 controls for the therapy and 2 controls for the virus;
  • 98 parameters (Table 3 and Table 4).

3. Open Loop Nash Equilibrium

In line with [6], we look for an Open Loop Nash equilibrium of the HIV/therapy game, considering the Pontryagin Maximum Principle necessary conditions. The model is difficult to analyze due to the number of variables and constraints and to the fact that state and co-state equations are “coupled”. For the sake of simplicity, let us relax the non negativity constraints Equation (31) on the state variables and verify them ex post, assuming that:
T i ( t f ) , M h ( t f ) , I p ( t f ) , I e ( t f ) , T i , j l ( t f ) , T i , j a ( t f ) , V j ( t f ) R
Theorem 1. 
[Therapy equilibrium necessary conditions] The optimal Therapy control satisfies the following necessary conditions:
u 1 * ( t ) = u c o * ( t ) , u f i * ( t ) , u 1 i b * ( t ) , u i i * ( t ) , u p i * ( t ) , u r t * ( t ) { 0 , 1 , u ^ 1 ( t ) }
where u ^ 1 ( t ) is the unique stationary point of the Therapy Hamiltonian Equation (A2). The co-states functions λ 1 ( t ) satisfy the co-state equations:
λ ˙ 1 , k ( t ) = - H T h e r x k x , u 1 , λ 1 , 1 , λ 1 , 2 , , λ 1 , 27 , t , k = 1 , , 27
with transversality conditions
λ 1 , k ( t f ) = 0 , k = 1 , , 27
The Therapy Hamiltonian is not necessarily concave w.r.t. its controls, nevertheless, some specific values of the model parameters lead to concavity, as the following Theorem states.
Theorem 2. 
[Therapy Hamiltonian Concavity sufficient conditions] If:
A c o A f i - h 1 , 2 2 ( t ) > 0
and:
A i b · A r t A p i A i i - h 4 , 6 2 ( t ) A p i - h 5 , 6 2 ( t ) A i i > 0
where h 1 , 2 ( t ) , h 4 , 6 ( t ) and h 5 , 6 ( t ) are given by Equations (A10)–(A12) respectively, then the Therapy Hamiltonian function is concave in u 1 .
Proof. 
See Appendix A. ☐
If concavity holds, the necessary conditions for the Therapy control Equation (35) takes the form:
u 1 * ( t ) = min 1 , m a x { 0 , u ^ 1 ( t ) }
where u ^ 1 ( t ) is the unique stationary point of the Therapy player.
Something more specific can be said for the HIV equilibrium strategies.
Theorem 3. 
[HIV equilibrium necessary conditions] The HIV mutation rates satisfy the following necessary conditions u 2 * ( t ) = ( u s , s v * ( t ) , u r , r v * ( t ) ) , where:
u s , s v * ( t ) = min { 1 , max { 0 , u ^ s , s v ( t ) } } , u r , r v * ( t ) = min { 1 , max { 0 , u ^ r , r v ( t ) } }
and:
u ^ s , s v ( t ) = B s , r - π s M h , s a ( t ) + T i , s a ( t ) λ 2 , 27 ( t ) + λ 2 , 26 ( t ) - 1 + u p i ( t ) B s , r
u ^ r , r v ( t ) = B r , s - π r M h , r a ( t ) + T i , r a ( t ) λ 2 , 27 ( t ) + λ 2 , 26 ( t ) - 1 + u p i ( t ) B r , s
The functions λ 2 , 26 ( t ) and λ 2 , 27 ( t ) are the co-states associated to the HIV virus variables V s and V r respectively. More generally the co-state equations are:
λ ˙ 2 , k ( t ) = - H H I V x k , x , u 2 , λ 2 , 1 , λ 2 , 2 , , λ 2 , 27 , t , k = 1 , , 27 ,
with transversality conditions:
λ 2 , k ( t f ) = 0 , k = 1 , , 27 .
Proof. 
See Appendix A. ☐
Observe that, due to the transversality conditions, it straight follows that u s , s v * ( t f ) = u r , r v * ( t f ) = 1 , so that by Equation (33) we obtain that in equilibrium condition the virus mutation rates from sensitive viruses to resistant viruses and vice versa vanish, that is u s , r v * ( t f ) = u r , s v * ( t f ) = 0 . Although it might appear as negative the fact that resistant viruses do not turn to sensitive any more, we must underline that this happens at the final time, when the number of viruses is sensibly decreased and a consequent long life expectancy has been obtained. In fact, our aim is to determine a therapy capable to prevent viral replication with its related side effects.
Substituting the optimal controls in the state Equations (3)–(29) and in the co-state Equation (36) and Equation (43), together with Equations (30), (37) and (44), we obtain a system of 81 autonomous strongly non linear first order differential equations, more precisely
  • 27 equations on the states;
  • 54 equations on the adjoint variables;
  • 27 boundary conditions for the states at time t 0 ;
  • 54 boundary conditions for the adjoint variables at time t f .

4. Numerical Simulations

An analytical solution of the state/co-state system can be obtained through a symbolic software, like for example, Wolfram Mathematica, nevertheless such a solution turns out to be as much long as uninformative. So that a numerical resolution turns out to be crucial in order to get some enlightening results. We formulate and apply an efficient and intuitive algorithm by combining the Forward-Backward Sweep algorithm [18,19,20,21] and a 4 th order Runge-Kutta routine (see e.g., [22]). We performed several simulations and we always got a unique solution to the necessary conditions. In what follows we present two particular instances, using the most common parameter values cited in the related literature, (see Appendix B). Please observe that differently from the theoretical analysis, here the HIV mutation rates u r , s v and u s , r v are explicitly shown, simply for an interpretation point of view.

4.1. Concave Instance

We first present a situation where the therapy Hamiltonian is concave, that is if conditions Equation (38) and Equation (39) hold. Such instance is characterized by very effective although highly toxic drugs which significantly reduce the virus mutation rates u r , s v and u s , r v .
In Figure 1, the first six graphs represent therapy dosages, the last two graphs show the HIV mutations, and the table shows the optimal payoffs related to immune system-therapy and to HIV, respectively.
Coreceptor antagonist ( u c o ) and fusion inhibitors ( u f i ) are both used for interfering with binding, fusion and entry of HIV to the host cell. Their dosages are the same for all therapy period (1000 days), more precisely they are administered at maximum doses for a few days and then they are no longer employed.
Immune boosting drugs ( u i b ) were introduced into this type of dynamics by Joshi in 2006 [12]. They represent the drugs that increase immune defenses enhancing the number of cells of immune system. In Joshi’s studies these drugs are used at a very low constant level (around 0 . 02 ) for the first 50 days. In our simulations, such a value turns out to be constantly equal to zero and this may have two explanations: first of all, it is not necessary to administrate a high toxicity drug if the lymphocyte counts remains high enough (around 800–900 c e l l s / m m 3 ), (see Figure 2). Secondarily, increasing the number of T lymphocytes means increasing the number of target cells for HIV virus, thus risking an increase in viral replication.
Figure 1. Equilibrium strategies and payoffs (concave instance).
Figure 1. Equilibrium strategies and payoffs (concave instance).
Mathematics 03 01139 g001
Figure 2. Lymphocytes, macrophage, and CTL dynamics (concave instance).
Figure 2. Lymphocytes, macrophage, and CTL dynamics (concave instance).
Mathematics 03 01139 g002
As we expected, Integrase ( u i i ) and Reverse Transcriptase Inhibitors ( u r t ) have very similar optimal dosages, because they are both used in the same phase of virions creation. To be more accurate, Reverse Transcriptase Inhibitors are used at the highest level for a few more days, in fact the interaction between the two enzymes, integrase (IN) and reverse transcriptase (RT) is functional and plays an important role in replication of the virus. Evidence has been accumulating that IN not only interacts closely with RT during reverse transcription but also remains associated with RT during the formation of the preintegrative complex.
Protease inhibitors drugs ( u p i ) block the viral protease enzyme necessary to produce mature virions upon budding from the host membrane. These drugs are used at maximum level for about the first 350 days, then the dosage decreases gradually becoming zero the last 100 days.
In fact it is well known that the resistance to protease inhibitors is a stepwise process. Briefly, primary mutations often have only a small effect on resistance, so that the PI therapy persists longer. Conversely, during continuous PI therapy, additional mutations emerge in the protease (commonly referred to as secondary mutations) leading to high-level PI resistance. These secondary mutations often have only a minor effect on drug susceptibility if they appear all alone but, in the presence of primary mutations, they can lead to a drastic increase in resistance, with the consequent decrease of the PI therapy, [23].
For what concerns HIV, the resistant–to–sensitive mutation rate from ( u r , s v ) is maximum for the first 350 days, then it drastically decreases, thanks to drug effects. At the same instant, the sensitive-to-resistant virus mutation rate ( u s , r v ), which is initially null, moves to the maximum level, because of the virus survival tendency. Nevertheless it does not remain at the maximum level, it begins to decrease, (around day 550) thanks to the protease inhibitors effect ( u p i ). At the end of the planning interval both the virus mutation rates vanish, as observed in Section 3 after Theorem 3. The Therapy and HIV payoffs are listed in the lower part of Figure 1. Observe that the high toxicity of drugs causes high costs for the immune system, and leads to low payoff values.
The white blood cell count is one of the most important tests to verify the presence of the HIV virus. It is an indicator of the health status of a patient and is one of the key factors in determining when to start HAART and prophylaxis for opportunistic infections. The tests usually measure the white blood cells present in the sample and calculate the percentage of CD4+ T lymphocytes. This value can be variable among different individuals and influenced by several factors. In general, the percentage of CD4+ T-cells is stable and it can be a parameter for evaluating the immune functions of a patient. Usually, it has a level between 500 and 1000 c e l l s / c m 3 in a healthy patient, and if it goes below this threshold, a possible infection could appear and damage the immune system. Observe that, despite the high toxicity of drugs, the lymphocyte count is sufficiently high (around 800–900 c e l l s / m m 3 , ) for the entire time interval (see Figure 2).
In Figure 2, Figure 3 and Figure 4 are shown the dynamics of all the 27 state variables of the model. In particular, in Figure 2 we find the dynamics of various types of healthy immune cells and in Figure 3 we find the infected cells dynamics. Finally in Figure 4 we find the dynamics of sensitive ( V s ) and resistant viruses ( V r ). Note that all optimal state functions are positive, so that the non-negativity constraints are fully satisfied. Sensitive viruses population decreases right from the start, under the effects of the therapy, and it becomes almost zero around day 400. Resistant viruses have the same behaviour and they become zero around day 800.
Figure 3. Infected cells dynamics (concave instance).
Figure 3. Infected cells dynamics (concave instance).
Mathematics 03 01139 g003
Figure 4. Sensitive and resistant viruses dynamics (concave instance).
Figure 4. Sensitive and resistant viruses dynamics (concave instance).
Mathematics 03 01139 g004

4.2. Nonconcave Instance

Let us now consider a nonconcave therapy instance, characterized by drugs with low toxicity but not very effective. In Figure 5, we can observe a prolonged intake over time is required. The performed algorithm is stable, in fact, results do not vary with the initial values given to adjoint variables at time t 0 (as it happens for other studies of the related literature, see e.g., [19,20]). Computational time is about 15–60 min in an Intel(R) core(TM) i-7 4.700 HQ CPU, 2.40 GHz.
HIV mutation rate u r , s v is maximum for the first 20 days then it decreases thanks to the combined action of therapy controls u c o , u f i and u p i . The reduction of u r , s v involves the immediate increase of the sensitive-to-resistant viruses mutation rate ( u s , r v ) to the maximum level. This causes the necessity to use the controls u i i and u r t at the highest level till day 250 in order to reduce to zero mutations u s , r v . Thanks to the prolonged therapy the latter remain at that null level until the end of the considered time interval. HIV attempts a further sting in the tail, by increasing the rate of mutation u r , s v to the highest level. In the meantime controls u i i and u r t are brought to zero for a small time interval, while the controls u c o , u f i and u p i continue to be used at their maximum level. Until the control u r , s v remains at a high level, the controls u i i and u r t perform a gradual increase.
We can observe how the combined action of 5 drugs leads to a decrease of the mutation rate u r , s v , around day 700. It continues to decrease for the remaining 300 days, after which both the viruses mutation rates are null. Immune boosting ( u i b ) continues to have a null dosage, for the same reasons seen in the concave instance.
The payoff for therapy and HIV are listed in the lower part of Figure 5. Due to the low toxicity of the drugs, the optimal immune system-therapy payoff turns to greater w.r.t. the one obtained in the concave instance, whilst the HIV payoff turns to be lower.
Figure 5. Equilibrium strategies and payoffs (nonconcave instance).
Figure 5. Equilibrium strategies and payoffs (nonconcave instance).
Mathematics 03 01139 g005
In Figure 6, Figure 7 and Figure 8 we report the dynamics of all 27 variable states. Note again that all optimal state functions are positive, so that the non-negativity constraints are fully satisfied. Although the drugs are not so effective, the lymphocytes count is high (around 800–900 c e l l s / m m 3 ). Sensitive viruses population ( V s ) becomes null around day 300, while resistant viruses population ( V r ) goes to zero after 900 days.
Figure 6. Lymphocytes, macrophage and CTL dynamics (nonconcave instance).
Figure 6. Lymphocytes, macrophage and CTL dynamics (nonconcave instance).
Mathematics 03 01139 g006
Figure 7. Infected cells dynamics (nonconcave instance).
Figure 7. Infected cells dynamics (nonconcave instance).
Mathematics 03 01139 g007
Figure 8. Sensitive and resistant viruses dynamics (nonconcave instance).
Figure 8. Sensitive and resistant viruses dynamics (nonconcave instance).
Mathematics 03 01139 g008

5. Conclusions

We have considered the game theory approach in a medical-therapeutic context for HIV treatment. We have formulated a differential game which describes the interaction between the immune system of a patient and the virus action. In order to model the efficient HAART therapy, we have considered 27 state functions and eight controls (six for the therapy and two for HIV). We look for an open Loop Nash equilibrium in order to find a stable long-term coexistence between the two players.
The optimality necessary conditions lead to a system of 81 first order strongly non-linear differential equations. Once proved the existence and uniqueness of the solution of such a system, we have performed some numerical simulations in order to get enlightening results. Just for explanatory examples, we considered two specific instances. Mathematically speaking the two instances differ in the concavity properties of the Hamiltonians w.r.t. the related controls. From the medical point of view, the concave instance is characterized by very effective although highly toxic drugs, whereas the nonconcave instance is represented by low effective and moderately toxic drugs. We observe quite different behaviours in the therapy strategies according to the type of drugs used. And this, in our opinion, is very interesting, as our aim is to find a tool which can support the clinician when combining the different drugs in the HAART therapy. Having an idea of the dynamical reactions to a given biological effect can help to make the optimal therapeutic decision.
An immediate, although not trivial, prospect could be the analysis of the problem following the Hamilton-Jacobi-Bellman approach to find a time consistent Markovian Nash equilibrium [24].
Another possible further extension of this model could be obtained in matching the Game Theory approach with the Fuzzy system issue [25,26,27,28], paying attention to the unavoidable complexity increase.

Acknowledgments

This research has been supported by the Italian Ministry of University and Research and the University of Padova.

Author Contributions

Alessandra Buratto and Rudy Cesaretto made the theoretical analysis; Rudy Cesaretto performed the numerical simulations; Rita Zamarchi wrote the introduction and took care of the immunological correctness of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix

A. Nash Equilibria

In this Appendix we report the proofs of Theorems 1 and 2. We first introduce the Hamiltonian functions for the two players in order to make the proofs easier to be read, see [29]. The Therapy Hamiltonian is:
H T h e r x , u 1 , λ 1 , 1 , λ 1 , 2 , , λ 1 , 27 , t = λ 1 , 0 M h ( t ) + I p ( t ) + I e ( t ) + i Γ T i ( t ) - 1 2 k Ω A k u k 2 + λ 1 , 1 { s 1 , T h θ T h θ T h + V s + V r + r T h 1 - T i + T i , j l + T i , j a T m a x + u i b T h - k 1 , s 1 - u f i 1 - u c o 1 - Q T h , V s V s T h - μ T h T h - k 1 , r 1 - Q T h , V r V r T h - k V s + V r T h - k V s + V r T h - k 3 M h , s a + M h , r a T h } + λ 1 , 2 { k V s T h + u i b T s - μ T s T s - k 1 , s 1 - u f i 1 - u c o 1 - Q T s , V s V s T s - k 1 , r 1 - Q T s , V r V r T s - k 3 M h , s a ( t ) + M h , r a T s } + λ 1 , 3 { k V r T h + u i b T r - μ T r T r - k 1 , s 1 - u f i 1 - u c o 1 - Q T r , V s V s T r - k 1 , r 1 - Q T r , V r V r T r - k 3 M h , s a + M h , r a T r } + λ 1 , 4 k V s + V r T h + u i b T g - μ T g T g + λ 1 , 5 s 1 , M h - μ M h M h - k 1 , s M 1 - u r t 1 - u p i V s M h - k 1 , r M V r M h + λ 1 , 6 c I p T i T i , j a - c q I p T i , j a - μ I p I p + λ 1 , 7 c q I p T i , j a - μ I e I e + λ 1 , 8 k 1 , s 1 - u f i 1 - u c o 1 - Q T h , V s V s T h + k 3 M h , s a T h ( t ) - μ T h , s l T h , s l - k 2 1 - u r t 1 - u i i T h , s l + λ 1 , 9 k 2 1 - u r t 1 - u i i T h , s l - μ T h , s a T h , s a - i s I e T h , s a + λ 1 , 10 k 1 , r 1 - Q T h , V r V r T h + k 3 M h , r a T h - μ T h , r l T h , r l - k 2 T h , r l + λ 1 , 11 k 2 T h , r l - μ T h , r a T h , r a - i r I e T h , r a + λ 1 , 12 k 1 , s 1 - u f i 1 - u c o 1 - Q T s , V s V s T s + k 3 M h , s a T s - μ T s , s l T s , s l - k 2 1 - u r t 1 - u i i T s , s l + λ 1 , 13 k 2 1 - u r t 1 - u i i T s , s l - μ T s , s a T s , s a - i s I e T s , s a + λ 1 , 14 k 1 , r 1 - Q T s , V r V r T s + k 3 M h , r a T s - μ T s , r l T s , r l - k 2 T s , r l + λ 1 , 15 k 2 T s , r l - μ T s , r a T s , r a - i r I e T s , r a + λ 1 , 16 k 1 , s 1 - u f i 1 - u c o 1 - Q T r , V s V s T r + k 3 M h , s a T r - μ T r , s l T r , s l - k 2 1 - u r t 1 - u i i T r , s l + λ 1 , 17 k 2 1 - u r t 1 - u i i T r , s l - μ T r , s a T r , s a - i s I e T r , s a + λ 1 , 18 k 1 , r 1 - Q T r , V r V r T r + k 3 M h , r a T r - μ T r , r l T r , r l - k 2 T r , r l + λ 1 , 19 k 2 T r , r l - μ T r , r a T r , r a - i r I e T r , r a + λ 1 , 20 k 1 , s 1 - u f i 1 - u c o 1 - Q T g , V s V s T g + k 3 M h , s a T g - μ T g , s l T g , s l - k 2 1 - u r t 1 - u i i T g , s l + λ 1 , 21 k 2 1 - u r t 1 - u i i T g , s l - μ T g , s a T g , s a - i s I e T g , s a + λ 1 , 22 k 1 , r 1 - Q T g , V r V r T g + k 3 M h , r a T g - μ T g , r l T g , r l - k 2 T g , r l + λ 1 , 23 k 2 T g , r l - μ T g , r a T g , r a - i r I e T g , r a + λ 1 , 24 k 1 , s M 1 - u r t 1 - u p i V s M h - μ M h , s a M h , s a + λ 1 , 25 k 1 , r M V r M h - μ M h , r a M h , r a + λ 1 , 26 { g 1 - u p i V s b 2 + V s + π s 1 - u p i u s , s v M h , s a + T i , s a + π r 1 - u p i 1 - u r , r v M h , r a + T i , r a - k 1 , s ρ s 1 - u f i 1 - u c o T i 1 - Q T i , V s V s - k 1 , s M ρ s 1 - u r t 1 - u p i M h V s - μ V s V s - s s T g + p s T s V s } + λ 1 , 27 { g V r b 2 + V r + π r u r , r v M h , r a + T i , r a + π s 1 - u s , s v M h , s a + T i , s a - k 1 , r ρ r T i 1 - Q T i , V r V r - k 1 , r M ρ r M h V r - μ V r V r - s r T g + p r T r V r }
The HIV Hamiltonian is
H H I V x , u 2 , λ 2 , 1 , λ 2 , 2 , , λ 2 , 27 , t = λ 2 , 0 V s + V r + i Γ j Θ T i , j l + T i , j a - 1 2 B s , r 1 - u s , s v 2 + B r , s 1 - u r , r v 2 + λ 2 , 1 { s 1 , T h θ T h θ T h + V s + V r + r T h 1 - T i + T i , j l + T i , j a T m a x + u i b T h - k 1 , s 1 - u f i 1 - u c o 1 - Q T h , V s V s T h - μ T h T h - k 1 , r 1 - Q T h , V r V r T h - k V s + V r T h - k V s + V r T h - k 3 M h , s a + M h , r a T h } + λ 2 , 2 { k V s T h + u i b T s - μ T s T s - k 1 , s 1 - u f i 1 - u c o 1 - Q T s , V s V s T s - k 1 , r 1 - Q T s , V r V r T s - k 3 M h , s a ( t ) + M h , r a T s } + λ 2 , 3 { k V r T h + u i b T r - μ T r T r - k 1 , s 1 - u f i 1 - u c o 1 - Q T r , V s V s T r - k 1 , r 1 - Q T r , V r V r T r - k 3 M h , s a + M h , r a T r } + λ 2 , 4 k V s + V r T h + u i b T g - μ T g T g + λ 2 , 5 s 1 , M h - μ M h M h - k 1 , s M 1 - u r t 1 - u p i V s M h - k 1 , r M V r M h + λ 2 , 6 c I p T i T i , j a - c q I p T i , j a - μ I p I p + λ 2 , 7 c q I p T i , j a - μ I e I e + λ 2 , 8 k 1 , s 1 - u f i 1 - u c o 1 - Q T h , V s V s T h + k 3 M h , s a T h ( t ) - μ T h , s l T h , s l - k 2 1 - u r t 1 - u i i T h , s l + λ 2 , 9 k 2 1 - u r t 1 - u i i T h , s l - μ T h , s a T h , s a - i s I e T h , s a + λ 2 , 10 k 1 , r 1 - Q T h , V r V r T h + k 3 M h , r a T h - μ T h , r l T h , r l - k 2 T h , r l + λ 2 , 11 k 2 T h , r l - μ T h , r a T h , r a - i r I e T h , r a + λ 2 , 12 k 1 , s 1 - u f i 1 - u c o 1 - Q T s , V s V s T s + k 3 M h , s a T s - μ T s , s l T s , s l - k 2 1 - u r t 1 - u i i T s , s l + λ 2 , 13 k 2 1 - u r t 1 - u i i T s , s l - μ T s , s a T s , s a - i s I e T s , s a + λ 2 , 14 k 1 , r 1 - Q T s , V r V r T s + k 3 M h , r a T s - μ T s , r l T s , r l - k 2 T s , r l + λ 2 , 15 k 2 T s , r l - μ T s , r a T s , r a - i r I e T s , r a + λ 2 , 16 k 1 , s 1 - u f i 1 - u c o 1 - Q T r , V s V s T r + k 3 M h , s a T r - μ T r , s l T r , s l - k 2 1 - u r t 1 - u i i T r , s l + λ 2 , 17 k 2 1 - u r t 1 - u i i T r , s l - μ T r , s a T r , s a - i s I e T r , s a + λ 2 , 18 k 1 , r 1 - Q T r , V r V r T r + k 3 M h , r a T r - μ T r , r l T r , r l - k 2 T r , r l + λ 2 , 19 k 2 T r , r l - μ T r , r a T r , r a - i r I e T r , r a + λ 2 , 20 k 1 , s 1 - u f i 1 - u c o 1 - Q T g , V s V s T g + k 3 M h , s a T g - μ T g , s l T g , s l - k 2 1 - u r t 1 - u i i T g , s l + λ 2 , 21 k 2 1 - u r t 1 - u i i T g , s l - μ T g , s a T g , s a - i s I e T g , s a + λ 2 , 22 k 1 , r 1 - Q T g , V r V r T g + k 3 M h , r a T g - μ T g , r l T g , r l - k 2 T g , r l + λ 2 , 23 k 2 T g , r l - μ T g , r a T g , r a - i r I e T g , r a + λ 2 , 24 k 1 , s M 1 - u r t 1 - u p i V s M h - μ M h , s a M h , s a + λ 2 , 25 k 1 , r M V r M h - μ M h , r a M h , r a + λ 2 , 26 { g 1 - u p i V s b 2 + V s + π s 1 - u p i u s , s v M h , s a + T i , s a + π r 1 - u p i 1 - u r , r v M h , r a + T i , r a - k 1 , s ρ s 1 - u f i 1 - u c o T i 1 - Q T i , V s V s - k 1 , s M ρ s 1 - u r t 1 - u p i M h V s - μ V s V s - s s T g + p s T s V s } + λ 2 , 27 { g V r b 2 + V r + π r u r , r v M h , r a + T i , r a + π s 1 - u s , s v M h , s a + T i , s a - k 1 , r ρ r T i 1 - Q T i , V r V r - k 1 , r M ρ r M h V r - μ V r V r - s r T g + p r T r V r }
The Hamiltonian functions of the two players are both almost twice differentiable with respect to the related controls and admit a unique stationary point. The maximization of the Hamiltonians is guarantee by Weierstrass Theorem, thanks to their continuity and the compactness of the feasible sets for the controls.
Proof of Theorem 1. 
The partial derivatives of the Therapy Hamiltonian w.r.t. the related controls are
H T h e r u c o = - A c o u c o - k 1 , s - 1 + u f i V s · [ ( - 1 + Q T g , V s ) λ 1 , 20 T g - ( - 1 + Q T g , V s ) λ 1 , 4 T g - ( - 1 + Q T h , V s ) λ 1 , 1 - λ 1 , 8 T h + ( - 1 + Q T r , V s ) λ 1 , 16 - λ 1 , 3 T s + ( - 1 + Q T s , V s ) λ 1 , 12 - λ 1 , 2 T s + ρ s λ 1 , 26 T g - Q T g , V s T g + T h - Q T h , V s T h + T r - Q T r , V s T r + T s - Q T s , V s T s ]
H T h e r u f i = - A f i u f i - k 1 , s - 1 + u c o V s · [ ( - 1 + Q T g , V s ) λ 1 , 20 T g - ( - 1 + Q T g , V s ) λ 1 , 4 T g - ( - 1 + Q T h , V s ) λ 1 , 1 - λ 1 , 8 T h + ( - 1 + Q T r , V s ) λ 1 , 16 - λ 1 , 3 T s + ( - 1 + Q T s , V s ) λ 1 , 12 - λ 1 , 2 T s + ρ s λ 1 , 26 T g - Q T g , V s T g + T h - Q T h , V s T h + T r - Q T r , V s T r + T s - Q T s , V s T s ]
H T h e r u i b = λ 1 , 4 T g + λ 1 , 1 T h + λ 1 , 3 T r + λ 1 , 2 T s - A i b u i b
H T h e r u i i = - A i i u i i - k 2 - 1 + u r t · λ 1 , 20 - λ 1 , 21 T g , s l + λ 1 , 8 - λ 1 , 9 T h , s l + λ 1 , 16 - λ 1 , 17 T r , s l + λ 1 , 12 - λ 1 , 13 T s , s l
H T h e r u p i = A p i u p i + k 1 , s M λ 1 , 24 - λ 1 , 5 M h - 1 + u r t V s + λ 1 , 26 { π r M h , r a + T g , r a + T h , r a + T r , r a + T s , r a · - 1 + u r , r v - π s u s , s v M h , r a + T g , r a + T h , r a + T r , r a + T s , r a + V s - k 1 , s M ρ s M h - 1 + u r t - g b 2 + V s }
H T h e r u r t = k 2 - 1 + u i i λ 1 , 20 - λ 1 , 21 T g , s l + λ 1 , 8 - λ 1 , 9 T h , s l + λ 1 , 16 - λ 1 , 17 T r , s l + λ 1 , 12 - λ 1 , 13 T s , s l - A r t u r t + k 1 , s M M h V s - 1 + u p i λ 1 , 24 - π s λ 1 , 26 - λ 1 , 5
There exists a unique Therapy stationary point, nevertheless we do not report it because of the length of its analytical formula.
The 6 × 6 Therapy Hessian matrix w.r.t. its controls is
H ˜ T h e r = - A c o h 1 , 2 ( t ) 0 0 0 0 h 1 , 2 ( t ) - A f i 0 0 0 0 0 0 - A i b 0 0 0 0 0 0 - A i i 0 h 4 , 6 ( t ) 0 0 0 0 - A p i h 5 , 6 ( t ) 0 0 0 h 4 , 6 ( t ) h 5 , 6 ( t ) - A r t
where
h 1 , 2 ( t ) = k 1 , s V s ( t ) { ( Q T h , V s - 1 ) λ 1 , 1 ( t ) - λ 1 , 8 ( t ) T h ( t ) - ( Q T s , V s - 1 ) λ 1 , 12 ( t ) - λ 1 , 2 ( t ) T s ( t ) - ( Q T r , V s - 1 ) λ 1 , 16 ( t ) - λ 1 , 3 ( t ) T r ( t ) - ( Q T g , V s - 1 ) λ 1 , 20 ( t ) - λ 1 , 4 ( t ) T g ( t ) + ρ s λ 1 , 26 ( t ) ( Q T h , V s - 1 ) T h ( t ) + ( Q T s , V s - 1 ) T s ( t ) + ( Q T r , V s - 1 ) T r ( t ) + ( Q T g , V s - 1 ) T g ( t ) }
h 4 , 6 ( t ) = k 2 [ - λ 1 , 12 ( t ) + λ 1 , 13 ( t ) T s , s l ( t ) + - λ 1 , 16 ( t ) + λ 1 , 17 ( t ) T r , s l ( t ) + - λ 1 , 20 ( t ) + λ 1 , 21 ( t ) T g , s l ( t ) + - λ 1 , 8 ( t ) + λ 1 , 9 ( t ) T h , s l ( t ) ]
h 5 , 6 ( t ) = k 1 , s M λ 1 , 24 ( t ) - ρ s λ 1 , 26 ( t ) - λ 1 , 5 ( t ) V s ( t ) M h ( t )
As the first leading principal minor ( - A c o ) is negative, the matrix cannot be in any case positively or semi-positively defined, so that the unique stationary point can be either a a saddle point or a local maximum. In the former case, one of the two extremes of the control feasible set is the maximum point. In the latter case, we must evaluate the Hamiltonian function not only in the extremes of the feasible control set but also in the stationary point, getting the point with the maximum objective value. In both case Equation (35) holds. The function evaluation necessary to determine the optimal Therapy control can be performed by mean of the numerical quadrature formulas with three different methods of numerical integration: the Composite rectangle method, the Composite trapezoid method; and the Composite Cavalieri-Simpson method, see [22]. ☐
Proof of Theorem 2. 
Let Equation (38) and Equation (39) hold, then the leading principal minors of the Therapy Hessian matrix are the following
2 th order
det H ˜ T h e r , 2 = A c o A f i - h 1 , 2 2 ( t ) > 0
3 rd order
- A i b · det H ˜ T h e r , 2 < 0
4 th order
A i i · A i b · det H ˜ T h e r , 2 > 0
5 th order
- A p i · A i i · A i b · det H ˜ T h e r , 2 < 0
Finally the determinant of the Hessian matrix, (principal minor of 6 th order), is
det ( H ˜ T h e r ) = det H ˜ T h e r , 2 · A i b · A r t A p i A i i - h 4 , 6 2 ( t ) A p i - h 5 , 6 2 ( t ) A i i > 0
The Hessian matrix is negative definite, viz the Therapy Hamiltonian function is concave in u 1 .
Proof of Theorem 3. 
The partial derivatives of the HIV Hamiltonian w.r.t. the related controls are
H H I V u s , s v = π s M h , s a ( t ) + T g , s a ( t ) + T h , s a ( t ) + T r , s a ( t ) + T s , s a ( t ) · λ 2 , 27 ( t ) + λ 2 , 26 ( t ) - 1 + u p i ( t ) - B s , r - 1 + u s , s v ( t )
H H I V u r , r v = π r M h , r a + T g , r a + T h , r a + T r , r a + T s , r a · λ 2 , 27 + λ 2 , 26 - 1 + u p i - B r , s - 1 + u r , r v
they admit a unique HIV stationary point, namely Equation (42). The HIV Hessian matrix is
H ˜ H I V = - B s , r 0 0 - B r , s
it is negative definite, because B s , r , B r , s > 0 , therefore the HIV Hamiltonian function is strictly concave in its controls u s , s v ( t ) and u r , r v ( t ) , and so the necessary conditions for the optimal controls are given by Equation (41). ☐

B. Parameters

The sets of parameter values we presented in the numerical Section 4 are very similar to the ones most commonly used in the related literature [5]. The two instances mainly differ from each other for the drugs toxicity parameters and HIV mutation costs, to consider concave and nonconcave properties. Table B1 shows the values we assigned to drugs toxicity and mutation costs for each instance. Table B2 and Table B3 show all common parameter values for both instances. Finally Table B4 contains boundary conditions.
Table B1. Weights.
Table B1. Weights.
SymbolDescriptionConcave InstanceNonconcave Instance
A c o Coreceptor antagonists toxicity5 10 - 2 . 7
A f i Fusion inhibitors toxicity5 10 - 2 . 7
A i b Immune boosting toxicity5 10 - 2 . 7
A i i Integrase inhibitors toxicity5 10 - 2 . 7
A p i Protease inhibitors toxicity5 10 - 2 . 7
A r t Reverse transcriptase inhibitors toxicity5 10 - 2 . 7
B s , r Mutation costs from sensitive to resistant viruses 10 - 10 10 - 13 . 7
B r , s Mutation costs from resistant to sensitive viruses 10 - 11 10 - 10 . 3
Table B2. Common parameters at both instances (Part 1).
Table B2. Common parameters at both instances (Part 1).
SymbolDescriptionValue
b 2 Half saturation constant 1 m m 3
cCytotoxic T Lymphocyte (CTL) activation rate 0 . 03 1 d a y
gInput rate of external viral source 25 c e l l s m m 3 × d a y
i s Death rate of cells infected by V s due to immune response 10 - 2 1 d a y
i r Death rate of cells infected by V r due to immune response 5 × 10 - 4 1 d a y
kRate T h convert to specific immune reaction cells ( T s or T r ) 10 - 6 1 d a y
k Rate T h convert to unspecific immune reaction cells ( T g ) 10 - 6 1 d a y
k 1 , s Rate CD4 T cells becomes infected by V s 0 . 0003 m m 3 v i r i o n s × d a y
k 1 , r Rate CD4 T cells becomes infected by V r 0 . 00005 m m 3 v i r i o n s × d a y
k 1 , s M Rate macrophages becomes infected by V s 0 . 00467 m m 3 v i r i o n s × d a y
k 1 , r M Rate macrophages becomes infected by V r 0 . 001 m m 3 v i r i o n s × d a y
k 2 Rate latently infected cells convert to actively infected cells 0 . 003 1 d a y
k 3 Rate infected macrophages infects CD4+ T cells 10 - 6 m m 3 c e l l s × d a y
μ T h Death rate of uninfected CD4+ T h cells 0 . 02 1 d a y
μ T s Death rate of uninfected CD4+ T s cells 0 . 02 1 d a y
μ T r Death rate of uninfected CD4+ T r cells 0 . 02 1 d a y
μ T g Death rate of uninfected CD4+ T g cells 0 . 02 1 d a y
μ M h Death rate of uninfected M h cells 0 . 005 1 d a y
μ I p Death rate of uninfected I p cells 0 . 001 1 d a y
μ I e Death rate of uninfected I e cells 0 . 1 1 d a y
μ T i , s l Death rate of latently infected T i , s l cells 0 . 1 1 d a y
μ T i , s a Death rate of actively infected T i , s a cells 0 . 24 1 d a y
μ T i , r l Death rate of latently infected T i , r l cells 0 . 2 1 d a y
μ T i , r a Death rate of actively infected T i , r a cells 0 . 5 1 d a y
μ M h , s a Death rate of actively infected M h , s a cells 0 . 04 1 d a y
μ M h , r a Death rate of actively infected M h , r a cells 0 . 07 1 d a y
μ V s Death rate of sensitive viruses V s 2 . 4 1 d a y
μ V r Death rate of resistant viruses V r 5 1 d a y
θ i Scaling parameter for type i cells1
Table B3. Common parameters at both instances (Part 2)
Table B3. Common parameters at both instances (Part 2)
SymbolDescriptionValue
π s Growth rate of sensitive viruses ( V s ) 8 1 d a y
π r Growth rate of resistant viruses ( V r ) 5 1 d a y
p s Specific immune response rate against V s 0 . 020 c e l l s m m 3 × d a y
p r Specific immune response rate against V r 0 . 001 c e l l s m m 3 × d a y
qGrowth rate of I e due to infected cells and I p 0 . 5 1 d a y
Q T h , V s Percentage of T h cells that recognize the virus V s 0
Q T h , V r Percentage of T h cells that recognize the virus V r 0
Q T s , V s Percentage of T s cells that recognize the virus V s 0.90
Q T s , V r Percentage of T s cells that recognize the virus V r 0
Q T r , V s Percentage of T r cells that recognize the virus V s 0
Q T r , V r Percentage of T r cells that recognize the virus V r 0.90
Q T g , V s Percentage of T g cells that recognize the virus V s 0
Q T g , V r Percentage of T g cells that recognize the virus V r 0
rrate of growth for CD4+ T Helper cell population 0 . 03 1 d a y
ρ s Average number of virions V s infecting a cell 1 v i r i o n s c e l l
ρ r Average number of virions V r infecting a cell 1 v i r i o n s c e l l
s 1 , T h Source/production T h cells 1 c e l l s m m 3 × d a y
s 1 , M h Source/production of macrophages 1 c e l l s m m 3 × d a y
s s Unspecific immune response rate against V s 9 . 7 c e l l s m m 3 × d a y
s r Unspecific immune response rate against V r 0 . 0009 c e l l s m m 3 × d a y
T m a x Maximum CD4+ T Helper cell population level 1700 c e l l s m m 3
Table B4. Common boundary conditions at both instances
Table B4. Common boundary conditions at both instances
SymbolDescriptionValue c e l l s m m 3
I 0 e Initial value of I e cells 10
I 0 p Initial value of I p cells 10
M h 0 Initial value of M h cells 1
T h 0 Initial value of T h cells 900
T s 0 Initial value of T s cells 100
T r 0 Initial value of T r cells 100
T g 0 Initial value of T g cells 1
T i , s l , 0 Initial value of latently T i cells infected by viruses V s 10
T i , s a , 0 Initial value of actively T i cells infected by viruses V s 0
T i , r l , 0 Initial value of latently T i cells infected by viruses V s 10
T i , r a , 0 Initial value of actively T i cells infected by viruses V r 0
M h , s a , 0 Initial value of actively M h cells infected by viruses V s 10
M h , r a , 0 Initial value of actively M h cells infected by viruses V r 10
V s 0 Initial value of viruses V s 500
V r 0 Initial value of viruses V r 300

References

  1. Moroni, M. AIDS, Se lo Conosci lo Eviti, 25 Anni di AIDS; Treccani: Firenze, Italy, 2007. [Google Scholar]
  2. Dybul, M.; Connors, M.; Fauci, A.; Paul, W. Immunology of HIV infection. In Fundamental Immunology, 5th ed.; Lippincott Williams & Wilkins: Philadelphia, PA, USA, 2003; pp. 1285–1318. [Google Scholar]
  3. Nachega, J.B.; Marconi, V.C.; van Zyl, G.U.; Gardner, E.M.; Preiser, W.; Hong, S.Y.; Mills, E.J.; Gross, R. HIV treatment adherence, drug resistance, virologic failure: Evolving concepts. Infect. Disord. Drug Targets 2011, 11, 167–174. [Google Scholar] [CrossRef] [PubMed]
  4. Joly, M.; Pinto, J.M. Role of mathematical modeling on the optimal control of HIV-1 pathogenesis. AIChE J. 2006, 52, 856–884. [Google Scholar] [CrossRef]
  5. Buratto, A.; Cesaretto, R.; Zamarchi, R. HIV Inside the Body: A Review of Mathematical Modeling; Padua University: Padua, Italy, 2015. [Google Scholar]
  6. Wu, J.; Zhang, M. A game theoretical approach to optimal control of dual drug delivery for HIV infection treatment. IEEE Trans. Syst. Man Cybern. B Cybern. 2010, 40, 694–702. [Google Scholar] [PubMed]
  7. Coppola, V.; Camus, O. Preventing without stigmatizing: The complex stakes of information on AIDS. Patient Educ. Couns. 2007, 67, 255–260. [Google Scholar] [CrossRef] [PubMed]
  8. Shonkwiler, R.W.; Herod, J. Mathematical Biology: An Introduction with Maple and Matlab; Springer Science & Business Media: New York, NY, USA, 2009; Volume 73. [Google Scholar]
  9. Nowak, M.A. Evolutionary Dynamics; Harvard University Press: Cambridge, MA, USA, 2006. [Google Scholar]
  10. Kirschner, D.E.; Webb, G. Immunotherapy of HIV-1 infection. J. Biol. Syst. 1998, 6, 71–83. [Google Scholar] [CrossRef]
  11. Caetano, M.A.; Yoneyama, T. Short and long period optimization of drug doses in the treatment of AIDS. An. Acad. Bras. Cienc. 2002, 74, 379–392. [Google Scholar] [CrossRef] [PubMed]
  12. Joshi, H.R. Optimal control of an HIV immunology model. Optim. Control Appl. Methods 2002, 23, 199–213. [Google Scholar] [CrossRef]
  13. Wodarz, D.; Nowak, M.A. Specific therapy regimes could lead to long-term immunological control of HIV. Proc. Natl. Acad. Sci. USA 1999, 96, 14464–14469. [Google Scholar] [CrossRef] [PubMed]
  14. Kirschner, D.; Perelson, A. A model for the immune system response to HIV: AZT treatment studies. Math. Popul. Dyn. Anal. Heterog. 1995, 1, 295–310. [Google Scholar]
  15. Wein, L.M.; Zenios, S.A.; Nowak, M.A. Dynamic multidrug therapies for HIV: A control theoretic approach. J. Theor. Biol. 1997, 185, 15–29. [Google Scholar] [CrossRef] [PubMed]
  16. Wein, L.M.; D’Amato, R.M.; Perelson, A.S. Mathematical analysis of antiretroviral therapy aimed at HIV-1 eradication or maintenance of low viral loads. J. Theor. Biol. 1998, 192, 81–98. [Google Scholar] [CrossRef] [PubMed]
  17. Haase, A.T.; Henry, K.; Zupancic, M.; Sedgewick, G.; Faust, R.A.; Melroe, H.; Cavert, W.; Gebhard, K.; Staskus, K.; Zhang, Z.Q.; et al. Quantitative image analysis of HIV-1 infection in lymphoid tissue. Science 1996, 274, 985–989. [Google Scholar] [CrossRef] [PubMed]
  18. Lenhart, S.; Workman, J.T. Optimal Control Applied to Biological Models; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  19. Velichenko, V. A numerical method for solving optimal control problems. USSR Comput. Math. Math. Phys. 1966, 6, 34–50. [Google Scholar] [CrossRef]
  20. Velichenko, V.; Pritykin, D. Control of the medical treatment of AIDS. Autom. Remote Control 2006, 67, 493–511. [Google Scholar] [CrossRef]
  21. Velichenko, V.; Pritykin, D. Numerical methods of optimal control of the HIV-infection dynamics. J. Comput. Syst. Sci. Int. 2006, 45, 894–905. [Google Scholar] [CrossRef]
  22. Quarteroni, F.S.A.; Saleri, F. Calcolo Scientifico; Springer: New York, NY, USA, 2008. [Google Scholar]
  23. Niu, P.; Lyonsl, N.; Hsuz, A.; Grannemanz, G.R.; Hos, D.D. Ordered accumulation of mutations in HIV protease confers resistance to ritonavir. Nat. Med. 1996, 2, 760–766. [Google Scholar]
  24. Dockner, E. Differential Games in Economics and Management Science; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  25. Dalal, N.; Greenhalgh, D.; Mao, X. A stochastic model for internal HIV dynamics. J. Math. Anal. Appl. 2008, 341, 1084–1101. [Google Scholar] [CrossRef] [Green Version]
  26. Grégio, J.M.; Caetano, M.A.; Yoneyama, T. State estimation and optimal long period clinical treatment of HIV seropositive patients. An. Acad. Bras. Cienc. 2009, 81, 3–12. [Google Scholar] [CrossRef] [PubMed]
  27. Tan, W.Y.; Xiang, Z. Some state space models of HIV pathogenesis under treatment by anti-viral drugs in HIV-infected individuals. Math. Biosci. 1999, 156, 69–94. [Google Scholar] [CrossRef]
  28. Tuckwell, H.C.; Le Corfec, E. A stochastic model for early HIV-1 population dynamics. J. Theor. Biol. 1998, 195, 451–463. [Google Scholar] [CrossRef] [PubMed]
  29. Seierstad, A.; Sydsaeter, K. Optimal Control Theory with Economic Applications; North-Holland: Amsterdam, The Netherlands, 1987; Volume 20. [Google Scholar]

Share and Cite

MDPI and ACS Style

Buratto, A.; Cesaretto, R.; Zamarchi, R. HIV vs. the Immune System: A Differential Game. Mathematics 2015, 3, 1139-1170. https://doi.org/10.3390/math3041139

AMA Style

Buratto A, Cesaretto R, Zamarchi R. HIV vs. the Immune System: A Differential Game. Mathematics. 2015; 3(4):1139-1170. https://doi.org/10.3390/math3041139

Chicago/Turabian Style

Buratto, Alessandra, Rudy Cesaretto, and Rita Zamarchi. 2015. "HIV vs. the Immune System: A Differential Game" Mathematics 3, no. 4: 1139-1170. https://doi.org/10.3390/math3041139

APA Style

Buratto, A., Cesaretto, R., & Zamarchi, R. (2015). HIV vs. the Immune System: A Differential Game. Mathematics, 3(4), 1139-1170. https://doi.org/10.3390/math3041139

Article Metrics

Back to TopTop