Knowledge-based structure prediction of MHC class I bound peptides: a study of 23 complexes

Physicochemical Heuristics for Identifying High Fidelity, Near-Native Structural Models of Peptide/MHC Complexes

There is long-standing interest in accurately modeling the structural features of peptides bound and presented by class I MHC proteins. This interest has grown with the advent of rapid genome sequencing and the prospect of personalized, peptide-based cancer vaccines, as well as the development of molecular and cellular therapeutics based on T cell receptor recognition of peptide-MHC. However, while the speed and accessibility of peptide-MHC modeling has improved substantially over the years, improvements in accuracy have been modest. Accuracy is crucial in peptide-MHC modeling, as T cell receptors are highly sensitive to peptide conformation and capturing fine details is therefore necessary for useful models. Studying nonameric peptides presented by the common class I MHC protein HLA-A*02:01, here we addressed a key question common to modern modeling efforts: from a set of models (or decoys) generated through conformational sampling, which is best? We found that the common strategy of decoy selection by lowest energy can lead to substantial errors in predicted structures. We therefore adopted a data-driven approach and trained functions capable of predicting near native decoys with exceptionally high accuracy. Although our implementation is limited to nonamer/HLA-A*02:01 complexes, our results serve as an important proof of concept from which improvements can be made and, given the significance of HLA-A*02:01 and its preference for nonameric peptides, should have immediate utility in select immunotherapeutic and other efforts for which structural information would be advantageous.

Structural Prediction of Peptide-MHC Binding Modes

The immune system is constantly protecting its host from the invasion of pathogens and the development of cancer cells. The specific CD8⁺ T-cell immune response against virus-infected cells and tumor cells is based on the T-cell receptor recognition of antigenic peptides bound to class I major histocompatibility complexes (MHC) at the surface of antigen presenting cells. Consequently, the peptide binding specificities of the highly polymorphic MHC have important implications for the design of vaccines, for the treatment of autoimmune diseases, and for personalized cancer immunotherapy. Evidence-based machine-learning approaches have been successfully used for the prediction of peptide binders and are currently being developed for the prediction of peptide immunogenicity. However, understanding and modeling the structural details of peptide/MHC binding is crucial for a better understanding of the molecular mechanisms triggering the immunological processes, estimating peptide/MHC affinity using universal physics-based approaches, and driving the design of novel peptide ligands. Unfortunately, due to the large diversity of MHC allotypes and possible peptides, the growing number of 3D structures of peptide/MHC (pMHC) complexes in the Protein Data Bank only covers a small fraction of the possibilities. Consequently, there is a growing need for rapid and efficient approaches to predict 3D structures of pMHC complexes. Here, we review the key characteristics of the 3D structure of pMHC complexes before listing databases and other sources of information on pMHC structures and MHC specificities. Finally, we discuss some of the most prominent pMHC docking software.

Structure Based Prediction of Neoantigen Immunogenicity

The development of immunological therapies that incorporate peptide antigens presented to T cells by MHC proteins is a long sought-after goal, particularly for cancer, where mutated neoantigens are being explored as personalized cancer vaccines. Although neoantigens can be identified through sequencing, bioinformatics and mass spectrometry, identifying those which are immunogenic and able to promote tumor rejection remains a significant challenge. Here we examined the potential of high-resolution structural modeling followed by energetic scoring of structural features for predicting neoantigen immunogenicity. After developing a strategy to rapidly and accurately model nonameric peptides bound to the common class I MHC protein HLA-A2, we trained a neural network on structural features that influence T cell receptor (TCR) and peptide binding energies. The resulting structurally-parameterized neural network outperformed methods that do not incorporate explicit structural or energetic properties in predicting CD8+ T cell responses of HLA-A2 presented nonameric peptides, while also providing insight into the underlying structural and biophysical mechanisms governing immunogenicity. Our proof-of-concept study demonstrates the potential for structure-based immunogenicity predictions in the development of personalized peptide-based vaccines.

A stochastic individual-based model to explore the role of spatial interactions and antigen recognition in the immune response against solid tumours

