Abstract
It was recently shown that the expansion of CD4+ T cells during a primary immune reaction to a peptide from cytochrome c decreases ∼0.5 log for every log increase in the number of cognate precursor cells, and that this remains valid over more than four orders of magnitude (Quiel et al. 2011. Proc. Natl. Acad. Sci. USA. 108: 3312–3317). This observed “power law” was explained by a mechanism where nondividing mature T cells inhibit the proliferation of less-differentiated cells of the same specificity. In this article, we interpret the same data by a mechanism where CD4+ T cells acquire cognate peptide-MHC (pMHC) complexes from the surface of APCs, thereby increasing the loss rate of pMHC. We show that a mathematical model implementing this “T cell grazing” mechanism, and having a T cell proliferation rate that is determined by the concentration of pMHC, explains the data equally well. As a consequence, the data no longer unequivocally support the previous explanation, and the increased loss of pMHC complexes on APCs at high T cell densities is an equally valid interpretation of this striking data.
Introduction
T cell responses to Ags are characterized by an expansion phase due to the Ag-specific activation and clonal expansion of cognate T cells, and a contraction phase during which cell densities approach a stable level that is called the memory phase (1). Several studies have demonstrated that the level of expansion during the expansion phase depends on the number of precursor cells recruited into the immune response (2–8). The most quantitative of these studies was recently performed by Quiel et al. (7), who varied the initial precursor density by >4 orders of magnitude, that is, from 1 to >30,000 cells/mouse, and measured the fold expansion at day 7 of the response (i.e., somewhere in the contraction phase). Their most striking result was that in a log–log representation, the fold expansion declines linearly with the precursor density over this wide range (see Fig. 1A); that is, the fold expansion drops ∼0.5 log for every log increase in the number of precursors (7). Because physiologically measured precursor frequencies cover a similar range, the tight relationship between the fold expansion and the precursor density was taken to reflect a fundamental property of Ag-mediated immune responses (7). Phenomenologically, the population density after 1 wk was very well described by cells (see Fig. 1B), where N(0) is number of naive precursors in the mouse (8), reflecting the same 0.5 log decrease in the fold expansion per log increase in the precursor density. Time courses of responses starting with 300 or 30,000 cells revealed that the peak of the response occurred earlier when the initial precursor density was higher (Fig. 1C).
The experimental system involved the adoptive transfer of CD4+ 5C.C7 TCR transgenic T cells that are specific for a peptide from cytochrome c peptide presented on I-Ek. The magnitude of the response was measured by an extremely sensitive PCR analysis specifically amplifying the CDR3 region of the Vβ-chain of the 5C.C7 TCR, allowing the detection of the progeny of a single cell in a mouse (7). Importantly, the effect shown in Fig. 1 was not due to a reduced recruitment of precursor cells at precursor high densities, because the difference in clonal expansion came about at a late stage in the response. Measuring the fold expansion at day 2, that is, before the release of proliferated cells from the lymphoid tissue, during responses starting with 300 or 30,000 precursors, it was found that both populations expanded ∼10-fold (7). Recording cell division by sacrificing mice at days 3.5, 4.5, 5.5, 6.5, and 7.5, after pulsing them with BrdU for their last 6 h, it was found that ∼60% of the cells were BrdU+ at days 3.5 and 4.5 when the response started with 300 cells, and that these percentages decline to 40 and 23%, when the response started with 30,000 cells (see Fig. 2D) (7). Thus the fraction of dividing cells is smaller and declines earlier when responses start with more precursor cells, and such a difference in the expansion rates can account for the observed effect (7).
The effect did not depend on regulatory T cells, because: 1) there were <1% Foxp3+ among the 5C.C7 cells, and this percentage did not depend on the precursor density; 2) the percentage Foxp3+ cells among the endogenous CD4+ T cells was normal (i.e., ∼10%) and did not depend on the precursor density; and 3) similar fold expansions were observed in mice receiving OT-II cells from wild type donors or from Scurfy donor mice that have no functional regulatory T cells (7). Several cytokines, like IFN-γ, IL-2, IL-7, and IL-15, and various inhibitory molecules, like Fas, TNF-α receptor, and CTLA-4, had hardly any effect on the differences in the fold expansion. Importantly, the effect was Ag specific because adding large numbers of T cells of another specificity hardly affected the fold expansion (7). Increasing the Ag dose, or increasing the density of CD11+ dendritic cells, increased the fold expansion, but for all precursor densities in a similar manner, implying that the observed reduction of the fold expansion cannot be caused by limiting Ag presentation at high precursor numbers only (7).
The data were interpreted by developing a mathematical model that allowed cognate T cells, after their priming by Ag, to differentiate from slowly dividing cells, to rapidly proliferating cells, to nondividing cells, to nondivided mature cells (8). The mature cells in this developmental cascade downregulate the differentiation of slowly dividing cells into rapidly proliferating cells. Intuitively, one can understand that this allows for an initial expansion that remains fairly independent of the initial precursor frequency, because the regulatory mature cells only appear later. When they appear, they readily account for the observed earlier and increased downregulation of the expansion at large population sizes (Figs. 1C, 2D). Combining the data and the modeling, it was suggested that the data are best explained by downregulatory feedback from differentiated (memory) cells on the clonal expansion of less-differentiated cells (7, 8).
The mathematical model defined a cascade from proliferating to mature nondividing cells, and consisted of two populations of dividing cells and two populations of mature cells, that is, of four differential equations (8). Because of the inclusion of several feedback loops, and proliferation functions allowing for competitive effects of cell density on the division rate, the model contained 15 parameters that were all estimated by fitting the model to the data. The model described the data well, and it was found that the fit remained of similar quality when one of the feedback loops was eliminated from the model. This simplified model had 10 parameters that were also estimated by fitting. The mere fact that the models describe the data well demonstrates that the interpretation in terms of downregulatory effects from differentiated cells on less-differentiated dividing cells is a possible explanation. Because both models are complicated, and have a large number of parameters, the excellent agreement between model and data fails to prove, however, that this is the only interpretation of the data. The authors indeed acknowledge that they cannot exclude other explanations like competition between cognate cells for presented peptides (7).
One alternative explanation for the observed reduction in the fold expansion is an increased loss of peptide-MHC (pMHC) complexes on APCs at high densities of precursor cells. During cognate interactions with APCs, CD4+ T cells tend to acquire a variety of cell-surface molecules, including the specific pMHC complexes binding the TCRs in the immunological synapse (9–14). This leads to a form of Ag-specific competition between T cells binding the same pMHC on the same APCs (10, 15). This effect would be perfectly consistent with the observations of Quiel et al. (7) because T cells of another specificity hardly affected the fold expansion. Their observations that increasing the Ag concentration, or the density of CD11+ dendritic cells, increased the fold expansion at all precursor densities in a similar manner (7) indeed suggest that pMHCs on APCs become a limiting resource at all precursor densities tested. We therefore investigate whether these striking data are also consistent with competition for presented peptides in a system where the pMHC are being grazed by other cells of the same specificity. Our approach is to develop a caricature model that only implements the essentials of clonal expansion and contraction, and the loss of pMHC by cognate T cells binding APCs. Because this simple model with five or seven parameters describes the data equally well, we conclude that the data equally support this alternative explanation.
Materials and Methods
To investigate whether the data presented in Fig. 1 can also be explained in terms of a model where T cells divide in response to pMHC, and where the availability of pMHC declines as a function of the number of T cells that form conjugates with the APCs, we developed a simple model describing the number of 5C.C7 CD4+ T cells in a mouse, N, and the concentration of pMHC, C, on an average APC. The concentration of pMHC is scaled by its maximal value, that is, upon immunization with cytochrome c, we set C(0) = 1, and let it decline afterward. pMHC decline by normal turnover and are consumed by the 5C.C7 CD4+ T cells. Because APCs can only bind a limited number of T cells at one time, the rate at which T cells consume pMHC should saturate at large T cell densities, and we model this by a Hill function with saturation constant h and shape parameter n. We model proliferation by assuming that individual T cells divide at a rate proportional to their cognate pMHC concentration.
Data were presented by Quiel et al. (7). (A) The fold expansion at day 7, N(7)/N(0), as a function of the number of precursors, N(0). Dashed line is the regression line 2,512 × N(0)−0.54 fitted by Bocharov et al. (8). (B) The population size at day 7 as a function of the number of precursors, N(0). Dashed line is the equivalent regression line 2,512 × N(0)0.46 fitted by Bocharov et al. (8). Data were taken from Supplemental Table VI in Bocharov et al. (8) and by digitizing the data from Fig. 1d in Quiel et al. (7). (C) The time-course data of responses starting with N(0) = 300 (bullets) or N(0) = 3,0000 (triangles) 5C.C7 T cells, and the best fit to these data by the simplified model of Bocharov et al. (8) (solid and dashed lines). Data were taken from Supplemental Table III in Bocharov et al. (8).
The pMHC concentration depends on the de novo generation of peptides from the depot of cytochrome c that was injected at day zero, and which we assume declines exponentially at rate kc. For the average concentration of pMHC on an APC, we write(1)where P is the rate at which pMHC are produced from cytochrome c, kd is the dissociation rate of peptides from MHC molecules, and kg is the maximum rate at which T cells graze pMHC. The dissociation rate, kd, typically ranges between several hours and >10 h (16–19). Because we have no information on the concentration of cytochrome c over time, or on the rate at which pMHC are produced from the full protein, we scale the pMHC concentration by its maximum concentration P/kd, so that C is at most 1. In terms of this scaled concentration, we obtain in the absence of T cells
which for C(0) = 1 obeys the solution
(2)which defines a biphasic decline from the initial value C(0) = 1 to zero. The addition of T cell grazing can only increase the rate at which C(t) declines, and we indeed obtain that 0 < C(t) ≤ 1. This has the immediate advantage that we can define the maximum proliferation rate of the 5C.C7 CD4+ T cells as a proportionality constant, p, reflecting the initial division rate of a 5C.C7 CD4+ T cell (in days). The full model can be written as the following two differential equations:
(3)where N is the total number of cognate CD4+ T cells in a mouse, C is the scaled concentration of pMHC on an APC, p and d are proliferation and death rates, respectively, and kc, kd, and kg are loss rates of pMHC (all rates are per day; see Table I). We initialize the model with varying numbers of precursors cells by setting N(0) to various values between 1 and
cells/mouse. A standard dose of cytochrome c corresponds to an initial pMHC concentration of C(0) = 1, and an increased dose of Ag is modeled by setting C(0) = 1.1.
In the Results, we will show that Eq. 3 describes the data well even if we set the parameter ; that is, when the loss of pMHC concentration from its initial value C(0) = 1 is completely because of the grazing by the T cells. This means that at high T cell densities (i.e., whenever N>>h), we obtain that dC/dt ≃ − kgC. This yields an exponential decay of pMHC,
. Because the peak of the T cell response occurs when pC = d, we can solve for the time of the peak from
[when C(0) = 1]. Thus, large T cell populations are predicted to start with an expansion phase and peak at day t = ln[p/d]/kg. The number of T cells at the peak, Nmax, can then be obtained by substituting the time of the peak into the full solution of dN/dt, that is:
(4)The model was simulated using the variable time step Runge–Kutta integrator provided by Press et al. (20), and was fitted to the data using the Levenberg–Marquardt algorithm (21). Confidence intervals were obtained by bootstrapping the residuals 50 times (22).
Results
The model, Eq. 3, was fitted to the data from the two published time courses, starting with either 300 or 30,000 T cells (Fig. 2A), and to the two published data sets on the measured fold expansion at day 7 for various different precursor densities (Fig. 2B). The new model describes the data well. One salient feature of the data were that for both initial conditions, the population size expanded 10-fold during the first 2 days of the response (7). Comparing the corresponding initial expansions in Fig. 2A reveals that the new model captures this very well. Fitting the data, it turned out that most of the net loss of pMHC was estimated to be due to the uptake by T cells, and that the quality of the fit was hardly reduced when the normal production and turnover of pMHC was ignored by setting kd = 0 (see Eq. 3). Setting kd implies that the kc parameter becomes irrelevant and that we can fit the data with just five parameters (see Table I). Conversely, the model fails to describe the data if we do not allow for consumption of pMHC by T cells, by setting kg = 0. Even the five remaining parameters have large 95% confidence limits because a high maximum rate, p, of T cell proliferation can be compensated by a high saturation constant, h, and/or a high death rate, d (see Table I). The best estimate of the proliferation rate of p = 2.7 per day is a reasonable value, however, as similar values are obtained when CD4+ T cells proliferate vigorously during an acute lymphocytic choriomeningitis virus infection (23). The fact that Eq. 3 with only five parameters describes the data well (see Fig. 2A, 2B) demonstrates that the loss of pMHC by the uptake of cognate T cells is a sufficient explanation for the striking observations made by Quiel et al. (7).
Comparing the data with the best fit of the new model. (A) The time course of the total number of T cells starting with N(0) = 300 (bullets and solid line) or N(0) = 30,000 (triangles and dashed line) CD4+ 5C.C7 T cells. Because the time courses run fairly parallel until day 3, this is in good agreement with the data showing a similar fold expansion at day 3 (7). (B) The fold expansion, N(7)/N(0), as a function of the number of precursors, N(0). The heavy line is the prediction of the new model, and the dashed line is the regression line 2,512 × T(0)−0.54 fitted by Bocharov et al. (8). (C) The pMHC concentration, C, as a measure of the fraction cells entering division. (D) The measured fraction of BrdU+ T cells after a short pulse of 6 h, which should also be proportional to the fraction of cells entering division. The BrdU data were taken from Supplemental Table IV in Bocharov et al. (8). (E) The effect of 10% increase in the initial Ag concentration. (F) The maximum clone size during the response predicted by the new model (solid line). The dashed line is the regression line 3,981 × N(0)−0.48 fitted to the prediction of the Bocharov et al. (8) model. Parameters were estimated by fitting Eq. 3 to the data in (A) and (B). Their values and confidence intervals are given in Table I. Becaue setting kd = kc = 0 did not affect the quality of the fit, we have fitted the data with five free parameters.
To test whether our alternative explanation is also consistent with the other data in the original publication (7), we compare the activation rate of the T cells in our model with the BrdU data taken at various time points in the mice (Fig. 2D). Because the fraction of cells entering division is proportional to the concentration of pMHC, we depict in Fig. 2C the predicted time course of the pMHC decay after starting with 300 or 30,000 precursor cells. Although the actual values of the pMHC concentration in the model (Fig. 2C) cannot be compared directly with the fraction of cells becoming labeled during a 6-h pulse with BrdU (Fig. 2D), the model and the data are in general agreement as the predicted division rate, pC, is somewhat lower and declines earlier when the immune response starts with 30,000 cells, as is seen in the data (Fig. 2D).
The model also correctly predicts that the fold expansion increases similarly for all precursor densities if the Ag concentration is increased. Increasing the initial pMHC concentration by 10%, we find the same increase in the fold expansion for a wide range of precursor densities (Fig. 2E), such that the dependence of the fold expansion on the initial density is not affected. Like Bocharov et al. (8), we can use the model to predict the fold expansion at the peak of the response (see Fig. 2F), rather than that at day 7 (which is somewhere in the contraction phase; see Fig. 2A). However, because there is very little data on the maximum size of the response, these two predictions cannot be compared with data. An interesting difference between the two models is that at very large precursor densities, the regression model of Bocharov et al. (8) predicts that the maximum clone size will be smaller than the number of precursor cells, whereas Eq. 3 predicts a 13-fold expansion at the peak of the response when the precursor density is very large (Fig. 2F). The fold expansion for large precursor densities can be solved analytically (see Eq. 4), and for the estimated parameters, the peak corresponds to a 13-fold expansion occurring around day 4. Close inspection of the fit of the regression line of Bocharov et al. (8) to the predicted data (their Fig. 2D) indeed reveals a deviation at large precursor numbers, and solving their simplified mechanistic model numerically confirmed that their mechanism also predicts that there should always be an initial expansion phase.
Discussion
A simple model where the loss of pMHC complexes increases with the density of cognate CD4+ T cells is capable of explaining the observed 0.5 log reduction in the fold expansion of T cells with every log increase in the precursor density over more than four orders of magnitude (7, 8). Because T cells are known to specifically pick up cognate pMHC from APC (9, 11–14), we provide a simple and realistic interpretation of the data. Hence, our main result is that consumption of pMHC by cognate T cells is a natural alternative interpretation of the data.
The saturation function that describes the net loss rate of pMHC as a function of the number of grazing T cells estimated by the fitting procedure is depicted in Fig. 3. Unfortunately, this function remains phenomenological, and the exponent of the Hill function, n = 0.31, has no mechanistic underpinning. Because we have been unable to estimate the production rate of pMHC and the normal dissociation rate of peptides from the data, the quantitative interpretation of the rate at which T cells graze pMHC also remains uncertain. One could make the model more realistic by including an initial time lag between the onset of CD4+ T cell proliferation, and/or explicitly account for the uptake of cytochrome c by Ag processing cells, its fragmentation into short peptides, the transportation of peptide-loaded MHC molecules to the cell surface, and the full activation of the precursor cells. Finally, it is likely that these CD4+ T cells complete several divisions after their initial activation, and only need to see pMHC again later in their clonal expansion phase (24–28). Such a partially programmed response would make it easier to explain that the initial expansion hardly depends on the precursor density. Given that estimating just five parameters from the current data already gives rise to large confidence intervals (see Table I) illustrates that the data do not allow us to reliably fit more complicated models to the data. It would be interesting to extend our model with a more mechanistic description of the increased loss of pMHC generated by cognate T cells, if data on the loss rate of pMHC in the system become available.
The saturation function fitted by the model on a logarithmic (A) and a linear scale (B). Panels depict with parameters h = 1307 cells and n = 0.31.
Because the two very different models behave so similarly, we need new experimental data to establish which of the interpretations of the data are correct. Previously, several authors have suggested that T cells compete for access to pMHC on APCs (4, 17, 29–34), and that T cells pick up pMHC from APCs (6, 9, 11–14). However, to test which mechanism is responsible for the reduction in T cell expansion over more than four orders of magnitude in precursor densities, it would be best to directly test the proposed mechanisms. This seems possible for the new interpretation involving the increased loss of cognate pMHC on APC at high densities of T cells, because pMHC densities can, in principle, be measured directly ex vivo by activating cognate T cells with APC taken out at various time points (35, 36). Directly testing the maturation and feedback model is more difficult because which T cell subsets are responsible for the negative feedback is not known (8), nor are the molecular interactions underlying the proposed feedback.
Disclosures
The authors have no financial conflicts of interest.
Acknowledgments
Most of this article was written at the Santa Fe Institute and it was finished at the Kavli Institute for Theoretical Physics at University of California Santa Barbara (Sanata Barbara, CA).
Footnotes
This work was supported by U.S. Department of Energy Contract DE-AC52-06NA25396; National Institutes of Health Grants AI028433, OD011095, P01-AI071195, and P20-GM103452, and Contract HHSN272201000055C (to A.S.P.); and National Science Foundation Grant NSF PHY11-25915.
Abbreviation used in this article:
- pMHC
- peptide-MHC.
- Received January 3, 2013.
- Accepted March 15, 2013.