MHCPred: A server for quantitative prediction of peptide-MHC binding

Accurate T-cell epitope prediction is a principal objective of computational vaccinology. As a service to the immunology and vaccinology communities at large, we have implemented, as a server on the World Wide Web, a partial least squares-based multivariate statistical approach to the quantitative prediction of peptide binding to major histocom- patibility complexes (MHC), the key checkpoint on the antigen presentation pathway within adaptive cellular immunity. MHCPred implements robust statistical models for both Class I alleles (HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3301, HLA-A*6801, HLA-A*6802 and HLA-B*3501) and Class II alleles (HLA-DRB*0401, HLA-DRB*0401 and HLA-DRB*0701). MHCPred is available from the URL: http://www.jenner.ac.uk/MHCPred.     

Prediction of promiscuous peptides that bind HLA class I molecules

Promiscuous T-cell epitopes make ideal targets for vaccine development. We report here a computational system, MULTIPRED, for the prediction of peptide binding to the HLA-A2 supertype. It combines a novel representation of peptide/MHC interactions with a hidden Markov model as the prediction algorithm. MULTIPREDis both sensitive and specific, and demonstrates high accuracy of peptide-binding predictions for HLA-A*0201, *0204, and *0205 alleles, good accuracy for *0206 allele, and marginal accuracy for *0203 allele. MULTIPREDreplaces earlier requirements for individual prediction models for each HLA allelic variant and simplifies computational aspects of peptide-binding prediction. Preliminary testing indicates that MULTIPRED can predict peptide binding to HLA-A2 supertype molecules with high accuracy, including those allelic variants for which no experimental binding data are currently available.     

Transporter associated with antigen processing preselection of peptides binding to the MHC: a bioinformatic evaluation

TAP is responsible for the transit of peptides from the cytosol to the lumen of the endoplasmic reticulum. In an immunological context, this event is followed by the binding of peptides to MHC molecules before export to the cell surface and recognition by T cells. Because TAP transport precedes MHC binding, TAP preferences may make a significant contribution to epitope selection. To assess the impact of this preselection, we have developed a scoring function for TAP affinity prediction using the additive method, have used it to analyze and extend the TAP binding motif, and have evaluated how well this model acts as a preselection step in predicting MHC binding peptides. To distinguish between MHC alleles that are exclusively dependent on TAP and those exhibiting only a partial dependence on TAP, two sets of MHC binding peptides were examined: HLA-A*0201 was selected as a representative of partially TAP-dependent HLA alleles, and HLA-A*0301 represented fully TAP-dependent HLA alleles. TAP preselection has a greater impact on TAP-dependent alleles than on TAP-independent alleles. The reduction in the number of nonbinders varied from 10% (TAP-independent) to 33% (TAP-dependent), suggesting that TAP preselection is an important component in the successful in silico prediction of T cell epitopes.     

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.     

Improving MHC binding peptide prediction by incorporating binding data of auxiliary MHC molecules

MOTIVATION: Various computational methods have been proposed to tackle the problem of predicting the peptide binding ability for a specific MHC molecule. These methods are based on known binding peptide sequences. However, current available peptide databases do not have very abundant amounts of examples and are highly redundant. Existing studies show that MHC molecules can be classified into supertypes in terms of peptide-binding specificities. Therefore, we first give a method for reducing the redundancy in a given dataset based on information entropy, then present a novel approach for prediction by learning a predictive model from a dataset of binders for not only the molecule of interest but also for other MHC molecules. RESULTS: We experimented on the HLA-A family with the binding nonamers of A1 supertype (HLA-A*0101, A*2601, A*2902, A*3002), A2 supertype (A*0201, A*0202, A*0203, A*0206, A*6802), A3 supertype (A*0301, A*1101, A*3101, A*3301, A*6801) and A24 supertype (A*2301 and A*2402), whose data were collected from six publicly available peptide databases and two private sources. The results show that our approach significantly improves the prediction accuracy of peptides that bind a specific HLA molecule when we combine binding data of HLA molecules in the same supertype. Our approach can thus be used to help find new binders for MHC molecules.     

Integrated modeling of the major events in the MHC class I antigen processing pathway

