Previous Article | Next Article ![]()
Journal of Virology, November 2005, p. 13579-13586, Vol. 79, No. 21
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.21.13579-13586.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Statistical Center for HIV and AIDS Research and Prevention, Fred Hutchinson Cancer Research Center, Seattle, Washington 98109,1 UCLA AIDS Institute, Department of Medicine, Department of Microbiology, Immunology and Molecular Genetics, David Greffen School of Medicine, UCLA, Los Angeles, California 90095,2 Department of Medicine, University of Washington, Seattle, Washington 98101,3 Program in Infectious Disease, Fred Hutchinson Cancer Research Center, Seattle, Washington 981094
Received 27 April 2005/ Accepted 26 July 2005
|
|
|---|
|
|
|---|
These data have prompted many vaccinologists to focus on cellular immunity. Yet, despite clear evidence that CTLs have a crucial antiviral function, their impact in vivo remains obscure. This lack of understanding presents an obstacle to defining the requirements for a successful CTL-based vaccine. Mathematical modeling of the interaction of CTLs and HIV-1 is one approach to elucidating these requirements.
For complicated biological processes, we often rely on mathematical models to explore mechanisms that are beyond direct experimental measurement. Most models require the representation of contributing processes by rate constants, to allow the evaluation of a mechanism of interest. Typically, some of these rate constants are estimated from the data used to build the model. But as model complexity increases, varying a large number of poorly defined parameters risks identifying the best-fitting wrong model. To avoid this pitfall, ideally rate constants should be derived from dedicated experiments, so that the "full model," when finally assembled, is not internally circular.
A key parameter in such a model is the efficiency of CTLs in eliminating infected target cells. A direct estimate is derived from observing the impact of patient-derived HIV-1-specific CTLs on replication in HIV-1-infected cell lines in vitro. From titrations of CTL density against viral growth, we can derive the killing rate by fitting a simplified model. To confirm the relevance of these in vitro measurements using cell lines, we estimated the same parameter by analyzing data from an in vivo "adoptive transfer" experiment conducted by Brodie et al. in 1999 (6, 7). Our results provide similar independent estimates of this key parameter for modeling the impact of CTLs in HIV-1 pathogenesis.
|
|
|---|
HIV-1-specific cytotoxic T-lymphocyte clones. HIV-1-specific CTL clones were obtained by limiting dilution cloning from peripheral blood mononuclear cells of infected individuals, characterized for specificity and HLA restriction, and maintained as previously described (30). Clone 68A62 recognized A*02-restricted reverse transcriptase (RT) epitope ILKEPVHGV (RT amino acids [aa] 476 to 484); 18030D23 recognized the A*02-restricted Gag epitope SLYNTVATL (Gag aa 77 to 85, p17); and clones LWC8 and 115M21 recognized the B*14-restricted Env epitope ERYLKDQQL (Env aa 584 to 592, gp41).
Inhibition assay. Assays for inhibition of HIV-1 replication by CTLs were performed as previously described (41). Briefly, target cells were infected at a multiplicity of 0.01 tissue culture infectious dose per cell and cultured at 5 x 105 target cells per well in 24-well plates with CTL clones at the indicated ratios. Supernatant p24 antigen concentrations were assayed by enzyme-linked immunosorbent assay at 2- to 4-day intervals (Dupont, Boston, MA).
Data from an in vivo study of CTL adoptive transfer. Published data from a 1999 study of infusions of HIV-1-specific CTLs were chosen for analysis because they provided a clear scenario where CTL levels changed with corresponding alterations in the concentrations of infected cells. Brodie et al. (6) derived Gag (p17 or p24)-specific CTL clones from three HIV-infected subjects on antiretroviral therapy with stable CD4 counts (224 to 261 per mm3). The CTLs were genetically modified to express the neomycin phosphotransferase gene (neo), and infusions (1 x 109 to 3 x 109 cells per m2 of body surface area) were administered a week apart. After the second dose, neo-modified cells constituted 2 to 3.5% of the patients' CD8 compartment. A few days after each infusion, the concentration of HIV RNA-positive cells decreased dramatically; but, in a few weeks, the neo-CTL disappeared and the infected-cell count rebounded. The data consist of baseline CD4 and the percentages of neo-marked CD8 cells and HIV-producing CD4 cells (productively infected targets [PITs]) at up to 14 time-points. Brodie et al. made several prebaseline measurements of PITs, which, except for patient 1, who seems to have had a spike at baseline, can be averaged to yield a steady-state infection rate. For patient 1, we used the value at day zero. Several CTL measurements were missing for patient 2, and we also omitted the final observation at day 72, which we expect is also an infection spike.
Mathematical models and fitting methods. For the in vivo analysis, we used a rate-equation (ODE) model with five compartments: the two observed, denoted C and P; uninfected, quiescent CD4+ T cells, denoted Q; uninfected, activated CD4+ T cells (targets) denoted T; and infected CD4+ T cells in the "eclipse" phase before transcribing HIV genes, denoted I. Rapid loss of CTLs after infusion allowed autonomous modeling of their effect, independent of the entire CTL activation circuit (involving antigen-presenting cells and CD4+ T-cell "help," for which few rate constants are available).
The model's kinetic equations express the rates at which cells enter or leave a compartment. Starting with the dynamics of the uninfected immune system: quiescent CD4+ T cells are created at rate
, die at rate
Q, and are activated at rate
. Activated cells are eliminated at rate
T, with a fraction,
, reverting to resting and the others undergoing programmed cell death (apoptosis). (We ignored the naïve versus memory distinction to restrict model complexity, but we might have included it as in reference 32.) Evidence from another experiment (see Discussion) indicated that the CTLs were unable to divide in vivo. Hence we assumed that the initial dynamics reflected cell dispersal and modeled it by a constant growth rate, from zero to the peak. For after-peak dynamics, we assumed a simple exponential decay with rate
C. On the infection side, target cells are infected at rate
P and progress to the productive phase at rate
. (We ignored the latently infected T-cell compartment, because the drop in viral load occurred too quickly to be explained by latent production.)
PITs are deleted by apoptosis or immune system killing, except by neo+ CTLs, at rate
P; killing by neo+ CTLs occurred at rate
C.
![]() | (1) |
![]() | (2) |
![]() | (3) |
![]() | (4) |
![]() | (5) |
t
1 and 7
t
8 and is zero otherwise; the function noinf(t) is zero on these intervals and otherwise 1. We included the unobserved compartments in order to exploit the baseline CD4 and viral load data. We omitted free virus, V, since viral dynamics is fast relative to cellular dynamics (free-virion lifetime in vivo may be less than a half-hour (11); effectively, V is proportional to P.
In this model, some rate constants are known approximately from other studies and are not likely altered by HIV infection or drug therapy, while others are either unknown or will likely be affected. For the former, we adopted standard values (see Table 1). For the latter, we cannot assume the thymopoesis rate, the activation rate, the infection coefficient, or the PIT lifetime is the same in Brodie et al.'s (6) infected, treated subjects as in healthy or untreated ones. The target of this investigation is the killing rate,
, while
C, the death rate of the neo-CTLs, can also clearly be estimated. Hence
,
,
, and
P, as well as the initial values of Q, T, and I, remain as unknown parameters. We resolved these issues as follows. We chose to estimate
P, the death rate (inverse lifetime) of PITs, and
, the thymopoesis rate, in order to compare with known estimates.
|
View this table: [in a new window] |
TABLE 1. Adopted standard values for the rate constants used in this study
|
and
, we derived by solving the steady-state equations of the model.
A steady state exists for equations 1 to 4, omitting C, with (observed) values of Pss and CD4 = Qss + Tss + Iss + Pss, provided 
A
P B > 0, where A = Pss and B =
Q (CD4 Pss) are fixed. Hence we adopted as estimable parameters
(
P,
,
C,
), yielding a rectangular search region. For each vector
, we derived the other parameters and the unobserved steady-state compartment values (used as initial conditions) from the equations
![]() | (6) |
![]() | (7) |
![]() | (8) |
![]() | (9) |
![]() | (10) |
![]() | (11) |
P,
,
C,
).
We assumed Gaussian measurement errors and minimized the sum of squares of differences between model predictions and data points (nonlinear regression). We fit the model by minimizing the sum of squares
![]() | (12) |
We used the fourth-order Runge-Kutta algorithm or an implicit solver to solve the system (25). For the minimization, we used a derivative-free, hill-climbing routine with random restarts (details available from the authors upon request).
To make confidence intervals, we computed the Fisher information matrix F:
![]() | (13) |
/
p X(t); they were obtained by solving the system
![]() | (14) |
/
p X(0). (Recall that the initial conditions depended on a parameter to be estimated.) Vec stands for the right-hand sides of equations 1 to 5, and J (for Jacobian) is the matrix obtained by differentiating Vec with respect to the coordinates. We proceeded by standard likelihood theory to invert F and use the diagonal entries to construct the Wald-type 95% confidence intervals for
P,
,
C, and 90% confidence intervals for
(which is a linear combination of
P and
). Finally, we assessed convexity of the likelihood by starting the function minimizer at random points within the search region, to ascertain if the same minimum was always reached (it was on 10 restarts); also, we checked that the information matrix evaluated at the minimum was positive definite.
In order to discuss the steady state with a given level of neo-CTL killing, we solved equations 1 to 4 with C fixed to obtain Pnew:
![]() | (15) |
=
Q +
and
=
P +
C. (Steady states in the model are neutrally stable, an artifact of omitting free virus; nevertheless they attract a perturbed state.)
For the in vitro analysis, we reduced the basic model to three compartments: uninfected target cells, infected cells in eclipse phase, and PITs. That is, we omitted Q and equations 1 and 5, assumed C = rT(0) is constant, set
T = 0, and added a term,
T, in equation 2 representing the growth rate of target cells. In equation 4, we replaced
C with
C/(1 +
C); parameter
allows for saturation at a high effector/target ratio. Substituting rT(0) for C, we performed a nonlinear regression of observed log inhibition (LI): LI = log[V(0; r)/V(8; r)], where V(t; r) is virus on day t measured in any units, e.g., nanograms/ml of p24 antigen, versus log[P(0; r)/P(8; r)], where P(t; r) is the predicted density of PITs. We estimated
and
while keeping other rate constants fixed.
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. Dose-response and model fit to the in vitro data for two CTL clones. Numbers on the x axis are 100 times the CTL/target ratio in the well, reflecting a multiplicity of infection of 0.01; numbers on the y axis are the log10 inhibitions of p24 antigen, relative to controls, after 8 days.
|
As simple regression models are inappropriate here, we employed for the statistical estimation a three-compartment ODE model of uninfected target cells, infected cells in the eclipse period, and PITs. In this analysis, sensitivity to unknown parameters mattered more than statistical error. (Indeed, estimating both killing and saturation coefficients, the "best" curves nearly passed through the observations [Fig. 1].)
The point estimates were not sensitive to varying the basic reproductive number (R0) but were to varying PIT lifetime (parameter
P) and eclipse period (parameter
); neither was known accurately in these cell lines. So we performed a sensitivity analysis, varying PIT lifetime from 2 to 3 days and the eclipse period from 1 to 3 days. We exhibit point estimates and sensitivity ranges in Table 2; Fig. 1 shows several model fits.
|
View this table: [in a new window] |
TABLE 2. Point estimates and sensitivity ranges
|
Analysis of adoptive transfer experiments demonstrates similar killing rates in vivo.
While analysis of the antiviral activity of CTLs against HIV-1 in vitro affords a controlled setting for defining the killing rate, this model is a highly contrived experimental system. We therefore obtained independent estimates of this parameter from experiments reported by Brodie et al. (6). The estimated parameters and confidence intervals are reported in Table 3; we display one model fit in Fig. 2. (The other two were similar.) The technical definition of the killing parameter in vivo is the following. Let CTLs specific for an HIV antigen, with density denoted "C," kill PITs, with density denoted "P," at a rate proportional to the product P · C. The sought after coefficient
is simply the constant of proportionality.
|
View this table: [in a new window] |
TABLE 3. In vivo estimated parameters
|
![]() View larger version (11K): [in a new window] |
FIG. 2. PIT and neo-marked CTL densities and model fit for subject 1, from the study by Brodie et al (6).
|
0.129 µl cell1 day1 (mean of the point estimates in Table 3; the "meta-analytic" 95% confidence interval of the mean was 0.19, 0.44).If total body cell counts are used instead of densities,
2 x 1010 day1.
For a more concrete representation, consider the killing rate per CTL per microliter. The subjects of this study had about 250 CD4+ T lymphocytes per µl, up to 2% productively infected, providing up to 5 PITs per µl in peripheral blood. From the approximate equation
P/C =
P
t, it is apparent that each CTL was able to kill an average of about 0.65 PIT per day.
|
|
|---|
![]() View larger version (7K): [in a new window] |
FIG. 3. Theoretical prediction from the model for the impact on PITs in subject 1, if the infused CTLs had persisted at peak level (i.e., setting C = 0).
|
Despite the observed effects on the lifetime of PITs, the transferred CTLs in the study by Brodie et al. (6) did not perturb the established steady state sufficiently to abolish viremia. According to our model, even if the infused CTLs persisted at their peak level, the impact on PITs would have been about a 10-fold-lowered density (Fig. 3). This is considerably less than the effect required to drive HIV-1 from steady state to extinction (disregarding the latent reservoir). That goal would require about a 30% frequency of these CTLs in the CD8+ T-lymphocyte compartment. The observation that viremia was not appreciably reduced in this study is consistent with the modest reduction of PITs, only to 0.5% of the CD4+ compartment, still in the range observed in chronic infection. In the context of the above hypotheses regarding the threshold of CTLs needed to prevent acute infection, this likely reflects a greater difficulty in destabilizing an already existing steady state.
Many of these conclusions about steady states and eradication follow simply from contemplating the reproductive number (number of daughter PITs produced by one mother PIT in its lifetime) in the simplest infection model (i.e., ignoring eclipse phase):
![]() | (16) |
P is replaced by
P +
C, which pushes R below 1. The PIT density will subsequently fall. However, that does not imply that the possibility of a steady state is abolished; for as infection lessens T will increase. In fact, in the infected steady state T was small, about 0.13 cell per µl; i.e., most activated CD4+ T cells were infected. A new steady state may form at a lower PIT density and a higher target cell density. That is in fact what the model predicts; solving the steady-state equations with C fixed at the peak level gave, e.g., for subject 1, a new solution with P about 0.3% of CD4.
For natural infection in untreated individuals,
and T will be larger and
P smaller than in the patients of Brodie et al. The reason that
., the initial PIT death rate in natural infection, will be smaller is that endogenous CTLs will not yet be activated or expanded. All three changes increase R, yielding the basic reproductive number, denoted R0; as we mentioned, R0
3 is reasonable. If we imagine having stimulated Gag-specific CTLs in uninfected individuals to the peak neo-CTL level, with the same killing potential, then R0 will be multiplied by the factor
![]() | (17) |
Caveats to our analyses include issues of statistical power, data selection, and model choice. Concerning statistical issues, we note that, in the in vivo study, the thymopoiesis rate,
, was particularly hard to estimate; the stability of the point estimates is surprising given the wide confidence intervals. Observed T-cell replacement rates range from, depending on health and age of the subject, 0.5 to 5 cells per µl per day, (18, 22, 23), so the estimates are high but reasonable. (Some theories of HIV pathogenesis even support an increased T-cell production [10].) Activation rates in the steady state fit to the baseline data (not included in Table 2) came out to be about 3% of CD4, somewhat elevated over the normal 1% (26), as is usually observed in HIV-infected patients (4, 24). In the in vitro study, the killing rate estimates, which reflect the slopes near zero (e.g., in Fig. 1), were essentially based on 3 data points per clone. Concerning the variation of estimated killing rates, a crude ranking, based on maximum height of the inhibition curves, would give 18030D23 > LWC8 > 68A62 > 115M21, but the ranking based on estimated
is 18030D23 > 115M21 > LWC8 > 68A62. The reason that 115M21 rose in this rankings is the single observation at an effector/target ratio of 1.56 (dilution factor, 64:1).
Two other objections to our mathematical techniques are that we assumed only passive measurement errors and that we should have described infection by a stochastic, branching-type model. To the former, we note that, for fitting nonlinear, higher-dimensional, active-noise models to data, adequate methods are lacking at this time (but we have proposed one [38]). To the latter, we remark that, with billions of infected targets and HIV-recognizing CTLs in the study by Brodie et al. and hundreds of thousands in the in vitro experiments, laws of mass action should be adequate. In contrast, for primary HIV, where both infection and response begin at low frequencies, a stochastic process is more appropriate (e.g., see reference 33).
Concerning data selection, we can compare our approach with prior work. The investigators in reference 19 depleted CD8 cells from macaques with monoclonal antibodies; when infected with an SIV-HIV chimera, the monkeys suffered elevated viral loads relative to controls. In reference 27, the investigators similarly perturbed the course of events in primary infection with SIVmac, but also depleted CD8 cells from chronically infected animals. Virus replication peaked again in these animals, but decreased when virus-specific CTLs reappeared. Unfortunately, neither work presented sufficient longitudinal data in the form we required for fitting models; counts of PITs and virus-specific CTLs were the crucial elements for our in vivo analysis. The authors of reference 12 reported a similar experiment and also addressed the theoretical mechanisms in terms of models. However, they concluded that the rapid rise in viremia could not be explained in their model by decreased killing, but rather by increased production. Another possible confounding factor in the CD8-depletion experiments is that the administration of the monoclonal antibody or the destruction of CD8 cells might have produced a general inflammationresulting in increased CD4 activation and viral production.
In reference 8, the authors studied CTL impact on SHIV-89.6P infection in vaccinated and unvaccinated macaques and concluded that the slopes after peak viremia were statistically indistinguishable between the two groups. Even granting that the vaccine acted by amplifying SHIV-specific CTLs, which moreover were potent killers, their method differs sufficiently from ours (exponential growth or decay equations fitted separately to prior- or post-peak viral load and CTL data, as opposed to nonlinear models fitted jointly to perturbed levels of both populations from an existing steady state), that we cannot assess the reason for the disparity.
Concerning model choice, a key assumption made in order to analyze the data of Brodie et al. was that the CTLs did not divide in vivo. This assumption is supported by an experiment in the humanized mouse system (20). CTLs were labeled with the dye CFSE, which binds to structural proteins and hence is evenly divided between daughter cells at division (17). The investigators observed that transferred CTLs did not divide in vivo. Moreover, the CTLs were rapidly lost in the HIV-infected host, probably due to apoptosis. This experiment should be repeated in the SIVmac system. (In reference 35, it is argued that performing the experiment in monkeys would also help resolve the origin of a putative CTL "memory defect.")
An additional concern in the adoptive transfer experiments is that the engineered CTLs might have had impaired function, either due to the manipulations ex vivo or lack of CD4+ T-cell helpso we have probably underestimated the natural killing efficiency of CTLs in vivo. Indeed, CTLs specific for HIV and active in chronic infection are probably not representative of CTLs in primary infection (16) or CTLs active against other viruses, due to a memory or killing defect (modeled in reference 35).
In summary, taken together with the observation that CTLs appear simultaneously with the initial drop in viremia (15), our results provide strong positive evidence that control of viremia in HIV infection is a consequence of CTL activity. (In reference 34, the authors provided a complementary negative argument that control is not due to target cell limitation.) Further studies are needed to elucidate the contributions of CTLs with various specificities and to explore whether enhancing the most potent killers is the key to a successful T-cell vaccine.
We thank S. Brodie and P. Greenberg for valuable discussions about the adoptive transfer experiments and B. Walker and S. Kalams for providing the CTL clones for the in vitro assays.
|
|
|---|
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»