Spatial interactions between cancer and immune cells, as well as the recognition of tumour antigens by cells of the immune system, play a key role in the immune response against solid tumours. The existing mathematical models generally focus only on one of these key aspects. We present here a spatial stochastic individual-based model that explicitly captures antigen expression and recognition. In our model, each cancer cell is characterised by an antigen profile which can change over time due to either epimutations or mutations. The immune response against the cancer cells is initiated by the dendritic cells that recognise the tumour antigens and present them to the cytotoxic T cells. Consequently, T cells become activated against the tumour cells expressing such antigens. Moreover, the differences in movement between inactive and active immune cells are explicitly taken into account by the model. Computational simulations of our model clarify the conditions for the emergence of tumour clearance, dormancy or escape, and allow us to assess the impact of antigenic heterogeneity of cancer cells on the efficacy of immune action. Ultimately, our results highlight the complex interplay between spatial interactions and adaptive mechanisms that underpins the immune response against solid tumours, and suggest how this may be exploited to further develop cancer immunotherapies.

Structure-based Methods for Binding Mode and Binding Affinity Prediction for Peptide-MHC Complexes

Understanding the mechanisms involved in the activation of an immune response is essential to many fields in human health, including vaccine development and personalized cancer immunotherapy. A central step in the activation of the adaptive immune response is the recognition, by T-cell lymphocytes, of peptides displayed by a special type of receptor known as Major Histocompatibility Complex (MHC). Considering the key role of MHC receptors in T-cell activation, the computational prediction of peptide binding to MHC has been an important goal for many immunological applications. Sequence- based methods have become the gold standard for peptide-MHC binding affinity prediction, but structure-based methods are expected to provide more general predictions (i.e., predictions applicable to all types of MHC receptors). In addition, structural modeling of peptide-MHC complexes has the potential to uncover yet unknown drivers of T-cell activation, thus allowing for the development of better and safer therapies. In this review, we discuss the use of computational methods for the structural modeling of peptide-MHC complexes (i.e., binding mode prediction) and for the structure-based prediction of binding affinity.

General Prediction of Peptide-MHC Binding Modes Using Incremental Docking: A Proof of Concept

The class I major histocompatibility complex (MHC) is capable of binding peptides derived from intracellular proteins and displaying them at the cell surface. The recognition of these peptide-MHC (pMHC) complexes by T-cells is the cornerstone of cellular immunity, enabling the elimination of infected or tumoral cells. T-cell-based immunotherapies against cancer, which leverage this mechanism, can greatly benefit from structural analyses of pMHC complexes. Several attempts have been made to use molecular docking for such analyses, but pMHC structure remains too challenging for even state-of-the-art docking tools. To overcome these limitations, we describe the use of an incremental meta-docking approach for structural prediction of pMHC complexes. Previous methods applied in this context used specific constraints to reduce the complexity of this prediction problem, at the expense of generality. Our strategy makes no assumption and can potentially be used to predict binding modes for any pMHC complex. Our method has been tested in a re-docking experiment, reproducing the binding modes of 25 pMHC complexes whose crystal structures are available. This study is a proof of concept that incremental docking strategies can lead to general geometry prediction of pMHC complexes, with potential applications for immunotherapy against cancer or infectious diseases.

Emerging Concepts in TCR Specificity: Rationalizing and (Maybe) Predicting Outcomes

T cell specificity emerges from a myriad of processes, ranging from the biological pathways that control T cell signaling to the structural and physical mechanisms that influence how TCRs bind peptides and MHC proteins. Of these processes, the binding specificity of the TCR is a key component. However, TCR specificity is enigmatic: TCRs are at once specific but also cross-reactive. Although long appreciated, this duality continues to puzzle immunologists and has implications for the development of TCR-based therapeutics. In this review, we discuss TCR specificity, emphasizing results that have emerged from structural and physical studies of TCR binding. We show how the TCR specificity/cross-reactivity duality can be rationalized from structural and biophysical principles. There is excellent agreement between predictions from these principles and classic predictions about the scope of TCR cross-reactivity. We demonstrate how these same principles can also explain amino acid preferences in immunogenic epitopes and highlight opportunities for structural considerations in predictive immunology.

Computational Modeling of T Cell Receptor Complexes

