Toward more accurate pan-specific MHC-peptide binding prediction: a review of current methods and tools

Development and Clinical Applications of Therapeutic Cancer Vaccines with Individualized and Shared Neoantigens

Neoantigens, presented as peptides on the surfaces of cancer cells, have recently been proposed as optimal targets for immunotherapy in clinical practice. The promising outcomes of neoantigen-based cancer vaccines have inspired enthusiasm for their broader clinical applications. However, the individualized tumor-specific antigens (TSA) entail considerable costs and time due to the variable immunogenicity and response rates of these neoantigens-based vaccines, influenced by factors such as neoantigen response, vaccine types, and combination therapy. Given the crucial role of neoantigen efficacy, a number of bioinformatics algorithms and pipelines have been developed to improve the accuracy rate of prediction through considering a series of factors involving in HLA-peptide-TCR complex formation, including peptide presentation, HLA-peptide affinity, and TCR recognition. On the other hand, shared neoantigens, originating from driver mutations at hot mutation spots (e.g., KRASG12D), offer a promising and ideal target for the development of therapeutic cancer vaccines. A series of clinical practices have established the efficacy of these vaccines in patients with distinct HLA haplotypes. Moreover, increasing evidence demonstrated that a combination of tumor associated antigens (TAAs) and neoantigens can also improve the prognosis, thus expand the repertoire of shared neoantigens for cancer vaccines. In this review, we provide an overview of the complex process involved in identifying personalized neoantigens, their clinical applications, advances in vaccine technology, and explore the therapeutic potential of shared neoantigen strategies.

ConvNeXt-MHC: improving MHC-peptide affinity prediction by structure-derived degenerate coding and the ConvNeXt model

Peptide binding to major histocompatibility complex (MHC) proteins plays a critical role in T-cell recognition and the specificity of the immune response. Experimental validation such peptides is extremely resource-intensive. As a result, accurate computational prediction of binding peptides is highly important, particularly in the context of cancer immunotherapy applications, such as the identification of neoantigens. In recent years, there is a significant need to continually improve the existing prediction methods to meet the demands of this field. We developed ConvNeXt-MHC, a method for predicting MHC-I-peptide binding affinity. It introduces a degenerate encoding approach to enhance well-established panspecific methods and integrates transfer learning and semi-supervised learning methods into the cutting-edge deep learning framework ConvNeXt. Comprehensive benchmark results demonstrate that ConvNeXt-MHC outperforms state-of-the-art methods in terms of accuracy. We expect that ConvNeXt-MHC will help us foster new discoveries in the field of immunoinformatics in the distant future. We constructed a user-friendly website at http://www.combio-lezhang.online/predict/, where users can access our data and application.

MHCII-peptide presentation: an assessment of the state-of-the-art prediction methods

Major histocompatibility complex Class II (MHCII) proteins initiate and regulate immune responses by presentation of antigenic peptides to CD4+ T-cells and self-restriction. The interactions between MHCII and peptides determine the specificity of the immune response and are crucial in immunotherapy and cancer vaccine design. With the ever-increasing amount of MHCII-peptide binding data available, many computational approaches have been developed for MHCII-peptide interaction prediction over the last decade. There is thus an urgent need to provide an up-to-date overview and assessment of these newly developed computational methods. To benchmark the prediction performance of these methods, we constructed an independent dataset containing binding and non-binding peptides to 20 human MHCII protein allotypes from the Immune Epitope Database, covering DP, DR and DQ alleles. After collecting 11 known predictors up to January 2022, we evaluated those available through a webserver or standalone packages on this independent dataset. The benchmarking results show that MixMHC2pred and NetMHCIIpan-4.1 achieve the best performance among all predictors. In general, newly developed methods perform better than older ones due to the rapid expansion of data on which they are trained and the development of deep learning algorithms. Our manuscript not only draws a full picture of the state-of-art of MHCII-peptide binding prediction, but also guides researchers in the choice among the different predictors. More importantly, it will inspire biomedical researchers in both academia and industry for the future developments in this field.

Historical perspective and future directions: computational science in immuno-oncology

Immuno-oncology holds promise for transforming patient care having achieved durable clinical response rates across a variety of advanced and metastatic cancers. Despite these achievements, only a minority of patients respond to immunotherapy, underscoring the importance of elucidating molecular mechanisms responsible for response and resistance to inform the development and selection of treatments. Breakthroughs in molecular sequencing technologies have led to the generation of an immense amount of genomic and transcriptomic sequencing data that can be mined to uncover complex tumor-immune interactions using computational tools. In this review, we discuss existing and emerging computational methods that contextualize the composition and functional state of the tumor microenvironment, infer the reactivity and clonal dynamics from reconstructed immune cell receptor repertoires, and predict the antigenic landscape for immune cell recognition. We further describe the advantage of multi-omics analyses for capturing multidimensional relationships and artificial intelligence techniques for integrating omics data with histopathological and radiological images to encapsulate patterns of treatment response and tumor-immune biology. Finally, we discuss key challenges impeding their widespread use and clinical application and conclude with future perspectives. We are hopeful that this review will both serve as a guide for prospective researchers seeking to use existing tools for scientific discoveries and inspire the optimization or development of novel tools to enhance precision, ultimately expediting advancements in immunotherapy that improve patient survival and quality of life.

RPEMHC: improved prediction of MHC-peptide binding affinity by a deep learning approach based on residue-residue pair encoding

MOTIVATION

Binding of peptides to major histocompatibility complex (MHC) molecules plays a crucial role in triggering T cell recognition mechanisms essential for immune response. Accurate prediction of MHC-peptide binding is vital for the development of cancer therapeutic vaccines. While recent deep learning-based methods have achieved significant performance in predicting MHC-peptide binding affinity, most of them separately encode MHC molecules and peptides as inputs, potentially overlooking critical interaction information between the two.

RESULTS

In this work, we propose RPEMHC, a new deep learning approach based on residue-residue pair encoding to predict the binding affinity between peptides and MHC, which encode an MHC molecule and a peptide as a residue-residue pair map. We evaluate the performance of RPEMHC on various MHC-II-related datasets for MHC-peptide binding prediction, demonstrating that RPEMHC achieves better or comparable performance against other state-of-the-art baselines. Moreover, we further construct experiments on MHC-I-related datasets, and experimental results demonstrate that our method can work on both two MHC classes. These extensive validations have manifested that RPEMHC is an effective tool for studying MHC-peptide interactions and can potentially facilitate the vaccine development.

AVAILABILITY

The source code of the method along with trained models is freely available at https://github.com/lennylv/RPEMHC.

DeepHLAPred: a deep learning-based method for non-classical HLA binder prediction

Human leukocyte antigen (HLA) is closely involved in regulating the human immune system. Despite great advance in detecting classical HLA Class I binders, there are few methods or toolkits for recognizing non-classical HLA Class I binders. To fill in this gap, we have developed a deep learning-based tool called DeepHLAPred. The DeepHLAPred used electron-ion interaction pseudo potential, integer numerical mapping and accumulated amino acid frequency as initial representation of non-classical HLA binder sequence. The deep learning module was used to further refine high-level representations. The deep learning module comprised two parallel convolutional neural networks, each followed by maximum pooling layer, dropout layer, and bi-directional long short-term memory network. The experimental results showed that the DeepHLAPred reached the state-of-the-art performanceson the cross-validation test and the independent test. The extensive test demonstrated the rationality of the DeepHLAPred. We further analyzed sequence pattern of non-classical HLA class I binders by information entropy. The information entropy of non-classical HLA binder sequence implied sequence pattern to a certain extent. In addition, we have developed a user-friendly webserver for convenient use, which is available at http://www.biolscience.cn/DeepHLApred/ . The tool and the analysis is helpful to detect non-classical HLA Class I binder. The source code and data is available at https://github.com/tangxingyu0/DeepHLApred .

DeepMHCI: an anchor position-aware deep interaction model for accurate MHC-I peptide binding affinity prediction

MOTIVATION

