Presentation of antigenic peptides by products of the major histocompatibility complex

Molecules encoded by the major histocompatibility complex (MHC) are polymorphic integral membrane proteins adapted to the presentation of peptide fragments of foreign antigens to antigen-specific T-cells. The diversity of infectious agents to which an immune response must be mounted poses a unique problem for receptor–ligand interactions; how can proteins whose polymorphism is necessarily limited bind an array of peptides almost infinite in its complexity? Both MHC class I and class II determinants have achieved this goal by harnessing a limited number of peptide side chains to anchor the epitope in place while exploiting conserved features of peptide structure, independent of their primary sequence. While class I molecules interact predominantly with the N- and C-termini of peptides, class II determinants form an extensive hydrogen bonding network along the length of the peptide backbone. Such a strategy ensures high-affinity binding, while selectively exposing the unique features of each ligand for recognition by the T-cell receptor. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

Sub‐angstrom modeling of complexes between flexible peptides and globular proteins

A wide range of regulatory processes in the cell are mediated by flexible peptides that fold upon binding to globular proteins. Computational efforts to model these interactions are hindered by the large number of rotatable bonds in flexible peptides relative to typical ligand molecules, and the fact that different peptides assume different backbone conformations within the same binding site. In this study, we present Rosetta FlexPepDock, a novel tool for refining coarse peptide–protein models that allows significant changes in both peptide backbone and side chains. We obtain high resolution models, often of sub‐angstrom backbone quality, over an extensive and general benchmark that is based on a large nonredundant dataset of 89 peptide–protein interactions. Importantly, side chains of known binding motifs are modeled particularly well, typically with atomic accuracy. In addition, our protocol has improved modeling quality for the important application of cross docking to PDZ domains. We anticipate that the ability to create high resolution models for a wide range of peptide–protein complexes will have significant impact on structure‐based functional characterization, controlled manipulation of peptide interactions, and on peptide‐based drug design. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Design of enhanced agonists through the use of a new virtual screening method: application to peptides that bind class I major histocompatibility complex (MHC) molecules

A new screening procedure is described that uses docking calculations to design enhanced agonist peptides that bind to major histocompatibility complex (MHC) class I receptors. The screening process proceeds via single mutations of one amino acid at the positions that directly interact with the MHC receptor. The energetic and structural effects of these mutations have been studied using fragments of the original ligand that vary in length. The results of these docking studies indicate that the mutant affinity ranking of long peptides can be practically reproduced with a screening approach performed using fragments of six residues. Fragments of four and five residues could mimic, in some cases, the structural arrangement of the side chains of the full-length peptide. We have compared the structural and energetic results of the docking calculations with experimental data using three unrelated ligand peptides that differ greatly in their affinity for the MHC complex. Analysis of the affinity of the fragments led to the identification of three important parameters in the construction of fragments that mimic the structural and energetic properties of the full-length ligand: the length of the fragment; its intermolecular energy; and the number and localization, internal or terminal, of the anchor residues. The results of this new peptide-design methodology have been applied to suggest new peptides derived from the MUC1-8 peptide that could be used as murine vaccines that trigger the immune response through the MHC class I protein H-2K(b).     

Direct transport of newly synthesized HLA-DR from the trans-Golgi network to major histocompatibility complex class II containing compartments (MIICS) demonstrated using a novel tyrosine-sulfated chimera.