T-cell receptor (TCR) binding to peptide/MHC determines specificity and initiates signaling in antigen-specific cellular immune responses. Structures of TCR-pMHC complexes have provided enormous insight to cellular immune functions, permitted a rational understanding of processes such as pathogen escape, and led to the development of novel approaches for the design of vaccines and other therapeutics. As production, crystallization, and structure determination of TCR-pMHC complexes can be challenging, there is considerable interest in modeling new complexes. Here we describe a rapid approach to TCR-pMHC modeling that takes advantage of structural features conserved in known complexes, such as the restricted TCR binding site and the generally conserved diagonal docking mode. The approach relies on the powerful Rosetta suite and is implemented using the PyRosetta scripting environment. We show how the approach can recapitulate changes in TCR binding angles and other structural details, and highlight areas where careful evaluation of parameters is needed and alternative choices might be made. As TCRs are highly sensitive to subtle structural perturbations, there is room for improvement. Our method nonetheless generates high-quality models that can be foundational for structure-based hypotheses regarding TCR recognition.

Quantitative prediction of peptide binding to HLA-DP1 protein

The exogenous proteins are processed by the host antigen-processing cells. Peptidic fragments of them are presented on the cell surface bound to the major hystocompatibility complex (MHC) molecules class II and recognized by the CD4+ T lymphocytes. The MHC binding is considered as the crucial prerequisite for T-cell recognition. Only peptides able to form stable complexes with the MHC proteins are recognized by the T-cells. These peptides are known as T-cell epitopes. All T-cell epitopes are MHC binders, but not all MHC binders are T-cell epitopes. The T-cell epitope prediction is one of the main priorities of immunoinformatics. In the present study, three chemometric techniques are combined to derive a model for in silico prediction of peptide binding to the human MHC class II protein HLA-DP1. The structures of a set of known peptide binders are described by amino acid z-descriptors. Data are processed by an iterative self-consisted algorithm using the method of partial least squares, and a quantitative matrix (QM) for peptide binding prediction to HLA-DP1 is derived. The QM is validated by two sets of proteins and showed an average accuracy of 86 percent.

A combination of epitope prediction and molecular docking allows for good identification of MHC class I restricted T-cell epitopes

•

Combining epitope prediction methods with molecular docking techniques to identify MHC class I restricted T-cell epitopes.

•

Based on available experimental data, the prediction accuracy is up to 90%.

•

Providing a valuable step forward for the design of better vaccines.

•

Better understanding the activation of T-cell epitopes by MHC binding peptides.

In silico identification of T-cell epitopes is emerging as a new methodology for the study of epitope-based vaccines against viruses and cancer. In order to improve accuracy of prediction, we designed a novel approach, using epitope prediction methods in combination with molecular docking techniques, to identify MHC class I restricted T-cell epitopes. Analysis of the HIV-1 p24 protein and influenza virus matrix protein revealed that the present approach is effective, yielding prediction accuracy of over 80% with respect to experimental data. Subsequently, we applied such a method for prediction of T-cell epitopes in SARS coronavirus (SARS-CoV) S, N and M proteins. Based on available experimental data, the prediction accuracy is up to 90% for S protein. We suggest the use of epitope prediction methods in combination with 3D structural modelling of peptide-MHC-TCR complex to identify MHC class I restricted T-cell epitopes for use in epitope based vaccines like HIV and human cancers, which should provide a valuable step forward for the design of better vaccines and may provide in depth understanding about activation of T-cell epitopes by MHC binding peptides.

Computational Resources for MHC Ligand Identification

Advances in the high-throughput determination of functional modulators of major histocompatibility complex (MHC) and improved computational predictions of MHC ligands have rendered the rational design of immunomodulatory peptides feasible. Proteome-derived peptides and 'reverse vaccinology' by computational means will play a driving role in future vaccine design. Here we review the molecular mechanisms of the MHC mediated immune response, present the computational approaches that have emerged in this area of biotechnology, and provide an overview of publicly available computational resources for predicting and designing new peptidic MHC ligands.

T-cell epitope vaccine design by immunoinformatics