Computationally predicting major histocompatibility complex class I (MHC-I) peptide binding affinity is an important problem in immunological bioinformatics, which is also crucial for the identification of neoantigens for personalized therapeutic cancer vaccines. Recent cutting-edge deep learning-based methods for this problem cannot achieve satisfactory performance, especially for non-9-mer peptides. This is because such methods generate the input by simply concatenating the two given sequences: a peptide and (the pseudo sequence of) an MHC class I molecule, which cannot precisely capture the anchor positions of the MHC binding motif for the peptides with variable lengths. We thus developed an anchor position-aware and high-performance deep model, DeepMHCI, with a position-wise gated layer and a residual binding interaction convolution layer. This allows the model to control the information flow in peptides to be aware of anchor positions and model the interactions between peptides and the MHC pseudo (binding) sequence directly with multiple convolutional kernels.

RESULTS

The performance of DeepMHCI has been thoroughly validated by extensive experiments on four benchmark datasets under various settings, such as 5-fold cross-validation, validation with the independent testing set, external HPV vaccine identification, and external CD8+ epitope identification. Experimental results with visualization of binding motifs demonstrate that DeepMHCI outperformed all competing methods, especially on non-9-mer peptides binding prediction.

AVAILABILITY AND IMPLEMENTATION

DeepMHCI is publicly available at https://github.com/ZhuLab-Fudan/DeepMHCI.

Computational Methods in Immunology and Vaccinology: Design and Development of Antibodies and Immunogens

The design of new biomolecules able to harness immune mechanisms for the treatment of diseases is a prime challenge for computational and simulative approaches. For instance, in recent years, antibodies have emerged as an important class of therapeutics against a spectrum of pathologies. In cancer, immune-inspired approaches are witnessing a surge thanks to a better understanding of tumor-associated antigens and the mechanisms of their engagement or evasion from the human immune system. Here, we provide a summary of the main state-of-the-art computational approaches that are used to design antibodies and antigens, and in parallel, we review key methodologies for epitope identification for both B- and T-cell mediated responses. A special focus is devoted to the description of structure- and physics-based models, privileged over purely sequence-based approaches. We discuss the implications of novel methods in engineering biomolecules with tailored immunological properties for possible therapeutic uses. Finally, we highlight the extraordinary challenges and opportunities presented by the possible integration of structure- and physics-based methods with emerging Artificial Intelligence technologies for the prediction and design of novel antigens, epitopes, and antibodies.

CrossDome: an interactive R package to predict cross-reactivity risk using immunopeptidomics databases

T-cell-based immunotherapies hold tremendous potential in the fight against cancer, thanks to their capacity to specifically targeting diseased cells. Nevertheless, this potential has been tempered with safety concerns regarding the possible recognition of unknown off-targets displayed by healthy cells. In a notorious example, engineered T-cells specific to MAGEA3 (EVDPIGHLY) also recognized a TITIN-derived peptide (ESDPIVAQY) expressed by cardiac cells, inducing lethal damage in melanoma patients. Such off-target toxicity has been related to T-cell cross-reactivity induced by molecular mimicry. In this context, there is growing interest in developing the means to avoid off-target toxicity, and to provide safer immunotherapy products. To this end, we present CrossDome, a multi-omics suite to predict the off-target toxicity risk of T-cell-based immunotherapies. Our suite provides two alternative protocols, i) a peptide-centered prediction, or ii) a TCR-centered prediction. As proof-of-principle, we evaluate our approach using 16 well-known cross-reactivity cases involving cancer-associated antigens. With CrossDome, the TITIN-derived peptide was predicted at the 99+ percentile rank among 36,000 scored candidates (p-value < 0.001). In addition, off-targets for all the 16 known cases were predicted within the top ranges of relatedness score on a Monte Carlo simulation with over 5 million putative peptide pairs, allowing us to determine a cut-off p-value for off-target toxicity risk. We also implemented a penalty system based on TCR hotspots, named contact map (CM). This TCR-centered approach improved upon the peptide-centered prediction on the MAGEA3-TITIN screening (e.g., from 27th to 6th, out of 36,000 ranked peptides). Next, we used an extended dataset of experimentally-determined cross-reactive peptides to evaluate alternative CrossDome protocols. The level of enrichment of validated cases among top 50 best-scored peptides was 63% for the peptide-centered protocol, and up to 82% for the TCR-centered protocol. Finally, we performed functional characterization of top ranking candidates, by integrating expression data, HLA binding, and immunogenicity predictions. CrossDome was designed as an R package for easy integration with antigen discovery pipelines, and an interactive web interface for users without coding experience. CrossDome is under active development, and it is available at https://github.com/AntunesLab/crossdome.

Advancing our knowledge of antigen processing with computational modelling, structural biology, and immunology

Antigen processing is an immunological mechanism by which intracellular peptides are transported to the cell surface while bound to Major Histocompatibility Complex molecules, where they can be surveyed by circulating CD8+ or CD4+ T-cells, potentially triggering an immunological response. The antigen processing pathway is a complex multistage filter that refines a huge pool of potential peptide ligands derived from protein degradation into a smaller ensemble for surface presentation. Each stage presents unique challenges due to the number of ligands, the polymorphic nature of MHC and other protein constituents of the pathway and the nature of the interactions between them. Predicting the ensemble of displayed peptide antigens, as well as their immunogenicity, is critical for improving T cell vaccines against pathogens and cancer. Our predictive abilities have always been hindered by an incomplete empirical understanding of the antigen processing pathway. In this review, we highlight the role of computational and structural approaches in improving our understanding of antigen processing, including structural biology, computer simulation, and machine learning techniques, with a particular focus on the MHC-I pathway.

Artificial intelligence applied in neoantigen identification facilitates personalized cancer immunotherapy

The field of cancer neoantigen investigation has developed swiftly in the past decade. Predicting novel and true neoantigens derived from large multi-omics data became difficult but critical challenges. The rise of Artificial Intelligence (AI) or Machine Learning (ML) in biomedicine application has brought benefits to strengthen the current computational pipeline for neoantigen prediction. ML algorithms offer powerful tools to recognize the multidimensional nature of the omics data and therefore extract the key neoantigen features enabling a successful discovery of new neoantigens. The present review aims to outline the significant technology progress of machine learning approaches, especially the newly deep learning tools and pipelines, that were recently applied in neoantigen prediction. In this review article, we summarize the current state-of-the-art tools developed to predict neoantigens. The standard workflow includes calling genetic variants in paired tumor and blood samples, and rating the binding affinity between mutated peptide, MHC (I and II) and T cell receptor (TCR), followed by characterizing the immunogenicity of tumor epitopes. More specifically, we highlight the outstanding feature extraction tools and multi-layer neural network architectures in typical ML models. It is noted that more integrated neoantigen-predicting pipelines are constructed with hybrid or combined ML algorithms instead of conventional machine learning models. In addition, the trends and challenges in further optimizing and integrating the existing pipelines are discussed.

Immunoinformatics-Aided Design of a Peptide Based Multiepitope Vaccine Targeting Glycoproteins and Membrane Proteins against Monkeypox Virus

Monkeypox is a self-limiting zoonotic viral disease and causes smallpox-like symptoms. The disease has a case fatality ratio of 3-6% and, recently, a multi-country outbreak of the disease has occurred. The currently available vaccines that have provided immunization against monkeypox are classified as live attenuated vaccinia virus-based vaccines, which pose challenges of safety and efficacy in chronic infections. In this study, we have used an immunoinformatics-aided design of a multi-epitope vaccine (MEV) candidate by targeting monkeypox virus (MPXV) glycoproteins and membrane proteins. From these proteins, seven epitopes (two T-helper cell epitopes, four T-cytotoxic cell epitopes and one linear B cell epitopes) were finally selected and predicted as antigenic, non-allergic, interferon-γ activating and non-toxic. These epitopes were linked to adjuvants to design a non-allergic and antigenic candidate MPXV-MEV. Further, molecular docking and molecular dynamics simulations predicted stable interactions between predicted MEV and human receptor TLR5. Finally, the immune-simulation analysis showed that the candidate MPXV-MEV could elicit a human immune response. The results obtained from these in silico experiments are promising but require further validation through additional in vivo experiments.

Design of a multi-epitope vaccine against the pathogenic fungi Candida tropicalis using an in silico approach

Background

Candida tropicalis causes tropical invasive fungal infections, with a high mortality. This fungus has been found to be resistant to antifungal classes such as azoles, echinocandins, and polyenes in several studies. As a result, it is vital to identify novel approaches to prevent and treat C. tropicalis infections. In this study, an in silico technique was utilized to deduce and evaluate a powerful multivalent epitope-based vaccine against C. tropicalis, which targets the secreted aspartic protease 2 (SAP2) protein. This protein is implicated in virulence and host invasion.

