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

Immunogenicity of Trypanosoma cruzi Multi-Epitope Recombinant Protein as an Antigen Candidate for Chagas Disease Vaccine in Humans

Chagas disease, caused by the protozoan Trypanosoma cruzi (T. cruzi), is the most significant neglected tropical disease affecting individuals in the Americas. Currently, available drugs, such as nifurtimox and benznidazole (BZN), are both toxic and ineffective in the chronic phase of the disease. A promising alternative is the development of a Chagas disease vaccine, although this effort is hampered by the complexity of the parasite and HLA polymorphisms. In addition, the activation of epitope-specific CD8+ T cells is critical to conferring a robust cell-mediated immune response and protection by producing IFN-γ and perforin. Thus, the antigen (s) for the development of a Chagas vaccine or immunotherapy must include CD8+ T cell epitopes. In this study, we aimed to develop a multi-epitope recombinant protein as a novel human vaccine for Chagas disease. Sixteen database programs were used to predict de novo 40 potential epitopes for the HLA-A*02:01 allele. Nine out of the 40 predicted epitopes were able to elicit IFN-γ production in Peripheral Blood Mononuclear Cells (PBMCs) from Chagas patients. Molecular docking revealed a good binding affinity among the epitopes with diverse HLA molecules. Therefore, a recombinant multi-epitope protein including these nine T. cruzi CD8+ epitopes was expressed and demonstrated to recall an antigen-specific immune response in ex-vivo assays using PBMCs from Chagas patients with the HLA-A*02 allele. These findings support the development of this multi-epitope protein as a promising candidate human vaccine against Chagas disease.

Predicting MHC class I binder: existing approaches and a novel recurrent neural network solution

Major histocompatibility complex (MHC) possesses important research value in the treatment of complex human diseases. A plethora of computational tools has been developed to predict MHC class I binders. Here, we comprehensively reviewed 27 up-to-date MHC I binding prediction tools developed over the last decade, thoroughly evaluating feature representation methods, prediction algorithms and model training strategies on a benchmark dataset from Immune Epitope Database. A common limitation was identified during the review that all existing tools can only handle a fixed peptide sequence length. To overcome this limitation, we developed a bilateral and variable long short-term memory (BVLSTM)-based approach, named BVLSTM-MHC. It is the first variable-length MHC class I binding predictor. In comparison to the 10 mainstream prediction tools on an independent validation dataset, BVLSTM-MHC achieved the best performance in six out of eight evaluated metrics. A web server based on the BVLSTM-MHC model was developed to enable accurate and efficient MHC class I binder prediction in human, mouse, macaque and chimpanzee.

Advanced strategies for development of vaccines against human bacterial pathogens

Infectious diseases are one of the main grounds of death and disabilities in human beings globally. Lack of effective treatment and immunization for many deadly infectious diseases and emerging drug resistance in pathogens underlines the need to either develop new vaccines or sufficiently improve the effectiveness of currently available drugs and vaccines. In this review, we discuss the application of advanced tools like bioinformatics, genomics, proteomics and associated techniques for a rational vaccine design.

The role of immunoinformatics in the development of T-cell peptide-based vaccines against Mycobacterium tuberculosis

INTRODUCTION

Tuberculosis (TB) is a major health problem worldwide. The BCG, the only authorized vaccine to fight TB, shows a variable protection in the adult population highlighting the need of a new vaccine. Immunoinformatics offers a variety of tools that can predict immunogenic T-cell peptides of Mycobacterium tuberculosis (Mtb) that can be used to create a new vaccine. Immunoinformatics has made possible the identification of immunogenic T-cell peptides of Mtb that have been tested in vitro showing a potential for using these molecules as part of a new TB vaccine.

AREAS COVERED

This review summarizes the most common immunoinformatics tools to identify immunogenic T-cell peptides and presents a compilation about research studies that have identified T-cell peptides of Mtb by using immunoinformatics. Also, it is provided a summary of the TB vaccines undergoing clinical trials.

EXPERT OPINION

In the next few years, the field of peptide-based vaccines will keep growing along with the development of more efficient and sophisticated immunoinformatic tools to identify immunogenic peptides with a greater accuracy.

Main Strategies for the Identification of Neoantigens

Genetic instability of tumors leads to the appearance of numerous tumor-specific somatic mutations that could potentially result in the production of mutated peptides that are presented on the cell surface by the MHC molecules. Peptides of this kind are commonly called neoantigens. Their presence on the cell surface specifically distinguishes tumors from healthy tissues. This feature makes neoantigens a promising target for immunotherapy. The rapid evolution of high-throughput genomics and proteomics makes it possible to implement these techniques in clinical practice. In particular, they provide useful tools for the investigation of neoantigens. The most valuable genomic approach to this problem is whole-exome sequencing coupled with RNA-seq. High-throughput mass-spectrometry is another option for direct identification of MHC-bound peptides, which is capable of revealing the entire MHC-bound peptidome. Finally, structure-based predictions could significantly improve the understanding of physicochemical and structural features that affect the immunogenicity of peptides. The development of pipelines combining such tools could improve the accuracy of the peptide selection process and decrease the required time. Here we present a review of the main existing approaches to investigating the neoantigens and suggest a possible ideal pipeline that takes into account all modern trends in the context of neoantigen discovery.

Combining Three-Dimensional Modeling with Artificial Intelligence to Increase Specificity and Precision in Peptide-MHC Binding Predictions

The reliable prediction of the affinity of candidate peptides for the MHC is important for predicting their potential antigenicity and thus influences medical applications, such as decisions on their inclusion in T cell-based vaccines. In this study, we present a rapid, predictive computational approach that combines a popular, sequence-based artificial neural network method, NetMHCpan 4.0, with three-dimensional structural modeling. We find that the ensembles of bound peptide conformations generated by the programs MODELLER and Rosetta FlexPepDock are less variable in geometry for strong binders than for low-affinity peptides. In tests on 1271 peptide sequences for which the experimental dissociation constants of binding to the well-characterized murine MHC allele H-2Db are known, by applying thresholds for geometric fluctuations the structure-based approach in a standalone manner drastically improves the statistical specificity, reducing the number of false positives. Furthermore, filtering candidates generated with NetMHCpan 4.0 with the structure-based predictor led to an increase in the positive predictive value (PPV) of the peptides correctly predicted to bind very strongly (i.e., K d < 100 nM) from 40 to 52% (p = 0.027). The combined method also significantly improved the PPV when tested on five human alleles, including some with limited data for training. Overall, an average increase of 10% in the PPV was found over the standalone sequence-based method. The combined method should be useful in the rapid design of effective T cell-based vaccines.

