Peptide Modulation of Class I Major Histocompatibility Complex Protein Molecular Flexibility and the Implications for Immune Recognition*

Background: Peptide modulation of MHC flexibility can affect recognition by immune receptors. Results: Different peptides alter the flexibility of the MHC protein at sites around the peptide-binding groove. Conclusion: Peptide modulation of MHC flexibility is not limited to specific peptides or isolated regions. Significance: Peptide modulation of MHC flexibility indicates an extension of antigenicity from the peptide to the MHC. T cells use the αβ T cell receptor (TCR) to recognize antigenic peptides presented by class I major histocompatibility complex proteins (pMHCs) on the surfaces of antigen-presenting cells. Flexibility in both TCRs and peptides plays an important role in antigen recognition and discrimination. Less clear is the role of flexibility in the MHC protein; although recent observations have indicated that mobility in the MHC can impact TCR recognition in a peptide-dependent fashion, the extent of this behavior is unknown. Here, using hydrogen/deuterium exchange, fluorescence anisotropy, and structural analyses, we show that the flexibility of the peptide binding groove of the class I MHC protein HLA-A*0201 varies significantly with different peptides. The variations extend throughout the binding groove, impacting regions contacted by TCRs as well as other activating and inhibitory receptors of the immune system. Our results are consistent with statistical mechanical models of protein structure and dynamics, in which the binding of different peptides alters the populations and exchange kinetics of substates in the MHC conformational ensemble. Altered MHC flexibility will influence receptor engagement, impacting conformational adaptations, entropic penalties associated with receptor recognition, and the populations of binding-competent states. Our results highlight a previously unrecognized aspect of the “altered self” mechanism of immune recognition and have implications for specificity, cross-reactivity, and antigenicity in cellular immunity.


T cells use the ␣␤ T cell receptor (TCR) to recognize antigenic peptides presented by class I major histocompatibility complex proteins (pMHCs) on the surfaces of antigen-presenting cells. Flexibility in both TCRs and peptides plays an important role in antigen recognition and discrimination. Less clear is the role of flexibility in the MHC protein; although recent observations have indicated that mobility in the MHC can impact TCR recognition in a peptide-dependent fashion, the extent of this behavior is unknown.
Here, using hydrogen/deuterium exchange, fluorescence anisotropy, and structural analyses, we show that the flexibility of the peptide binding groove of the class I MHC protein HLA-A*0201 varies significantly with different peptides. The variations extend throughout the binding groove, impacting regions contacted by TCRs as well as other activating and inhibitory receptors of the immune system. Our results are consistent with statistical mechanical models of protein structure and dynamics, in which the binding of different peptides alters the populations and exchange kinetics of substates in the MHC conformational ensemble. Altered MHC flexibility will influence receptor engagement, impacting conformational adaptations, entropic penalties associated with receptor recognition, and the populations of binding-competent states. Our results highlight a previously unrecognized aspect of the "altered self" mechanism of immune recognition and have implications for specificity, cross-reactivity, and antigenicity in cellular immunity.
␣␤ T cell receptors (TCRs) 4 on the surfaces of CD8ϩ T cells recognize antigenic peptides bound and presented by class I major histocompatibility complex proteins (class I pMHC complexes) to initiate and propagate an antigen-dependent immune response. Both TCR and peptide conformational properties play key roles in antigen recognition and discrimination. Conformational changes in the TCR complementaritydetermining region binding loops often accompany binding (1), and measurements of complementarity-determining region loop dynamics have directly linked loop flexibility to TCR cross-reactivity (2,3). Likewise, conformational changes in peptides often accompany TCR binding (e.g. Refs. 4 and 5), and peptide flexibility can be altered through peptide modifications or MHC polymorphisms (6 -8).
Less appreciated, however, is the role that MHC flexibility plays in TCR recognition. Examples of conformational changes in MHC peptide-binding domains occurring upon TCR binding are mostly limited to small shifts in the shapes and positions of the ␣1 or ␣2 helices (e.g. Ref. 9). In one telling case, however, the conformation of the ␣2 helix of the human class I MHC protein HLA-A*0201 (HLA-A2) was shown to undergo a large reorganization upon binding of the A6 TCR (10). Notably, this reorganization is dependent on the peptide as it has not been observed upon recognition of other HLA-A2-presented peptides by the same TCR (11,12).
In protein binding, crystallographically observed conformational differences have been shown to occur in regions with enhanced motions as the lower energy barriers that facilitate conformational changes translate into faster rates of conformational exchange (13). Indeed, with HLA-A2, the mobility of the region that undergoes a structural change upon binding of the A6 TCR was shown to vary with different peptides (10). This may not be a unique observation as similar effects are believed to contribute to differential T cell recognition of native and modified MART-1 peptides (6), and MHC flexibility can be altered by micropolymorphisms (7,14,15). As demonstrated by the A6 TCR, peptide-dependent MHC dynamics can impact TCR recognition by altering the barriers for conformational adjustments, influencing the entropic costs for receptor binding and shifting the populations of binding-competent states. Although altered MHC flexibility has not been apparent structurally, differences in flexibility are generally difficult to ascertain from examination of structure alone.
Here, we explore the extent to which different peptides modify the flexibility of the peptide-binding domain of the class I MHC protein in detail. Using hydrogen/deuterium exchange, fluorescence anisotropy, and an analysis of the class I pMHC structural database, we found a surprising degree of peptide dependence to the motional properties of the HLA-A2 peptidebinding groove. The results were not limited to one region, but extended throughout both ␣-helices as well as the ␤-sheet floor of the groove. Regions affected included those that interact with TCRs as well as other activating and inhibitory receptors of the immune system, such as the CD8 coreceptor and natural killer receptors. Our observations emphasize that association of peptides with HLA-A2 proteins physically alters the properties of the protein. Indeed, our results encompass yet broaden the "altered self" model of immune recognition, which emerged from the finding that receptors of the cellular immune system recognize MHC proteins altered through the binding and presentation of antigens (16). Our findings have implications for immunological specificity, cross-reactivity, and ultimately, the determinants of antigenicity for class I MHC-presented peptides.