Rational design of epitope-driven vaccines is a key goal of immunoinformatics. Typically, candidate selection relies on the prediction of MHC-peptide binding only, as this is known to be the most selective step in the MHC class I antigen processing pathway. However, proteasomal cleavage and transport by the transporter associated with antigen processing (TAP) are essential steps in antigen processing as well. While prediction methods exist for the individual steps, no method has yet offered an integrated prediction of all three major processing events. Here we present WAPP, a method combining prediction of proteasomal cleavage, TAP transport, and MHC binding into a single prediction system. The proteasomal cleavage site prediction employs a new matrix-based method that is based on experimentally verified proteasomal cleavage sites. Support vector regression is used for predicting peptides transported by TAP. MHC binding is the last step in the antigen processing pathway and was predicted using a support vector machine method, SVMHC. The individual methods are combined in a filtering approach mimicking the natural processing pathway. WAPP thus predicts peptides that are cleaved by the proteasome at the C terminus, transported by TAP, and show significant affinity to MHC class I molecules. This results in a decrease in false positive rates compared to MHC binding prediction alone. Compared to prediction of MHC binding only, we report an increased overall accuracy and a lower rate of false positive predictions for the HLA-A*0201, HLA-B*2705, HLA-A*01, and HLA-A*03 alleles using WAPP. The method is available online through our prediction server at http://www-bs.informatik.uni-tuebingen.de/WAPP     

MHCPred 2.0

The accurate computational prediction of T-cell epitopes can greatly reduce the experimental overhead implicit in candidate epitope identification within genomic sequences. In this article we present MHCPred 2.0, an enhanced version of our online, quantitative T-cell epitope prediction server. The previous version of MHCPred included mostly alleles from the human leukocyte antigen A (HLA-A) locus. In MHCPred 2.0, mouse models are added and computational constraints removed. Currently the server includes 11 human HLA class I, three human HLA class II, and three mouse class I models. Additionally, a binding model for the human transporter associated with antigen processing (TAP) is incorporated into the new MHCPred. A tool for the design of heteroclitic peptides is also included within the server. To refine the veracity of binding affinities prediction, a confidence percentage is also now calculated for each peptide predicted. AVAILABILITY: As previously, MHCPred 2.0 is freely available at the URL http://www.jenner.ac.uk/MHCPred/ CONTACT: Darren R. Flower (darren.flower@jenner.ac.uk).     

Random screening of proteins for HLA-A*0201-binding nine-amino acid peptides is not sufficient for identifying CD8 T cell epitopes recognized in the context of HLA-A*0201

HLA-A2 is the most frequent HLA molecule in Caucasians with HLA-A*0201 representing the most frequent allele; it was also the first human HLA allele for which peptide binding prediction was developed. The Bioinformatics and Molecular Analysis Section of the National Institutes of Health (BIMAS) and the University of Tübingen (Syfpeithi) provide the most popular prediction algorithms of peptide/MHC interaction on the World Wide Web. To test these predictions, HLA-A*0201-binding nine-amino acid peptides were searched by both algorithms in 19 structural CMV proteins. According to Syfpeithi, the top 2% of predicted peptides should contain the naturally presented epitopes in 80% of predictions (www.syfpeithi.de). Because of the high number of predicted peptides, the analysis was limited to 10 randomly chosen proteins. The top 2% of peptides predicted by both algorithms were synthesized corresponding to 261 peptides in total. PBMC from 10 HLA-A*0201-positive and CMV-seropositive healthy blood donors were tested by ex vivo stimulation with all 261 peptides using crossover peptide pools. IFN-gamma production in T cells measured by CFC was used as readout. However, only one peptide was found to be stimulating in one single donor. As a result of this work, we report a potential new T cell target protein, one previously unknown CD8-T cell-stimulating peptide, and an extensive list of CMV-derived potentially strong HLA-A*0201-binding peptides that are not recognized by T cells of HLA-A*0201-positive CMV-seropositive donors. We conclude that MHC/peptide binding predictions are helpful for locating epitopes in known target proteins but not necessarily for screening epitopes in proteins not known to be T cell targets.     

Toward the quantitative prediction of T-cell epitopes: QSAR studies on peptides having affinity with the class I MHC molecular HLA-A*0201.

It would be useful for vaccine development to develop a method of rapidly identifying peptide epitopes. In this paper, the empirical three-dimensional quantitative structure-affinity relationship (3D-QSAR) methods were used to study the relationship between the three dimensional structural parameters (the isotropic surface area, ISA, and the electronic charge index, ECI) of the HLA-A*0201 binding peptide and the HLA-A*0201/peptide binding affinities. A set of 102 peptides having affinity with the class I MHC HLA-A*0201 molecule was used as training set. A test set of 40 peptides was used to determine the predictive value of the models. The 3D-QSAR models yielded a q2 = 0.5724 and a high rpred2 = 0.6955. The standard regression coefficients indicated that the hydrophobic interactions played an important role in peptide-MHC molecule binding and predicted the specific amino acid residue essential at a certain position of the peptide. The approach tested in the current paper is highly complementary to many of the methods described in references and possesses good predictability. It is a rapid and convenient method to detect high affinity peptide epitopes.     

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.     

Epitope-Based Peptides Prediction from Proteome of Nipah Virus

