- Research article
- Open Access
- Open Peer Review
- Published:

# Low-level HIV-1 replication and the dynamics of the resting CD4^{+}T cell reservoir for HIV-1 in the setting of HAART

*BMC Infectious Diseases***volume 8**, Article number: 2 (2008)

## Abstract

### Background

In the setting of highly active antiretroviral therapy (HAART), plasma levels of human immunodeficiency type-1 (HIV-1) rapidly decay to below the limit of detection of standard clinical assays. However, reactivation of remaining latently infected memory CD4^{+} T cells is a source of continued virus production, forcing patients to remain on HAART despite clinically undetectable viral loads. Unfortunately, the latent reservoir decays slowly, with a half-life of up to 44 months, making it the major known obstacle to the eradication of HIV-1 infection. However, the mechanism underlying the long half-life of the latent reservoir is unknown. The most likely potential mechanisms are low-level viral replication and the intrinsic stability of latently infected cells.

### Methods

Here we use a mathematical model of T cell dynamics in the setting of HIV-1 infection to probe the decay characteristics of the latent reservoir upon initiation of HAART. We compare the behavior of this model to patient derived data in order to gain insight into the role of low-level viral replication in the setting of HAART.

### Results

By comparing the behavior of our model to patient derived data, we find that the viral dynamics observed in patients on HAART could be consistent with low-level viral replication but that this replication would not significantly affect the decay rate of the latent reservoir. Rather than low-level replication, the intrinsic stability of latently infected cells and the rate at which they are reactivated primarily determine the observed reservoir decay rate according to the predictions of our model.

### Conclusion

The intrinsic stability of the latent reservoir has important implications for efforts to eradicate HIV-1 infection and suggests that intensified HAART would not accelerate the decay of the latent reservoir.

## Background

The latent reservoir for HIV-1 in resting CD4^{+} T cells is generated when productively infected CD4^{+} T lymphoblasts revert back to the resting state, becoming memory T cells, instead of succumbing to viral cytopathic effects or host cytolytic effector mechanisms [1–4]. The result is a state of viral latency in resting memory CD4^{+} T cells, cells that are extremely quiescent, with little to no transcription of HIV-1 genes [5–7]. Given that memory T cells form the basis for lifelong immunity to recall antigens, it is not surprising that the average half-life of the latent reservoir in the setting of HAART can be as long as four years [8, 9]. However, the basis for the remarkable stability of the latent reservoir has remained controversial.

The two most reasonable mechanisms for maintenance of the latent reservoir in the setting of HAART are 1) replenishment by low-level viral replication [10–20] and 2) the intrinsic stability of latently infected cells (i.e. memory T cells) [8, 9, 21–23]. While some studies have suggested that low-level viral replication confers stability by continuously reseeding the latent reservoir despite HAART [10, 19, 20], other studies have provided experimental evidence at odds with a major role for viral replication in maintaining the latent reservoir [24, 25]. These studies have shown that in many patients responding well to HAART, there is no evolution of drug resistance, suggesting a lack of viral replication [26]. We have previously shown that the maximal rate at which new cells enter the reservoir in the setting of HAART is extremely low [27]. These studies provide indirect evidence that intrinsic stability of memory T cells and not replenishment by ongoing viral replication is the major reason for the stability of the latent reservoir.

Mathematical models have proven useful for the analysis of several aspects of HIV-1 infection including the dynamics of viral replication [28–31], the effects of immune responses [32–35], and the mechanism of CD4 depletion [32, 36–38]. We present here a mathematical analysis of CD4^{+} T cell dynamics in the setting of HIV-1 infection in order to explore the dynamics of the latent resting CD4^{+} T cell reservoir. We extend elegant models of HIV-1 and CD4^{+} T cell dynamics previously described by Alan Perelson and Martin Nowak [28, 32] to explore how low-level viral replication influences the observed decay of the latent reservoir in patients on HAART. A recent study [39] analyzed the persistence of the latent reservoir in the setting of HAART with a model similar to ours. However, this study [39] did not focus on the decay properties of latently infected cells in relation to low-level viral replication. Also, because the authors did not constrain the maximum amount of viral replication compatible with available experimental data from patients on HAART, this study [39] was unable to answer the clinically significant question of whether realistic levels of residual viral replication in the setting HAART affect the experimentally observed decay rate of the latent reservoir. In this study, we calculate the well-known replication threshold below which infection cannot be sustained [40] for our model, and discuss latent reservoir replenishment above and below this threshold. Having explicitly illustrated the primary factors involved in establishing and maintaining the latent reservoir, we offer the first explicit analysis of the relationship between low-level viral replication and the decay rate of the latent reservoir. Our results indicate that the impact of viral replication on the decay rate of the latent reservoir rapidly diminishes with increasing inhibition by HAART. For levels of viral replication likely to occur in the setting of HAART, our model predicts that the decay of the latent reservoir is determined by intrinsic (e.g. death rate) as well as extrinsic properties (e.g. activation rate) of latently infected cells. We further apply our theoretical analysis to patient-derived data and show that any flow of new cells into the latent reservoir by viral replication is unlikely to impact the observed, *in vivo* decay rate. Thus theoretical predictions and experimental data both suggest that the long half-life of the latent reservoir is likely attributable to the intrinsic stability and reactivation rate of latently infected T cells. As such, any possible shortening of the half-life of latently infected CD4^{+} T cells would depend on increasing the death rate or reactivation rate of these cells. These results have important implications for eradication of HIV-1 infection by conventional HAART and suggest the necessity of developing strategies to target the latent reservoir specifically by approaches either apart from or in addition to HAART intensification.

## Methods

### Mathematical model of T cell dynamics

To explore the possibility that low-level viral replication reseeds the latent reservoir and thereby slows the decay of the reservoir in the setting of HAART, we used a slightly modified version of previously described mathematical models of T cell dynamics in HIV-1 infection [28, 32]. By virtue of its foundation in previously described models [28, 32], our model makes the same assumptions as these models: constant model parameters (e.g. viral infectivity) as well as a well-stirred, homogeneous virus and cell populations. Our modifications to these models make relatively few additional assumptions about the physiology of CD4^{+} T cell dynamics, as described below. Our model (Figure 1) considers uninfected (activated) target cells (*U*), productively infected (activated) cells (*P*) and latently infected cells (*L*). We consider infected cells as those that contain a stably integrated HIV-1 genome. Uninfected cells are introduced into the system according to a zero order rate constant, *λ*, which represents production of target cells (e.g. by activation of resting CD4^{+} T cells as well as thymic or homeostatic production). Uninfected cells are infected by virions to become productively infected cells at a rate proportional to a constant, *β*, which represents the virus infectivity and is a reflection of viral replication. Because free virions have a very high turnover rate and are made by productively infected cells (infected activated CD4^{+} T cells), their dynamics closely resemble those of productively infected cells. Therefore free virions are not explicitly represented but rather are represented by productively infected cells. Thus the rate of infection is represented as a non-linear term (*βUP*) proportional to both the target (uninfected) cells and the productively infected cells [32]. Only productively infected cells replenish the latent reservoir [41]. We model the reversion of productively infected cells to latency at a rate proportional to *α*
_{
R
}, while latently infected cells may reactivate at a rate proportional to *α*
_{
Q
}.

Productively infected and latently infected cells each have intrinsic death rates represented by *d*
_{
P
}, and *d*
_{
L
}, respectively. Uninfected activated CD4+ T cells can be lost by conversion into resting T cells or by cell death. However, since resting T cells cannot be productively infected by HIV, they are irrelevant for our model, and thus we subsume any kind of loss of activated T cells under the rate constant *d*
_{
U
}. The model is then defined by a system of three ordinary differential equations:

We do not incorporate other permissive cell types (e.g. macrophages) since they comprise a small fraction of all infected cells [31] and would not significantly change the overall qualitative as well as quantitative findings.

