Peptide binding by class I and class II MHC molecules

Exploring gene regulatory interaction networks and predicting therapeutic molecules for hypopharyngeal cancer and EGFR-mutated lung adenocarcinoma

Hypopharyngeal cancer gene regulatory networks and therapeutic molecules are a disease that is associated with EGFR-mutated lung adenocarcinoma. Here we utilized a bioinformatics approach to identify genetic commonalities between these two diseases. To this end, we examined microarray datasets from GEO (Gene Expression Omnibus) to identify differentially expressed genes, common genes, and hub genes between the selected two diseases. Our analyses identified potential therapeutic molecules for the selected diseases based on 10 hub genes with the highest interactions according to the degree topology method and the maximum clique centrality (MCC). These therapeutic molecules may have the potential for simultaneous treatment of these diseases.

Evaluation of potential MHC-I allele-specific epitopes in Zika virus proteins and the effects of mutations on peptide-MHC-I interaction studied using in silico approaches

Zika virus (ZIKV) infection is a global health concern due to its association with microcephaly and neurological complications. The development of a T-cell vaccine is important to combat this disease. In this study, we propose ZIKV major histocompatibility complex I (MHC-I) epitopes based on in silico screening consensus followed by molecular docking, PRODIGY, and molecular dynamics (MD) simulation analyses. The effects of the reported mutations on peptide-MHC-I (pMHC-I) complexes were also evaluated. In general, our data indicate an allele-specific peptide-binding human leukocyte antigen (HLA) and potential epitopes. For HLA-B44, we showed that the absence of acidic residue Glu at P2, due to the loss of the electrostatic interaction with Lys45, has a negative impact on the pMHC-I complex stability and explains the low free energy estimated for the immunodominant peptide E-4 (IGVSNRDFV). Our MD data also suggest the deleterious effects of acidic residue Asp at P1 on the pMHC-I stability of HLA-B8 due to destabilization of the α-helix and β-strand. Free energy estimation further indicated that the mutation from Val to Ala at P9 of peptide E-247 (DAHAKRQTV), which was found exclusively in microcephaly samples, did not reduce HLA-B8 affinity. In contrast, the mutation from Thr to Pro at P2 of the peptide NS5-832 (VTKWTDIPY) decreased the interaction energy, number of intermolecular interactions, and adversely affected its binding mode with HLA-A1. Overall, our findings are important with regard to the design of T-cell peptide vaccines and for understanding how ZIKV escapes recognition by CD8 + T-cells.

Major Histocompatibility Complex Class I (FLA-E*01801) Molecular Structure in Domestic Cats Demonstrates Species-Specific Characteristics in Presenting Viral Antigen Peptides

Feline immunodeficiency virus (FIV) infection in domestic cats is the smallest usable natural model for lentiviral infection studies. FLA-E*01801 was applied to FIV AIDS vaccine research. We determined the crystal structure of FLA-E*01801 complexed with a peptide derived from FIV (gag positions 40 to 48; RMANVSTGR [RMA9]). The A pocket of the FLA-E*01801 complex plays a valuable restrictive role in peptide binding. Mutation experiments and circular-dichroism (CD) spectroscopy revealed that peptides with Asp at the first position (P1) could not bind to FLA-E*01801. The crystal structure and in vitro refolding of the mutant FLA-E*01801 complex demonstrated that Glu⁶³ and Trp¹⁶⁷ in the A pocket play important roles in restricting P1D. The B pocket of the FLA-E*01801 complex accommodates M/T/A/V/I/L/S residues, whereas the negatively charged F pocket prefers R/K residues. Based on the peptide binding motif, 125 FLA-E*01801-restricted FIV nonapeptides (San Diego isolate) were identified. Our results provide the structural basis for peptide presentation by the FLA-E*01801 molecule, especially A pocket restriction on peptide binding, and identify the potential cytotoxic T lymphocyte (CTL) epitope peptides of FIV presented by FLA-E*01801. These results will benefit both the reasonable design of FLA-E*01801-restricted CTL epitopes and the further development of the AIDS vaccine.IMPORTANCE Feline immunodeficiency virus (FIV) is a viral pathogen in cats, and this infection is the smallest usable natural model for lentivirus infection studies. To examine how FLA I presents FIV epitope peptides, we crystallized and solved the first classic feline major histocompatibility complex class I (MHC-I) molecular structure. Surprisingly, pocket A restricts peptide binding. Trp¹⁶⁷ blocks the left side of pocket A, causing P1D to conflict with Glu⁶³ We also identified the FLA-E*01801 binding motif X (except D)-(M/T/A/V/I/L/S)-X-X-X-X-X-X-(R/K) based on structural and biochemical experiments. We identified 125 FLA-E*01801-restricted nonapeptides from FIV. These results are valuable for developing peptide-based FIV and human immunodeficiency virus (HIV) vaccines and for studying how MHC-I molecules present peptides.

An N-end rule pathway that recognizes proline and destroys gluconeogenic enzymes

Cells synthesize glucose if deprived of it, and destroy gluconeogenic enzymes upon return to glucose-replete conditions. We found that the Gid4 subunit of the ubiquitin ligase GID in the yeast Saccharomyces cerevisiae targeted the gluconeogenic enzymes Fbp1, Icl1, and Mdh2 for degradation. Gid4 recognized the N-terminal proline (Pro) residue and the ~5-residue-long adjacent sequence motifs. Pck1, the fourth gluconeogenic enzyme, contains Pro at position 2; Gid4 directly or indirectly recognized Pro at position 2 of Pck1, contributing to its targeting. These and related results identified Gid4 as the recognition component of the GID-based proteolytic system termed the Pro/N-end rule pathway. Substrates of this pathway include gluconeogenic enzymes that bear either the N-terminal Pro residue or a Pro at position 2, together with adjacent sequence motifs.