The identification of MHC class II epitope-based peptides are urgently needed for appropriate vaccination against Nipah virus (NiV) because there are currently no approved vaccines for human NiV infection. In the present study, prediction and modeling of T cell epitopes of NiV antigenic proteins nucleocapsid, phosphoprotein, matrix, fusion, glycoprotein, L protein, W protein, V protein and C protein followed by the binding simulation studies of predicted highest binding scores with their corresponding MHC class II alleles were done. Immunoinformatic tool ProPred was used to predict the promiscuous MHC class II epitopes of viral antigenic proteins. PEPstr server did the 3D structure models of the epitopes and Modeller 9.10 did alleles. We docked epitope with allele structure using the AutoDock 4.2 Tool. The docked peptide–allele complex structure was optimized using molecular dynamics simulation for 5 ps with the CHARMM-22 force field using NAnoscale Molecular Dynamics program incorporated in visual molecular dynamics (VMD 1.9.2) and then evaluating the stability of complex structure by calculating RMSD values. Epitope MKLQFSLGS of Matrix protein has considerable binding energy and score with DRBI*0421 MHC class II allele. This predicted peptide has potential to induce T cell-mediated immune response and is expected to useful in designing epitope-based vaccines against NiV after further testing by wet lab studies.     

Stepwise identification of HLA-A*0201-restricted CD8+ T-cell epitope peptides from herpes simplex virus type 1 genome boosted by a StepRank scheme

Identification of immunodominant epitopes is the first step in the rational design of peptide vaccines aimed at T-cell immunity. To date, however, it is yet a great challenge for accurately predicting the potent epitope peptides from a pool of large-scale candidates with an efficient manner. In this study, a method that we named StepRank has been developed for the reliable and rapid prediction of binding capabilities/affinities between proteins and genome-wide peptides. In this procedure, instead of single strategy used in most traditional epitope identification algorithms, four steps with different purposes and thus different computational demands are employed in turn to screen the large-scale peptide candidates that are normally generated from, for example, pathogenic genome. The steps 1 and 2 aim at qualitative exclusion of typical nonbinders by using empirical rule and linear statistical approach, while the steps 3 and 4 focus on quantitative examination and prediction of the interaction energy profile and binding affinity of peptide to target protein via quantitative structure-activity relationship (QSAR) and structure-based free energy analysis. We exemplify this method through its application to binding predictions of the peptide segments derived from the 76 known open-reading frames (ORFs) of herpes simplex virus type 1 (HSV-1) genome with or without affinity to human major histocompatibility complex class I (MHC I) molecule HLA-A*0201, and find that the predictive results are well compatible with the classical anchor residue theory and perfectly match for the extended motif pattern of MHC I-binding peptides. The putative epitopes are further confirmed by comparisons with 11 experimentally measured HLA-A*0201-restrcited peptides from the HSV-1 glycoproteins D and K. We expect that this well-designed scheme can be applied in the computational screening of other viral genomes as well.     

HLA-A3 supermotif defined by quantitative structure-activity relationship analysis

Activation of a cytotoxic T cell requires specific binding of antigenic peptides to major histocompatibility complex (MHC) molecules. This paper reports a study of peptides binding to members of the HLA-A3 superfamily using a recently developed 2D-QSAR method, called the additive method. Four alleles with high phenotype frequency were included in the study: A*0301, A*1101, A*3101 and A*6801. The influence of each of the 20 amino acids at each position of the peptide on binding was studied. A refined A3 supertype motif was defined in the study.     

MHC binding peptides for designing of vaccines against Japanese encephalitis virus: A computational approach

Japanese encephalitis (JE), a viral disease has seen a drastic and fatal enlargement in the northern states of India in the current decade. The better and exact cure for the disease is still in waiting. For the cause an in silico strategy in the development of the peptide vaccine has been taken here for the study. A computational approach to find out the Major Histocompatibility Complex (MHC) binding peptide has been implemented. The prediction analysis identified MHC class I (using propred I) and MHC class II (using propred) binding peptides at an expectable percent predicted IC (50) threshold values. These predicted Human leukocyte antigen [HLA] allele binding peptides were further analyzed for potential conserved region using an Immune Epitope Database and Analysis Resource (IEDB). This analysis shows that HLA-DRB1*0101, HLA-DRB3*0101, HLA-DRB1*0401, HLA-DRB1*0102 and HLA-DRB1*07:01% of class II (in genotype 2) and HLA-A*0101, HLA-A*02, HLA-A*0301, HLA-A*2402, HLA-B*0702 and HLA-B*4402% of HLA I (in genotype 3) bound peptides are conserved. The predicted peptides MHC class I are ILDSNGDIIGLY, FVMDEAHFTDPA, KTRKILPQIIK, RLMSPNRVPNYNLF, APTRVVAAEMAEAL, YENVFHTLW and MHC class II molecule are TTGVYRIMARGILGT, NYNLFVMDEAHFTDP, AAAIFMTATPPGTTD, GDTTTGVYRIMARGI and FGEVGAVSL found to be top ranking with potential super antigenic property by binding to all HLA. Out of these the predicted peptide FVMDEAHFTDPA for allele HLA-A*02:01 in MHC class I and NYNLFVMDEAHFTDP for allele HLA-DRB3*01:01 in MHC class II was observed to be most potent and can be further proposed as a significant vaccine in the process. The reported results revealed that the immune-informatics techniques implemented in the development of small size peptide is useful in the development of vaccines against the Japanese encephalitis virus (JEV).     