Binding of antigenic peptides to major histocompatibility complex (MHC) class II glycoproteins occurs in specialized endocytic compartments of antigen-presenting cells, which in man are termed MIICs. Newly synthesized MHC class II molecules are transported from the trans-Golgi network to MIICs, but previous studies of this important step in antigen processing have failed to conclusively determine whether most immature MHC class II complexes are transported directly to the processing compartments or are first transiently exposed at the cell surface. To attempt to resolve this question, I constructed a chimeric HLA-DRalpha chain containing two optimal tyrosine sulfation motifs. When expressed in a human B lymphoblastoid cell line lacking functional DRalpha chains, the chimera was correctly incorporated into complexes containing endogenous beta and invariant chains, transported to the trans-Golgi network, and efficiently sulfated. Pulse-chase experiments showed that the sulfated complexes were rapidly transported to processing compartments with kinetics consistent with direct transport from the trans-Golgi network. The rate of maturation was not significantly altered in cells expressing a temperature-sensitive mutant of dynamin under conditions where the endocytosis of transferrin was inhibited by 95%, confirming that endocytosis was not required for delivery to MIICs. Maturation of MHC class II-containing complexes was inhibited by aluminum fluoride and brefeldin A, indicating the involvement of heterotrimeric G-proteins and ADP-ribosylation factor in the transport event(s). The procedure described provides a unique mechanism to study critical events in antigen processing and presentation.     

Functional evidence that conserved TCR CDR alpha 3 loop docking governs the cross-recognition of closely related peptide:class I complexes.

The TCR recognizes its peptide:MHC (pMHC) ligand by assuming a diagonal orientation relative to the MHC helices, but it is unclear whether and to what degree individual TCRs exhibit docking variations when contacting similar pMHC complexes. We analyzed monospecific and cross-reactive recognition by diverse TCRs of an immunodominant HVH-1 glycoprotein B epitope (HSV-8p) bound to two closely related MHC class I molecules, H-2K(b) and H-2K(bm8). Previous studies indicated that the pMHC portion likely to vary in conformation between the two complexes resided at the N-terminal part of the complex, adjacent to peptide residues 2-4 and the neighboring MHC side chains. We found that CTL clones sharing TCR beta-chains exhibited disparate recognition patterns, whereas those with drastically different TCRbeta-chains but sharing identical TCRalpha CDR3 loops displayed identical functional specificity. This suggested that the CDRalpha3 loop determines the TCR specificity in our model, the conclusion supported by modeling of the TCR over the actual HSV-8:K(b) crystal structure. Importantly, these results indicate a remarkable conservation in CDRalpha3 positioning, and, therefore, in docking of diverse TCRalphabeta heterodimers onto variant peptide:class I complexes, implying a high degree of determinism in thymic selection and T cell activation.     

C-terminal anchoring of a peptide to class II MHC via the P10 residue is compatible with a peptide bulge.

The binding of antigenic peptide to class II MHC is mediated by hydrogen bonds between the MHC and the peptide, by salt bridges, and by hydrophobic interactions. The latter are confined to a number of deeper pockets within the peptide binding groove, and peptide side chains that interact with these pockets are referred to as anchor residues. T cell recognition involves solvent-accessible peptide residues along with minor changes in MHC helical pitch induced by the anchor residues. In class I MHC there is an added level of epitope complexity that results from binding of longer peptides that bulge out into the solvent-accessible, T cell contact area. Unlike class I MHC, class II MHC does not bind peptides of discrete length, and the possibility of peptide bulging has not been clearly addressed. A peptide derived from position 24-37 of integrin beta(3) can either bind or not bind to the class II MHC molecule HLA DRB3*0101 based on a polymorphism at the P9 anchor. We show that the loss of binding can be compensated by changes at the P10 position. We propose that this could be an example of a class II peptide bulge. Although not as efficient as P9 anchoring, the use of P10 as an anchor adds another possible mechanism by which T cell epitopes can be generated in the class II presentation system.     

The three-dimensional structure of HLA-B27 at 2.1 A resolution suggests a general mechanism for tight peptide binding to MHC.

