BioMed Research International

Volume 2014 (2014), Article ID 472869, 9 pages

http://dx.doi.org/10.1155/2014/472869

Research Article

Computational Study to Determine When to Initiate and Alternate Therapy in HIV Infection

Matthias Haering,¹ Andreas Hördt,² Michael Meyer-Hermann,1,3 and Esteban A. Hernandez-Vargas¹

¹Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Inhoffenstraße 7, 38124 Braunschweig, Germany

²Institut für Geophysik und extraterrestrische Physik, Technische Universität Braunschweig, Mendelssohnstraße 3, 38106 Braunschweig, Germany

³Institute for Biochemistry, Biotechnology and Bioinformatics, Braunschweig University of Technology, Braunschweig, 38106 Braunschweig, Germany

Received 15 February 2014; Revised 7 April 2014; Accepted 10 April 2014; Published 11 May 2014

Academic Editor: Filippo Castiglione

Copyright © 2014 Matthias Haering et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

HIV is a widespread viral infection without cure. Drug treatment has transformed HIV disease into a treatable long-term infection. However, the appearance of mutations within the viral genome reduces the susceptibility of HIV to drugs. Therefore, a key goal is to extend the time until patients exhibit resistance to all existing drugs. Current HIV treatment guidelines seem poorly supported as practitioners have not achieved a consensus on the optimal time to initiate and to switch antiretroviral treatments. We contribute to this discussion with predictions derived from a mathematical model of HIV dynamics. Our results indicate that early therapy initiation (within 2 years postinfection) is critical to delay AIDS progression. For patients who have not received any therapy during the first 3 years postinfection, switch in response to virological failure may outperform proactive switching strategies. In case that proactive switching is opted, the switching time between therapies should not be larger than 100 days. Further clinical trials are needed to either confirm or falsify these predictions.

1. Introduction

According to the last report from UNAIDS in 2011 [1], 34 million people live with human immunodeficiency virus (HIV). Although the number of new cases of HIV infection is declining, the number of people living with HIV is increasing; therefore the problems of continuing treatment of chronic infection are of major importance in today's social medicine.

Nowadays, the drugs to treat HIV type 1 (HIV-1) infection belong to four distinct classes [2, 3]: reverse transcriptase inhibitors, protease inhibitors, integrase inhibitors, and fusion inhibitors (Figure 1). Currently, highly active antiretroviral therapies (HAART) generally comprise three different drugs to combat different parts of the HIV cycle. These therapies prevent immune deterioration, reduce morbidity and mortality, and prolong the life expectancy of people infected with HIV.

Figure 1: HIV infection cycle affected by the four distinct drug classes. The drug classes are shown in red boxes.

Nevertheless, HAART is not always successful. Many patients have long-term complications while others experience virological failure (inability to maintain HIV RNA levels below 50 copies/mL) [3]. In most cases, viral rebound is associated with the emergence of resistance-conferring mutations within the viral genome, resulting in reduced viral susceptibility to one or more of the drugs. This is related to the reverse transcription process of viral RNA into DNA, which is highly prone to errors, introducing on average one mutation for each viral genome transcribed [2].

The primary goal of the initial regimen proposed in the guidelines for the use of antiretroviral agents in HIV-1 infected adults and adolescents by Department of Health and Human Services (DHHS) [3] is to suppress viral replication to the maximum and to sustain this level of suppression as long as possible. Even though HAART can reduce the viral load in the blood by at least five orders of magnitude, ongoing low-level replication still occurs; hence, the risk of developing resistance is always present. Complete virus eradication by HAART does not appear to be achievable in the foreseeable future. Furthermore, antiretroviral guidelines [3] have not achieved a consensus on two fundamental questions in HIV treatment: (i) when to start antiretroviral treatment and (ii) when to change for a new antiretroviral treatment.

The standard procedure for therapy initiation is still a point of discussion in the expert panel over the last 20 years [3]. However, there is a general consensus that antiretroviral therapy should be initiated in all patients with a history of an AIDS-defining illness or when CD4+ T cell counts are less than 350 cells/mm³. It is strongly recommended to start therapy if the CD4+ T cell count is between 350 and 500 cells/mm³ and there is a recommendation with moderate urgency for patients with CD4+ T cell counts > 500 cells/mm³ (Table 1).

