help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wofsy, C.
Right arrow Articles by Goldstein, B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Wofsy, C.
Right arrow Articles by Goldstein, B.

Biophys J, February 2001, p. 606-612, Vol. 80, No. 2

Calculations Show Substantial Serial Engagement of T Cell Receptors

Carla Wofsy,* Daniel Coombs,dagger and Byron GoldsteinDagger

 *Department of Mathematics and Statistics, University of New Mexico, Albuquerque, New Mexico 87131;  dagger Program in Applied Mathematics, University of Arizona, Tucson, Arizona 85721; and  Dagger Theoretical Biology and Biophysics Group, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 USA


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

The serial engagement model provides an attractive and plausible explanation for how a typical antigen presenting cell, exhibiting a low density of peptides recognized by a T cell, can initiate T cell responses. If a single peptide displayed by a major histocompatibility complex (MHC) can bind, sequentially, to different T cell receptors (TCR), then a few peptides can activate many receptors. To date, arguments supporting and questioning the prevalence of serial engagement have centered on the down-regulation of TCR after contact of T cells with antigen presenting cells. Recently, the existence of serial engagement has been challenged by the demonstration that engagement of TCR can down-regulate nonengaged bystander TCR. Here we show that for binding and dissociation rates that characterize interactions between T cell receptors and peptide-MHC, substantial serial engagement occurs. The result is independent of mechanisms and measurements of receptor down-regulation. The conclusion that single peptide-MHC engage many TCR, before diffusing out of the contact region between the antigen-presenting cell and the T cell, is based on a general first passage time calculation for a particle alternating between states in which different diffusion coefficients govern its transport.


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

The activation of a T cell begins with the formation of an "immunological synapse" (Shaw and Dustin, 1997; Grakoui et al., 1999), a region of contact between a T cell and an antigen presenting cell (APC). Complementary adhesion molecules first form bonds between the T cell and APC, holding the two surfaces that constitute the contact area within about 40 nm of each other. These adhesion bonds rapidly migrate to the outer region of the contact area. The surface of the APC also displays a heterogeneous population of peptide-major histocompatibility complexes (MHC). If activation is to proceed, the homogeneous population of T cell receptors (TCR) must bind to a subpopulation of the peptides on the APC. When TCR bind to peptide-MHC, this brings the surfaces closer, within about 15 nm. The TCR-peptide-MHC bonds are found, after a short transition, in the inner region of the contact area, surrounded by a ring of adhesion bonds.

The TCR-peptide bond is weak, characterized by a rapid dissociation rate constant and a half-life in the range of 1 to 20 s (reviewed in Davis et al., 1998). Recruitment of CD4 or CD8 coreceptors on T cells that bind to MHC molecules helps stabilize the TCR-peptide-MHC bond and increase the half-life (Garcia et al., 1996), but even T cells deficient in coreceptors or with blocked coreceptors can be activated (Hampl et al., 1997). T cell activation has been triggered with APCs having fewer than 100 specific peptide-MHC on their surfaces. Further, over a period of hours during which the immunological synapse is maintained, thousands of TCR are internalized. This observation led Valitutti et al. (1995) to propose that a single peptide-MHC can interact with many TCR, often causing the TCR to be internalized. They estimated that when there was a low density of specific peptide-MHC on the APC, a single peptide could serially engage 200 TCR, i.e., the ratio of TCR internalized to peptide-MHC on an APC was approximately 200. In different experiments, Itoh et al. (1999) observed down-regulation of 80 to 100 TCR per peptide-MHC. However, serial engagement has been challenged by San José et al. (2000), who showed that TCR that have not been engaged by peptide-MHC can be internalized in a peptide-MHC-dependent manner. This observation raises the possibility that the internalization of more TCR than peptides on the APC may be due to a bystander effect rather than serial engagement.

We present here a way to decide whether serial engagement occurs under physiological conditions that are independent of TCR internalization. Our approach is to determine how serial engagement depends on the parameters governing the T cell-APC interaction, i.e., the forward and reverse rate constants for the TCR-peptide bond, the surface densities of peptide-MHC and TCR, the radius of the contact area, and the diffusion coefficients of the TCR, peptide-MHC, and the bound complex that they form with each other. We derive analytic expressions, in terms of the fundamental parameters, for the total number of TCR bonds formed, serially, by a single peptide-MHC, before the complex diffuses out of the contact area, and for the rate of TCR bond formation per peptide-MHC. Applying these expressions to data, we obtain bounds on the number of TCR encounters per peptide-MHC, for peptides with distinct binding properties and signaling behavior.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Mean residence time for a particle alternating between two diffusing states