Cell surface complexes of class I MHC molecules and bound peptide antigens serve as specific recognition elements controlling the cytotoxic immune response. The 2.1 A structure of the human class I MHC molecule HLA-B27 provides a detailed composite image of a co-crystallized collection of HLA-B27-bound peptides, indicating that they share a common main-chain structure and length. It also permits direct visualization of the conservation of arginine as an "anchor" side chain at the second peptide position, which is bound in a potentially HLA-B27-specific pocket and may therefore have a role in the association of HLA-B27 with several diseases. Tight peptide binding to class I MHC molecules appears to result from the extensive contacts found at the ends of the cleft between peptide main-chain atoms and conserved MHC side chains, which also involve the peptide in stabilizing the three-dimensional fold of HLA-B27. The concentration of binding interactions at the peptide termini permits extensive sequence (and probably some length) variability in the center of the peptide, where it is exposed for T cell recognition.     

Mimetics of a T cell epitope based on poly-N-acylated amine backbone structures induce T cells in vitro and in vivo

Peptidomimetics of the major histocompatibility complex (MHC) class I-restricted ovalbumin-derived T cell epitope SIINFEKL were generated by replacing parts of the peptide backbone by a poly-N-acylated amine (PAA) backbone with aromatic, heteroaromatic, and pseudoaromatic side chains that branch off of the main chain at the amine nitrogen. The structure of the PAAs was designed to position this side chain in the central epitope anchor pocket of the MHC molecule. A number of biologically active PAAs were found that induced cytolysis by the mouse cytotoxic T cell clone 4G3. Competition experiments with independent peptides that are known to bind to the restricting MHC molecule H-2K(b) suggest that the PAAs are bound by the MHC molecules at the same site as conventional peptide epitopes. The PAAs were active also in vivo and induced primary cytotoxic T cell responses in mice.     

Conformational and structural characteristics of peptides binding to HLA-DR molecules.

A fundamental characteristic of MHC class I and class II proteins is their unusual capacity to form stable complexes with a wide spectrum of peptide ligands. In this study, sets of peptide analogues containing long chain-biotinylated lysine individually substituted for each amino acid in the sequence have been used to explore the structural requirements for the formation of peptide-MHC class II protein complexes. Based on the ability of the analogs to bind both the MHC protein and fluorescent streptavidin, receptor contact residues were identified and from their spacing the conformation of the bound peptides could be inferred. Six separate peptides were studied; three defined by HLA-DR1Dw1-restricted T cells, and three identified by T cells restricted through alleles other than HLA-DR1Dw1. The similar patterns of fluorescent signals observed when the former three peptides were studied indicated that they shared conformational features when bound to HLA-DR1Dw1. In contrast when the latter three peptides were examined, the data indicated that they shared some but not all of the conformational features characteristic of the peptides known to elicit HLA-DR1Dw1-restricted T cells. When the peptide sequences were aligned based on the critical contact residues, two positions of structural homology were apparent. In each sequence, an amino acid with a bulky hydrophobic side chain could be identified separated by four residues from a small amino acid. These minimal structural requirements were consistent with recent experiments demonstrating that only a small number of side chains in the peptide were necessary for binding to the MHC protein.     

Spatial structure of isoleucine pentapeptides Glu-Phe-Leu-Arg-Ile-NH2 and Pro-Phe-Tyr-Arg-Ile-NH2

The spatial structure of two cardioactive isoleucine pentapeptides Glu-Phe-Leu-Arg-Ile-NH2 (I) and Pro-Phe-Tyr-Arg-Ile-NH2 (II) have been investigated using the theoretical conformational analysis. The low-energy conformations of these molecules were found, the values of dihedral angles of the backbone and side chains of the amino acid residues constituting these peptides were determined, and the energies of intra- and interresidual interactions were estimated. It was revealed that the spatial structure of molecule I can exist as five and that of molecule II as seven stable backbone forms.     

Ab initio prediction of peptide-MHC binding geometry for diverse class I MHC allotypes