Table 1: Panel recommendations for therapy initiation.

Urgency of initiation

CD4+ T cell count
<350 cells/mm³
Strong I

350–500 cells/mm³
Strong II

>500 cells/mm³
Moderate III

Transmission risk
Perinatal
Strong I

Heterosexual
Strong I

Other
Strong III

The urgency is rated by validity: I: data from randomized trials, II: data from well-designed nonrandomized trials or observational cohort studies with long-term clinical outcomes, III: expert opinion [3].

The expert panel [3] pointed out the absence of cohort studies that conclusively demonstrate a clinical benefit of HAART in patients with CD4+ T count > 350 cells/mm³. For some patients, the potential risks of short- or long-term drug-related complications and nonadherence to long-term therapy may offset possible benefits of earlier therapy initiation.

Computational biology may play an important role in evaluating the impact of the initiation time of HAART during HIV infection. Using a mathematical model, [4] suggested to initiate treatment in stages of the infection when the viral load can be easily controlled: in the acute phase of the infection when the viral load peaks and moderately in the asymptomatic phase. The very early phase and the AIDS phase are considered hardly controllable. Thus, authors in [4] argued in favour of an early but not immediate treatment initiation. The main drawback of this study is the lack of a mathematical model that is able to reproduce the whole disease trajectory, which limits the long-term assessments of treatment strategies. The observation that HAART timing has a strong impact on the disease outcome is supported by computational results by [5, 6]. However, these studies have similar drawbacks, as they do not envisage long-term dynamics (more than 8 years when AIDS may appear) and different treatment protocols.

For the second fundamental question in HIV treatment, when treatment should be alternated in response to virological failure, the answer is more complex [3, 7, 8]. One clinical goal is to delay the time until patients exhibit strains resistant to all existing regimens. There is a crucial trade-off between switching therapies. On the one hand, switching early carries the risk of poor adherence to a new drug regimen and prematurely exhausting the limited number of remaining salvage therapies. On the other hand, switching drugs too late allows the accumulation of mutations that leads to multidrug resistance [2, 3]. The most aggressive approach would be to change therapies for any repeated detectable viremia (e.g., two consecutive HIV RNA > 50 copies/mL after suppression). The most conservative strategy has been to allow detectable viremia up to an arbitrary level (e.g., 1000–500 copies/mL). This latter approach is called switch on virological failure [3].

Recently, Martinez-Picado et al. [9] suggested that switching between therapies can decrease the likelihood of viral resistance and prolong the pre-AIDS period. This hypothesis was supported by clinical trials [9] called Switching Antiviral Therapy Combination against HIV (SWATCH), which consists of two HAART regimens that are periodically alternated every 3 months. The alternation of 2 drug regimens could inhibit the emergence of highly resistant genotypes.

In this paper, we focus on the trade-off between the inhibition of long-term reservoirs and the promotion of resistant genotypes and thus help to precisely work out the advantages of early treatment. We also investigate how frequently antiretroviral regimens should be alternated in order to attain optimal proactive treatments.

2. Material and Methods

2.1. Mathematical Model

A typical HIV infection response consists of three stages: an initial acute infection, a long asymptomatic period, and a final increase in viral load with simultaneous collapse in healthy CD4+ T cell counts. The majority of existing mathematical models [4, 11–20] give a good representation of either the first two stages [11, 13, 21–24] or the last stage of the infection [25] but do not describe the three stages observed in HIV infection in a single framework. A mathematical model that is able to represent the typical HIV infection dynamics including all three stages was suggested by [26]. The model includes the possibility of full parameter variations, without losing the capability to describe the three stages.

The results in [26] indicate that HIV infection can be considered as two feedback systems (Figure 2). One provides the fast dynamics presented in the early stages of infection as a result of the fast infection process of CD4+ T cells. The second feedback sustains a constant slow infection process in macrophages over the years accompanied by the long time survival conditions of macrophages.

Figure 2: Nonlinear HIV model. represents the uninfected CD4+ T cells, represents the infected CD4+ T cells with the th strain, represents uninfected macrophages, represents infected macrophages with the th strain, represents the th strain, and is the sum of all strains.