Computational investigation of peptide binding stabilities of HLA-B*27 and HLA-B*44 alleles

Major Histocompatibility Complex (MHC) is a cell surface glycoprotein that binds to foreign antigens and presents them to T lymphocyte cells on the surface of Antigen Presenting Cells (APCs) for appropriate immune recognition. Recently, studies focusing on peptide-based vaccine design have allowed a better understanding of peptide immunogenicity mechanisms, which is defined as the ability of a peptide to stimulate CTL-mediated immune response. Peptide immunogenicity is also known to be related to the stability of peptide-loaded MHC (pMHC) complex. In this study, ENCoM server was used for structure-based estimation of the impact of single point mutations on pMHC complex stabilities. For this purpose, two human MHC molecules from the HLA-B*27 group (HLA-B*27:05 and HLA-B*27:09) in complex with four different peptides (GRFAAAIAK, RRKWRRWHL, RRRWRRLTV and IRAAPPPLF) and three HLA-B*44 molecules (HLA-B*44:02, HLA-B*44:03 and HLA-B*44:05) in complex with two different peptides (EEYLQAFTY and EEYLKAWTF) were analyzed. We found that the stability of pMHC complexes is dependent on both peptide sequence and MHC allele. Furthermore, we demonstrate that allele-specific peptide-binding preferences can be accurately revealed using structure-based computational methods predicting the effect of mutations on protein stability.

Postgenomic Approaches and Bioinformatics Tools to Advance the Development of Vaccines against Bacteria of the Burkholderia cepacia Complex

Bacteria of the Burkholderia cepacia complex (Bcc) remain an important cause of morbidity and mortality among patients suffering from cystic fibrosis. Eradication of these pathogens by antimicrobial therapy often fails, highlighting the need to develop novel strategies to eradicate infections. Vaccines are attractive since they can confer protection to particularly vulnerable patients, as is the case of cystic fibrosis patients. Several studies have identified specific virulence factors and proteins as potential subunit vaccine candidates. So far, no vaccine is available to protect from Bcc infections. In the present work, we review the most promising postgenomic approaches and selected web tools available to speed up the identification of immunogenic proteins with the potential of conferring protection against Bcc infections.

PSSMHCpan: a novel PSSM-based software for predicting class I peptide-HLA binding affinity

Predicting peptide binding affinity with human leukocyte antigen (HLA) is a crucial step in developing powerful antitumor vaccine for cancer immunotherapy. Currently available methods work quite well in predicting peptide binding affinity with HLA alleles such as HLA-A*0201, HLA-A*0101, and HLA-B*0702 in terms of sensitivity and specificity. However, quite a few types of HLA alleles that are present in the majority of human populations including HLA-A*0202, HLA-A*0203, HLA-A*6802, HLA-B*5101, HLA-B*5301, HLA-B*5401, and HLA-B*5701 still cannot be predicted with satisfactory accuracy using currently available methods. Furthermore, currently the most popularly used methods for predicting peptide binding affinity are inefficient in identifying neoantigens from a large quantity of whole genome and transcriptome sequencing data. Here we present a Position Specific Scoring Matrix (PSSM)-based software called PSSMHCpan to accurately and efficiently predict peptide binding affinity with a broad coverage of HLA class I alleles. We evaluated the performance of PSSMHCpan by analyzing 10-fold cross-validation on a training database containing 87 HLA alleles and obtained an average area under receiver operating characteristic curve (AUC) of 0.94 and accuracy (ACC) of 0.85. In an independent dataset (Peptide Database of Cancer Immunity) evaluation, PSSMHCpan is substantially better than the popularly used NetMHC-4.0, NetMHCpan-3.0, PickPocket, Nebula, and SMM with a sensitivity of 0.90, as compared to 0.74, 0.81, 0.77, 0.24, and 0.79. In addition, PSSMHCpan is more than 197 times faster than NetMHC-4.0, NetMHCpan-3.0, PickPocket, sNebula, and SMM when predicting neoantigens from 661 263 peptides from a breast tumor sample. Finally, we built a neoantigen prediction pipeline and identified 117 017 neoantigens from 467 cancer samples of various cancers from TCGA. PSSMHCpan is superior to the currently available methods in predicting peptide binding affinity with a broad coverage of HLA class I alleles.

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.

Ab-initio conformational epitope structure prediction using genetic algorithm and SVM for vaccine design

BACKGROUND AND OBJECTIVE

T-cell epitope structure identification is a significant challenging immunoinformatic problem within epitope-based vaccine design. Epitopes or antigenic peptides are a set of amino acids that bind with the Major Histocompatibility Complex (MHC) molecules. The aim of this process is presented by Antigen Presenting Cells to be inspected by T-cells. MHC-molecule-binding epitopes are responsible for triggering the immune response to antigens. The epitope's three-dimensional (3D) molecular structure (i.e., tertiary structure) reflects its proper function. Therefore, the identification of MHC class-II epitopes structure is a significant step towards epitope-based vaccine design and understanding of the immune system.

METHODS

In this paper, we propose a new technique using a Genetic Algorithm for Predicting the Epitope Structure (GAPES), to predict the structure of MHC class-II epitopes based on their sequence. The proposed Elitist-based genetic algorithm for predicting the epitope's tertiary structure is based on Ab-Initio Empirical Conformational Energy Program for Peptides (ECEPP) Force Field Model. The developed secondary structure prediction technique relies on Ramachandran Plot. We used two alignment algorithms: the ROSS alignment and TM-Score alignment. We applied four different alignment approaches to calculate the similarity scores of the dataset under test. We utilized the support vector machine (SVM) classifier as an evaluation of the prediction performance.

RESULTS