EXPERIMENTAL PROCEDURES
Proteins and Peptides-pMHC complexes were produced by refolding bacterially expressed inclusion bodies of the HLA-A2 heavy chain and the ␤ 2 -microglobulin subunit in the presence of excess peptide as described previously (17). Refolded complexes were purified by a combination of ion exchange and size exclusion chromatography. Peptides were synthesized commercially (GenScript). Cysteine mutations in HLA-A2 were performed by PCR mutagenesis and confirmed by sequencing.
Hydrogen/Deuterium Exchange-Mass Spectrometry-Hydrogen/deuterium exchange-mass spectrometry (HDX-MS) was performed as described previously (18). Briefly, exchange was initiated by diluting protein samples at 25°C in 25 mM HEPES, 50 mM NaCl (pH 7.4) 10-fold with the same buffer made with 99.9% 2 H 2 O. pMHC concentrations were 30 M after dilution. Excess free peptide was present at a concentration of 30 M. Exchange was performed within 1 h after purification over a size exclusion column. After dilution into 2 H 2 O, 5-l aliquots were removed at various time points. For each aliquot, exchange was quenched at 0°C with 100 l of 0.1% trifluoroacetic acid (TFA) at pH 2.4. Pepsin (8 l at 1 mg/ml) was added and digestion performed for 5 min before freezing in liquid nitrogen.
The identities of pepsin-generated fragments of HLA-A2 were previously established in our laboratory by LC/MS/MS (18). The analysis of the exchange behavior of these fragments was performed by MALDI mass spectrometry as described previously (18). Briefly, a matrix solution of 2,5-dihydroxybenzoic acid in 50:50 0.1% TFA (pH 2.4) and acetonitrile was mixed 1:1 with the acidified protein digest on a chilled MALDI plate. The plate was dried under vacuum to minimize back exchange; experiments with fully deuterated fragments indicated that back exchange never exceeded 5%. Spectra were acquired on a Bruker Autoflex III Smartbeam MALDI-TOF mass spectrom-eter. Data were collected out to hours, although only data to 10 min were selected for comparison as shown in Fig. 1D. This timescale was chosen as it shows clear differences among the samples, and each sample reached an apparent plateau after 10 min, permitting an exponential fit and an estimate of error.
Fluorescence Anisotropy-Protein labeling and analysis by fluorescence anisotropy were performed as described previously (10). Briefly, purified single-cysteine mutants in 20 mM Na 2 HPO 4 , 75 mM NaCl (pH 7.4) were mixed with a 10-fold molar excess of Alexa Fluor 488 C5-maleimide or BODIPY-FL N-(2-aminoethyl) maleimide along with 10-fold excess tris(2carboxyethyl)phosphine. Parallel reactions with wild-type protein were performed to investigate nonspecific labeling. After mixing protein and label for 1 h at room temperature, excess label was removed by dialysis followed by chromatography. Labeling efficiencies measured spectrophotometrically varied between 92 and 97%. UV images of reduced and nonreduced SDS-PAGE gels verified that fluorescence emanated only from the chain with the cysteine mutant, and fluorescence measurements confirmed that wild-type protein had insignificant amounts of nonspecifically incorporated label.
Steady-state anisotropy measurements were performed in a Beacon 2000 instrument (Invitrogen) in 20 mM Na 2 HPO 4 , 75 mM NaCl (pH 7.4) buffer with excess free peptide (between 20and 100-fold molar excess). Anisotropy measurements on individual samples were repeated 50 times, and samples were measured in duplicate utilizing two independently prepared samples. The values in Fig. 2 and Tables 1 and 2 were collected at 100 nM protein concentration and 25°C. Additional data for most samples were also collected at concentrations of 10 nM, 500 nM, and 1 M and at 4 and 35°C.
Measurements on denatured samples were performed by diluting protein in 8 M urea (pH 10) (19). Refolding of fluorescently labeled heavy chain cysteine mutants with new peptides was performed by denaturing a 1-ml solution of fully assembled, labeled Tax⅐HLA-A2 as above, dialyzing against 500 ml of denaturant over 24 h with three buffer changes to remove peptide, and then using the labeled heavy chain in new small scale refolding reactions with different peptides.
Structural Analyses-Solvent-accessible surface areas were calculated from peptide⅐HLA-A2 crystallographic structures using Discovery Studio 3.5 (Accelrys), using a grid-based algorithm with 960 points/Å 2 and a probe radius of 1.4 Å. Analyses of helix geometry for peptide⅐HLA-A2 and peptide⅐H-2K b structures were performed with HELANALϩ (20). Two-dimensional plots of peptide-protein interactions were generated with LigPlotϩ (21). When multiple molecules were present in the asymmetric unit, calculations were performed on only the first molecule.

