Homology modeling and molecular dynamics simulations of MUC1-9/H-2Kb complex suggest novel binding interactions

Human MUC1 is over-expressed in human adenocarcinomas and has been used as a target for immunotherapy studies. The 9-mer MUC1-9 peptide has been identified as one of the peptides which binds to murine MHC class I H-2K^b. The structure of MUC1-9 in complex with H-2K^b has been modeled and simulated with classical molecular dynamics, based on the x-ray structure of the SEV9 peptide/H-2K^b complex. Two independent trajectories with the solvated complex (10 ns in length) were produced. Approximately 12 hydrogen bonds were identified during both trajectories to contribute to peptide/MHC complex, as well as 1-2 water mediated hydrogen bonds. Stability of the complex was also confirmed by buried surface area analysis, although the corresponding values were about 20% lower than those of the original x-ray structure. Interestingly, a bulged conformation of the peptide’s central region, partially characterized as a β-turn, was found exposed form the binding groove. In addition, P1 and P9 residues remained bound in the A and F binding pockets, even though there was a suggestion that P9 was more flexible. The complex lacked numerous water mediated hydrogen bonds that were present in the reference peptide x-ray structure. Moreover, local displacements of residues Asp4, Thr5 and Pro9 resulted in loss of some key interactions with the MHC molecule. This might explain the reduced affinity of the MUC1-9 peptide, relatively to SEV9, for the MHC class I H-2K^b.     

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).     

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.     

Crystal structure of a non-canonical high affinity peptide complexed with MHC class I: a novel use of alternative anchors

The crystal structure of a non-standard peptide, YEA9, in complex with H-2Kb, at 1.5 A resolution demonstrates how YEA9 peptide can bind with surprisingly high affinity through insertion of alternative, long, non-canonical anchors into the B and E pockets. The use of "alternative pockets" represents a new mode of high affinity peptide binding, that should be considered when predicting peptide epitopes for MHC class I. These novel interactions encountered in this non-canonical high affinity peptide-MHC complex should help predict additional binding peptides from primary protein sequences and aid in the design of alternative approaches for peptide-based vaccines.     

Insights into the structure of the LC13 TCR/HLA-B8-EBV peptide complex with molecular dynamics simulations

One key step in the immune response against infected or tumor cells is the recognition of the T-cell receptor (TCR) by class I major histocompatibility complexes. The complex between the HLA-B8 molecule and the immunodominant peptide with sequence FLRGRAYGL, derived from the Epstein-Barr virus, with the LC13 TCR has been determined by X-ray diffraction. The complex has been used as a starting point in a molecular dynamics study in order to investigate the dynamics of the complex association and to explore the specific interactions of the complex formation. The analyzed structures provided evidence that the peptide adopts an open type β-turn conformation close to C-terminal part, which dominates peptide/TCR interactions. Conformational energy landscape analysis indicated the presence of two conformational clusters in the peptide’s structure, underlying the backbone flexibility of the peptide despite being surrounded by two receptors. The peptide/MHC/TCR interface was found to hold significant number of solvent molecules, more specifically the peptide has been found to have approximately seventeen hydrogen bonds with water molecules. The molecular dynamics simulation indicated the disruption of some MHC/TCR contacts, mainly with the CDR1α loop. However, several other interactions emerged that resulted in a stable association during the 20 ns trajectory, as revealed by the buried surface area analysis.     

Cooperativity during the formation of peptide/MHC class II complexes