Results

By focusing on the Sap2 protein, 11 highly antigenic, non-allergic, non-toxic, and conserved epitopes were identified. These were subsequently paired with RS09 and flagellin adjuvants, as well as a pan HLA DR-binding epitope (PADRE) sequence to create a vaccine candidate that elicited both cell-mediated and humoral immune responses. It was projected that the vaccine design would be soluble, stable, antigenic, and non-allergic. Ramachandran plot analysis was applied to validate the vaccine construct’s 3-dimensional model. The vaccine construct was tested (at 100 ns) using molecular docking and molecular dynamics simulations, which demonstrated that it can stably connect with MHC-I and Toll-like receptor molecules. Based on in silico studies, we have shown that the vaccine construct can be expressed in E. coli. We surmise that the vaccine design is unrelated to any human proteins, indicating that it is safe to use.

Conclusions

The vaccine design looks to be an effective option for preventing C. tropicalis infections, based on the outcomes of the studies. A fungal vaccine can be proposed as prophylactic medicine and could provide initial protection as sometimes diagnosis of infection could be challenging. However, more in vitro and in vivo research is needed to prove the efficacy and safety of the proposed vaccine design.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43141-022-00415-3.

SILVI, an open-source pipeline for T-cell epitope selection

High-throughput screening of available genomic data and identification of potential antigenic candidates have promoted the development of epitope-based vaccines and therapeutics. Several immunoinformatic tools are available to predict potential epitopes and other immunogenicity-related features, yet it is still challenging and time-consuming to compare and integrate results from different algorithms. We developed the R script SILVI (short for: from in silico to in vivo), to assist in the selection of the potentially most immunogenic T-cell epitopes from Human Leukocyte Antigen (HLA)-binding prediction data. SILVI merges and compares data from available HLA-binding prediction servers, and integrates additional relevant information of predicted epitopes, namely BLASTp alignments with host proteins and physical-chemical properties. The two default criteria applied by SILVI and additional filtering allow the fast selection of the most conserved, promiscuous, strong binding T-cell epitopes. Users may adapt the script at their discretion as it is written in open-source R language. To demonstrate the workflow and present selection options, SILVI was used to integrate HLA-binding prediction results of three example proteins, from viral, bacterial and parasitic microorganisms, containing validated epitopes included in the Immune Epitope Database (IEDB), plus the Human Papillomavirus (HPV) proteome. Applying different filters on predicted IC50, hydrophobicity and mismatches with host proteins allows to significantly reduce the epitope lists with favourable sensitivity and specificity to select immunogenic epitopes. We contemplate SILVI will assist T-cell epitope selections and can be continuously refined in a community-driven manner, helping the improvement and design of peptide-based vaccines or immunotherapies. SILVI development version is available at: github.com/JoanaPissarra/SILVI2020 and https://doi.org/10.5281/zenodo.6865909.

CD8TCEI-EukPath: A Novel Predictor to Rapidly Identify CD8+ T-Cell Epitopes of Eukaryotic Pathogens Using a Hybrid Feature Selection Approach

Computational prediction to screen potential vaccine candidates has been proven to be a reliable way to provide guarantees for vaccine discovery in infectious diseases. As an important class of organisms causing infectious diseases, pathogenic eukaryotes (such as parasitic protozoans) have evolved the ability to colonize a wide range of hosts, including humans and animals; meanwhile, protective vaccines are urgently needed. Inspired by the immunological idea that pathogen-derived epitopes are able to mediate the CD8⁺ T-cell-related host adaptive immune response and with the available positive and negative CD8⁺ T-cell epitopes (TCEs), we proposed a novel predictor called CD8TCEI-EukPath to detect CD8⁺ TCEs of eukaryotic pathogens. Our method integrated multiple amino acid sequence-based hybrid features, employed a well-established feature selection technique, and eventually built an efficient machine learning classifier to differentiate CD8⁺ TCEs from non-CD8⁺ TCEs. Based on the feature selection results, 520 optimal hybrid features were used for modeling by utilizing the LightGBM algorithm. CD8TCEI-EukPath achieved impressive performance, with an accuracy of 79.255% in ten-fold cross-validation and an accuracy of 78.169% in the independent test. Collectively, CD8TCEI-EukPath will contribute to rapidly screening epitope-based vaccine candidates, particularly from large peptide-coding datasets. To conduct the prediction of CD8⁺ TCEs conveniently, an online web server is freely accessible (http://lab.malab.cn/∼hrs/CD8TCEI-EukPath/).

Molecular mimicry between tumor associated antigens and microbiota-derived epitopes

Background

The gut microbiota profile is unique for each individual and are composed by different bacteria species according to individual birth-to-infant transitions. In the last years, the local and systemic effects of microbiota on cancer onset, progression and response to treatments, such as immunotherapies, has been extensively described. Here we offer a new perspective, proposing a role for the microbiota based on the molecular mimicry of tumor associated antigens by microbiome-associated antigens.

Methods

In the present study we looked for homology between published TAAs and non-self microbiota-derived epitopes. Blast search for sequence homology was combined with extensive bioinformatics analyses.

Results

Several evidences for homology between TAAs and microbiota-derived antigens have been found. Strikingly, three cases of 100% homology between the paired sequences has been identified. The predicted average affinity to HLA molecules of microbiota-derived antigens is very high (< 100 nM). The structural conformation of the microbiota-derived epitopes is, in general, highly similar to the corresponding TAA. In some cases, it is identical and contact areas with both HLA and TCR chains are indistinguishable. Moreover, the spatial conformation of TCR-facing residues can be identical in paired TAA and microbiota-derived epitopes, with exactly the same values of planar as well as dihedral angles.

Conclusions

The data reported in the present study show for the first time the high homology in the linear sequence as well as in structure and conformation between TAAs and peptides derived from microbiota species of the Firmicutes and the Bacteroidetes phyla, which together account for 90% of gut microbiota. Cross-reacting CD8⁺ T cell responses are very likely induced. Therefore, the anti-microbiota T cell memory may turn out to be an anti-cancer T cell memory, able to control the growth of a cancer developed during the lifetime if the expressed TAA is similar to the microbiota epitope. This may ultimately represent a relevant selective advantage for cancer patients and may lead to a novel preventive anti-cancer vaccine strategy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03512-6.

DeepSeqPanII: An Interpretable Recurrent Neural Network Model With Attention Mechanism for Peptide-HLA Class II Binding Prediction

Human leukocyte antigen (HLA) complex molecules play an essential role in immune interactions by presenting peptides on the cell surface to T cells. With significant deep learning progress, a series of neural network-based models have been proposed and demonstrated with their excellent performances for peptide-HLA class I binding prediction. However, there is still a lack of effective binding prediction models for HLA class II protein binding with peptides due to its inherent challenges. We present a novel sequence-based pan-specific neural network structure, DeepSeaPanII, for peptide-HLA class II binding prediction in this work. Our model is an end-to-end neural network model without the need for pre-or post-processing on input samples compared with existing pan-specific models. Besides state-of-the-art performance in binding affinity prediction, DeepSeqPanII can also extract biological insight on the binding mechanism over the peptide by its attention mechanism-based binding core prediction capability. The leave-one-allele-out cross-validation and benchmark evaluation results show that our proposed network model achieved state-of-the-art performance in HLA-II peptide binding. The source code and trained models are freely available at https://github.com/pcpLiu/DeepSeqPanII.

Prediction and characterization of the T cell epitopes for the major soybean protein allergens using bioinformatics approaches

Protein allergens is a health risk for consumption of soybeans. To understand allerginicity mechanism, T cell epitopes of 7 soybean allergens were predicted and screened by abilities to induce cytokine interleukin (IL) 4. The relationships among amino acid composition, properties, allergenicity, and pepsin hydrolysis sites were analyzed. Among the 138 T cell epitopes identified, YIKDVFRVIPSEVLS, KDVFRVIPSEVLSNS, DVFRVIPSEVLSNSY of Gly m 6.0501 (P04347), and AKADALFKAIEAYLL, ADALFKAIEAYLLAH of Gly m 4.0101 (P26987) were the most possible epitope candidates. In T cell epitopes pattern, the frequencies of amino acids Q, D, E, P, and G decreased, while F, I, N, V, K, H, A, L, and S increased. Hydrophobic residues at positions p1 and p2 and positively charged residues in positions p13 might contribute to allergenicity. Most of epitopes could be hydrolyzed by pepsin into small polypeptides within 12 residues length, and the anti-digestive epitope regions contained I, V, S, N, and Q residues. T cell epitopes EEQRQQEGVIVELSK from Gly m 5.03 (P25974) showed resistance to pepsin hydrolysis and would cause a higher Th2 cell response. This research provides basis for the development of hypoallergenic soybean products in the soybean industry as well as for the immunotherapy design for protein allergy.

Peptide-Based Vaccines for Tuberculosis

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis. As a result of the coronavirus disease 2019 (COVID-19) pandemic, the global TB mortality rate in 2020 is rising, making TB prevention and control more challenging. Vaccination has been considered the best approach to reduce the TB burden. Unfortunately, BCG, the only TB vaccine currently approved for use, offers some protection against childhood TB but is less effective in adults. Therefore, it is urgent to develop new TB vaccines that are more effective than BCG. Accumulating data indicated that peptides or epitopes play essential roles in bridging innate and adaptive immunity and triggering adaptive immunity. Furthermore, innovations in bioinformatics, immunoinformatics, synthetic technologies, new materials, and transgenic animal models have put wings on the research of peptide-based vaccines for TB. Hence, this review seeks to give an overview of current tools that can be used to design a peptide-based vaccine, the research status of peptide-based vaccines for TB, protein-based bacterial vaccine delivery systems, and animal models for the peptide-based vaccines. These explorations will provide approaches and strategies for developing safer and more effective peptide-based vaccines and contribute to achieving the WHO's End TB Strategy.

Machine Learning Techniques for the Prediction of B-Cell and T-Cell Epitopes as Potential Vaccine Targets with a Specific Focus on SARS-CoV-2 Pathogen: A Review

The only part of an antigen (a protein molecule found on the surface of a pathogen) that is composed of epitopes specific to T and B cells is recognized by the human immune system (HIS). Identification of epitopes is considered critical for designing an epitope-based peptide vaccine (EBPV). Although there are a number of vaccine types, EBPVs have received less attention thus far. It is important to mention that EBPVs have a great deal of untapped potential for boosting vaccination safety-they are less expensive and take a short time to produce. Thus, in order to quickly contain global pandemics such as the ongoing outbreak of coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), as well as epidemics and endemics, EBPVs are considered promising vaccine types. The high mutation rate of SARS-CoV-2 has posed a great challenge to public health worldwide because either the composition of existing vaccines has to be changed or a new vaccine has to be developed to protect against its different variants. In such scenarios, time being the critical factor, EBPVs can be a promising alternative. To design an effective and viable EBPV against different strains of a pathogen, it is important to identify the putative T- and B-cell epitopes. Using the wet-lab experimental approach to identify these epitopes is time-consuming and costly because the experimental screening of a vast number of potential epitope candidates is required. Fortunately, various available machine learning (ML)-based prediction methods have reduced the burden related to the epitope mapping process by decreasing the potential epitope candidate list for experimental trials. Moreover, these methods are also cost-effective, scalable, and fast. This paper presents a systematic review of various state-of-the-art and relevant ML-based methods and tools for predicting T- and B-cell epitopes. Special emphasis is placed on highlighting and analyzing various models for predicting epitopes of SARS-CoV-2, the causative agent of COVID-19. Based on the various methods and tools discussed, future research directions for epitope prediction are presented.