Since determining the crystallographic structure of all peptide-MHC complexes is infeasible, an accurate prediction of the conformation is a critical computational problem. These models can be useful for determining binding energetics, predicting the structures of specific ternary complexes with T-cell receptors, and designing new molecules interacting with these complexes. The main difficulties are (1) adequate sampling of the large number of conformational degrees of freedom for the flexible peptide, (2) predicting subtle changes in the MHC interface geometry upon binding, and (3) building models for numerous MHC allotypes without known structures. Whereas previous studies have approached the sampling problem by dividing the conformational variables into different sets and predicting them separately, we have refined the Biased-Probability Monte Carlo docking protocol in internal coordinates to optimize a physical energy function for all peptide variables simultaneously. We also imitated the induced fit by docking into a more permissive smooth grid representation of the MHC followed by refinement and reranking using an all-atom MHC model. Our method was tested by a comparison of the results of cross-docking 14 peptides into HLA-A*0201 and 9 peptides into H-2K(b) as well as docking peptides into homology models for five different HLA allotypes with a comprehensive set of experimental structures. The surprisingly accurate prediction (0.75 A backbone RMSD) for cross-docking of a highly flexible decapeptide, dissimilar to the original bound peptide, as well as docking predictions using homology models for two allotypes with low average backbone RMSDs of less than 1.0 A illustrate the method's effectiveness. Finally, energy terms calculated using the predicted structures were combined with supervised learning on a large data set to classify peptides as either HLA-A*0201 binders or nonbinders. In contrast with sequence-based prediction methods, this model was also able to predict the binding affinity for peptides to a different MHC allotype (H-2K(b)), not used for training, with comparable prediction accuracy.     

Peptide-binding motifs for the I-Ad MHC class II molecule: alternate pH-dependent binding behavior

The ability of peptides to form stable complexes with MHC class II molecules expressed in the host determines their ability to recruit CD4 T cells during an immune response. In this study, we sought to define the features of the antigenic peptides that control their kinetic stability with I-A(d) because of the diversity of peptides that this molecule is known to present. Peptide dissociation assays indicated that each pocket of I-A(d) displays exquisite sensitivity to side chain structure, size, and charge. Most surprising were results related to the P1 pocket, which has been difficult to define by conventional competition assays. Our studies revealed a considerable degree of specificity in the P1 pocket but also an unexpected degree of structural flexibility. Amino acids with neutral side chains such as Met and the alternatively negatively charged Glu are both highly favored at P1. Interestingly, these two options at the P1 pocket in I-A(d) display dramatically different pH-dependent interactions with the class II molecule. These findings are discussed in the context of a structural model to explain these data and in light of the immunological implications of pH-dependent behavior of class II-peptide complexes in acidic endosomal compartments, where DM-catalyzed loading of class II molecules takes place, and at the neutral pH of the APC cell surface, where class II-peptide complexes promote activation of CD4 T cells.     

Hydrogen bond integrity between MHC class II molecules and bound peptide determines the intracellular fate of MHC class II molecules.

MHC class II molecules associate with peptides through pocket interactions and the formation of hydrogen bonds. The current paradigm suggests that the interaction of side chains of the peptide with pockets in the class II molecule is responsible for the formation of stable class II-peptide complexes. However, recent evidence has shown that the formation of hydrogen bonds between genetically conserved residues of the class II molecule and the main chain of the peptide contributes profoundly to peptide stability. In this study, we have used I-A(k), a class II molecule known to form strong pocket interactions with bound peptides, to probe the general importance of hydrogen bond integrity in peptide acquisition. Our studies have revealed that abolishing hydrogen bonds contributed by positions 81 or 82 in the beta-chain of I-A(k) results in class II molecules that are internally degraded when trafficked through proteolytic endosomal compartments. The presence of high-affinity peptides derived from either endogenous or exogenous sources protects the hydrogen bond-deficient variant from intracellular degradation. Together, these data indicate that disruption of the potential to form a complete hydrogen bond network between MHC class II molecules and bound peptides greatly diminishes the ability of class II molecules to bind peptides. The subsequent failure to stably acquire peptides leads to protease sensitivity of empty class II molecules, and thus to proteolytic degradation before export to the surface of APCs.     