To generate an effective immune response, class II major histocompatibility complex molecules (MHCII) must present a diverse array of peptide ligands for recognition by T lymphocytes. Peptide/MHCII complexes are stabilized by hydrophobic anchoring of peptide side chains to pockets in the MHCII protein and the formation of hydrogen bonds to the peptide backbone. Many current models of peptide/MHCII association assume an additive and independent contribution of the interactions between major MHCII pockets and corresponding side chains in the peptide. However, significant conformational rearrangements occur in both the peptide and MHCII during binding. Therefore, we hypothesize that peptide binding to MHCII could be viewed as a folding process in which both molecules cooperate to produce the final conformation. To directly test this hypothesis, we adapt a serial mutagenesis strategy to study cooperativity in the interaction of the human MHCII HLA-DR1 and a peptide derived from influenza hemagglutinin. Substitutions in either the peptide or HLA-DR1 that are predicted to interfere with hydrogen bond formation show cooperative effects on complex stability and affinity. Substitution of a peptide side chain that provides a hydrophobic contact also contributes to the cooperative effect, suggesting a role for all energetic sources in the folding process. We propose that cooperativity throughout the peptide-binding groove reflects the folding of segments of the MHCII molecule into helices around the peptide with a concomitant folding of the peptide into a polyproline helix. The implications of cooperativity for peptide/MHCII structure and epitope selection are discussed.     

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.     

Peptidic termini play a significant role in TCR recognition

TCR recognition of class I MHC is dependent on the composition of the antigenic peptide and the MHC. Single amino acid substitutions in either the MHC or the peptide may dramatically alter recognition. While the major interactions between TCR and the peptide/MHC complex appear to be focused on the complementarity-determining region (CDR)3, it is also clear from the cocrystal structure of class I MHC and TCR that the amino and carboxyl ends of the peptide may play a role through interactions with the CDR1. In this work we show that gp33 variants substituted at the peptidic termini at the putative CDR1 contact regions show improved recognition in B6 mice. The rank order of recognition is different using the P14 transgenic T cells, suggesting that one reason for improved recognition is a change in the TCR repertoire that recognizes the peptide. However, the affinity of the TCR by some of the peptide/MHC complex with increased recognition is improved, as shown by increased tetramer binding to P14 T cells. These substitutions at the termini of the peptide-binding cleft cause localized conformational changes as seen by changes in mAb binding and crystallographic structures. The different peptide structures also show different conformations in the center of the peptide, but these are shown to be energetically similar and thus most likely have no significance with respect to TCR recognition. Therefore, small conformational changes, localized to the CDR1 contact regions, may play a significant role in TCR recognition.     

Structure of a human autoimmune TCR bound to a myelin basic protein self-peptide and a multiple sclerosis-associated MHC class II molecule

Multiple sclerosis is mediated by T-cell responses to central nervous system antigens such as myelin basic protein (MBP). To investigate self-peptide/major histocompatibility complex (MHC) recognition and T-cell receptor (TCR) degeneracy, we determined the crystal structure, at 2.8 A resolution, of an autoimmune TCR (3A6) bound to an MBP self-peptide and the multiple sclerosis-associated MHC class II molecule, human leukocyte antigen (HLA)-DR2a. The complex reveals that 3A6 primarily recognizes the N-terminal portion of MBP, in contrast with antimicrobial and alloreactive TCRs, which focus on the peptide center. Moreover, this binding mode, which may be frequent among autoimmune TCRs, is compatible with a wide range of orientation angles of TCR to peptide/MHC. The interface is characterized by a scarcity of hydrogen bonds between TCR and peptide, and TCR-induced conformational changes in MBP/HLA-DR2a, which likely explain the low observed affinity. Degeneracy of 3A6, manifested by recognition of superagonist peptides bearing substitutions at nearly all TCR-contacting positions, results from the few specific interactions between 3A6 and MBP, allowing optimization of interface complementarity through variations in the peptide.     

Main chain hydrogen bond interactions in the binding of proline-rich gluten peptides to the celiac disease-associated HLA-DQ2 molecule

