- 1Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN, United States
- 2Harper Cancer Research Institute, University of Notre Dame, South Bend, IN, United States
Structural biology of peptides presented by class I and class II MHC proteins has transformed immunology, impacting our understanding of fundamental immune mechanisms and allowing researchers to rationalize immunogenicity and design novel vaccines. However, proteins are not static structures as often inferred from crystallographic structures. Their components move and breathe individually and collectively over a range of timescales. Peptides bound within MHC peptide-binding grooves are no exception and their motions have been shown to impact recognition by T cell and other receptors in ways that influence function. Furthermore, peptides tune the motions of MHC proteins themselves, which impacts recognition of peptide/MHC complexes by other proteins. Here, we review the motional properties of peptides in MHC binding grooves and discuss how peptide properties can influence MHC motions. We briefly review theoretical concepts about protein motion and highlight key data that illustrate immunological consequences. We focus primarily on class I systems due to greater availability of data, but segue into class II systems as the concepts and consequences overlap. We suggest that characterization of the dynamic “energy landscapes” of peptide/MHC complexes and the resulting functional consequences is one of the next frontiers in structural immunology.
Introduction
Presentation of peptides by class I MHC (MHC-I) proteins to T cell receptors (TCRs) is a key component of cellular immunity. The first structures of peptide/MHC-I structures illustrated how peptides were presented by MHC-I proteins, lying in an extended form embedded within a binding groove formed by two flanking α helices and a β sheet floor (1). The general architecture of peptide/MHC-I complexes is widely recognized, with representations found within thousands of reviews, research publications, and textbooks ranging from general biology to advanced immunology. The solution of peptide/MHC-I structures answered fundamental questions in immunology, including how a single receptor can simultaneously recognize both self (MHC) and non-self (peptide) in antigen recognition (2). The subsequent solution of structures of complexes between TCRs and peptide/MHC-I complexes illustrated at an atomic level how simultaneous recognition of self/non-self occurs (3, 4).
The first high-resolution crystal structure of a TCR-peptide/MHC-I complex was that of the A6 TCR bound to the HTLV-1 Tax11–19 antigen presented by HLA-A2 (3). Comparison of the TCR-bound complex with that of the free peptide/MHC-I revealed that the peptide undergoes a conformational change upon TCR binding, centered around the central core—the peptide is essentially “squished” into the binding groove by the receptor (5) (Figure 1A). Many subsequent structures demonstrated that peptide conformational changes frequently occur upon TCR binding (Figure 1B). The regularity of this occurrence is likely related to the fact that in MHC-I complexes peptides are usually not flat within the binding groove, but bulge due to the tethering of the N- and C-terminal residues, with the degree of bulging increasing with peptide length (6–9). Peptides can also be “pulled” away from the binding groove in response to TCR binding, an occurrence which could be related to the identity of the primary anchor residues (10, 11).
Figure 1. Conformational changes in peptides and heterogeneity in peptides bound to class I MHC (MHC-I). (A) Illustration of the change in the backbone of the Tax11–19 peptide bound to HLA-A2 upon binding of the A6 T cell receptors (TCR). The conformational change is centered upon amino acids 6 and 7, with a maximal displacement of 3 Å occurring at the α carbon of position 6. The root mean square (RMS) deviation for contiguous peptide backbone atoms between TCR-free and TCR-bound is 1.3 Å. (B) Statistics of peptide conformational changes that occur upon TCR binding for all peptide/MHC-I complexes for which TCR-free/bound structures exist in the Protein Data Bank. The figure shows a Box plot of peptide backbone RMS deviations between free and bound structures. Individual values are indicated by red dots, and the interquartile range between the 75th and 25th percentiles indicated in yellow. Whiskers extend to the furthest values that lie within the 75th and 25th percentile value ±1.5× the interquartile range. RMS deviations are binned according to the Freedman–Diaconis rule for a total of 12 bins (12). The blue dashed curve indicates the population distribution. Complexes and TCR-free/bound PDB codes are provided for the 25th, 50th, and 75th percentile values, as well as for the points that demarcate the whiskers and for the system which displayed the largest peptide conformational change upon binding. (C) Weak electron density for the triply modified 10-mer peptide GP2 bound to HLA-A2 (top) and the 11-mer peptide BZLF1 bound to HLA-B35 (bottom) (13, 14). Gaps in the density show regions where conformational heterogeneity is likely to exist. Density is calculated from a 2Fo−Fc map and contoured at 1σ. (D) Multiple conformations of the anchor modified MART-127–35 ALG nonameric peptide bound to HLA-A2 (15). The electron density was sufficiently clear to allow refinement of the backbone at positions 4 and 5 in two different conformations (blue and gold regions).
In some crystal structures of peptide/MHC-I complexes, peptides are poorly refined in the binding groove, with side chains and even backbones lacking electron density (Figure 1C). There can be multiple reasons for weak or missing electron density in protein structures, such as poor crystal morphology or even the existence of multiple peptides in one crystal, as was the case in the very first structure of HLA-A2 (16). Another reason for poor electron density is structural heterogeneity, stemming from the existence of multiple peptide conformations or the interconversion between different conformations on the timescale of the X-ray diffraction experiment. Such heterogeneity was demonstrated in an early experiment with the GP2 HER-2/neu epitope, which had missing density in the peptide center when bound to HLA-A2 (17). In other cases, electron density is clear enough to identify peptides in multiple conformations (15) (Figure 1D).
What do binding-induced conformational changes and structural heterogeneity have in common? They are both indicative of molecular motion, or protein dynamics. Proteins are not static molecules, and their atoms move independently and collectively over a wide range of timescales. Peptide/MHC complexes are no exception, and peptide conformational changes, weak electron density, and structural heterogeneity give some insight into the influence that motion can have in antigen recognition. Below, we outline how peptide motion can be important in influencing antigenicity and suggest how it may be considered in efforts to predict and even manipulate immune recognition.
Theoretical Considerations of Protein Motion
Although not long considered in structural and molecular immunology, the study of protein dynamics has a long and rich history. The first computational studies of protein motion were performed in the 1970s [e.g., Ref. (18)] and experimental studies using various forms of spectroscopy [fluorescence, nuclear magnetic resonance (NMR)] date to the 1960s (19). The question of how protein motion drives biological function was considered as early as 1972 in discussions about conformational changes occurring in hemoglobin as it performs its physiological functions (20).
Often, protein dynamics that influence molecular recognition are described using terms such as induced fit, preexisting equilibrium, or conformational selection. Induced fit is commonly used to describe structural differences that are apparent between free and bound structures, whereas preexisting equilibrium or conformational selection is often used to describe rapid motion that occurs in an un-bound protein. Although used to describe distinct scenarios, each of these terms ultimately reflects the influence of motion (21, 22). The distinctions boil down to timescales, or the rates at which proteins exhibit flexibility and undergo conformational changes. Timescales and rates in return relate to the heights of free energy barriers between conformations and the frequencies with which these barriers are overcome. In that sense, classical induced fit and preexisting equilibria/conformational selection mechanisms reflect extremes on a continuum—slow, or low frequency motion with high barriers for the former, and rapid, high frequency motion with low barriers for the latter (Figure 2A). The barriers between different structural states can be modulated by changes in the molecular environment, contributing to “induced fit” changes occurring during binding—but nonetheless, induced fit inherently reflects a propensity to move. In some cases, assignment of structural changes to either induced fit or conformational selection can be clear, such as seen with recognition of peptides by the A6 and G10 TCRs (23, 24). More often though, when induced fit, structural heterogeneity or poor electron density is indicated from crystal structures, the heights of the barriers can only be guessed at, so where on the continuum the associated protein motions lie is often unknown.
Figure 2. Free energy landscapes and peptide tuning of the class I MHC (MHC-I) energy landscape. (A) Schematic showing the variation in protein-free energy with conformation. Atomic coordinates are along the x axis and free energy is along the y axis. Wells along the conformation axis represent different structural substates separated by barriers. The heights of the barriers yield the rates at which the protein moves between (or samples) different conformations. Low barriers translate into rapid, high-frequency motions, whereas high barriers translate into slower, low-frequency motions. The number and energy levels of the structural substates separated by the barriers gives the protein entropy. (B) Peptide tuning of the MHC-I energy landscape. Ligplot analysis (25) of the interactions the Tax and Flu M1 peptides make with HLA-A2. Peptides are shown with purple bonds and contacting MHC residues are indicated. Interactions formed by hydrophobic atoms are shown with red hashes. Hydrogen bonds are shown with green lines with distances indicated. HLA-A2 residues making significant interactions in one complex but not the other are circled in red. (C) Illustration of how peptides can tune the MHC-I energy landscape. The image shows a traditional folding funnel, with the native peptide/MHC-I architecture at the bottom of the funnel. Zooming into the tip of the funnel reveals the energy landscape of the assembled complex, as diagrammed in panel (A). Due to the different interactions formed by different peptides as shown in panel (B), the energy landscape is altered, changing MHC-I protein dynamics. Figure adapted from Ref. (26) and used by permission.
How can protein motions—fast, slow, or intermediate—influence molecular recognition, such as occurs between a TCR and peptide/MHC complex? To form the most stable complexes (i.e., those with the lowest free energy), proteins need to optimize shape and chemical complementarity within the protein–protein interface. If motions are required to achieve this, then these motions can limit the rate at which complexes form, in turn limiting the overall affinity of the complex. An example is found in the antibody maturation process, during which mutations are introduced that remove conformational heterogeneity and its associated motion from the antibody binding site (27). Accordingly, one signature of antibody maturation is an increase in the rate of association rates as maturation proceeds, strengthening binding (28).
A related mechanism is associated with the population of multiple conformational states—the various “wells” between the barriers in Figure 2A. Conformational heterogeneity is directly related to entropy: the more states that are populated, the greater the entropy. If heterogeneity is reduced upon binding as is usually expected, then structural heterogeneity will increase the entropic cost for binding, and affinity will weaken. Again using antibodies as an example, affinity maturation is associated with a reduction in antibody conformational heterogeneity and lower entropic costs for binding (27).
Consider TCR recognition of a peptide/MHC-I complex in the context of the above discussion. A peptide that must move to optimize fit with a TCR, such as the Tax11–19 antigen presented by HLA-A2 in Figure 1A, will result in slower binding compared to a more rigid peptide presented in an optimal conformation. Indeed, recognition of the Tax11–19/HLA-A2 complex by the A6 and B7 TCRs proceeds with association rates in the range of ~0.5 × 106 M−1 s−1, orders of magnitude slower than the diffusion controlled limit (29). Other TCR binding reactions associated with peptide conformational changes occur with similar rates, or even slower (although it is important to note that TCR motional properties can also influence association rates, confounding the assignment of slow association rates purely to peptide conformational changes) (30). Clear demonstration of entropic penalties associated with peptide conformational changes or motion are more difficult to find given the varied contributions to binding entropy changes, but many TCRs bind with unfavorable entropy changes (31).
Assessing Peptide Motions in MHC-I Binding Grooves
As discussed above, indications of peptide motions within the MHC-I binding groove can come from crystallographic structures. This can include the presence of multiple peptide conformations, missing electron density, or significant conformational changes occurring upon TCR binding. However, moving from structural indications of motion to more detailed, actual assessments of motion requires additional experiments.
The simplest approach is to use crystal structures as the basis for molecular dynamics (MD) simulations. MD simulations take the initial set of atomic coordinates and compute the time-dependent variations in structure using classical laws of motion and a “force field” that describes the interactions between atoms. With ever-growing improvements in computer hardware, MD simulations have advanced considerably, such that we are now seeing simulations of entire protein folding reactions on timescales extending beyond microseconds (32). MD simulations have been used extensively to study the motions of peptides bound to MHC-I proteins (15, 33–51), as reviewed recently by Freund and colleagues (52). MD simulations have the advantage of tracking actual motion in atomic detail, capturing different conformations as a function of time. In most cases, however, MD simulations are limited to capturing motions that occur on fast timescales, and of course are virtual, with confidence in the results depending on the quality of the initial structure and a wide range of parameters that can be tuned to optimize speed versus accuracy.
Direct experimental measurements of peptide motions are less common than MD simulations, owing to the challenges of producing and working with peptide/MHC-I complexes. NMR spectroscopy is the gold standard for experimentally monitoring protein motion in atomic detail. NMR can provide motional detail across fast (picosecond to nanosecond) and slow (millisecond or greater) timescales and in some cases, can give insight into structural aspects of motion (e.g., amplitudes and directions). However, NMR on proteins the size and complexity of peptide/MHC complexes is technically challenging and of low throughput. Correspondingly, NMR has seen less use in studying peptide/MHC-I complexes (53–58). In two recent NMR studies, peptides bound to MHC-I proteins were shown to interconvert between conformations seen in different crystallographic structures (53, 56), confirming the expectation that crystallographically observed conformational changes reflect the propensity to move.
Fluorescence spectroscopy is commonly used to study protein motion, and it has seen some use in studying peptide/MHC-I complexes (35, 45, 59–62). In one notable study, fluorescence polarization was used along with MD simulations to show that allelic variations in HLA-B27 led to differences in peptide motions (35). A limitation of fluorescence is a reliance on bulky fluorescent labels, which can perturb structure and dynamics (63) [though in at least one case, intrinsic tryptophan fluorescence has been used to study MHC-I behavior (62)]. Infrared spectroscopy is another approach that may prove to be of use in studying the motions of peptides bound to MHC proteins, as it can report on fast motions using approaches analogous to NMR, but with reporters that are less perturbative—and sometimes non-perturbative in the case of C-D probes—to the native structure and dynamics of interest (64, 65).
Peptide Motions and Antigenicity: Modified Peptides and Tumor Neoepitopes
From the discussion thus far, it is evident that motion of peptides in MHC-I binding grooves will impact TCR binding, influencing association rates, entropic penalties, or both. Although exceptions are known, in general peptide antigenicity scales with the affinity of the TCR for the peptide/MHC complex, with thresholds at the high and low limits (30, 66, 67). All other factors being equal, increasingly mobile peptides will be recognized more weakly by TCRs, with correspondingly reduced antigenicity.
The clearest example of the influence of peptide motion on antigenicity comes from the family of MART-1 tumor antigens and associated variants. The MART-127–35 nonamer (sequence AAGIGILTV) is weakly immunogenic, attributed to its weak binding to HLA-A2 (68, 69). Weak binding stems from the presence of a suboptimal alanine at the first primary anchor position (peptide position 2). Modification of peptide primary anchors is a well-known strategy for improving peptide binding, and in some cases, antigenicity (70). Curiously, replacing the second alanine of the MART-127–35 nonamer with leucine eliminates antigenicity with multiple T cell clones, despite strengthening MHC binding. The crystal structure of the anchor modified variant bound to HLA-A2 showed multiple conformations of the peptide in the binding groove, compared to a single conformation with the native, unmodified peptide (15). NMR confirmed the crystallographic data, showing the modified peptide indeed sampled multiple conformations (38). Exploration of peptide motions using MD simulations suggested the multiple conformations were attributable to rapid, nanosecond motions in the backbone of the modified peptide that were absent from the native peptide (comfortingly, the MD simulations recapitulated the conformations that were seen crystallographically). Binding experiments with recombinant TCRs confirmed the more dynamic modified peptide was in fact recognized more weakly, explaining its reduced antigenicity.
The results with the MART-1 nonamer suggest that screening for changes in peptide motion may be useful in helping predict antigenicity. One area where this may be particularly helpful is in predicting the antigenicity of “neoepitopes” arising from mutations present in tumor genomes. Neoepitopes are of considerable interest in cancer immunotherapy (71). However, predictive algorithms historically used in identifying immunodominant epitopes from viral genomes have performed poorly in predicting neoepitope antigenicity (37, 72). Toward this end, we used structural modeling and MD simulations to investigate a small set of neoepitopes and their wild-type counterparts bound to the mouse MHC-I protein H-2Kd. Although only a small number of peptides were examined, there was a positive relationship between mutations which led to reduced peptide motions and antigenicity (37). The relationship between peptide rigidity and immunogenicity is consistent with the results seen with the MART-1 nonamer, as well as the theoretical considerations noted above: mutations that reduce motion should enhance T cell recognition by increasing TCR association rates and decreasing entropic costs for binding. Other factors are undoubtedly important to consider in neoantigen immunogenicity (peptide processing, MHC binding, amino acid composition, etc.), but incorporating structural modeling and predictions of how mutations alter fast peptide dynamics may provide another layer of sophistication for prediction methods and increase the likelihood of identifying those neoepitopes most likely to induce tumor rejection.
Extension of Peptide Dynamics into MHC-I Proteins
Despite the common illustrations that render peptides and MHC-I proteins as distinct components (e.g., Figure 1), peptides are usually closely packed within MHC-I binding grooves [excluding those that are unusually long and bulge extensively (7, 8, 73–75)]. Therefore, it should not be surprising that peptide features can influence features of the peptide-binding groove. The possibilities for peptide-dependent structural shifts in MHC-I α helices were first noted in 1996 (76). More recently, we performed a comprehensive analysis of 51 structures of peptide/HLA-A2 complexes and found systematic deviations in the width of the peptide-binding groove and the bends and positions of both the α1 and α2 helices (26).
If there are peptide-dependent structural MHC-I shifts, there must be peptide-dependent MHC-I motions. Indeed, this was observed in 2009 when recognition of the Tax and Tel1p peptides by the A6 TCR was shown to yield different conformations in the α2 helix of HLA-A2 (34). It was shown that the α2 helix moved differently on the nanosecond timescale with the two peptides bound, which contributed to the structural shift seen with TCR binding to Tel1p versus Tax. The motions underlying the helix conformational change fit were found to be dependent on peptide structural features, and were correlated with differential peptide dynamics. Importantly, the different positioning of the MHC α2 helix was necessary for optimal interactions with the TCR, demonstrating that the “tuning” of HLA-A2 motions by the Tax and Tel1p proteins was indeed functionally important. Peptide-dependent alterations in the properties of the MHC-I α2 helix have been seen in other cases (77), suggesting that this region of the protein may be particularly sensitive to features of different peptides.
Returning to the MART-1 peptide system, via MD simulations, alteration of the MART-127–35 nonamer was shown to alter the fluctuations of the α1 and α2 helices (38). More recently, fluorescence anisotropy and hydrogen/deuterium exchange was used to show that peptide-dependent tuning of MHC-I protein dynamics is a general phenomenon, not limited to particular peptides or a single site in the peptide-binding domain (26). Further evidence has been seen in a NMR study of HLA-B35, where it was shown that different peptides altered MHC sampling of minor conformations that had not been observed crystallographically (58). A recent analysis of temperature factors in peptide/MHC-I structures provided additional data and highlighted regions which might be particularly susceptible to include the α1 and α2 helices, as well as regions in the α3 domain and even β2-microglobulin (52). Thus, there is a structural and dynamic “extension of antigenicity” from the peptide to the MHC protein, with significant potential to influence TCR and T cell recognition. There is also evidence that allele-specific variations in different MHC-I proteins can modulate the susceptibility for influencing protein motions (35, 36, 44, 47–49, 78), suggesting a previously unrecognized complexity in how MHC genetics can influence antigenicity.
How can different peptides modulate MHC protein dynamics? Although they may be structurally similar, at the atomic level different peptides will form many different interactions with an MHC protein, as shown in Figure 2B. The different interatomic interactions will be of disparate energies; even a hydrogen bond between two identical atoms will be of different strengths if the geometries differ (26, 76). Thus, the free energies of structural substates will vary with different peptides bound. This is diagrammed in Figure 2C, which shows a peptide/MHC “folding funnel” leading from unfolded/unassembled protein to the assembled peptide/MHC-I native state. Different peptides will tune the energy landscape of the assembled complex as discussed above. In other words, the levels of the wells describing different structural states will move up or down depending on the identity of the peptide. This will change the heights of the barriers between them, altering the rates of conformational exchange, i.e., changing protein motions. Moreover, with different energies for different substates, the overall entropy of the protein will change, potentially leading to different entropy changes upon TCR binding and impacting binding affinity. Note that as crystallographic structures do not show energy and do not provide clear insights into structural substates, alterations in peptide/MHC dynamics can occur even in the absence of any apparent crystallographic differences (58). This same phenomenon explains how allelic variations can change peptide–protein dynamics: amino acid differences between MHC variants result in altered protein–peptide contacts, affecting the free energies of structural substates and with corresponding impacts on protein dynamics.
Peptide-dependent protein motion could also explain observations of peptide-dependent binding of other proteins to MHC-I, such as NK receptors (79). Indeed, it is believed that alterations in protein fluctuations is a key mechanistic component of peptide loading and exchange, perhaps in conjunction with partial peptide dissociation and the corresponding alterations in MHC dynamics (42, 44–46, 48). In classic biochemical terms, peptides are acting as allosteric effectors, modifying the binding of these other proteins and subsequent immunological functions. The underlying mechanisms by which peptides allosterically alter motions at remote sites is not clear, but could involve either discrete pathways of motion or more global alterations of the protein energy landscape and subsequent dynamics (80).
Dynamics in Class II MHC (MHC-II) Proteins
The discussion above has centered on peptide motion within MHC-I proteins. Are there parallels with class II MHC proteins? The theoretical possibility of course exists, and the biochemical and functional implications are the same: peptide motion in MHC-II binding grooves can influence the binding of TCRs and other receptors, and peptides and MHC-II allelic variations have potential to influence MHC-II protein motions in functionally significant ways. MHC-I and MHC-II proteins are structurally homologous. However, in MHC-II, the peptide termini extend from the binding groove and, therefore, MHC-II-bound peptides are more extended, lack the bulges seen with peptides bound to MHC-I, and are more hydrogen-bonded to the protein (1). Likely for these reasons, large-scale peptide conformational changes upon TCR binding are seen less frequently with MHC-II systems (although there is less structural data for MHC-II). However, TCRs also bind peptide/MHC-II complexes with slow kinetics and occasionally unfavorable entropy changes (30, 31). Hairpin-style secondary structures have been seen in class II-presented peptides (81) and flanking regions of class II-presented peptides that lie outside of the groove can impact TCR binding, possibly due to dynamic effects (82). Therefore, although the data are less apparent, motions of peptides bound to MHC-II systems should also be considered when taking stock of the physical influences on TCR binding and peptide antigenicity. Furthermore, there is clear evidence for peptide influences on the motion of the MHC-II molecule, particularly in regions that interact with the peptide-exchange catalyst HLA-DM, as well as a fundamental role for protein motion in MHC-II peptide-exchange itself (83–86). Thus, in addition to having converged on fundamentally similar structures and functions, MHC-I and MHC-II protein functions may be similarly influenced by—and take advantage of—peptide and protein motional properties.
Conclusion
Structural biology of peptides presented by MHC proteins has transformed our understanding of immunology, with impacts ranging from our understanding of fundamental immune mechanisms to the design and optimization of vaccines. However, considering peptide/MHC complexes not as static structures but as molecules that move and breathe has opened new avenues of investigation and shed new light on factors that impact immunogenicity. The mobility of peptides within MHC protein binding grooves can impact antigen immunogenicity, regardless of whether mobility or conformational variability is apparent from crystallographic structures. This is relevant in vaccine design, as well as predicting immunogenic epitopes from pathogen and tumor genomes. Efforts to predict peptide mobility may be particularly helpful in screening cancer neoepitopes, which have proven particularly challenging for predictive immunology. Peptides also “tune” the mobility of MHC proteins themselves, contributing to antigenicity and affecting interactions with other proteins. The latter can include the machinery of peptide loading and exchange as well as other activating and inhibitory immune receptors. Broader characterization of the energy landscapes of peptide/MHC complexes and the resulting functional consequences is clearly one of the next frontiers in structural immunology.
Author Contributions
The manuscript was developed and written by all three authors.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Funding
This study was supported by grant R35GM118166 from the National Institutes of General Medical Sciences, United States National Institutes of Health.
References
1. Madden DR. The three-dimensional structure of peptide-MHC complexes. Annu Rev Immunol (1995) 13(1):587–622. doi:10.1146/annurev.iy.13.040195.003103
2. Bevan MJ. Altered self, altered world. J Immunol (2004) 173(5):2897–8. doi:10.4049/jimmunol.173.5.2897
3. Garboczi DN, Ghosh P, Utz U, Fan QR, Biddison WE, Wiley DC. Structure of the complex between human T-cell receptor, viral peptide and HLA-A2. Nature (1996) 384(6605):134–41. doi:10.1038/384134a0
4. Garcia KC, Degano M, Stanfield RL, Brunmark A, Jackson MR, Peterson PA, et al. An alphabeta T cell receptor structure at 2.5 A and its orientation in the TCR-MHC complex [see comments]. Science (1996) 274(5285):209–19. doi:10.1126/science.274.5285.209
5. Madden DR, Garboczi DN, Wiley DC. The antigenic identity of peptide-MHC complexes: a comparison of the conformations of five viral peptides presented by HLA-A2 [published erratum appears in Cell 1994 Jan 28;76(2):following 410]. Cell (1993) 75(4):693–708. doi:10.1016/0092-8674(93)90490-H
6. Tynan FE, Burrows SR, Buckle AM, Clements CS, Borg NA, Miles JJ, et al. T cell receptor recognition of a ‘super-bulged’ major histocompatibility complex class I-bound peptide. Nat Immunol (2005) 6(11):1114–22. doi:10.1038/ni1257
7. Speir JA, Stevens J, Joly E, Butcher GW, Wilson IA. Two different, highly exposed, bulged structures for an unusually long peptide bound to rat MHC class I RT1-Aa. Immunity (2001) 14(1):81–92. doi:10.1016/S1074-7613(01)00091-7
8. Ekeruche-Makinde J, Miles JJ, van den Berg HA, Skowera A, Cole DK, Dolton G, et al. Peptide length determines the outcome of TCR/peptide-MHCI engagement. Blood (2013) 121(7):1112–23. doi:10.1182/blood-2012-06-437202
9. Guo HC, Jardetzky TS, Garrett TP, Lane WS, Strominger JL, Wiley DC. Different length peptides bind to HLA-Aw68 similarly at their ends but bulge out in the middle. Nature (1992) 360:364–6. doi:10.1038/360364a0
10. Wang Y, Singh NK, Spear TT, Hellman LM, Piepenbrink KH, McMahan RH, et al. How an alloreactive T-cell receptor achieves peptide and MHC specificity. Proc Natl Acad Sci U S A (2017) 114(24):E4792–801. doi:10.1073/pnas.1700459114
11. Madura F, Rizkallah PJ, Holland CJ, Fuller A, Bulek A, Godkin AJ, et al. Structural basis for ineffective T-cell responses to MHC anchor residue-improved “heteroclitic” peptides. Eur J Immunol (2015) 45(2):584–91. doi:10.1002/eji.201445114
12. Freedman D, Diaconis P. On the histogram as a density estimator: L2 theory. Probab Theory Relat Fields (1981) 57(4):453–76.
13. Sharma AK, Kuhns JJ, Yan S, Friedline RH, Long B, Tisch R, et al. Class I major histocompatibility complex anchor substitutions alter the conformation of T cell receptor contacts. J Biol Chem (2001) 276(24):21443–9. doi:10.1074/jbc.M010791200
14. Miles JJ, Elhassen D, Borg NA, Silins SL, Tynan FE, Burrows JM, et al. CTL recognition of a bulged viral peptide involves biased TCR selection. J Immunol (2005) 175(6):3826–34. doi:10.4049/jimmunol.175.6.3826
15. Borbulevych OY, Insaidoo FK, Baxter TK, Powell DJ Jr, Johnson LA, Restifo NP, et al. Structures of MART-1(26/27-35) peptide/HLA-A2 complexes reveal a remarkable disconnect between antigen structural homology and T cell recognition. J Mol Biol (2007) 372(5):1123–36. doi:10.1016/j.jmb.2007.07.025
16. Bjorkman PJ, Saper MA, Samraoui B, Bennett WS, Strominger JL, Wiley DC. Structure of the human class I histocompatibility antigen, HLA-A2. Nature (1987) 329(6139):506–12. doi:10.1038/329506a0
17. Kuhns JJ, Batalia MA, Yan S, Collins EJ. Poor binding of a HER-2/neu epitope (GP2) to HLA-A2.1 is due to a lack of interactions with the center of the peptide. J Biol Chem (1999) 274(51):36422–7. doi:10.1074/jbc.274.51.36422
18. Karplus M, McCammon JA. Protein structural fluctuations during a period of 100 ps. Nature (1979) 277(5697):578–578. doi:10.1038/277578a0
19. Stryer L. Fluorescence spectroscopy of proteins. Science (1968) 162(3853):526–33. doi:10.1126/science.162.3853.526
20. Moon RB, Richards JH. Conformational studies of various hemoglobins by natural-abundance (13)C NMR spectroscopy. Proc Natl Acad Sci U S A (1972) 69(8):2193–7. doi:10.1073/pnas.69.8.2193
21. Csermely P, Palotai R, Nussinov R. Induced fit, conformational selection and independent dynamic segments: an extended view of binding events. Trends Biochem Sci (2010) 35(10):539–46. doi:10.1016/j.tibs.2010.04.009
22. Boehr DD, Nussinov R, Wright PE. The role of dynamic conformational ensembles in biomolecular recognition. Nat Chem Biol (2009) 5(11):789–96. doi:10.1038/nchembio.232
23. Scott DR, Borbulevych OY, Piepenbrink KH, Corcelli SA, Baker BM. Disparate degrees of hypervariable loop flexibility control T-cell receptor cross-reactivity, specificity, and binding mechanism. J Mol Biol (2011) 414(3):385–400. doi:10.1016/j.jmb.2011.10.006
24. Lee JK, Stewart-Jones G, Dong T, Harlos K, Di Gleria K, Dorrell L, et al. T cell cross-reactivity and conformational changes during TCR engagement. J Exp Med (2004) 200(11):1455–66. doi:10.1084/jem.20041251
25. Laskowski RA, Swindells MB. LigPlot+: multiple ligand–protein interaction diagrams for drug discovery. J Chem Inf Model (2011) 51(10):2778–86. doi:10.1021/ci200227u
26. Hawse WF, Gloor BE, Ayres CM, Kho K, Nuter E, Baker BM. Peptide modulation of class I major histocompatibility complex protein molecular flexibility and the implications for immune recognition. J Biol Chem (2013) 288(34):24372–81. doi:10.1074/jbc.M113.490664
27. Schmidt AG, Xu H, Khan AR, O’Donnell T, Khurana S, King LR, et al. Preconfiguration of the antigen-binding site during affinity maturation of a broadly neutralizing influenza virus antibody. Proc Natl Acad Sci U S A (2013) 110(1):264–9. doi:10.1073/pnas.1218256109
28. Foote J, Milstein C. Kinetic maturation of an immune response. Nature (1991) 352(6335):530–2. doi:10.1038/352530a0
29. Davis-Harrison RL, Armstrong KM, Baker BM. Two different T cell receptors use different thermodynamic strategies to recognize the same peptide/MHC ligand. J Mol Biol (2005) 346(2):533–50. doi:10.1016/j.jmb.2004.11.063
30. Stone JD, Chervin AS, Kranz DM. T-cell receptor binding affinities and kinetics: impact on T-cell activity and specificity. Immunology (2009) 126(2):165–76. doi:10.1111/j.1365-2567.2008.03015.x
31. Armstrong KM, Insaidoo FK, Baker BM. Thermodynamics of T-cell receptor-peptide/MHC interactions: progress and opportunities. J Mol Recognit (2008) 21:275–87. doi:10.1002/jmr.896
32. Lindorff-Larsen K, Piana S, Dror RO, Shaw DE. How fast-folding proteins fold. Science (2011) 334(6055):517–20. doi:10.1126/science.1208351
33. Zacharias M, Springer S. Conformational flexibility of the MHC class I {alpha}1-{alpha}2 domain in peptide bound and free states: a molecular dynamics simulation study. Biophys J (2004) 87(4):2203–14. doi:10.1529/biophysj.104.044743
34. Borbulevych OY, Piepenbrink KH, Gloor BE, Scott DR, Sommese RF, Cole DK, et al. T cell receptor cross-reactivity directed by antigen-dependent tuning of peptide-MHC molecular flexibility. Immunity (2009) 31(6):885–96. doi:10.1016/j.immuni.2009.11.003
35. Pohlmann T, Bockmann RA, Grubmuller H, Uchanska-Ziegler B, Ziegler A, Alexiev U. Differential peptide dynamics is linked to MHC polymorphism. J Biol Chem (2004) 279:28197–201. doi:10.1074/jbc.C400128200
36. Narzi D, Becker CM, Fiorillo MT, Uchanska-Ziegler B, Ziegler A, Böckmann RA. Dynamical characterization of two differentially disease associated MHC class I proteins in complex with viral and self-peptides. J Mol Biol (2012) 415(2):429–42. doi:10.1016/j.jmb.2011.11.021
37. Duan F, Duitama J, Al Seesi S, Ayres CM, Corcelli SA, Pawashe AP, et al. Genomic and bioinformatic profiling of mutational neoepitopes reveals new rules to predict anticancer immunogenicity. J Exp Med (2014) 211(11):2231–48. doi:10.1084/jem.20141308
38. Insaidoo FK, Borbulevych OY, Hossain M, Santhanagopolan SM, Baxter TK, Baker BM. Loss of T cell antigen recognition arising from changes in peptide and major histocompatibility complex protein flexibility: implications for vaccine design. J Biol Chem (2011) 286(46):40163–73. doi:10.1074/jbc.M111.283564
39. Sieker F, Straatsma TP, Springer S, Zacharias M. Differential tapasin dependence of MHC class I molecules correlates with conformational changes upon peptide dissociation: a molecular dynamics simulation study. Mol Immunol (2008) 45(14):3714–22. doi:10.1016/j.molimm.2008.06.009
40. Saini SK, Ostermeir K, Ramnarayan VR, Schuster H, Zacharias M, Springer S. Dipeptides promote folding and peptide binding of MHC class I molecules. Proc Natl Acad Sci U S A (2013) 110(38):15383–8. doi:10.1073/pnas.1308672110
41. Abualrous ET, Saini SK, Ramnarayan VR, Ilca FT, Zacharias M, Springer S. The carboxy terminus of the ligand peptide determines the stability of the MHC class I molecule H-2Kb: a combined molecular dynamics and experimental study. PLoS One (2015) 10(8):e0135421. doi:10.1371/journal.pone.0135421
42. Garstka MA, Fritzsche S, Lenart I, Hein Z, Jankevicius G, Boyle LH, et al. Tapasin dependence of major histocompatibility complex class I molecules correlates with their conformational flexibility. FASEB J (2011) 25(11):3989–98. doi:10.1096/fj.11-190249
43. Hein Z, Uchtenhagen H, Abualrous ET, Saini SK, Janßen L, Van Hateren A, et al. Peptide-independent stabilization of MHC class I molecules breaches cellular quality control. J Cell Sci (2014) 127(13):2885–97. doi:10.1242/jcs.145334
44. Abualrous ET, Fritzsche S, Hein Z, Al-Balushi MS, Reinink P, Boyle LH, et al. F pocket flexibility influences the tapasin dependence of two differentially disease-associated MHC class I proteins. Eur J Immunol (2015) 45(4):1248–57. doi:10.1002/eji.201445307
45. Fleischmann G, Fisette O, Thomas C, Wieneke R, Tumulka F, Schneeweiss C, et al. Mechanistic basis for epitope proofreading in the peptide-loading complex. J Immunol (2015) 195(9):4503–13. doi:10.4049/jimmunol.1501515
46. Fisette O, Wingbermühle S, Tampé R, Schäfer LV. Molecular mechanism of peptide editing in the tapasin–MHC I complex. Sci Rep (2016) 6:19085. doi:10.1038/srep19085
47. Bailey A, van Hateren A, Elliott T, Werner JM. Two polymorphisms facilitate differences in plasticity between two chicken major histocompatibility complex class I proteins. PLoS One (2014) 9(2):e89657. doi:10.1371/journal.pone.0089657
48. Bailey A, Dalchau N, Carter R, Emmott S, Phillips A, Werner JM, et al. Selector function of MHC I molecules is determined by protein plasticity. Sci Rep (2015) 5:14928. doi:10.1038/srep14928
49. Nurzia E, Narzi D, Cauli A, Mathieu A, Tedeschi V, Caristi S, et al. Interaction pattern of Arg 62 in the A-pocket of differentially disease-associated HLA-B27 subtypes suggests distinct TCR binding modes. PLoS One (2012) 7(3):e32865. doi:10.1371/journal.pone.0032865
50. Mage MG, Dolan MA, Wang R, Boyd LF, Revilleza MJ, Robinson H, et al. The peptide-receptive transition state of MHC class I molecules: insight from structure and molecular dynamics. J Immunol (2012) 189(3):1391–9. doi:10.4049/jimmunol.1200831
51. Fisette O, Wingbermühle S, Schäfer LV. Partial dissociation of truncated peptides influences the structural dynamics of the MHCI binding groove. Front Immunol (2017) 8:408. doi:10.3389/fimmu.2017.00408
52. Wieczorek M, Abualrous ET, Sticht J, Álvaro-Benito M, Stolzenberg S, Noé F, et al. Major histocompatibility complex (MHC) class I and MHC class II proteins: conformational plasticity in antigen presentation. Front Immunol (2017) 8:292. doi:10.3389/fimmu.2017.00292
53. Insaidoo FK, Zajicek J, Baker BM. A general and efficient approach for NMR studies of peptide dynamics in class I MHC peptide binding grooves. Biochemistry (2009) 48(41):9708–10. doi:10.1021/bi9008787
54. Nakagawa M, Chiba-Kamoshida K, Udaka K, Nakanishi H. NMR study on the interaction between MHC class I protein and its antigen peptide. Biochem Biophys Res Commun (2000) 278(3):609–13. doi:10.1006/bbrc.2000.3844
55. Beerbaum M, Ballaschk M, Erdmann N, Schnick C, Diehl A, Uchanska-Ziegler B, et al. NMR spectroscopy reveals unexpected structural variation at the protein–protein interface in MHC class I molecules. J Biomol NMR (2013) 57(2):167–78. doi:10.1007/s10858-013-9777-z
56. Hawse WF, De S, Greenwood AI, Nicholson LK, Zajicek J, Kovrigin EL, et al. TCR scanning of peptide/MHC through complementary matching of receptor and ligand molecular flexibility. J Immunol (2014) 192(6):2885–91. doi:10.4049/jimmunol.1302953
57. Kurimoto E, Kuroki K, Yamaguchi Y, Yagi-Utsumi M, Igaki T, Iguchi T, et al. Structural and functional mosaic nature of MHC class I molecules in their peptide-free form. Mol Immunol (2013) 55(3–4):393–9. doi:10.1016/j.molimm.2013.03.014
58. Yanaka S, Ueno T, Shi Y, Qi J, Gao GF, Tsumoto K, et al. Peptide-dependent conformational fluctuation determines the stability of the human leukocyte antigen class I complex. J Biol Chem (2014) 289(35):24680–90. doi:10.1074/jbc.M114.566174
59. Binz AK, Rodriguez RC, Biddison WE, Baker BM. Thermodynamic and kinetic analysis of a peptide-class I MHC interaction highlights the noncovalent nature and conformational dynamics of the class I heterotrimer. Biochemistry (2003) 42(17):4954–61. doi:10.1021/bi034077m
60. Gakamsky DM, Lewitzki E, Grell E, Saulquin X, Malissen B, Montero-Julian F, et al. Kinetic evidence for a ligand-binding-induced conformational transition in the T cell receptor. Proc Natl Acad Sci U S A (2007) 104(42):16639–44. doi:10.1073/pnas.0707061104
61. Gakamsky DM, Davis DM, Strominger JL, Pecht I. Assembly and dissociation of human leukocyte antigen (HLA)-A2 studied by real-time fluorescence resonance energy transfer. Biochemistry (2000) 39(36):11163–9. doi:10.1021/bi000763z
62. Saini SK, Abualrous ET, Tigan AS, Covella K, Wellbrock U, Springer S. Not all empty MHC class I molecules are molten globules: tryptophan fluorescence reveals a two-step mechanism of thermal denaturation. Mol Immunol (2013) 54(3–4):386–96. doi:10.1016/j.molimm.2013.01.004
63. Scott Daniel R, Vardeman CF II, Corcelli Steven A, Baker Brian M. Limitations of time-resolved fluorescence suggested by molecular simulations: assessing the dynamics of T cell receptor binding loops. Biophys J (2012) 103(12):2532–40. doi:10.1016/j.bpj.2012.10.037
64. Fabian H, Huser H, Narzi D, Misselwitz R, Loll B, Ziegler A, et al. HLA-B27 subtypes differentially associated with disease exhibit conformational differences in solution. J Mol Biol (2008) 376(3):798–810. doi:10.1016/j.jmb.2007.12.009
65. Le Sueur AL, Horness RE, Thielges MC. Applications of two-dimensional infrared spectroscopy. Analyst (2015) 140(13):4336–49. doi:10.1039/C5AN00558B
66. Stone JD, Kranz D. Role of T cell receptor affinity in the efficacy and specificity of adoptive T cell therapies. Front Immunol (2013) 4:244. doi:10.3389/fimmu.2013.00244
67. Hebeisen M, Oberle SG, Presotto D, Speiser DE, Zehn D, Rufer N. Molecular insights for optimizing t cell receptor specificity against cancer. Front Immunol (2013) 4:154. doi:10.3389/fimmu.2013.00154
68. Romero P, Gervois N, Schneider J, Escobar P, Valmori D, Pannetier C, et al. Cytolytic T lymphocyte recognition of the immunodominant HLA-A*0201-restricted Melan-A/MART-1 antigenic peptide in melanoma. J Immunol (1997) 159(5):2366–74.
69. Valmori D, Fonteneau JF, Lizana CM, Gervois N, Liénard D, Rimoldi D, et al. Enhanced generation of specific tumor-reactive CTL in vitro by selected melan-A/MART-1 immunodominant peptide analogues. J Immunol (1998) 160(4):1750–8.
70. Parkhurst MR, Salgaller ML, Southwood S, Robbins PF, Sette A, Rosenberg SA, et al. Improved induction of melanoma-reactive CTL with peptides from the melanoma antigen gp100 modified at HLA-A*0201-binding residues. J Immunol (1996) 157(6):2539–48.
71. Gubin MM, Artyomov MN, Mardis ER, Schreiber RD. Tumor neoantigens: building a framework for personalized cancer immunotherapy. J Clin Invest (2015) 125(9):3413–21. doi:10.1172/JCI80008
72. Schmidt J, Guillaume P, Dojcinovic D, Karbach J, Coukos G, Luescher I. In silico and cell-based analyses reveal strong divergence between prediction and observation of T cell recognized tumor antigen T cell epitopes. J Biol Chem (2017) 292(28):11840–9. doi:10.1074/jbc.M117.789511
73. Tynan FE, Borg NA, Miles JJ, Beddoe T, El-Hassen D, Silins SL, et al. High resolution structures of highly bulged viral epitopes bound to major histocompatibility complex class I: implications for T-cell receptor engagement and T-cell immunodominance. J Biol Chem (2005) 280(25):23900–9. doi:10.1074/jbc.M503060200
74. Probst-Kepper M, Hecht HJ, Herrmann H, Janke V, Ocklenburg F, Klempnauer J, et al. Conformational restraints and flexibility of 14-meric peptides in complex with HLA-B*3501. J Immunol (2004) 173(9):5610–6. doi:10.4049/jimmunol.173.9.5610
75. Hassan C, Chabrol E, Jahn L, Kester MG, de Ru AH, Drijfhout JW, et al. Naturally processed non-canonical HLA-A*02:01 presented peptides. J Biol Chem (2015) 290(5):2593–603. doi:10.1074/jbc.M114.607028
76. Smith KJ, Reid SW, Stuart DI, McMichael AJ, Jones EY, Bell JI. An altered position of the [alpha]2 helix of MHC class i is revealed by the crystal structure of HLA-B*3501. Immunity (1996) 4(3):203–13. doi:10.1016/S1074-7613(00)80429-X
77. Borbulevych OY, Piepenbrink KH, Baker BM. Conformational melding permits a conserved binding geometry in TCR recognition of foreign and self molecular mimics. J Immunol (2011) 186(5):2950–8. doi:10.4049/jimmunol.1003150
78. Hulsmeyer M, Fiorillo MT, Bettosini F, Sorrentino R, Saenger W, Ziegler A, et al. Dual, HLA-B27 subtype-dependent conformation of a self-peptide. J Exp Med (2004) 199(2):271–81. doi:10.1084/jem.20031690
79. Dam J, Guan R, Natarajan K, Dimasi N, Chlewicki LK, Kranz DM, et al. Variable MHC class I engagement by Ly49 natural killer cell receptors demonstrated by the crystal structure of Ly49C bound to H-2Kb. Nat Immunol (2003) 4(12):1213–22. doi:10.1038/ni1006
80. Motlagh HN, Wrabl JO, Li J, Hilser VJ. The ensemble nature of allostery. Nature (2014) 508(7496):331–9. doi:10.1038/nature13001
81. Zavala-Ruiz Z, Strug I, Walker BD, Norris PJ, Stern LJ. A hairpin turn in a class II MHC-bound peptide orients residues outside the binding groove for T cell recognition. Proc Natl Acad Sci U S A (2004) 101(36):13279–84. doi:10.1073/pnas.0403371101
82. Cole DK, Gallagher K, Lemercier B, Holland CJ, Junaid S, Hindley JP, et al. Modification of the carboxy-terminal flanking region of a universal influenza epitope alters CD4+ T-cell repertoire selection. Nat Commun (2012) 3:665. doi:10.1038/ncomms1665
83. Painter CA, Stern LJ. Conformational variation in structures of classical and non-classical MHCII proteins and functional implications. Immunol Rev (2012) 250(1):144–57. doi:10.1111/imr.12003
84. Ferrante A, Templeton M, Hoffman M, Castellini MJ. The thermodynamic mechanism of peptide–MHC class II complex formation is a determinant of susceptibility to HLA-DM. J Immunol (2015) 195(3):1251–61. doi:10.4049/jimmunol.1402367
85. Wieczorek M, Sticht J, Stolzenberg S, Günther S, Wehmeyer C, El Habre Z, et al. MHC class II complexes sample intermediate states along the peptide exchange pathway. Nat Commun (2016) 7:13224. doi:10.1038/ncomms13224
Keywords: MHC, peptide, flexibility, dynamics, antigenicity
Citation: Ayres CM, Corcelli SA and Baker BM (2017) Peptide and Peptide-Dependent Motions in MHC Proteins: Immunological Implications and Biophysical Underpinnings. Front. Immunol. 8:935. doi: 10.3389/fimmu.2017.00935
Received: 05 May 2017; Accepted: 21 July 2017;
Published: 07 August 2017
Edited by:
Tim Elliott, University of Southampton, United KingdomReviewed by:
Dhruv Sethi, South Asian University, IndiaShouxiong Huang, University of Cincinnati, United States
Andrew K. Sewell, Cardiff University, United Kingdom
Copyright: © 2017 Ayres, Corcelli and Baker. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Brian M. Baker, YnJpYW4tYmFrZXImI3gwMDA0MDtuZC5lZHU=