Here, following the work by [8], we extend the mutation tree from [10] into a nonlinear model with mutations (1)–(5). Unlike existing models [4, 8, 11–20], our model is able to adequately represent the three stages of HIV infection and the dynamics of resistant genotypes when HAART treatment is introduced. The model is defined by the following set of differential equations:

where represents the uninfected CD4+ T cells, represents the infected CD4+ T cells with the th strain (genetic variant or subtype of the virus), represents uninfected macrophages, represents infected macrophages with the th strain, represents the th strain, and is the sum of all strains.

Parameters and represent the source terms of new CD4+ T cells and macrophages, respectively. HIV, as other pathogens, triggers the proliferation of immune cells. Homeostatic proliferation is modeled using a logistic growth model limited by viral load. This would allow convergence to a high percentage of the reservoirs (macrophages) being infected without allowing the total population to expand at unrealistic growth rates. Parameters and are the maximum proliferation rate for CD4+ T cells and macrophages, respectively. and represent the respective half-velocity constants.

The infection rate constant is represented with for CD4+ T cells and for macrophages. Viral proliferation is achieved in infected CD4+ T cells and infected macrophages with rate constants and , respectively. These parameters depend on the fitness of the genotype () and the therapy () that is being used.

We consider two therapies () composed of reverse transcriptase inhibitors with effectiveness for CD4+ T cells and for macrophages and protease inhibitors with effectiveness for CD4+ T cells and for macrophages. Note that, based on clinical evidence [27], inhibitors are more effective in CD4+ T cells than in macrophages; this is considered by and . The therapy parameter values used in this study are given in Figure 3.

Figure 3: Four variant mutation trees. The wild type (WT) is susceptible to both therapies. Genotype 1 (G1) is susceptible to therapy 2 and genotype 2 (G2) is susceptible to therapy 1. The highly resistant genotype (HRG) is not affected by any therapy. Parameter values were taken from [10].

The mutation rate is expressed by , and represents the genetic connections between genotypes. The degradation rates for the relevant species are , , , , and . Parameter values are presented in Table 2.

Table 2: Parameter values for (1)–(5).

Parameter
Units
Nominal value
Source

Cells/mm³ day
10
[28]

Cells/mm³ day
0.15
[28]

mm³/day copies
4.57 × 10⁻⁵
[28]

mm³/day copies
4.33 × 10⁻⁸
[26]

Copies/cell day
38
[29]

Copies/cell day
35
[29]

Day⁻¹
0.01
[29]

Day⁻¹
0.4
[29]

Day⁻¹
1 × 10⁻³
[26]

Day⁻¹
1 × 10⁻³
[26]

Day⁻¹
2.4
[29]

Day⁻¹
0.01
[26]

Day⁻¹
0.003
[26]

Copies/mm³
300
[26]

Copies/mm³
220
[26]

The meaning of the different terms and parameters is illustrated in a block diagram (Figure 2). For example, the left-hand side of (1) defines the rate of change of the uninfected CD4+ T cells. The right-hand side defines how CD4+ T cells are altered during infection: a constant production of CD4+ T cells is expressed by . CD4+ T cell depletion is proportional to the death rate and the current number of healthy CD4+ T cells. The cell proliferation is proportional to the number of healthy cells and the total viral load , which represents the sum over all strains (). The infection is described by the sum term on the right-hand side of (1) and is also proportional to the existing number of uninfected cells and the strains . Each strain has its own infection rate constant , which is why we have to sum over index . The other equations may be understood in a similar fashion.

2.2. Viral Mutation Tree

Nowadays, it is considered critical to take viral mutation into account during the development of treatment strategies. The process of reverse transcription is extremely error-prone and it is during this step that mutations can occur. High levels of resistance can be produced by substitutions of a single amino acid [3]. For instance, when lamivudine is used as a single agent, resistant strains will appear in a few weeks [2]. This is the reason why monotherapy has been discontinued and HAART is composed of at least three different drugs. As a result, multiple mutations are required for resistance to occur to all drugs in one regimen.

As a simple motivating example, we consider the mutation tree with 4 variants and 2 possible antiretroviral treatments proposed in [10]. The wild type genotype (WT) would be the most prolific variant in the absence of any drugs (Figure 3). However, it is also the variant that all drug combinations have been designed to combat and therefore is susceptible to all therapies. The wild type genotype (WT) can mutate to either genotype 1 (G1) which is susceptible to therapy 2 or genotype 2 (G2) which is susceptible to therapy 1. After mutations, the highly resistant genotype (HRG) is a genotype with a low proliferation rate, but resistant to all drug therapies.