Binding of peptide epitopes to major histocompatibility complex proteins involves multiple hydrogen bond interactions between the peptide main chain and major histocompatibility complex residues. The crystal structure of HLA-DQ2 complexed with the alphaI-gliadin epitope (LQPFPQPELPY) revealed four hydrogen bonds between DQ2 and peptide main chain amides. This is remarkable, given that four of the nine core residues in this peptide are proline residues that cannot engage in amide hydrogen bonding. Preserving main chain hydrogen bond interactions despite the presence of multiple proline residues in gluten peptides is a key element for the HLA-DQ2 association of celiac disease. We have investigated the relative contribution of each main chain hydrogen bond interaction by preparing a series of N-methylated alphaI epitope analogues and measuring their binding affinity and off-rate constants to DQ2. Additionally, we measured the binding of alphaI-gliadin peptide analogues in which norvaline, which contains a backbone amide hydrogen bond donor, was substituted for each proline. Our results demonstrate that hydrogen bonds at P4 and P2 positions are most important for binding, whereas the hydrogen bonds at P9 and P6 make smaller contributions to the overall binding affinity. There is no evidence for a hydrogen bond between DQ2 and the P1 amide nitrogen in peptides without proline at this position. This is a unique feature of DQ2 and is likely a key parameter for preferential binding of proline-rich gluten peptides and development of celiac disease.     

Structure of a covalently stabilized complex of a human alphabeta T-cell receptor, influenza HA peptide and MHC class II molecule, HLA-DR1

An alphabeta T-cell receptor (alphabetaTCR)/hemagglutinin (HA) peptide/human leukocyte antigen (HLA)-DR1 complex was stabilized by flexibly linking the HA peptide with the human HA1.7 alphabetaTCR, to increase the local concentration of the interacting proteins once the peptide has been loaded onto the major histocompatibility complex (MHC) molecule. The structure of the complex, determined by X-ray crystallography, has a binding mode similar to that of the human B7 alphabetaTCR on a pMHCI molecule. Twelve of the 15 MHC residues contacted are at the same positions observed earlier in class I MHC/peptide/TCR complexes. One contact, to an MHC loop outside the peptide-binding site, is conserved and specific to pMHCII complexes. TCR gene usage in the response to HA/HLA-DR appears to conserve charged interactions between three lysines of the peptide and acidic residues on the TCR.     

Crystallographic structure of the human leukocyte antigen DRA, DRB3*0101: models of a directional alloimmune response and autoimmunity

We describe structural studies of the human leukocyte antigen DR52a, HLA-DRA/DRB3*0101, in complex with an N-terminal human platelet integrin alphaII(B)betaIII glycoprotein peptide which contains a Leu/Pro dimorphism. The 33:Leu dimorphism is the epitope for the T cell directed response in neonatal alloimmune thrombocytopenia and post-transfusion purpura in individuals with the alphaII(B)betaIII 33:Pro allele, and defines the unidirectional alloimmune response. This condition is always associated with DR52a. The crystallographic structure has been refined to 2.25 A. There are two alphabeta heterodimers to the asymmetric unit in space group P4(1)2(1)2. The molecule is characterized by two prominent hydrophobic pockets at either end of the peptide binding cleft and a deep, narrower and highly charged P4 opening underneath the beta 1 chain. Further, the peptide in the second molecule displays a sharp upward turn after pocket P9. The structure reveals the role of pockets and the distinctive basic P4 pocket, shared by DR52a and DR3, in selecting their respective binding peptide repertoire. We observe an interesting switch in a residue from the canonically assigned pocket 6 seen in prior class II structures to pocket 4. This occludes the P6 pocket helping to explain the distinctive "1-4-9" peptide binding motif. A beta57 Asp-->Val substitution abrogates the salt-bridge to alpha76 Arg and along with a hydrophobic beta37 is important in shaping the P9 pocket. DRB3*0101 and DRB1*0301 belong to an ancestral haplotype and are associated with many autoimmune diseases linked to antigen presentation, but whereas DR3 is susceptible to type 1 diabetes DR52a is not. This dichotomy is explored for clues to the disease.     

Promiscuous antigen presentation by the nonclassical MHC Ib Qa-2 is enabled by a shallow, hydrophobic groove and self-stabilized peptide conformation.