An Update on "Reverse Vaccinology": The Pathway from Genomes and Epitope Predictions to Tailored, Recombinant Vaccines

In this chapter, we review the computational approaches that have led to a new generation of vaccines in recent years. There are many alternative routes to develop vaccines based on the concept of reverse vaccinology. They all follow the same basic principles-mining available genome and proteome information for antigen candidates, and recombinantly expressing them for vaccine production. Some of the same principles have been used successfully for cancer therapy approaches. In this review, we focus on infectious diseases, describing the general workflow from bioinformatic predictions of antigens and epitopes down to examples where such predictions have been used successfully for vaccine development.

Identification of Immunogenic MHC Class II Human HER3 Peptides that Mediate Anti-HER3 CD4+ Th1 Responses and Potential Use as a Cancer Vaccine

The HER3/ERBB3 receptor is an oncogenic receptor tyrosine kinase that forms heterodimers with EGFR family members and is overexpressed in numerous cancers. HER3 overexpression associates with reduced survival and acquired resistance to targeted therapies, making it a potential therapeutic target in multiple cancer types. Here, we report on immunogenic, promiscuous MHC class II-binding HER3 peptides, which can generate HER3-specific CD4⁺ Th1 antitumor immune responses. Using an overlapping peptide screening methodology, we identified nine MHC class II-binding HER3 epitopes that elicited specific Th1 immune response in both healthy donors and breast cancer patients. Most of these peptides were not identified by current binding algorithms. Homology assessment of amino acid sequence BLAST showed >90% sequence similarity between human and murine HER3/ERBB3 peptide sequences. HER3 peptide-pulsed dendritic cell vaccination resulted in anti-HER3 CD4⁺ Th1 responses that prevented tumor development, significantly delayed tumor growth in prevention models, and caused regression in multiple therapeutic models of HER3-expressing murine tumors, including mammary carcinoma and melanoma. Tumors were robustly infiltrated with CD4⁺ T cells, suggesting their key role in tumor rejection. Our data demonstrate that class II HER3 promiscuous peptides are effective at inducing HER3-specific CD4⁺ Th1 responses and suggest their applicability in immunotherapies for human HER3-overexpressing tumors.

OUP accepted manuscript

Motivation

Computationally predicting major histocompatibility complex (MHC)-peptide binding affinity is an important problem in immunological bioinformatics. Recent cutting-edge deep learning-based methods for this problem are unable to achieve satisfactory performance for MHC class II molecules. This is because such methods generate the input by simply concatenating the two given sequences: (the estimated binding core of) a peptide and (the pseudo sequence of) an MHC class II molecule, ignoring biological knowledge behind the interactions of the two molecules. We thus propose a binding core-aware deep learning-based model, DeepMHCII, with a binding interaction convolution layer, which allows to integrate all potential binding cores (in a given peptide) with the MHC pseudo (binding) sequence, through modeling the interaction with multiple convolutional kernels.

Results

Extensive empirical experiments with four large-scale datasets demonstrate that DeepMHCII significantly outperformed four state-of-the-art methods under numerous settings, such as 5-fold cross-validation, leave one molecule out, validation with independent testing sets and binding core prediction. All these results and visualization of the predicted binding cores indicate the effectiveness of our model, DeepMHCII, and the importance of properly modeling biological facts in deep learning for high predictive performance and efficient knowledge discovery.

Availability and implementation

DeepMHCII is publicly available at https://github.com/yourh/DeepMHCII.

Supplementary information

Supplementary data are available at Bioinformatics online.

Long-term memory T cells as preventive anticancer immunity elicited by TuA-derived heteroclitic peptides

The host's immune system may be primed against antigens during the lifetime (e.g. microorganisms antigens-MoAs), and swiftly recalled upon growth of a tumor expressing antigens similar in sequence and structure. C57BL/6 mice were immunized in a preventive setting with tumor antigens (TuAs) or corresponding heteroclitic peptides specific for TC-1 and B16 cell lines. Immediately or 2-months after the end of the vaccination protocol, animals were implanted with cell lines. The specific anti-vaccine immune response as well as tumor growth were regularly evaluated for 2 months post-implantation. The preventive vaccination with TuA or their heteroclitic peptides (hPep) was able to delay (B16) or completely suppress (TC-1) tumor growth when cancer cells were implanted immediately after the end of the vaccination. More importantly, TC-1 tumor growth was significantly delayed, and suppressed in 6/8 animals, also when cells were implanted 2-months after the end of the vaccination. The vaccine-specific T cell response provided a strong immune correlate to the pattern of tumor growth. A preventive immunization with heteroclitic peptides resembling a TuA is able to strongly delay or even suppress tumor growth in a mouse model. More importantly, the same effect is observed also when tumor cells are implanted 2 months after the end of vaccination, which corresponds to 8 - 10 years in human life. The observed potent tumor control indicates that a memory T cell immunity elicited during the lifetime by a antigens similar to a TuA, i.e. viral antigens, may ultimately represent a great advantage for cancer patients and may lead to a novel preventive anti-cancer vaccine strategy.

