Derivation of an amino acid similarity matrix for peptide: MHC binding and its application as a Bayesian prior

RLpMIEC: High-Affinity Peptide Generation Targeting Major Histocompatibility Complex-I Guided and Interpreted by Interaction Spectrum-Navigated Reinforcement Learning

Major histocompatibility complex (MHC) plays a vital role in presenting epitopes (short peptides from pathogenic proteins) to T-cell receptors (TCRs) to trigger the subsequent immune responses. Vaccine design targeting MHC generally aims to find epitopes with a high binding affinity for MHC presentation. Nevertheless, to find novel epitopes usually requires high-throughput screening of bulk peptide database, which is time-consuming, labor-intensive, more unaffordable, and very expensive. Excitingly, the past several years have witnessed the great success of artificial intelligence (AI) in various fields, such as natural language processing (NLP, e.g., GPT-4), protein structure prediction and engineering (e.g., AlphaFold2), and so on. Therefore, herein, we propose a deep reinforcement-learning (RL)-based generative algorithm, RLpMIEC, to quantitatively design peptide targeting MHC-I systems. Specifically, RLpMIEC combines the energetic spectrum (namely, the molecular interaction energy component, MIEC) based on the peptide-MHC interaction and the sequence information to generate peptides with strong binding affinity and precise MIEC spectra to accelerate the discovery of candidate peptide vaccines. RLpMIEC performs well in all the generative capability evaluations and can generate peptides with strong binding affinities and precise MIECs and, moreover, with high interpretability, demonstrating its powerful capability in participation for accelerating peptide-based vaccine development.

MHC-I-presented non-canonical antigens expand the cancer immunotherapy targets in acute myeloid leukemia

Identification of tumor neoantigens is indispensable for the development of cancer immunotherapies. However, we are still lacking knowledge about the potential neoantigens derived from sequences outside protein-coding regions. Here, we comprehensively characterized the immunopeptidome landscape by integrating multi-omics data in acute myeloid leukemia (AML). Both canonical and non-canonical MHC-associated peptides (MAPs) in AML were identified. We found that the quality and characteristics of ncMAPs are comparable or superior to cMAPs, suggesting ncMAPs are indispensable sources for tumor neoantigens. We further proposed a computational framework to prioritize the neoantigens by integrating additional transcriptome and immunopeptidome in normal tissues. Notably, 6 of prioritized 13 neoantigens were derived from ncMAPs. The expressions of corresponding source genes are highly related to infiltrations of immune cells. Finally, a risk model was developed, which exhibited good performance for clinical prognosis in AML. Our findings expand potential cancer immunotherapy targets and provide in-depth insights into AML treatment, laying a new foundation for precision therapies in AML.

APEX-pHLA: A novel method for accurate prediction of the binding between exogenous short peptides and HLA class I molecules

Human leukocyte antigen (HLA) molecules play critically significant role within the realm of immunotherapy due to their capacities to recognize and bind exogenous antigens such as peptides, subsequently delivering them to immune cells. Predicting the binding between peptides and HLA molecules (pHLA) can expedite the screening of immunogenic peptides and facilitate vaccine design. However, traditional experimental methods are time-consuming and inefficient. In this study, an efficient method based on deep learning was developed for predicting peptide-HLA binding, which treated peptide sequences as linguistic entities. It combined the architectures of textCNN and BiLSTM to create a deep neural network model called APEX-pHLA. This model operated without limitations related to HLA class I allele variants and peptide segment lengths, enabling efficient encoding of sequence features for both HLA and peptide segments. On the independent test set, the model achieved Accuracy, ROC_AUC, F1, and MCC is 0.9449, 0.9850, 0.9453, and 0.8899, respectively. Similarly, on an external test set, the results were 0.9803, 0.9574, 0.8835, and 0.7863, respectively. These findings outperformed fifteen methods previously reported in the literature. The accurate prediction capability of the APEX-pHLA model in peptide-HLA binding might provide valuable insights for future HLA vaccine design.

How antisense transcripts can evolve to encode novel proteins

Protein coding features can emerge de novo in non coding transcripts, resulting in emergence of new protein coding genes. Studies across many species show that a large fraction of evolutionarily novel non-coding RNAs have an antisense overlap with protein coding genes. The open reading frames (ORFs) in these antisense RNAs could also overlap with existing ORFs. In this study, we investigate how the evolution an ORF could be constrained by its overlap with an existing ORF in three different reading frames. Using a combination of mathematical modeling and genome/transcriptome data analysis in two different model organisms, we show that antisense overlap can increase the likelihood of ORF emergence and reduce the likelihood of ORF loss, especially in one of the three reading frames. In addition to rationalising the repeatedly reported prevalence of de novo emerged genes in antisense transcripts, our work also provides a generic modeling and an analytical framework that can be used to understand evolution of antisense genes.

GIHP: Graph convolutional neural network based interpretable pan-specific HLA-peptide binding affinity prediction

Accurately predicting the binding affinities between Human Leukocyte Antigen (HLA) molecules and peptides is a crucial step in understanding the adaptive immune response. This knowledge can have important implications for the development of effective vaccines and the design of targeted immunotherapies. Existing sequence-based methods are insufficient to capture the structure information. Besides, the current methods lack model interpretability, which hinder revealing the key binding amino acids between the two molecules. To address these limitations, we proposed an interpretable graph convolutional neural network (GCNN) based prediction method named GIHP. Considering the size differences between HLA and short peptides, GIHP represent HLA structure as amino acid-level graph while represent peptide SMILE string as atom-level graph. For interpretation, we design a novel visual explanation method, gradient weighted activation mapping (Grad-WAM), for identifying key binding residues. GIHP achieved better prediction accuracy than state-of-the-art methods across various datasets. According to current research findings, key HLA-peptide binding residues mutations directly impact immunotherapy efficacy. Therefore, we verified those highlighted key residues to see whether they can significantly distinguish immunotherapy patient groups. We have verified that the identified functional residues can successfully separate patient survival groups across breast, bladder, and pan-cancer datasets. Results demonstrate that GIHP improves the accuracy and interpretation capabilities of HLA-peptide prediction, and the findings of this study can be used to guide personalized cancer immunotherapy treatment. Codes and datasets are publicly accessible at: https://github.com/sdustSu/GIHP.

HLAIImaster: a deep learning method with adaptive domain knowledge predicts HLA II neoepitope immunogenic responses