The prediction accuracy and the Area Under Receiver Operating Characteristic (ROC) Curve (AUC) were calculated as measures of performance. The calculations are performed on twelve similarity-reduced datasets of the Immune Epitope Data Base (IEDB) and a large dataset of peptide-binding affinities to HLA-DRB1*0101. The results showed that GAPES was reliable and very accurate. We achieved an average prediction accuracy of 93.50% and an average AUC of 0.974 in the IEDB dataset. Also, we achieved an accuracy of 95.125% and an AUC of 0.987 on the HLA-DRB1*0101 allele of the Wang benchmark dataset.

CONCLUSIONS

The results indicate that the proposed prediction technique "GAPES" is a promising technique that will help researchers and scientists to predict the protein structure and it will assist them in the intelligent design of new epitope-based vaccines.

DynaDom: structure-based prediction of T cell receptor inter-domain and T cell receptor-peptide-MHC (class I) association angles

Background

T cell receptor (TCR) molecules are involved in the adaptive immune response as they distinguish between self- and foreign-peptides, presented in major histocompatibility complex molecules (pMHC). Former studies showed that the association angles of the TCR variable domains (Vα/Vβ) can differ significantly and change upon binding to the pMHC complex. These changes can be described as a rotation of the domains around a general Center of Rotation, characterized by the interaction of two highly conserved glutamine residues.

Methods

We developed a computational method, DynaDom, for the prediction of TCR Vα/Vβ inter-domain and TCR/pMHC orientations in TCRpMHC complexes, which allows predicting the orientation of multiple protein-domains. In addition, we implemented a new approach to predict the correct orientation of the carboxamide endgroups in glutamine and asparagine residues, which can also be used as an external, independent tool.

Results

The approach was evaluated for the remodeling of 75 and 53 experimental structures of TCR and TCRpMHC (class I) complexes, respectively. We show that the DynaDom method predicts the correct orientation of the TCR Vα/Vβ angles in 96 and 89% of the cases, for the poses with the best RMSD and best interaction energy, respectively. For the concurrent prediction of the TCR Vα/Vβ and pMHC orientations, the respective rates reached 74 and 72%. Through an exhaustive analysis, we could show that the pMHC placement can be further improved by a straightforward, yet very time intensive extension of the current approach.

Conclusions

The results obtained in the present remodeling study prove the suitability of our approach for interdomain-angle optimization. In addition, the high prediction rate obtained specifically for the energetically highest ranked poses further demonstrates that our method is a powerful candidate for blind prediction. Therefore it should be well suited as part of any accurate atomistic modeling pipeline for TCRpMHC complexes and potentially other large molecular assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-016-0071-7) contains supplementary material, which is available to authorized users.

An Overview of Bioinformatics Tools and Resources in Allergy

The rapidly increasing number of characterized allergens has created huge demands for advanced information storage, retrieval, and analysis. Bioinformatics and machine learning approaches provide useful tools for the study of allergens and epitopes prediction, which greatly complement traditional laboratory techniques. The specific applications mainly include identification of B- and T-cell epitopes, and assessment of allergenicity and cross-reactivity. In order to facilitate the work of clinical and basic researchers who are not familiar with bioinformatics, we review in this chapter the most important databases, bioinformatic tools, and methods with relevance to the study of allergens.

Prediction of peptide binding to a major histocompatibility complex class I molecule based on docking simulation

Binding between major histocompatibility complex (MHC) class I molecules and immunogenic epitopes is one of the most important processes for cell-mediated immunity. Consequently, computational prediction of amino acid sequences of MHC class I binding peptides from a given sequence may lead to important biomedical advances. In this study, an efficient structure-based method for predicting peptide binding to MHC class I molecules was developed, in which the binding free energy of the peptide was evaluated by two individual docking simulations. An original penalty function and restriction of degrees of freedom were determined by analysis of 361 published X-ray structures of the complex and were then introduced into the docking simulations. To validate the method, calculations using a 50-amino acid sequence as a prediction target were performed. In 27 calculations, the binding free energy of the known peptide was within the top 5 of 166 peptides generated from the 50-amino acid sequence. Finally, demonstrative calculations using a whole sequence of a protein as a prediction target were performed. These data clearly demonstrate high potential of this method for predicting peptide binding to MHC class I molecules.

modPDZpep: a web resource for structure based analysis of human PDZ-mediated interaction networks

Background

PDZ domains recognize short sequence stretches usually present in C-terminal of their interaction partners. Because of the involvement of PDZ domains in many important biological processes, several attempts have been made for developing bioinformatics tools for genome-wide identification of PDZ interaction networks. Currently available tools for prediction of interaction partners of PDZ domains utilize machine learning approach. Since, they have been trained using experimental substrate specificity data for specific PDZ families, their applicability is limited to PDZ families closely related to the training set. These tools also do not allow analysis of PDZ-peptide interaction interfaces.

Results

We have used a structure based approach to develop modPDZpep, a program to predict the interaction partners of human PDZ domains and analyze structural details of PDZ interaction interfaces. modPDZpep predicts interaction partners by using structural models of PDZ-peptide complexes and evaluating binding energy scores using residue based statistical pair potentials. Since, it does not require training using experimental data on peptide binding affinity, it can predict substrates for diverse PDZ families. Because of the use of simple scoring function for binding energy, it is also fast enough for genome scale structure based analysis of PDZ interaction networks. Benchmarking using artificial as well as real negative datasets indicates good predictive power with ROC-AUC values in the range of 0.7 to 0.9 for a large number of human PDZ domains. Another novel feature of modPDZpep is its ability to map novel PDZ mediated interactions in human protein-protein interaction networks, either by utilizing available experimental phage display data or by structure based predictions.

Conclusions

In summary, we have developed modPDZpep, a web-server for structure based analysis of human PDZ domains. It is freely available at http://www.nii.ac.in/modPDZpep.html or http://202.54.226.235/modPDZpep.html.

Reviewers

This article was reviewed by Michael Gromiha and Zoltán Gáspári.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0151-4) contains supplementary material, which is available to authorized users.

Chemical Modification of Influenza CD8+ T-Cell Epitopes Enhances Their Immunogenicity Regardless of Immunodominance