Vaccination is generally considered to be the most effective method of preventing infectious diseases. All vaccinations work by presenting a foreign antigen to the immune system in order to evoke an immune response. The active agent of a vaccine may be intact but inactivated ('attenuated') forms of the causative pathogens (bacteria or viruses), or purified components of the pathogen that have been found to be highly immunogenic. The increased understanding of antigen recognition at molecular level has resulted in the development of rationally designed peptide vaccines. The concept of peptide vaccines is based on identification and chemical synthesis of B-cell and T-cell epitopes which are immunodominant and can induce specific immune responses. The accelerating growth of bioinformatics techniques and applications along with the substantial amount of experimental data has given rise to a new field, called immunoinformatics. Immunoinformatics is a branch of bioinformatics dealing with in silico analysis and modelling of immunological data and problems. Different sequence- and structure-based immunoinformatics methods are reviewed in the paper.

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

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

Ranking of binding and nonbinding peptides to MHC class I molecules using inverse folding approach: implications for vaccine design

T cell recognition of the peptide-MHC complex initiates a cascade of immunological events necessary for immune responses. Accurate T-cell epitope prediction is an important part of the vaccine designing. Development of predictive algorithms based on sequence profile requires a very large number of experimental binding peptide data to major histocompatibility complex (MHC) molecules. Here we used inverse folding approach to study the peptide specificity of MHC Class-I molecule with the aim of obtaining a better differentiation between binding and nonbinding sequence. Overlapping peptides, spanning the entire protein sequence, are threaded through the backbone coordinates of a known peptide fold in the MHC groove, and their interaction energies are evaluated using statistical pairwise contact potentials. We used the Miyazawa & Jernigan and Betancourt & Thirumalai tables for pairwise contact potentials, and two distance criteria (Nearest atom >> 4.0 A & C-beta >> 7.0 A) for ranking the peptides in an ascending order according to their energy values, and in most cases, known antigenic peptides are highly ranked. The predictions from threading improved when used multiple templates and average scoring scheme. In general, when structural information about a protein-peptide complex is available, the current application of the threading approach can be used to screen a large library of peptides for selection of the best binders to the target protein. The proposed scheme may significantly reduce the number of peptides to be tested in wet laboratory for epitope based vaccine design.

MODPROPEP: a program for knowledge-based modeling of protein-peptide complexes

MODPROPEP is a web server for knowledge-based modeling of protein–peptide complexes, specifically peptides in complex with major histocompatibility complex (MHC) proteins and kinases. The available crystal structures of protein–peptide complexes in PDB are used as templates for modeling peptides of desired sequence in the substrate-binding pocket of MHCs or protein kinases. The substrate peptides are modeled using the same backbone conformation as in the template and the side-chain conformations are obtained by the program SCWRL. MODPROPEP provides a number of user-friendly interfaces for visualizing the structure of the modeled protein–peptide complexes and analyzing the contacts made by the modeled peptide ligand in the substrate-binding pocket of the MHC or protein kinase. Analysis of these specific inter-molecular contacts is crucial for understanding structural basis of the substrate specificity of these two protein families. This software also provides appropriate interfaces for identifying, putative MHC-binding peptides in the sequence of an antigen or phosphorylation sites on the substrate protein of a kinase, by scoring these inter-molecular contacts using residue-based statistical pair potentials. MODPROPEP would complement various available sequence-based programs (SYFPEITHI, SCANSITE, etc.) for predicting substrates of MHCs and protein kinases. The program is available at http://www.nii.res.in/modpropep.html

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.

In silico grouping of peptide/HLA class I complexes using structural interaction characteristics

MOTIVATION

Classification of human leukocyte antigen (HLA) proteins into supertypes underpins the development of epitope-based vaccines with wide population coverage. Current methods for HLA supertype definition, based on common structural features of HLA proteins and/or their functional binding specificities, leave structural interaction characteristics among different HLA supertypes with antigenic peptides unexplored.

METHODS

We describe the use of structural interaction descriptors for the analysis of 68 peptide/HLA class I crystallographic structures. Interaction parameters computed include the number of intermolecular hydrogen bonds between each HLA protein and its corresponding bound peptide, solvent accessibility, gap volume and gap index.

RESULTS