While significant strides have been made in predicting neoepitopes that trigger autologous CD4+ T cell responses, accurately identifying the antigen presentation by human leukocyte antigen (HLA) class II molecules remains a challenge. This identification is critical for developing vaccines and cancer immunotherapies. Current prediction methods are limited, primarily due to a lack of high-quality training epitope datasets and algorithmic constraints. To predict the exogenous HLA class II-restricted peptides across most of the human population, we utilized the mass spectrometry data to profile >223 000 eluted ligands over HLA-DR, -DQ, and -DP alleles. Here, by integrating these data with peptide processing and gene expression, we introduce HLAIImaster, an attention-based deep learning framework with adaptive domain knowledge for predicting neoepitope immunogenicity. Leveraging diverse biological characteristics and our enhanced deep learning framework, HLAIImaster is significantly improved against existing tools in terms of positive predictive value across various neoantigen studies. Robust domain knowledge learning accurately identifies neoepitope immunogenicity, bridging the gap between neoantigen biology and the clinical setting and paving the way for future neoantigen-based therapies to provide greater clinical benefit. In summary, we present a comprehensive exploitation of the immunogenic neoepitope repertoire of cancers, facilitating the effective development of "just-in-time" personalized vaccines.

Comprehensive profiling of cancer neoantigens from aberrant RNA splicing

BACKGROUND

Cancer neoantigens arise from protein-altering somatic mutations in tumor and rank among the most promising next-generation immuno-oncology agents when used in combination with immune checkpoint inhibitors. We previously developed a computational framework, REAL-neo, for identification, quality control, and prioritization of both class-I and class-II human leucocyte antigen (HLA)-presented neoantigens resulting from somatic single-nucleotide mutations, small insertions and deletions, and gene fusions. In this study, we developed a new module, SPLICE-neo, to identify neoantigens from aberrant RNA transcripts from two distinct sources: (1) DNA mutations within splice sites and (2) de novo RNA aberrant splicings.

METHODS

First, SPLICE-neo was used to profile all DNA splice-site mutations in 11,892 tumors from The Cancer Genome Atlas (TCGA) and identified 11 profiles of splicing donor or acceptor site gains or losses. Transcript isoforms resulting from the top seven most frequent profiles were computed using novel logic models. Second, SPLICE-neo identified de novo RNA splicing events using RNA sequencing reads mapped to novel exon junctions from either single, double, or multiple exon-skipping events. The aberrant transcripts from both sources were then ranked based on isoform expression levels and z-scores assuming that individual aberrant splicing events are rare. Finally, top-ranked novel isoforms were translated into protein, and the resulting neoepitopes were evaluated for neoantigen potential using REAL-neo. The top splicing neoantigen candidates binding to HLA-A*02:01 were validated using in vitro T2 binding assays.

RESULTS

We identified abundant splicing neoantigens in four representative TCGA cancers: BRCA, LUAD, LUSC, and LIHC. In addition to their substantial contribution to neoantigen load, several splicing neoantigens were potent tumor antigens with stronger bindings to HLA compared with the positive control of antigens from influenza virus.

CONCLUSIONS

SPLICE-neo is the first tool to comprehensively identify and prioritize splicing neoantigens from both DNA splice-site mutations and de novo RNA aberrant splicings. There are two major advances of SPLICE-neo. First, we developed novel logic models that assemble and prioritize full-length aberrant transcripts from DNA splice-site mutations. Second, SPLICE-neo can identify exon-skipping events involving more than two exons, which account for a quarter to one-third of all skipping events.

In Silico Analysis of a Candidate Multi-epitope Peptide Vaccine Against Human Brucellosis

Brucellosis is one of the neglected endemic zoonoses in the world. Vaccination appears to be a promising health strategy to prevent it. This study used advanced computational techniques to develop a potent multi-epitope vaccine for human brucellosis. Seven epitopes from four main brucella species that infect humans were selected. They had significant potential to induce cellular and humoral responses. They showed high antigenic ability without the allergenic characteristic. In order to improve its immunogenicity, suitable adjuvants were also added to the structure of the vaccine. The physicochemical and immunological properties of the vaccine were evaluated. Then its two and three-dimensional structure was predicted. The vaccine was docked with toll-like receptor4 to assess its ability to stimulate innate immune responses. For successful expression of the vaccine protein in Escherichia coli, in silico cloning, codon optimization, and mRNA stability were evaluated. The immune simulation was performed to reveal the immune response profile of the vaccine after injection. The designed vaccine showed the high ability to induce immune response, especially cellular responses to human brucellosis. It showed the appropriate physicochemical properties, a high-quality structure, and a high potential for expression in a prokaryotic system.

TripHLApan: predicting HLA molecules binding peptides based on triple coding matrix and transfer learning

Human leukocyte antigen (HLA) recognizes foreign threats and triggers immune responses by presenting peptides to T cells. Computationally modeling the binding patterns between peptide and HLA is very important for the development of tumor vaccines. However, it is still a big challenge to accurately predict HLA molecules binding peptides. In this paper, we develop a new model TripHLApan for predicting HLA molecules binding peptides by integrating triple coding matrix, BiGRU + Attention models, and transfer learning strategy. We have found the main interaction site regions between HLA molecules and peptides, as well as the correlation between HLA encoding and binding motifs. Based on the discovery, we make the preprocessing and coding closer to the natural biological process. Besides, due to the input being based on multiple types of features and the attention module focused on the BiGRU hidden layer, TripHLApan has learned more sequence level binding information. The application of transfer learning strategies ensures the accuracy of prediction results under special lengths (peptides in length 8) and model scalability with the data explosion. Compared with the current optimal models, TripHLApan exhibits strong predictive performance in various prediction environments with different positive and negative sample ratios. In addition, we validate the superiority and scalability of TripHLApan's predictive performance using additional latest data sets, ablation experiments and binding reconstitution ability in the samples of a melanoma patient. The results show that TripHLApan is a powerful tool for predicting the binding of HLA-I and HLA-II molecular peptides for the synthesis of tumor vaccines. TripHLApan is publicly available at https://github.com/CSUBioGroup/TripHLApan.git.

APE-Gen2.0: Expanding Rapid Class I Peptide-Major Histocompatibility Complex Modeling to Post-Translational Modifications and Noncanonical Peptide Geometries

The recognition of peptides bound to class I major histocompatibility complex (MHC-I) receptors by T-cell receptors (TCRs) is a determinant of triggering the adaptive immune response. While the exact molecular features that drive the TCR recognition are still unknown, studies have suggested that the geometry of the joint peptide-MHC (pMHC) structure plays an important role. As such, there is a definite need for methods and tools that accurately predict the structure of the peptide bound to the MHC-I receptor. In the past few years, many pMHC structural modeling tools have emerged that provide high-quality modeled structures in the general case. However, there are numerous instances of non-canonical cases in the immunopeptidome that the majority of pMHC modeling tools do not attend to, most notably, peptides that exhibit non-standard amino acids and post-translational modifications (PTMs) or peptides that assume non-canonical geometries in the MHC binding cleft. Such chemical and structural properties have been shown to be present in neoantigens; therefore, accurate structural modeling of these instances can be vital for cancer immunotherapy. To this end, we have developed APE-Gen2.0, a tool that improves upon its predecessor and other pMHC modeling tools, both in terms of modeling accuracy and the available modeling range of non-canonical peptide cases. Some of the improvements include (i) the ability to model peptides that have different types of PTMs such as phosphorylation, nitration, and citrullination; (ii) a new and improved anchor identification routine in order to identify and model peptides that exhibit a non-canonical anchor conformation; and (iii) a web server that provides a platform for easy and accessible pMHC modeling. We further show that structures predicted by APE-Gen2.0 can be used to assess the effects that PTMs have in binding affinity in a more accurate manner than just using solely the sequence of the peptide. APE-Gen2.0 is freely available at https://apegen.kavrakilab.org.