T cells are essential players in the defense against infection. By targeting the MHC class I antigen-presenting pathway with peptide-based vaccines, antigen-specific T cells can be induced. However, low immunogenicity of peptides poses a challenge. Here, we set out to increase immunogenicity of influenza-specific CD8+ T cell epitopes. By substituting amino acids in wild type sequences with non-proteogenic amino acids, affinity for MHC can be increased, which may ultimately enhance cytotoxic CD8+ T cell responses. Since preventive vaccines against viruses should induce a broad immune response, we used this method to optimize influenza-specific epitopes of varying dominance. For this purpose, HLA-A*0201 epitopes GILGFVFTL, FMYSDFHFI and NMLSTVLGV were selected in order of decreasing MHC-affinity and dominance. For all epitopes, we designed chemically enhanced altered peptide ligands (CPLs) that exhibited greater binding affinity than their WT counterparts; even binding scores of the high affinity GILGFVFTL epitope could be improved. When HLA-A*0201 transgenic mice were vaccinated with selected CPLs, at least 2 out of 4 CPLs of each epitope showed an increase in IFN-γ responses of splenocytes. Moreover, modification of the low affinity epitope NMLSTVLGV led to an increase in the number of mice that responded. By optimizing three additional influenza epitopes specific for HLA-A*0301, we show that this strategy can be extended to other alleles. Thus, enhancing binding affinity of peptides provides a valuable tool to improve the immunogenicity and range of preventive T cell-targeted peptide vaccines.

HLaffy: estimating peptide affinities for Class-1 HLA molecules by learning position-specific pair potentials

MOTIVATION

T-cell epitopes serve as molecular keys to initiate adaptive immune responses. Identification of T-cell epitopes is also a key step in rational vaccine design. Most available methods are driven by informatics and are critically dependent on experimentally obtained training data. Analysis of a training set from Immune Epitope Database (IEDB) for several alleles indicates that the sampling of the peptide space is extremely sparse covering a tiny fraction of the possible nonamer space, and also heavily skewed, thus restricting the range of epitope prediction.

RESULTS

We present a new epitope prediction method that has four distinct computational modules: (i) structural modelling, estimating statistical pair-potentials and constraint derivation, (ii) implicit modelling and interaction profiling, (iii) feature representation and binding affinity prediction and (iv) use of graphical models to extract peptide sequence signatures to predict epitopes for HLA class I alleles.

CONCLUSIONS

HLaffy is a novel and efficient epitope prediction method that predicts epitopes for any Class-1 HLA allele, by estimating the binding strengths of peptide-HLA complexes which is achieved through learning pair-potentials important for peptide binding. It relies on the strength of the mechanistic understanding of peptide-HLA recognition and provides an estimate of the total ligand space for each allele. The performance of HLaffy is seen to be superior to the currently available methods.

AVAILABILITY AND IMPLEMENTATION

The method is made accessible through a webserver http://proline.biochem.iisc.ernet.in/HLaffy

CONTACT

: nchandra@biochem.iisc.ernet.in

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Third-Kind Encounters in Biomedicine: Immunology Meets Mathematics and Informatics to Become Quantitative and Predictive

The understanding of the immune response is right now at the center of biomedical research. There are growing expectations that immune-based interventions will in the midterm provide new, personalized, and targeted therapeutic options for many severe and highly prevalent diseases, from aggressive cancers to infectious and autoimmune diseases. To this end, immunology should surpass its current descriptive and phenomenological nature, and become quantitative, and thereby predictive.Immunology is an ideal field for deploying the tools, methodologies, and philosophy of systems biology, an approach that combines quantitative experimental data, computational biology, and mathematical modeling. This is because, from an organism-wide perspective, the immunity is a biological system of systems, a paradigmatic instance of a multi-scale system. At the molecular scale, the critical phenotypic responses of immune cells are governed by large biochemical networks, enriched in nested regulatory motifs such as feedback and feedforward loops. This network complexity confers them the ability of highly nonlinear behavior, including remarkable examples of homeostasis, ultra-sensitivity, hysteresis, and bistability. Moving from the cellular level, different immune cell populations communicate with each other by direct physical contact or receiving and secreting signaling molecules such as cytokines. Moreover, the interaction of the immune system with its potential targets (e.g., pathogens or tumor cells) is far from simple, as it involves a number of attack and counterattack mechanisms that ultimately constitute a tightly regulated multi-feedback loop system. From a more practical perspective, this leads to the consequence that today's immunologists are facing an ever-increasing challenge of integrating massive quantities from multi-platforms.In this chapter, we support the idea that the analysis of the immune system demands the use of systems-level approaches to ensure the success in the search for more effective and personalized immune-based therapies.

Major histocompatibility complex linked databases and prediction tools for designing vaccines

Presently, the major histocompatibility complex (MHC) is receiving considerable interest owing to its remarkable role in antigen presentation and vaccine design. The specific databases and prediction approaches related to MHC sequences, structures and binding/nonbinding peptides have been aggressively developed in the past two decades with their own benchmarks and standards. Before using these databases and prediction tools, it is important to analyze why and how the tools are constructed along with their strengths and limitations. The current review presents insights into web-based immunological bioinformatics resources that include searchable databases of MHC sequences, epitopes and prediction tools that are linked to MHC based vaccine design, including population coverage analysis. In T cell epitope forecasts, MHC class I binding predictions are very accurate for most of the identified MHC alleles. However, these predictions could be further improved by integrating proteasome cleavage (in conjugation with transporter associated with antigen processing (TAP) binding) prediction, as well as T cell receptor binding prediction. On the other hand, MHC class II restricted epitope predictions display relatively low accuracy compared to MHC class I. To date, pan-specific tools have been developed, which not only deliver significantly improved predictions in terms of accuracy, but also in terms of the coverage of MHC alleles and supertypes. In addition, structural modeling and simulation systems for peptide-MHC complexes enable the molecular-level investigation of immune processes. Finally, epitope prediction tools, and their assessments and guidelines, have been presented to immunologist for the design of novel vaccine and diagnostics.