2.3. Therapy

The DHHS panel [3] recommends therapies with two nucleoside analogues and either protease inhibitors or nucleoside reverse transcriptase inhibitors. The combination of these drugs is crucial in controlling the development of resistance. For simulation purposes, we consider the common clinical strategy suggested by the DHHS panel in [3], which recommends change therapy when virological failure is presented, that is, switch on virological failure (SVF): introduce a new regimen if there is detectable viremia (HIV RNA > 1,000 copies/mL) and drug-resistant genotypes are identified. We also consider the SWATCH strategy proposed in [9]. The rationale behind this strategy is that one could preempt virologic rebound and reduce accumulating drug-resistant genotypes by alternating treatments. This strategy is implemented as follows:

 SWATCH: alternate between two regimens every 3 months.

2.4. Monte Carlo Simulations

As the interpatient variability and the extreme error sensitivity of the HIV replication process are nonnegligible in HIV treatment, we perform Monte Carlo (MC) simulations to analyse treatment strategies in stochastic environments [30]. Our MC algorithm consists of repeated random sampling, providing numerical results that can be interpreted with statistical methods. To this end, we consider 1000 repeated simulations with randomly perturbed parameters (normal distribution and 30% deviation from the nominal parameter values in Table 2). The described MC simulations without therapy, with SVF, and SWATCH treatments were carried out at six different initiation times (0.5, 1, 1.5, 2, 3, and 4 years). The differential equations (1)–(5) are solved with the toolbox from the MATLAB library. We record the year when immunological failure appears (), that is, when CD4+ T cells sink under 200 copies/mm³. represents the number of cases where no immunological failure occurs during 30 years. The MC procedure is illustrated as a flow diagram in Figure 4. The analysis of the MC simulation results is based on a two-way ANOVA test. The significance is identified as value ; the data are further analyzed by a two-way -test and a Bonferroni posttest.

Figure 4: Monte Carlo simulation algorithm. The Simulation box includes the simulation of the four-genotype mutation model under a certain treatment. For each of the six times of treatment initiation, the simulation will be repeated 1000 times and the data can be analysed.

3. Results and Discussion

For the scenario with no treatment, we obtain an average time of immunological failure of 8.5 years postinfection from our MC simulations. In addition, 99% of the simulated cases may progress to AIDS at different time scales. These results are consistent with clinical observations in [31], suggesting that model (1)–(5) can adequately represent the basic features of HIV infection.

3.1. Treatment Initiation

The treatment strategies SVF and SWATCH (Figure 5) were implemented in model (1)–(5) using different initiation times (0.5, 1, 1.5, 2, 3, and 4 years) and parameter values as presented in Table 2. During the first 5 years, we observe that both SVF and SWATCH satisfy the levels required for healthy immunological responses (CD4+ T > 500 cells/mm³) and decrease the viral load to undetectable levels (<50 copies/mL) as suggested in [3].

Figure 5: SVF and SWATCH treatment strategies. CD4+ T cells and viral load are displayed for different initiation times (). The top two panels show simulation results for the SVF strategy and the bottom panels show the results for the SWATCH strategy. Deterministic simulations are based on nominal parameter values from Table 2.

For both treatment strategies, our results support the hypothesis that an early treatment initiation is more beneficial than a late start of the therapy. However, simulation results for the SVF strategy (Figure 5) suggest that there is only little improvement for a very early treatment initiation (<1 year).

If a second treatment is introduced, a second virological failure may appear, faster than the first one, progressing to AIDS (Figure 5). This is consistent with the observation in [32] that persistent low-level viremia and long-term reservoirs promote a second virological failure.

The clinical trial SWATCH [9] suggested an alternative solution to minimize HIV resistance mutation by alternating between two regimens every three months while viral load is suppressed. In our simulations, the alternation is visible in form of the high-frequency oscillations. For early initiation treatment (<2 years), the results in Figure 5 indicate that the SWATCH approach is clearly superior to SVF because the time for immunological failure may occur 10 years later compared to the SVF strategy. However, SWATCH is more sensitive to the initiation time. For instance, when years, SWATCH and SVF may show similar performance. We conclude that SWATCH treatment may provide significant extension of the time to virological failure only if the treatment is initiated before the third year postinfection.