Energetics and cooperativity of the hydrogen bonding and anchor interactions that bind peptides to MHC class II protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).     

MHC Class I/peptide interactions: Binding specificity and kinetics

Recent developments in the preparation of soluble analogues of the major histocompatibility complex (MHC) class l molecules as well as in the applications of real time biosensor technology have permitted the direct analysis of the binding of MHC class l molecules to antigenic peptides. Using synthetic peptide analogues with cysteine substitutions at appropriate positions, peptides can be immobilized on a dextran-modified gold biosensor surface with a specific spatial orientation. A full set of such substituted peptides (known as ‘pepsicles’, as they are peptides on a stick) representing antigenic or self peptides can be used in the functional mapping of the MHC class l peptide binding site. Scans of sets of peptide analogues reveal that some amino acid side chains of the peptide are critical to stable binding to the MHC molecule, while others are not. This is consistent with functional experiments using substituted peptides and three-dimensional molecular models of MHC/peptide complexes. Details analysis of the kinetic dissociation rates (k_d) of the MHC molecules from the specifically coupled solid phase peptides revels that the stability of the complex is a function of the particular peptide, its coupling position, and the MHC molecule. Measured k_d values for antigenic peptide/class I interactions at 25°C are in the range of ca 10^?4–10^?6/s. Biosensor methodology for the analysis of the binding of MHC class I molecules to solid-phase peptides using real time surface plasmon resonance offers a rational approach to the general analysis of protein/peptide interactions.     

Towards universal structure-based prediction of class II MHC epitopes for diverse allotypes

The binding of peptide fragments of antigens to class II MHC proteins is a crucial step in initiating a helper T cell immune response. The discovery of these peptide epitopes is important for understanding the normal immune response and its misregulation in autoimmunity and allergies and also for vaccine design. In spite of their biomedical importance, the high diversity of class II MHC proteins combined with the large number of possible peptide sequences make comprehensive experimental determination of epitopes for all MHC allotypes infeasible. Computational methods can address this need by predicting epitopes for a particular MHC allotype. We present a structure-based method for predicting class II epitopes that combines molecular mechanics docking of a fully flexible peptide into the MHC binding cleft followed by binding affinity prediction using a machine learning classifier trained on interaction energy components calculated from the docking solution. Although the primary advantage of structure-based prediction methods over the commonly employed sequence-based methods is their applicability to essentially any MHC allotype, this has not yet been convincingly demonstrated. In order to test the transferability of the prediction method to different MHC proteins, we trained the scoring method on binding data for DRB1*0101 and used it to make predictions for multiple MHC allotypes with distinct peptide binding specificities including representatives from the other human class II MHC loci, HLA-DP and HLA-DQ, as well as for two murine allotypes. The results showed that the prediction method was able to achieve significant discrimination between epitope and non-epitope peptides for all MHC allotypes examined, based on AUC values in the range 0.632-0.821. We also discuss how accounting for peptide binding in multiple registers to class II MHC largely explains the systematically worse performance of prediction methods for class II MHC compared with those for class I MHC based on quantitative prediction performance estimates for peptide binding to class II MHC in a fixed register.     

Automatic sequence design of major histocompatibility complex class I binding peptides impairing CD8+ T cell recognition