Dipeptides promote folding and peptide binding of MHC class I molecules

MHC class I molecules bind only those peptides with high affinity that conform to stringent length and sequence requirements. We have now investigated which peptides can aid the in vitro folding of class I molecules, and we find that the dipeptide glycyl-leucine efficiently supports the folding of HLA-A*02:01 and H-2K(b) into a peptide-receptive conformation that rapidly binds high-affinity peptides. Treatment of cells with glycyl-leucine induces accumulation of peptide-receptive H-2K(b) and HLA-A*02:01 at the surface of cells. Other dipeptides with a hydrophobic second amino acid show similar enhancement effects. Our data suggest that the dipeptides bind into the F pocket like the C-terminal amino acids of a high-affinity peptide.

Retroviral infections in sheep and goats: small ruminant lentiviruses and host interaction

Small ruminant lentiviruses (SRLV) are members of the Retrovirus family comprising the closely related Visna/Maedi Virus (VMV) and the Caprine Arthritis-Encephalitis Virus (CAEV), which infect sheep and goats. Both infect cells of the monocyte/macrophage lineage and cause lifelong infections. Infection by VMV and CAEV can lead to Visna/Maedi (VM) and Caprine Arthritis-Encephalitis (CAE) respectively, slow progressive inflammatory diseases primarily affecting the lungs, nervous system, joints and mammary glands. VM and CAE are distributed worldwide and develop over a period of months or years, always leading to the death of the host, with the consequent economic and welfare implications. Currently, the control of VM and CAE relies on the control of transmission and culling of infected animals. However, there is evidence that host genetics play an important role in determining Susceptibility/Resistance to SRLV infection and disease progression, but little work has been performed in small ruminants. More research is necessary to understand the host-SRLV interaction.

Molecular cloning and characterization of sea bass (Dicentrarchus labrax, L.) MHC class I heavy chain and β2-microglobulin

In this work, the gene and cDNA of sea bass (Dicentrarchus labrax) β2-microglobulin (Dila-β2m) and several cDNAs of MHC class I heavy chain (Dila-UA) were characterized. While Dila-β2m is single-copy, numerous Dila-UA transcripts were identified per individual with variability at the peptide-binding domain (PBD), but also with unexpected diversity from the connective peptide (CP) through the 3' untranslated region (UTR). Phylogenetic analysis segregates Dila-β2m and Dila-UA into each subfamily cluster, placing them in the fish class and branching Dila-MHC-I with lineage U. The α1 domains resemble those of the recently proposed L1 trans-species lineage. Although no Dila-specific α1, α2 or α3 sub-lineages could be observed, two highly distinct sub-lineages were identified at the CP/TM/CYT regions. The three-dimensional homology model of sea bass MHC-I complex is consistent with other characterized vertebrate structures. Furthermore, basal tissue-specific expression profiles were determined for both molecules, and expression of β2m was evaluated after poly I:C stimulus. Results suggest these molecules are orthologues of other β2m and teleost classical MHC-I and their basic structure is evolutionarily conserved, providing relevant information for further studies on antigen presentation in this fish species.

TCR-MHC docking orientation: natural selection, or thymic selection?

T cell receptors (TCR) dock on their peptide-major histocompatibility complex (pMHC) targets in a conserved orientation. Since amino acid sidechains are the foundation of specific protein-protein interactions, a simple explanation for the conserved docking orientation is that key amino acids encoded by the TCR and MHC genes have been selected and maintained through evolution in order to preserve TCR/pMHC binding. Expectations that follow from the hypothesis that TCR and MHC evolved to interact are discussed in light of the data that both support and refute them. Finally, an alternative and equally simple explanation for the driving force behind the conserved docking orientation is described.

What guides MHC-restricted TCR recognition?

MHC-encoded molecules govern adaptive immune responses by presenting peptides to T cell receptors (TCRs). Based on TCR-MHC crystal structures, we revisit the extent of TCR binding degeneracy, a property with important biological consequences because the diversity of TCR ligands that can be encountered exceeds the number of T cell clones present in a person at any one time. We also discuss whether the approximate diagonal binding of TCR on MHC molecules is due to an intrinsic property of the TCR variable regions, or results from the action of the CD4 and CD8 coreceptors during intrathymic T cell selection. Finally, we discuss how MHC restriction of antigen recognition might have emerged during evolution.

HLA-peptide binding prediction using structural and modeling principles

Short peptides binding to specific human leukocyte antigen (HLA) alleles elicit immune response. These candidate peptides have potential utility in peptide vaccine design and development. The binding of peptides to allele-specific HLA molecule is estimated using competitive binding assay and biochemical binding constants. Application of this method for proteome-wide screening in parasites, viruses, and virulent bacterial strains is laborious and expensive. However, short listing of candidate peptides using prediction approaches have been realized lately. Prediction of peptide binding to HLA alleles using structural and modeling principles has gained momentum in recent years. Here, we discuss the current status of such prediction.

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.

Modeling the Peptide−T Cell Receptor Interaction by the Comparative Molecular Similarity Indices Analysis−Soft Independent Modeling of Class Analogy Technique