RESULTS
Peptides Alter the Hydrogen/Deuterium Exchange Behavior of the HLA-A2 Peptide-binding Domain-Although NMR is the most powerful approach for characterizing protein flexibility, the size and complexity of pMHC complexes has made their characterization by NMR challenging. We thus began examining whether and how different peptides alter HLA-A2 flexibility by monitoring hydrogen/deuterium exchange via proteolysis followed by mass spectrometry (HDX-MS). In HDX-MS, protein prepared in aqueous solution is incubated in 2 H 2 O and digested into fragments and deuteron incorporation for each fragment is determined via mass spectrometry (22). Increased levels of deuteration indicate that backbone amides are more accessible to exchange through greater sampling of more open, solvent-exposed states, whereas lower levels of incorporation indicate that backbone amides are protected from exchange through less frequent excursions into solvent-exposed states. HDX-MS has been utilized in a number of cases to explore protein motions, including a recent study of conformational lability in class II MHC proteins (23) as well as our recent study of the effects of binding on the properties pMHC complexes and TCRs (18). Changes in hydrogen/deuterium exchange upon events such as ligand binding result from alterations in exchange kinetics for backbone amides, which manifest as a different number of deuterons incorporated as a function of incubation time.
We performed HDX-MS on four different peptide⅐HLA-A2 complexes. The peptides chosen were 9-or 10-mers for which crystallographic structures of the pMHC complexes were available ( Fig. 1A) (10, 24 -26). Care was taken to ensure that the data were not influenced by peptide dissociation. Each peptide chosen was a well characterized peptide with optimal anchors for HLA-A2, and exchange was performed on freshly purified samples at micromolar concentrations well above the nanomolar K D values believed to characterize the interaction of optimal peptides with class I MHC proteins. (Note that accurate K D values are difficult to measure for class I MHC proteins as the peptide-free molecule is unstable, rapidly aggregates, and cannot be directly purified. However, the half-life of the Tax⅐HLA-A2 complex is 8 h at 25°C with an estimated K D of 18 nM (19), and the T m values of the Tax and Flu M1 complexes measured by thermal denaturation are identical within error (24,27). Further, as indicated below, identical results were obtained with and without the presence of excess peptide in the exchange reactions. Due to its instability, measurements on the peptide-free protein could not be performed.) Our analysis focused on fragments that encompass the HLA-A2 peptide-binding domain, including three fragments that span the ␣1 helix, two that span the ␣2 helix, and one that forms a strand of the C-terminal half of the ␤-sheet floor of the binding groove (Fig. 1B). These regions were chosen as they are the primary sites of interaction with TCRs, and in the case of the ␤-sheet region, they interact with other activating and inhibitory receptors of the immune system. These regions are not glycosylated in class I MHC proteins and are not likely to be impacted by the potential formation of TCR-pMHC clusters in an immunological synapse.
The exchange behavior of each fragment followed single or double exponential kinetics as typically seen for HDX-MS data (Fig. 1, C and D). Exchange in each fragment varied with peptide, which we quantified as the degree of deuteration after 10 min of exchange (Fig. 1E). The results were reproducible with repeated measurements on separate samples (supplemental Fig. S1a). The values reported in Fig. 1 included 30 M excess peptide, although experiments without additional peptide yielded identical results (supplemental Fig. S1b).
The variation in exchange behavior with different peptides is consistent with a peptide dependence on the motions of the HLA-A2-binding domain. However, the results could also be explained by structural variations in the different complexes. The latter seems unlikely as the peptide-binding domains of the four peptide⅐HLA-A2 complexes superimpose with an average root mean square deviation of 1.2 Å. However, an influence from the different structures and compositions of the peptides is possible. To explore both possibilities, we calculated the solvent-accessible surface areas of the seven fragments in the four pMHC structures (Fig. 1F). The surface areas of each fragment were very similar in all four structures, and there was no correlation between the degree of deuteration and solvent accessibility. Thus, the peptide-dependent variation in the HDX-MS data cannot easily be explained by differences in the crystallographic structures or in the compositions of the peptides.
Fluorescence Anisotropy Reveals Substantial Local Variation in HLA-A2 Flexibility with Different Peptides-To more directly assess the peptide dependence of HLA-A2 flexibility, we examined protein motions using fluorescence anisotropy. When a fluorophore is excited with polarized light, emission is depolarized by the molecular motion that occurs over the fluorescence lifetime. In a protein covalently labeled with a fluorescent probe, motional modes leading to depolarization are protein tumbling, rotation of the probe around its tether, and the influence of backbone and neighboring side chain motion on probe dynamics. As different peptides do not alter the size and shape of the pMHC complex, protein tumbling will be invariant with peptide. Thus, provided a probe attached to the MHC does not interact directly with a peptide, differences in fluorescence depolarization can be attributed to the influence that peptides have on MHC solution structure and dynamics.
We studied the peptide-dependent variation of MHC flexibility by introducing cysteine mutations at six positions within the ␣-helices of the HLA-A2 heavy chain groove and two positions within the ␤-sheet that forms the floor of the groove. The six positions on the ␣-helices encompass regions contacted by TCRs and were within the fragments studied by HDX-MS. The positions on the ␤-sheet were also in a fragment identified by HDX-MS and are in proximity to the binding sites for the ␤ 2 -microglobulin subunit, the CD8 coreceptor, and, for murine class I MHC proteins, the Ly49 class of inhibitory receptors (28,29). The positions chosen for labeling have side chains in the wild-type protein that are polar or charged and, in the structures of the various pMHC complexes, point away from the protein and make no contacts to peptide ( Fig. 2A and  supplemental Fig. S2).
Each mutant pMHC complex was independently labeled with Alexa Fluor 488 C5-maleimide. We examined the four peptides studied by HDX-MS as well as gp100 T2M , another tight binding peptide whose structure bound to HLA-A2 is available (30). Single-cysteine mutants were refolded with the five peptides, fluorescently labeled, and repurified prior to collecting data. Excess peptide was included to ensure stability of the complex, and each set of measurements was performed on two independently prepared samples. To help benchmark the measurements, we also collected data for chemically denatured samples as well as unconjugated fluorescein.
We observed a dramatic variation in anisotropy with the position of the probe (Fig. 2B). For example, with the Flu M1 peptide, the anisotropy ranged from 128 to 201 mA for probes on the ␣1 helix and from 88 to 112 mA for probes on the ␣2 helix. To help place the numbers in perspective, the theoretical maximum anisotropy for a fully constrained molecule is 400 mA, and we measured a value of 9 mA for unconjugated Alexa Fluor 488. Measurements on the denatured samples ranged from 40 to 64 mA. The measurements in Fig. 2 were reproducible and independent of concentration over a 100-fold range (10 nM to 1 M). The general variations in anisotropy with peptide and position were maintained at 4 and 35°C. The independence of the results with both concentration and temperature further rules out any influence of complex instability.
More striking than the positional variation in the anisotropy was how it varied with different peptides. These data are also in Fig. 2B and recast as percentages of deviation from the mean in Fig. 2C. Although in general the ␣2 helix was more flexible than the ␣1 helix, the flexibilities at sites within both helices varied considerably. Large (ϳ2-fold) variations were seen at the center of the ␣1 helix, as well as the end of the ␣2 helix. Variations of 20 -40% were seen at the other positions. Our previous studies demonstrated that a 25-30% shift was associated with a large MHC conformational rearrangement and higher entropic penalty upon TCR binding (10). Thus, the variation in flexibility within the HLA-A2 peptide-binding groove is significant. The variation in the region of the ␤-sheet floor was also significant; of particular interest is the wide variation for position 115, an amino acid critical for recognition of HLA-A2 by the CD8 coreceptor (31).
We examined in detail the possibility that nonspecific labeling could influence the results. No nonspecific covalent labeling of wild-type HLA-A2 heavy chains was observed in parallel labeling reactions. To rule out any noncovalent nonspecific labeling, a sample of purified protein labeled at position 151 and refolded with the Tax peptide was dialyzed under denaturing conditions to remove peptide. After extensive dialysis, the sample was recovered and split into five aliquots, and the labeled HLA-A2 heavy chain in each aliquot was refolded a second time with each of the peptides studied. Measurements on these samples were essentially identical to the measurements on samples that had not undergone this process (Table 1).
We also repeated select measurements using BODIPY-FL, a fluorophore with structural and chemical properties very different from Alexa Fluor 488. Although the absolute anisotropy numbers differed as expected due to the different fluorescence lifetimes and linker arms, the variation seen with different peptides was very similar with the two probes (Table 2), ruling out any probe-specific effects.
Evidence for Differential MHC Dynamics from Crystallographic Structures-Although our HDX-MS and fluorescence anisotropy data show that different peptides alter the motional properties of HLA-A2 protein, as indicated above, the overall architecture of class I MHC proteins is relatively invariant with peptides. Nonetheless, some variation in the position and geometry of the ␣1 and ␣2 helices with peptide has been noted (e.g. Refs. 10, 12, and 32). These variations are similar to those that have been observed in comparative molecular dynamics simulations of different peptides bound to the same class I MHC protein (6,7,10).
To more precisely quantify the structural differences in the same class I MHC protein crystallized with different peptides,  we assembled a database of 51 nonredundant peptide⅐HLA-A2 structures (supplemental Table S1A) and examined how the geometry of the HLA-A2 ␣1 and ␣2 helices varies with peptide.
We computed the average bending angle of the helices (defined as the angle formed by least squared-fit lines between the (i-3)-i and i-(iϩ3) ␣-carbons of a helix) (20) as well as the width of the peptide-binding groove. The results indicated small but significant peptide-dependent variations (Fig. 3, A and B). We also compared binding groove geometry in all cases in which two peptide⅐HLA-A2 complexes were present in the crystallographic asymmetric unit. Although the database was smaller (26 complexes), the differences between the two molecules were similar to those observed in the larger dataset (supplemental Fig. S3a). Lastly, we examined cases in which structures of the same peptide⅐HLA-A2 complex were available both free and bound to a TCR. Although the size of the bound dataset was even smaller (only 12 different TCR-peptide⅐HLA-A2 structures were available, with three pairs of TCRs bound to the same pMHC), differences of the same magnitude were again observed between free and bound (supplemental Fig. S3b).
To further examine the structures, we examined the variation in crystallographic B-factors (or temperature factors). Although B-factors are impacted by model error and lattice imperfections and report more on static disorder within the crystalline environment as opposed to motions in solution (33,34), we reasoned that an analysis of a large number of pMHC complexes would allow some of these limitations to be circumvented.
The average and standard deviations for normalized B-factors of ␣-carbons of the ␣1/␣2 helices and the C-terminal half of the ␤-sheet floor of the groove are shown in Fig. 3C. In general, the average normalized B-factors fit with expected constraints on protein secondary structure, increasing toward the termini of the helices and in the turns of the ␤-sheet. Importantly however, for all three elements of structure, the variation in the average values was considerable. The standard deviations   for the values in the ␣1 helix for example ranged from 19 to 30% of the mean values, with an average of 24%. Similar results were seen for the ␣2 helix and the ␤-sheet. These results are consistent with the peptide dependence of flexibility detected by HDX-MS and fluorescence anisotropy. To assess whether our results are unique to HLA-A2, we repeated the B-factor analysis on structures of different peptides bound to the murine H-2K b class I MHC protein (supplemental Table S1B). Although the number of complexes was smaller (28 nonredundant H-2K b structures when compared with 51 for HLA-A2), the same general results were seen; the average standard deviations of the B-factors of the ␣-carbons of the H-2K b ␣1 helix, ␣2 helix, and ␤-sheet were 30, 32, and 24%, respectively. Thus, the results observed with HLA-A2 may be generalizable to other class I MHC proteins.