An automatic protein design procedure was used to compute amino acid sequences of peptides likely to bind the HLA-A2 major histocompatibility complex (MHC) class I allele. The only information used by the procedure are a structural template, a rotamer library, and a well established classical empirical force field. The calculations are performed on six different templates from x-ray structures of HLA-A0201-peptide complexes. Each template consists of the bound peptide backbone and the full atomic coordinates of the MHC protein. Sequences within 2 kcal/mol of the minimum energy sequence are computed for each template, and the sequences from all the templates are combined and ranked by their energies. The five lowest energy peptide sequences and five other low energy sequences re-ranked on the basis of their similarity to peptides known to bind the same MHC allele are chemically synthesized and tested for their ability to bind and form stable complexes with the HLA-A2 molecule. The most efficient binders are also tested for inhibition of the T cell receptor recognition of two known CD8(+) T effectors. Results show that all 10 peptides bind the expected MHC protein. The six strongest binders also form stable HLA-A2-peptide complexes, albeit to varying degrees, and three peptides display significant inhibition of CD8(+) T cell recognition. These results are rationalized in light of our knowledge of the three-dimensional structures of the HLA-A2-peptide and HLA-A2-peptide-T cell receptor complexes.     

Changes at the floor of the peptide-binding groove induce a strong preference for proline at position 3 of the bound peptide: molecular dynamics simulations of HLA-A*0217

We report on molecular dynamics simulations of major histocompatibility complex (MHC)-peptide complexes. Class I MHC molecules play an important role in cellular immunity by presenting antigenic peptides to cytotoxic T cells. Pockets in the peptide-binding groove of MHC molecules accommodate anchor side chains of the bound peptide. Amino acid substitutions in MHC affect differences in the peptide-anchor motifs. HLA-A*0217, human MHC class I molecule, differs from HLA-A*0201 only by three amino acid residues substitutions (positions 95, 97, and 99) at the floor of the peptide-binding groove. A*0217 showed a strong preference for Pro at position 3 (p3) and accepted Phe at p9 of its peptide ligands, but these preferences have not been found in other HLA-A2 ligands. To reveal the structural mechanism of these observations, the A*0217-peptide complexes were simulated by 1000 ps molecular dynamics at 300 K with explicit solvent molecules and compared with those of the A*0201-peptide complexes. We examined the distances between the anchor side chain of the bound peptide and the pocket, and the rms fluctuations of the bound peptides and the HLA molecules. On the basis of the results from our simulations, we propose that Pro at p3 serves as an optimum residue to lock the dominant anchor residue (p9) tightly into pocket F and to hold the peptide in the binding groove, rather than a secondary anchor residue fitting optimally the complementary pocket. We also found that Phe at p9 is used to occupy the space created by replacements of three amino acid residues at the floor within the groove. These findings would provide a novel understanding in the peptide-binding motifs of class I MHC molecules.     

Structure-based identification of MHC binding peptides: Benchmarking of prediction accuracy

Identification of MHC binding peptides is essential for understanding the molecular mechanism of immune response. However, most of the prediction methods use motifs/profiles derived from experimental peptide binding data for specific MHC alleles, thus limiting their applicability only to those alleles for which such data is available. In this work we have developed a structure-based method which does not require experimental peptide binding data for training. Our method models MHC-peptide complexes using crystal structures of 170 MHC-peptide complexes and evaluates the binding energies using two well known residue based statistical pair potentials, namely Betancourt-Thirumalai (BT) and Miyazawa-Jernigan (MJ) matrices. Extensive benchmarking of prediction accuracy on a data set of 1654 epitopes from class I and class II alleles available in the SYFPEITHI database indicate that BT pair-potential can predict more than 60% of the known binders in case of 14 MHC alleles with AUC values for ROC curves ranging from 0.6 to 0.9. Similar benchmarking on 29,522 class I and class II MHC binding peptides with known IC(50) values in the IEDB database showed AUC values higher than 0.6 for 10 class I alleles and 9 class II alleles in predictions involving classification of a peptide to be binder or non-binder. Comparison with recently available benchmarking studies indicated that, the prediction accuracy of our method for many of the class I and class II MHC alleles was comparable to the sequence based methods, even if it does not use any experimental data for training. It is also encouraging to note that the ranks of true binding peptides could further be improved, when high scoring peptides obtained from pair potential were re-ranked using all atom forcefield and MM/PBSA method.