While the results above were obtained for one particular parameter set, we now consider the nonlinearity of the problem and the corresponding sensitivity to parameter variations by using an MC approach with 1000 random samples (Table 3). In comparison to the MC simulations without treatment, we note in Table 3 that the number of cases without immunological failure () increases substantially (approximately 30%). For instance, when , MC simulations suggest that SVF may achieve 29.8% cases without immunological failure while SWATCH may achieve up to 46.6%. Furthermore, the appearance of virological failure is prolonged approximately by 4–7 years ( value ).

Table 3: Monte Carlo simulation outcome for SVF and SWATCH strategies.

(a) SVF

(years)
(years)

0.5
14.01 ± 4.34
298

1
13.93 ± 4.39
292

1.5
13.82 ± 4.43
284

2
13.61 ± 4.38
277

3*
13.06 ± 4.08
259

4*
12.83 ± 4.55
253

(b) SWATCH

(years)
(years)

0.5
15.54 ± 4.79
466

1
15.45 ± 4.81
454

1.5
15.13 ± 4.79
442

2*
14.69 ± 4.68
427

3*
13.67 ± 4.38
386

4*
12.59 ± 3.84
336

is the average year when immunological failure occurs for the 1000 random samples and the respective standard deviation . is the number of cases that do not experience immunological failure after 30 years postinfection. *Star represents a statistically significant difference with respect to the group of years (). There is a significant difference between SWATCH and SVF for all the groups when .

Our studies support the hypothesis that an early intervention has significant positive impact on postponing the progression to AIDS. Figure 6 shows the average time when immunological failure occurs () as a function of the initial time (). We can see that SVF performance decreases almost linearly with respect to the initial time, while SWATCH shows approximately a parabolic behaviour. Note that the SWATCH strategy can outperform the SVF strategy only when therapy is initiated before the second year postinfection ( value ). Therefore, it is not recommended to use the SWATCH strategy on patients who have not initiated treatment during 3 years postinfection.

Figure 6: Monte Carlo simulation outcome for SVF (+) and SWATCH (□) strategies. Averaged time of immunological failure () based on 1000 random samples is plotted against different initiation times (). A parabolic function was fitted to the SWATCH data (solid black line) and a linear function to the SVF data (solid blue line).

Our results lead us to recommend avoiding treatment initiation after 2 years postinfection. In addition, for patients who have not received any treatment within 3 years or more postinfection, the SVF strategy is a better alternative since the advantage of SWATCH is fading (Figure 6) while the risk of long-term drug toxication could be smaller with the SVF strategy.

3.2. Alternating between Treatments

Computational studies [7, 8, 10] and clinical trials [9] suggest that a proactive alternating strategy like SWATCH may yield promising results. However, not enough work has been done to analyse the synergistic effect of initiation time and period of alternation between treatments.

Here, we ran model (1)–(5) with different switching and different initiation times. The first observation is that fast switching times between regimens yield longer periods without immunological failure (Figure 7). In addition, for switching times larger than 150 days, simulation trajectories reached a steady state, meaning that long periods between treatment switches should be avoided. An aspect that is not covered here, however, is that fast switching between regimens (<60 days) could increase drug toxicity and may lead to bad adherence to therapy. Further studies (pharmacokinetics/pharmacodynamics) are needed to include such issues to develop a comprehensive management concept.

Figure 7: Optimal switching between regimens. Time of immunological failure () with respect to switching time for different initiation times ().

A second important conclusion from Figure 7 is that the SWATCH strategy seems to perform poorly if the treatment was initiated 2 years postinfection. To attain a high effectiveness of the treatment, the switching frequency between antiretroviral regimens may need to be increased depending on the delay of therapy initiation. Considering or years, our results suggest that very short switching times about 10–60 days would be required to obtain the same performance as if the therapy was started earlier. Consistent with the results in Figure 6, this investigation leads to the conclusion that the SWATCH strategy should be initiated within 2 years postinfection to achieve responses superior to common clinical approaches (SVF).