A set of 38 epitopes and 183 non-epitopes, which bind to alleles of the HLA-A3 supertype, was subjected to a combination of comparative molecular similarity indices analysis (CoMSIA) and soft independent modeling of class analogy (SIMCA). During the process of T cell recognition, T cell receptors (TCR) interact with the central section of the bound nonamer peptide; thus only positions 4-8 were considered in the study. The derived model distinguished 82% of the epitopes and 73% of the non-epitopes after cross-validation in five groups. The overall preference from the model is for polar amino acids with high electron density and the ability to form hydrogen bonds. These so-called "aggressive" amino acids are flanked by small-sized residues, which enable such residues to protrude from the binding cleft and take an active role in TCR-mediated T cell recognition. Combinations of "aggressive" and "passive" amino acids in the middle part of epitopes constitute a putative TCR binding motif.

Anchor profiles of HLA-specific peptides: analysis by a novel affinity scoring method and experimental validation

The study of intermolecular interactions is a fundamental research subject in biology. Here we report on the development of a quantitative structure-based affinity scoring method for peptide-protein complexes, named PepScope. The method operates on the basis of a highly specific force field function (CHARMM) that is applied to all-atom structural representations of peptide-receptor complexes. Peptide side-chain contributions to total affinity are scored after detailed rotameric sampling followed by controlled energy refinement. A de novo approach to estimate dehydration energies was developed, based on the simulation of individual amino acids in a solvent box filled with explicit water molecules. Transferability of the method was demonstrated by its application to the hydrophobic HLA-A2 and -A24 receptors, the polar HLA-A1, and the sterically ruled HLA-B7 receptor. A combined theoretical and experimental study on 39 anchor substitutions in FxSKQYMTx/HLA-A2 and -A24 complexes indicated a prediction accuracy of about two thirds of a log-unit in Kd. Analysis of free energy contributions identified a great role of desolvation and conformational strain effects in establishing a given specificity profile. Interestingly, the method rightly predicted that most anchor profiles are less specific than so far assumed. This suggests that many potential T-cell epitopes could be missed with current prediction methods. The results presented in this work may therefore significantly affect T-cell epitope discovery programs applied in the field of peptide vaccine development.

The H-2KkMHC Peptide-Binding Groove Anchors the Backbone of an Octameric Antigenic Peptide in an Unprecedented Mode

A wealth of data has accumulated on the structure of mouse MHC class I (MHCI) molecules encoded by the H-2(b) and H-2(d) haplotypes. In contrast, there is a dearth of structural data regarding H-2(k)-encoded molecules. Therefore, the structures of H-2K(k) complexed to an octameric peptide from influenza A virus (HA(259-266)) and to a nonameric peptide from SV40 (SV40(560-568)) have been determined by x-ray crystallography at 2.5 and 3.0 A resolutions, respectively. The structure of the H-2K(k)-HA(259-266) complex reveals that residues located on the floor of the peptide-binding groove contact directly the backbone of the octameric peptide and force it to lie deep within the H-2K(k) groove. This unprecedented mode of peptide binding occurs despite the presence of bulky residues in the middle of the floor of the H-2K(k) peptide-binding groove. As a result, the Calpha atoms of peptide residues P5 and P6 are more buried than the corresponding residues of H-2K(b)-bound octapeptides, making them even less accessible to TCR contact. When bound to H-2K(k), the backbone of the SV40(560-568) nonapeptide bulges out of the peptide-binding groove and adopts a conformation reminiscent of that observed for peptides bound to H-2L(d). This structural convergence occurs despite the totally different architectures of the H-2L(d) and H-2K(k) peptide-binding grooves. Therefore, these two H-2K(k)-peptide complexes provide insights into the mechanisms through which MHC polymorphism outside primary peptide pockets influences the conformation of the bound peptides and have implications for TCR recognition and vaccine design.

Structure of a pheromone receptor-associated MHC molecule with an open and empty groove

Neurons in the murine vomeronasal organ (VNO) express a family of class Ib major histocompatibility complex (MHC) proteins (M10s) that interact with the V2R class of VNO receptors. This interaction may play a direct role in the detection of pheromonal cues that initiate reproductive and territorial behaviors. The crystal structure of M10.5, an M10 family member, is similar to that of classical MHC molecules. However, the M10.5 counterpart of the MHC peptide-binding groove is open and unoccupied, revealing the first structure of an empty class I MHC molecule. Similar to empty MHC molecules, but unlike peptide-filled MHC proteins and non-peptide–binding MHC homologs, M10.5 is thermally unstable, suggesting that its groove is normally occupied. However, M10.5 does not bind endogenous peptides when expressed in mammalian cells or when offered a mixture of class I–binding peptides. The F pocket side of the M10.5 groove is open, suggesting that ligands larger than 8–10-mer class I–binding peptides could fit by extending out of the groove. Moreover, variable residues point up from the groove helices, rather than toward the groove as in classical MHC structures. These data suggest that M10s are unlikely to provide specific recognition of class I MHC–binding peptides, but are consistent with binding to other ligands, including proteins such as the V2Rs.

MHC-like protein M10.5 is expressed in the vomeronasal organ. The structure does not bind endogenous class I-binding peptides, but is thought to interact with a larger V2R pheromone receptor.

Determination of Structural Principles Underlying Three Different Modes of Lymphocytic Choriomeningitis Virus Escape from CTL Recognition

