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 Anderson, T. G.
Right arrow Articles by McConnell, H. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Anderson, T. G.
Right arrow Articles by McConnell, H. M.

Biophys J, November 1999, p. 2451-2461, Vol. 77, No. 5

Interpretation of Biphasic Dissociation Kinetics for Isomeric Class II Major Histocompatibility Complex-Peptide Complexes

Thomas G. Anderson and Harden M. McConnell *

Department of Chemistry, Stanford University, Stanford, California 94305 USA

    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

Antigenic peptides bound to class II major histocompatibility complex (MHC) proteins play a key role in the distinction between "self" and "nonself" by the cellular immune system. Although the formation and dissociation of these complexes are often thought of in terms of the simple mechanism MHC + P right-left-harpoons  MHC-P, studies of MHC-peptide dissociation kinetics suggest that multiple interconverting forms of the bound MHC-peptide complex can be formed. However, the precise relationship between observed dissociation data and proposed multiple-complex mechanisms has not been systematically examined. Here we provide a mathematical analysis to fill this gap and attempt to clarify the kinetic behavior that is expected to result from the proposed mechanisms. We also examine multiple-complex dynamics that can be "hidden" in conventional experiments. Although we focus on MHC-peptide interactions, the analysis provided here is fully general and applies to any ligand-receptor system having two distinct bound states.

    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

Proteins of the class II major histocompatibility complex (MHC) are expressed on the surface of antigen-presenting cells, where they display peptide fragments to receptors on CD4+ T helper cells. Like antibodies, these T-cell receptors (TCRs) exhibit great diversity at their recognition site (Davis and Bjorkman, 1988; Davis, 1990; Garboczi et al., 1996; Garcia et al., 1996; Jardetzky, 1997), enabling them to distinguish between different peptides in the context of the MHC protein, effecting the specific cellular immune response (reviewed in Abbas et al., 1994). It is generally assumed that the immunologically active form of an MHC-peptide complex adopts a single unique structure that can be recognized by the TCR. However, a growing body of evidence suggests that multiple isomeric complexes can be formed in MHC-peptide systems (Dornmair et al., 1989; Sadegh-Nasseri and McConnell, 1989; Beeson and McConnell, 1994; Beeson et al., 1996; Schmitt et al., 1998a,b). In several cases, these isomeric complexes can be distinguished by T cells (Viner et al., 1996; Rabinowitz et al., 1997).

The structural changes involved in the formation of isomeric complexes are not well understood. Although the static structures of several MHC-peptide complexes from mice as well as humans have been characterized by x-ray crystallography (Brown et al., 1993; Fremont et al., 1996, 1998; Jardetzky et al., 1996; Scott et al., 1998), such structures provide little information about the dynamics of the complexes. Rather than rely upon structural data alone, this laboratory has investigated MHC-peptide isomers indirectly by studying the kinetics of formation and dissociation of MHC-peptide complexes. Kinetic experiments have been instrumental in elucidating the mechanism of MHC-peptide binding (Fig. 1 A), which is thought to involve numerous reactions. Initially present complexes of MHC with endogenous peptides (MPe) can dissociate to form empty MHC (M) (Witt and McConnell, 1991, 1992). This "active" unbound MHC can reversibly convert to a "dormant" form (MI) (Rabinowitz et al., 1998; Natarajan et al., 1999) or be irreversibly inactivated (forming MX) (Mason and McConnell, 1994). We shall not discuss these steps in detail here. In studies of binding to labeled peptides (P*), it has been shown in several cases that MHC can form multiple isomeric complexes ({MP*}1 and {MP*}2) (Sadegh-Nasseri and McConnell, 1989; Witt and McConnell, 1992; Sadegh-Nasseri et al., 1994). Reactions of these complexes, particularly interconversion and dissociation, will be the focus of our discussion. Additional reactions not shown in Fig. 1 have also been proposed, for example, the displacement of one peptide by another via a transient two-peptide intermediate (Tampé and McConnell, 1991; de Kroon and McConnell, 1993, 1994; Witt and McConnell, 1994). We will not discuss these ancillary reactions further, restricting our discussion to the subset of reactions shown in Fig. 1 B.



View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 1   Dynamics of the peptide-MHC system. (A) Reactions that have been inferred from kinetics experiments to occur in MHC-peptide interaction. The species involved are M, "active" unbound MHC; P*, labeled peptide; MI, "dormant" inactive MHC; MX, permanently inactivated MHC; {MPe}, MHC bound to endogenous peptide; {MP*Pe}, ternary complex of labeled and endogenous peptides with MHC; {MP*}1 and {MP*}2, isomeric forms of the labeled peptide-MHC complex. (B) The subset of reactions involved in the dissociation of the MHC-peptide complexes.

To assess the biological significance of isomeric MHC-peptide complexes, it is important to know how and under what conditions the different isomers are formed. If one isomer is formed in vanishingly small amounts, or if the two complexes interconvert very rapidly, the complexes might not be distinguished by T cells. However, if isomeric MHC-peptide complexes interconvert slowly, they might give rise to distinct immune responses because of the short (~10 s) duration of MHC-peptide/TCR interactions under physiological conditions (Matsui et al., 1991, 1994). This could be particularly important if different MHC-peptide isomers are preferentially formed in different in vivo environments. To address these issues, it is necessary to solve the biophysical problem of determining the proportions of the isomeric complexes that are present under various conditions, as well as the rates at which they form and interconvert.

Unfortunately, the MHC-peptide system presents a number of obstacles to conventional kinetic analysis. Side reactions of the empty MHC molecule and formation of empty MHC by the dissociation of prebound endogenous peptides make it difficult to interpret measurements of the binding rates and equilibrium binding constants of labeled peptides. Consequently, in the analysis presented here we restrict our discussion to MHC-peptide dissociation reactions (Fig. 1 B). Although multiple forms of peptide-MHC complexes are known (Boniface et al., 1996; Dadaglio et al., 1997; Runnels et al., 1996), isomeric complexes of the labeled MHC-peptide complex are not directly observed in kinetic experiments, so the relative concentrations of these complexes are generally unknown. The absence of well-defined initial conditions greatly complicates the analysis of kinetic data from MHC-peptide reactions (Steinfeld et al., 1989).