The structural interactions patterns of peptide/HLA class I complexes investigated herein vary among individual alleles and may be grouped in a supertype dependent manner. Using the proposed methodology, eight HLA class I supertypes were defined based on existing experimental crystallographic structures which largely overlaps (77% consensus) with the definitions by binding motifs. This mode of classification, which considers conformational information of both peptide and HLA proteins, provides an alternative to the characterization of supertypes using either peptide or HLA protein information alone.

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.

Structural prediction of peptides binding to MHC class I molecules

Peptide binding to class I major histocompatibility complex (MHCI) molecules is a key step in the immune response and the structural details of this interaction are of importance in the design of peptide vaccines. Algorithms based on primary sequence have had success in predicting potential antigenic peptides for MHCI, but such algorithms have limited accuracy and provide no structural information. Here, we present an algorithm, PePSSI (peptide-MHC prediction of structure through solvated interfaces), for the prediction of peptide structure when bound to the MHCI molecule, HLA-A2. The algorithm combines sampling of peptide backbone conformations and flexible movement of MHC side chains and is unique among other prediction algorithms in its incorporation of explicit water molecules at the peptide-MHC interface. In an initial test of the algorithm, PePSSI was used to predict the conformation of eight peptides bound to HLA-A2, for which X-ray data are available. Comparison of the predicted and X-ray conformations of these peptides gave RMSD values between 1.301 and 2.475 A. Binding conformations of 266 peptides with known binding affinities for HLA-A2 were then predicted using PePSSI. Structural analyses of these peptide-HLA-A2 conformations showed that peptide binding affinity is positively correlated with the number of peptide-MHC contacts and negatively correlated with the number of interfacial water molecules. These results are consistent with the relatively hydrophobic binding nature of the HLA-A2 peptide binding interface. In summary, PePSSI is capable of rapid and accurate prediction of peptide-MHC binding conformations, which may in turn allow estimation of MHCI-peptide binding affinity.

Computational methods for prediction of T-cell epitopes--a framework for modelling, testing, and applications

Computational models complement laboratory experimentation for efficient identification of MHC-binding peptides and T-cell epitopes. Methods for prediction of MHC-binding peptides include binding motifs, quantitative matrices, artificial neural networks, hidden Markov models, and molecular modelling. Models derived by these methods have been successfully used for prediction of T-cell epitopes in cancer, autoimmunity, infectious disease, and allergy. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures and performed according to strict standards. This requires careful selection of data for model building, and adequate testing and validation. A range of web-based databases and MHC-binding prediction programs are available. Although some available prediction programs for particular MHC alleles have reasonable accuracy, there is no guarantee that all models produce good quality predictions. In this article, we present and discuss a framework for modelling, testing, and applications of computational methods used in predictions of T-cell epitopes.

A novel MHCp binding prediction model

Many statistical and molecular mechanics models have been developed and tested for major histocompatibility complex peptide (MHCp) binding predictions during the last decade. The statistical model prediction using pooled peptide sequence data and three-dimensional modeling prediction by molecular mechanics calculations have been assessed for efficiency and human leukocyte antigen diversity coverage. We describe a novel predictive model using information gleaned from 29 human MHCp crystal structures. The validation for the new model is performed using four different sets of data: (1) MHCp crystal structures, (2) peptides with known IC(50) binding values, (3) peptides tested positive by tetramer staining, (4) peptides with known binding information at the MHCBN database. The model produces high prediction efficiencies (average 60 %) with good sensitivity (approximately 50%-73%) and specificity (52%-58%) values. The average positive predictive value of the model is 89%, while the average negative predictive value is only 18%. The efficiency is very high in predicting binders and very low in predicting nonbinders. This model is superior to many existing methods because of its potential application to any given MHC allele whose sequence is clearly defined.

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.

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.

Efficient discovery of immune response targets by cyclical refinement of QSAR models of peptide binding

Peptides that induce and recall T-cell responses are called T-cell epitopes. T-cell epitopes may be useful in a subunit vaccine against malaria. Computer models that simulate peptide binding to MHC are useful for selecting candidate T-cell epitopes since they minimize the number of experiments required for their identification. We applied a combination of computational and immunological strategies to select candidate T-cell epitopes. A total of 86 experimental binding assays were performed in three rounds of identification of HLA-A11 binding peptides from the six preerythrocytic malaria antigens. Thirty-six peptides were experimentally confirmed as binders. We show that the cyclical refinement of the ANN models results in a significant improvement of the efficiency of identifying potential T-cell epitopes.