DISCUSSION
Studying recognition of the Tel1p and Tax peptides by the A6 TCR, we previously observed that peptide-dependent motions of the HLA-A2 peptide-binding groove could influence TCR recognition, specificity, and cross-reactivity (10). Here, we provide concrete evidence that the "tuning" of the flexibility of class I MHC proteins by different peptides is a general phenomenon and not localized to certain peptides or a single region of the protein.
The alteration of protein motions by ligand binding can be interpreted in the context of energy landscapes, which relate protein conformation and energy and define the populations and transition kinetics for the conformational substates that make up the structural ensemble of a protein (35,36). In altering flexibility, the binding of a ligand to a particular set of substates stabilizes those and destabilizes others, altering the energy landscape and thus entropy and protein dynamics (37). Importantly, the binding of different ligands can stabilize or destabilize different sets of substates, resulting in ligand-dependent changes in protein motions (38 -40). In addition to our results here, this behavior has been demonstrated for a variety of proteins, with differential motions impacting processes as diverse as signaling, enzyme activity, and allostery (reviewed in Ref. 41). In structural terms, the stabilization/destabilization of different states by various ligands can arise from formation of different types, strengths, and numbers of noncovalent interactions, as well as different steric limitations on what conformations can be sampled as the protein moves. Given that peptide binding is coupled to class I MHC folding, with the peptide deeply integrated into the structure of the protein, opportunities exist for peptides to impact the stability and thus populations of numerous states (as a simple example, when the structures of the Tax and Flu M1 peptides bound to HLA-A2 are compared, there are 11 amino acids that are uniquely contacted in the two structures (supplemental Fig. S4)). Our results are interpreted in this context in Fig. 4, which shows a stylized class I MHC folding funnel and stylized variations in the folded state energy landscape with different peptides.
It may be notable that beyond showing peptide-specific differences, the correlation between the differences measured by HDX-MS and fluorescence anisotropy is not strong. However, it is important to note that the two experiments probe very different time and length scales. The HDX-MS experiments probe motions indirectly on the millisecond timescale and more slowly, and cannot localize the motions to specific amino acids within the peptidic fragments, which range from 7 to 20 amino acids. Conversely, the fluorescence anisotropy experiments probe motions on the fast, nanosecond timescale directly at the labeled positions. Thus, correlations between the HDX-MS and fluorescence anisotropy results need not be expected. The variation across timescales, however, is important as it speaks to the potential for altered MHC flexibility to impact TCR recognition through multiple ways, including binding entropies (fast timescales), conformational changes (fast and slow timescales), and altered populations of bindingcompetent states (slow timescales).
A limitation of our experiments is that they do not provide information about what conformational states are sampled. At the nanosecond timescale, some insight can be gained from prior molecular dynamics simulations of the same class I MHC protein with different peptides bound, each of which showed peptide-dependent variations in the overall shape, positions, and stabilities of the ␣-helices of the peptide-binding groove (6,7,10). Although these studies did not address variation in the ␤-sheet floor of the binding groove, given the tight connectivity between peptides, the helices, and the ␤-sheet, a peptide dependence to motion within individual ␤-strands should not be surprising (indeed, the differences in peptide-HLA-A2 contacts in different complexes include multiple residues of the ␤-sheet (supplemental Fig. S4)). A more direct comparison between our results and these simulations is limited by the finite simulation times of the simulations (50 -400 ns); advanced sampling methods that extend the range of molecular simulations should be expected to provide further insight.
Our findings have implications for the determinants of antigenicity and the functional distinctions that are often made between a peptide and a presenting class I MHC protein. When considered alongside our previous findings with the Tax and Tel1p peptides (10), our results imply that there can be an "extension of antigenicity" from the peptide to the MHC and that distinctions commonly made between the two weaken at the atomic level (42). For example, patterns of similar TCR-MHC interatomic interactions in structures of TCRs sharing variable gene segments bound to the same MHC have been used to help argue that TCRs maintain an intrinsic affinity toward MHC proteins (43). However, if the strength of TCR-MHC contacts can vary due to altered MHC dynamics, this limits the extent to which evolution can have imparted an MHC bias onto TCR gene segments. This could help explain why clear patterns of conserved TCR-MHC contacts have been difficult to identify, even in structures where the same TCR is engaged to the same MHC protein presenting different peptides (4).
Our results are also relevant to the surprising degree of specificity seen in many TCR-pMHC interactions. Although TCR cross-reactivity is well appreciated (44), in many cases, TCRs have been shown to be highly sensitive to changes in peptide structure or composition, sometimes with no clear structural explanation. Such fine specificity has been implicated in T cell antagonism as well as the poor performance of modified peptides as vaccine candidates (45,46). Our findings suggest that in some cases, TCR specificity could arise from dynamical changes resulting from peptide modification.
We focused largely on a single human class I MHC allele, HLA-A*0201. Although the structural analysis suggests similar behavior with the murine class I MHC H-2K b , an important question to resolve is whether and how MHC micropolymorphisms influence the peptide dependence of dynamics. This may be particularly likely for polymorphic positions whose side chains are embedded within the peptide-binding groove (15). Less clear is the extent to which peptide-dependent dynamics occur in class II MHC proteins. Although peptide association with class II molecules is less of a folding reaction and peptide termini in class II are not embedded in the protein as with class I, a recent examination of a large database of class II pMHC structures revealed structural variations resembling those seen with class I MHC proteins (47), possibly reflecting intrinsic class II MHC dynamics that could be modulated by peptide.
In addition to T cell receptors, a variety of other molecules interact with MHC proteins. These include the CD4/CD8 coreceptors expressed by T cells, and for class I MHC molecules, a range of activating and inhibitory receptors expressed by natural killer and other cell types (48). Structural and binding data indicate that these molecules interact with class I MHC molecules in a way distinct from T cell receptors, making few, if any contacts to peptides (28, 29, 49 -51). Key contact regions include the short arm of the ␣2 helix and the N-terminal end of the ␣1 helix, as well as the sides or bottoms of the ␤-sheet floor of the peptide-binding groove. Our data indicate that peptides can modulate the flexibility of each of these regions, suggesting a mechanism for the paradoxical peptide dependence that has been reported for some of these interactions (e.g. Ref. 51). Thus it may be possible that the tuning of MHC dynamics by different peptides reflects a fundamental mechanism for modulating antigenicity.