Lymphocytic choriomeningitis virus infection of H-2(b) mice generates a strong CD8(+) CTL response mainly directed toward three immunodominant epitopes, one of which, gp33, is presented by both H-2D(b) and H-2K(b) MHC class I molecules. This CTL response acts as a selective agent for the emergence of viral escape variants. These variants generate altered peptide ligands (APLs) that, when presented by class I MHC molecules, antagonize CTL recognition and ultimately allow the virus to evade the cellular immune response. The emergence of APLs of the gp33 epitope is particularly advantageous for LCMV, as it allows viral escape in the context of both H-2D(b) and H-2K(b) MHC class I molecules. We have determined crystal structures of three different APLs of gp33 in complex with both H-2D(b) and H-2K(b). Comparison between these APL/MHC structures and those of the index gp33 peptide/MHC reveals the structural basis for three different strategies used by LCMV viral escape mutations: 1) conformational changes in peptide and MHC residues that are potential TCR contacts, 2) impairment of APL binding to the MHC peptide binding cleft, and 3) introduction of subtle changes at the TCR/pMHC interface, such as the removal of a single hydroxyl group.

Predicting peptide binding to MHC pockets via molecular modeling, implicit solvation, and global optimization

Development of a computational prediction method based on molecular modeling, global optimization, and implicit solvation has produced accurate structure and relative binding affinity predictions for peptide amino acids binding to five pockets of the MHC molecule HLA-DRB1*0101. Because peptide binding to MHC molecules is essential to many immune responses, development of such a method for understanding and predicting the forces that drive binding is crucial for pharmaceutical design and disease treatment. Underlying the development of this prediction method are two hypotheses. The first is that pockets formed by the peptide binding groove of MHC molecules are independent, separating the prediction of peptide amino acids that bind within individual pockets from those that bind between pockets. The second hypothesis is that the native state of a system composed of an amino acid bound to a protein pocket corresponds to the system's lowest free energy. The prediction method developed from these hypotheses uses atomistic-level modeling, deterministic global optimization, and three methods of implicit solvation: solvent-accessible area, solvent-accessible volume, and Poisson-Boltzmann electrostatics. The method predicts relative binding affinities of peptide amino acids for pockets of HLA-DRB1*0101 by determining computationally an amino acid's global minimum energy conformation. Prediction results from the method are in agreement with X-ray crystallography data and experimental binding assays.

Dynamic flexibility of a peptide-binding groove of human HLA-DR1 class II MHC molecules: normal mode analysis of the antigen peptide-class II MHC complex

Class II major histocompatibility complex (MHC) has tolerance for binding longer antigen peptides than those bound by class I MHC. In this paper, a normal mode analysis on HLA-DR1 class II MHC involving an antigen peptide indicated that the peptide-binding groove had some different dynamic characteristics from that of HLA-A2 class I MHC. The dynamic changes in the class I groove with removal of the bound peptide were limited primarily to the central region and the C-terminal side (corresponding to the C-terminal side of the bound peptide) of the groove, while the dynamic changes in the class II groove with removal of the bound peptide extended to the whole of the groove, and were especially remarkable around a strand located in the N-terminal side (corresponding to the N-terminal side of the bound peptide) of the groove. These results suggest that the N-terminal side of the class II groove is more flexible than the same side of the class I groove, and this flexibility may allow some N-terminal residues of the bound peptide to extend outside the class II groove. Definite anti-correlative motions with removal of the bound peptide appeared between two alpha-helical regions of class II MHC as in the case of class I MHC. These motions of the class II groove may play an important role in obtaining "a flexible dynamic fit" against diverse longer peptides both of whose terminals extend outside the groove.

Associations between measles antibody levels and the protein structure of class II human leukocyte antigens

In this study, we describe associations between variation in human leukocyte antigen (HLA) DP and DQ amino acid sequences and low measles antibody levels after measles immunization. We tested serum samples from 242 children for measles immunoglobulin G antibodies. We performed class II HLA typing and examined associations between DQ and DP exon 2 amino acid sequences and antibody levels. No DP amino acid variants were associated with seronegativity. However, 11 DQA and 6 DQB amino acid variants were associated with seronegativity (p<0.005). These amino acid variants were highly correlated, and the significant DQA amino acids were only found in alleles *0201, *0301, *0401, *0501, and *0601. Two of the amino acids associated with measles seronegativity were located in predicted binding pockets of the DQ molecule; one was present in the leader sequence. Among the DQB alleles, all of the amino acid variants associated with seronegativity were present only in the DQB*0201 allele. Two of the amino acids associated with seronegativity were located in predicted binding pockets of the DQ molecule; two were located in the leader sequence. Our data suggest that specific DQA and DQB amino acid variations are associated with measles seronegativity after vaccination.

Amino acid dimorphism and parasite immune evasion: cellular immune responses to a promiscuous epitope of Plasmodium falciparum merozoite surface protein 1 displaying dimorphic amino acid polymorphism are highly constrained