A comparative molecular similarity indices (CoMSIA) study of peptide binding to the HLA-A3 superfamily

Epitope identification is the basis of modern vaccine design. The present paper studied the supermotif of the HLA-A3 superfamily, using comparative molecular similarity indices analysis (CoMSIA). Four alleles with high phenotype frequencies were used: A*1101, A*0301, A*3101 and A*6801. Five physicochemical properties-steric bulk, electrostatic potential, local hydrophobicity, hydrogen-bond donor and acceptor abilities-were considered and 'all fields' models were produced for each of the alleles. The models have a moderate level of predictivity and there is a good correlation between the data. A revised HLA-A3 supermotif was defined based on the comparison of favoured and disfavoured properties for each position of the MHC bound peptide. The present study demonstrated that CoMSIA is an effective tool for studying peptide-MHC interactions.     

Synthesis and ex vivo profiling of chemically modified cytomegalovirus CMVpp65 epitopes.

The effect of substituting unnatural hydrophobic amino acids into the critical MHC binding residues of an HLA-A*0201-restricted cytomegalovirus CMVpp65 epitope, NLVPMVATV, has been investigated. A new set of peptides containing the amino acids tert-butyl glycine (Tgl), cyclohexyl glycine (Chg), neo-pentyl glycine (Npg), cyclohexyl alanine (Cha) and cyclo leucine (Cyl), at either position 2, to mimic Leu, or position 9, to mimic Val, have been synthesised. Immunological profiling using class I MHC stabilisation assays to assess MHC binding affinity, and enzyme-linked immunospot (ELISPOT) assays to assess the ability of the modified peptides to re-stimulate a specific cytotoxic T-lymphocyte (CTL) response, compared to the native epitope, have been performed. It was found that the majority of the unnatural substitutions resulted in a decrease in either HLA-A*0201 binding affinity or cytotoxic T-cell activity. However, the HLA-A*0201 binding affinity was unrelated to the ability to re-stimulate a T-cell response. Minimisation and molecular dynamics studies proved helpful in dissecting the ELISPOT responses. Two principal peptide binding modes were found by minimisation, designated kinked and straight. Peptides that bound in a kinked conformation were poor at re-stimulating a T-cell response. Of the peptides that bound in a straight conformation, molecular dynamics (MD) simulations revealed that those capable of re-stimulating the strongest responses had the greatest degree of flexibility (as determined by RMSD values across the MD simulation) around the P6 residue, one of the residues important for T-cell receptor recognition.     

HLA class I-associated diseases with a suspected autoimmune etiology: HLA-B27 subtypes as a model system

Although most autoimmune diseases are connected to major histocompatibility complex (MHC) class II alleles, a small number of these disorders exhibit a variable degree of association with selected MHC class I genes, like certain human HLA-A and HLA-B alleles. The basis for these associations, however, has so far remained elusive. An understanding might be obtained by comparing functional, biochemical, and biophysical properties of alleles that are minimally distinct from each other, but are nevertheless differentially associated to a given disease, like the HLA-B*27:05 and HLA-B*27:09 antigens, which differ only by a single amino acid residue (Asp116His) that is deeply buried within the binding groove. We have employed a number of approaches, including X-ray crystallography and isotope-edited infrared spectroscopy, to investigate biophysical characteristics of the two HLA-B27 subtypes complexed with up to ten different peptides. Our findings demonstrate that the binding of these peptides as well as the conformational flexibility of the subtypes is greatly influenced by interactions of the C-terminal peptide residue. In particular, a basic C-terminal peptide residue is favoured by the disease-associated subtype HLA-B*27:05, but not by HLA-B*27:09. This property appears also as the only common denominator of distinct HLA class I alleles, among them HLA-B*27:05, HLA-A*03:01 or HLA-A*11:01, that are associated with diseases suspected to have an autoimmune etiology. We postulate here that the products of these alleles, due to their unusual ability to bind with high affinity to a particular peptide set during positive T cell selection in the thymus, are involved in shaping an abnormal T cell repertoire which predisposes to the acquisition of autoimmune diseases.