Predicting Hemagglutinin MHC-II Ligand Analogues in Anti-TNFα Biologics: Implications for Immunogenicity of Pharmaceutical Proteins

The purpose of this study was to evaluate the extent of overlapping immunogenic peptides between three pharmaceutical biologics and influenza viruses. Clinical studies have shown that subsets of patients with rheumatoid arthritis (RA) develop anti-drug antibodies towards anti-TNFα biologics. We postulate that common infectious pathogens, including influenza viruses, may sensitize RA patients toward recombinant proteins. We hypothesize that embedded within infliximab (IFX), adalimumab (ADA), and etanercept (ETN) are ligands of class II major histocompatibility complex (MHC-II) that mimic T cell epitopes derived from influenza hemagglutinin (HA). The rationale is that repeated administration of the biologics would reactivate HA-primed CD4 T cells, stimulating B cells to produce cross-reactive antibodies. Custom scripts were constructed using MATLAB to compare MHC-II ligands of HA and the biologics; all ligands were predicted using tools in Immune Epitope Database and Resources (IEDB). We analyzed three HLA-DR1 alleles (0101, 0401 and 1001) that are prominent in RA patients, and two alleles (0103 and 1502) that are not associated with RA. The results indicate that 0401 would present more analogues of HA ligands in the three anti-TNFα biologics compared to the other alleles. The approach led to identification of potential ligands in IFX and ADA that shares sequence homology with a known HA-specific CD4 T cell epitope. We also discovered a peptide in the complementarity-determining region 3 (CDR-3) of ADA that encompasses both a potential CD4 T cell epitope and a known B cell epitope in HA. The results may help generate new hypotheses for interrogating patient variability of immunogenicity of the anti-TNFα drugs. The approach would aid development of new recombinant biologics by identifying analogues of CD4 T cell epitopes of common pathogens at the preclinical stage.

Analysis of HLA A*02 association with vaccine efficacy in the RV144 HIV-1 vaccine trial

The RV144 HIV-1 vaccine trial demonstrated partial efficacy of 31% against HIV-1 infection. Studies into possible correlates of protection found that antibodies specific to the V1 and V2 (V1/V2) region of envelope correlated inversely with infection risk and that viruses isolated from trial participants contained genetic signatures of vaccine-induced pressure in the V1/V2 region. We explored the hypothesis that the genetic signatures in V1 and V2 could be partly attributed to selection by vaccine-primed T cells. We performed a T-cell-based sieve analysis of breakthrough viruses in the RV144 trial and found evidence of predicted HLA binding escape that was greater in vaccine versus placebo recipients. The predicted escape depended on class I HLA A*02- and A*11-restricted epitopes in the MN strain rgp120 vaccine immunogen. Though we hypothesized that this was indicative of postacquisition selection pressure, we also found that vaccine efficacy (VE) was greater in A*02-positive (A*02(+)) participants than in A*02(-) participants (VE = 54% versus 3%, P = 0.05). Vaccine efficacy against viruses with a lysine residue at site 169, important to antibody binding and implicated in vaccine-induced immune pressure, was also greater in A*02(+) participants (VE = 74% versus 15%, P = 0.02). Additionally, a reanalysis of vaccine-induced immune responses that focused on those that were shown to correlate with infection risk suggested that the humoral responses may have differed in A*02(+) participants. These exploratory and hypothesis-generating analyses indicate there may be an association between a class I HLA allele and vaccine efficacy, highlighting the importance of considering HLA alleles and host immune genetics in HIV vaccine trials.

IMPORTANCE

The RV144 trial was the first to show efficacy against HIV-1 infection. Subsequently, much effort has been directed toward understanding the mechanisms of protection. Here, we conducted a T-cell-based sieve analysis, which compared the genetic sequences of viruses isolated from infected vaccine and placebo recipients. Though we hypothesized that the observed sieve effect indicated postacquisition T-cell selection, we also found that vaccine efficacy was greater for participants who expressed HLA A*02, an allele implicated in the sieve analysis. Though HLA alleles have been associated with disease progression and viral load in HIV-1 infection, these data are the first to suggest the association of a class I HLA allele and vaccine efficacy. While these statistical analyses do not provide mechanistic evidence of protection in RV144, they generate testable hypotheses for the HIV vaccine community and they highlight the importance of assessing the impact of host immune genetics in vaccine-induced immunity and protection. (This study has been registered at ClinicalTrials.gov under registration no. NCT00223080.).

Immunoinformatics: a brief review

A large volume of data relevant to immunology research has accumulated due to sequencing of genomes of the human and other model organisms. At the same time, huge amounts of clinical and epidemiologic data are being deposited in various scientific literature and clinical records. This accumulation of the information is like a goldmine for researchers looking for mechanisms of immune function and disease pathogenesis. Thus the need to handle this rapidly growing immunological resource has given rise to the field known as immunoinformatics. Immunoinformatics, otherwise known as computational immunology, is the interface between computer science and experimental immunology. It represents the use of computational methods and resources for the understanding of immunological information. It not only helps in dealing with huge amount of data but also plays a great role in defining new hypotheses related to immune responses. This chapter reviews classical immunology, different databases, and prediction tool. Further, it briefly describes applications of immunoinformatics in reverse vaccinology, immune system modeling, and cancer diagnosis and therapy. It also explores the idea of integrating immunoinformatics with systems biology for the development of personalized medicine. All these efforts save time and cost to a great extent.

Thermodynamics of Peptide-MHC Class II Interactions: Not all Complexes are Created Equal

The adaptive immune response begins when CD4+ T cells recognize antigenic peptides bound to class II molecules of the Major Histocompatibility Complex (MHCII). The interaction between peptides and MHCII has been historically interpreted as a rigid docking event. However, this model has been challenged by the evidence that conformational flexibility plays an important role in peptide-MHCII complex formation. Thermodynamic analysis of the binding reaction suggests a model of complexation in which the physical-chemical nature of the peptide determines the variability in flexibility of the substates in the peptide-MHC conformational ensemble. This review discusses our understanding of the correlation between thermodynamics of peptide binding and structural features of the resulting complex as well as their impact on HLA-DM activity and on our ability to predict MHCII-restricted epitopes.

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.