Like most other surface-exposed antigens of Plasmodium falciparum, the leading malaria vaccine candidate merozoite surface protein (MSP)-1 contains a large number of dimorphic amino acid positions. This type of diversity is presumed to be associated with parasite immune evasion and represents one major obstacle to malaria subunit vaccine development. To understand the precise role of antigen dimorphism in immune evasion, we have analyzed the flexibility of CD4 T cell immune responses against a semi-conserved sequence stretch of the N-terminal block of MSP-1. While this sequence contains overlapping promiscuous T cell epitopes and is a target for growth inhibitory antibodies, three dimorphic amino acid positions may limit its suitability as component of a multi-epitope malaria vaccine. We have analyzed the CD4 T cell responses in a group of human volunteers immunized with a synthetic malaria peptide vaccine containing a single MSP-143-53 sequence variant. All human T cell lines and HLA-DR- or -DP-restricted T cell clones studied were exclusively specific for the sequence variant used for immunization. Competition peptide binding assays with affinity-purified HLA-DR molecules indicated that dimorphism does not primarily affect HLA binding. Modeling studies of the dominant restricting HLA-DRB1*0801 molecule showed that the dimorphic amino acids represent potential TCR contact residues. Lack of productive triggering of the TCR by MHC/variant peptide ligand complexes thus seems to be the characteristic feature of parasite immune evasion associated with antigen dimorphism.

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.

Dynamic characteristics of a peptide-binding groove of human HLA-A2 class I MHC molecules: normal mode analysis of the antigen peptide-class I MHC complex

Class I major histocompatibility complex (MHC) binds antigen peptides with various sequences. We performed a normal mode analysis of HLA-A2 MHC that binds three peptides with different affinity. HLA-A2 MHC has a peptide-binding groove composed of two alpha-helices (residue 49-84, residue 140-179). Some residues in the center of the groove showed an increase in fluctuations and some residue pairs between two helix groups showed a negative change in correlations by removing the antigen peptide. The extent of the fluctuation and correlation changes correlated well with the experimental ranking of the three peptides in binding affinity. Some definite anti-correlative motions were found between two helix groups in low frequency modes (<50 cm(-1)) by removing the antigen peptide. We propose that the above anti-correlative motions play an important role to bind the antigen peptide, especially in obtaining a "dynamic fit."

Types of inter-atomic interactions at the MHC-peptide interface: identifying commonality from accumulated data

BACKGROUND

Quantitative information on the types of inter-atomic interactions at the MHC-peptide interface will provide insights to backbone/sidechain atom preference during binding. Qualitative descriptions of such interactions in each complex have been documented by protein crystallographers. However, no comprehensive report is available to account for the common types of inter-atomic interactions in a set of MHC-peptide complexes characterized by variation in MHC allele and peptide sequence. The available x-ray crystallography data for these complexes in the Protein Databank (PDB) provides an opportunity to identify the prevalent types of such interactions at the binding interface.

RESULTS

We calculated the percentage distributions of four types of interactions at varying inter-atomic distances. The mean percentage distribution for these interactions and their standard deviation about the mean distribution is presented. The prevalence of SS and SB interactions at the MHC-peptide interface is shown in this study. SB is clearly dominant at an inter-atomic distance of 3A.

CONCLUSION

The prevalently dominant SB interactions at the interface suggest the importance of peptide backbone conformation during MHC-peptide binding. Currently, available algorithms are developed for protein sidechain prediction upon fixed backbone template. This study shows the preference of backbone atoms in MHC-peptide binding and hence emphasizes the need for accurate peptide backbone prediction in quantitative MHC-peptide binding calculations.

Recovery of known T-cell epitopes by computational scanning of a viral genome

A new computational method (EpiDock) is proposed for predicting peptide binding to class I MHC proteins, from the amino acid sequence of any protein of immunological interest. Starting from the primary structure of the target protein, individual three-dimensional structures of all possible MHC-peptide (8-, 9- and 10-mers) complexes are obtained by homology modelling. A free energy scoring function (Fresno) is then used to predict the absolute binding free energy of all possible peptides to the class I MHC restriction protein. Assuming that immunodominant epitopes are usually found among the top MHC binders, the method can thus be applied to predict the location of immunogenic peptides on the sequence of the protein target. When applied to the prediction of HLA-A*0201-restricted T-cell epitopes from the Hepatitis B virus, EpiDock was able to recover 92% of known high affinity binders and 80% of known epitopes within a filtered subset of all possible nonapeptides corresponding to about one tenth of the full theoretical list. The proposed method is fully automated and fast enough to scan a viral genome in less than an hour on a parallel computing architecture. As it requires very few starting experimental data, EpiDock can be used: (i) to predict potential T-cell epitopes from viral genomes (ii) to roughly predict still unknown peptide binding motifs for novel class I MHC alleles.

T cell activity correlates with oligomeric peptide-major histocompatibility complex binding on T cell surface

Recognition of virally infected cells by CD8+ T cells requires differentiation between self and nonself peptide-class I major histocompatibility complexes (pMHC). Recognition of foreign pMHC by host T cells is a major factor in the rejection of transplanted organs from the same species (allotransplant) or different species (xenotransplant). AHIII12.2 is a murine T cell clone that recognizes the xenogeneic (human) class I MHC HLA-A2.1 molecule (A2) and the syngeneic murine class I MHC H-2 D(b) molecule (D(b)). Recognition of both A2 and D(b) are peptide-dependent, and the sequences of the peptides recognized have been determined. Alterations in the antigenic peptides bound to A2 cause large changes in AHIII12.2 T cell responsiveness. Crystal structures of three representative peptides (agonist, null, and antagonist) bound to A2 partially explain the changes in AHIII12.2 responsiveness. Using class I pMHC octamers, a strong correlation is seen between T cell activity and the affinity of pMHC complexes for the T cell receptor. However, contrary to previous studies, we see similar half-lives for the pMHC multimers bound to the AHIII12.2 cell surface.