Computational cancer neoantigen prediction: current status and recent advances

Over the last few decades, immunotherapy has shown significant therapeutic efficacy in a broad range of cancer types. Antitumor immune responses are contingent on the recognition of tumor-specific antigens, which are termed neoantigens. Tumor neoantigens are ideal targets for immunotherapy since they can be recognized as non-self antigens by the host immune system and thus are able to elicit an antitumor T-cell response. There are an increasing number of studies that highlight the importance of tumor neoantigens in immunoediting and in the sensitivity to immune checkpoint blockade. Therefore, one of the most fundamental tasks in the field of immuno-oncology research is the identification of patient-specific neoantigens. To this end, a plethora of computational approaches have been developed in order to predict tumor-specific aberrant peptides and quantify their likelihood of binding to patients' human leukocyte antigen molecules in order to be recognized by T cells. In this review, we systematically summarize and present the most recent advances in computational neoantigen prediction, and discuss the challenges and novel methods that are being developed to resolve them.

•

Tumors have the ability to acquire immune escape mechanisms.

•

Tumor-specific aberrant peptides (neoantigens) can elicit an immune response by the host immune system.

•

The identification of neoantigens is one of the most fundamental tasks in the field of immuno-oncology research.

•

A plethora of computational approaches have been developed in order to predict patient-specificneoantigens.

Screening and Analysis of Anaplasma marginale Tunisian Isolates Reveal the Diversity of lipA Phylogeographic Marker and the Conservation of OmpA Protein Vaccine Candidate

Bovine anaplasmosis caused by Anaplasma marginale is a disease responsible for serious animal health problems and great economic losses all over the world. Thereby, the identification of A. marginale isolates from various bioclimatic areas in each country, the phylogeographic analysis of these isolates based on the most informative markers, and the evaluation of the most promising candidate antigens are crucial steps in developing effective vaccines against a wide range of A. marginale strains. In order to contribute to this challenge, a total of 791 bovine samples from various bioclimatic areas of Tunisia were tested for the occurrence of A. marginale DNA through msp4 gene fragment amplification. Phylogeographic analysis was performed by using lipA and sucB gene analyses, and the genetic relationship with previously characterized A. marginale isolates and strains was analyzed by applying similarity comparison and phylogenetic analysis. To evaluate the conservation of OmpA protein vaccine candidate, almost complete ompA nucleotide sequences were also obtained from Tunisian isolates, and various bioinformatics software were used in order to analyze the physicochemical properties and the secondary and tertiary structures of their deduced proteins and to predict their immunodominant epitopes of B and T cells. A. marginale DNA was detected in 19 bovine samples (2.4%). Risk factor analysis shows that cattle derived from subhumid bioclimatic area were more infected than those that originated from other areas. The analysis of lipA phylogeographic marker indicated a higher diversity of Tunisian A. marginale isolates compared with other available worldwide isolates and strains. Molecular, phylogenetic, and immuno-informatics analyses of the vaccine candidate OmpA protein demonstrated that this antigen and its predicted immunodominant epitopes of B and T cells appear to be highly conserved between Tunisian isolates and compared with isolates from other countries, suggesting that the minimal intraspecific modifications will not affect the potential cross-protective capacity of humoral and cell-mediated immune responses against multiple A. marginale worldwide strains.

Connecting MHC-I-binding motifs with HLA alleles via deep learning

The selection of peptides presented by MHC molecules is crucial for antigen discovery. Previously, several predictors have shown impressive performance on binding affinity. However, the decisive MHC residues and their relation to the selection of binding peptides are still unrevealed. Here, we connected HLA alleles with binding motifs via our deep learning-based framework, MHCfovea. MHCfovea expanded the knowledge of MHC-I-binding motifs from 150 to 13,008 alleles. After clustering N-terminal and C-terminal sub-motifs on both observed and unobserved alleles, MHCfovea calculated the hyper-motifs and the corresponding allele signatures on the important positions to disclose the relation between binding motifs and MHC-I sequences. MHCfovea delivered 32 pairs of hyper-motifs and allele signatures (HLA-A: 13, HLA-B: 12, and HLA-C: 7). The paired hyper-motifs and allele signatures disclosed the critical polymorphic residues that determine the binding preference, which are believed to be valuable for antigen discovery and vaccine design when allele specificity is concerned.

Ko-Han Lee et al. develop MHCfovea, a machine-learning method for predicting peptide-binding by MHC molecules and inferring peptide motifs and MHC allele signatures. They demonstrate that MHCfovea is capable of detecting meaningful hyper-motifs and allele signatures, making it a useful resource for the community.

MHCVision: estimation of global and local false discovery rate for MHC class I peptide binding prediction

Motivation

MHC-peptide binding prediction has been widely used for understanding the immune response of individuals or populations, each carrying different MHC molecules as well as for the development of immunotherapeutics. The results from MHC-peptide binding prediction tools are mostly reported as a predicted binding affinity (IC₅₀) and the percentile rank score, and global thresholds e.g. IC₅₀ value < 500 nM or percentile rank < 2% are generally recommended for distinguishing binding peptides from non-binding peptides. However, it is difficult to evaluate statistically the probability of an individual peptide binding prediction to be true or false solely considering predicted scores. Therefore, statistics describing the overall global false discovery rate (FDR) and local FDR, also called posterior error probability (PEP) are required to give statistical context to the natively produced scores.

Result

We have developed an algorithm and code implementation, called MHCVision, for estimation of FDR and PEP values for the predicted results of MHC-peptide binding prediction from the NetMHCpan tool. MHCVision performs parameter estimation using a modified expectation maximization framework for a two-component beta mixture model, representing the distribution of true and false scores of the predicted dataset. We can then estimate the PEP of an individual peptide’s predicted score, and conversely the probability that it is true. We demonstrate that the use of global FDR and PEP estimation can provide a better trade-off between sensitivity and precision over using currently recommended thresholds from tools.

Availability and implementation

https://github.com/PGB-LIV/MHCVision.

Supplementary information

Supplementary data are available at Bioinformatics online.

Unraveling the Allosteric Communication Mechanisms in T-Cell Receptor-Peptide-Loaded Major Histocompatibility Complex Dynamics Using Molecular Dynamics Simulations: An Approach Based on Dynamic Cross Correlation Maps and Residue Interaction Energy Calculations

Antigen presentation by major histocompatibility complex (MHC) proteins to T-cell receptors (TCRs) plays a crucial role in triggering the adaptive immune response. Most of our knowledge on TCR-peptide-loaded major histocompatibility complex (pMHC) interaction stemmed from experiments yielding static structures, yet the dynamic aspects of this molecular interaction are equally important to understand the underlying molecular mechanisms and to develop treatment strategies against diseases such as cancer and autoimmune diseases. To this end, computational biophysics studies including all-atom molecular dynamics simulations have provided useful insights; however, we still lack a basic understanding of an overall allosteric mechanism that results in conformational changes in the TCR and subsequent T-cell activation. Previous hydrogen-deuterium exchange and nuclear magnetic resonance studies provided clues regarding these molecular mechanisms, including global rigidification and allosteric effects on the constant domain of TCRs away from the pMHC interaction site. Here, we show that molecular dynamics simulations can be used to identify how this overall rigidification may be related to the allosteric communication within TCRs upon pMHC interaction via essential dynamics and nonbonded residue-residue interaction energy analyses. The residues taking part in the rigidification effect are highlighted with an intricate analysis on residue interaction changes, which lead to a detailed outline of the complex formation event. Our results indicate that residues of the Cβ domain of TCRs show significant differences in their nonbonded interactions upon complex formation. Moreover, the dynamic cross correlations between these residues are also increased, in line with their nonbonded interaction energy changes. Altogether, our approach may be valuable for elucidating intramolecular allosteric changes in the TCR structure upon pMHC interaction in molecular dynamics simulations.

Bioinformatic prediction of immunodominant regions in spike protein for early diagnosis of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

Background

To contain the pandemics caused by SARS-CoV-2, early detection approaches with high accuracy and accessibility are critical. Generating an antigen-capture based detection system would be an ideal strategy complementing the current methods based on nucleic acids and antibody detection. The spike protein is found on the outside of virus particles and appropriate for antigen detection.