BACKGROUND: Qa-2 is a nonclassical MHC Ib antigen, which has been implicated in both innate and adaptive immune responses, as well as embryonic development. Qa-2 has an unusual peptide binding specificity in that it requires two dominant C-terminal anchor residues and is capable of associating with a substantially more diverse array of peptide sequences than other nonclassical MHC. RESULTS: We have determined the crystal structure, to 2.3 A, of the Q9 gene of murine Qa-2 complexed with a self-peptide derived from the L19 ribosomal protein, which is abundant in the pool of peptides eluted from the Q9 groove. The 9 amino acid peptide is bound high in a shallow, hydrophobic binding groove of Q9, which is missing a C pocket. The peptide makes few specific contacts and exhibits extremely poor shape complementarity to the MHC groove, which facilitates the presentation of a degenerate array of sequences. The L19 peptide is in a centrally bulged conformation that is stabilized by intramolecular interactions from the invariant P7 histidine anchor residue and by a hydrophobic core of preferred secondary anchor residues that have minimal interaction with the MHC. CONCLUSIONS: Unexpectedly, the preferred secondary peptide residues that exhibit tenuous contact with Q9 contribute significantly to the overall stability of the peptide-MHC complex. The structure of this complex implies a "conformational" selection by Q9 for peptide residues that optimally stabilize the large bulge rather than making intimate contact with the MHC pockets.     

Disruption of Hydrogen Bonds between Major Histocompatibility Complex Class II and the Peptide N-Terminus Is Not Sufficient to Form a Human Leukocyte Antigen-DM Receptive State of Major Histocompatibility Complex Class II

Peptide presentation by MHC class II is of critical importance to the function of CD4+ T cells. HLA-DM resides in the endosomal pathway and edits the peptide repertoire of newly synthesized MHC class II molecules before they are exported to the cell surface. HLA-DM ensures MHC class II molecules bind high affinity peptides by targeting unstable MHC class II:peptide complexes for peptide exchange. Research over the past decade has implicated the peptide N-terminus in modulating the ability of HLA-DM to target a given MHC class II:peptide combination. In particular, attention has been focused on both the hydrogen bonds between MHC class II and peptide, and the occupancy of the P1 anchor pocket. We sought to solve the crystal structure of a HLA-DR1 molecule containing a truncated hemagglutinin peptide missing three N-terminal residues compared to the full-length sequence (residues 306–318) to determine the nature of the MHC class II:peptide species that binds HLA-DM. Here we present structural evidence that HLA-DR1 that is loaded with a peptide truncated to the P1 anchor residue such that it cannot make select hydrogen bonds with the peptide N-terminus, adopts the same conformation as molecules loaded with full-length peptide. HLA-DR1:peptide combinations that were unable to engage up to four key hydrogen bonds were also unable to bind HLA-DM, while those truncated to the P2 residue bound well. These results indicate that the conformational changes in MHC class II molecules that are recognized by HLA-DM occur after disengagement of the P1 anchor residue.     

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

Major histocompatibility complex (MHC) class I molecules (proteins) bind peptides of eight to ten amino acids to present them at the cell surface to cytotoxic T cells. The class I binding groove binds the peptide via hydrogen bonds with the peptide termini and via diverse interactions with the anchor residue side chains of the peptide. To elucidate which of these interactions is most important for the thermodynamic and kinetic stability of the peptide-bound state, we have combined molecular dynamics simulations and experimental approaches in an investigation of the conformational dynamics and binding parameters of a murine class I molecule (H-2K^b) with optimal and truncated natural peptide epitopes. We show that the F pocket region dominates the conformational and thermodynamic properties of the binding groove, and that therefore the binding of the C terminus of the peptide to the F pocket region plays a crucial role in bringing about the peptide-bound state of MHC class I.     

The structure and stability of an HLA-A*0201/octameric tax peptide complex with an empty conserved peptide-N-terminal binding site