The analysis starts with a model in which a particle diffuses in a circular region of radius a, alternating between state 1, where it diffuses with diffusion coefficient D1, and state 2, where its diffusion coefficient is D2 (Fig. 1). The transition from state 1 to state 2 occurs at a rate lambda 1 and the transition from state 2 to state 1 occurs at a rate lambda 2. We calculate the mean residence time for the particle in the region, or equivalently, the mean time until the particle crosses the boundary of the region. For the case where the diffusing particle is a peptide-MHC in an immunological synapse, alternating between unbound and TCR-bound states, the calculation gives the mean time the complex remains in the contact area between the APC and the T cell.



View larger version (65K):
[in this window]
[in a new window]
 
FIGURE 1   Peptide-MHC diffuse on an antigen-presenting cell (APC) and bind to T cell receptors (TCR) on a T cell, in the region of contact of the two cells. In the model used to calculate the number and rate of serial TCR encounters by a peptide-MHC, the contact area is taken to be circular, with radius a. The initial position of the complex is specified by the distance r from the center of the region. The diffusion coefficient governing movement of the peptide-MHC is D1 when the complex is not bound to a TCR (pathway 1 to 2; then binding occurs at position 2). The diffusion coefficient is D2 when the complex is TCR-bound (pathway 2 to 3 to 4; dissociation occurs at position 4).

In the Appendix, we show that the mean residence time, < t1> , for a particle starting in state 1 at a random location within the circular region, is given by
⟨t<SUB>1</SUB>⟩=<FENCE><FR><NU>a<SUP>2</SUP></NU><DE>8</DE></FR></FENCE><FENCE><FR><NU>&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB></NU><DE>&lgr;<SUB>1</SUB>D<SUB>2</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB></DE></FR></FENCE> (1)

+<FR><NU>&lgr;<SUB>1</SUB>D<SUB>2</SUB>(D<SUB>2</SUB>−D<SUB>1</SUB>)</NU><DE>(&lgr;<SUB>1</SUB>D<SUB>2</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB>)<SUP>2</SUP></DE></FR> <FENCE>1−<FR><NU>2I<SUB>1</SUB>(&agr;a)</NU><DE>&agr;aI<SUB>0</SUB>(&agr;a)</DE></FR></FENCE>
where I0 and I1 are modified Bessel functions (Abramowitz and Stegun, 1964) and alpha  = ((lambda 1/D1) + (lambda 2/D2))1/2.

Application to free and TCR-bound peptide-MHC

For a peptide-MHC that is either free (state 1) or bound to a TCR (state 2), the transition rates are
&lgr;<SUB>1</SUB>=<A><AC>k</AC><AC>&cjs1171;</AC></A><SUB><UP>on</UP></SUB>T &lgr;<SUB>2</SUB>=k<SUB><UP>off</UP></SUB> (2)
where <A><AC>k</AC><AC>&cjs1171;</AC></A>on is the effective two-dimensional forward rate constant (Dustin et al., 1996; Shaw and Dustin, 1997) for the binding of peptide-MHC to TCR, within the contact area between an APC and T cell, koff is the reverse rate constant, and T denotes the concentration of unbound TCR in the contact area. We make the approximation that T is constant and uniform in the contact area over the period of interest. The estimates we present are for low densities of peptide-MHC, so that binding does not significantly deplete the pool of unbound TCR, nor do peptide-MHC compete among themselves for unbound TCR. Further, the estimates hold only for short times, before internalization and other processes cause extensive down-regulation of TCR or spatial redistribution of TCR and peptide-MHC in the contact area.

There are three diffusion coefficients of interest, DP for the free peptide-MHC on the APC, DT for the free TCR on the T cell, and DB for the bound complex. If, when a bridge forms between a TCR and a peptide-MHC, there are no induced cytoskeletal interactions, then from the Einstein relation it follows that the diffusion coefficient of the bound complex is
D<SUB><UP>B</UP></SUB>=<FR><NU>D<SUB><UP>P</UP></SUB>D<SUB><UP>T</UP></SUB></NU><DE>D<SUB><UP>P</UP></SUB>+D<SUB><UP>T</UP></SUB></DE></FR> (3)
From Eq. 3, we see that DB <=  DP. If the cell-cell bridge formed by a bound TCR-peptide-MHC induces cytoskeletal constraints beyond those that act separately on unbound mobile TCR and peptide-MHC, the bound complex may be essentially immobile; then DB approx  0.

In the notation of the general model, D1 = DP and D2 = DB. The expression for the mean residence time in the contact area (Eq. 1) simplifies in the two extreme cases:
⟨t<SUB>1</SUB>⟩=<FR><NU>a<SUP>2</SUP></NU><DE>8D<SUB><UP>P</UP></SUB></DE></FR> <UP>if</UP> D<SUB><UP>B</UP></SUB>=D<SUB><UP>P</UP></SUB> (4)