PeptiVAX: A new adaptable peptides-delivery platform for development of CTL-based, SARS-CoV-2 vaccines

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) posed a threat to public health and the global economy, necessitating the development of various vaccination strategies. Mutations in the SPIKE protein gene, a crucial component of mRNA and adenovirus-based vaccines, raised concerns about vaccine efficacy, prompting the need for rapid vaccine updates. To address this, we leveraged PeptiCRAd, an oncolytic vaccine based on tumor antigen decorated oncolytic adenoviruses, creating a vaccine platform called PeptiVAX. First, we identified multiple CD8 T-cell epitopes from highly conserved regions across coronaviruses, expanding the range of T-cell responses to non-SPIKE proteins. We designed short segments containing the predicted epitopes presented by common HLA-Is in the global population. Testing the immunogenicity, we characterized T-cell responses to candidate peptides in peripheral blood mononuclear cells (PBMCs) from pre-pandemic healthy donors and ICU patients. As a proof of concept in mice, we selected a peptide with epitopes predicted to bind to murine MHC-I haplotypes. Our technology successfully elicited peptide-specific T-cell responses, unaffected by the use of unarmed adenoviral vectors or adeno-based vaccines encoding SPIKE. In conclusion, PeptiVAX represents a fast and adaptable SARS-CoV-2 vaccine delivery system that broadens T-cell responses beyond the SPIKE protein, offering potential benefits for vaccine effectiveness.

Stimulating T cell responses against patient-derived breast cancer cells with neoantigen peptide-loaded peripheral blood mononuclear cells

Breast cancer stands as a formidable global health challenge for women. While neoantigens exhibit efficacy in activating T cells specific to cancer and instigating anti-tumor immune responses, the accuracy of neoantigen prediction remains suboptimal. In this study, we identified neoantigens from the patient-derived breast cancer cells, PC-B-142CA and PC-B-148CA cells, utilizing whole-genome and RNA sequencing. The pVAC-Seq pipeline was employed, with minor modification incorporating criteria (1) binding affinity of mutant (MT) peptide with HLA (IC50 MT) ≤ 500 nm in 3 of 5 algorithms and (2) IC50 wild type (WT)/MT > 1. Sequencing results unveiled 2513 and 3490 somatic mutations, and 646 and 652 non-synonymous mutations in PC-B-142CA and PC-B-148CA, respectively. We selected the top 3 neoantigens to perform molecular dynamic simulation and synthesized 9-12 amino acid neoantigen peptides, which were then pulsed onto healthy donor peripheral blood mononuclear cells (PBMCs). Results demonstrated that T cells activated by ADGRL1E274K, PARP1E619K, and SEC14L2R43Q peptides identified from PC-B-142CA exhibited significantly increased production of interferon-gamma (IFN-γ), while PARP1E619K and SEC14L2R43Q peptides induced the expression of CD107a on T cells. The % tumor cell lysis was notably enhanced by T cells activated with MT peptides across all three healthy donors. Moreover, ALKBH6V83M and GAAI823T peptides from PC-B-148CA remarkably stimulated IFN-γ- and CD107a-positive T cells, displaying high cell-killing activity against target cancer cells. In summary, our findings underscore the successful identification of neoantigens with anti-tumor T cell functions and highlight the potential of personalized neoantigens as a promising avenue for breast cancer treatment.

Comparison of the somatic genomic landscape between central- and peripheral-type non-small cell lung cancer

OBJECTIVE

Lung cancer is classified into central and peripheral types based on the anatomic location. The present study aimed to explore the distinct patterns of genomic alterations between central- and peripheral-type non-small cell lung cancers (NSCLCs) with negative driver genes and identify potential driver genes and biomarkers to improve therapy strategies for NSCLC.

METHODS

Whole-exome sequencing (WES) was performed with 182 tumor/control pairs of samples from 145 Chinese NSCLC patients without EGFR, ALK, or ROS1 alterations. Significantly mutated genes (SMGs) and somatic copy number alterations (SCNAs) were identified. Subsequently, tumor mutation burden (TMB), weighted genome integrity index (wGII), copy number alteration (CNA) burden, Shannon diversity index (SDI), intratumor heterogeneity (ITH), neoantigen load (NAL), and clonal variations were evaluated in central- and peripheral-type NSCLCs. Furthermore, mutational signature analysis and survival analysis were performed.

RESULTS

TP53 was the most frequently mutated gene in NSCLC and more frequently mutated in central-type NSCLC. Higher wGII, ITH, and SDI were found in central-type lung adenocarcinoma (LUAD) than in peripheral-type LUAD. The NAL of central-type lung squamous cell carcinoma (LUSC) with stage III/IV was significantly higher than that of peripheral-type LUSC. Mutational signature analysis revealed that SBS10b, SBS24, and ID7 were significantly different in central- and peripheral-type NSCLCs. Furthermore, central-type NSCLC was found to evolve at a higher level with fewer clones and more subclones, particularly in central-type LUSC. Survival analysis revealed that TMB, CNA burden, NAL, subclonal driver mutations, and subclonal mutations were negatively related to the overall survival (OS) and the progression-free survival (PFS) of central-type LUAD.

CONCLUSIONS

Central-type NSCLC tended to evolve at a higher level and might suggest a favorable response to immunotherapy. Our study also identified several potential driver genes and promising biomarkers for the prognosis and prediction of chemotherapy responses in NSCLC.

Accurate TCR-pMHC interaction prediction using a BERT-based transfer learning method

Accurate prediction of TCR-pMHC binding is important for the development of cancer immunotherapies, especially TCR-based agents. Existing algorithms often experience diminished performance when dealing with unseen epitopes, primarily due to the complexity in TCR-pMHC recognition patterns and the scarcity of available data for training. We have developed a novel deep learning model, 'TCR Antigen Binding Recognition' based on BERT, named as TABR-BERT. Leveraging BERT's potent representation learning capabilities, TABR-BERT effectively captures essential information regarding TCR-pMHC interactions from TCR sequences, antigen epitope sequences and epitope-MHC binding. By transferring this knowledge to predict TCR-pMHC recognition, TABR-BERT demonstrated better results in benchmark tests than existing methods, particularly for unseen epitopes.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptides bound to the human MHC, HLA, has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within our curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer pHLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work may be applied towards predicting antigen immunogenicity, and receptor cross-reactivity.