The relationships between the observed dissociation rate constants and the microscopic rate constants are complicated for all but the simplest cases. To illustrate our analysis of MHC-peptide dissociation kinetics, we shall discuss kinetic data for a simulated two-complex MHC-peptide system with arbitrarily chosen microscopic rate constants. In addition, we will briefly apply the analysis to previously published kinetic experiments to demonstrate the application of this analysis to real experimental data. Although we present our discussion in terms of MHC-peptide interactions, the analysis provided here is fully general and applies to dissociation studies of any ligand-receptor system having two kinetically distinct bound states.

    EXPERIMENTAL TECHNIQUES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

Preparation and dissociation of MHC-peptide complexes

The preparation of MHC-peptide complexes involves incubating MHC proteins with an excess of peptide that is labeled with a radioactive tag (Sadegh-Nasseri et al., 1994), a fluorescent moiety (Tampé and McConnell, 1991; Witt and McConnell, 1994), or some other detectable group (Jensen, 1992). This incubation is typically carried out under quasiphysiological conditions of pH 5.3, 150 mM sodium chloride, and 37°C, to mimic the conditions under which MHC loading occurs within the endosomal compartments of antigen-presenting cells (Tulp et al., 1994). After a period of incubation, which generally lasts 10-20 h, the complex formed is separated from unbound peptide, using a size exclusion column (Witt and McConnell, 1991; de Kroon and McConnell, 1993). Samples of the incubation mixture may be taken at various time points during the course of the binding reaction; the amount of labeled peptide bound to the MHC can then be measured to provide a profile of the binding kinetics (Witt and McConnell, 1991; de Kroon and McConnell, 1993; Liang et al., 1995).

In a dissociation experiment, labeled MHC-peptide complex is prepared and isolated as described above. The complex is then incubated under the dissociation conditions of interest, which may or may not be the same as the binding conditions. The amount of labeled complex is then measured over time as the labeled peptide dissociates from the MHC. An unlabeled competitor peptide is sometimes added to the incubation to inhibit rebinding of dissociated peptide. The collected data are normalized and fit to either a mono- or biexponential decay curve (Witt and McConnell, 1994), as illustrated by Eqs. 1 and 2:
<FR><NU>[{<UP>MP*</UP>}](t)</NU><DE>[{<UP>MP*</UP>}]<SUB>0</SUB></DE></FR>=e<SUP><UP>−k</UP><SUP><UP>o</UP></SUP><SUB><UP>off</UP></SUB><UP>t</UP></SUP> (1)

<FR><NU>[{<UP>MP*</UP>}](t)</NU><DE>[{<UP>MP*</UP>}]<SUB>0</SUB></DE></FR>=F<SUP><UP>o</UP></SUP>e<SUP><UP>−k</UP><SUP><UP>o</UP></SUP><SUB><UP>fast</UP></SUB><UP>t</UP></SUP>+S<SUP><UP>o</UP></SUP>e<SUP><UP>−k</UP><SUP><UP>o</UP></SUP><SUB><UP>slow</UP></SUB><UP>t</UP></SUP> (2)
The superscript o indicates the observed kinetic parameters. For the biphasic curve, the magnitudes of the fast and slow exponential phases are denoted by Fo and So, respectively. For the biphasic dissociation curve, Fo + So = 1. Such a dissociation curve is shown in Fig. 2.



View larger version (16K):
[in this window]
[in a new window]
 
FIGURE 2   A typical biphasic dissociation curve for the simulated MHC-peptide system discussed in the text. The time dependence of the total concentration of complex is fit by [{MP*}]/[{MP*}]0 = 0.18e-1.2t + 0.82e-0.04t. The observed parameters are the magnitude of fast phase, Fo = 0.18; the rate constant of the fast phase, kfasto = 1.2 h-1; and the rate constant of the slow phase, kslowo = 0.04 h-1.

Numerical simulations of experimental data

Dissociation curves and other experimental data for the hypothetical MHC-peptide system discussed in this paper were generated by the numerical evaluation of rate equations. Differential equations describing the rates of reaction for the species shown in Fig. 1 A were integrated numerically using Mathematica 3.0 (Wolfram Research) to determine the concentrations of the species as functions of time. Magnitudes of the fast and slow exponential phases of dissociation were taken from fits of the biexponential equation (Eq. 2) to the simulated dissociation data. Except as otherwise noted, values for the rate constants and initial concentrations were taken from previous kinetic simulations (Beeson et al., 1996; Rabinowitz et al., 1998). Rate constants: dissociation of endogenous peptide, ke = 2 h-1; reversible inactivation of MHC, kai = 15 h-1; activation of dormant MHC, kia = 0.25 h-1; irreversible inactivation of MHC, kx = 0.05 h-1; formation of {MP*}1, kon,1 = 360 µM-1 h-1; dissociation of {MP*}1, koff,1 = 0.7 h-1; conversion of {MP*}1 to {MP*}2, k12 = 0.48 h-1; conversion of {MP*}2 to {MP*}1, k21 = 0.048 h-1; formation of {MP*}2, kon,2 = 61.7 µM-1 h-1; dissociation of {MP*}2, koff,2 = 0.012 h-1. Initial concentrations used for simulation of binding reactions: endogenous complex, [{MPe}]0 = 1 µM; all other MHC species = 0; labeled peptide, [P*]0 = 100 µM. For dissociation reactions, the initial concentration of labeled peptide was set at zero, and the initial concentration of unlabeled competitor peptide was set at 100 µM.

    RESULTS AND DISCUSSION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

Ambiguous dissociation kinetics