⟨t<SUB>1</SUB>⟩=<FENCE><FR><NU>a<SUP>2</SUP></NU><DE>8D<SUB><UP>P</UP></SUB></DE></FR></FENCE><FENCE>1+<FR><NU>&lgr;<SUB>1</SUB></NU><DE>&lgr;<SUB>2</SUB></DE></FR></FENCE>=<FR><NU>a<SUP>2</SUP></NU><DE>8D<SUB><UP>P</UP></SUB></DE></FR> (1+<A><AC>K</AC><AC>&cjs1171;</AC></A>T) <UP>if</UP> D<SUB><UP>B</UP></SUB>=0 (5)
where K = <A><AC>k</AC><AC>&cjs1171;</AC></A>on/koff is the two-dimensional equilibrium constant.

Number and rate of TCR hits per peptide-MHC

We can now estimate the number and rate of encounters ("hits") between a peptide-MHC and T cell receptors in an immunological synapse. The number of TCR hits by a single peptide-MHC, while it remains in the contact area, is approximately equal to the mean time in the contact area, < t1> , divided by the mean cycle time, i.e., the sum of the mean time for a free peptide-MHC to bind to a TCR, 1/(<A><AC>k</AC><AC>&cjs1171;</AC></A>onT), and the mean lifetime of the resulting bond, 1/koff. Then
<UP>hits</UP>≈<FR><NU>1</NU><DE>(1/(<A><AC>k</AC><AC>&cjs1171;</AC></A><SUB><UP>on</UP></SUB>T))+(1/k<SUB><UP>off</UP></SUB>)</DE></FR> ⟨t<SUB>1</SUB>⟩=<FR><NU>k<SUB><UP>off</UP></SUB><A><AC>K</AC><AC>&cjs1171;</AC></A>T</NU><DE>1+<A><AC>K</AC><AC>&cjs1171;</AC></A>T</DE></FR> ⟨t<SUB>1</SUB>⟩ (6)
The hitting rate per peptide is the number of hits divided by the time a peptide spends in the contact area
<UP>hits/s</UP>≈<FR><NU>k<SUB><UP>off</UP></SUB><A><AC>K</AC><AC>&cjs1171;</AC></A>T</NU><DE>1+<A><AC>K</AC><AC>&cjs1171;</AC></A>T</DE></FR> (7)
From Eq. 7, we see that the hitting rate depends only on the two-dimensional equilibrium constant and the reverse rate constant, or equivalently, on the forward and reverse rate constants. The hitting rate does not depend on the diffusion coefficients or the radius of the contact area. An alternative way to obtain the hitting rate is simply to note that there is one hit per mean cycle time.

Limiting cases

Two limits are helpful in understanding what determines the rate of serial engagement, for different peptides
<UP>hits/s</UP>≈<FENCE><AR><R><C><UP>k</UP><SUB><UP>off</UP></SUB></C><C><UP>if</UP></C><C><UP><A><AC>K</AC><AC>&cjs1171;</AC></A>T &z.Gt; 1</UP></C></R><R><C><OVL><UP>k</UP><SUB><UP>on</UP></SUB></OVL>T</C><C><UP>if</UP></C><C><UP><A><AC>K</AC><AC>&cjs1171;</AC></A>T &z.Lt; 1</UP></C></R></AR></FENCE> (8)
When KT 1, so that a peptide-MHC complex spends a large fraction of its time bound to TCR, the dissociation rate determines the hitting rate. In this limit, increasing the TCR concentration no longer increases the hitting rate. The TCR concentration is already sufficiently high that when a peptide-MHC complex dissociates, it will immediately find and bind to a new TCR. In the other limit, when KT 1, the rate of hitting is determined by the rate of binding. In this case, dissociation is so rapid relative to binding that the time between TCR hits is essentially the time it takes a free peptide-MHC to bind to a TCR.

If we consider peptide-MHC complexes with similar rates of binding to TCR in an immunological synapse, but with a wide range of dissociation rate constants, the hitting rates can range from near 0 (when koff is so small that KT 1 and the hitting rate is approximately koff) up to <A><AC>k</AC><AC>&cjs1171;</AC></A>onT (when koff is so large that KT 1).

Potentially, the total number of TCR hits while a peptide-MHC remains in the contact area (Eq. 6) depends on all of the system parameters, through < t1> (Eq. 1). However, if the TCR-bound complex is immobile, then < t1> is given by Eq. 5 and
<UP>hits</UP>≈<FR><NU><A><AC>k</AC><AC>&cjs1171;</AC></A><SUB><UP>on</UP></SUB>Ta<SUP>2</SUP></NU><DE>8D<SUB><UP>P</UP></SUB></DE></FR> <UP>if</UP> D<SUB><UP>B</UP></SUB>=0 (9)
Equation 9 shows that in the special case where diffusion of the bound complex is negligible, the total number of hits is independent of the reverse rate constant koff. The hitting rate, however, depends on koff.