Protein-peptide complex prediction through fragment interaction patterns

The number of protein-peptide interactions in a cell is so large that experimental determination of all these complex structures would be a daunting task. Although homology modeling and refinement protocols have vastly improved the number and quality of predicted structural models, ab initio methods are still challenged by both the large number of possible docking sites and the conformational space accessible to flexible peptides. We present a method that addresses these challenges by sampling the entire accessible surface of a protein with a reduced conformational space of interacting backbone fragment pairs from unrelated structures. We demonstrate its potential by predicting ab initio the bound structure for a variety of protein-peptide complexes. In addition, we show the potential of our method for the discovery of domain interaction sites and domain-domain docking.

Using databases and data mining in vaccinology

Throughout time functional immunology has accumulated vast amounts of quantitative and qualitative data relevant to the design and discovery of vaccines. Such data includes, but is not limited to, components of the host and pathogen genome (including antigens and virulence factors), T- and B-cell epitopes and other components of the antigen presentation pathway and allergens. In this review the authors discuss a range of databases that archive such data. Built on such information, increasingly sophisticated data mining techniques have developed that create predictive models of utilitarian value. With special reference to epitope data, the authors discuss the strengths and weaknesses of the available techniques and how they can aid computer-aided vaccine design deliver added value for vaccinology.

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 prediction of Major Histocompatibility Complex (MHC) epitopes

Because of the enormous diversity of both MHC proteins and peptide epitopes, computational epitope prediction methods are needed in order to supplement limited experimental data. These prediction methods are useful for guiding experiments and have many potential biomedical applications. Unlike popular sequence-based methods, structure-based epitope prediction methods can predict epitopes for multiple MHC types with highly distinct peptide binding propensities. In this chapter, we describe in detail our previously developed structure-based epitope prediction methods for both class I and class II MHC proteins. We also discuss the relative advantages and disadvantages of sequence-based versus structure-based methods and how to evaluate prediction performance.

EpicCapo: epitope prediction using combined information of amino acid pairwise contact potentials and HLA-peptide contact site information

Background

Epitope identification is an essential step toward synthetic vaccine development since epitopes play an important role in activating immune response. Classical experimental approaches are laborious and time-consuming, and therefore computational methods for generating epitope candidates have been actively studied. Most of these methods, however, are based on sophisticated nonlinear techniques for achieving higher predictive performance. The use of these techniques tend to diminish their interpretability with respect to binding potential: that is, they do not provide much insight into binding mechanisms.

Results

We have developed a novel epitope prediction method named EpicCapo and its variants, EpicCapo⁺ and EpicCapo^+REF. Nonapeptides were encoded numerically using a novel peptide-encoding scheme for machine learning algorithms by utilizing 40 amino acid pairwise contact potentials (referred to as AAPPs throughout this paper). The predictive performances of EpicCapo⁺ and EpicCapo^+REF outperformed other state-of-the-art methods without losing interpretability. Interestingly, the most informative AAPPs estimated by our study were those developed by Micheletti and Simons while previous studies utilized two AAPPs developed by Miyazawa & Jernigan and Betancourt & Thirumalai. In addition, we found that all amino acid positions in nonapeptides could effect on performances of the predictive models including non-anchor positions. Finally, EpicCapo^+REF was applied to identify candidates of promiscuous epitopes. As a result, 67.1% of the predicted nonapeptides epitopes were consistent with preceding studies based on immunological experiments.

Conclusions

Our method achieved high performance in testing with benchmark datasets. In addition, our study identified a number of candidates of promiscuous CTL epitopes consistent with previously reported immunological experiments. We speculate that our techniques may be useful in the development of new vaccines. The R implementation of EpicCapo^+REF is available at http://pirun.ku.ac.th/~fsciiok/EpicCapoREF.zip. Datasets are available at http://pirun.ku.ac.th/~fsciiok/Datasets.zip.

PAAQD: Predicting immunogenicity of MHC class I binding peptides using amino acid pairwise contact potentials and quantum topological molecular similarity descriptors

Prediction of peptide immunogenicity is a promising approach for novel vaccine discovery. Conventionally, epitope prediction methods have been developed to accelerate the process of vaccine production by searching for candidate peptides from pathogenic proteins. However, recent studies revealed that peptides with high binding affinity to major histocompatibility complex molecules (MHCs) do not always result in high immunogenicity. Therefore, it is promising to predict the peptide immunogenicity rather than epitopes in order to discover new vaccines more effectively. To this end, we developed a novel T-cell reactivity predictor which we call PAAQD. Nonapeptides were encoded numerically, using combining information of amino acid pairwise contact potentials (AAPPs) and quantum topological molecular similarity (QTMS) descriptors. Encoded data were used in the construction of our classification model. Our numerical experiments suggested that the predictive performance of PAAQD is at least comparable with POPISK, one of the pioneering techniques for T-cell reactivity prediction. Also, our experiment suggested that the first and eighth positions of nonapeptides are the most important for immunogenicity and most of the anchor residues in epitope prediction were not important in T-cell reactivity prediction. The R implementation of PAAQD is available at http://pirun.ku.ac.th/~fsciiok/PAAQD.rar.

NEURAL MODELS FOR PREDICTING VIRAL VACCINE TARGETS

We applied artificial neural networks (ANN) for the prediction of targets of immune responses that are useful for study of vaccine formulations against viral infections. Using a novel data representation, we developed a system termed MULTIPRED that can predict peptide binding to multiple related human leukocyte antigens (HLA). This implementation showed high accuracy in the prediction of the promiscuous peptides that bind to five HLA-A2 allelic variants. MULTIPRED is useful for the identification of peptides that bind multiple HLA-A2 variants as a group. By implementing ANN as a classification engine, we enabled both the prediction of peptides binding to multiple individual HLA-A2 molecules and the prediction of promiscuous binders using a single model. The ANN MULTIPRED predicts peptide binding to HLA-A*0205 with excellent accuracy (area under the receiver operating characteristic curve — AROC > 0.90), and to HLA-A*0201, HLA-A*0204 and HLA-A*0206 with high accuracy (AROC > 0.85). Antigenic regions with high density of binders ("antigenic hot-spots") represent best targets for vaccine design. MULTIPRED not only predicts individual 9-mer binders but also predicts antigenic hot spots. Two HLA-A2 hot-spots in Severe Acute Respiratory Syndrome Coronavirus (SARS–CoV) membrane protein were predicted by using MULTIPRED.