The crystal structure of the human class I MHC molecule HLA-A2 complexed with of an octameric peptide, Tax8 (LFGYPVYV), from human T cell lymphotrophic virus-1 (HTLV-1) has been determined. This structure is compared with a newly refined, higher resolution (1.8 A) structure of HLA-A2 complexed with the nonameric Tax9 peptide (LLFGYPVYV) with one more N-terminal residue. Despite the absence of a peptide residue (P1) bound in the conserved N-terminal peptide-binding pocket of the Tax8/HLA-A2 complex, the structures of the two complexes are essentially identical. Water molecules in the Tax8 complex replace the terminal amino group of the Tax9 peptide and mediate a network of hydrogen bonds among the secondary structural elements at that end of the peptide-binding groove. Thermal denaturation measurements indicate that the Tax8 complex is much less stable, DeltaTm = 16 degrees C, than the Tax9 complex, but both can sensitize target cells for lysis by some Tax-specific CTL from HTLV-1 infected individuals. The absence of a P1 peptide residue is thus not enough to prevent formation of a "closed conformation" of the peptide-binding site. TCR affinity measurements and cytotoxic T cell assays indicate that the Tax8/HLA-A2 complex does not functionally cross-react with the A6-TCR-bearing T cell clone specific for Tax9/HLA-A2 complexes.     

Structural definition of the H-2Kd peptide-binding motif

Classic major histocompatibility complex (MHC) proteins associate with antigen- and self-derived peptides in an allele-specific manner. Herein we present the crystal structure of the MHC class I protein H-2K(d) (K(d)) expressed by BALB/c mice in complex with an antigenic peptide derived from influenza A/PR/8/34 nucleoprotein (Flu, residues 147-155, TYQRTRALV). Analysis of our structure in conjunction with the sequences of naturally processed epitopes provides a comprehensive understanding of the dominant K(d) peptide-binding motif. We find that Flu residues Tyr(P2), Thr(P5), and Val(P9) are sequestered into the B, C, and F pockets of the K(d) groove, respectively. The shape and chemistry of the polymorphic B pocket make it an optimal binding site for the side chain of Tyr(P2) as the dominant anchoring residue of nonameric peptides. The non-polar F pocket limits the amino acid repertoire at P9 to hydrophobic residues such as Ile, Leu, or Val, whereas the C pocket restricts the size of the P5-anchoring side chain. We also show that Flu is accommodated in the complex through an unfavorable kink in the otherwise extended peptide backbone due to the presence of a prominent ridge in the K(d) groove. Surprisingly, this backbone conformation is strikingly similar to D(b)-presented peptides despite the fact that these proteins employ distinct motif-anchoring strategies. The results presented in this study provide a solid foundation for the understanding of K(d)-restricted antigen presentation and recognition events.     

MPID: MHC-Peptide Interaction Database for sequence-structure-function information on peptides binding to MHC molecules

SUMMARY: Binding of short antigenic peptides to Major histocompatibility complex (MHC) proteins is the first step in T-cell mediated immune response. To understand the structural principles governing MHC-specific peptide recognition and binding, we have developed the MHC-Peptide Interaction Database (MPID), containing sequence-structure-function information. MPID (version 1.2) contains curated x-ray crystallographic data on 86 MHC peptide complexes, with precomputed interaction parameters (solvent accessibility, hydrogen bonds, gap volume and gap index). A user-friendly web interface and query tools will facilitate the development of predictive algorithms for MHC-peptide binding from a structural viewpoint. AVAILABILITY: Freely accessible from http://surya.bic.nus.edu.sg/mpid.     

Cutting edge: a single, essential hydrogen bond controls the stability of peptide-MHC class II complexes.

The binding of peptides to MHC class II molecules is mediated in part by a conserved array of intermolecular hydrogen bonds. We have evaluated the consequences of disrupting the hydrogen bond between beta-His-81 of the class II molecule and bound peptide. These studies revealed that peptide dissociation rates were accelerated by factors ranging to 200-fold. The sensitivity of a peptide to loss of the hydrogen bond is inversely correlated with the inherent kinetic stability of the peptide-MHC complex. The same relationship has been observed between inherent kinetic stability and the susceptibility to DM. Given that the rate enhancement observed for MHC class II I-Ad protein mutated at position 81 in the beta-chain is comparable with DM-catalyzed rates for other class II molecules, we suggest that DM could function by stabilizing a peptide-MHC intermediate in which one or more hydrogen bonds between the peptide and MHC, such as that contributed by the beta-His-81 hydrogen bond, are disrupted.