Zooming in on the hydrophobic ridge of H-2D(b): implications for the conformational variability of bound peptides

Class I major histocompatibility complex (MHC) molecules, which display intracellularly processed peptides on the cell surface for scanning by T-cell receptors (TCRs), are extraordinarily polymorphic. MHC polymorphism is believed to result from natural selection, since individuals heterozygous at the corresponding loci can cope with a larger number of pathogens. Here, we present the crystal structures of the murine MHC molecule H-2D(b) in complex with the peptides gp276 and np396 from the lymphocytic choriomeningitis virus (LCMV), solved at 2.18 A and 2.20 A resolution, respectively. The most prominent feature of H-2D(b) is a hydrophobic ridge that cuts across its antigen-binding site, which is conserved in the L(d)-like family of class I MHC molecules. The comparison with previously solved crystal structures of peptide/H-2D(b) complexes shows that the hydrophobic ridge focuses the conformational variability of the bound peptides in a "hot-spot", which could allow optimal TCR interaction and discrimination. This finding suggests a functional reason for the conservation of this structural element.

Class I major histocompatibility complex anchor substitutions alter the conformation of T cell receptor contacts

An immunogenic peptide (GP2) derived from HER-2/neu binds to HLA-A2.1 very poorly. Some altered-peptide ligands (APL) of GP2 have increased binding affinity and generate improved cytotoxic T lymphocyte recognition of GP2-presenting tumor cells, but most do not. Increases in binding affinity of single-substitution APL are not additive in double-substitution APL. A common first assumption about peptide binding to class I major histocompatibility complex is that each residue binds independently. In addition, immunologists interested in immunotherapy frequently assume that anchor substitutions do not affect T cell receptor contact residues. However, the crystal structures of two GP2 APL show that the central residues change position depending on the identity of the anchor residue(s). Thus, it is clear that subtle changes in the identity of anchor residues may have significant effects on the positions of the T cell receptor contact residues.

Towards the MHC-peptide combinatorics

The exponentially increased sequence information on major histocompatibility complex (MHC) alleles points to the existence of a high degree of polymorphism within them. To understand the functional consequences of MHC alleles, 36 nonredundant MHC-peptide complexes in the protein data bank (PDB) were examined. Induced fit molecular recognition patterns such as those in MHC-peptide complexes are governed by numerous rules. The 36 complexes were clustered into 19 subgroups based on allele specificity and peptide length. The subgroups were further analyzed for identifying common features in MHC-peptide binding pattern. The four major observations made during the investigation were: (1) the positional preference of peptide residues defined by percentage burial upon complex formation is shown for all the 19 subgroups and the burial profiles within entries in a given subgroup are found to be similar; (2) in class I specific 8- and 9-mer peptides, the fourth residue is consistently solvent exposed, however this observation is not consistent in class I specific 10-mer peptides; (3) an anchor-shift in positional preference is observed towards the C terminal as the peptide length increases in class II specific peptides; and (4) peptide backbone atoms are proportionately dominant at the MHC-peptide interface.

Specificity in protein-protein interactions: the structural basis for dual recognition in endonuclease colicin-immunity protein complexes

Bacteria producing endonuclease colicins are protected against their cytotoxic activity by virtue of a small immunity protein that binds with high affinity and specificity to inactivate the endonuclease. DNase binding by the immunity protein occurs through a "dual recognition" mechanism in which conserved residues from helix III act as the binding-site anchor, while variable residues from helix II define specificity. We now report the 1.7 A crystal structure of the 24.5 kDa complex formed between the endonuclease domain of colicin E9 and its cognate immunity protein Im9, which provides a molecular rationale for this mechanism. Conserved residues of Im9 form a binding-energy hotspot through a combination of backbone hydrogen bonds to the endonuclease, many via buried solvent molecules, and hydrophobic interactions at the core of the interface, while the specificity-determining residues interact with corresponding specificity side-chains on the enzyme. Comparison between the present structure and that reported recently for the colicin E7 endonuclease domain in complex with Im7 highlights how specificity is achieved by very different interactions in the two complexes, predominantly hydrophobic in nature in the E9-Im9 complex but charged in the E7-Im7 complex. A key feature of both complexes is the contact between a conserved tyrosine residue from the immunity proteins (Im9 Tyr54) with a specificity residue on the endonuclease directing it toward the specificity sites of the immunity protein. Remarkably, this tyrosine residue and its neighbour (Im9 Tyr55) are the pivots of a 19 degrees rigid-body rotation that relates the positions of Im7 and Im9 in the two complexes. This rotation does not affect conserved immunity protein interactions with the endonuclease but results in different regions of the specificity helix being presented to the enzyme.

Modeling the interactions of a peptide-major histocompatibility class I ligand with its receptors. I. Recognition by two alpha beta T cell receptors

A three-dimensional model of the complex between an Influenza Hemagglutinin peptide, Ha255-262, and its restricting element, the mouse major histocompatibility complex (MHC) class I molecule, Kk, was built by homology modeling and subsequently refined by simulated annealing and restrained molecular dynamics. Next, three-dimensional models of two different T cell receptors (TCRs) both specific for the Ha255-262/Kk complex were generated based on previously published TCR X-ray structures. Finally, guided by the recently published X-ray structures of ternary TCR/peptide/MHC-I complexes, the TCR models were successfully docked into the Ha255-262/Kk model. We have previously used a systematic and exhaustive panel of 144 single amino acid substituted analogs to analyze both MHC binding and T cell recognition of the parental viral peptide. This large body of experimental data was used to evaluate the models. They were found to account well for the experimentally obtained data, lending considerable support to the proposed models and suggesting a universal docking mode for alpha beta TCRs to MHC-peptide complexes. Such models may also be useful in guiding future rational experimentation.

The structural basis for the increased immunogenicity of two HIV-reverse transcriptase peptide variant/class I major histocompatibility complexes

Designing altered peptide ligands to generate specific immunological reactivity when bound to class I major histocompatibility complexes is important for both therapeutic and prophylactic reasons. We have previously shown that two altered peptides, derived from human immunodeficiency virus (HIV)-reverse transcriptase (RT) residues 309-317, are more immunogenic in vitro than the wild-type peptide. One peptide variant, I1Y, was able to stimulate RT-specific cytotoxic T cells from the blood of three HIV-infected individuals better than the wild-type RT peptide. Both I1Y and I1F peptide variants increase the cell surface half-life of the peptide-class I complex approximately 3-fold over that of the RT peptide but have different immunological activities. These peptides are candidates for the design of vaccines for HIV due to their increased immunogenicity. To understand the basis for the increased cell surface stability compared with wild-type peptide and to understand the differences in T cell recognition between I1Y and I1F, we determined the x-ray crystal structures of the two class I MHC-peptide complexes. These structures indicate that the increased cell surface half-life is due to pi-pi stacking interactions between Trp-167 of HLA-A2.1 and the aromatic P1 residues of I1F and I1Y. Comparison of the structures and modeling potential T cell receptor (TCR) interactions suggests that T cell interactions and immunogenicity are different between I1Y and I1F for two reasons. First, subtle changes in the steric and polar properties of the I1Y peptide affect TCR engagement. Second, water-mediated hydrogen bond interactions between the P1-Tyr and the P4-Glu peptide residues increase peptide side chain rigidity of residues critical for TCR engagement.

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

Class I major histocompatibility complex (MHC) molecules bind short peptides derived from proteins synthesized within the cell. These complexes of peptide and class I MHC (pMHC) are transported from the endoplasmic reticulum to the cell surface. If a clonotypic T cell receptor expressed on a circulating T cell binds to the pMHC complex, the cell presenting the pMHC is killed. In this manner, some tumor cells expressing aberrant proteins are recognized and removed by the immune system. However, not all tumors are recognized efficiently. One reason hypothesized for poor T cell recognition of tumor-associated peptides is poor binding of those peptides to class I MHC molecules. Many peptides, derived from the proto-oncogene HER-2/neu have been shown to be recognized by cytotoxic T cells derived from HLA-A2(+) patients with breast cancer and other adenocarcinomas. Seven of these peptides were found to bind with intermediate to poor affinity. In particular, GP2 (HER-2/neu residues 654-662) binds very poorly even though it is predicted to bind well based upon the presence of the correct HLA-A2.1 peptide-binding motif. Altering the anchor residues to those most favored by HLA-A2.1 did not significantly improve binding affinity. The crystallographic structure shows that unlike other class I-peptide structures, the center of the peptide does not assume one specific conformation and does not make stabilizing contacts with the peptide-binding cleft.

Quantitative analysis of peptide-MHC class II interaction

The tremendous progress in the field of basic immunology and immunochemistry made in the last decade has significantly advanced our understanding of antigen processing and presentation by MHC class I and II proteins. In this review different techniques to study peptide interaction with MHC class II molecules are summarized and their impact on the elucidation of quantitative parameters, like affinities or kinetic data, is discussed. A recently introduced method based on synthetic combinatorial peptide libraries allows to quantify the binding contribution of each amino acid residue in a class II ligand and is presented in more detail. As this knowledge is fundamental for current investigations in modern medicine, e.g. for novel immune system based therapy concepts, further aspects like the design of new high affinity MHC class II ligands and the prediction of peptide antigens are discussed.

Predicting binding affinities of protein ligands from three-dimensional models: application to peptide binding to class I major histocompatibility proteins

A simple and fast free energy scoring function (Fresno) has been developed to predict the binding free energy of peptides to class I major histocompatibility (MHC) proteins. It differs from existing scoring functions mainly by the explicit treatment of ligand desolvation and of unfavorable protein-ligand contacts. Thus, it may be particularly useful in predicting binding affinities from three-dimensional models of protein-ligand complexes. The Fresno function was independently calibrated for two different training sets: (a) five HLA-A0201-peptide structures, which had been determined by X-ray crystallography, and (b) three-dimensional models of 37 H-2K(k)-peptide structures, which had been obtained by knowledge-based homology modeling. For both training sets, a good cross-validated fit to experimental binding free energies was obtained with predictive errors of 3-3.5 kJ/mol. As expected, lipophilic interactions were found to contribute the most to HLA-A0201-peptide interactions, whereas H-bonding predominates in H-2K(k) recognition. Both cross-validated models were afterward used to predict the binding affinity of a test set of 26 peptides to HLA-A0204 (an HLA allele closely related to HLA-A0201) and of a series of 16 peptides to H-2K(k). Predictions were more accurate for HLA-A2-binding peptides as the training set had been built from experimentally determined structures. The average error in predicting the binding free energy of the test peptides was 3.1 kJ/mol. For the homology model-derived equation, the average error in predicting the binding free energy of peptides to K(k) was significantly higher (5.4 kJ/mol) but still very acceptable. The present scoring function is thus able to predict with a good accuracy binding free energies from three-dimensional models, at the condition that the backbone coordinates of the MHC-bound peptide have first been determined with an accuracy of about 1-1.5 A. Furthermore, it may be easily recalibrated for any protein-ligand complex.

Long-range effects in protein--ligand interactions mediate peptide specificity in the human major histocompatibilty antigen HLA-B27 (B*2701)

B*2701 differs from all other HLA-B27 subtypes of known peptide specificity in that, among its natural peptide ligands, arginine is not the only allowed residue at peptide position 2. Indeed, B*2701 is unique in binding many peptides with Gln2 in vivo. However, the mutation (Asp74Tyr) responsible for altered selectivity is far away from the B pocket of the peptide binding site to which Gln/Arg2 binds. Here, we present a model that explains this effect. It is proposed that a new rotameric state of the conserved Lys70 is responsible for the unique B*2701 binding motif. This side chain should be either kept away from pocket B through its interaction with Asp74 in most HLA-B27 subtypes, or switched to this pocket if residue 74 is Tyr as in B*2701. Involvement of Lys70 in pocket B would thus allow binding of peptides with Gln2. Binding of Arg2-containing peptides to B*2701 is also possible because Lys70 could adopt another conformation, H-bonded to Asn97, which preserves the same binding mode of Arg2 as in B*2705. This model was experimentally validated by mutating Lys70 into Ala in B*2701. Edman sequencing of the B*2701(K70A) peptide pool showed only Arg2, characteristic of HLA-B27-bound peptides, and no evidence for Gln2. This supports the computational model and demonstrates that allowance of B*2701 for peptides with Gln2 is due to the long-range effect of the polymorphic residue 74 of HLA-B27, by inducing a conformational switch of the conserved Lys70.

Nonapeptide analogues containing (R)-3-hydroxybutanoate and beta-homoalanine oligomers: synthesis and binding affinity to a class I major histocompatibility complex protein

Crystal structures of antigenic peptides bound to class I MHC proteins suggest that chemical modifications of the central part of the bound peptide should not alter binding affinity to the MHC restriction protein but could perturb the T-cell response to the parent epitope. In our effort in designing nonpeptidic high-affinity ligands for class I MHC proteins, oligomers of (R)-3-hydroxybutanoate and(or) beta-homoalanine have been substituted for the central part of a HLA-B27-restricted T-cell epitope of viral origin. The affinity of six modified peptides to the B2705 allele was determined by an in vitro stabilization assay. Four out of the six designed analogues presented an affinity similar to that of the parent peptide. Two compounds, sharing the same stereochemistry (R,R,S,S) at the four stereogenic centers of the nonpeptidic spacer, bound to B2705 with a 5-6-fold decreased affinity. Although the chiral spacers do not strongly interact with the protein active site, there are configurations which are not accepted by the MHC binding groove, probably because of improper orientation of some lateral substituents in the bound state and different conformational behavior in the free state. However we demonstrate that beta-amino acids can be incorporated in the sequence of viral T-cell epitopes without impairing MHC binding. The presented structure-activity relationships open the door to the rational design of peptide-based vaccines and of nonnatural T-cell receptor antagonists aimed at blocking peptide-specific T-cell responses in MHC-associated autoimmune diseases.

Structural evidence of T cell xeno-reactivity in the absence of molecular mimicry

The T cell receptor (TCR), from a xeno-reactive murine cytotoxic T lymphocyte clone AHIII12.2, recognizes murine H-2Db complexed with peptide p1027 (FAPGVFPYM), as well as human HLA-A2.1 complexed with peptide p1049 (ALWGFFPVL). A commonly proposed model (the molecular mimicry model) used to explain TCR cross-reactivity suggests that the molecular surfaces of the recognized complexes are similar in shape, charge, or both, in spite of the primary sequence differences. To examine the mechanism of xeno-reactivity of AHIII12.2, we have determined the crystal structures of A2/p1049 and Db/p1027 to 2.5 A and 2.8 A resolution, respectively. The crystal structures show that the TCR footprint regions of the two class I complexes are significantly different in shape and charge. We propose that rather than simple molecular mimicry, unpredictable arrays of common and differential contacts on the two class I complexes are used for their recognition by the same TCR.

From peptides to peptidomimetics: design of nonpeptide ligands for major histocompatibility proteins

The ever increasing data available on antigen presentation by class I or class II histocompatibility proteins have made these glycoproteins highly interesting pharmaceutical targets for either vaccination or immunosuppressive therapy of autoimmune diseases and cancers. Herewith, we review the design and biological properties of the very first nonpeptide ligands of major histocompatibility proteins as well as their potential application in vaccination, Major Histocompatibility Complex (MHC) blockade or T cell receptor antagonism.

Altered peptide ligand design: altering immune responses to class I MHC/peptide complexes

Class I proteins are responsible for binding proteins from endogenously synthesized proteins and displaying them on the cell surface. Our understanding of this process has reached the point where we can manipulate the biochemical properties of peptide/class I binding and determine the effects of this alteration on subsequent immune responses. In this article, we will review the biochemistry of peptide/class I binding, and the effects of structure on this interaction between class I proteins and their peptide ligands. We will review the data which suggest that the major relevant biochemical parameter of class I peptide binding is the off-rate. We will show that the design of altered ligands with improved binding, thermostability and immunogenicity is possible.