Methods

In this study, we utilized bioinformatics approaches to explore the immunodominant fragments on spike protein of SARS-CoV-2.

Results

The S1 subunit of spike protein was identified with higher sequence specificity. Three immunodominant fragments, Spike_56-94, Spike_199-264, and Spike_577-612, located at the S1 subunit were finally selected via bioinformatics analysis. The glycosylation sites and high-frequency mutation sites on spike protein were circumvented in the antigen design. All the identified fragments present qualified antigenicity, hydrophilicity, and surface accessibility. A recombinant antigen with a length of 194 amino acids (aa) consisting of the selected immunodominant fragments as well as a universal Th epitope was finally constructed.

Conclusion

The recombinant peptide encoded by the construct contains multiple immunodominant epitopes, which is expected to stimulate a strong immune response in mice and generate qualified antibodies for SARS-CoV-2 detection.

A Proteome-Wide Immunoinformatics Tool to Accelerate T-Cell Epitope Discovery and Vaccine Design in the Context of Emerging Infectious Diseases: An Ethnicity-Oriented Approach

Emerging infectious diseases (EIDs) caused by viruses are increasing in frequency, causing a high disease burden and mortality world-wide. The COVID-19 pandemic caused by the novel SARS-like coronavirus (SARS-CoV-2) underscores the need to innovate and accelerate the development of effective vaccination strategies against EIDs. Human leukocyte antigen (HLA) molecules play a central role in the immune system by determining the peptide repertoire displayed to the T-cell compartment. Genetic polymorphisms of the HLA system thus confer a strong variability in vaccine-induced immune responses and may complicate the selection of vaccine candidates, because the distribution and frequencies of HLA alleles are highly variable among different ethnic groups. Herein, we build on the emerging paradigm of rational epitope-based vaccine design, by describing an immunoinformatics tool (Predivac-3.0) for proteome-wide T-cell epitope discovery that accounts for ethnic-level variations in immune responsiveness. Predivac-3.0 implements both CD8+ and CD4+ T-cell epitope predictions based on HLA allele frequencies retrieved from the Allele Frequency Net Database. The tool was thoroughly assessed, proving comparable performances (AUC ~0.9) against four state-of-the-art pan-specific immunoinformatics methods capable of population-level analysis (NetMHCPan-4.0, Pickpocket, PSSMHCPan and SMM), as well as a strong accuracy on proteome-wide T-cell epitope predictions for HIV-specific immune responses in the Japanese population. The utility of the method was investigated for the COVID-19 pandemic, by performing in silico T-cell epitope mapping of the SARS-CoV-2 spike glycoprotein according to the ethnic context of the countries where the ChAdOx1 vaccine is currently initiating phase III clinical trials. Potentially immunodominant CD8+ and CD4+ T-cell epitopes and population coverages were predicted for each population (the Epitope Discovery mode), along with optimized sets of broadly recognized (promiscuous) T-cell epitopes maximizing coverage in the target populations (the Epitope Optimization mode). Population-specific epitope-rich regions (T-cell epitope clusters) were further predicted in protein antigens based on combined criteria of epitope density and population coverage. Overall, we conclude that Predivac-3.0 holds potential to contribute in the understanding of ethnic-level variations of vaccine-induced immune responsiveness and to guide the development of epitope-based next-generation vaccines against emerging pathogens, whose geographic distributions and populations in need of vaccinations are often well-defined for regional epidemics.

Anthem: a user customised tool for fast and accurate prediction of binding between peptides and HLA class I molecules

Neopeptide-based immunotherapy has been recognised as a promising approach for the treatment of cancers. For neopeptides to be recognised by CD8+ T cells and induce an immune response, their binding to human leukocyte antigen class I (HLA-I) molecules is a necessary first step. Most epitope prediction tools thus rely on the prediction of such binding. With the use of mass spectrometry, the scale of naturally presented HLA ligands that could be used to develop such predictors has been expanded. However, there are rarely efforts that focus on the integration of these experimental data with computational algorithms to efficiently develop up-to-date predictors. Here, we present Anthem for accurate HLA-I binding prediction. In particular, we have developed a user-friendly framework to support the development of customisable HLA-I binding prediction models to meet challenges associated with the rapidly increasing availability of large amounts of immunopeptidomic data. Our extensive evaluation, using both independent and experimental datasets shows that Anthem achieves an overall similar or higher area under curve value compared with other contemporary tools. It is anticipated that Anthem will provide a unique opportunity for the non-expert user to analyse and interpret their own in-house or publicly deposited datasets.

An Immunoinformatics Study to Predict Epitopes in the Envelope Protein of SARS-CoV-2

COVID-19 is a new viral emergent disease caused by a novel strain of coronavirus. This virus has caused a huge problem in the world as millions of people are affected by this disease. We aimed at designing a peptide vaccine for COVID-19 particularly for the envelope protein using computational methods to predict epitopes inducing the immune system. The envelope protein sequence of SARS-CoV-2 has been retrieved from the NCBI database. The bioinformatics analysis was carried out by using the Immune Epitope Database (IEDB) to predict B- and T-cell epitopes. The predicted HTL and CTL epitopes were docked with HLA alleles and binding energies were evaluated. The allergenicity of predicted epitopes was analyzed, the conservancy analysis was performed, and the population coverage was determined throughout the world. Some overlapped CTL, HTL, and B-cell epitopes were suggested to become a universal candidate for peptide-based vaccine against COVID-19. This vaccine peptide could simultaneously elicit humoral and cell-mediated immune responses. We hope to confirm our findings by adding complementary steps of both in vitro and in vivo studies to support this new universal predicted candidate.

Elucidating the Efficacy of Vaccination against Vibriosis in Lates calcarifer Using Two Recombinant Protein Vaccines Containing the Outer Membrane Protein K (r-OmpK) of Vibrio alginolyticus and the DNA Chaperone J (r-DnaJ) of Vibrio harveyi

Recombinant cell vaccines expressing the OmpK and DnaJ of Vibrio were developed and subsequently, a vaccination efficacy trial was carried out on juvenile seabass (~5 cm; ~20 g). The fish were divided into 5 groups of 50 fish per group, kept in triplicate. Groups 1 and 2 were injected with 10⁷ CFU/mL of the inactivated recombinant cells vaccines, the pET-32/LIC-OmpK and pET-32/LIC-DnaJ, respectively. Group 3 was similarly injected with 10⁷ CFU/mL of inactivated E. coli BL21 (DE3), Group 4 with 10⁷ CFU/mL of formalin killed whole cells V. harveyi, and Group 5 with PBS solution. Serum, mucus, and gut lavage were used to determine the antibody levels before all fish were challenged with V. harveyi, V. alginolyticus, and V. parahemolyticus, respectively on day 15 post-vaccination. There was significant increase in the serum and gut lavage antibody titers in the juvenile seabass vaccinated with r-OmpK vaccine. In addition, there was an up-regulation for TLR2, MyD88, and MHCI genes in the kidney and intestinal tissues of r-OmpK vaccinated fish. At the same time, r-OmpK triggered higher expression level of interleukin IL-10, IL-8, IL-1ß in the spleen, intestine, and kidney compared to r-DnaJ. Overall, r-OmpK and r-DnaJ triggered protection by curbing inflammation and strengthening the adaptive immune response. Vaccinated fish also demonstrated strong cross protection against heterologous of Vibrio isolates, the V. harveyi, V. alginolyticus, and V. parahaemolyticus. The fish vaccinated with r-OmpK protein were completely protected with a relative per cent of survival (RPS) of 90 percent against V. harveyi and 100 percent against V. alginolyticus and V. parahaemolyticus. A semi-quantitative PCR detection of Vibrio spp. from the seawater containing the seabass also revealed that vaccination resulted in reduction of pathogen shedding. In conclusion, our results suggest r-OmpK as a candidate vaccine molecule against multiple Vibrio strain to prevent vibriosis in marine fish.

Alkali Metal Cationization of Tumor-associated Antigen Peptides for Improved Dissociation and Measurement by Differential Ion Mobility-Mass Spectrometry