Predicting sequences and structures of MHC-binding peptides: a computational combinatorial approach

Peptides bound to MHC molecules on the surface of cells convey critical information about the cellular milieu to immune system T cells. Predicting which peptides can bind an MHC molecule, and understanding their modes of binding, are important in order to design better diagnostic and therapeutic agents for infectious and autoimmune diseases. Due to the difficulty of obtaining sufficient experimental binding data for each human MHC molecule, computational modeling of MHC peptide-binding properties is necessary. This paper describes a computational combinatorial design approach to the prediction of peptides that bind an MHC molecule of known X-ray crystallographic or NMR-determined structure. The procedure uses chemical fragments as models for amino acid residues and produces a set of sequences for peptides predicted to bind in the MHC peptide-binding groove. The probabilities for specific amino acids occurring at each position of the peptide are calculated based on these sequences, and these probabilities show a good agreement with amino acid distributions derived from a MHC-binding peptide database. The method also enables prediction of the three-dimensional structure of MHC-peptide complexes. Docking, linking, and optimization procedures were performed with the XPLOR program [1].

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.

Structure-based prediction of binding peptides to MHC class I molecules: application to a broad range of MHC alleles

Specific binding of antigenic peptides to major histocompatibility complex (MHC) class I molecules is a prerequisite for their recognition by cytotoxic T-cells. Prediction of MHC-binding peptides must therefore be incorporated in any predictive algorithm attempting to identify immunodominant T-cell epitopes, based on the amino acid sequence of the protein antigen. Development of predictive algorithms based on experimental binding data requires experimental testing of a very large number of peptides. A complementary approach relies on the structural conservation observed in crystallographically solved peptide-MHC complexes. By this approach, the peptide structure in the MHC groove is used as a template upon which peptide candidates are threaded, and their compatibility to bind is evaluated by statistical pairwise potentials. Our original algorithm based on this approach used the pairwise potential table of Miyazawa and Jernigan (Miyazawa S, Jernigan RL, 1996, J Mol Biol 256:623-644) and succeeded to correctly identify good binders only for MHC molecules with hydrophobic binding pockets, probably because of the high emphasis of hydrophobic interactions in this table. A recently developed pairwise potential table by Betancourt and Thirumalai (Betancourt MR, Thirumalai D, 1999, Protein Sci 8:361-369) that is based on the Miyazawa and Jernigan table describes the hydrophilic interactions more appropriately. In this paper, we demonstrate how the use of this table, together with a new definition of MHC contact residues by which only residues that contribute exclusively to sequence specific binding are included, allows the development of an improved algorithm that can be applied to a wide range of MHC class I alleles.

Knowledge-based grouping of modeled HLA peptide complexes

Human leukocyte antigens are the most polymorphic of human genes and multiple sequence alignment shows that such polymorphisms are clustered in the functional peptide binding domains. Because of such polymorphism among the peptide binding residues, the prediction of peptides that bind to specific HLA molecules is very difficult. In recent years two different types of computer based prediction methods have been developed and both the methods have their own advantages and disadvantages. The nonavailability of allele specific binding data restricts the use of knowledge-based prediction methods for a wide range of HLA alleles. Alternatively, the modeling scheme appears to be a promising predictive tool for the selection of peptides that bind to specific HLA molecules. The scoring of the modeled HLA-peptide complexes is a major concern. The use of knowledge based rules (van der Waals clashes and solvent exposed hydrophobic residues) to distinguish binders from nonbinders is applied in the present study. The rules based on (1) number of observed atomic clashes between the modeled peptide and the HLA structure, and (2) number of solvent exposed hydrophobic residues on the modeled peptide effectively discriminate experimentally known binders from poor/nonbinders. Solved crystal complexes show no vdW Clash (vdWC) in 95% cases and no solvent exposed hydrophobic peptide residues (SEHPR) were seen in 86% cases. In our attempt to compare experimental binding data with the predicted scores by this scoring scheme, 77% of the peptides are correctly grouped as good binders with a sensitivity of 71%.

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.