Previous Article | Next Article ![]()
Journal of Virology, September 2002, p. 8963-8965, Vol. 76, No. 17
0022-538X/02/$04.00+0 DOI: 10.1128/JVI.76.17.8963-8965.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Ecology and Evolution, ETH Zürich, 8092 Zürich, Switzerland,1 INRIA-Lorraine, 54600 Villers-les-Nancy, France2
Received 12 February 2002/ Accepted 20 May 2002
|
|
|---|
|
|
|---|
Latently infected CD4+ T lymphocytes have been shown to harbor replication-competent integrated provirus even after years of suppressive therapy (8, 12, 25). The pool of latently infected cells is established early during primary human immunodeficiency virus type 1 (HIV-1) infection (6), presumably by the reversion of productively infected memory CD4+ T cells into the resting state (20). In vitro these cells can be induced to produce infectious virus, and activation has been shown to occur also in vivo upon encounter with the appropriate recall antigens (14) or by the administration of activating lymphokines (7). In patients with apparently complete suppression of HIV-1 replication by HAART, the source of rebounding virus after the cessation of therapy was identified as the latently infected cell pool (26). However, the exact contribution of this pool to the residual virus production during fully suppressive HAART has not been clarified. Here we present a mathematical model on the dynamics of latently infected cells and their impact on residual virus replication. We show that a decelerating decay might be an inherent characteristic of this system.
Our goal was to model virus production by reactivated latently infected T cells after prolonged fully suppressive therapy when the initial major virus-producing population has already dwindled. This restricted goal allowed us to neglect much of the complexity of HIV-1 infection. We considered latently infected cells, L, which die at rate
LL and which are reactivated at rate
L, and cells actively producing virus, T*, which arise by the reactivation of latently infected cells and which die at rate
T*. The generation of latently or productively infected cells by new rounds of infection is assumed to be fully suppressed by the antiviral therapy. The reversion of cells actively producing virus to the resting state can also be neglected, as this cell pool plummets after the start of therapy. Virus particles, V, are produced at rate pT* and decay at rate cV. Thus we write
![]() | ((1a)) |
![]() | ((1b)) |
![]() | ((1c)) |
The clearance of latently infected cells reflects the sum of activation and death. Over time, one obtains a simple exponential decay as L(t) = L0
, where L0 is the initial size of the latently infected cell pool. The dynamics of latently infected cells are much slower than those of cells actively producing virus, which allowed us to assume a quasi-steady state for the latter. By setting equation 1b to zero we derived T*(t) = (
/
)L(t). Similarly, as virus particles turn over much faster than virus-producing cells (16, 18, 22), the virus level will follow the latter at a constant ratio: V(t) = (p/c)T*(t). In all, the virus level maintained by production from reactivated latently infected cells can be obtained as
![]() | ((2)) |
This form of the model predicts that virus production by this cell population declines parallel to the number of latently infected cells at a rate reflecting the sum of reactivation and death. However, the pool of latently infected cells is likely to be heterogeneous in the activation rate. Some cells are specific to common antigens and hence have a high probability to become activated, while other cells are specific for rare antigens and might persist for a long time without activation. This situation can be modeled by employing a continuous range of activation rates instead of a single fixed value. This transforms equation 2 into the form
![]() |
) defines the initial distribution of the latently infected cell pool with respect to the activation rate, which has a maximum value of
max. The total initial pool size of around 106 cells (4) can be obtained as L(0) =
L0(
)d
. It can be shown that for all biologically plausible distributions (L0[
] should be a smooth function of
and nonzero when
is close to zero; among others, uniform, normal, exponential, log normal, and gamma distributions fulfill these criteria) the behavior of the virus level is decelerating decay over time. The total number of latently infected cells declines according to the equation
![]() |
Figure 1a shows an example with a normal distribution of activation rates for the initial cell pool. The decline of the virus level maintained by reactivated cells and that of the latently infected cell pool decelerate over time. This is consistent with previous observations on the apparent biphasic decay of latently infected cells (2, 19). Figure 1a also indicates the slopes reported for the first few weeks (half-life around 6 to 8 days [2, 17]) and the following months (half-life around 6 months [23, 27]).
![]() View larger version (12K): [in a new window] |
FIG. 1. Total body load of virus and latently infected cells during HAART. (a) Number of latently infected cells (dashed line) and the level of virus maintained by these cells (heavy solid line). Light lines, decays with half-lives reported for the first few weeks of treatment (6 to 8 days) and for the following months (6 months). (b) Overall virus level (heavy solid line) including virus produced by the initial virus-producing population (dash-dotted line) and by long-lived persistently infected cells (dashed line). The first two phases of decay are dominated by these two cell populations, while the third phase is dominated by production from reactivated latently infected cells (light line). In both panels the initial population of latently infected cells is normally distributed, with parameters µ = 0.1d-1 and = 0.1d-1, and is scaled to yield a total of 106 cells. The distribution has been truncated at zero and renormalized to exclude negative activation rates. Further parameters and initial values were max = 1d-1, p/c = 100, L = 0.001d-1, pMM0 = 3 x 109d-1, M = 0.066d-1, T0 = 5 x 108, and = 0.7d-1.
|
L, which might be extremely slow (half-life of 43.9 months [11]). Note that our model predicts that the virus level decays somewhat faster than the latently infected cell pool (Fig. 1). The reason for this is illustrated by equation 2, which describes the simple case with a single fixed activation rate. The steady-state ratio of the virus level to the number of latently infected cells is linearly proportional to the activation rate. With a distributed activation rate, this implies that as the distribution shifts toward lower values, the ratio of virus level to the number of latently infected cells also declines over time, which results in a faster decay of the former.
To model the total virus level from the initiation of therapy, we needed to incorporate also the initial fast-producing cell populations into the model. In this we adopted the results of Perelson et al. (17) to obtain the total virus level as V(t) = VL(t) + (p/c)T0e-
t + (pM/c)M0
, where T0 is the initial population of cells actively producing virus, M0 is the initial population of long-lived persistently infected cells, and pM and
M are the virus production rate and death rate, respectively, of the latter cells. We set the initial values and the additional parameters to correspond approximately to the averages obtained by Perelson et al. (17). The graph of overall virus production was plotted in Fig. 1b. The result is an apparently triphasic decline in which the third phase reflects the long tail of the decelerating production by reactivated latently infected cells. Deceleration at this stage is already very slow, which results in an apparently linear third phase. However, the slope of third-phase virus decay always remains steeper than the innate death rate of latently infected cells due to the decreasing ratio of virus to cells.
A decelerating decay of virus production by latently infected cells is not the only possible explanation for the observed persistence and long-term dynamics of HIV-1 under HAART. It has been shown that low steady-state viral loads can be explained in terms of density-dependent infected-cell death rates or compartment models including a drug sanctuary or a privileged cell type where drugs have a reduced effect (3). The waning of the HIV-specific cytotoxic immune response may also result in a decelerating virus decay if the response has a strong effect on the life span of virus-producing cells (1). A predecessor to our model is the work by Ferguson et al. (10), who assumed a distributed activation state (in terms of division times and susceptibility to infection) of the primary virus-producing cell population. Finally, Grossman et al. (13) attributed the slower second phase to the decreasing frequency of local virus bursts and also showed that under imperfect therapy new infectious bursts initiated by reactivated latently infected cells may continuously reseed the reservoir and maintain a stable steady-state virus level.
To summarize, our model predicts that the source of residual virus production after long-term sustained suppressive treatment might indeed be the pool of latently infected cells. By its decelerating decay, virus production from this compartment eventually dominates production from any other compartment with a fixed turnover faster than the innate decay of latently infected cells. Importantly, however, the deceleration of virus decay can be explained by production from this single compartment without invoking multiple virus-producing cell types with graded turnover rates. We have also demonstrated that heterogeneity in the activation rate can explain the observed behavior of the latently infected cells themselves: fast initial decay, which eventually decelerates to the slow innate turnover of this compartment. Finally, testable predictions of our model are that (i) during treatment the distribution of the specificity of latently infected cells is gradually shifted toward rare antigens as the cells specific for common antigens are expected to be depleted faster and (ii) the ratio of virus level to latently infected cell level declines over time.
V.M. was partially supported by the Hungarian Scientific Research Fund (OTKA), and S.B. gratefully acknowledges financial support from the Swiss National Science Foundation.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»