Tandem mass spectrometry (MS/MS) is a highly sensitive and selective method for the detection of tumor-associated peptide antigens. These short, nontryptic sequences may lack basic residues, resulting in the formation of predominantly [peptide + H]+ ions in electrospray. These singly charged ions tend to undergo inefficient dissociation, leading to issues in sequence determination. Addition of alkali metal salts to the electrospray solvent can drive the formation of [peptide + H + metal]2+ ions that have enhanced dissociation characteristics relative to [peptide + H]+ ions. Both previously identified tumor-associated antigens and predicted neoantigen sequences were investigated. The previously reported rearrangement mechanism in MS/MS of sodium-cationized peptides is applied here to demonstrate complete C-terminal sequencing of tumor-associated peptide antigens. Differential ion mobility spectrometry (DIMS) is shown to selectively enrich [peptide + H + metal]2+ species by filtering out singly charged interferences at relatively low field strengths, offsetting the decrease in signal intensity associated with the use of alkali metal cations.

Immunoinformatics: Predicting Peptide-MHC Binding

Immunoinformatics is a discipline that applies methods of computer science to study and model the immune system. A fundamental question addressed by immunoinformatics is how to understand the rules of antigen presentation by MHC molecules to T cells, a process that is central to adaptive immune responses to infections and cancer. In the modern era of personalized medicine, the ability to model and predict which antigens can be presented by MHC is key to manipulating the immune system and designing strategies for therapeutic intervention. Since the MHC is both polygenic and extremely polymorphic, each individual possesses a personalized set of MHC molecules with different peptide-binding specificities, and collectively they present a unique individualized peptide imprint of the ongoing protein metabolism. Mapping all MHC allotypes is an enormous undertaking that cannot be achieved without a strong bioinformatics component. Computational tools for the prediction of peptide-MHC binding have thus become essential in most pipelines for T cell epitope discovery and an inescapable component of vaccine and cancer research. Here, we describe the development of several such tools, from pioneering efforts to the current state-of-the-art methods, that have allowed for accurate predictions of peptide binding of all MHC molecules, even including those that have not yet been characterized experimentally.

3CL hydrolase-based multiepitope peptide vaccine against SARS-CoV-2 using immunoinformatics

The present study provides the first multiepitope vaccine construct using the 3CL hydrolase protein of SARS-CoV-2. The coronavirus 3CL hydrolase (Mpro) enzyme is essential for proteolytic maturation of the virus. This study was based on immunoinformatics and structural vaccinology strategies. The design of the multiepitope vaccine was built using helper T-cell and cytotoxic T-cell epitopes from the 3CL hydrolase protein along with an adjuvant to enhance immune response; these are joined to each other by short peptide linkers. The vaccine also carries potential B-cell linear epitope regions, B-cell discontinuous epitopes, and interferon-γ-inducing epitopes. Epitopes of the constructed multiepitope vaccine were found to be antigenic, nonallergic, nontoxic, and covering large human populations worldwide. The vaccine construct was modeled, validated, and refined by different programs to achieve a high-quality three-dimensional structure. The resulting high-quality model was applied for conformational B-cell epitope selection and docking analyses with toll-like receptor-3 for understanding the capability of the vaccine to elicit an immune response. In silico cloning and codon adaptation were also performed with the pET-19b plasmid vector. The designed multiepitope peptide vaccine may prompt the development of a vaccine to control SARS-CoV-2 infection.

A comprehensive review and performance evaluation of bioinformatics tools for HLA class I peptide-binding prediction

Human leukocyte antigen class I (HLA-I) molecules are encoded by major histocompatibility complex (MHC) class I loci in humans. The binding and interaction between HLA-I molecules and intracellular peptides derived from a variety of proteolytic mechanisms play a crucial role in subsequent T-cell recognition of target cells and the specificity of the immune response. In this context, tools that predict the likelihood for a peptide to bind to specific HLA class I allotypes are important for selecting the most promising antigenic targets for immunotherapy. In this article, we comprehensively review a variety of currently available tools for predicting the binding of peptides to a selection of HLA-I allomorphs. Specifically, we compare their calculation methods for the prediction score, employed algorithms, evaluation strategies and software functionalities. In addition, we have evaluated the prediction performance of the reviewed tools based on an independent validation data set, containing 21 101 experimentally verified ligands across 19 HLA-I allotypes. The benchmarking results show that MixMHCpred 2.0.1 achieves the best performance for predicting peptides binding to most of the HLA-I allomorphs studied, while NetMHCpan 4.0 and NetMHCcons 1.1 outperform the other machine learning-based and consensus-based tools, respectively. Importantly, it should be noted that a peptide predicted with a higher binding score for a specific HLA allotype does not necessarily imply it will be immunogenic. That said, peptide-binding predictors are still very useful in that they can help to significantly reduce the large number of epitope candidates that need to be experimentally verified. Several other factors, including susceptibility to proteasome cleavage, peptide transport into the endoplasmic reticulum and T-cell receptor repertoire, also contribute to the immunogenicity of peptide antigens, and some of them can be considered by some predictors. Therefore, integrating features derived from these additional factors together with HLA-binding properties by using machine-learning algorithms may increase the prediction accuracy of immunogenic peptides. As such, we anticipate that this review and benchmarking survey will assist researchers in selecting appropriate prediction tools that best suit their purposes and provide useful guidelines for the development of improved antigen predictors in the future.

T Cell Epitope Predictions

Throughout the body, T cells monitor MHC-bound ligands expressed on the surface of essentially all cell types. MHC ligands that trigger a T cell immune response are referred to as T cell epitopes. Identifying such epitopes enables tracking, phenotyping, and stimulating T cells involved in immune responses in infectious disease, allergy, autoimmunity, transplantation, and cancer. The specific T cell epitopes recognized in an individual are determined by genetic factors such as the MHC molecules the individual expresses, in parallel to the individual's environmental exposure history. The complexity and importance of T cell epitope mapping have motivated the development of computational approaches that predict what T cell epitopes are likely to be recognized in a given individual or in a broader population. Such predictions guide experimental epitope mapping studies and enable computational analysis of the immunogenic potential of a given protein sequence region.

A comprehensive analysis of amino-peptidase N1 protein (APN) from Anopheles culicifacies for epitope design using Immuno-informatics models

Analysis of the Amino-peptidase N (APN) protein from Anopheles culicifacies as a vector based Transmission Blocking Vaccines (TBV) target has been considered for malaria vaccine development. Short peptides as potential epitopes for B cells and cytotoxic T cells and/or helper T cells were identified using prediction models provided by NetCTL and IEDB servers. Antigenicity determination, allergenicity, immunogenicity, epitope conservancy analysis, atomic interaction with HLA allele specific structure models and population coverage were investigated in this study. The analysis of the target protein helped to identify conserved regions as potential epitopes of APN in various Anopheles species. The T cell epitopes like peptides were further analyzed by using molecular docking to check interactions against the allele specific HLA models. Thus, we report the predicted B cell (VDERYRL) and T cell (RRYLATTQF for HLA class I and LKATFTVSI for HLA class II) epitopes like peptides from APN protein of Anopheles culicifacies (Diptera: Culicidae) for further consideration as vaccine candidates subsequent to in vitro and in vivo analysis.

Transfer String Kernel for Cross-Context DNA-Protein Binding Prediction

Through sequence-based classification, this paper tries to accurately predict the DNA binding sites of transcription factors (TFs) in an unannotated cellular context. Related methods in the literature fail to perform such predictions accurately, since they do not consider sample distribution shift of sequence segments from an annotated (source) context to an unannotated (target) context. We, therefore, propose a method called "Transfer String Kernel" (TSK) that achieves improved prediction of transcription factor binding site (TFBS) using knowledge transfer via cross-context sample adaptation. TSK maps sequence segments to a high-dimensional feature space using a discriminative mismatch string kernel framework. In this high-dimensional space, labeled examples of the source context are re-weighted so that the revised sample distribution matches the target context more closely. We have experimentally verified TSK for TFBS identifications on 14 different TFs under a cross-organism setting. We find that TSK consistently outperforms the state-of-the-art TFBS tools, especially when working with TFs whose binding sequences are not conserved across contexts. We also demonstrate the generalizability of TSK by showing its cutting-edge performance on a different set of cross-context tasks for the MHC peptide binding predictions.

Moving from Empirical to Rational Vaccine Design in the 'Omics' Era