Our results reveal the importance of early therapy to delay the AIDS progression. There are hints that an early intervention could improve the patient healing process. The recent clinical results by [33] showed the possibility of a mechanistic cure for patients who are treated immediately after the infection. A mechanistic cure means a permanent viral suppression in the absence of therapy to levels that prevent immunodeficiency and transmission. However, it was not possible to simulate this behaviour using model (1)–(5). This is likely because the model is based on ODEs, implying that populations arbitrarily close to zero can recover. Discrete approaches could bring new insights into this field.

4. Conclusions

From our studies, we conclude that antiretroviral treatment strategies initiated after 2 years postinfection are not beneficial to extend the time to progression to AIDS.

Another significant result is that the SWATCH strategy outperforms SVF only when therapies are initiated within 2 years postinfection and switching periods for SWATCH strategy are less than 90 days. Large switching periods between regimens (>100 days) should be avoided during the application of SWATCH.

This work is a step forward for defining criteria for when to initiate and alternate therapy. Future work will be directed to the experimental evaluation of the presented results.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by BMBF-GerontoSys initiative (GerontoShield), portfolio HGF (metabolic dysfunction), and BMBF e-Med project SYSIMIT.

References

UNAIDS, Global Report: UNAIDS Report on the Global AIDS Epidemic 2010, 2010.

F. Clavel and A. J. Hance, “HIV Drug Resistance,” The New England Journal of Medicine, vol. 350, no. 10, pp. 1023–1035, 2004.

DHHS, Panel on Antiretroviral Guidelines for Adults and AdoLescents. Guidelines for the Use of Antiretroviral Agents in HIV-1-Infected Adults and Adolescents, Department of Health and Human Services (DHHS), 2013.

A. M. Jeffrey, X. Xia, and I. K. Craig, “When to initiate HIV therapy: a control theoretic approach,” IEEE Transactions on Biomedical Engineering, vol. 50, no. 11, pp. 1213–1220, 2003.

F. Castiglione and P. Paci, “Criticality of timing for anti-HIV therapy initiation,” PLoS ONE, vol. 5, no. 12, Article ID e15294, 2010.

P. Paci, F. Martini, M. Bernaschi, G. D'Offizi, and F. Castiglione, “Timely HAART initiation may pave the way for a better viral control,” BMC Infectious Diseases, vol. 11, article 56, 2011.

E. A. Hernandez-Vargas, P. Colaneri, and R. H. Middleton, “Optimal therapy scheduling for a simplified HIV infection model,” Automatica, vol. 49, pp. 2874–2880, 2013.

E. A. Hernandez-Vargas, P. Colaneri, and R. H. Middleton, “Switching strategies to mitigate HIV mutation,” IEEE Transactions on Control Systems Technology, 2014.

J. Martinez-Picado, E. Negredo, L. Ruiz et al., “Alternation of antiretroviral drug regimens for HIV infection: a randomized, controlled trial,” Annals of Internal Medicine, vol. 139, no. 2, pp. 81–89, 2003.

E. Hernandez-Vargas, P. Colaneri, R. Middleton, and F. Blanchini, “Discrete-time control for switched positive systems with application to mitigating viral escape,” International Journal of Robust and Nonlinear Control, vol. 21, no. 10, pp. 1093–1111, 2011.

M. A. Nowak, R. M. May, and R. M. Anderson, “The evolutionary dynamics of HIV-1 quasispecies and the development of immunodeficiency disease,” AIDS, vol. 4, no. 11, pp. 1095–1103, 1990.

B. Adams, H. Banks, H. Kwon, and H. Tran, “Dynamic multidrug therapies for HIV: optimal and STI control approaches,” Mathematical Biosciences and Engineering, vol. 1, pp. 223–241, 2004.

A. S. Perelson, D. E. Kirschner, and R. de Boer, “Dynamics of HIV infection of CD4+ T cells,” Mathematical Biosciences, vol. 114, no. 1, pp. 81–125, 1993.

S. H. Bajaria, G. Webb, M. Cloyd, and D. Kirschner, “Dynamics of naive and memory CD4+ T lymphocytes in HIV-1 disease progression,” Journal of Acquired Immune Deficiency Syndromes, vol. 30, no. 1, pp. 41–58, 2002.