Fig. 2 shows a simulated dissociation curve for a complex of MHC with peptide P*. The distinct biexponential shape of the curve is not consistent with the simple reaction MP* right-arrow M + P*. Two distinct MHC-peptide complexes must be involved in the dissociation, as illustrated by Fig. 1 B (Sadegh-Nasseri and McConnell, 1989). If the two complexes were distinguishable by spectroscopic or other means, their individual concentrations could be measured over time to provide a complete picture of the kinetics of the system (Fig. 3). However, for kinetic studies of MHC-peptide systems, isomeric complexes of the MHC-peptide complex generally appear as a single signal (see Experimental Techniques), providing no information about the individual concentrations of the complexes. As a consequence, many different MHC-peptide systems with disparate microscopic behaviors could give rise to the same observed dissociation curve (Fig. 3).



View larger version (13K):
[in this window]
[in a new window]
 
FIGURE 3   Many different two-complex MHC-peptides produce the observed dissociation curve shown in Fig. 2. Shown are dissociation curves for (A) Scheme 1, (B) Scheme 2, (C) Scheme 3, and (D) Scheme 4, discussed in the text. The behavior of Scheme 5 is very similar to that of B. Dashed lines indicate the concentrations of complexes {MP*}1 (- - -) and {MP*}2 (--- - ---); the total concentration of complex is indicated by the solid line. Note that although the behavior of the individual complexes differs from one reaction scheme to another, the observed total concentration is identical for all of them.

There are four microscopic rate constants involved in the dissociation of a two-complex MHC-peptide system (Fig. 1 B), yet there are only two macroscopic rate constants observed (Fig. 2). Without concentration data, there are two degrees of freedom in the "solution space" of systems that are consistent with our observed dissociation curve (Fig. 4). We want to determine where our system lies in this solution space. That is, we wish to know the microscopic kinetic parameters (rate constants and initial concentrations) of the system that produced the dissociation curve shown in Fig. 2. Before addressing this issue, we shall first examine several "candidate" systems that are consistent with these dissociation data.



View larger version (52K):
[in this window]
[in a new window]
 
FIGURE 4   The set of systems that produce the dissociation curve shown in Fig. 2 lies on a two-dimensional surface in "kinetic space." The three axes correspond to the values of the microscopic rate constants k12, k21, and koff,2; the fourth microscopic rate constant, koff,1, is related to the others by Eq. 20 (Table 1). The fraction of complex in the form {MP*}1 is indicated by color. Each point on the surface shown represents an MHC-peptide system that produces the observed dissociation curve; the points corresponding to Schemes 1-5 are indicated. The dotted line indicates the subset of this solution space, which is also consistent with the data shown in Fig. 4. The gray dot indicates the system that was used to numerically simulate the data in Figs. 2 and 4.

The parallel system