An ideal vaccine provides long lasting protection against a pathogen by eliciting a well-rounded immune response which engages both innate and adaptive immunity. However, we have a limited understanding of how components of innate immunity, antibody and cell-mediated adaptive immunity interact and function together at a systems level. With advances in high-throughput ‘Omics’ methodologies it has become possible to capture global changes in the host, at a cellular and molecular level, that are induced by vaccination and infection. Analysis of these datasets has shown the promise of discovering mechanisms behind vaccine mediated protection, immunological memory, adverse effects as well as development of more efficient antigens and adjuvants. In this review, we will discuss how systems vaccinology takes advantage of new technology platforms and big data analysis, to enable the rational development of better vaccines.

Molecular Design of Peptide-Fc Fusion Drugs

Background:

Peptide-Fc fusion drugs, also known as peptibodies, are a category of biological therapeutics in which the Fc region of an antibody is genetically fused to a peptide of interest. However, to develop such kind of drugs is laborious and expensive. Rational design is urgently needed.

Methods:

We summarized the key steps in peptide-Fc fusion technology and stressed the main computational resources, tools, and methods that had been used in the rational design of peptide-Fc fusion drugs. We also raised open questions about the computer-aided molecular design of peptide-Fc.

Results:

The design of peptibody consists of four steps. First, identify peptide leads from native ligands, biopanning, and computational design or prediction. Second, select the proper Fc region from different classes or subclasses of immunoglobulin. Third, fuse the peptide leads and Fc together properly. At last, evaluate the immunogenicity of the constructs. At each step, there are quite a few useful resources and computational tools.

Conclusion:

Reviewing the molecular design of peptibody will certainly help make the transition from peptide leads to drugs on the market quicker and cheaper.

In silico prediction of nonpermissive HLA-DPB1 mismatches in unrelated HCT by functional distance

In silico prediction of high-risk donor-recipient HLA mismatches after unrelated donor (UD) hematopoietic cell transplantation (HCT) is an attractive, yet elusive, objective. Nonpermissive T-cell epitope (TCE) group mismatches were defined by alloreactive T-cell cross-reactivity for 52/80 HLA-DPB1 alleles (TCE-X). More recently, a numerical functional distance (FD) scoring system for in silico prediction of TCE groups based on the median impact of exon 2-encoded amino acid polymorphism on T-cell alloreactivity was developed for all DPB1 alleles (TCE-FD), including the 28/80 common alleles not assigned by TCE-X. We compared clinical outcome associations of nonpermissive DPB1 mismatches defined by TCE-X or TCE-FD in 8/8 HLA-matched UD-HCT for acute leukemia, myelodysplastic syndrome, and chronic myelogenous leukemia between 1999 and 2011 (N = 2730). Concordance between the 2 models was 92.3%, with most differences arising from DPB1*06:01 and DPB1*19:01 being differently assigned by TCE-X and TCE-FD. In both models, nonpermissive mismatches were associated with reduced overall survival (hazard ratio [HR], 1.15, P < .006 and HR, 1.12, P < .03), increased transplant-related mortality (HR, 1.31, P < .001 and HR, 1.26, P < .001) as well as acute (HR, 1.16, P < .02 and HR, 1.22, P < .001) and chronic (HR, 1.20, P < .003 and HR, 1.22, P < .001) graft-versus-host disease (GVHD). We show that in silico prediction of nonpermissive DPB1 mismatches significantly associated with major transplant outcomes is feasible for any DPB1 allele with known exon 2 sequence based on experimentally elaborated FD scores. This proof-of-principle observation opens new avenues for developing HLA risk-prediction models in HCT and has practical implications for UD searches.

DeepSeqPan, a novel deep convolutional neural network model for pan-specific class I HLA-peptide binding affinity prediction

Interactions between human leukocyte antigens (HLAs) and peptides play a critical role in the human immune system. Accurate computational prediction of HLA-binding peptides can be used for peptide drug discovery. Currently, the best prediction algorithms are neural network-based pan-specific models, which take advantage of the large amount of data across HLA alleles. However, current pan-specific models are all based on the pseudo sequence encoding for modeling the binding context, which is based on 34 positions identified from the HLA protein-peptide bound structures in early works. In this work, we proposed a novel deep convolutional neural network model (DCNN) for HLA-peptide binding prediction, in which the encoding of the HLA sequence and the binding context are both learned by the network itself without requiring the HLA-peptide bound structure information. Our DCNN model is also characterized by its binding context extraction layer and dual outputs with both binding affinity output and binding probability outputs. Evaluation on public benchmark datasets shows that our DeepSeqPan model without HLA structural information in training achieves state-of-the-art performance on a large number of HLA alleles with good generalization capability. Since our model only needs raw sequences from the HLA-peptide binding pairs, it can be applied to binding predictions of HLAs without structure information and can also be applied to other protein binding problems such as protein-DNA and protein-RNA bindings. The implementation code and trained models are freely available at https://github.com/pcpLiu/DeepSeqPan .

The Utility of Supertype Clustering in Prediction for Class II MHC-Peptide Binding

MOTIVATION

Extensive efforts have been devoted to understanding the antigenic peptides binding to MHC class I and II molecules since they play a fundamental role in controlling immune responses and due their involvement in vaccination, transplantation, and autoimmunity. The genes coding for the MHC molecules are highly polymorphic, and it is difficult to build computational models for MHC molecules with few know binders. On the other hand, previous studies demonstrated that some MHC molecules share overlapping peptide binding repertoires and attempted to group them into supertypes. Herein, we present a framework of the utility of supertype clustering to gain more information about the data to improve the prediction accuracy of class II MHC-peptide binding.

RESULTS

We developed a new method, called superMHC, for class II MHC-peptide binding prediction, including three MHC isotypes of HLA-DR, HLA-DP, and HLA-DQ, by using supertype clustering in conjunction with RLS regression. The supertypes were identified by using a novel repertoire dissimilarity index to quantify the difference in MHC binding specificities. The superMHC method achieves the state-of-the-art performance and is demonstrated to predict binding affinities to a series of MHC molecules with few binders accurately. These results have implications for understanding receptor-ligand interactions involved in MHC-peptide binding.

Improved methods for predicting peptide binding affinity to MHC class II molecules

Major histocompatibility complex class II (MHC-II) molecules are expressed on the surface of professional antigen-presenting cells where they display peptides to T helper cells, which orchestrate the onset and outcome of many host immune responses. Understanding which peptides will be presented by the MHC-II molecule is therefore important for understanding the activation of T helper cells and can be used to identify T-cell epitopes. We here present updated versions of two MHC-II-peptide binding affinity prediction methods, NetMHCII and NetMHCIIpan. These were constructed using an extended data set of quantitative MHC-peptide binding affinity data obtained from the Immune Epitope Database covering HLA-DR, HLA-DQ, HLA-DP and H-2 mouse molecules. We show that training with this extended data set improved the performance for peptide binding predictions for both methods. Both methods are publicly available at www.cbs.dtu.dk/services/NetMHCII-2.3 and www.cbs.dtu.dk/services/NetMHCIIpan-3.2.

Trans-Allelic Model for Prediction of Peptide:MHC-II Interactions

Major histocompatibility complex class two (MHC-II) molecules are trans-membrane proteins and key components of the cellular immune system. Upon recognition of foreign peptides expressed on the MHC-II binding groove, CD4⁺ T cells mount an immune response against invading pathogens. Therefore, mechanistic identification and knowledge of physicochemical features that govern interactions between peptides and MHC-II molecules is useful for the design of effective epitope-based vaccines, as well as for understanding of immune responses. In this article, we present a comprehensive trans-allelic prediction model, a generalized version of our previous biophysical model, that can predict peptide interactions for all three human MHC-II loci (HLA-DR, HLA-DP, and HLA-DQ), using both peptide sequence data and structural information of MHC-II molecules. The advantage of this approach over other machine learning models is that it offers a simple and plausible physical explanation for peptide–MHC-II interactions. We train the model using a benchmark experimental dataset and measure its predictive performance using novel data. Despite its relative simplicity, we find that the model has comparable performance to the state-of-the-art method, the NetMHCIIpan method. Focusing on the physical basis of peptide–MHC binding, we find support for previous theoretical predictions about the contributions of certain binding pockets to the binding energy. In addition, we find that binding pocket P5 of HLA-DP, which was not previously considered as a primary anchor, does make strong contribution to the binding energy. Together, the results indicate that our model can serve as a useful complement to alternative approaches to predicting peptide–MHC interactions.