D. E. Kirschner, R. Mehr, and A. S. Perelson, “Role of the thymus in pediatric HIV-1 infection,” Journal of Acquired Immune Deficiency Syndromes and Human Retrovirology, vol. 18, no. 2, pp. 95–109, 1998.

P. Ye, D. E. Kirschner, and A. P. Kourtis, “The thymus during HIV disease: role in pathogenesis and in immune recovery,” Current HIV Research, vol. 2, no. 2, pp. 177–183, 2004.

D. Kirschner, G. F. Webb, and M. Cloyd, “Model of HIV-1 disease progression based on virus-induced lymph node homing and homing-induced apoptosis of CD4+ lymphocytes,” Journal of Acquired Immune Deficiency Syndromes, vol. 24, no. 4, pp. 352–362, 2000.

L. Rong and A. S. Perelson, “Modeling HIV persistence, the latent reservoir, and viral blips,” Journal of Theoretical Biology, vol. 260, no. 2, pp. 308–331, 2009.

Z. Shu, E. A. Hernandez-Vargas, and R. H. Middleton, “A math ematical study on immune activation and related dynamics in HIV infection,” in Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC '11), vol. 2, pp. 3682–3687, 2011.

I. B. Hogue, S. H. Bajaria, B. A. Fallert, S. Qin, T. A. Reinhart, and D. E. Kirschner, “The dual role of dendritic cells in the immune response to human immunodeficiency virus type 1 infection,” Journal of General Virology, vol. 89, no. 9, pp. 2228–2239, 2008.

A. S. Perelson, A. U. Neumann, M. Markowitz, J. M. Leonard, and D. D. Ho, “HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time,” Science, vol. 271, no. 5255, pp. 1582–1586, 1996.

D. S. Callaway and A. S. Perelson, “HIV-1 infection and low steady state viral loads,” Bulletin of Mathematical Biology, vol. 64, no. 1, pp. 29–64, 2002.

A. Yates, J. Stark, N. Klein, R. Antia, and R. Callard, “Understanding the slow depletion of memory CD4+ T cells in HIV infection,” PLoS Medicine, vol. 4, no. 5, Article ID e177, 2007.

N. Dalal, D. Greenhalgh, and X. Mao, “A stochastic model for internal HIV dynamics,” Journal of Mathematical Analysis and Applications, vol. 341, no. 2, pp. 1084–1101, 2008.

N. I. Stilianakis, K. Dietz, and D. Schenzle, “Analysis of a model for the pathogenesis of AIDS,” Mathematical Biosciences, vol. 145, no. 1, pp. 27–46, 1997.

E. A. Hernandez-Vargas and R. H. Middleton, “Modeling the three stages in HIV infection,” Journal of Theoretical Biology, vol. 320, pp. 33–40, 2013.

J. M. Orenstein, “The macrophage in HIV infection,” Immunobiology, vol. 204, no. 5, pp. 598–602, 2001.

D. Kirschner and A. Perelson, “A model for the immune system response to HIV: AZT treatment studies,” in Mathematical Population Dynamics, pp. 295–310, Wuerz, Winnepeg, Canada, 1995.

M. Hadjiandreou, R. Conejeros, and V. S. Vassiliadis, “Towards a long-term model construction for the dynamic simulation of HIV infection,” Mathematical Biosciences and Engineering, vol. 4, no. 3, pp. 489–504, 2007.

C. J. Mode, Applications of Monte Carlo Methods in Biology, Medicine and Other Fields of Science, InTech, 2011.

A. S. Fauci, G. Pantaleo, S. Stanley, and D. Weissman, “Immunopathogenic mechanisms of HIV infection,” Annals of Internal Medicine, vol. 124, no. 7, pp. 654–663, 1996.

S. Sungkanuparph, R. K. Groger, E. T. Overton, V. J. Fraser, and W. G. Powderly, “Persistent low-level viraemia and virological failure in HIV-1-infected patients treated with highly active antiretroviral therapy,” HIV Medicine, vol. 7, no. 7, pp. 437–441, 2006.

A. Sáez-Cirión, C. Bacchus, L. Hocqueloux et al., “Post-treatment HIV-1 controllers with a long-term Virological remission after the interruption of early initiated antiretroviral therapy ANRS VISCONTI study,” PLoS Pathogens, vol. 9, Article ID e1003211, 2013.