Mismatch repair deficiency is not sufficient to elicit tumor immunogenicity

DNA mismatch repair deficiency (MMRd) is associated with a high tumor mutational burden (TMB) and sensitivity to immune checkpoint blockade (ICB) therapy. Nevertheless, most MMRd tumors do not durably respond to ICB and critical questions remain about immunosurveillance and TMB in these tumors. In the present study, we developed autochthonous mouse models of MMRd lung and colon cancer. Surprisingly, these models did not display increased T cell infiltration or ICB response, which we showed to be the result of substantial intratumor heterogeneity of mutations. Furthermore, we found that immunosurveillance shapes the clonal architecture but not the overall burden of neoantigens, and T cell responses against subclonal neoantigens are blunted. Finally, we showed that clonal, but not subclonal, neoantigen burden predicts ICB response in clinical trials of MMRd gastric and colorectal cancer. These results provide important context for understanding immune evasion in cancers with a high TMB and have major implications for therapies aimed at increasing TMB.

Deep neural networks predict class I major histocompatibility complex epitope presentation and transfer learn neoepitope immunogenicity

Identifying neoepitopes that elicit an adaptive immune response is a major bottleneck to developing personalized cancer vaccines. Experimental validation of candidate neoepitopes is extremely resource intensive and the vast majority of candidates are non-immunogenic, creating a needle-in-a-haystack problem. Here we address this challenge, presenting computational methods for predicting class I major histocompatibility complex (MHC-I) epitopes and identifying immunogenic neoepitopes with improved precision. The BigMHC method comprises an ensemble of seven pan-allelic deep neural networks trained on peptide-MHC eluted ligand data from mass spectrometry assays and transfer learned on data from assays of antigen-specific immune response. Compared with four state-of-the-art classifiers, BigMHC significantly improves the prediction of epitope presentation on a test set of 45,409 MHC ligands among 900,592 random negatives (area under the receiver operating characteristic = 0.9733; area under the precision-recall curve = 0.8779). After transfer learning on immunogenicity data, BigMHC yields significantly higher precision than seven state-of-the-art models in identifying immunogenic neoepitopes, making BigMHC effective in clinical settings.

Modeling of MT. P495, an mRNA-based vaccine against the phosphate-binding protein PstS1 of Mycobacterium tuberculosis

Tuberculosis (TB) is a contagious disease that predominantly affects the lungs, but can also spread to other organs via the bloodstream. TB affects about one-fourth population of the world. With age, the effectiveness of Bacillus Calmette-Guérin (BCG), the only authorized TB vaccine, decreases. In the quest for a prophylactic and immunotherapeutic vaccine, in this study, a hypothetical mRNA vaccine is delineated, named MT. P495, implementing in silico and immunoinformatics approaches to evaluate key aspects and immunogenic epitopes across the PstS1, a highly conserved periplasmic protein of Mycobacterium tuberculosis (Mtb). PstS1 elicited the potential to generate 99.9% population coverage worldwide. The presence of T- and B-cell epitopes across the PstS1 protein were validated using several computational prediction tools. Molecular docking and dynamics simulation confirmed stable epitope-allele interaction. Immune cell response to the antigen clearance rate was verified by the in silico analysis of immune simulation. Codon optimization confirmed the efficient translation of the mRNA in the host cell. With Toll-like receptors, the vaccine exhibited stable and strong interactions. Findings suggest that the MT. P495 vaccine probably will elicit specific immune responses against Mtb. This mRNA vaccine model is a ready source for further wet-lab validation to confirm the efficacy of this proposed vaccine candidate.

High-throughput, targeted MHC class I immunopeptidomics using a functional genetics screening platform

Identification of CD8+ T cell epitopes is critical for the development of immunotherapeutics. Existing methods for major histocompatibility complex class I (MHC class I) ligand discovery are time intensive, specialized and unable to interrogate specific proteins on a large scale. Here, we present EpiScan, which uses surface MHC class I levels as a readout for whether a genetically encoded peptide is an MHC class I ligand. Predetermined starting pools composed of >100,000 peptides can be designed using oligonucleotide synthesis, permitting large-scale MHC class I screening. We exploit this programmability of EpiScan to uncover an unappreciated role for cysteine that increases the number of predicted ligands by 9-21%, reveal affinity hierarchies by analysis of biased anchor peptide libraries and screen viral proteomes for MHC class I ligands. Using these data, we generate and iteratively refine peptide binding predictions to create EpiScan Predictor. EpiScan Predictor performs comparably to other state-of-the-art MHC class I peptide binding prediction algorithms without suffering from underrepresentation of cysteine-containing peptides. Thus, targeted immunopeptidomics using EpiScan will accelerate CD8+ T cell epitope discovery toward the goal of individual-specific immunotherapeutics.

Integrated Clinical, Molecular and Immunological Characterization of Pulmonary Sarcomatoid Carcinomas Reveals an Immune Escape Mechanism That May Influence Therapeutic Strategies

Pulmonary sarcomatoid carcinoma (PSC) has highly aggressive biological behaviour and poor clinical outcomes, raising expectations for new therapeutic strategies. We characterized 179 PSC by immunohistochemistry, next-generation sequencing and in silico analysis using a deep learning algorithm with respect to clinical, immunological and molecular features. PSC was more common in men, older ages and smokers. Surgery was an independent factor (p < 0.01) of overall survival (OS). PD-L1 expression was detected in 82.1% of all patients. PSC patients displaying altered epitopes due to processing mutations showed another PD-L1-independent immune escape mechanism, which also significantly influenced OS (p < 0.02). The effect was also maintained when only advanced tumour stages were considered (p < 0.01). These patients also showed improved survival with a significant correlation for immunotherapy (p < 0.05) when few or no processing mutations were detected, although this should be interpreted with caution due to the small number of patients studied. Genomic alterations for which there are already approved drugs were present in 35.4% of patients. Met exon 14 skipping was found more frequently (13.7%) and EGFR mutations less frequently (1.7%) than in other NSCLC. In summary, in addition to the divergent genomic landscape of PSC, the specific immunological features of this prognostically poor subtype should be considered in therapy stratification.

Entropic analysis of antigen-specific CDR3 domains identifies essential binding motifs shared by CDR3s with different antigen specificities

Antigen-specific T cell receptor (TCR) sequences can have prognostic, predictive, and therapeutic value, but decoding the specificity of TCR recognition remains challenging. Unlike DNA strands that base pair, TCRs bind to their targets with different orientations and different lengths, which complicates comparisons. We present scanning parametrized by normalized TCR length (SPAN-TCR) to analyze antigen-specific TCR CDR3 sequences and identify patterns driving TCR-pMHC specificity. Using entropic analysis, SPAN-TCR identifies 2-mer motifs that decrease the diversity (entropy) of CDR3s. These motifs are the most common patterns that can predict CDR3 composition, and we identify "essential" motifs that decrease entropy in the same CDR3 α or β chain containing the 2-mer, and "super-essential" motifs that decrease entropy in both chains. Molecular dynamics analysis further suggests that these motifs may play important roles in binding. We then employ SPAN-TCR to resolve similarities in TCR repertoires against different antigens using public databases of TCR sequences.

