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
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 |
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 |
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:
|
(1)
|
|
(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 |
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*
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.
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.
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.
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.
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.
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
K1
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} ( ) 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:
|
(3)
|
Consequently, the molar fraction of complex in the form
{MP*}1 extrapolated to zero binding time can be
expressed as
|
(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.
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
|
(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
|
(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
|
(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
|
(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.
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 ( ) 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
(
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:
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 |
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
).
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
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.
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.