PeptX: using genetic algorithms to optimize peptides for MHC binding

Background

The binding between the major histocompatibility complex and the presented peptide is an indispensable prerequisite for the adaptive immune response. There is a plethora of different in silico techniques for the prediction of the peptide binding affinity to major histocompatibility complexes. Most studies screen a set of peptides for promising candidates to predict possible T cell epitopes. In this study we ask the question vice versa: Which peptides do have highest binding affinities to a given major histocompatibility complex according to certain in silico scoring functions?

Results

Since a full screening of all possible peptides is not feasible in reasonable runtime, we introduce a heuristic approach. We developed a framework for Genetic Algorithms to optimize peptides for the binding to major histocompatibility complexes. In an extensive benchmark we tested various operator combinations. We found that (1) selection operators have a strong influence on the convergence of the population while recombination operators have minor influence and (2) that five different binding prediction methods lead to five different sets of "optimal" peptides for the same major histocompatibility complex. The consensus peptides were experimentally verified as high affinity binders.

Conclusion

We provide a generalized framework to calculate sets of high affinity binders based on different previously published scoring functions in reasonable runtime. Furthermore we give insight into the different behaviours of operators and scoring functions of the Genetic Algorithm.

Emerging vaccine informatics

Vaccine informatics is an emerging research area that focuses on development and applications of bioinformatics methods that can be used to facilitate every aspect of the preclinical, clinical, and postlicensure vaccine enterprises. Many immunoinformatics algorithms and resources have been developed to predict T- and B-cell immune epitopes for epitope vaccine development and protective immunity analysis. Vaccine protein candidates are predictable in silico from genome sequences using reverse vaccinology. Systematic transcriptomics and proteomics gene expression analyses facilitate rational vaccine design and identification of gene responses that are correlates of protection in vivo. Mathematical simulations have been used to model host-pathogen interactions and improve vaccine production and vaccination protocols. Computational methods have also been used for development of immunization registries or immunization information systems, assessment of vaccine safety and efficacy, and immunization modeling. Computational literature mining and databases effectively process, mine, and store large amounts of vaccine literature and data. Vaccine Ontology (VO) has been initiated to integrate various vaccine data and support automated reasoning.

Large-scale characterization of peptide-MHC binding landscapes with structural simulations

Class I major histocompatibility complex proteins play a critical role in the adaptive immune system by binding to peptides derived from cytosolic proteins and presenting them on the cell surface for surveillance by T cells. The varied peptide binding specificity of these highly polymorphic molecules has important consequences for vaccine design, transplantation, autoimmunity, and cancer development. Here, we describe a molecular modeling study of MHC-peptide interactions that integrates sampling techniques from protein-protein docking, loop modeling, de novo structure prediction, and protein design in order to construct atomically detailed peptide binding landscapes for a diverse set of MHC proteins. Specificity profiles derived from these landscapes recover key features of experimental binding profiles and can be used to predict peptide binding with reasonable accuracy. Family wide comparison of the predicted binding landscapes recapitulates previously reported patterns of specificity divergence and peptide-repertoire diversity while providing a structural basis for observed specificity patterns. The size and sequence diversity of these structure-based binding landscapes enable us to identify subtle patterns of covariation between peptide sequence positions; analysis of the associated structural models suggests physical interactions that may mediate these sequence correlations.

Predicting peptide binding to Major Histocompatibility Complex molecules

The Major Histocompatibility Complex (MHC) constitutes an important part of the human immune system. During infection, pathogenic proteins are processed into peptide fragments by the antigen processing machinery. These peptides bind to MHC molecules and the MHC-peptide complex is then transported to the cell membrane where it elicits an immune response via T-cell binding. Understanding the molecular mechanism of this process will greatly assist in determining the aetiology of various diseases and in the design of effective drugs. One of the most challenging aspects of this area of research is understanding the specificity and sensitivity of the binding process. An empirical approach to the problem is unfeasible as there are over 512 billion potential binding peptides for each MHC molecule. Computational approaches offer the promise of predicting peptide binding, thus dramatically reducing the number of peptides proceeding to experimental verification. Various bioinformatic approaches have been developed to predict whether or not a particular peptide will bind to a particular MHC allele. Currently, peptide binding prediction methods can be categorised into three major groups: motif- and scoring matrix-based methods, artificial intelligence- (AI-) based methods, and structure-based methods. The first two are sequence-based approaches and are generally based on common sequence motifs in peptides known to bind to MHC molecules. The structure-based approach concerns the structural features and the distribution of energy between the binding peptide and the MHC molecule. Although knowledge of the molecular structure of the MHC molecules is expected to lead to better predictions of peptide binding, the development of structure-based methods has been relatively slow compared to sequence-based methods. Comparisons of various methods showed that the best sequence-based methods significantly outperform structure-based methods. This may be improved by producing more structures and binding data desperately needed by many alleles, especially class II molecules. On the other hand, the large number of verification methods and indicators used by structure-based studies hinders critical evaluation of the methods. Adopting commonly used assessment procedures can demonstrate the relative performance of structure-based methods in a straightforward comparison with other methods. This review provides an overview of current methods for predicting peptide binding to the MHC, with a focus on structure-based methods, and explores the potential for future development in this area.

Comparative analysis of epitope predictions: proposed library of putative vaccine candidates for HIV

Designing a vaccine for a disease is one of the crucial tasks that involve millions and billions of dollars, several decades and yet there is no guarantee of successful results. Several pharmaceutical companies are investing their money and time in such activities. Computational biology could be of great help in these activities by proving a library of plausible candidates that might actually show some positive responses. MHC binding peptide prediction is one such area where the immense power of computers could be used to get a breakthrough. In this direction several databases and servers have been developed by many labs to predict the MHC binding peptides. These short peptides on the antigen surface are recognized by the MHC molecule and are presented to the receptors of T-cells for further immune response. Peptides that bind to a given MHC molecule share sequence similarity. Here we present a comparative study of servers that can predict the MHC binding peptides in a given protein sequence of the antigen. Based on this comparative analysis on HIV data, we are able to propose a library of putative vaccine candidates for the env GP-160 protein of HIV-1.