Estimates from data

We can now use the general expression we derived for the hitting rate (Eq. 7) to see if, under typical experimental conditions, serial engagement occurs. The two-dimensional dissociation constant, KD = 1/K = koff/<A><AC>k</AC><AC>&cjs1171;</AC></A>on, has been determined to be 1 × 109 cm-2 (K = 1 × 10-9 cm2) for the binding of TCR to the peptide MCC88-103 complexed with mobile MHC molecules on a planar bilayer (Grakoui et al., 1999). For this peptide, koff = 0.057 s-1. For a T cell with 3 × 104 TCR per cell (Shaw and Dustin, 1997) and a surface area S = 5 × 10-6 cm2 (Grakoui et al., 1999), the TCR concentration at the start of the experiment T = 6 × 109 cm-2, and KT = 6. From Eq. 7, the hitting rate therefore equals 0.05 s-1.

We can estimate how many hits a single peptide makes while in the contact area by multiplying this hitting rate by the average time a peptide spends in the contact area, < t1> . The radius of the contact area is about 5 × 10-4 cm (Grakoui et al., 1999). We assume the diffusion coefficient of the MHC on an APC is typical of transmembrane proteins (reviewed by Saxton and Jacobson, 1997) and take DP = 3 × 10-10 cm2/s. For major histocompatibility antigens on various cell lines DP approx  1 - 4 × 10-10 cm2/s (Wade et al., 1989; Edidin et al., 1991; Qui et al., 1996; Munnelly et al., 2000). From Eqs. 4 and 5 we know that a2/(8DP<=  < t1>  <=  (1 + KT)a2/(8DP) and for the case being considered, 100 s <=  < t1>  <=  700 s. (The lower bound corresponds to both the free and bound peptide diffusing at the same rate, and the upper bound corresponds to the bound complex being immobile.) Thus, we estimate that at the start of an experiment, when APC and T cell have attached and a contact area has formed, a peptide initially in the contact area will engage 5 to 35 TCR in the period of 100 to 700 s (2 to 12 min) before leaving the contact area. Over a 5-h period, as in the experiments of Valitutti et al. (1995), a peptide will enter and leave the contact area numerous times, engaging additional TCR. For the set of parameters we have considered, we therefore predict that a significant number of serial engagements occur. (From the results we have derived, we cannot estimate how many serial encounters occur over long periods of time, since T will decrease with time and become nonuniform as TCR internalization occurs, and the mean time a peptide-MHC spends in the contact area will change as complexes enter the region from the periphery.)

Effect of internalization

Internalization of TCR, and other processes that may contribute to TCR down-regulation on a time scale consistent with published data (Valitutti et al., 1995; San José et al., 2000; Niedergang et al., 1997) do not affect our estimates of the extent of serial engagement in the initial period of cell-cell contact. It is easy to apply our general model to the case where bound TCR are subject to internalization. In this case, there are two ways for a peptide-MHC to make a transition from state 2 (bound) to state 1 (free), i.e., by internalization of the bound TCR (at rate lambda ) or by dissociation of the bond between the TCR and the peptide-MHC (at rate koff). Then lambda 2 = koff + lambda  and the hitting rate is
<UP>hits/s</UP>≈<FR><NU>1</NU><DE>(1/(<A><AC>k</AC><AC>&cjs1171;</AC></A><SUB><UP>on</UP></SUB>T))+(1/(k<SUB><UP>off</UP></SUB>+&lgr;))</DE></FR>=<FR><NU>k<SUB><UP>off</UP></SUB><A><AC>K</AC><AC>&cjs1171;</AC></A>T</NU><DE>1+(<A><AC>K</AC><AC>&cjs1171;</AC></A>T/(1+(&lgr;/k<SUB><UP>off</UP></SUB>))</DE></FR> (10)
To obtain a rough estimate of lambda , in Fig. 2 we fit a simple exponential decay to the data of Valitutti et al. (1995) and obtained lambda  approx  .028/min or 4.7 × 10-4/s when the peptide density is 25 nM and lambda  approx  .038/min or 6.3 × 10-4/s when the peptide density is 20 µM. For most TCR-peptide-MHC, koff >=  0.01 s-1 (Davis et al., 1998). Then lambda /koff 1, and Eq. 10 indicates that internalization does not alter the previously estimated rate (Eq. 7) of serial TCR encounters by peptide-MHC in the immunological synapse.



View larger version (20K):
[in this window]
[in a new window]
 
FIGURE 2   The data, from Fig. 1b of Valitutti et al. (1995), show the time course of CD3 down-regulation when T cells are conjugated with peptide-pulsed APC at peptide concentrations of 25 nM () and 20 µM (), corresponding to approximately 50 and 1500 peptides per APC (see Valitutti et al., 1995, for details). The solid lines are separate nonlinear least squares fits to these data of the function: exp(-lambda t) + fd, where lambda  is the apparent internalization rate and fd is the fraction of TCR remaining after down-regulation is complete. The best fit values of the parameters are: lambda  = 0.028 ± 0.005, fd = 0.50 ± 0.02 (); lambda  = 0.038 ± 0.003, fd = 0.14 ± 0.02 ().

Serial engagement of TCR by peptide agonists, weak agonists, and antagonists

To estimate the hitting rate (Eq. 7), we must know koff and the two-dimensional equilibrium binding constant, K. As discussed, for the one peptide (MC88-103) TCR interaction for which K has been determined (Grakoui et al., 1999), K = 1 × 10-9 cm2, koff = 0.057 s-1 and, therefore, <A><AC>k</AC><AC>&cjs1171;</AC></A>on = 5.7 × 10-11 cm2/s. Grakoui et al. (1999) determined from a biosensor study that the three-dimensional forward rate constant kon = 900 M-1 s-1.

If we assume that for all the peptides they studied, <A><AC>k</AC><AC>&cjs1171;</AC></A>on is proportional to kon, we can estimate hitting rates and upper and lower bounds on the number of peptide-TCR serial encounters. This is done in Table 1, where we see that serial engagement is predicted to occur for agonists, weak agonists, and antagonists.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 1   Estimates of hitting rates and initial number of hits for peptide-MHC interacting with TCR

In reaching this conclusion, we assumed that the mean time for a peptide-MHC to dissociate from a TCR is the same as that obtained from measurements where one of the reactants is in solution, such as in a BIACORE experiment. This assumption will break down if there is a significant probability that when the bond between a TCR and a peptide-MHC dissociates, the peptide-MHC will rebind to the same TCR. Such rebinding will increase the effective lifetime of the bond, i.e., reduce the value of koff, and, from Eq. 7, reduce the hitting rate. Dustin (1997) investigated this rebinding question for the binding of the glycoprotein CD2 on T cells to CD58 on glass-supported planar bilayers and showed that dissociation of CD2-CD58 bonds led to the creation of new partners rather than reformation of the same pairs. Apparently, the relative diffusion between a CD2 and CD58 pair upon dissociation was such that the two molecules became well separated. Is this true as well for TCR/peptide-MHC pairs in the contact region? In the experiments of Grakoui et al. (1999), 80% of the TCR population appeared to be immobile in a fluorescence photobleaching recovery experiment. Let us estimate what the diffusion coefficient of the peptide-MHC, Dp, must be to achieve TCR/peptide-MHC separation after dissociation, in the extreme case when all TCR are immobile. Upon dissociation, in a time t the peptide-MHC diffuses a mean square distance r2 = 4Dpt. Assuming re-formation of the same pair does not occur, the peptide-MHC will diffuse a mean time Delta t = 1/(<A><AC>k</AC><AC>&cjs1171;</AC></A>onT) before binding another TCR. (Recall that <A><AC>k</AC><AC>&cjs1171;</AC></A>on is the effective two-dimensional forward rate constant and T is the concentration of unbound TCR in the contact area.) If the peptide-MHC diffuses a distance that is less than the average distance between TCR, then rebinding to the same TCR may become significant. For randomly distributed TCR on a surface, the mean square distance between TCR is < s2>  = 1/(pi T). Thus we expect rebinding of the same pairs to be negligible and koff in the contact area to be the intrinsic off rate constant when Dp > < s2> /(4Delta t) = <A><AC>k</AC><AC>&cjs1171;</AC></A>on/(4pi ). The peptide in Table 1 with the largest forward rate constant, <A><AC>k</AC><AC>&cjs1171;</AC></A>on = 9.7 × 10-10 cm2/s, is the weak agonist N72T. The inequality predicts that if Dp > 8 × 10-11 cm2/s, re-formation of the same TCR-peptide-MHC bond will be negligible and koff in the contact region will be the same as determined from solution measurements. Although the measured value, DP approx  - 4 × 10-10 cm2/s (Wade et al., 1989; Edidin et al., 1991; Qui et al., 1996; Munnelly et al., 2000), is close to the predicted value for re-formation of a TCR-peptide-MHC bond, indirect evidence suggests that koff is not reduced in the contact area. The off rate constant or, equivalently, the mean lifetime of the TCR-peptide MHC bond, is the dominant factor in determining whether a peptide is an agonist, weak agonist, or antagonist (Matsui et al., 1994; Lyons et al., 1996; Kersh et al., 1998). N72T is a weak agonist with an intrinsic koff = 0.14 s-1. If koff was reduced by a factor of just 2 or 3 in the contact area, we would expect N72T to be a full agonist (see Table 1). Since N72T is not a full agonist, we conclude that its intrinsic koff characterizes bond dissociation in the contact area as well as in solution.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

We have derived an expression, Eq. 6, for the number of TCR a peptide-MHC binds, sequentially, before diffusing out of the immunological synapse, assuming the peptide-MHC starts from a random position within the contact area. From this expression we obtained the per peptide hitting rate, Eq. 7. We also obtained expressions for the shortest and longest times a peptide-MHC spends in the contact area, Eqs. 4 and 5. These times correspond to the two extremes, where the motion of the peptide-MHC is not influenced by binding to the TCR and where the peptide-MHC becomes immobile upon forming a bond. We used these results to see if serial engagement occurs at the start of an experiment, when an APC and T cell have come into contact, and the peptide-MHC and TCR concentrations are still uniform in the contact area. We estimated, for each of a series of peptides, lower and upper bounds on the number of hits (TCR engagements) a peptide-MHC makes before it leaves the contact area (Table 1).

The estimates of serial engagement rates, per peptide-MHC, depend on the TCR concentration being approximately uniform in the contact area during the first few minutes of the experiment but are independent of how TCR are transported on the T cell surface. The predictions are the same whether active (Valitutti et al., 1995; Wülfing and Davis, 1998) or passive (diffusive) transport mechanisms dominate the movement of TCR into the immunological synapse. The estimates show that for agonists (i.e., peptides that trigger T cell responses), weak agonists, and antagonists, serial engagement occurs.

Previous arguments supporting serial engagement (Valitutti et al., 1995; Itoh et al., 1999; Lanzavecchia and Sallusto, 2000) have been based on the observation of extensive TCR down-regulation, apparently triggered by a relatively small number of peptide-MHC. If bystander effects lead to the internalization of more than one TCR per peptide-TCR encounter (San José et al., 2000; Niedergang et al., 1997), published estimates of the number of serial engagements per peptide (Valitutti et al., 1995) may be high. Nevertheless, for the parameters that characterize peptide-TCR interactions, serial engagement is expected to be robust.

If all peptide-MHC that interact with TCR, whether agonists, weak agonists or antagonists, undergo serial engagement, what role does serial engagement play in T cell activation? The simple formation of bonds between TCR and peptide-MHC within the immunological synapse is probably not sufficient to initiate a T cell response. As with all other multisubunit immune recognition receptors (MIRR), there is considerable evidence that receptor aggregation must occur before a TCR can initiate a signaling cascade (Boniface et al., 1998; Bachmann et al., 1998; Bachmann and Ohashi, 1999). Oligomerization of TCR bound to peptide-MHC is followed rapidly by a series of TCR modifications. First, specific tyrosines residing in immunoreceptor tyrosine-based activation motifs (ITAMs) on subunits of the TCR-CD3 complex are phosphorylated by the protein tyrosine kinase (PTK) Lck. This results in recruitment of the PTK ZAP-70 from the cytosol to the phosphorylated TCR zeta  chain. Zap-70 in turn becomes phosphorylated and, as long as the signaling complex remains intact, the signaling cascade proceeds (reviewed in Germain and Stefanova, 1999; Lanzavecchia et al., 1999).

The mean lifetime of the TCR-peptide-MHC bond (1/koff) appears to be the dominant factor in determining whether a peptide is an agonist, weak agonist, or antagonist (Matsui et al., 1994; Lyons et al., 1996; Kersh et al., 1998). From Table 1 we see that there is a strong correlation between koff, the hitting rate, and a peptide's ability to activate T cells. In Table 1, agonists have the smallest koff values and slowest hitting rates and the antagonist has the highest koff and the fastest hitting rate. However, a second mechanism, kinetic proofreading, has been invoked to explain the correlation between T cell activation and the lifetime of the TCR-peptide-MHC bond (McKeithan, 1995). The kinetic proofreading model postulates that for a particular cellular response to occur, the TCR must complete a series of modifications (e.g., phosphorylations, associations with enzymes, adaptors). If the lifetime of the TCR-peptide-MHC bond is too short, then almost always, a bound TCR will dissociate from the peptide-MHC, become disengaged from any signaling molecules it has associated with, and be dephosphorylated to its basal level before it undergoes the necessary number of modifications to achieve activation. As observed, the kinetic proofreading model predicts that the best agonist peptides will have the longest TCR-peptide-MHC bond lifetimes.

If kinetic proofreading is all that matters, then the longer the lifetime (the smaller the koff) of the TCR-peptide-MHC bond, the better the peptide will be at activating T cells. If that is so, why have no agonists been reported with half-lives longer than about 30 s (Corr et al., 1994; Lyons et al., 1996; Grakoui et al., 1999; Davis et al., 1999)? Lanzavecchia et al. (1999) have argued that in addition to the lifetime of the TCR-peptide-MHC bond, of key importance in T cell activation is the number of triggered TCR. For a fixed kon, as koff decreases, the hitting rate (Eq. 7), decreases. Thus, what appears to place an upper limit on the lifetime of an agonist peptide is the balance between the lifetime of the TCR-peptide-MHC bond and the hitting rate.

We noted that in Table 1, the antagonist has the highest hitting rate. One way antagonism can occur is if a critical initiating kinase is limiting, as is the initiating kinase Lyn for the MIRR Fcepsilon RI, in rat basophilic leukemia cells (Torigoe et al., 1997; Wofsy et al., 1997). As proposed by Torigoe et al. (1998), the initiating kinase can be tied up in unproductive associations with TCR that are repeatedly forming short-lived bonds with antagonist peptide-MHC. The result is that less initiating kinase is available for association with TCR that form long-lived bonds with agonist peptide-MHC. A second way the agonist may be inhibited is through the formation of mixed oligomers composed of TCR-agonist bonds and TCR-antagonist bonds. These heterogeneous oligomers will have shorter lifetimes than oligomers formed solely by TCR-agonist bonds and be less effective at signaling (Davis et al., 1998). Of course, if the lifetime of the TCR-antagonist bond becomes too short, even the initiating kinase will not have time to associate. This cannot be compensated for by raising the hitting rate and creating more oligomers, because in the limit of large koff, the hitting rate becomes independent of koff and approaches its maximum value, <A><AC>k</AC><AC>&cjs1171;</AC></A>on T (Eq. 9). The balance between hitting rates and the TCR-peptide-MHC bond lifetime, i.e., between serial engagement and kinetic proofreading, plays a dominant role in determining the range of koff values over which peptides are active, whether as agonists, weak agonists, or antagonists.


    APPENDIX
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Derivation of mean time in the contact area

For a particle diffusing in a circular region of radius a, with an initial position r <=  a and initial state i (i = 1, 2), we define ti as the mean time to reach the boundary of the region. The mean residence times t1 and t2 satisfy the partial differential equations
D<SUB>1</SUB>∇<SUP>2</SUP>t<SUB>1</SUB>−&lgr;<SUB>1</SUB>t<SUB>1</SUB>+&lgr;<SUB>1</SUB>t<SUB>2</SUB>+1=0 (11)

D<SUB>2</SUB>∇<SUP>2</SUP>t<SUB>2</SUB>−&lgr;<SUB>2</SUB>t<SUB>2</SUB>+&lgr;<SUB>2</SUB>t<SUB>1</SUB>+1=0 (12)
where D1, D2, lambda 1, and lambda 2 are the diffusion coefficients for the two states, and rates of transition between the states, defined previously. A heuristic derivation of Eqs. 11 and 12, based on a two-dimensional random walk, is analogous to the derivation of Eqs. 10a, b presented in Appendix A of Goldstein et al. (1984). The boundary conditions for the application we consider here are ti(a) = 0 and ti(0) is finite, for i = 1, 2.

Multiplying Eq. 11 by lambda 2 and Eq. 12 by lambda 1 and adding the resulting equations gives Poisson's equation for an average t
D∇<SUP>2</SUP>t=<UP>−</UP>1 (13)
where
D=<FR><NU>&lgr;<SUB>2</SUB></NU><DE>&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB></DE></FR> D<SUB>1</SUB>+<FR><NU>&lgr;<SUB>1</SUB></NU><DE>&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB></DE></FR> D<SUB>2</SUB> (14a)

t=<FR><NU>&lgr;<SUB>2</SUB>D<SUB>1</SUB></NU><DE>&lgr;<SUB>1</SUB>D<SUB>2</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB></DE></FR> t<SUB>1</SUB>+<FR><NU>&lgr;<SUB>1</SUB>D<SUB>2</SUB></NU><DE>&lgr;<SUB>1</SUB>D<SUB>2</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB></DE></FR> t<SUB>2</SUB> (14b)
The solution to Eq. 13, subject to the boundary conditions, is t(r) = (a2 - r2)/(4D). Using the definition of t (Eq. 14b) to eliminate t2 in Eq. 11 gives this equation for t1:
D<SUB>1</SUB>∇<SUP>2</SUP>t<SUB>1</SUB>−(&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB>/D<SUB>2</SUB>)t<SUB>1</SUB>=<UP>−</UP>1 (15)

−(&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB>D<SUB>1</SUB>/D<SUB>2</SUB>)(a<SUP>2</SUP>−r<SUP>2</SUP>)/(4D)
The general solution is the following sum of the general solution to the corresponding homogeneous equation (Abramowitz and Stegun, 1964) and a particular solution to Eq. 15:
t<SUB>1</SUB>(r)=AI<SUB>0</SUB>(&agr;r)+BK<SUB>0</SUB>(&agr;r) (16)

+(a<SUP>2</SUP>−r<SUP>2</SUP>)/(4D)+<FR><NU>(1−D<SUB>1</SUB>/D)(D<SUB>2</SUB>/D)</NU><DE>&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB></DE></FR>
where I0 and K0 are modified Bessel functions and alpha  = ((lambda 1/D1) + (lambda 2/D2))1/2. The function K0 becomes infinite as r tends to 0, so to maintain a finite solution, we need B = 0. Applying the other boundary condition, t1 (a) = 0, gives
t<SUB>1</SUB>(r)=<FR><NU>a<SUP>2</SUP>−r<SUP>2</SUP></NU><DE>4D</DE></FR>+<FR><NU>(1−D<SUB>1</SUB>/D)(D<SUB>2</SUB>/D)</NU><DE>&lgr;<SUB>1</SUB>+&lgr;<SUB>2</SUB></DE></FR><FENCE>1−<FR><NU>I<SUB>0</SUB>(&agr;r)</NU><DE>I<SUB>0</SUB>(&agr;a)</DE></FR></FENCE> (17)
Averaging over r, we obtain Eq. 1.

    ACKNOWLEDGMENTS

This work was supported by National Institutes of Health grant GM35556 and National Science Foundation grant MCB9723897 and performed in part under the auspices of the U.S. Department of Energy.

    FOOTNOTES

Received for publication 15 June 2000 and in final form 21 November 2000.

Address reprint requests to Dr. Byron Goldstein, Theoretical Biology and Biophysics Group, Theoretical Division, T-10, MS K710, Los Alamos National Laboratory, Los Alamos, NM 87545. Tel.: 505-667-6538; Fax: 505-665-3493; E-mail: bxg{at}lanl.gov.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
RESULTS
DISCUSSION
APPENDIX
REFERENCES

Biophys J, February 2001, p. 606-612, Vol. 80, No. 2
© 2001 by the Biophysical Society   0006-3495/01/02/606/07  $2.00



This article has been cited by other articles:


Home page
Biophys. JHome page
O. Dushek and D. Coombs
Analysis of Serial Engagement and Peptide-MHC Transport in T Cell Receptor Microclusters
Biophys. J., May 1, 2008; 94(9): 3447 - 3460.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
S. N. Arkhipov and I. V. Maly
Quantitative Analysis of the Role of Receptor Recycling in T Cell Polarization
Biophys. J., December 1, 2006; 91(11): 4306 - 4316.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
J. R. Wedagedera and N. J. Burroughs
T-Cell Activation: A Queuing Theory Analysis at Low Agonist Density
Biophys. J., September 1, 2006; 91(5): 1604 - 1618.
[Abstract] [Full Text] [PDF]


Home page
J. Immunol.Home page
A. Jansson, E. Barnes, P. Klenerman, M. Harlen, P. Sorensen, S. J. Davis, and P. Nilsson
A Theoretical Framework for Quantitative Analysis of the Molecular Basis of Costimulation
J. Immunol., August 1, 2005; 175(3): 1575 - 1585.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
C. Utzny, M. Faroudi, and S. Valitutti
Frequency Encoding of T-Cell Receptor Engagement Dynamics in Calcium Time Series
Biophys. J., January 1, 2005; 88(1): 1 - 14.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
D. Coombs and B. Goldstein
Effects of the Geometry of the Immunological Synapse on the Delivery of Effector Molecules
Biophys. J., October 1, 2004; 87(4): 2215 - 2220.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
D. Coombs, M. Dembo, C. Wofsy, and B. Goldstein
Equilibrium Thermodynamics of Cell-Cell Adhesion Mediated by Multiple Ligand-Receptor Pairs
Biophys. J., March 1, 2004; 86(3): 1408 - 1423.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
S.-J. E. Lee, Y. Hori, and A. K. Chakraborty
Low T cell receptor expression and thermal fluctuations contribute to formation of dynamic multifocal synapses in thymocytes
PNAS, April 15, 2003; 100(8): 4383 - 4388.
[Abstract] [Full Text] [PDF]


Home page
J. Immunol.Home page
J. R. Faeder, W. S. Hlavacek, I. Reischl, M. L. Blinov, H. Metzger, A. Redondo, C. Wofsy, and B. Goldstein
Investigation of Early Events in Fc{epsilon}RI-Mediated Signaling Using a Detailed Mathematical Model
J. Immunol., April 1, 2003; 170(7): 3769 - 3781.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Wofsy, C.
Right arrow Articles by Goldstein, B.