Computational prediction of MHC anchor locations guides neoantigen identification and prioritization

Neoantigens are tumor-specific peptide sequences resulting from sources such as somatic DNA mutations. Upon loading onto major histocompatibility complex (MHC) molecules, they can trigger recognition by T cells. Accurate neoantigen identification is thus critical for both designing cancer vaccines and predicting response to immunotherapies. Neoantigen identification and prioritization relies on correctly predicting whether the presenting peptide sequence can successfully induce an immune response. Because most somatic mutations are single-nucleotide variants, changes between wild-type and mutated peptides are typically subtle and require cautious interpretation. A potentially underappreciated variable in neoantigen prediction pipelines is the mutation position within the peptide relative to its anchor positions for the patient's specific MHC molecules. Whereas a subset of peptide positions are presented to the T cell receptor for recognition, others are responsible for anchoring to the MHC, making these positional considerations critical for predicting T cell responses. We computationally predicted anchor positions for different peptide lengths for 328 common HLA alleles and identified unique anchoring patterns among them. Analysis of 923 tumor samples shows that 6 to 38% of neoantigen candidates are potentially misclassified and can be rescued using allele-specific knowledge of anchor positions. A subset of anchor results were orthogonally validated using protein crystallography structures. Representative anchor trends were experimentally validated using peptide-MHC stability assays and competition binding assays. By incorporating our anchor prediction results into neoantigen prediction pipelines, we hope to formalize, streamline, and improve the identification process for relevant clinical studies.

Recent trends in next generation immunoinformatics harnessed for universal coronavirus vaccine design

The ongoing pandemic of coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has globally devastated public health, the economies of many countries and quality of life universally. The recent emergence of immune-escaped variants and scenario of vaccinated individuals being infected has raised the global concerns about the effectiveness of the current available vaccines in transmission control and disease prevention. Given the high rate mutation of SARS-CoV-2, an efficacious vaccine targeting against multiple variants that contains virus-specific epitopes is desperately needed. An immunoinformatics approach is gaining traction in vaccine design and development due to the significant reduction in time and cost of immunogenicity studies and increasing reliability of the generated results. It can underpin the development of novel therapeutic methods and accelerate the design and production of peptide vaccines for infectious diseases. Structural proteins, particularly spike protein (S), along with other proteins have been studied intensively as promising coronavirus vaccine targets. Numbers of promising online immunological databases, tools and web servers have widely been employed for the design and development of next generation COVID-19 vaccines. This review highlights the role of immunoinformatics in identifying immunogenic peptides as potential vaccine targets, involving databases, and prediction and characterization of epitopes which can be harnessed for designing future coronavirus vaccines.

Challenges in neoantigen-directed therapeutics

A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.

A Sample Guideline for Reverse Vaccinology Approach for the Development of Subunit Vaccine Using Varicella Zoster as a Model Disease

For the development of multi-peptide vaccine, identification of antigenic epitopes is crucial. If it is done using wet lab techniques, the identification process can be time-consuming, laborious, and cost-intensive. In silico tools, on the other hand, enable researchers to predict potential epitopes with little to no cost for further in vivo and in vitro testing. The rapid identification process using in silico tools helps in responding to health emergencies faster. Developing an efficient and high coverage vaccine is one of the ways to reduce morbidity and mortality rates of the diseases and protect the affected populations. In this chapter, we introduce the necessary tools and methodology for the identification and characterization of antigenic epitopes to design a multi-epitope vaccine using varicella-zoster virus as an example vector model.

HLA Class I Supertype Classification Based on Structural Similarity

HLA class I proteins, a critical component in adaptive immunity, bind and present intracellular Ags to CD8+ T cells. The extreme polymorphism of HLA genes and associated peptide binding specificities leads to challenges in various endeavors, including neoantigen vaccine development, disease association studies, and HLA typing. Supertype classification, defined by clustering functionally similar HLA alleles, has proven helpful in reducing the complexity of distinguishing alleles. However, determining supertypes via experiments is impractical, and current in silico classification methods exhibit limitations in stability and functional relevance. In this study, by incorporating three-dimensional structures we present a method for classifying HLA class I molecules with improved breadth, accuracy, stability, and flexibility. Critical for these advances is our finding that structural similarity highly correlates with peptide binding specificity. The new classification should be broadly useful in peptide-based vaccine development and HLA-disease association studies.

Biophysical and biochemical insights in the design of immunoassays

BACKGROUND

Rapid antigen assays have been attractive for decentralized, point of care diagnostics because of their low cost, robustness, and ease of use. The development of a diagnostic assay for a newly emerging infectious disease needs to take into account the progression of a disease, whether there is human to human transmission, and patient biomarker levels with time, and these all impact the choice of antigen targets and affinity agents.

SCOPE OF REVIEW

The factors involved in the biophysical design of rapid antigen immunoassays are discussed, focusing on antigen selection and designing for cross-reactivity. State of the art in the biophysical characterization of protein-ligand or antigen-antibody interactions, the different types of affinity agents used in immunoassays, and biochemical conjugation strategies are described.

MAJOR CONCLUSIONS

Antigen choice is a critical factor in immunoassay diagnostic development, and should account for the properties of the virion, virus, and disease progression. Biophysical and biochemical aspects of immunoassays are critical for performance.

GENERAL SIGNIFICANCE

This review can serve as an instructive guide to aid in diagnostic development for future emerging diseases.

Tumor microenvironment: barrier or opportunity towards effective cancer therapy

Tumor microenvironment (TME) is a specialized ecosystem of host components, designed by tumor cells for successful development and metastasis of tumor. With the advent of 3D culture and advanced bioinformatic methodologies, it is now possible to study TME’s individual components and their interplay at higher resolution. Deeper understanding of the immune cell’s diversity, stromal constituents, repertoire profiling, neoantigen prediction of TMEs has provided the opportunity to explore the spatial and temporal regulation of immune therapeutic interventions. The variation of TME composition among patients plays an important role in determining responders and non-responders towards cancer immunotherapy. Therefore, there could be a possibility of reprogramming of TME components to overcome the widely prevailing issue of immunotherapeutic resistance. The focus of the present review is to understand the complexity of TME and comprehending future perspective of its components as potential therapeutic targets. The later part of the review describes the sophisticated 3D models emerging as valuable means to study TME components and an extensive account of advanced bioinformatic tools to profile TME components and predict neoantigens. Overall, this review provides a comprehensive account of the current knowledge available to target TME.

Not all T cell epitopes are equally desired: a review of in silico tools for the prediction of cytokine-inducing potential of T-cell epitopes

Assessment of protective or harmful T cell response induced by any antigenic epitope is important in designing any immunotherapeutic molecule. The understanding of cytokine induction potential also helps us to monitor antigen-specific cellular immune responses and rational vaccine design. The classical immunoinformatics tools served well for prediction of B cell and T cell epitopes. However, in the last decade, the prediction algorithms for T cell epitope inducing specific cytokines have also been developed and appreciated in the scientific community. This review summarizes the current status of such tools, their applications, background algorithms, their use in experimental setup and functionalities available in the tools/web servers.

Quantitative immunopeptidomics reveals a tumor stroma-specific target for T cell therapy

T cell receptor (TCR)-based immunotherapy has emerged as a promising therapeutic approach for the treatment of patients with solid cancers. Identifying peptide-human leukocyte antigen (pHLA) complexes highly presented on tumors and rarely expressed on healthy tissue in combination with high-affinity TCRs that when introduced into T cells can redirect T cells to eliminate tumor but not healthy tissue is a key requirement for safe and efficacious TCR-based therapies. To discover promising shared tumor antigens that could be targeted via TCR-based adoptive T cell therapy, we employed population-scale immunopeptidomics using quantitative mass spectrometry across ~1500 tumor and normal tissue samples. We identified an HLA-A*02:01-restricted pan-cancer epitope within the collagen type VI α-3 (COL6A3) gene that is highly presented on tumor stroma across multiple solid cancers due to a tumor-specific alternative splicing event that rarely occurs outside the tumor microenvironment. T cells expressing natural COL6A3-specific TCRs demonstrated only modest activity against cells presenting high copy numbers of COL6A3 pHLAs. One of these TCRs was affinity-enhanced, enabling transduced T cells to specifically eliminate tumors in vivo that expressed similar copy numbers of pHLAs as primary tumor specimens. The enhanced TCR variants exhibited a favorable safety profile with no detectable off-target reactivity, paving the way to initiate clinical trials using COL6A3-specific TCRs to target an array of solid tumors.

A comprehensive analysis of the IEDB MHC class-I automated benchmark

In 2014, the Immune Epitope Database automated benchmark was created to compare the performance of the MHC class I binding predictors. However, this is not a straightforward process due to the different and non-standardized outputs of the methods. Additionally, some methods are more restrictive regarding the HLA alleles and epitope sizes for which they predict binding affinities, while others are more comprehensive. To address how these problems impacted the ranking of the predictors, we developed an approach to assess the reliability of different metrics. We found that using percentile-ranked results improved the stability of the ranks and allowed the predictors to be reliably ranked despite not being evaluated on the same data. We also found that given the rate new data are incorporated into the benchmark, a new method must wait for at least 4 years to be ranked against the pre-existing methods. The best-performing tools with statistically indistinguishable scores in this benchmark were NetMHCcons, NetMHCpan4.0, ANN3.4, NetMHCpan3.0 and NetMHCpan2.8. The results of this study will be used to improve the evaluation and display of benchmark performance. We highly encourage anyone working on MHC binding predictions to participate in this benchmark to get an unbiased evaluation of their predictors.

MHCRoBERTa: pan-specific peptide-MHC class I binding prediction through transfer learning with label-agnostic protein sequences

Predicting the binding of peptide and major histocompatibility complex (MHC) plays a vital role in immunotherapy for cancer. The success of Alphafold of applying natural language processing (NLP) algorithms in protein secondary struction prediction has inspired us to explore the possibility of NLP methods in predicting peptide-MHC class I binding. Based on the above motivations, we propose the MHCRoBERTa method, RoBERTa pre-training approach, for predicting the binding affinity between type I MHC and peptides. Analysis of the results on benchmark dataset demonstrates that MHCRoBERTa can outperform other state-of-art prediction methods with an increase of the Spearman rank correlation coefficient (SRCC) value. Notably, our model gave a significant improvement on IC50 value. Our method has achieved SRCC value and AUC value as 0.785 and 0.817, respectively. Our SRCC value is 14.3% higher than NetMHCpan3.0 (the second highest SRCC value on pan-specific) and is 3% higher than MHCflurry (the second highest SRCC value on all methods). The AUC value is also better than any other pan-specific methods. Moreover, we visualize the multi-head self-attention for the token representation across the layers and heads by this method. Through the analysis of the representation of each layer and head, we can show whether the model has learned the syntax and semantics necessary to perform the prediction task well. All these results demonstrate that our model can accurately predict the peptide-MHC class I binding affinity and that MHCRoBERTa is a powerful tool for screening potential neoantigens for cancer immunotherapy. MHCRoBERTa is available as an open source software at github (https://github.com/FuxuWang/MHCRoBERTa).

Comparative Immunoreactivity Analyses of Hantaan Virus Glycoprotein-Derived MHC-I Epitopes in Vaccination

MHC-I antigen processes and presentation trigger host-specific anti-viral cellular responses during infection, in which epitope-recognizing cytotoxic T lymphocytes eliminate infected cells and contribute to viral clearance through a cytolytic killing effect. In this study, Hantaan virus (HTNV) GP-derived 9-mer dominant epitopes were obtained with high affinity to major HLA-I and H-2 superfamilies. Further immunogenicity and conservation analyses selected 11 promising candidates, and molecule docking (MD) was then simulated with the corresponding MHC-I alleles. Two-way hierarchical clustering revealed the interactions between GP peptides and MHC-I haplotypes. Briefly, epitope hotspots sharing good affinity to a wide spectrum of MHC-I molecules highlighted the biomedical practice for vaccination, and haplotype clusters represented the similarities among individuals during T-cell response establishment. Cross-validation proved the patterns observed through both MD simulation and public data integration. Lastly, 148 HTNV variants yielded six types of major amino acid residue replacements involving four in nine hotspots, which minimally influenced the general potential of MHC-I superfamily presentation. Altogether, our work comprehensively evaluates the pan-MHC-I immunoreactivity of HTNV GP through a state-of-the-art workflow in light of comparative immunology, acknowledges present discoveries, and offers guidance for ongoing HTNV vaccine pursuit.

Cancer Neoantigens: Challenges and Future Directions for Prediction, Prioritization, and Validation

Prioritization of immunogenic neoantigens is key to enhancing cancer immunotherapy through the development of personalized vaccines, adoptive T cell therapy, and the prediction of response to immune checkpoint inhibition. Neoantigens are tumor-specific proteins that allow the immune system to recognize and destroy a tumor. Cancer immunotherapies, such as personalized cancer vaccines, adoptive T cell therapy, and immune checkpoint inhibition, rely on an understanding of the patient-specific neoantigen profile in order to guide personalized therapeutic strategies. Genomic approaches to predicting and prioritizing immunogenic neoantigens are rapidly expanding, raising new opportunities to advance these tools and enhance their clinical relevance. Predicting neoantigens requires acquisition of high-quality samples and sequencing data, followed by variant calling and variant annotation. Subsequently, prioritizing which of these neoantigens may elicit a tumor-specific immune response requires application and integration of tools to predict the expression, processing, binding, and recognition potentials of the neoantigen. Finally, improvement of the computational tools is held in constant tension with the availability of datasets with validated immunogenic neoantigens. The goal of this review article is to summarize the current knowledge and limitations in neoantigen prediction, prioritization, and validation and propose future directions that will improve personalized cancer treatment.

Identification of alternative splicing-derived cancer neoantigens for mRNA vaccine development

Messenger RNA (mRNA) vaccines have shown great potential for anti-tumor therapy due to the advantages in safety, efficacy and industrial production. However, it remains a challenge to identify suitable cancer neoantigens that can be targeted for mRNA vaccines. Abnormal alternative splicing occurs in a variety of tumors, which may result in the translation of abnormal transcripts into tumor-specific proteins. High-throughput technologies make it possible for systematic characterization of alternative splicing as a source of suitable target neoantigens for mRNA vaccine development. Here, we summarized difficulties and challenges for identifying alternative splicing-derived cancer neoantigens from RNA-seq data and proposed a conceptual framework for designing personalized mRNA vaccines based on alternative splicing-derived cancer neoantigens. In addition, several points were presented to spark further discussion toward improving the identification of alternative splicing-derived cancer neoantigens.

GroEL/S overexpression helps to purge deleterious mutations and reduce genetic diversity during adaptive protein evolution

Chaperones are proteins that help other proteins fold. They also affect the adaptive evolution of their client proteins by buffering the effect of deleterious mutations and increasing the genetic diversity of evolving proteins. We study how the bacterial chaperone GroE (GroEL + GroES) affects the evolution of green fluorescent protein (GFP). To this end we subjected GFP to multiple rounds of mutation and selection for its color phenotype in four replicate E. coli populations, and studied its evolutionary dynamics through high-throughput sequencing and mutant engineering. We evolved GFP both under stabilizing selection for its ancestral (green) phenotype, and to directional selection for a new (cyan) phenotype. We did so both under low and high expression of the chaperone GroE. In contrast to previous work, we observe that GroE does not just buffer but also helps purge deleterious (fluorescence reducing) mutations from evolving populations. In doing so, GroE helps reduce the genetic diversity of evolving populations. In addition, it causes phenotypic heterogeneity in mutants with the same genotype, helping to enhance their fluorescence in some cells, and reducing it in others. Our observations show that chaperones can affect adaptive evolution in more than one way.

CRESSP: a comprehensive pipeline for prediction of immunopathogenic SARS-CoV-2 epitopes using structural properties of proteins

The development of autoimmune diseases following SARS-CoV-2 infection, including multisystem inflammatory syndrome, has been reported, and several mechanisms have been suggested, including molecular mimicry. We developed a scalable, comparative immunoinformatics pipeline called cross-reactive-epitope-search-using-structural-properties-of-proteins (CRESSP) to identify cross-reactive epitopes between a collection of SARS-CoV-2 proteomes and the human proteome using the structural properties of the proteins. Overall, by searching 4 911 245 proteins from 196 352 SARS-CoV-2 genomes, we identified 133 and 648 human proteins harboring potential cross-reactive B-cell and CD8+ T-cell epitopes, respectively. To demonstrate the robustness of our pipeline, we predicted the cross-reactive epitopes of coronavirus spike proteins, which were recognized by known cross-neutralizing antibodies. Using single-cell expression data, we identified PARP14 as a potential target of intermolecular epitope spreading between the virus and human proteins. Finally, we developed a web application (https://ahs2202.github.io/3M/) to interactively visualize our results. We also made our pipeline available as an open-source CRESSP package (https://pypi.org/project/cressp/), which can analyze any two proteomes of interest to identify potentially cross-reactive epitopes between the proteomes. Overall, our immunoinformatic resources provide a foundation for the investigation of molecular mimicry in the pathogenesis of autoimmune and chronic inflammatory diseases following COVID-19.

Identification of Neoantigens in Cancer Cells as Targets for Immunotherapy

The clinical benefits of immune checkpoint blockage (ICB) therapy have been widely reported. In patients with cancer, researchers have demonstrated the clinical potential of antitumor cytotoxic T cells that can be reinvigorated or enhanced by ICB. Compared to self-antigens, neoantigens derived from tumor somatic mutations are believed to be ideal immune targets in tumors. Candidate tumor neoantigens can be identified through immunogenomic or immunopeptidomic approaches. Identification of neoantigens has revealed several points of the clinical relevance. For instance, tumor mutation burden (TMB) may be an indicator of immunotherapy. In various cancers, mutation rates accompanying neoantigen loads may be indicative of immunotherapy. Furthermore, mismatch repair-deficient tumors can be eradicated by T cells in ICB treatment. Hence, immunotherapies using vaccines or adoptive T-cell transfer targeting neoantigens are potential innovative strategies. However, significant efforts are required to identify the optimal epitopes. In this review, we summarize the recent progress in the identification of neoantigens and discussed preclinical and clinical studies based on neoantigens. We also discuss the issues remaining to be addressed before clinical applications of these new therapeutic strategies can be materialized.

In silico analyses and experimental validation of the MHC class-I restricted epitopes of Ebolavirus GP

Ebolavirus (EBOV) causes an extremely high mortality and prevalence disease called Ebola virus disease (EVD). There is only one glycoprotein (GP) on the virus particle surface, which mediates entry into the host cell. MHC class-I restricted CD8 + T cell responses are important antiviral immune responses. Therefore, it is of great importance to understand EBOV GP-specific MHC class-I restricted epitopes within immunogenicity. In this study, computational approaches were employed to predict the dominant MHC class-I molecule epitopes of EBOV GP for mouse H2 and major alleles of HLA class-I supertypes. Our results yielded 42 dominant epitopes in H2 haplotypes and 301 dominant epitopes in HLA class-I haplotypes. After validation by ELISpot assay, in-depth analyses to ascertain their nature of conservation, immunogenicity, and docking with the corresponding MHC class-I molecules were undertaken. Our study predicted MHC class-I restricted epitopes that may aid the advancement of anti-EBOV immune responses. And the integrated strategy of epitope prediction, validation, and comparative analyses were postulated, promising for epitope-based immunotherapy development and application to viral epidemics.

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.

In silico tools and databases for designing cancer immunotherapy

Immunotherapy is a rapidly growing therapy for cancer which have numerous benefits over conventional treatments like surgery, chemotherapy, and radiation. Overall survival of cancer patients has improved significantly due to the use of immunotherapy. It acts as a novel pillar for treating different malignancies from their primary to the metastatic stage. Recent preferments in high-throughput sequencing and computational immunology leads to the development of targeted immunotherapy for precision oncology. In the last few decades, several computational methods and resources have been developed for designing immunotherapy against cancer. In this review, we have summarized cancer-associated genomic, transcriptomic, and mutation profile repositories. We have also enlisted in silico methods for the prediction of vaccine candidates, HLA binders, cytokines inducing peptides, and potential neoepitopes. Of note, we have incorporated the most important bioinformatics pipelines and resources for the designing of cancer immunotherapy. Moreover, to facilitate the scientific community, we have developed a web portal entitled ImmCancer (https://webs.iiitd.edu.in/raghava/immcancer/), comprises cancer immunotherapy tools and repositories.

Low neoantigen expression and poor T-cell priming underlie early immune escape in colorectal cancer

Immune evasion is a hallmark of cancer, and therapies that restore immune surveillance have proven highly effective in cancers with high tumor mutation burden (TMB) (e.g., those with microsatellite instability (MSI)). Whether low TMB cancers, which are largely refractory to immunotherapy, harbor potentially immunogenic neoantigens remains unclear. Here, we show that tumors from all patients with microsatellite stable (MSS) colorectal cancer (CRC) express clonal predicted neoantigens despite low TMB. Unexpectedly, these neoantigens are broadly expressed at lower levels compared to those in MSI CRC. Using a versatile platform for modulating neoantigen expression in CRC organoids and transplantation into the distal colon of mice, we show that low expression precludes productive cross priming and drives immediate T cell dysfunction. Strikingly, experimental or therapeutic rescue of priming rendered T cells capable of controlling tumors with low neoantigen expression. These findings underscore a critical role of neoantigen expression level in immune evasion and therapy response.

Infectious disease mRNA vaccines and a review on epitope prediction for vaccine design

Messenger RNA (mRNA) vaccines have recently emerged as a new type of vaccine technology, showing strong potential to combat the COVID-19 pandemic. In addition to SARS-CoV-2 which caused the pandemic, mRNA vaccines have been developed and tested to prevent infectious diseases caused by other viruses such as Zika virus, the dengue virus, the respiratory syncytial virus, influenza H7N9 and Flavivirus. Interestingly, mRNA vaccines may also be useful for preventing non-infectious diseases such as diabetes and cancer. This review summarises the current progresses of mRNA vaccines designed for a range of diseases including COVID-19. As epitope study is a primary component in the in silico design of mRNA vaccines, we also survey on advanced bioinformatics and machine learning algorithms which have been used for epitope prediction, and review on user-friendly software tools available for this purpose. Finally, we discuss some of the unanswered concerns about mRNA vaccines, such as unknown long-term side effects, and present with our perspectives on future developments in this exciting area.

In-silico design of envelope based multi-epitope vaccine candidate against Kyasanur forest disease virus

Kyasanur forest disease virus (KFDV) causing tick-borne hemorrhagic fever which was earlier endemic to western Ghats, southern India, it is now encroaching into new geographic regions, but there is no approved medicine or effective vaccine against this deadly disease. In this study, we did in-silico design of multi-epitope subunit vaccine for KFDV. B-cell and T-cell epitopes were predicted from conserved regions of KFDV envelope protein and two vaccine candidates (VC1 and VC2) were constructed, those were found to be non-allergic and possess good antigenic properties, also gives cross-protection against Alkhurma hemorrhagic fever virus. The 3D structures of vaccine candidates were built and validated. Docking analysis of vaccine candidates with toll-like receptor-2 (TLR-2) by Cluspro and PatchDock revealed strong affinity between VC1 and TLR2. Ligplot tool was identified the intermolecular hydrogen bonds between vaccine candidates and TLR-2, iMOD server confirmed the stability of the docking complexes. JCAT sever ensured cloning efficiency of both vaccine constructs and in-silico cloning into pET30a (+) vector by SnapGene showed successful translation of epitope region. IMMSIM server was identified increased immunological responses. Finally, multi-epitope vaccine candidates were designed and validated their efficiency, it may pave the way for up-coming vaccine and diagnostic kit development.

The Cancer Epitope Database and Analysis Resource: A Blueprint for the Establishment of a New Bioinformatics Resource for Use by the Cancer Immunology Community

Recent years have witnessed a dramatic rise in interest towards cancer epitopes in general and particularly neoepitopes, antigens that are encoded by somatic mutations that arise as a consequence of tumorigenesis. There is also an interest in the specific T cell and B cell receptors recognizing these epitopes, as they have therapeutic applications. They can also aid in basic studies to infer the specificity of T cells or B cells characterized in bulk and single-cell sequencing data. The resurgence of interest in T cell and B cell epitopes emphasizes the need to catalog all cancer epitope-related data linked to the biological, immunological, and clinical contexts, and most importantly, making this information freely available to the scientific community in a user-friendly format. In parallel, there is also a need to develop resources for epitope prediction and analysis tools that provide researchers access to predictive strategies and provide objective evaluations of their performance. For example, such tools should enable researchers to identify epitopes that can be effectively used for immunotherapy or in defining biomarkers to predict the outcome of checkpoint blockade therapies. We present here a detailed vision, blueprint, and work plan for the development of a new resource, the Cancer Epitope Database and Analysis Resource (CEDAR). CEDAR will provide a freely accessible, comprehensive collection of cancer epitope and receptor data curated from the literature and provide easily accessible epitope and T cell/B cell target prediction and analysis tools. The curated cancer epitope data will provide a transparent benchmark dataset that can be used to assess how well prediction tools perform and to develop new prediction tools relevant to the cancer research community.

Personalized Neoantigen-Pulsed DC Vaccines: Advances in Clinical Applications

In the past few decades, great progress has been made in the clinical application of dendritic cell (DC) vaccines loaded with personalized neoantigens. Personalized neoantigens are antigens arising from somatic mutations in cancers, with specificity to each patient. DC vaccines work based on the fundamental characteristics of DCs, which are professional antigen-presenting cells (APCs), responsible for the uptake, processing, and presentation of antigens to T cells to activate immune responses. Neoantigens can exert their antitumor effects only after they are taken up by APCs and presented to T cells. In recent years, neoantigen-based personalized tumor therapeutic vaccines have proven to be safe, immunogenic and feasible treatment strategies in patients with melanoma and glioblastoma that provide new hope in the treatment of cancer patients and a new approach to cure cancer. In addition, according to ClinicalTrials.gov, hundreds of registered DC vaccine trials are either completed or ongoing worldwide, of which 9 are in early phase I, 191 in phase I, 166 in phase II and 8 in phase III. Hundreds of clinical studies on therapeutic tumor vaccines globally have proven that DC vaccines are stable, reliable and very safe. However, in this process, many other factors still limit the effectiveness of the vaccine. This review will focus on the current research progress on personalized neoantigen-pulsed DC vaccines, their limitations and future research directions of DC vaccines loaded with neoantigens. This review aims to provide a better understanding of DCs biology and manipulation of activated DCs for DCs researchers to produce the next generation of highly efficient cancer vaccines for patients.

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.