Parameter values were chosen based on previously published reports, where possible. Most productively infected cells have been shown to have extremely short half-lives [30, 31], while latently infected cells have remarkable stability [8, 9]. These characteristics are reflected in the death rates of productively infected cells, *δ*
_{
P
}, and latently infected cells, *δ*
_{
L
}. The rate at which target cells (uninfected activated CD4^{+} T cells) turn over, *δ*
_{
U
}, represents a balance between death, proliferation and reversion to resting. In the case of no infection (i.e. in an HIV^{-} individual), the steady-state number of target cells, $\overline{U}$ is *λ*/*δ*
_{
U
}or *λ* = *δ*
_{
U
}
*U*. Therefore *δ*
_{
U
}also represents fraction of target cells that are replaced daily. Previously published studies estimate that at most 2% of target cells are replaced daily in an uninfected individual, therefore we set *δ*
_{
U
}= 0.02 day^{-1} (based on direct measurements in [42, 43]). Furthermore, we find that for *δ*
_{
U
}= 0.02 day^{-1}, our model reproduces the clinically observed, smooth rebound [29, 44, 45] of the viral load (reflected here by the number of productively infected cells) after interruption of HAART with little to no oscillation (data not shown). Because activated cells both proliferate rapidly and die rapidly, it is possible that the net death rate (death rate less proliferation rate) could be very small and on the order of 0.001 day^{-1}. However, we do not expect that *δ*
_{
U
}could be much larger than our chosen value of 0.02 day^{-1}. A larger value for *δ*
_{
U
}would be contradictory to our expectation that *δ*
_{
U
}<<*δ*
_{
P
}, which follows from the highly cytopathic nature of HIV. Regardless, while we choose *δ*
_{
U
}to best fit experimentally observed data, we nonetheless test values of *δ*
_{
U
}as low as 0.001 day^{-1} and as large as 1.0 day^{-1} and find our results to be robust through this entire range.

Infection of an activated CD4^{+} T cell by HIV leads to significant dysregulation of cellular processes as the virus subverts the host cellular machinery. This dysfunction may also lead to decreased conversion to the resting state, which is reflected by the substantially lower (by several orders of magnitude) number of latently infected cells compared to productively infected cells [3]. We therefore chose a rate for conversion of productively infected cells to latency, *α*
_{
R
}, to reflect this experimentally observed balance. For the numerical calculations and simulations performed in this report, we used the following parameter choices (unless specified otherwise): *δ*
_{
P
}= 0.5 day^{-1}, *δ*
_{
L
}= 0.0001 day^{-1}, *α*
_{
Q
}= 0.0005 day^{-1}, *α*
_{
R
}= *α*
_{
Q
}/100, and *δ*
_{
U
}= 0.02 day^{-1} with *λ* constrained to 2 × 10^{9} cells/day by the choice of *δ*
_{
U
}and the steady-state value of uninfected cells in the setting of no infection, which is approximately 10^{11} cells (Table 1). However, we tested a wide range of parameter values in order to check the validity of our results (data not shown).

All simulations and calculations were performed with MATLAB version 7.2.0.232. **Patient Data.** Patients' viral load records were obtained with informed consent.

## Results

### A threshold of viral replication is necessary for existence of the latent reservoir

Analysis of our model (equations 1–3 in *Methods*) leads to several predictions about the nature of HIV-1 infection and its impact on CD4^{+} T cell dynamics. The model has two steady-state solutions: the trivial steady-state reflects conditions under which infection does not successfully occur and the system remains at the uninfected steady-state $(\overline{U}=\frac{\lambda}{{\delta}_{U}},\overline{P}=\overline{L}=0)$. The non-trivial steady-state values of all three state variables $(\overline{U},\overline{P},\overline{L})$ are given by the following equations:

where *A* = *δ*
_{
L
}
*δ*
_{
P
}+ *δ*
_{
P
}
*α*
_{
Q
}+ *δ*
_{
L
}
*α*
_{
R
}and *B* = *δ*
_{
L
}+ *α*
_{
Q
}. These solutions indicate that the steady-state number of latently infected cells depends on *β*, the infectivity, as does the steady state number of productively infected cells. Inspection of equation 6 immediately reveals that the non-trivial steady-state solutions (equations 4–6) are not biological (*P*, *L* ≤ 0) for *β* ≤ *β*
_{
crit
}, where

Because *δ*
_{
P
}>> *δ*
_{
L
}, *α*
_{
Q
}>> *α*
_{
R
}, we can approximate *A* ≈ *δ*
_{
P
}(*δ*
_{
L
}+ *α*
_{
Q
}) = *δ*
_{
P
}
*B* and therefore

A successful infection is not established (the trivial steady-state occurs) for *β* ≤ *β*
_{
crit
}. A critical threshold dependence of the steady-state for infected cells on the infectivity was previously described for models of HIV viral and T cell dynamics [40]. Our model demonstrates a similar property and so from the condition that *L* > 0, it follows that *β* > *β*
_{
crit
}for a successful infection. Stability analysis shows that the non-trivial steady-state is stable (and trivial steady-state is unstable) when *β* > *β*
_{
crit
}, and that the converse is true when *β* ≤ *β*
_{
crit
}.

Based on equation 8, it is apparent that *β*
_{
crit
}is almost entirely determined by *δ*
_{
U
}, *δ*
_{
P
}and *λ*. However, the trivial steady-state solution enforces that the number of target cells in the HIV^{-}individual is *λ*/*δ*
_{
U
}. Because the quantity *λ*/*δ*
_{
U
}is roughly constant across all uninfected individuals, the primary factor influencing *β*
_{
crit
}is *δ*
_{
P
}, the death rate of productively infected cells. We find a stable, physiologically non-trivial steady state for *β* > *β*
_{
crit
}that reflects the need for *β* to be large enough (to generate enough productively infected cells) to overcome the large death rate of productively infected cells. Likewise, the maintenance of latently infected cells requires at least a minimum of productively infected cells sufficient for regular transitions to latency despite rapid cell death. For rate constants consistent with previously published studies (Table 1), the predicted steady-state levels of both $\overline{P}$ and $\overline{L}$ change minimally when small perturbations are made to *β*
_{
Untreated
}. However, both $\overline{P}$ and $\overline{L}$ rapidly drop to zero as *β* approaches *β*
_{
crit
}. If *β* ≤ *β*
_{
crit
}, then the level of ongoing viral replication is unable to maintain infection and the infection subsequently decays.

In order to understand how the size of the latent reservoir relates to viral replication, we now examine the analytic relationship between infectivity and the ratio $\overline{L}/\overline{U}$, as obtained by dividing equation 6 by equation 4. From $\overline{L}/\overline{U}$ we solve for the viral infectivity

in an untreated patient at steady-state. Therefore viral replication cannot replenish the latent reservoir if *β* is reduced from *β*
_{
Untreated
}to *β*
_{
crit
}where

Thus even in the setting of ongoing viral replication and HAART, there should be no net entry of new cells into the latent reservoir if HAART has attenuated the viral replication such that *β* ≤ *β*
_{
crit
}(where *β*
_{
crit
}represents the minimum amount of viral replication necessary to maintain a latent reservoir). Below, we will argue that the experimental data constrain *β* to below *β*
_{
crit
}in the setting of HAART.

### Low-level viral replication does not significantly affect the decay rate of the latent reservoir

By preventing new entry into the latent reservoir, HAART induces the latent reservoir to decay at a rate that depends on the degree of inhibition of viral replication. With our model, T cell dynamics in the setting of HAART can be simulated by decreasing the infectivity, *β*, from *β*
_{
Untreated
}to *rβ*
_{
Untreated
}, where residual viral replication is represented by *r* = *β*
_{
HAART
}/*β*
_{
Untreated
}, which reflects the fraction of *β*
_{
Untreated
}to which viral infectivity is reduced in the setting of HAART. If we assume that *U* changes slowly relative to *P* and *L* (in the setting of HAART, the number of uninfected cells changes very little) [28], we can substitute $\overline{U}$ for *U* in equations 1 and 2 – obtaining a system of 2 linear, ordinary differential equations (for *P* and *L*). This system is characterized by two eigenvalues, where

is the asymptotic (least negative) eigenvalue and reflects the decay rate of the latent reservoir. The maximum value of Λ, Λ_{
r = 0}, is the fastest decay rate of the latent reservoir and occurs when *r* = 0, reflecting a complete inhibition of viral replication. By expanding equation 11 to zeroeth order in *α*
_{
R
}for *r* = 0, the maximum reservoir decay rate may be approximated as Λ_{
r= 0 }≈ *B* = *δ*
_{
L
}+ *α*
_{
Q
}. The decay rate of the latent reservoir is therefore determined primarily by the activation rate of the latently infected cells (*α*
_{
Q
}) as well as the net death rate of latently infected cells (*δ*
_{
L
}), which reflects the difference between the death rate and proliferation rate of latently infected cells. Thus the faster latently infected cells die, the faster they decay (larger *δ*
_{
L
}). Conversely, the faster latently infected cells divide (e.g. homeostatic proliferation), the slower they decay (smaller *δ*
_{
L
}). In principle, however, because resting CD4^{+} T cells both divide and die very slowly, it is likely that *α*
_{
Q
}> *δ*
_{
L
}. Therefore, the reactivation rate of latently infected cells makes the dominant contribution to Λ, the decay rate of latently infected CD4^{+} T cells.

Further, we observe that Λ increases as *A*(1 - *r*)/(*α*
_{
R
}+ *B*) until it approximately reaches its maximum (Λ_{
r = 0}) when *r* = *r** ≈ (*δ*
_{
P
}- *δ*
_{
L
}- *α*
_{
R
}- *α*
_{
Q
})/*δ*
_{
P
}, at which point the reservoir decay rate is virtually independent of *β* (Figure 2). This relationship between Λ and 1 - *r*, a consequence of the nonlinear structure of the model, reflects the threshold of viral replication. Once viral replication has been reduced to a level so low that the virus population cannot sustain itself anymore, the replenishment of the latent CD4^{+} T cell reservoir through ongoing replication becomes negligible and has virtually no effect on the reservoir decay rate. This relationship between Λ and 1 - *r* also holds for other parameter combinations and the predicted reservoir decay rates closely match simulation results. (Note that for *β* > *β*
_{
crit
}, Λ describes the decay of the reservoir towards the new steady state, not towards zero.) In deriving an expression for the decay rate, we assume that the model parameters remain constant from when HAART is started. However, we also consider the consequences of linearizing around the post-HAART value of $\overline{U}$. At long time scales, the system might approach the post-HAART steady-state; therefore we must consider this possibility. For the case of *β* ≤ *β*
_{
crit
}, we find that linearization around the post-HAART value of $\overline{U}$ does not significantly change our previous calculations – the reservoir continues to decay at a constant rate determined by equation 11 except that *r* = *β*
_{
HAART
}/*β*
_{
crit
}, which does not significantly change the numerical value of Λ. For the case of *β* > *β*
_{
crit
}, linearization around the post-HAART value of $\overline{U}$ leads to Λ ≈ 0 (because *r* → 1 in equation 11), suggesting a very slow (almost infinitely slow) decay of the latent reservoir near the post-HAART steady-state. These calculations predict that for the case of *β* > *β*
_{
crit
}, the reservoir decay rate (as defined by the asymptotic eigenvalue of the equations defining our model) slows down as the post-HAART steady-state is achieved. Consistent with results that we have described in a previous study (regimes 1–3 of reference [27]), the decay of the resting CD4^{+} T cell reservoir cannot be accurately described as an exponential decay near the post-HAART steady-state in the setting of viral replication sufficient for continual replenishment of the reservoir (*β* > *β*
_{
crit
}).

### Comparison with previous models of latent reservoir persistence

The decay of the latent reservoir was previously modeled by Muller *et al.* [46], who made the simplification that HAART eliminates the flow of cells into latent reservoir from the productively infected cell compartment. This simplification reduced the dynamics of the latent reservoir in the setting of HAART to a strict exponential decay dependent on only the death rate and reactivation rate of latently infected cells, and allowed Muller *et al.* to study how the reservoir decays when the cells of the latent reservoir have a distribution of reactivation rates rather than one constant reactivation rate. While quite insightful, the study by Muller *et al.* does not address the effect of viral replication on the decay of the latent reservoir by virtue of their simplifying assumption that viral replication does replenish the latent reservoir.

The persistence of low-level viremia and the latent HIV-1 reservoir in the setting of HAART was also addressed by Kim and Perelson [39], who use a simplified variation of a previous model of viral dynamics [30, 31]. Kim and Perelson perform a thorough analysis of their model's parameter space to understand how various model parameters contribute to persistence of low-level viremia and a latent reservoir. Considering the effects of viral replication, Kim and Perelson find different regimes of behavior for their model (with respect to latently infected cells and viral load) that depend upon the degree of HAART efficacy (*ε*
_{
HAART
}) relative to a critical drug efficacy (*ε*
_{
crit
}) (where Kim and Perelson's efficacy, *ε*, corresponds to 1 - *r* in our model so that a larger efficacy indicates greater inhibition of viral replication). Kim and Perelson find that only when *ε*
_{
HAART
}> *ε*
_{
crit
}do both the latent reservoir and viral load decay towards zero. This finding is robust for all physiologic ranges of their model parameter values.

While offering insight into the persistence of low-level viremia and the latent reservoir, the analysis presented in Kim and Perelson does not directly address two key clinical issues: (1) whether *ε*
_{
HAART
}> *ε*
_{
crit
}in the typical HAART-treated patient and (2) how the decay rate of the latent reservoir is affected by further increasing the efficacy of HAART (e.g. by intensification of the HAART regimen) when *ε*
_{
HAART
}is already greater than *ε*
_{
crit
}. Because it has been previously suggested that intensification of HAART can increase the decay rate of the latent reservoir [20], we repeated Kim and Perelson's analysis in order to determine whether increasing *ε*
_{
HAART
}despite the fact that *ε*
_{
HAART
}> *ε*
_{
crit
}would predict an increase in the decay rate of the latent reservoir [20]. We find that the decay dynamics of the latent reservoir remain approximately the same regardless of *ε*
_{
HAART
}, as long as *ε*
_{
HAART
}> *ε*
_{
crit
}. This finding is consistent with their steady-state analysis, which found that the eigenvalue describing the long-term decay characteristics of the latent reservoir is independent of *ε* (i.e. it is unaffected by ongoing viral replication). This eigenvalue, solved for by Kim and Perelson, is consistent with the simplification of reservoir dynamics made by Muller *et al.* [46] as well as equal to our derived approximation of -Λ_{
r= 0 }(≈ -*B* = -(*δ*
_{
L
}+ *α*
_{
Q
})) and is therefore dependent on only the net death and reactivation rate of latently infected cells. Kim and Perelson's larger model (defined by equations 1–7 of [39]) also demonstrates a threshold effect for the decay of the latent reservoir, whereby increasing HAART efficacy beyond a certain point does not accelerate the decay rate of the latent reservoir. Whether the average HAART-treated patient has reached this point of viral inhibition remains to be addressed.

### Maximal levels of viral replication consistent with patient data do not affect the decay rate of the latent reservoir

In the previous sections we have shown that a threshold of viral replication is necessary for existence of the latent reservoir and that the decay rate of the reservoir is rapidly maximized with decreasing replication. To constrain viral replication in patients on HAART, we compare the behavior of the model with data from patients on HAART. If HAART reduces *β* (equal to *β*
_{
untreated
}) by a factor *r* = *β*
_{
HAART
}/*β*
_{
untreated
}<*β*
_{
crit
}/*β*
_{
untreated
}, our model predicts that the number of productively infected cells (and by extension, viral load) decays to zero. We find that when r is sufficiently less than *β*
_{
crit
}/*β*
_{
untreated
}(in our system *r* ≲ 0.95 *β*
_{
crit
}/*β*
_{
untreated
}), the behavior of the decay of productively infected cells can be described as undergoing a rapid and approximately exponential decay, reflecting the elimination of productively infected cells produced by new infections. This initial rapid decay is followed by a slower, also approximately exponential decay, reflecting reactivation of the decaying latently infected cell pool (Figure 3A), towards zero. This result is consistent with clinically observed decay of viral load in HIV^{+} patients as they start HAART with the exception that our model does not predict a second phase decay. The latter is absent because we include only CD4^{+} T cell compartments in our model. This decay behavior of the productively infected cell pool demonstrated by our model occurs regardless of the parameter choices tested, as long as *r* <*β*
_{
crit
}/*β*
_{
untreated
}. If HAART reduces *β* such that *r* ≳ *β*
_{
crit
}/*β*
_{
untreated
}(indicating a suboptimal HAART regimen that cannot fully suppress the infection), then we must consider two cases. In our model, we find two regimes of behavior for the decay of productively infected cells that depend on the value of *δ*
_{
U
}(the rate at which target cells are removed from the system).

The dependence of the system's qualitative behavior on *δ*
_{
U
}is robust to choice of other model parameters. With the other model parameter values that we use, if *δ*
_{
U
}≲ 0.1035 day^{-1} then suboptimal HAART causes the productively infected population to rapidly decay but subsequently rebound, undergoing a dampened oscillation (see Appendix) towards a final steady-state that is greater than zero if *r* > *β*
_{
crit
}/*β*
_{
untreated
}or a zero steady-state if *r* = *β*
_{
crit
}/*β*
_{
untreated
}(Figure 3B). In contrast, if *δ*
_{
U
}≳ 0.1035 day^{-1} then in the setting of a suboptimal HAART regime productively infected cells undergo a smooth but slow sub-exponential decay (see Appendix) toward a greater-than-zero final steady state if *r* > *β*
_{
crit
}/*β*
_{
untreated
}or a zero steady-state if *r* = *β*
_{
crit
}/*β*
_{
untreated
}(Figure 3C).

Previous longitudinal studies of viral load in HIV^{+} patients starting HAART have observed a characteristic decay of the viral load [4, 18, 30, 31] consistent with the patients we have observed (Figure 3D). In these patients, who are typical of patients responding well to HAART, the viral load is observed to make a rapid, approximately exponential decay, which reflects decay of productively infected cells. In these patients, who are typical of patients responding well to HAART, the viral load is observed to make a rapid, approximately exponential decay, which reflects decay of productively infected cells. The viral load decays to below the limit of detection, reaching a quasi-steady-state that has recently been observed experimentally [47] and is consistent with a much slower phase of decay (half-life of >67 weeks). These observed dynamics of low-level viremia most likely represents reactivation from the decay of long-lived reservoirs for HIV-1 (e.g. resting CD4^{+} T cells). Clinically, the viral load of a patient who responds well to HAART (i.e. a successful and optimal HAART regimen) never makes a sustained rebound after the initial, rapid drop. In our simulations (Figure 3A), we observe this characteristic decay of productively infected cells only for values of *β* up to *β* ≈ *β*
_{
crit
}. Our model therefore conservatively constrains *β*
_{
HAART
}≲ *β*
_{
crit
}based on the available experimental data, but in fact *β*
_{
HAART
}may be as low as zero.

For almost all values of *r* that we tested, the decay of latently infected cells, *L*, towards a new steady-state (that depends on the value of *r*) or eradication was the same (Figure 3E). For the case where *r* > *β*
_{
crit
}/*β*
_{
untreated
}, the decay of the reservoir can be described as approximately exponential only initially. However, at some point when the reservoir size approaches the non-zero post-HAART steady-state, the reservoir decay becomes sub-exponential and can no longer be described with the asymptotic eigenvalue. This is apparent by the curves in Figure 3E for *r* > *β*
_{
crit
}/*β*
_{
untreated
}, which initially decay along the exponential decay curves where *r* ≤ *β*
_{
crit
}/*β*
_{
untreated
}but subsequently separate as sub-exponential decays towards the non-zero post-HAART steady-state.

Based on the constraint that $r=\frac{{\beta}_{HAART}}{{\beta}_{Untreated}}\lesssim \frac{{\beta}_{crit}}{{\beta}_{Untreated}}$ predicted by our model, *r* may range between 0.1428 to 0.9940 for values of *δ*
_{
U
}ranging from 0.001 day^{-1} to 1.0 day^{-1}. For these values of *r*, we calculate that the decay of the latent reservoir observed in patients who are successfully treated with HAART is essentially at its maximum value Λ_{
r= 0 }(which occurs when there is no viral replication) and is therefore independent of any ongoing, low-level viral replication.

For these calculations, we used parameter values consistent with previously reported values (Table 1) and a 44 month half-life for the latent reservoir [8, 9]. Previous studies have estimated the fraction of resting CD4^{+} T cells that consists of latently infected cells to be on the order of 1 in 10^{4} [2–4]. Although *U* represents target cells (i.e. activated CD4^{+} T cells) in our model, we estimate $\overline{L}/\overline{U}$ = 10^{-4} because the number of activated CD4^{+} T cells (e.g. CD69^{+}, CD25^{+} and/or HLA-DR^{+}) is on the same order of magnitude as the number of resting CD4^{+} T cells in patients on HAART [48, 49]. Alternative parameter choices also consistent with previously reported values, but for which the latent reservoir decays with a half-life not equal to 44 months (e.g. 6 months, 12 months) or alternative values of $\overline{L}/\overline{U}$ (e.g. 10^{-5}–10^{-6}) show that Λ is nonetheless very close to Λ_{
r = 0}.

### Implications for eradication of the latent reservoir

For values of *β* that are consistent with clinically observed viral load decay characteristics in the setting of HAART, the reservoir decay rate remains essentially constant at Λ_{
r= 0 }Because Λ_{
r= 0 }≈ *B* = *α*
_{
Q
}+ *δ*
_{
L
}, the major determinants of the reservoir decay rate are the activation rate and the intrinsic death rate of latently infected cells. We used our model to compare the potential for increasing the decay rate of the latent reservoir by enhanced inhibition of viral infectivity, which would occur in the setting of intensified HAART, to enhanced reactivation of latently infected cells, a therapeutic option that is currently being developed [50–52]. Figure 4A describes the case of viral inhibition by HAART. For non-suppressive reductions in the infectivity (*r* > *β*
_{
crit
}/*β*), the decay dynamics of the latent reservoir diverge from an exponential decay and asymptotically move toward a non-zero steady-state number of latently infected cells. However, for a suppressive reduction in infectivity (*r* <*β*
_{
crit
}/*β*), the latent reservoir decays exponentially with a rate constant equal to *B* = *α*
_{
Q
}+ *δ*
_{
L
}towards zero latently infected cells (Figure 4A). Our model therefore predicts that any drug or drug combination reducing HIV infectivity by a factor *r* <*β*
_{
crit
}/*β* will cause the latent reservoir to decay with the same dynamics. We do not know the value of *r* for current HAART regimens, although our analysis suggests that *r* <*β*
_{
crit
}/*β*. Nonetheless, we must consider the possibility that *r* > 0, in which case intensified HAART may further reduce *r* without hastening the decay of the latent reservoir. It is also possible that standard HAART regimens reduce infectivity such that *r* = 0 (i.e. *β*
_{
HAART
}= 0), in which case intensification of standard HAART regimens would have no conceivable benefit.

Unlike the case of decreasing infectivity, we find that increasing reactivation of latently infected cells will greatly increase the decay of the reservoir (Figure 4B). This result suggests that while the full potential of HAART on the reservoir decay rate has been achieved, enhanced reactivation of latently infected cell may offer another approach to accelerating the decay rate of the latent reservoir.

## Discussion

Despite rapid decay of free virus and of most productively infected cells after the initiation of HAART, the latent reservoir for HIV-1 in resting CD4^{+} T cells remains as an inducible source of viremia that can potentially lead to rebound in plasma virus levels upon discontinuation of HAART. The latent reservoir is extremely stable, with a half-life of 44 months [8, 9, 53]. Because the latent reservoir represents a major barrier to eradication, many studies have focused on developing methods to purge the reservoir [50–52]. Successful eradication of the latent reservoir, however, depends on understanding the mechanism underlying its remarkable stability.

The two factors believed to be most important in determining the decay rate of the latent reservoir are 1) the degree of replenishment of the reservoir by ongoing viral replication that continues despite HAART and 2) the intrinsic stability of resting memory T cells, the cells that harbor latent HIV-1. How each factor contributes to stabilizing the reservoir has important implications in development of methods to eradicate the reservoir. Evidence for replenishment of the reservoir in the setting of HAART comes from clinical studies showing accelerated reservoir decay for patients on an intensified HAART regimen over that observed in patients on standard HAART regimen [10, 19, 20]. The results of these studies suggest that standard HAART does not stop all replication, raising the possibility that ongoing replication could play a role in maintaining the latent reservoir. However, there are several lines of evidence that argue against ongoing viral replication in patients responding well to HAART. First, a meta-analysis of a large number of clinical trials of antiretroviral drugs has shown that the most successful HAART regimens are those that include two drugs with an extremely low genetic barrier to resistance (efavirenz and lamivudine or emtricitabine) (John A. Bartlett, *12*
^{th}
*Conference on Retroviruses and Opportunistic Infections*, 2005). Second, patients who demonstrate consistent suppression of viremia on HAART do not develop drug-resistant virus [9, 25, 26, 54]. Third, direct examination of the residual viremia has shown a static evolutionary pattern, with no evidence of evolution even during blips [24, 26]. Fourth, in many patients on HAART, most of the residual viremia is comprised of a small number of viral clones that are released into the plasma for long periods of time without evolution, suggesting a mechanism of persistence other than viral replication [55]. Using patients who demonstrate these predominant plasma clones, we have quantitatively constrained the rate at which new cells are introduced in the latent reservoir [27]. This study confirmed the minimal contribution of new entrants into the latent reservoirs of these patients. The lack of evidence for replication-driven replenishment of the latent reservoir suggests that intrinsic stability of the reservoir is the most significant factor in determining the observed reservoir decay rate [8, 9, 21–23].

There is no direct evidence, however, against the ability of low-level replication in the setting of HAART to stabilize the latent reservoir. In this study, we extend a previously described mathematical model of HIV dynamics by incorporating latently infected cells in order to understand the effect of low-level viral replication on the decay of the latent reservoir. Previously, more complicated models of T cell dynamics in the setting of HIV infection have been published [33, 56]. However, these models often incorporate physiologic processes that have not yet been well characterized in the experimental literature. Because the kinetics of these processes have not been quantified experimentally, it is difficult to infer much information from these models. By contrast, we use a more concise model of T cell dynamics that incorporates widely accepted physiologic relationships between different T cell populations. While all aspects of T cell dynamics are not yet known, the kinetics of the T cell activation processes included in our model represent the best characterized in the experimental literature. Based on the previous track record of similar models and the preponderance of data supporting the structure as well as results of our model, we believe that our model provides insight into the effects of low-level viral replication in the setting of HAART.

Like previous models, in our model a threshold of viral replication must be surpassed to overcome the cytopathic effects of HIV-1 infection [40]. Our model predicts two regimes for how ongoing viral replication impacts the decay rate of the latent reservoir. In the first regime, high levels of ongoing viral replication slow the decay rate of the latent reservoir; in this regime, attenuation of viral replication would hasten the decay of the latent reservoir. In the second regime, lower levels of viral replication do not significantly affect the decay rate of the latent reservoir. The low frequency at which a productively infected cell transitions to latency requires that many productively infected cells be present before one latently infected cell is produced. In this regime, the reduced number of productively infected cells is insufficient to consistently replenish the latent reservoir and therefore further attenuation of viral replication would not increase the decay rate of the latent reservoir. Based on previously reported model parameter values, we numerically explored the boundaries of these regimes by calculating the predicted reservoir decay rate for various levels of viral replication in the setting of HAART. We found that the impact of viral replication on the reservoir decay rate is rapidly minimized with decreasing viral replication. It has been hypothesized that there exist other cellular reservoirs for HIV-1, most likely in the form of monocytes and tissue macrophages [57, 58], that are possibly maintained in the setting of HAART by ongoing viral replication. While these reservoirs may potentially be purged faster with intensification of HAART, our study shows that the decay rate of the resting CD4^{+} T cell reservoir cannot be increased in this manner.

Based on our calculations of *r** (the reduction in *β* necessary before the decay of the latent reservoir becomes independent of ongoing replication; Figure 2) and our model predictions for the dynamics of viral load (reflected by productively infected cells), we would intuitively expect standard HAART to reduce viral replication into the regime where residual replication would have an insignificant effect on the reservoir decay rate. To test this theory, we compared previously reported patient data to the behavior of our model for different levels of HAART efficacy. Our analysis found that the behavior of the model can quite clearly reproduce experimental findings observed in the setting of suppressive HAART. In this regime of low-level viral replication, our model predicts that viral replication does not significantly affect the decay rate of the latent reservoir.

If a small amount of viral replication is insufficient to affect the decay rate of the latent reservoir, then it may be asked why complete suppression should be the goal of HAART. Despite our analysis suggesting that low-level viral replication does not affect the decay rate of the latent reservoir, ongoing viral replication in the setting of intense selective pressure by HAART may lead to selection of drug resistant viruses, which could potentially return the infection back or close to the original steady-state achieved by wild-type virus in the absence of drugs. Thus achieving suppression of replication, while not necessarily accelerating reservoir decay, is critical to preventing the evolution of drug resistance [59]. While a lack of clinically observed drug resistance in patients on HAART has been thought to reflect a lack of viral replication, a lack of viral evolution may also be observed if a viable drug-resistant mutant requires too many mutations to occur in a reasonable amount of time. The rare and unlucky HIV^{+} patient may archive an HIV mutant that is one mutation away from becoming completely drug resistant, which may occur from possible low-level viral replication. However, there is a highly significant correlation between patient non-adherence and development of multi-drug resistance, and patient non-adherence remains by far the most likely cause for the emergence of multi-drug resistance. Furthermore, because patient adherence is self-reported, non-adherence remains the most likely cause of multi-drug resistance even in reportedly adherent patients. Therefore, prevention of multi-drug resistance is not a validated rationale for intensified HAART [60]. However, future work will be helpful in understanding whether low-level viral replication in the setting of HAART may lead to multi-drug resistance.

## Conclusion

Our results have several implications. Standard antiretroviral medications have some toxicities and side effects, which may motivate patient non-adherence leading to the subsequent development of drug resistance. Intensification of HAART may also lead to intensification of drug toxicities and adverse effects. Therefore implementation of an intensified HAART regimen can only be justified if the benefit would outweigh the cost in patient morbidity. If there is essentially no risk of drug resistance in compliant patients on standard HAART regimens, then intensified HAART would have little to no benefit towards maintaining suppression of viremia. If the decay rate of the latent reservoir in patients on standard HAART regimens has reached a maximum, then an intensified HAART regimen would have no benefit toward eradication of the latent reservoir – suggesting that other approaches would be necessary for eradication of the reservoir. In contrast, HAART simplification (reduction in the number of drugs) has been recently discussed as an option for reducing the burden of treatment on patients who have already achieved sustained virologic suppression [61]. Based on our findings, we expect that a simplified HAART regimen would maintain the decay rate of the latent reservoir as long as long virologic rebound is not observed. Yet whether such a simplified HAART regimen would provide effective protection against resistance mutations is a different question, which is beyond the scope of our analysis. Our results suggest that the intrinsic dynamic properties of the reservoir are the primary factors that determine the decay rate of the latent reservoir. In particular, we find that the activation rate and intrinsic death rate of the latently infected cells primarily determine the decay rate of the latent reservoir and may serve as alternative targets to be exploited for accelerating the decay of the reservoir [62–64].

Eradication of the HIV infection is and will remain an extremely difficult undertaking. The complex dynamic relationship between viral replication, the latent reservoir, and the evolution of drug resistance, as well as management of these issues in the context of minimizing patient morbidity poses a daunting task for health care providers. Insight into these dynamic relationships will provide a greater understanding of the underlying principles guiding patient care and will hopefully lead to new approaches for achieving the ultimate goal of HIV eradication.

## Appendix

We assume that the number of latently infected cells does not change significantly over the time scale of oscillations that occur immediately after initiation of HAART. Therefore we set *dL*/*dt* = 0 and our model defined by equations 1–3 reduces to