One candidate is the "parallel" system, with two noninterconverting complexes with the dissociation rate constants and initial populations shown in Scheme 1. 
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(18%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>∼0</SUB></LL><UL><SUB>∼0</SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(82%)</UL></LIM></C></R><R><C><SUB> 1.2 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>0.04 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 1</UP>
As shown, the rate of interconversion between the complexes in this system is negligible. This system produces the dissociation curve shown in Fig. 2. Interpreting the dissociation kinetics according to this scheme, the observed fast and slow exponential phases of the dissociation curve correspond to the dissociation of complexes {MP*}1 and {MP*}2, respectively. The microscopic observed rate constants koff,1 and koff,2 are equal to kfasto and kslowo, and the amplitudes of the fast and slow exponential phases in the dissociation curve faithfully reflect the initial populations of the two complexes. This is the simplest interpretation of the observed dissociation data.

Sequential systems

The parallel system is not the only solution consistent with our data, however. One could also account for the observed dissociation curve with a "sequential" system in which the two complexes can slowly interconvert, as shown in Scheme 2. 
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(21%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>0.04 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>«0.04 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(79%)</UL></LIM></C></R><R><C><SUB> 1.2 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>∼0</SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 2</UP>
Here, complex {MP*}2 does not dissociate directly at an appreciable rate. Rather, it dissociates through the "kinetic intermediate" {MP*}1. The initial populations of the two complexes are still close to the magnitudes of the observed fast and slow dissociation phases, and the microscopic rate constants koff,1 and k21 are equal to the observed rate constants kfasto and kslowo. Scheme 2 also produces the observed dissociation curve shown in Fig. 2.

Other schemes that are consistent with the observed dissociation curve can be proposed in which the interconversion of the complexes is more rapid. In such systems, the observed kinetics are distorted from the microscopic rates. For example, the system shown in Scheme 3 produces the MP* dissociation curve in Fig. 2, even though the microscopic rate constants do not match the observed rate constants and the initial populations of the two complexes are far from the magnitudes of the observed fast and slow dissociation phases.
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(42%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>0.08 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>0.56 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(58%)</UL></LIM></C></R><R><C><SUB> 0.60 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>∼0</SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 3</UP>
The kinetic distortions seen in this sequential scheme can be understood in terms of a steady-state approximation for the fast-dissociating complex {MP*}1. In the fast phase of dissociation, the population of {MP*}1 rapidly drops to its steady-state level, with a rate constant kfasto about equal to koff,1 + k12 = 0.60 h-1 + 0.56 h-1 = 1.16 h-1. After {MP*}1 reaches this level, the steady-state approximation tells us that the system decays with a rate constant kslowo of (k21 koff,1)/(koff,1 + k12) = 0.041 h-1. These are the values of kfasto and kslowo that are observed in Fig. 2. The steady-state approximation also helps rationalize the observed magnitude of the fast dissociation phase. Because the initial population of fast-dissociating {MP*}1 is partitioned over dissociation and interconversion reactions, the initial 42% population of {MP*}1 gives rise to only an 18% observed fast phase. It should be recognized that the steady-state approximation is not a precise treatment of this reaction scheme. A more detailed analysis, provided in the Appendix and summarized in Table 1, shows that the relationship between the observed kinetics and the microscopic rate constants is in fact much more complicated.

Kinetic partitioning distorts the rate constants and populations in any MHC-peptide system for which the rate of conversion of {MP*}1 to {MP*}2 is comparable to the other reaction rates. This means that, except in the case of the parallel system, the magnitude of the observed fast phase in a dissociation reaction underestimates the initial population of the complex {MP*}1, sometimes to a great degree. Consider the system shown in Scheme 4, which is also consistent with the dissociation curve in Fig. 2.
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(100%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>0.19 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>0.80 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(0%)</UL></LIM></C></R><R><C><SUB> 0.25 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>∼0</SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 4</UP>
In this case, all of the complex is initially present in the faster-dissociating form {MP*}1, yet the magnitude of the observed fast dissociation phase is only 18%.

Complete mechanisms

In the examples discussed so far, one or more microscopic reaction rates are negligible. However, most of the solution space for our dissociation curve consists of reaction systems in which all of the microscopic reactions proceed at significant rates. One example of such a system is shown in Scheme 5; like the other examples, it too generates the biphasic dissociation curve in Fig. 2.
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(20%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>0.015 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>0.060 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(80%)</UL></LIM></C></R><R><C><SUB> 1.14 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>0.026 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 5</UP>
This particular example is interesting because the two complexes are initially at equilibrium: the ratio of the initial concentrations of {MP*}2 to {MP*}1 is equal to K1right-left-harpoons 2 = k12/k21 = 4.

Although there is a diverse range of systems that are consistent with our observed dissociation curve, we can place limits on the values of the microscopic parameters of the actual MHC-peptide system. As shown in Fig. 4, the microscopic rate constant koff,1 must be less than or equal to the observed kfasto (1.2 h-1), and koff,2 must be less than or equal to kslowo (0.04 h-1). The rate constants for interconversion, k12 and k21, can range from 0 (in the case of the parallel scheme) to a value less than the observed kfasto. The fraction of complex initially in the form {MP*}1 must be greater than Fo, the observed magnitude of the fast dissociation phase, and can, in principle, be as high as 100%.

For simplicity, the observed magnitudes (Fo and So) and rate constants (kfasto and kslowo) of biphasic MHC-peptide dissociation curves in previous reports have often been taken as readouts of the initial concentrations and dissociation rate constants of the faster and slower dissociating complexes (Beeson and McConnell, 1994; Beeson et al., 1996; Rabinowitz et al., 1997). However, as the above examples demonstrate, this simple interpretation is valid only if the complexes interconvert very slowly relative to the dissociation reactions. Slow interconversion need not be the case, however. Experiments involving "regeneration" of a fast dissociation phase from partially dissociated MHC-peptide complexes show that the interconversion rate of isomeric complexes can be comparable to or faster than that of dissociation reactions (Schmitt et al., 1998b). In the absence of additional information, therefore, it is generally prudent to consider the full range of consistent systems when interpreting MHC-peptide dissociation kinetics.

Information from variation of binding time

Although it is generally not possible to determine the initial concentrations of the two isomeric complexes in a typical dissociation experiment, one can nevertheless manipulate these concentrations to some extent by varying the length of the binding incubation that precedes the dissociation. Because the two MHC-peptide complexes are almost certainly formed at different rates, the relative populations of the complexes should depend on the length of the binding incubation. Therefore, by changing the binding time, one can systematically vary the ratio of the complexes' concentrations for different dissociation experiments. This variation in concentrations shows up indirectly in the observed magnitude of fast dissociation phase, Fo, over a range of binding incubation times (Fig. 5).



View larger version (12K):
[in this window]
[in a new window]
 
FIGURE 5   The magnitude of observed fast phase, Fo, after binding incubations of 0.2, 0.5, 1, 2, 4, 8, 13, and 20 h for two simulations of a two-complex reaction system with the same observed dissociation rate constants as in Fig. 2. The simulations used the rate constants given in Experimental Techniques, with an initial peptide concentration of 100 µM and 1 µM MHC protein initially in the form {MPe} (open circle ) or MI (). (A) The value of Fo is greatest after short binding incubations and decreases for longer binding incubations. (B) The ratio of the observed slow and fast dissociation phases as a function of binding time. Note that the precise shapes of these curves depend upon reactions that precede complex formation, but the values extrapolated to binding times of 0 and infinity (1.074 and 32.4) are in the same ratio as the observed rate constants (kfasto/kslowo = 30) and are independent of the binding process. (C) A close-up of the curve from B, showing the ratio of the slow to the fast dissociation phase after short incubation times for the two simulations. The solid line illustrates the limiting value of the slope of this curve, as calculated from the kinetic model. The So/Fo ratio extrapolated to zero binding time is 1.074, corresponding to a magnitude of the fast phase F0o = 0.482.

Measurement of the magnitude of fast dissociation phase after different binding times provides additional information about the microscopic kinetics of the system. Fig. 5 A shows a plot of the observed magnitude of fast dissociation phase, Fo, for simulated dissociations after binding incubations of up to 20 h; this curve is an approximately exponential decay from a maximum of F0o = 0.48 at zero binding time, with a limiting value of 0.03 at long binding times. The data can also be plotted as the ratio of the observed slow and fast dissociation phases, producing the sigmoidal curve shown in Fig. 5 B. The precise shapes of these curves depend upon the rate of complex formation, which is limited by the rate at which MHC becomes available for binding the labeled peptide through dissociation of endogenous peptide or by activation of "dormant" empty MHC (Rabinowitz et al., 1998). However, the limits of the curves in zero and infinite binding time do not depend on these ancillary processes.

In the limit of zero binding time, the ratio of the isomeric complexes formed is entirely determined by the ratio of their binding rate constants kon,1 and kon,2. If the binding reaction is carried out under the same conditions of temperature and pH as the dissociation reactions, these binding rate constants are related to the dissociation and interconversion rate constants of the complexes through a thermodynamic cycle:
<FR><NU>k<SUB><UP>on,1</UP></SUB></NU><DE>k<SUB><UP>off,1</UP></SUB></DE></FR> <FR><NU>k<SUB>12</SUB></NU><DE>k<SUB>21</SUB></DE></FR> <FR><NU>k<SUB><UP>off,2</UP></SUB></NU><DE>k<SUB><UP>on,2</UP></SUB></DE></FR>=K<SUB><UP>M⇌MP1</UP></SUB>K<SUB><UP>MP1⇌MP2</UP></SUB>K<SUB><UP>MP2⇌M</UP></SUB>=1 (3)
Consequently, the molar fraction of complex in the form {MP*}1 extrapolated to zero binding time can be expressed as
X<SUB>{<UP>MP*</UP>}<SUB>1</SUB></SUB>(0)=<FR><NU>k<SUB><UP>off,1</UP></SUB>k<SUB>21</SUB></NU><DE>k<SUB><UP>off,1</UP></SUB>k<SUB>21</SUB>+k<SUB>12</SUB>k<SUB><UP>off,2</UP></SUB></DE></FR> (4)
Because the magnitude of the observed fast dissociation phase is related to the fraction of {MP*}1 through Eq. 18 (Table 1), measurement of the magnitude of fast dissociation phase in the limit of zero binding time, F0o, provides additional information about the values of the microscopic rate constants of the MHC-peptide system.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 1   Relationships between the observed and microscopic kinetic parameters of a two-complex MHC-peptide system

We can calculate the F0o values that would be observed for the candidate systems described by Schemes 1-5, using the rate constants provided. Of the five candidates discussed, Schemes 2, 4, and 5 are not consistent with the observed value F0o = 0.48 for the binding variation data in Fig. 5. In terms of the "kinetic space" of MHC-peptide systems shown in Fig. 4, our measurement of F0o restricts the range of possible candidates to the one-dimensional slice indicated by the dotted line.

Solving the equations in Table 1 for the parallel and sequential mechanisms gives general expressions for the limits of the microscopic parameters that are consistent with MHC-peptide dissociation experiments; these are shown in Table 2. Systems involving interconversion as well as direct dissociation of both complexes are characterized by values falling between the two extremes.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 2   Limits on the microscopic parameters consistent with observed dissociation curves

Like the magnitude of the fast dissociation phase at zero binding time, the value of Fo after infinite binding time is determined entirely by the system's microscopic interconversion and dissociation rate constants. In the limit of infinite binding time, the two isomeric MHC-peptide complexes are in equilibrium with each other; hence their populations are in the ratio of k12/k21, and the fraction of complex in the form {MP*}1 is
X<SUB>{<UP>MP*</UP>}<SUB>1</SUB></SUB>(<UP>eq</UP>)=<FR><NU>k<SUB>12</SUB></NU><DE>k<SUB>12</SUB>+k<SUB>21</SUB></DE></FR> (5)
The magnitude of the fast dissociation phase after binding to equilibrium is related to this fraction of {MP*}1 and is therefore also characteristic of the MHC-peptide system.

Interestingly, the magnitudes of the fast and slow exponential phases at zero and infinite binding time are related to one another in a fairly simple way. It can be shown that the kinetic equations in Table 1 lead to the relationship
<FENCE><FR><NU>S<SUP><UP>o</UP></SUP></NU><DE>F<SUP><UP>o</UP></SUP></DE></FR></FENCE><SUB><UP>eq</UP></SUB>=<FR><NU>k<SUP><UP>o</UP></SUP><SUB><UP>fast</UP></SUB></NU><DE>k<SUP><UP>o</UP></SUP><SUB><UP>slow</UP></SUB></DE></FR> <FENCE><FR><NU>S<SUP><UP>o</UP></SUP></NU><DE>F<SUP><UP>o</UP></SUP></DE></FR></FENCE><SUB>0</SUB> (6)
where So is the normalized magnitude of the slow phase. Because of this, it is useful to determine the value of Fo at both zero and infinite time as a check on the two-complex model. In a sigmoidal plot of slow/fast phase versus binding time, the extrapolated values at zero and infinity should be in the same ratio as the slow and fast observed dissociation rate constants.

Another useful feature of the So/Fo curve is illustrated in Fig. 5 C. At short incubation times, it can be shown that the curve approaches the line y = y0 + (1/2)y0(kfasto - kslowo)t. This allows one to estimate the value of the So/Fo ratio---and therefore the value of Fo---at zero time from the value at a binding time t close to zero, using the formula
<FENCE><FR><NU>S<SUP><UP>o</UP></SUP></NU><DE>F<SUP><UP>o</UP></SUP></DE></FR></FENCE><SUB>0</SUB>≅<FENCE><FR><NU>2</NU><DE>2+(k<SUP><UP>o</UP></SUP><SUB><UP>fast</UP></SUB>−k<SUP><UP>o</UP></SUP><SUB><UP>slow</UP></SUB>)t</DE></FR></FENCE><FENCE><FR><NU>S<SUP><UP>o</UP></SUP></NU><DE>F<SUP><UP>o</UP></SUP></DE></FR></FENCE><SUB><UP>t</UP></SUB> (7)
It is important to note that the foregoing discussion is valid only if the binding conditions used to prepare the MHC-peptide complex are identical to the conditions of the dissociation reactions (except for peptide concentration). If binding and dissociation are performed under different conditions, the magnitude of the fast dissociation phase will depend on the microscopic rate constants under both sets of conditions, which greatly complicates the analysis of dissociation behavior for various binding times.

Even after dissociation is measured following a range of binding times, there is still a continuum of candidate systems that are consistent with the kinetic data, indicated by the dotted line in Fig. 4. To determine where our system lies along this line, we need measurements of the concentrations of the two complexes; the indirect readout provided by the magnitude of observed dissociation fast phase is not sufficient. Conventional kinetic measurements of MHC-peptide dissociation do not distinguish between different isomeric forms of the MHC-peptide complex, but other techniques may be used to detect and measure the two complexes. For example, 19F NMR has been used to directly measure the relative populations of different complexes in several MHC-peptide systems (Schmitt et al., 1998a).

Consider the example of the peptide PCC, consisting of residues 89-104 of pigeon cytochrome c, bound to the class II MHC molecule I-Ek. Kinetic studies (Schmitt et al., 1998b) have shown that this system shows distinctly biphasic dissociation kinetics at pH 5.3 and 25°C, with observed rate constants of kfasto = 1.70 -0.29/+0.44 h-1 and kslowo = 0.01397 -0.00006/+0.00006 h-1. From the reported dissociation results following different binding times, we can estimate the ratio of the slow to the fast phase at zero binding time as (Fo/So)0 = 2.72 -0.40/+0.47, corresponding to a fast phase at zero binding time of F0o = 0.27 -0.03/+0.03. Separate NMR measurements (Schmitt et al., 1998a) have shown that when the complexes are isolated at pH 7.0 and 25°C, the fraction of complex in the form {MP*}1 is X{MP*}1 = 0.5 -0.045/+0.045, and the magnitude of the fast phase in a subsequent dissociation at pH 5.3 and 25°C is Fo(X1 = 0.5) = 0.141 -0.013/+0.013. Using these values, Eq. 4, and the equations in Table 1, we calculate that the microscopic rate constants of this system are koff,1 = 0.50 -0.15/+0.25 h-1; k12 = 1.17 -0.28/+0.40 h-1; k21 = 0.043 -0.027/+0.068 h-1; koff,2 = 0.0013 -0.0013/+0.0072 h-1. Based on these calculated rate constants, the equilibrium fraction of {MP*}1 at pH 5.3 is X{MP*}1(eq) = 0.036 -0.026/+0.078; this is very different from the value (~0.5) measured at pH 7.0. This example shows that some of the microscopic parameters, namely k21 and koff,2, are difficult to determine precisely; nevertheless the calculated ranges do provide useful information about the system.

Hidden biphasic kinetics

Observation of a biphasic dissociation curve proves that isomeric forms of the MHC-peptide complex are present, even though they may not be detected directly. The inverse of this statement is not true, however. That is, the absence of an observed fast dissociation phase does not rule out the presence of multiple complexes. Indeed, if the dynamics of the MHC-peptide complexes are dominated by interconversion rather than dissociation, it can be virtually impossible to observe a biphasic dissociation curve, even though two different complexes are initially present.

For complexes prepared under the same conditions as the dissociation measurement, the maximum magnitude of observable fast phase is F0o, as can be seen in Fig. 5. From Eq. 18, we can derive an upper limit for F0o in terms of the interconversion rate of the complexes. For a relative interconversion rate of
k<SUB><UP>i</UP></SUB>=<FR><NU>(k<SUB>12</SUB>+k<SUB>21</SUB>)</NU><DE>(k<SUB><UP>off,1</UP></SUB>+k<SUB><UP>off,2</UP></SUB>)</DE></FR> (8)
the value of F0o cannot be greater than 1/(1 + ki), and generally will be less than this, especially if the complex {MP*}2 has a significant dissociation rate. For example, if two MHC-peptide complexes have dissociation rate constants of 0.09 h-1 and 0.01 h-1 and interconversion rate constants of 6 h-1 and 4 h-1 (Scheme 5), then ki = 100, and so F0o must be less than 0.01. 
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(86%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>4 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>6 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(14%)</UL></LIM></C></R><R><C><SUB> 0.09 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>0.01 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 6</UP>
Simulation of this system shows that the observed dissociation rate constants are 10.1 h-1 and 0.042 h-1, but the value of F0o is only 0.0037 (Fig. 6, circles). This fast phase, with magnitude less than 0.4%, is impossible to detect using current methods.



View larger version (20K):
[in this window]
[in a new window]
 
FIGURE 6   "Hidden" biphasic kinetics in MHC-peptide dissociation reactions: dissociation curves for the systems described by Schemes 6 (open circle ) and 7 (), after a very short binding incubation to maximize the amount of fast phase observed. The solid lines are single-exponential curve fits to the data.

Can systems with such rapid interconversion rates be physiologically relevant? In the example described above, the two complexes have characteristic lifetimes with respect to dissociation (tau off = 1/koff) of 11 h and 100 h. By comparison, their lifetimes with respect to interconversion are 10 min and 15 min. Although the interconversion reactions are much more rapid than dissociation, interconversion is still slow enough to keep the two complexes distinct from a physiological point of view, even though they are not kinetically resolvable. In a physiological context, the important time scale is the lifetime of the ternary MHC-peptide/TCR complex, which is on the order of seconds (Matsui et al., 1991, 1994).

Significantly, the magnitude of a fast dissociation phase may be undetectably small, even though the system initially contains a large fraction of {MP*}1. In the example shown in Scheme 6, the fraction of complex {MP*}1 present after a very short binding incubation is 85.7%; nevertheless, only a 0.4% fast phase is observed. The relatively fast interconversion renders the fast-dissociating complex "kinetically invisible" in dissociation experiments.

Another way in which a fast-dissociating MHC-peptide complex can evade kinetic detection is to have similar rate constants for the observed fast and slow dissociation phases. If the rate of interconversion is comparable to the rate of dissociation and the equilibrium constant between the two complexes favors the faster-dissociating {MP*}1, then the two observed rate constants can be very similar. Consider the system shown in Scheme 7:
<AR><R><C><LIM><OP>{<UP>MP*</UP>}<SUB>1</SUB></OP><UL>(99.6%)</UL></LIM></C><C><LIM><OP><ARROW>⇌</ARROW></OP><LL><SUB>0.1 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></LL><UL><SUB>0.004 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></UL></LIM></C><C><LIM><OP>{<UP>MP*</UP>}<SUB>2</SUB></OP><UL>(0.4%)</UL></LIM></C></R><R><C><SUB> 0.09 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB>↘</C><C></C><C>↙<SUB>0.01 <UP>h</UP><SUP><UP>−1</UP></SUP></SUB></C></R><R><C></C><C><UP>M</UP>+<UP>P*</UP></C></R></AR>

<UP>Scheme 7</UP>
Here the microscopic dissociation rate constants are the same as in Scheme 6, but the relative rate of interconversion, ki, is only (0.004 h-1 + 0.1 h-1)/(0.09 h-1 + 0.01 h-1) = 1.04. Even though the maximum magnitude of the fast dissociation phase for this system is a relatively high 30.4%, the observed rate constants of this system, 0.124 h-1 and 0.0805 h-1, differ by only a factor of 1.5 (Fig. 6, squares). A double-exponential dissociation curve with such similar rate constants is virtually impossible to distinguish from a single-exponential curve, given the noise in current experiments (on the order of 1%). In general, two-complex systems like this one, with small equilibrium constants of interconversion (Ki = k12/k21), will have similar observed dissociation rate constants.

The examples of Schemes 6 and 7 illustrate that biphasic dissociation kinetics of a two-complex system can be difficult to resolve if either ki or 1/Ki is greater than 1. Many two-complex systems may fall into one or both of the "unresolvable" categories illustrated by these schemes. This suggests the possibility that isomeric MHC-peptide complexes may be much more common than previously believed.

Although biphasic kinetics of an MHC-peptide system may be unresolvable under typical dissociation conditions, other conditions may produce more clearly biphasic behavior. MHC-peptide dissociation reactions are intrinsically relaxation processes, in which a sample prepared in the presence of labeled peptide is perturbed by removing unbound peptide and/or adding an excess of unlabeled competitor peptide. But other conditions, such as temperature and pH, can also be changed when dissociation is initiated. For example, a bound MHC-peptide complex may be stored at a low temperature (4°C) and then warmed to carry out a dissociation experiment at 25°C or 37°C. Alternatively, samples may be prepared at pH 5.3 for dissociations carried out at pH 7.0.

As discussed earlier, if the MHC-peptide complex is prepared under the same conditions as the dissociation reaction, the maximum magnitude of the dissociation fast phase is observed in the limit of zero binding time; this is also shown in Fig. 7 A. However, if the conditions are not the same, it may be possible to prepare the system such that the initial fraction of fast-dissociating {MP*}1---and hence the observed fast dissociation phase---is enhanced. This enhancement may permit the resolution of biphasic behavior for some systems. Fig. 7 A also shows that an interesting situation can result when the complex is prepared so that the initial fraction of {MP*}1 is very small. If the amount of {MP*}1 is small enough, the magnitude of observed fast dissociation can be zero or even negative! Such a dissociation appears to have a "lag period" at early time points (Fig. 7 B).



View larger version (27K):
[in this window]
[in a new window]
 
FIGURE 7   The magnitude of the observed dissociation fast phase generally does not reflect the initial amount of fast-dissociating complex. (A) The relationship between the magnitude of the observed fast phase, Fo, and the initial fraction of the fast-dissociating complex {MP*}1 is linear and depends on the interconversion and dissociation rates of the two complexes. For a given set of observed dissociations, consistent mechanisms span the range from parallel (e.g., Scheme 1), indicated by the heavy line labeled par, to sequential (e.g., Scheme 3), indicated by line seq. A dashed line (--- - ---) corresponds to the actual mechanism used for the simulation of Figs. 2 and 4. The observed fast dissociation phases in the limits of zero and infinite binding time are indicated by the vertical dotted lines. (B) A dissociation curve with a negative fast phase. Shown is a simulated dissociation curve for the two-complex system described in Scheme 3, starting from an initial population of 0.8% {MP*}1 and 99.2% {MP*}2; for this curve, Fo = -0.034. Inset: The early portion of the dissociation curve deviates from a single-exponential curve (dotted line).

Another useful perturbation of the MHC-peptide system involves the addition of other chemical species to the reaction system. Dimethyl sulfoxide has been demonstrated to enhance peptide release from MHC (Schmitt et al., 1998a), as has the peptide dynorphin A (1-13) (de Kroon and McConnell, 1993; Schmitt, 1999). The chaperone molecule DM also catalyzes peptide release from MHC (Denzin and Cresswell, 1995). These substances appear to act disproportionately on the fast-dissociating isomer of the MHC-peptide complex, and as such may provide an aid in resolving two-complex dissociation kinetics.

    CONCLUSIONS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

Technical obstacles to measuring MHC-peptide reaction kinetics have largely been overcome (Witt and McConnell, 1993), but it is clear that such studies by themselves provide limited information about the isomers' dynamics. Although existing kinetic measurements cannot provide complete information about the reactions between MHC and peptides, the application of different techniques for characterizing MHC-peptide complexes and kinetics can help to focus the picture of these interactions.

This analysis also offers the prospect of resolving the mechanistic effects observed in changes in temperature (Witt and McConnell, 1994), pH (Witt and McConnell, 1991; Boniface et al., 1993; Schmitt et al., 1998b), and peptide sequence (Dornmair et al., 1991; Beeson et al., 1996; Liang et al., 1996; Schmitt et al., 1998b), which may be relevant to the biological function of these multiple-isomer MHC-peptide complexes. A further objective of this work is to extend our analysis to the binding reactions of peptides to MHC molecules. These reactions, which have obvious importance in the understanding of MHC-peptide interactions, have been demonstrated to show complicated kinetics (Tampé and McConnell, 1991; Witt and McConnell, 1991, 1992; de Kroon and McConnell, 1993; Mason and McConnell, 1994; Rabinowitz et al., 1998; Natarajan et al., 1999) and are not merely the "reverse" of the dissociation reactions.

Finally, we hope to use this work to guide single-molecule fluorescence studies of MHC-peptide systems. Structural differences between isomeric complexes may influence environment-sensitive fluorescent labels. Any difference in fluorescence between complexes could be exploited in single-molecule studies, for which the amount of fluorescence would provide a direct measure of the state of the bound peptide. Intensity correlation functions may also permit the determination of isomerization rates (Wennman et al., 1997). Single molecule techniques are discussed in a recent issue of Science (Gimzewski and Joachim, 1999; Mehta et al., 1999; Moerner and Orrit, 1999; Weiss, 1999).

    APPENDIX: ANALYTICAL SOLUTION OF THE RATE EQUATIONS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

In general, for a system of N chemical species interacting by first-order (A right-arrow B) reactions, the time-dependent concentrations of the species can be described by an (N - 1)-exponential curve plus a constant equilibrium value (Moore and Pearson, 1981). For the irreversible dissociation of isomeric MHC-peptide complexes shown in Fig. 1 B, N = 3 (complexes {MP*}1 and {MP*}2, and the dissociated state), and the equilibrium concentrations of the complexes are zero. Hence the total concentration of the two complexes during an irreversible dissociation reaction can be described by a biexponential curve.

Expressions for the macroscopic rate constants kfasto and kslowo of this biexponential curve can be found by integrating the rate equations for the reactions shown in Fig. 1 B. The concentrations of the three MHC species are described by the differential equations
<FR><NU><UP>d</UP>[{<UP>MP*</UP>}<SUB>1</SUB>]</NU><DE><UP>d</UP>t</DE></FR>=<UP>−</UP>(k<SUB><UP>off,1</UP></SUB>+k<SUB>12</SUB>)[{<UP>MP*</UP>}<SUB>1</SUB>]+k<SUB>21</SUB>[{<UP>MP*</UP>}<SUB>2</SUB>] (9)

<FR><NU><UP>d</UP>[{<UP>MP*</UP>}<SUB>2</SUB>]</NU><DE><UP>d</UP>t</DE></FR>=k<SUB>12</SUB>[{<UP>MP*</UP>}<SUB>1</SUB>]−(k<SUB>21</SUB>+k<SUB><UP>off,2</UP></SUB>)[{<UP>MP*</UP>}<SUB>2</SUB>] (10)

<FR><NU><UP>d</UP>[<UP>M</UP>]</NU><DE><UP>d</UP>t</DE></FR>=k<SUB><UP>off,1</UP></SUB>[{<UP>MP*</UP>}<SUB>1</SUB>]+k<SUB><UP>off,2</UP></SUB>[{<UP>MP*</UP>}<SUB>2</SUB>] (11)
Solving this system of differential equations gives the expressions for the macroscopic rate constants kfasto and kslowo shown in Table 1, as well as expressions for the coefficients Fo and So.

The observed magnitudes of fast and slow phase, Fo and So, correspond to the sums of the two complexes' fast and slow coefficients. Expressions for these magnitudes can also be derived by examining the total dissociation rate of the complexes in terms of the concentrations and rate constants involved. From the mechanism in Fig. 1 B, we can express the initial rate of the dissociation reaction of an MHC-peptide system as
<UP>−</UP><FR><NU><UP>d</UP>[{<UP>MP*</UP>}]</NU><DE><UP>d</UP>t</DE></FR>(0)=[{<UP>MP*</UP>}<SUB>1</SUB>]k<SUB><UP>off,1</UP></SUB>+[{<UP>MP*</UP>}<SUB>2</SUB>]k<SUB><UP>off,2</UP></SUB> (12)
Normalized with respect to the total initial concentration of the MHC-peptide complex, this becomes
<UP>−</UP><FR><NU><UP>dMP*</UP></NU><DE><UP>d</UP>t</DE></FR>(0)=X<SUB>{<UP>MP*</UP>}<SUB>1</SUB></SUB>k<SUB><UP>off,1</UP></SUB>+X<SUB>{<UP>MP*</UP>}<SUB>2</SUB></SUB>k<SUB><UP>off,2</UP></SUB> (13)
where X{MP*}1 and X{MP*}2 are the initial fractions of the two complexes.

A similar rate equation can be written based on the empirical biphasic concentration curve (Eq. 2) that is fit to the dissociation data. Differentiating this (already normalized) function with respect to time gives the initial dissociation rate as
<UP>−</UP><FR><NU><UP>dMP*</UP></NU><DE><UP>d</UP>t</DE></FR>(0)=F<SUP><UP>o</UP></SUP>k<SUP><UP>o</UP></SUP><SUB><UP>fast</UP></SUB>+S<SUP><UP>o</UP></SUP>k<SUP><UP>o</UP></SUP><SUB><UP>slow</UP></SUB> (14)
Combining Eqs. 13 and 14 gives the relation
X<SUB>{<UP>MP*</UP>}<SUB>1</SUB></SUB>k<SUB><UP>off,1</UP></SUB>+X<SUB>{<UP>MP*</UP>}<SUB>2</SUB></SUB>k<SUB><UP>off,2</UP></SUB>=F<SUP><UP>o</UP></SUP>k<SUP><UP>o</UP></SUP><SUB><UP>fast</UP></SUB>+S<SUP><UP>o</UP></SUP>k<SUP><UP>o</UP></SUP><SUB><UP>slow</UP></SUB> (15)
Given that X{MP*}1 + X{MP*}2 = 1 and Fo + So = 1, Eq. 15 can be rearranged to give the expressions for Fo and So shown in Table 1.

    ACKNOWLEDGMENTS

The authors thank Lutz Schmitt for many helpful discussions.

This material is based upon work supported under a National Science Foundation Graduate Fellowship. This work was also supported by grant 5R37 AI13587-23 from the National Institutes of Health.

    FOOTNOTES

Received for publication 16 December 1998 and in final form 2 August 1999.

Address reprint requests to Dr. Harden M. McConnell, Department of Chemistry, Stanford University, Stanford, CA 94305-5080. Tel.: 650-723-4571; Fax: 650-723-4943; E-mail: harden{at}leland.stanford.edu.

    REFERENCES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL TECHNIQUES
RESULTS AND DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES