Coupling in silico and in vitro analysis of peptide-MHC binding: a bioinformatic approach enabling prediction of superbinding peptides and anchorless epitopes

The ability to define and manipulate the interaction of peptides with MHC molecules has immense immunological utility, with applications in epitope identification, vaccine design, and immunomodulation. However, the methods currently available for prediction of peptide-MHC binding are far from ideal. We recently described the application of a bioinformatic prediction method based on quantitative structure-affinity relationship methods to peptide-MHC binding. In this study we demonstrate the predictivity and utility of this approach. We determined the binding affinities of a set of 90 nonamer peptides for the MHC class I allele HLA-A*0201 using an in-house, FACS-based, MHC stabilization assay, and from these data we derived an additive quantitative structure-affinity relationship model for peptide interaction with the HLA-A*0201 molecule. Using this model we then designed a series of high affinity HLA-A2-binding peptides. Experimental analysis revealed that all these peptides showed high binding affinities to the HLA-A*0201 molecule, significantly higher than the highest previously recorded. In addition, by the use of systematic substitution at principal anchor positions 2 and 9, we showed that high binding peptides are tolerant to a wide range of nonpreferred amino acids. Our results support a model in which the affinity of peptide binding to MHC is determined by the interactions of amino acids at multiple positions with the MHC molecule and may be enhanced by enthalpic cooperativity between these component interactions.     

Exploiting the Peptide — MHC Water Interface in the Computer-Aided Design of Non-Natural Peptides that Bind to the Class I MHC Molecule HLA-A2

Abstract

Class I major histocompatibility complex (MHC) molecules bind peptides derived from intra-cellular proteins and present them to cytotoxic T cells. Certain human immunological diseases are associated with errors in this process. Here we describe an approach to the design of non-natural peptides that could potentially interfere with peptide presentation associated with autoimmune diseases. We have shown previously that the interaction of the peptide GILGFVFTL with the MHC molecule HLA-A2 is mediated by a network of water molecules. In principle, the addition of hydroxyl groups to the peptide could allow for an enhanced interaction of the modified peptide with this water network. Here we illustrate this approach using a peptide having the non-natural amino acid homoserine at position 3, GIhSGFVFTL, and also peptides in which the Cα(F5)—CO—NH1—Cα(V6) peptide bond is replaced by an ether. Cα(F5)—CH(X)—O—Cα(V6), to give the non-natural peptide GILGF—CH(X)—O—VFTL, where X = CH₂OH or CH₃. In a 200 ps solvated molecular dynamics simulation of the HLA-A2 complexes of each peptide for GIhSGFVFTL and GILGF—CH(CH₂OH)—O—VFTL the peptide conformation remained essentially unchanged from that of GILGFVFTL in the X-ray structure of its complex with HLA-A2. In contrast, for GILGF—CH(CH₃)—O—VFTL the peptide conformation deviated from the X-ray conformation, indicating the importance of the hydroxyl group.     

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

Candidate epitope identification using peptide property models: application to cancer immunotherapy

Peptides derived from pathogens or tumors are selectively presented by the major histocompatibility complex proteins (MHC) to the T lymphocytes. Antigenic peptide-MHC complexes on the cell surface are specifically recognized by T cells and, in conjunction with co-factor interactions, can activate the T cells to initiate the necessary immune response against the target cells. Peptides that are capable of binding to multiple MHC molecules are potential T cell epitopes for diverse human populations that may be useful in vaccine design. Bioinformatical approaches to predict MHC binding peptides can facilitate the resource-consuming effort of T cell epitope identification. We describe a new method for predicting MHC binding based on peptide property models constructed using biophysical parameters of the constituent amino acids and a training set of known binders. The models can be applied to development of anti-tumor vaccines by scanning proteins over-expressed in cancer cells for peptides that bind to a variety of MHC molecules. The complete algorithm is described and illustrated in the context of identifying candidate T cell epitopes for melanomas and breast cancers. We analyzed MART-1, S-100, MBP, and CD63 for melanoma and p53, MUC1, cyclin B1, HER-2/neu, and CEA for breast cancer. In general, proteins over-expressed in cancer cells may be identified using DNA microarray expression profiling. Comparisons of model predictions with available experimental data were assessed. The candidate epitopes identified by such a computational approach must be evaluated experimentally but the approach can provide an efficient and focused strategy for anti-cancer immunotherapy development.     

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.     

A single autoimmune T cell receptor recognizes more than a million different peptides

The T cell receptor (TCR) orchestrates immune responses by binding to foreign peptides presented at the cell surface in the context of major histocompatibility complex (MHC) molecules. Effective immunity requires that all possible foreign peptide-MHC molecules are recognized or risks leaving holes in immune coverage that pathogens could quickly evolve to exploit. It is unclear how a limited pool of <10(8) human TCRs can successfully provide immunity to the vast array of possible different peptides that could be produced from 20 proteogenic amino acids and presented by self-MHC molecules (>10(15) distinct peptide-MHCs). One possibility is that T cell immunity incorporates an extremely high level of receptor degeneracy, enabling each TCR to recognize multiple peptides. However, the extent of such TCR degeneracy has never been fully quantified. Here, we perform a comprehensive experimental and mathematical analysis to reveal that a single patient-derived autoimmune CD8(+) T cell clone of pathogenic relevance in human type I diabetes recognizes >one million distinct decamer peptides in the context of a single MHC class I molecule. A large number of peptides that acted as substantially better agonists than the wild-type "index" preproinsulin-derived peptide (ALWGPDPAAA) were identified. The RQFGPDFPTI peptide (sampled from >10(8) peptides) was >100-fold more potent than the index peptide despite differing from this sequence at 7 of 10 positions. Quantification of this previously unappreciated high level of CD8(+) T cell cross-reactivity represents an important step toward understanding the system requirements for adaptive immunity and highlights the enormous potential of TCR degeneracy to be the causative factor in autoimmune disease.     

Chicken major histocompatibility complex class II molecules of the B haplotype present self and foreign peptides

The chicken major histocompatibility complex (MHC), or B-complex, mediates genetic resistance and susceptibility to infectious disease. For example, the B19 haplotype is associated with susceptibility to Marek's disease. Here, we describe the sequencing and analysis of peptides presented by B19 MHC class II molecules. A B19/B19 B-cell line was used for the immunoaffinity purification of MHC class II molecules, which was followed by acid elution of the bound peptides. The eluted peptides were then analysed using tandem mass spectrometry. Thirty peptide sequences were obtained, ranging from 11 to 25 amino acids in length. Source protein cellular localization included the plasma membrane, cytosol and endosomal pathway. In addition, five peptides from the envelope glycoprotein of chicken syncytial virus (CSV) were identified. Chicken syncytial virus had been used as a helper virus along with reticuloendotheliosis virus strain T for transformation of B19/B19B cells. Alignment and analysis of the peptide sequence pool provided a putative peptide-binding motif for the B19 MHC class II.     

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.     

The dominantly expressed class II molecule from a resistant MHC haplotype presents only a few Marek’s disease virus peptides by using an unprecedented binding motif

Viral diseases pose major threats to humans and other animals, including the billions of chickens that are an important food source as well as a public health concern due to zoonotic pathogens. Unlike humans and other typical mammals, the major histocompatibility complex (MHC) of chickens can confer decisive resistance or susceptibility to many viral diseases. An iconic example is Marek’s disease, caused by an oncogenic herpesvirus with over 100 genes. Classical MHC class I and class II molecules present antigenic peptides to T lymphocytes, and it has been hard to understand how such MHC molecules could be involved in susceptibility to Marek’s disease, given the potential number of peptides from over 100 genes. We used a new in vitro infection system and immunopeptidomics to determine peptide motifs for the 2 class II molecules expressed by the MHC haplotype B2, which is known to confer resistance to Marek’s disease. Surprisingly, we found that the vast majority of viral peptide epitopes presented by chicken class II molecules arise from only 4 viral genes, nearly all having the peptide motif for BL2*02, the dominantly expressed class II molecule in chickens. We expressed BL2*02 linked to several Marek’s disease virus (MDV) peptides and determined one X-ray crystal structure, showing how a single small amino acid in the binding site causes a crinkle in the peptide, leading to a core binding peptide of 10 amino acids, compared to the 9 amino acids in all other reported class II molecules. The limited number of potential T cell epitopes from such a complex virus can explain the differential MHC-determined resistance to MDV, but raises questions of mechanism and opportunities for vaccine targets in this important food species, as well as providing a basis for understanding class II molecules in other species including humans.

This study shows that chicken MHC class II molecules present peptides from only a handful of the more than 100 genes of the oncogenic herpesvirus Marek’s disease virus, explaining the strong genetic association of chicken MHC with resistance and susceptibility to this and other economically-important pathogens.     

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.     

One-pot, mix-and-read peptide-MHC tetramers

Background

Cytotoxic T Lymphocytes (CTL) recognize complexes of peptide ligands and Major Histocompatibility Complex (MHC) class I molecules presented at the surface of Antigen Presenting Cells (APC). Detection and isolation of CTL''s are of importance for research on CTL immunity, and development of vaccines and adoptive immune therapy. Peptide-MHC tetramers have become important reagents for detection and enumeration of specific CTL''s. Conventional peptide-MHC-tetramer production involves recombinant MHC production, in vitro refolding, biotinylation and tetramerization; each step followed by various biochemical steps such as chromatographic purification, concentration etc. Such cumbersome production protocols have limited dissemination and restricted availability of peptide-MHC tetramers effectively precluding large-scale screening strategies involving many different peptide-MHC tetramers.

Methodology/Principal Findings

We have developed an approach whereby any given tetramer specificity can be produced within 2 days with very limited effort and hands-on time. The strategy is based on the isolation of correctly oxidized, in vivo biotinylated recombinant MHC I heavy chain (HC). Such biotinylated MHC I HC molecules can be refolded in vitro, tetramerized with streptavidin, and used for specific T cell staining-all in a one-pot reaction without any intervening purification steps.

Conclusions/Significance

We have developed an efficient “one-pot, mix-and-read” strategy for peptide-MHC tetramer generation, and demonstrated specific T cell straining comparable to a commercially available MHC-tetramer. Here, seven peptide-MHC tetramers representing four different human MHC (HLA) class I proteins have been generated. The technique should be readily extendable to any binding peptide and pre-biotinylated MHC (at this time we have over 40 different pre-biotinylated HLA proteins). It is simple, robust, and versatile technique with a very broad application potential as it can be adapted both to small- and large-scale production of one or many different peptide-MHC tetramers for T cell isolation, or epitope screening.     

An on-bead assay for the identification of non-natural peptides targeting the androgen receptor-cofactor interaction

An efficient and rapid on-bead screening method was established to identify non-natural peptides that target the Androgen Receptor-cofactor interaction. Binding of the Androgen Receptor ligand binding domain to peptide sequences displayed on beads in a One-Bead-One-Compound format could be screened using fluorescence microscopy. The method was applied to generate and screen both a focussed and a random peptide library. Resynthesis of the peptide hits allowed for the verification of the affinity of the selected peptides for the Androgen Receptor in a competitive fluorescence polarization assay. For both libraries strong Androgen Receptor binding peptides were found, both with non-natural and natural amino acids. The peptides identified with natural amino acids showed great similarity in terms of preferred amino acid sequence with peptides previously isolated from biological screens, thus validating the screening approach. The non-natural peptides featured important novel chemical transformations on the relevant hydrophobic amino acid positions interacting with the Androgen Receptor. This screening approach expands the molecular diversity of peptide inhibitors for nuclear receptors.     

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.     

Cutting edge: the conversion of arginine to citrulline allows for a high-affinity peptide interaction with the rheumatoid arthritis-associated HLA-DRB1*0401 MHC class II molecule

Rheumatoid arthritis (RA) is genetically associated with MHC class II molecules that contain the shared epitope. These MHC molecules may participate in disease pathogenesis by selectively binding arthritogenic peptides for presentation to autoreactive CD4(+) T cells. The nature of the arthritogenic Ag is not known, but recent work has identified posttranslationally modified proteins containing citrulline (deiminated arginine) as specific targets of the IgG Ab response in RA patients. To understand how citrulline might evoke an autoimmune reaction, we have studied T cell responses to citrulline-containing peptides in HLA-DRB1*0401 transgenic (DR4-IE tg) mice. In this study, we demonstrate that the conversion of arginine to citrulline at the peptide side-chain position interacting with the shared epitope significantly increases peptide-MHC affinity and leads to the activation CD4(+) T cells in DR4-IE tg mice. These results reveal how DRB1 alleles with the shared epitope could initiate an autoimmune response to citrullinated self-Ags in RA patients.     

A Comparative Approach Linking Molecular Dynamics of Altered Peptide Ligands and MHC with In Vivo Immune Responses

Background

The recognition of peptide in the context of MHC by T lymphocytes is a critical step in the initiation of an adaptive immune response. However, the molecular nature of the interaction between peptide and MHC and how it influences T cell responsiveness is not fully understood.

Results

We analyzed the immunological consequences of the interaction of MHC class II (I-A^u) restricted 11-mer peptides of myelin basic protein with amino acid substitutions at position 4. These mutant peptides differ in MHC binding affinity, CD4⁺ T cell priming, and alter the severity of peptide-induced experimental allergic encephalomyelitis. Using molecular dynamics, a computational method of quantifying intrinsic movements of proteins at high resolution, we investigated conformational changes in MHC upon peptide binding. We found that irrespective of peptide binding affinity, MHC deformation appears to influence costimulation, which then leads to effective T cell priming and disease induction. Although this study compares in vivo and molecular dynamics results for three altered peptide ligands, further investigation with similar complexes is essential to determine whether spatial rearrangement of peptide-MHC and costimulatory complexes is an additional level of T cell regulation.     

Contribution of individual amino acids within MHC molecule or antigenic peptide to TCR ligand potency

The TCR recognition of peptides bound to MHC class II molecules is highly flexible in some T cells. Although progress has been made in understanding the interactions within the trimolecular complex, to what extent the individual components and their amino acid composition contribute to ligand recognition by individual T cells is not completely understood. We investigated how single amino acid residues influence Ag recognition of T cells by combining several experimental approaches. We defined TCR motifs for CD4+ T cells using peptide synthetic combinatorial libraries in the positional scanning format (PS-SCL) and single amino acid-modified peptide analogues. The similarity of the TCR motifs defined by both methods and the identification of stimulatory antigenic peptides by the PS-SCL approach argue for a contribution of each amino acid residue to the overall potency of the antigenic peptide ligand. In some instances, however, motifs are formed by adjacent amino acids, and their combined influence is superimposed on the overall contribution of each amino acid within the peptide epitope. In contrast to the flexibility of the TCR to interact with different peptides, recognition was very sensitive toward modifications of the MHC-restriction element. Exchanges of just one amino acid of the MHC molecule drastically reduced the number of peptides recognized. The results indicate that a specific MHC molecule not only selects certain peptides, but also is crucial for setting an affinity threshold for TCR recognition, which determines the flexibility in peptide recognition for a given TCR.     

A structure-based approach to designing non-natural peptides that can activate anti-melanoma cytotoxic T cells

Tumor antigens presented by major histocompatibility complex (MHC) class I molecules and recognized by CD8(+) cytotoxic T lymphocytes (CTLs) may generate an efficient antitumor immune response after appropriate immunization. Antigenic peptides can be used in vivo to induce antitumor or antiviral immunity. The efficiency of naked peptides may be greatly limited by their degradation in the biological fluids. We present a rational, structure-based approach to design structurally modified, peptidase-resistant and biologically active analogues of human tumor antigen MAGE-1.A1. This approach is based on our understanding of the peptide interaction with the MHC and the T cell receptor and its precise degradation pathway. Knowledge of these mechanisms led to the design of a non-natural, minimally modified analogue of MAGE-1.A1, [Aib2, NMe-Ser8]MAGE-1.A1, which was highly peptidase-resistant and bound to MHC and activated MAGE-1.A1-specific anti-melanoma CTLs. Thus, we showed that it is possible to structurally modify peptide epitopes to obtain analogues that are still specifically recognized by CTLs. Such analogues may represent interesting leads for antitumor synthetic vaccines.     

Binding of rationally designed non-natural peptides to the human leukocyte antigen HLA-B*2705

High-affinity ligands of non-peptidic nature, binding to the class I major histocompatibility complex protein HLA B*2705 whose expression is strongly linked to the pathogenesis of the autoimmune disease ankylosing spondylitis, should give way to a selective immunotherapy by blocking or antagonising the interaction with autoreactive T cell clones. Here we present experimental data on the binding of modified peptides, designed to optimally bind to HLA-B*2705 by filling a hydrophobic binding pocket (pocket D) with nonencoded aromatic amino acids. Three peptides with altered side chains (alpha-naphthylalanine, beta-naphthylalanine and homophenylalanine) in position 3 were synthesised. The thermal denaturation profiles of the HLA protein in complex with the modified peptides, monitored by circular dichroism spectroscopy, showed a significant shift towards higher melting temperatures with respect to the parent T cell epitope. The proposed binding mode of the nonnatural peptides was checked by site-directed mutagenesis of the pocket D, hypothesised to accommodate the large hydrophobic side chains. Reducing the size and depth of the pocket by mutating Leu l56 into Trp only affects the binding of the non-natural ligands, thus providing experimental evidence that the nonnatural peptide amino acids bind as predicted to the host MHC protein. © 1998 European Peptide Society and John Wiley & Sons, Ltd.     

烟曲霉抗原Asp f16HLA-A*0201限制性的CD8+CTL抗原表位生物信息学预测与实验室鉴定

目的 预测与鉴定烟曲霉抗原Asp f16的HLA-A *0201限制性CD8+细胞毒性T细胞(CTL)抗原表位.方法 以国人常见的HLA-A*0201位点为靶点,依据生物信息学软件扫描烟曲霉特异性抗原Asp f16的全部427个氨基酸序列.使用HLA-A *0201转基因小鼠制备骨髓来源的树突状细胞(DC)和CTL.流式细胞仪技术检测DC表面MHC Ⅱ类抗原,CD80,CD86和CD11c的表达来验证其是否成熟.ELISPOT试验检测烟曲霉抗原多肽特异性CTL产生的细胞因子IFN-γ.四聚体(Tetramer)试验证实烟曲霉特异性CTL与抗原肽,HLA-A*0201分子复合体的亲和性.结果 根据与MHC I类分子结合的半衰期评分,选择了3个HLA-A*0201限制性抗原表位.流式细胞仪分析示成熟DC高表达HLA Ⅱ类抗原,CD80,CD86和CD11c.Tetramer试验证实烟曲霉特异性T细胞受体与抗原肽,HLA-A*0201分子复合体的高亲和性.ELISPOT实验结果 表明烟曲霉抗原肽体外可以活化CD8+CTL,被负载了抗原肽的DC刺激活化后可以产生IFN-γ.结论 本研究成功鉴定烟曲霉抗原Asp f16的HLA-A*0201限制性CD8+CTL表位,可作为疫苗设计的候选表位,为进一步研发新型抗烟曲霉疫苗提供参考.