We linearize equations A1 and A2 around the pre-HAART steady-state values $\overline{U}$ and $\overline{P}$ to obtain:

Based on the pre-HAART steady-state values $\overline{U}$ and $\overline{P}$ described in equations 4 and 5 and the approximation that *A*/*B* ≈ *δ*
_{
P
}(described in the main text), we can simplify equation A3 to

We determined the two eigenvalues of this system, which reflect the kinetic properties of the dynamics of *U* and *P*:

If *U* and *P* change without oscillations, then the eigenvalues of this system defined by equation A5 must have no imaginary parts. This restricts

and therefore

Based on numerical values of our model parameters, we determine that for approximately *δ*
_{
U
}≤ 0.1035, there will be no oscillations in the dynamics of *U* and *P*. The dominant factor that determines the lower bound on *δ*
_{
U
}in equation A7 is $2{\delta}_{P}\sqrt{\frac{\overline{L}}{\overline{U}}\frac{({\alpha}_{Q}+{\delta}_{L})}{{\alpha}_{R}}}$, which is several orders of magnitude larger than the other components of *δ*
_{
U
}. Examination of this expression reveals that the lower bound on *δ*
_{
U
}depends approximately linearly on *δ*
_{
P
}. The other major factor capable of producing large variation in the lower bound on *δ*
_{
U
}is the quantity $\frac{\overline{L}}{\overline{U}}$, which is variable from person to person and is also viral strain dependent. Over time it is also likely that the activation rate of latently infected cells decreases. Taking the limit as *α*
_{
Q
}goes to zero, we find that the lower bound on *δ*
_{
U
}is not significantly affected – decreasing by slightly more than a factor of two, from 0.1035 to 0.0437.

## References

- 1.
Chun TW, Finzi D, Margolick J, Chadwick K, Schwartz D, Siliciano RF: In vivo fate of HIV-1-infected T cells: quantitative analysis of the transition to stable latency. Nat Med. 1995, 1: 1284-1290. 10.1038/nm1295-1284.

- 2.
Chun TW, Stuyver L, Mizell SB, Ehler LA, Mican JA, Baseler M, Lloyd AL, Nowak MA, Fauci AS: Presence of an inducible HIV-1 latent reservoir during highly active antiretroviral therapy. Proc Natl Acad Sci USA. 1997, 94: 13193-13197. 10.1073/pnas.94.24.13193.

- 3.
Chun TW, Carruth L, Finzi D, Shen X, DiGiuseppe JA, Taylor H, Hermankova M, Chadwick K, Margolick J, Quinn TC, Kuo YH, Brookmeyer R, Zeiger MA, Barditch-Crovo P, Siliciano RF: Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature. 1997, 387: 183-188. 10.1038/387183a0.

- 4.
Finzi D, Hermankova M, Pierson T, Carruth LM, Buck C, Chaisson RE, Quinn TC, Chadwick K, Margolick J, Brookmeyer R, Gallant J, Markowitz M, Ho DD, Richman DD, Siliciano RF: Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. Science. 1997, 278: 1295-1300. 10.1126/science.278.5341.1295.

- 5.
Hermankova M, Siliciano JD, Zhou Y, Monie D, Chadwick K, Margolick JB, Quinn TC, Siliciano RF: Analysis of human immunodeficiency virus type 1 gene expression in latently infected resting CD4+ T lymphocytes in vivo. J Virol. 2003, 77: 7383-7392. 10.1128/JVI.77.13.7383-7392.2003.

- 6.
Lassen KG, Bailey JR, Siliciano RF: Analysis of human immunodeficiency virus type 1 transcriptional elongation in resting CD4+ T cells in vivo. J Virol. 2004, 78: 9105-9114. 10.1128/JVI.78.17.9105-9114.2004.

- 7.
Brooks DG, Hamer DH, Arlen PA, Gao L, Bristol G, Kitchen CM, Berger EA, Zack JA: Molecular characterization, reactivation, and depletion of latent HIV. Immunity. 2003, 19: 413-423. 10.1016/S1074-7613(03)00236-X.

- 8.
Finzi D, Blankson J, Siliciano JD, Margolick JB, Chadwick K, Pierson T, Smith K, Lisziewicz J, Lori F, Flexner C, Quinn TC, Chaisson RE, Rosenberg E, Walker B, Gange S, Gallant J, Siliciano RF: Latent infection of CD4+ T cells provides a mechanism for lifelong persistence of HIV-1, even in patients on effective combination therapy. Nat Med. 1999, 5: 512-517. 10.1038/8394.

- 9.
Siliciano JD, Kajdas J, Finzi D, Quinn TC, Chadwick K, Margolick JB, Kovacs C, Gange SJ, Siliciano RF: Long-term follow-up studies confirm the stability of the latent reservoir for HIV-1 in resting CD4+ T cells. Nat Med. 2003, 9: 727-728. 10.1038/nm880.

- 10.
Chun TW, Nickle DC, Justement JS, Large D, Semerjian A, Curlin ME, O'Shea MA, Hallahan CW, Daucher M, Ward DJ, Moir S, Mullins JI, Kovacs C, Fauci AS: HIV-infected individuals receiving effective antiviral therapy for extended periods of time continually replenish their viral reservoir. J Clin Invest. 2005, 115: 3250-3255. 10.1172/JCI26197.

- 11.
Dornadula G, Zhang H, VanUitert B, Stern J, Livornese L, Ingerman MJ, Witek J, Kedanis RJ, Natkin J, DeSimone J, Pomerantz RJ: Residual HIV-1 RNA in blood plasma of patients taking suppressive highly active antiretroviral therapy. JAMA. 1999, 282: 1627-1632. 10.1001/jama.282.17.1627.

- 12.
Frenkel LM, Wang Y, Learn GH, McKernan JL, Ellis GM, Mohan KM, Holte SE, De Vange SM, Pawluk DM, Melvin AJ, Lewis PF, Heath LM, Beck IA, Mahalanabis M, Naugler WE, Tobin NH, Mullins JI: Multiple viral genetic analyses detect low-level human immunodeficiency virus type 1 replication during effective highly active antiretroviral therapy. J Virol. 2003, 77: 5721-5730. 10.1128/JVI.77.10.5721-5730.2003.

- 13.
Furtado MR, Callaway DS, Phair JP, Kunstman KJ, Stanton JL, Macken CA, Perelson AS, Wolinsky SM: Persistence of HIV-1 transcription in peripheral-blood mononuclear cells in patients receiving potent antiretroviral therapy. N Engl J Med. 1999, 340: 1614-1622. 10.1056/NEJM199905273402102.

- 14.
Havlir DV, Strain MC, Clerici M, Ignacio C, Trabattoni D, Ferrante P, Wong JK: Productive infection maintains a dynamic steady state of residual viremia in human immunodeficiency virus type 1-infected persons treated with suppressive antiretroviral therapy for five years. J Virol. 2003, 77: 11212-11219. 10.1128/JVI.77.20.11212-11219.2003.

- 15.
Imamichi H, Crandall KA, Natarajan V, Jiang MK, Dewar RL, Berg S, Gaddam A, Bosche M, Metcalf JA, Davey RT, Lane HC: Human immunodeficiency virus type 1 quasi species that rebound after discontinuation of highly active antiretroviral therapy are similar to the viral quasi species present before initiation of therapy. J Infect Dis. 2001, 183: 36-50. 10.1086/317641.

- 16.
Natarajan V, Bosche M, Metcalf JA, Ward DJ, Lane HC, Kovacs JA: HIV-1 replication in patients with undetectable plasma virus receiving HAART. Highly active antiretroviral therapy. Lancet. 1999, 353: 119-120. 10.1016/S0140-6736(05)76156-0.

- 17.
Tobin NH, Learn GH, Holte SE, Wang Y, Melvin AJ, McKernan JL, Pawluk DM, Mohan KM, Lewis PF, Mullins JI, Frenkel LM: Evidence that low-level viremias during effective highly active antiretroviral therapy result from two processes: expression of archival virus and replication of virus. J Virol. 2005, 79: 9625-9634. 10.1128/JVI.79.15.9625-9634.2005.

- 18.
Zhang L, Ramratnam B, Tenner-Racz K, He Y, Vesanen M, Lewin S, Talal A, Racz P, Perelson AS, Korber BT, Markowitz M, Ho DD: Quantifying residual HIV-1 replication in patients receiving combination antiretroviral therapy. N Engl J Med. 1999, 340: 1605-1613. 10.1056/NEJM199905273402101.

- 19.
Ramratnam B, Mittler JE, Zhang L, Boden D, Hurley A, Fang F, Macken CA, Perelson AS, Markowitz M, Ho DD: The decay of the latent reservoir of replication-competent HIV-1 is inversely correlated with the extent of residual viral replication during prolonged anti-retroviral therapy. Nat Med. 2000, 6: 82-85. 10.1038/71577.

- 20.
Ramratnam B, Ribeiro R, He T, Chung C, Simon V, Vanderhoeven J, Hurley A, Zhang L, Perelson AS, Ho DD, Markowitz M: Intensification of antiretroviral therapy accelerates the decay of the HIV-1 latent reservoir and decreases, but does not eliminate, ongoing virus replication. J Acquir Immune Defic Syndr. 2004, 35: 33-37. 10.1097/00126334-200401010-00004.

- 21.
Hermankova M, Ray SC, Ruff C, Powell-Davis M, Ingersoll R, D'Aquila RT, Quinn TC, Siliciano JD, Siliciano RF, Persaud D: HIV-1 drug resistance profiles in children and adults with viral load of <50 copies/ml receiving combination therapy. JAMA. 2001, 286: 196-207. 10.1001/jama.286.2.196.

- 22.
Siliciano JD, Siliciano RF: A long-term latent reservoir for HIV-1: discovery and clinical implications. J Antimicrob Chemother. 2004, 54: 6-9. 10.1093/jac/dkh292.

- 23.
Strain MC, Gunthard HF, Havlir DV, Ignacio CC, Smith DM, Leigh-Brown AJ, Macaranas TR, Lam RY, Daly OA, Fischer M, Opravil M, Levine H, Bacheler L, Spina CA, Richman DD, Wong JK: Heterogeneous clearance rates of long-lived lymphocytes infected with HIV: intrinsic stability predicts lifelong persistence. Proc Natl Acad Sci USA. 2003, 100: 4819-4824. 10.1073/pnas.0736332100.

- 24.
Nettles RE, Kieffer TL, Kwon P, Monie D, Han Y, Parsons T, Cofrancesco J, Gallant JE, Quinn TC, Jackson B, Flexner C, Carson K, Ray S, Persaud D, Siliciano RF: Intermittent HIV-1 viremia (Blips) and drug resistance in patients receiving HAART. JAMA. 2005, 293: 817-829. 10.1001/jama.293.7.817.

- 25.
Persaud D, Siberry GK, Ahonkhai A, Kajdas J, Monie D, Hutton N, Watson DC, Quinn TC, Ray SC, Siliciano RF: Continued production of drug-sensitive human immunodeficiency virus type 1 in children on combination antiretroviral therapy who have undetectable viral loads. J Virol. 2004, 78: 968-979. 10.1128/JVI.78.2.968-979.2004.

- 26.
Kieffer TL, Finucane MM, Nettles RE, Quinn TC, Broman KW, Ray SC, Persaud D, Siliciano RF: Genotypic analysis of HIV-1 drug resistance at the limit of detection: virus production without evolution in treated adults with undetectable HIV loads. J Infect Dis. 2004, 189: 1452-1465. 10.1086/382488.

- 27.
Sedaghat AR, Siliciano JD, Brennan TP, Wilke CO, Siliciano RF: Limits on replenishment of the resting CD4+ T cell reservoir for HIV in patients on HAART. PLoS Pathog. 2007, 3: e122-10.1371/journal.ppat.0030122.

- 28.
Bonhoeffer S, May RM, Shaw GM, Nowak MA: Virus dynamics and drug therapy. Proc Natl Acad Sci USA. 1997, 94: 6971-6976. 10.1073/pnas.94.13.6971.

- 29.
Nowak MA, Bonhoeffer S, Shaw GM, May RM: Anti-viral drug treatment: dynamics of resistance in free virus and infected cell populations. J Theor Biol. 1997, 184: 203-217. 10.1006/jtbi.1996.0307.

- 30.
Perelson AS, Neumann AU, Markowitz M, Leonard JM, Ho DD: HIV-1 dynamics in vivo: virion clearance rate, infected cell life-span, and viral generation time. Science. 1996, 271: 1582-1586. 10.1126/science.271.5255.1582.

- 31.
Perelson AS, Essunger P, Cao Y, Vesanen M, Hurley A, Saksela K, Markowitz M, Ho DD: Decay characteristics of HIV-1-infected compartments during combination therapy. Nature. 1997, 387: 188-191. 10.1038/387188a0.

- 32.
De Boer RJ, Perelson AS: Target cell limited and immune control models of HIV infection: a comparison. J Theor Biol. 1998, 190: 201-214. 10.1006/jtbi.1997.0548.

- 33.
Rouzine IM, Sergeev RA, Glushtsov AI: Two types of cytotoxic lymphocyte regulation explain kinetics of immune response to human immunodeficiency virus. Proc Natl Acad Sci USA. 2006, 103: 666-671. 10.1073/pnas.0510016103.

- 34.
Wick D, Self SG: What's the matter with HIV-directed killer T cells?. J Theor Biol. 2002, 219: 19-31.

- 35.
Nowak MA, May RM, Anderson RM: The evolutionary dynamics of HIV-1 quasispecies and the development of immunodeficiency disease. AIDS. 1990, 4: 1095-1103. 10.1097/00002030-199011000-00007.

- 36.
Ribeiro RM, Mohri H, Ho DD, Perelson AS: In vivo dynamics of T cell activation, proliferation, and death in HIV-1 infection: why are CD4+ but not CD8+ T cells depleted?. Proc Natl Acad Sci USA. 2002, 99: 15572-15577. 10.1073/pnas.242358099.

- 37.
Grossman Z, Meier-Schellersheim M, Sousa AE, Victorino RM, Paul WE: CD4+ T-cell depletion in HIV infection: are we closer to understanding the cause?. Nat Med. 2002, 8: 319-323. 10.1038/nm0402-319.

- 38.
Krakauer DC, Nowak M: T-cell induced pathogenesis in HIV: bystander effects and latent infection. Proc Biol Sci. 1999, 266: 1069-1075. 10.1098/rspb.1999.0745.

- 39.
Kim H, Perelson AS: Viral and latent reservoir persistence in HIV-1-infected patients on therapy. PLoS Comput Biol. 2006, 2 (10): e135-10.1371/journal.pcbi.0020135.

- 40.
Bonhoeffer S, Coffin JM, Nowak MA: Human immunodeficiency virus drug therapy and virus load. J Virol. 1997, 71: 3275-3278.

- 41.
Pierson T, McArthur J, Siliciano RF: Reservoirs for HIV-1: mechanisms for viral persistence in the presence of antiviral immune responses and antiretroviral therapy. Annu Rev Immunol. 2000, 18: 665-708. 10.1146/annurev.immunol.18.1.665.

- 42.
Mohri H, Bonhoeffer S, Monard S, Perelson AS, Ho DD: Rapid turnover of T lymphocytes in SIV-infected rhesus macaques. Science. 1998, 279: 1223-1227. 10.1126/science.279.5354.1223.

- 43.
Kovacs JA, Lempicki RA, Sidorov IA, Adelsberger JW, Herpin B, Metcalf JA, Sereti I, Polis MA, Davey RT, Tavel J, Falloon J, Stevens R, Lambert L, Dewar R, Schwartzentruber DJ, Anver MR, Baseler MW, Masur H, Dimitrov DS, Lane HC: Identification of dynamically distinct subpopulations of T lymphocytes that are differentially affected by HIV. J Exp Med. 2001, 194: 1731-1741. 10.1084/jem.194.12.1731.

- 44.
Davey RT, Bhat N, Yoder C, Chun TW, Metcalf JA, Dewar R, Natarajan V, Lempicki RA, Adelsberger JW, Miller KD, Kovacs JA, Polis MA, Walker RE, Falloon J, Masur H, Gee D, Baseler M, Dimitrov DS, Fauci AS, Lane HC: HIV-1 and T cell dynamics after interruption of highly active antiretroviral therapy (HAART) in patients with a history of sustained viral suppression. Proc Natl Acad Sci USA. 1999, 96: 15109-15114. 10.1073/pnas.96.26.15109.

- 45.
Wei X, Ghosh SK, Taylor ME, Johnson VA, Emini EA, Deutsch P, Lifson JD, Bonhoeffer S, Nowak MA, Hahn BH: Viral dynamics in human immunodeficiency virus type 1 infection. Nature. 1995, 373: 117-122. 10.1038/373117a0.

- 46.
Muller V, Vigueras-Gomez JF, Bonhoeffer S: Decelerating decay of latently infected cells during prolonged therapy for human immunodeficiency virus type 1 infection. J Virol. 2002, 76: 8963-8965. 10.1128/JVI.76.17.8963-8965.2002.

- 47.
Maldarelli F, Palmer S, King MS, Wiegand A, Polis MA, Mican J, Kovacs JA, Davey RT, Rock-Kress D, Dewar R, Liu S, Metcalf JA, Rehm C, Brun SC, Hanna GJ, Kempf DJ, Coffin JM, Mellors JW: ART suppresses plasma HIV-1 RNA to a stable set point predicted by pretherapy viremia. PLoS Pathog. 2007, 3: e46-10.1371/journal.ppat.0030046.

- 48.
Sereti I, Herpin B, Metcalf JA, Stevens R, Baseler MW, Hallahan CW, Kovacs JA, Davey RT, Lane HC: CD4 T cell expansions are associated with increased apoptosis rates of T lymphocytes during IL-2 cycles in HIV infected patients. AIDS. 2001, 15: 1765-1775. 10.1097/00002030-200109280-00004.

- 49.
Kovacs JA, Vogel S, Metcalf JA, Baseler M, Stevens R, Adelsberger J, Lempicki R, Hengel RL, Sereti I, Lambert L, Dewar RL, Davey RT, Walker RE, Falloon J, Polis MA, Masur H, Lane HC: Interleukin-2 induced immune effects in human immunodeficiency virus-infected patients receiving intermittent interleukin-2 immunotherapy. Eur J Immunol. 2001, 31: 1351-1360. 10.1002/1521-4141(200105)31:5<1351::AID-IMMU1351>3.0.CO;2-9.

- 50.
Lehrman G, Hogue IB, Palmer S, Jennings C, Spina CA, Wiegand A, Landay AL, Coombs RW, Richman DD, Mellors JW, Coffin JM, Bosch RJ, Margolis DM: Depletion of latent HIV-1 infection in vivo: a proof-of-concept study. Lancet. 2005, 366: 549-555. 10.1016/S0140-6736(05)67098-5.

- 51.
van Praag RM, Prins JM, Roos MT, Schellekens PT, Ten B, Yong SL, Schuitemaker H, Eerenberg AJ, Jurriaans S, de WF, Fox CH, Goudsmit J, Miedema F, Lange JM: OKT3 and IL-2 treatment for purging of the latent HIV-1 reservoir in vivo results in selective long-lasting CD4+ T cell depletion. J Clin Immunol. 2001, 21: 218-226. 10.1023/A:1011091300321.

- 52.
Williams SA, Chen LF, Kwon H, Fenard D, Bisgrove D, Verdin E, Greene WC: Prostratin antagonizes HIV latency by activating NF-kappaB. J Biol Chem. 2004, 279: 42008-42017. 10.1074/jbc.M402124200.

- 53.
Persaud D, Pierson T, Ruff C, Finzi D, Chadwick KR, Margolick JB, Ruff A, Hutton N, Ray S, Siliciano RF: A stable latent reservoir for HIV-1 in resting CD4(+) T lymphocytes in infected children. J Clin Invest. 2000, 105: 995-1003.

- 54.
Ruff CT, Ray SC, Kwon P, Zinn R, Pendleton A, Hutton N, Ashworth R, Gange S, Quinn TC, Siliciano RF, Persaud D: Persistence of wild-type virus and lack of temporal structure in the latent reservoir for human immunodeficiency virus type 1 in pediatric patients with extensive antiretroviral exposure. J Virol. 2002, 76: 9481-9492. 10.1128/JVI.76.18.9481-9492.2002.

- 55.
Bailey JR, Sedaghat AR, Kieffer T, Brennan T, Lee PK, Wind-Rotolo M, Haggerty CM, Kamireddi AR, Liu Y, Lee J, Persaud D, Gallant JE, Cofrancesco J, Quinn TC, Wilke CO, Ray SC, Siliciano JD, Nettles RE, Siliciano RF: Residual Human Immunodeficiency Virus Type 1 Viremia in Some Patients on Antiretroviral Therapy Is Dominated by a Small Number of Invariant Clones Rarely Found in Circulating CD4+ T Cells. J Virol. 2006, 80: 6441-6457. 10.1128/JVI.00591-06.

- 56.
Bajaria SH, Webb G, Cloyd M, Kirschner D: Dynamics of naive and memory CD4+ T lymphocytes in HIV-1 disease progression. J Acquir Immune Defic Syndr. 2002, 30: 41-58. 10.1097/00126334-200205010-00006.

- 57.
Crowe S, Zhu T, Muller WA: The contribution of monocyte infection and trafficking to viral persistence, and maintenance of the viral reservoir in HIV infection. J Leukoc Biol. 2003, 74: 635-641. 10.1189/jlb.0503204.

- 58.
Schrager LK, D'Souza MP: Cellular and anatomical reservoirs of HIV-1 in patients receiving potent antiretroviral combination therapy. JAMA. 1998, 280: 67-71. 10.1001/jama.280.1.67.

- 59.
Coffin JM: HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science. 1995, 267: 483-489. 10.1126/science.7824947.

- 60.
Gulick RM, Ribaudo HJ, Shikuma CM, Lalama C, Schackman BR, Meyer WA, Acosta EP, Schouten J, Squires KE, Pilcher CD, Murphy RL, Koletar SL, Carlson M, Reichman RC, Bastow B, Klingman KL, Kuritzkes DR: Three- vs four-drug antiretroviral regimens for the initial treatment of HIV-1 infection: a randomized controlled trial. JAMA. 2006, 296: 769-781. 10.1001/jama.296.7.769.

- 61.
Swindells S, DiRienzo AG, Wilkin T, Fletcher CV, Margolis DM, Thal GD, Godfrey C, Bastow B, Ray MG, Wang H, Coombs RW, McKinnon J, Mellors JW: Regimen simplification to atazanavir-ritonavir alone as maintenance antiretroviral therapy after sustained virologic suppression. JAMA. 2006, 296: 806-814. 10.1001/jama.296.7.806.

- 62.
Nunnari G, Pomerantz RJ: IL-7 as a potential therapy for HIV-1-infected individuals. Expert Opin Biol Ther. 2005, 5: 1421-1426. 10.1517/14712598.5.11.1421.

- 63.
Wang FX, Xu Y, Sullivan J, Souder E, Argyris EG, Acheampong EA, Fisher J, Sierra M, Thomson MM, Najera R, Frank I, Kulkosky J, Pomerantz RJ, Nunnari G: IL-7 is a potent and proviral strain-specific inducer of latent HIV-1 cellular reservoirs of infected individuals on virally suppressive HAART. J Clin Invest. 2005, 115: 128-137. 10.1172/JCI200522574.

- 64.
Kulkosky J, Culnan DM, Roman J, Dornadula G, Schnell M, Boyd MR, Pomerantz RJ: Prostratin: activation of latent HIV-1 expression suggests a potential inductive adjuvant therapy for HAART. Blood. 2001, 98: 3006-3015. 10.1182/blood.V98.10.3006.

- 65.
Perelson AS: Modelling viral and immune system dynamics. Nat Rev Immunol. 2002, 2: 28-36. 10.1038/nri700.

- 66.
Nelson PW, Perelson AS: Mathematical analysis of delay differential equation models of HIV-1 infection. Math Biosci. 2002, 179: 73-94. 10.1016/S0025-5564(02)00099-8.

- 67.
Dixit NM, Perelson AS: HIV dynamics with multiple infections of target cells. Proc Natl Acad Sci USA. 2005, 102: 8198-8203. 10.1073/pnas.0407498102.

- 68.
Haase AT, Henry K, Zupancic M, Sedgewick G, Faust RA, Melroe H, Cavert W, Gebhard K, Staskus K, Zhang ZQ, Dailey PJ, Balfour HH, Erice A, Perelson AS: Quantitative image analysis of HIV-1 infection in lymphoid tissue. Science. 1996, 274: 985-989. 10.1126/science.274.5289.985.

- 69.
Cavert W, Notermans DW, Staskus K, Wietgrefe SW, Zupancic M, Gebhard K, Henry K, Zhang ZQ, Mills R, McDade H, Schuwirth CM, Goudsmit J, Danner SA, Haase AT: Kinetics of response in lymphoid tissues to antiretroviral therapy of HIV-1 infection. Science. 1997, 276: 960-964. 10.1126/science.276.5314.960.

- 70.
Ho DD, Neumann AU, Perelson AS, Chen W, Leonard JM, Markowitz M: Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature. 1995, 373: 123-126. 10.1038/373123a0.

### Pre-publication history

The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/8/2/prepub

## Acknowledgements

The authors would like to acknowledge Igor M. Rouzine for helpful and insightful discussion about the manuscript. C.O.W. was supported by NIH grant AI 065960. R.F.S. was supported by NIH grants AI43222 and AI51178, by a grant from the Doris Duke Charitable Foundation, and by the Howard Hughes Medical Institute.

## Author information

## Additional information

### Competing interests

The author(s) declare that they have no competing interests.

### Authors' contributions

ARS, RFS and COW conceived the model. ARS and COW designed and analyzed the model. ARS, RFS and COW interpreted the results of the model and wrote the manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

## About this article

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Viral Load
- Decay Rate
- Infected Cell
- Viral Replication
- Cell Dynamic