T-cell epitope prediction and immune complex simulation using molecular dynamics: state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

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.

Immunoinformatics: an integrated scenario

Genome sequencing of humans and other organisms has led to the accumulation of huge amounts of data, which include immunologically relevant data. A large volume of clinical data has been deposited in several immunological databases and as a result immunoinformatics has emerged as an important field which acts as an intersection between experimental immunology and computational approaches. It not only helps in dealing with the huge amount of data but also plays a role in defining new hypotheses related to immune responses. This article reviews classical immunology, different databases and prediction tools. It also describes applications of immunoinformatics in designing in silico vaccination and immune system modelling. All these efforts save time and reduce cost.

Immunoinformatics and modeling perspective of T cell epitope-based cancer immunotherapy: a holistic picture

Cancer immunotherapy is fast gaining global attention with its unique position as a potential therapy showing promise in cancer prevention and cure. It utilizes the natural system of immunity as opposed to chemotherapy and radiotherapy that utilize chemical drugs and radiation, respectively. Cancer immunotherapy essentially involves treatment and/or prevention with vaccines in the form of peptide vaccines (T and B cell epitopes), DNA vaccines and vaccination using whole tumor cells, dendritic cells, viral vectors, antibodies and adoptive transfer of T cells to harness the body's own immune system towards the targeting of cancer cells for destruction. Given the time, cost and labor involved in the vaccine discovery and development, researchers have evinced interest in the novel field of immunoinformatics to cut down the escalation of these critical resources. Immunoinformatics is a relatively new buzzword in the scientific circuit that is showing its potential and delivering on its promise in expediting the development of effective cancer immunotherapeutic agents. This review attempts to present a holistic picture of our race against cancer and time using the science and technology of immunoinformatics and molecular modeling in T cell epitope-based cancer immunotherapy. It also attempts to showcase some problem areas as well as novel ones waiting to be explored where development of novel immunoinformatics tools and simulations in the context of cancer immunotherapy would be highly welcome.

Towards in silico design of epitope-based vaccines

BACKGROUND

Epitope-based vaccines (EVs) make use of immunogenic peptides (epitopes) to trigger an immune response. Due to their manifold advantages, EVs have recently been attracting growing interest. The success of an EV is determined by the choice of epitopes used as a basis. However, the experimental discovery of candidate epitopes is expensive in terms of time and money. Furthermore, for the final choice of epitopes various immunological requirements have to be considered.

METHODS

Numerous in silico approaches exist that can guide the design of EVs. In particular, computational methods for MHC binding prediction have already become standard tools in immunology. Apart from binding prediction and prediction of antigen processing, methods for epitope design and selection have been suggested. We review these in silico approaches for epitope discovery and selection along with their strengths and weaknesses. Finally, we discuss some of the obvious problems in the design of EVs.

CONCLUSION

State-of-the-art in silico approaches to MHC binding prediction yield high accuracies. However, a more thorough understanding of the underlying biological processes and significant amounts of experimental data will be required for the validation and improvement of in silico approaches to the remaining aspects of EV design.

An introduction to epitope prediction methods and software

In this paper, current prediction methods and algorithms for both T- and B cell epitopes are reviewed, and a comprehensive summary of epitope prediction software and databases currently available online is also provided. This review can offer researchers in this field a sense of direction and insights for future work. However, our main purpose is to introduce clinical and basic biomedical researchers who are not familiar with these biological analysis tools and databases to these online resources and to provide guidance on how to use them effectively.

Prediction of MHC binding peptide using Gibbs motif sampler, weight matrix and artificial neural network

The identification of MHC restricted epitopes is an important goal in peptide based vaccine and diagnostic development. As wet lab experiments for identification of MHC binding peptide are expensive and time consuming, in silico tools have been developed as fast alternatives, however with low performance. In the present study, we used IEDB training and blind validation datasets for the prediction of peptide binding to fourteen human MHC class I and II molecules using Gibbs motif sampler, weight matrix and artificial neural network methods. As compare to MHC class I predictor based on sequence weighting (Aroc=0.95 and CC=0.56) and artificial neural network (Aroc=0.73 and CC=0.25), MHC class II predictor based on Gibbs sampler did not perform well (Aroc=0.62 and CC=0.19). The predictive accuracy of Gibbs motif sampler in identifying the 9-mer cores of a binding peptide to DRB1 alleles are also limited (40¢), however above the random prediction (14¢). Therefore, the size of dataset (training and validation) and the correct identification of the binding core are the two main factors limiting the performance of MHC class-II binding peptide prediction. Overall, these data suggest that there is substantial room to improve the quality of the core predictions using novel approaches that capture distinct features of MHC-peptide interactions than the current approaches.

Shift-invariant adaptive double threading: learning MHC II-peptide binding

The major histocompatibility complex (MHC) plays important roles in the workings of the human immune system. Specificity of MHC binding to peptide fragments from cellular and pathogens' proteins has been found to correlate with disease outcome and pathogen or cancer evolution. In this paper we propose a novel approach to predicting binding configurations and energies for MHC class II molecules, whose epitopes are generally predicted less well than the MHC I epitopes due in part to larger variation in bound peptide length. We treat the relative position of the peptide as a hidden variable, and model the ensemble of different binding configurations, rather than use a separate alignment procedure to narrow it down to one. Thus, our predictor infers a distribution over peptide positions from the MHC II and peptide sequences, and computes the total binding affinity. The training procedure iterates the predictions with re-estimation of the parameters of the binding groove model. For a given relative peptide position, any MHC class I prediction model can be used. Here we choose the physics based model of Jojic et al. (2006). We show that the parameters of the binding model can be learned efficiently from the training data and then used to estimate binding energies for previously untested peptides. Our technique performs on par with previous approaches to MHC II epitope prediction. Furthermore, our model choice allows generalization to new MHC class II alleles, which were not a part of the training set.

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.