NetMHCcons: a consensus method for the major histocompatibility complex class I predictions

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.

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.

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.

Artificial intelligence and neoantigens: paving the path for precision cancer immunotherapy

Cancer immunotherapy has witnessed rapid advancement in recent years, with a particular focus on neoantigens as promising targets for personalized treatments. The convergence of immunogenomics, bioinformatics, and artificial intelligence (AI) has propelled the development of innovative neoantigen discovery tools and pipelines. These tools have revolutionized our ability to identify tumor-specific antigens, providing the foundation for precision cancer immunotherapy. AI-driven algorithms can process extensive amounts of data, identify patterns, and make predictions that were once challenging to achieve. However, the integration of AI comes with its own set of challenges, leaving space for further research. With particular focus on the computational approaches, in this article we have explored the current landscape of neoantigen prediction, the fundamental concepts behind, the challenges and their potential solutions providing a comprehensive overview of this rapidly evolving field.

Systems approach to design multi-epitopic peptide vaccine candidate against fowl adenovirus structural proteins for Gallus gallus domesticus

Introduction

Fowl adenovirus (FAdV) is a significant pathogen in poultry, causing various diseases such as hepatitis-hydropericardium, inclusion body hepatitis, and gizzard erosion. Different serotypes of FAdV are associated with specific conditions, highlighting the need for targeted prevention strategies. Given the rising prevalence of FAdV-related diseases globally, effective vaccination and biosecurity measures are crucial. In this study, we explore the potential of structural proteins to design a multi-epitope vaccine targeting FAdV.

Methods

We employed an in silico approach to design the multi-epitope vaccine. Essential viral structural proteins, including hexon, penton, and fiber protein, were selected as vaccine targets. T-cell and B-cell epitopes binding to MHC-I and MHC-II molecules were predicted using computational methods. Molecular docking studies were conducted to validate the interaction of the multi-epitope vaccine candidate with chicken Toll-like receptors 2 and 5.

Results

Our in silico methodology successfully identified potential T-cell and B-cell epitopes within the selected viral structural proteins. Molecular docking studies revealed strong interactions between the multi-epitope vaccine candidate and chicken Toll-like receptors 2 and 5, indicating the structural integrity and immunogenic potential of the designed vaccine.

Discussion

The designed multi-epitope vaccine presents a promising approach for combating FAdV infections in chickens. By targeting essential viral structural proteins, the vaccine is expected to induce a robust immunological response. The in silico methodology utilized in this study provides a rapid and cost-effective means of vaccine design, offering insights into potential vaccine candidates before experimental validation. Future studies should focus on in vitro and in vivo evaluations to further assess the efficacy and safety of the proposed vaccine.

NeoAgDT: optimization of personal neoantigen vaccine composition by digital twin simulation of a cancer cell population

MOTIVATION

Neoantigen vaccines make use of tumor-specific mutations to enable the patient's immune system to recognize and eliminate cancer. Selecting vaccine elements, however, is a complex task which needs to take into account not only the underlying antigen presentation pathway but also tumor heterogeneity.

RESULTS

Here, we present NeoAgDT, a two-step approach consisting of: (i) simulating individual cancer cells to create a digital twin of the patient's tumor cell population and (ii) optimizing the vaccine composition by integer linear programming based on this digital twin. NeoAgDT shows improved selection of experimentally validated neoantigens over ranking-based approaches in a study of seven patients.

AVAILABILITY AND IMPLEMENTATION

The NeoAgDT code is published on Github: https://github.com/nec-research/neoagdt.

Prognostic impact of HLA supertype mismatch in single-unit cord blood transplantation

The "human leukocyte antigen (HLA) supertype" is a functional classification of HLA alleles, which was defined by structural features and peptide specificities, and has been reportedly associated with the clinical outcomes of viral infections and autoimmune diseases. Although the disparity in each HLA locus was reported to have no clinical significance in single-unit cord blood transplantation (sCBT), the clinical significance of the HLA supertype in sCBT remains unknown. Therefore, we retrospectively analyzed clinical data of 1603 patients who received sCBT in eight institutes in Japan between 2000 and 2017. Each HLA allele was categorized into 19 supertypes, and the prognostic effect of disparities was then assessed. An HLA-B supertype mismatch was identified as a poor prognostic factor (PFS: hazard ratio [HR] = 1.23, p = 0.00044) and was associated with a higher cumulative incidence (CI) of relapse (HR = 1.24, p = 0.013). However, an HLA-B supertype mismatch was not associated with the CI of acute and chronic graft-versus-host-disease. The multivariate analysis for relapse and PFS showed the significance of an HLA-B supertype mismatch independent of allelic mismatches, and other previously reported prognostic factors. HLA-B supertype-matched grafts should be selected in sCBT.

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.

DNAJB1-PRKACA fusion neoantigens elicit rare endogenous T cell responses that potentiate cell therapy for fibrolamellar carcinoma

Fibrolamellar carcinoma (FLC) is a liver tumor with a high mortality burden and few treatment options. A promising therapeutic vulnerability in FLC is its driver mutation, a conserved DNAJB1-PRKACA gene fusion that could be an ideal target neoantigen for immunotherapy. In this study, we aim to define endogenous CD8 T cell responses to this fusion in FLC patients and evaluate fusion-specific T cell receptors (TCRs) for use in cellular immunotherapies. We observe that fusion-specific CD8 T cells are rare and that FLC patient TCR repertoires lack large clusters of related TCR sequences characteristic of potent antigen-specific responses, potentially explaining why endogenous immune responses are insufficient to clear FLC tumors. Nevertheless, we define two functional fusion-specific TCRs, one of which has strong anti-tumor activity in vivo. Together, our results provide insights into the fragmented nature of neoantigen-specific repertoires in humans and indicate routes for clinical development of successful immunotherapies for FLC.

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.

Leveraging the immunoinformatics approach for designing the SARS-CoV-2 omicron-specific antigenic cassette of mRNA vaccine

Emergence of SARS-CoV-2 Omicron variant has presented a significant challenge to global health, demanding rapid development of mRNA-based vaccines. The mRNA-guided vaccine platforms offer various advantages over traditional vaccine platforms. The mRNA by nature is a short-lived molecule that guides the cells to manufacture antigenic proteins. In the present work, we have created an omicron spike antigenic protein sequence characterized by base composition analysis, modeling, and docking with the ACE-2 receptor. Further, we predicted the B-cell and T-cell epitopes followed by antigenicity, toxicity, and allergenicity. Finally, the protein was reverse translated, codon-optimized, and encoding mRNA sequence was checked for its stability by predicting the secondary structures. A comprehensive examination of in-silico data revealed 628.2 as a potent antigenic candidate that was finally used in Gemcovac®-OM, a heterologous booster mRNA vaccine for COVID-19.

Transfer learning improves pMHC kinetic stability and immunogenicity predictions

The cellular immune response comprises several processes, with the most notable ones being the binding of the peptide to the Major Histocompability Complex (MHC), the peptide-MHC (pMHC) presentation to the surface of the cell, and the recognition of the pMHC by the T-Cell Receptor. Identifying the most potent peptide targets for MHC binding, presentation and T-cell recognition is vital for developing peptide-based vaccines and T-cell-based immunotherapies. Data-driven tools that predict each of these steps have been developed, and the availability of mass spectrometry (MS) datasets has facilitated the development of accurate Machine Learning (ML) methods for class-I pMHC binding prediction. However, the accuracy of ML-based tools for pMHC kinetic stability prediction and peptide immunogenicity prediction is uncertain, as stability and immunogenicity datasets are not abundant. Here, we use transfer learning techniques to improve stability and immunogenicity predictions, by taking advantage of a large number of binding affinity and MS datasets. The resulting models, TLStab and TLImm, exhibit comparable or better performance than state-of-the-art approaches on different stability and immunogenicity test sets respectively. Our approach demonstrates the promise of learning from the task of peptide binding to improve predictions on downstream tasks. The source code of TLStab and TLImm is publicly available at https://github.com/KavrakiLab/TL-MHC.

Evaluating the genetic diversity of the Plasmodium vivax siap2 locus: A promising candidate for an effective malaria vaccine?

Malaria is the deadliest parasitic disease in the world. Traditional control measures have become less effective; hence, there is a need to explore alternative strategies, such as antimalarial vaccines. However, designing an anti-Plasmodium vivax vaccine is considered a challenge due to the complex parasite biology and the antigens' high genetic diversity. Recently, the sporozoite invasion-associated protein 2 (SIAP2) has been suggested as a potential antigen to be considered in vaccine design due to its significance during hepatocyte invasion. However, its use may be limited by the incomplete understanding of gene/protein diversity. Here, the genetic diversity of pvsiap2 using P. vivax DNA samples from Colombia was assessed. Through PCR amplification and sequencing, we compared the Colombian sequences with available worldwide sequences, revealing that pvsiap2 displays low genetic diversity. Molecular evolutionary analyses showed that pvsiap2 appears to be influenced by directional selection. Moreover, the haplotypes found differ by a few mutational steps and several of them were shared between different geographical areas. On the other hand, several conserved regions within PvSIAP2 were predicted as potential B-cell or T-cell epitopes. Considering these characteristics and its role in hepatocyte invasion, the PvSIAP2 protein emerges as a promising antigen to be considered in a multi-antigen-multi-stage (multivalent) fully effective vaccine against P. vivax malaria.

An Efficient Approach to the Accurate Prediction of Mutational Effects in Antigen Binding to the MHC1

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

HLAEquity: Examining biases in pan-allele peptide-HLA binding predictors

Peptide-HLA (pHLA) binding prediction is essential in screening peptide candidates for personalized peptide vaccines. Machine learning (ML) pHLA binding prediction tools are trained on vast amounts of data and are effective in screening peptide candidates. Most ML models report the ability to generalize to HLA alleles unseen during training ("pan-allele" models). However, the use of datasets with imbalanced allele content raises concerns about biased model performance. First, we examine the data bias of two ML-based pan-allele pHLA binding predictors. We find that the pHLA datasets overrepresent alleles from geographic populations of high-income countries. Second, we show that the identified data bias is perpetuated within ML models, leading to algorithmic bias and subpar performance for alleles expressed in low-income geographic populations. We draw attention to the potential therapeutic consequences of this bias, and we challenge the use of the term "pan-allele" to describe models trained with currently available public datasets.

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.

The SysteMHC Atlas v2.0, an updated resource for mass spectrometry-based immunopeptidomics

The SysteMHC Atlas v1.0 was the first public repository dedicated to mass spectrometry-based immunopeptidomics. Here we introduce a newly released version of the SysteMHC Atlas v2.0 (https://systemhc.sjtu.edu.cn), a comprehensive collection of 7190 MS files from 303 allotypes. We extended and optimized a computational pipeline that allows the identification of MHC-bound peptides carrying on unexpected post-translational modifications (PTMs), thereby resulting in 471K modified peptides identified over 60 distinct PTM types. In total, we identified approximately 1.0 million and 1.1 million unique peptides for MHC class I and class II immunopeptidomes, respectively, indicating a 6.8-fold increase and a 28-fold increase to those in v1.0. The SysteMHC Atlas v2.0 introduces several new features, including the inclusion of non-UniProt peptides, and the incorporation of several novel computational tools for FDR estimation, binding affinity prediction and motif deconvolution. Additionally, we enhanced the user interface, upgraded website framework, and provided external links to other resources related. Finally, we built and provided various spectral libraries as community resources for data mining and future immunopeptidomic and proteomic analysis. We believe that the SysteMHC Atlas v2.0 is a unique resource to provide key insights to the immunology and proteomics community and will accelerate the development of vaccines and immunotherapies.

DeepHLApan: A Deep Learning Approach for the Prediction of Peptide-HLA Binding and Immunogenicity

Neoantigens are crucial in distinguishing cancer cells from normal ones and play a significant role in cancer immunotherapy. The field of bioinformatics prediction for tumor neoantigens has rapidly developed, focusing on the prediction of peptide-HLA binding affinity. In this chapter, we introduce a user-friendly tool named DeepHLApan, which utilizes deep learning techniques to predict neoantigens by considering both peptide-HLA binding affinity and immunogenicity. We provide the application of DeepHLApan, along with the source code, docker version, and web-server. These resources are freely available at https://github.com/zjupgx/deephlapan and http://pgx.zju.edu.cn/deephlapan/ .

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.

Development and evaluation of a multi-epitope subunit vaccine against Mycoplasma synoviae infection

Mycoplasma synoviae is an extremely significant avian pathogen, causing substantial financial harm to poultry farmers worldwide, and impacting both chicken and turkey production. Multi-epitope vaccines offer higher immunity and lower allergenicity compared to conventional vaccines. In this study, our objective is to develop a multi-epitope vaccine for M. synoviae (MSMV) and to evaluate the immune responses and protective efficacy of MSMV in chickens. We successfully identified a total of 14 B-cell, 5 MHC-I, and 16 MHC-II binding epitopes from the immunodominant proteins RS01790, BMP, GrpE, RS00900, and RS00275. Subsequently, we synthesized the multi-epitope vaccine by connecting all conserved epitopes using appropriate linkers. The resulting MSMV demonstrated notable antigenicity, non-allergenic properties, and stability. Notably, the MSMV effectively stimulated high levels of antibody production in chickens. Furthermore, MSMV the vaccine elicited a robust cellular immune response in chickens, characterized by a well-balanced Th1/Th2-type cytokine profile and enhanced lymphocyte proliferation. In immune protection experiments, the vaccinated chickens exhibited reduced air sac lesion scores and tracheal mucosal thickness compared to their non-vaccinated chickens. Additionally, vaccinated chickens displayed lower M. synoviae loads in throat swabs. These findings collectively suggested that the MSMV holds significant potential as a promising vaccine candidate for managing M. synoviae infections.

Prioritization of infectious epitopes for translational investigation in type 1 diabetes etiology

Molecular mimicry is one mechanism by which infectious agents are thought to trigger islet autoimmunity in type 1 diabetes. With a growing number of reported infectious agents and islet antigens, strategies to prioritize the study of infectious agents are critically needed to expedite translational research into the etiology of type 1 diabetes. In this work, we developed an in-silico pipeline for assessing molecular mimicry in type 1 diabetes etiology based on sequence homology, empirical binding affinity to specific MHC molecules, and empirical potential for T-cell immunogenicity. We then assess whether potential molecular mimics were conserved across other pathogens known to infect humans. Overall, we identified 61 potentially high-impact molecular mimics showing sequence homology, strong empirical binding affinity, and empirical immunogenicity linked with specific MHC molecules. We further found that peptide sequences from 32 of these potential molecular mimics were conserved across several human pathogens. These findings facilitate translational evaluation of molecular mimicry in type 1 diabetes etiology by providing a curated and prioritized list of peptides from infectious agents for etiopathologic investigation. These results may also provide evidence for generation of infectious and HLA-specific preclinical models and inform future screening and preventative efforts in genetically susceptible populations.

Driver mutations in GNAQ and GNA11 genes as potential targets for precision immunotherapy in uveal melanoma patients

Uveal melanoma (UM) is the most common ocular malignancy in adults. Nearly 95% of UM patients carry the mutually exclusive mutations in the homologous genes GNAQ (amino acid change Q209L/Q209P) and GNA11 (aminoacid change Q209L). UM is located in an immunosuppressed organ and does not suffer immunoediting. Therefore, we hypothesize that driver mutations in GNAQ/11 genes could be recognized by the immune system. Genomic and transcriptomic data from primary uveal tumors were collected from the TCGA-UM dataset (n = 80) and used to assess the immunogenic potential for GNAQ/GNA11 Q209L/Q209P mutations using a variety of tools and HLA type information. All prediction tools showed stronger GNAQ/11 Q209L binding to HLA than GNAQ/11 Q209P. The immunogenicity analysis revealed that Q209L is likely to be presented by more than 73% of individuals in 1000 G databases whereas Q209P is only predicted to be presented in 24% of individuals. GNAQ/11 Q209L showed a higher likelihood to be presented by HLA-I molecules than almost all driver mutations analyzed. Finally, samples carrying Q209L had a higher immune-reactive phenotype. Regarding cancer risk, seven HLA genotypes with low Q209L affinity show higher frequency in uveal melanoma patients than in the general population. However, no clear association was found between any HLA genotype and survival. Results suggest a high potential immunogenicity of the GNAQ/11 Q209L variant that could allow the generation of novel therapeutic tools to treat UM like neoantigen vaccinations.

Bioinformatics analysis of Muscovy duck parvovirus REP and VP1 proteins

This article was aimed at analyzing the sequence, structure, and function of the two Muscovy duck parvovirus proteins, including REP and VP1. The antigenicity, physical and chemical properties, transmembrane regions, phosphorylation sites, glycosylation sites, three-dimensional structure, and linear epitope of VP1 and REP were predicted and analyzed through bioinformatics methods. A multi-epitope vaccine was also constructed based on the screened epitopes, and the vaccine was characterized, modeled, molecularly docked and molecularly cloned. The epitopes were screened according to the criteria of antigenicity, non-allergenicity and non-toxicity, and 12 epitope fragments were obtained. The B cell epitopes were analyzed according to four scales: β-turn, hydrophilicity, surface accessibility and antigenicity. Combined with the epitope prediction results based on structure, the final epitope prediction results were obtained. The multi-epitope vaccine used an EAAAK-linked adjuvant, a GPGPG-linked T-cell epitope, and a KK-linked B-cell epitope. The analysis showed that the vaccine was stable hydrophilic, antigenic, conserved and non-allergenic. Based on molecular docking it was shown that good interactions between the vaccine and the immune receptor were generated and were essential to generate an immune response. The final vaccine was reverse translated into cDNA and the DNA vaccine was designed by codon optimization and molecular cloning. Further trials are still needed to demonstrate the immunogenicity and other aspects of vaccine efficacy.Communicated by Ramaswamy H. Sarma.

An association study of HLA with levofloxacin-induced severe cutaneous adverse drug reactions in Han Chinese

Levofloxacin-induced severe cutaneous adverse drug reactions (LEV-SCARs) remain unexplored. An association study of human leukocyte antigen (HLA) alleles with LEV-SCARs among 12 patients, 806 healthy subjects, and 100 levofloxacin-tolerant individuals was performed. The carrier frequencies of HLA-B∗13:01 (odds ratio [OR]: 4.50; 95% confidence interval [CI]: 1.15-17.65; p = 0.043), HLA-B∗13:02 (OR: 6.14; 95% CI: 1.73-21.76; p = 0.0072), and serotype B13 (OR: 17.73; 95% CI: 3.61-86.95; p = 4.85 × 10-5) in patients with LEV-SCARs were significantly higher than those of levofloxacin-tolerant individuals. Molecular docking analysis suggested that levofloxacin formed more stable binding models with HLA-B∗13:01 and HLA-B∗13:02 than with non-risk HLA-B∗46:01. Mass spectrometry revealed that nonapeptides bound to HLA-B∗13:02 shifted at several positions after exposure to levofloxacin. Prospective screening for serotype B13 (sensitivity: 83%, specificity: 78%) and alternative drug treatment for carriers may significantly decrease the incidence of LEV-SCARs.

Identification of vaccine candidates against rhodococcus equi by combining pangenome analysis with a reverse vaccinology approach

Rhodococcus equi (R. equi) is a zoonotic opportunistic pathogen that can cause life-threatening infections. The rapid evolution of multidrug-resistant R. equi and the fact that there is no currently licensed effective vaccine against R. equi warrant the need for vaccine development. Reverse vaccinology (RV), which involves screening a pathogen's entire genome and proteome using various web-based prediction tools, is considered one of the most effective approaches for identifying vaccine candidates. Here, we performed a pangenome analysis to determine the core proteins of R. equi. We then used the RV approach to examine the subcellular localization, host and gut flora homology, antigenicity, transmembrane helices, physicochemical properties, and immunogenicity of the core proteins to select potential vaccine candidates. The vaccine candidates were then subjected to epitope mapping to predict the exposed antigenic epitopes that possess the ability to bind with major histocompatibility complex I/II (MHC I/II) molecules. These vaccine candidates and epitopes will form a library of elements for the development of a polyvalent or universal vaccine against R. equi. Sixteen R. equi complete proteomes were found to contain 6,238 protein families, and the core proteins consisted of 3,969 protein families (∼63.63% of the pangenome), reflecting a low degree of intraspecies genomic variability. From the pool of core proteins, 483 nonhost homologous membrane and extracellular proteins were screened, and 12 vaccine candidates were finally identified according to their antigenicity, physicochemical properties and other factors. These included four cell wall/membrane/envelope biogenesis proteins; four amino acid transport and metabolism proteins; one cell cycle control, cell division and chromosome partitioning protein; one carbohydrate transport and metabolism protein; one secondary metabolite biosynthesis, transport and catabolism protein; and one defense mechanism protein. All 12 vaccine candidates have an experimentally validated 3D structure available in the protein data bank (PDB). Epitope mapping of the candidates showed that 16 MHC I epitopes and 13 MHC II epitopes with the strongest immunogenicity were exposed on the protein surface, indicating that they could be used to develop a polypeptide vaccine. Thus, we utilized an analytical strategy that combines pangenome analysis and RV to generate a peptide antigen library that simplifies the development of multivalent or universal vaccines against R. equi and can be applied to the development of other vaccines.

Next-Generation TB Vaccines: Progress, Challenges, and Prospects

Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), is a prevalent global infectious disease and a leading cause of mortality worldwide. Currently, the only available vaccine for TB prevention is Bacillus Calmette-Guérin (BCG). However, BCG demonstrates limited efficacy, particularly in adults. Efforts to develop effective TB vaccines have been ongoing for nearly a century. In this review, we have examined the current obstacles in TB vaccine research and emphasized the significance of understanding the interaction mechanism between MTB and hosts in order to provide new avenues for research and establish a solid foundation for the development of novel vaccines. We have also assessed various TB vaccine candidates, including inactivated vaccines, attenuated live vaccines, subunit vaccines, viral vector vaccines, DNA vaccines, and the emerging mRNA vaccines as well as virus-like particle (VLP)-based vaccines, which are currently in preclinical stages or clinical trials. Furthermore, we have discussed the challenges and opportunities associated with developing different types of TB vaccines and outlined future directions for TB vaccine research, aiming to expedite the development of effective vaccines. This comprehensive review offers a summary of the progress made in the field of novel TB vaccines.

CapsNet-MHC predicts peptide-MHC class I binding based on capsule neural networks

The Major Histocompatibility Complex (MHC) binds to the derived peptides from pathogens to present them to killer T cells on the cell surface. Developing computational methods for accurate, fast, and explainable peptide-MHC binding prediction can facilitate immunotherapies and vaccine development. Various deep learning-based methods rely on separate feature extraction from the peptide and MHC sequences and ignore their pairwise binding information. This paper develops a capsule neural network-based method to efficiently capture the peptide-MHC complex features to predict the peptide-MHC class I binding. Various evaluations confirmed our method outperformance over the alternative methods, while it can provide accurate prediction over less available data. Moreover, for providing precise insights into the results, we explored the essential features that contributed to the prediction. Since the simulation results demonstrated consistency with the experimental studies, we concluded that our method can be utilized for the accurate, rapid, and interpretable peptide-MHC binding prediction to assist biological therapies.

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.

A novel transcriptional signature identifies T-cell infiltration in high-risk paediatric cancer

BACKGROUND

Molecular profiling of the tumour immune microenvironment (TIME) has enabled the rational choice of immunotherapies in some adult cancers. In contrast, the TIME of paediatric cancers is relatively unexplored. We speculated that a more refined appreciation of the TIME in childhood cancers, rather than a reliance on commonly used biomarkers such as tumour mutation burden (TMB), neoantigen load and PD-L1 expression, is an essential prerequisite for improved immunotherapies in childhood solid cancers.

METHODS

We combined immunohistochemistry (IHC) with RNA sequencing and whole-genome sequencing across a diverse spectrum of high-risk paediatric cancers to develop an alternative, expression-based signature associated with CD8⁺ T-cell infiltration of the TIME. Furthermore, we explored transcriptional features of immune archetypes and T-cell receptor sequencing diversity, assessed the relationship between CD8⁺ and CD4⁺ abundance by IHC and deconvolution predictions and assessed the common adult biomarkers such as neoantigen load and TMB.

RESULTS

A novel 15-gene immune signature, Immune Paediatric Signature Score (IPASS), was identified. Using this signature, we estimate up to 31% of high-risk cancers harbour infiltrating T-cells. In addition, we showed that PD-L1 protein expression is poorly correlated with PD-L1 RNA expression and TMB and neoantigen load are not predictive of T-cell infiltration in paediatrics. Furthermore, deconvolution algorithms are only weakly correlated with IHC measurements of T-cells.

CONCLUSIONS

Our data provides new insights into the variable immune-suppressive mechanisms dampening responses in paediatric solid cancers. Effective immune-based interventions in high-risk paediatric cancer will require individualised analysis of the TIME.

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.

Pervasiveness of HLA allele-specific expression loss across tumor types

BACKGROUND

Efficient presentation of mutant peptide fragments by the human leukocyte antigen class I (HLA-I) genes is necessary for immune-mediated killing of cancer cells. According to recent reports, patient HLA-I genotypes can impact the efficacy of cancer immunotherapy, and the somatic loss of HLA-I heterozygosity has been established as a factor in immune evasion. While global deregulated expression of HLA-I has also been reported in different tumor types, the role of HLA-I allele-specific expression loss - that is, the preferential RNA expression loss of specific HLA-I alleles - has not been fully characterized in cancer.

METHODS

Here, we use RNA and whole-exome sequencing data to quantify HLA-I allele-specific expression (ASE) in cancer using our novel method arcasHLA-quant.

RESULTS

We show that HLA-I ASE loss in at least one of the three HLA-I genes is a pervasive phenomenon across TCGA tumor types. In pancreatic adenocarcinoma, tumor-specific HLA-I ASE loss is associated with decreased overall survival specifically in the basal-like subtype, a finding that we validated in an independent cohort through laser-capture microdissection. Additionally, we show that HLA-I ASE loss is associated with poor immunotherapy outcomes in metastatic melanoma through retrospective analyses.

CONCLUSIONS

Together, our results highlight the prevalence of HLA-I ASE loss and provide initial evidence of its clinical significance in cancer prognosis and immunotherapy treatment.

Designing neoantigen cancer vaccines, trials, and outcomes

Neoantigen vaccines are based on epitopes of antigenic parts of mutant proteins expressed in cancer cells. These highly immunogenic antigens may trigger the immune system to combat cancer cells. Improvements in sequencing technology and computational tools have resulted in several clinical trials of neoantigen vaccines on cancer patients. In this review, we have looked into the design of the vaccines which are undergoing several clinical trials. We have discussed the criteria, processes, and challenges associated with the design of neoantigens. We searched different databases to track the ongoing clinical trials and their reported outcomes. We observed, in several trials, the vaccines boost the immune system to combat the cancer cells while maintaining a reasonable margin of safety. Detection of neoantigens has led to the development of several databases. Adjuvants also play a catalytic role in improving the efficacy of the vaccine. Through this review, we can conclude that the efficacy of vaccines can make it a potential treatment across different types of cancers.

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.

Diversity and evolution of computationally predicted T cell epitopes against human respiratory syncytial virus

Human respiratory syncytial virus (RSV) is a major cause of lower respiratory infection. Despite more than 60 years of research, there is no licensed vaccine. While B cell response is a major focus for vaccine design, the T cell epitope profile of RSV is also important for vaccine development. Here, we computationally predicted putative T cell epitopes in the Fusion protein (F) and Glycoprotein (G) of RSV wild circulating strains by predicting Major Histocompatibility Complex (MHC) class I and class II binding affinity. We limited our inferences to conserved epitopes in both F and G proteins that have been experimentally validated. We applied multidimensional scaling (MDS) to construct T cell epitope landscapes to investigate the diversity and evolution of T cell profiles across different RSV strains. We find the RSV strains are clustered into three RSV-A groups and two RSV-B groups on this T epitope landscape. These clusters represent divergent RSV strains with potentially different immunogenic profiles. In addition, our results show a greater proportion of F protein T cell epitope content conservation among recent epidemic strains, whereas the G protein T cell epitope content was decreased. Importantly, our results suggest that RSV-A and RSV-B have different patterns of epitope drift and replacement and that RSV-B vaccines may need more frequent updates. Our study provides a novel framework to study RSV T cell epitope evolution. Understanding the patterns of T cell epitope conservation and change may be valuable for vaccine design and assessment.

Epitope identification for p53R273C mutant

BACKGROUND

With the rise of immunotherapy based on cancer neoantigen, identification of neoepitopes has become an urgent problem to be solved. The TP53 R273C mutation is one of the hotspot mutations of TP53, however, the immunogenicity of this mutation is not yet clear. The aim of this study is to identify potential epitopes for p53R273C mutant.

METHODS

In this study, bioinformatic methods, peptide exchange assay, and peptide-immunized human leukocyte antigen (HLA) transgenic mouse model were used to explore the immunogenicity of this mutation.

RESULTS

Peptides with higher affinity to common HLA-A alleles (A*11:01, A*02:01) were discovered by computational prediction. All the 8-11 mer peptides contain the mutation site were synthesized and soluble peptides were used in the peptide exchange assay. However, the exchange efficiencies of these predicted peptides to HLAs were lower. Fortunately, other peptides with higher exchange efficiency were discovered. Then, the immunogenicity of these peptides was validated with the HLA-A2 transgenic mice model.

CONCLUSION

We identified three potential neoepitopes of p53^R273C for HLA-A*02:01, one potential neoepitope for HLA-A*11:01 and no neoepitope for HLA-A*24:02.

Shared 6mer Peptides of Human and Omicron (21K and 21L) at SARS-CoV-2 Mutation Sites

We investigated the short sequences involving Omicron 21K and Omicron 21L variants to reveal any possible molecular mimicry-associated autoimmunity risks and changes in those. We first identified common 6mers of the viral and human protein sequences present for both the mutant (Omicron) and nonmutant (SARS-CoV-2) versions of the same viral sequence and then predicted the binding affinities of those sequences to the HLA supertype representatives. We evaluated change in the potential autoimmunity risk, through comparative assessment of the nonmutant and mutant viral sequences and their similar human peptides with common 6mers and affinities to the same HLA allele. This change is the lost and the new, or de novo, autoimmunity risk, associated with the mutations in the Omicron 21K and Omicron 21L variants. Accordingly, e.g., the affinity of virus-similar sequences of the Ig heavy chain junction regions shifted from the HLA-B*15:01 to the HLA-A*01:01 allele at the mutant sequences. Additionally, peptides of different human proteins sharing 6mers with SARS-CoV-2 proteins at the mutation sites of interest and with affinities to the HLA-B*07:02 allele, such as the respective SARS-CoV-2 sequences, were lost. Among all, any possible molecular mimicry-associated novel risk appeared to be prominent in HLA-A*24:02 and HLA-B*27:05 serotypes upon infection with Omicron 21L. Associated disease, pathway, and tissue expression data supported possible new risks for the HLA-B*27:05 and HLA-A*01:01 serotypes, while the risks for the HLA-B*07:02 serotypes could have been lost or diminished, and those for the HLA-A*03:01 serotypes could have been retained, for the individuals infected with Omicron variants under study. These are likely to affect the complications related to cross-reactions influencing the relevant HLA serotypes upon infection with Omicron 21K and Omicron 21L.

Seq2Neo: A Comprehensive Pipeline for Cancer Neoantigen Immunogenicity Prediction

Neoantigens derived from somatic DNA alterations are ideal cancer-specific targets. In recent years, the combination therapy of PD-1/PD-L1 blockers and neoantigen vaccines has shown clinical efficacy in original PD-1/PD-L1 blocker non-responders. However, not all somatic DNA mutations result in immunogenicity among cancer cells and efficient tools to predict the immunogenicity of neoepitopes are still urgently needed. Here, we present the Seq2Neo pipeline, which provides a one-stop solution for neoepitope feature prediction using raw sequencing data. Neoantigens derived from different types of genome DNA alterations, including point mutations, insertion deletions and gene fusions, are all supported. Importantly, a convolutional neural network (CNN)-based model was trained to predict the immunogenicity of neoepitopes and this model showed an improved performance compared to the currently available tools in immunogenicity prediction using independent datasets. We anticipate that the Seq2Neo pipeline could become a useful tool in the prediction of neoantigen immunogenicity and cancer immunotherapy. Seq2Neo is open-source software under an academic free license (AFL) v3.0 and is freely available at Github.

Single-cell analysis reveals clonally expanded tumor-associated CD57+ CD8 T cells are enriched in the periphery of patients with metastatic urothelial cancer responding to PD-L1 blockade

BACKGROUND

A growing body of evidence suggests that T-cell responses against neoantigens are critical regulators of response to immune checkpoint blockade. We previously showed that circulating neoantigen-specific CD8 T cells in patients with lung cancer responding to anti-Programmed death-ligand 1 (PD-L1) (atezolizumab) exhibit a unique phenotype with high expression of CD57, CD244, and KLRG1. Here, we extended our analysis on neoantigen-specific CD8 T cells to patients with metastatic urothelial cancer (mUC) and further profiled total CD8 T cells to identify blood-based predictive biomarkers of response to atezolizumab.

METHODS

We identified tumor neoantigens from 20 patients with mUC and profiled their peripheral CD8 T cells using highly multiplexed combinatorial tetramer staining. Another set of patients with mUC treated with atezolizumab (n=30) or chemotherapy (n=40) were selected to profile peripheral CD8 T cells by mass cytometry. Using single-cell transcriptional analysis (single-cell RNA sequencing (scRNA-seq)), together with CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) and paired T-cell receptor (TCR) sequencing, we further characterized peripheral CD8 T cells in a subset of patients (n=16).

RESULTS

High frequency of CD57 was observed in neoantigen-specific CD8 T cells in patients with mUC responding to atezolizumab. Extending these findings to bulk CD8 T cells, we found higher frequency of CD57 expressing CD8 T cells before treatment in patients responding to atezolizumab (n=20, p<0.01) but not to chemotherapy. These findings were corroborated in a validation cohort (n=30, p<0.01) and notably were independent of known biomarkers of response. scRNA-seq analysis identified a clonally expanded cluster enriched within CD57⁺ CD8 T cells in responding patients characterized by higher expression of genes associated with activation, cytotoxicity, and tissue-resident memory markers. Furthermore, compared with CD57^- CD8 T cells, TCRs of CD57⁺ CD8 T cells showed increased overlap with the TCR repertoire of tumor-infiltrating T cells.

CONCLUSIONS

Collectively, we show high frequencies of CD57 among neoantigen-specific and bulk CD8 T cells in patients responding to atezolizumab. The TCR repertoire overlap between peripheral CD57⁺ CD8 T cells and tumor-infiltrating lymphocytes suggest that accumulation of peripheral CD57⁺ CD8 T cells is reflective of an ongoing antitumor T-cell response. Our findings provide evidence and rationale for using circulating CD8 T cells expressing CD57 as a readily accessible blood-based biomarker for selecting patients with mUC for atezolizumab therapy.

In silico SARS-CoV-2 vaccine development for Omicron strain using reverse vaccinology

BACKGROUND

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic began in 2019 but it remains as a serious threat today. To reduce and prevent spread of the virus, multiple vaccines have been developed. Despite the efforts in developing vaccines, Omicron strain of the virus has recently been designated as a variant of concern (VOC) by the World Health Organization (WHO).

OBJECTIVE

To develop a vaccine candidate against Omicron strain (B.1.1.529, BA.1) of the SARS-CoV-19.

METHODS

We applied reverse vaccinology methods for BA.1 and BA.2 as the vaccine target and a control, respectively. First, we predicted MHC I, MHC II and B cell epitopes based on their viral genome sequences. Second, after estimation of antigenicity, allergenicity and toxicity, a vaccine construct was assembled and tested for physicochemical properties and solubility. Third, AlphaFold2, RaptorX and RoseTTAfold servers were used to predict secondary structures and 3D structures of the vaccine construct. Fourth, molecular docking analysis was performed to test binding of our construct with angiotensin converting enzyme 2 (ACE2). Lastly, we compared mutation profiles on the epitopes between BA.1, BA.2, and wild type to estimate the efficacy of the vaccine.

RESULTS

We collected a total of 10 MHC I, 9 MHC II and 5 B cell epitopes for the final vaccine construct for Omicron strain. All epitopes were predicted to be antigenic, non-allergenic and non-toxic. The construct was estimated to have proper stability and solubility. The best modelled tertiary structures were selected for molecular docking analysis with ACE2 receptor.

CONCLUSIONS

These results suggest the potential efficacy of our newly developed vaccine construct as a novel vaccine candidate against Omicron strain of the coronavirus.

Viral immunogenic footprints conferring T cell cross-protection to SARS-CoV-2 and its variants

COVID-19 brought scenes from sci-fi movies into real life. Infected individuals include asymptomatic cases to severe disease leading to death, suggesting the involvement of the genetic constitution of populations and pathogens contributing to differential individuals' outcomes. To investigate shared immunogenic features between SARS-CoV-2 targets and other coronaviruses, we modeled their peptides in 3D structures of HLA-A*02:01 (pMHC), comparing their molecular surfaces These structures were also compared with a panel of epitopes from unrelated viruses, looking for potential triggers conferring cross-protection in uninfected individuals. As expected, SARS-CoV 1 and 2 peptides share molecular and physicochemical features, providing an explanation for the verified experimental immunogenicity among them. Surprisingly, even discordant sequences from human coronaviruses 229E, OC43 and epitopes from unrelated viruses involved in endemic human infections exhibit similar fingerprints of immunogenicity with SARS-CoV-2 peptides. The same approach indicates a conserved CD8+ T cell recognition between Wuhan SARS-CoV-2 sequences and altered peptides from Variants of Concern. Examination of structural data over epitope sequence analysis here could explain how previous infections may produce a heterologous immunity response in a global scale against emergent diseases such as Covid-19, mitigating its full lethal potential, and paves the way for the development of wide spectrum vaccine development.

Neoantigens in precision cancer immunotherapy: from identification to clinical applications

Immunotherapies targeting cancer neoantigens are safe, effective, and precise. Neoantigens can be identified mainly by genomic techniques such as next-generation sequencing and high-throughput single-cell sequencing; proteomic techniques such as mass spectrometry; and bioinformatics tools based on high-throughput sequencing data, mass spectrometry data, and biological databases. Neoantigen-related therapies are widely used in clinical practice and include neoantigen vaccines, neoantigen-specific CD8+ and CD4+ T cells, and neoantigen-pulsed dendritic cells. In addition, neoantigens can be used as biomarkers to assess immunotherapy response, resistance, and prognosis. Therapies based on neoantigens are an important and promising branch of cancer immunotherapy. Unremitting efforts are needed to unravel the comprehensive role of neoantigens in anti-tumor immunity and to extend their clinical application. This review aimed to summarize the progress in neoantigen research and to discuss its opportunities and challenges in precision cancer immunotherapy.

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.

Molecular correlates of clinical response and resistance to atezolizumab in combination with bevacizumab in advanced hepatocellular carcinoma

Atezolizumab (anti-programmed death-ligand 1 (PD-L1)) and bevacizumab (anti-vascular endothelial growth factor (VEGF)) combination therapy has become the new standard of care in patients with unresectable hepatocellular carcinoma. However, potential predictive biomarkers and mechanisms of response and resistance remain less well understood. We report integrated molecular analyses of tumor samples from 358 patients with hepatocellular carcinoma (HCC) enrolled in the GO30140 phase 1b or IMbrave150 phase 3 trial and treated with atezolizumab combined with bevacizumab, atezolizumab alone or sorafenib (multikinase inhibitor). Pre-existing immunity (high expression of CD274, T-effector signature and intratumoral CD8⁺ T cell density) was associated with better clinical outcomes with the combination. Reduced clinical benefit was associated with high regulatory T cell (Treg) to effector T cell (Teff) ratio and expression of oncofetal genes (GPC3, AFP). Improved outcomes from the combination versus atezolizumab alone were associated with high expression of VEGF Receptor 2 (KDR), Tregs and myeloid inflammation signatures. These findings were further validated by analyses of paired pre- and post-treatment biopsies, in situ analyses and in vivo mouse models. Our study identified key molecular correlates of the combination therapy and highlighted that anti-VEGF might synergize with anti-PD-L1 by targeting angiogenesis, Treg proliferation and myeloid cell inflammation.

Prediction of Major Histocompatibility Complex Binding with Bilateral and Variable Long Short Term Memory Networks

Simple Summary

Major histocompatibility complex molecules are of significant biological and clinical importance due to their utility in immunotherapy. The prediction of potential MHC binding peptides can estimate a T-cell immune response. The variable length of existing MHC binding peptides creates difficulty for MHC binding prediction algorithms. Thus, we utilized a bilateral and variable long-short term memory neural network to address this specific problem and developed a novel MHC binding prediction tool.

Abstract

As an important part of immune surveillance, major histocompatibility complex (MHC) is a set of proteins that recognize foreign molecules. Computational prediction methods for MHC binding peptides have been developed. However, existing methods share the limitation of fixed peptide sequence length, which necessitates the training of models by peptide length or prediction with a length reduction technique. Using a bidirectional long short-term memory neural network, we constructed BVMHC, an MHC class I and II binding prediction tool that is independent of peptide length. The performance of BVMHC was compared to seven MHC class I prediction tools and three MHC class II prediction tools using eight performance criteria independently. BVMHC attained the best performance in three of the eight criteria for MHC class I, and the best performance in four of the eight criteria for MHC class II, including accuracy and AUC. Furthermore, models for non-human species were also trained using the same strategy and made available for applications in mice, chimpanzees, macaques, and rats. BVMHC is composed of a series of peptide length independent MHC class I and II binding predictors. Models from this study have been implemented in an online web portal for easy access and use.

HLAB: learning the BiLSTM features from the ProtBert-encoded proteins for the class I HLA-peptide binding prediction

Human Leukocyte Antigen (HLA) is a type of molecule residing on the surfaces of most human cells and exerts an essential role in the immune system responding to the invasive items. The T cell antigen receptors may recognize the HLA-peptide complexes on the surfaces of cancer cells and destroy these cancer cells through toxic T lymphocytes. The computational determination of HLA-binding peptides will facilitate the rapid development of cancer immunotherapies. This study hypothesized that the natural language processing-encoded peptide features may be further enriched by another deep neural network. The hypothesis was tested with the Bi-directional Long Short-Term Memory-extracted features from the pretrained Protein Bidirectional Encoder Representations from Transformers-encoded features of the class I HLA (HLA-I)-binding peptides. The experimental data showed that our proposed HLAB feature engineering algorithm outperformed the existing ones in detecting the HLA-I-binding peptides. The extensive evaluation data show that the proposed HLAB algorithm outperforms all the seven existing studies on predicting the peptides binding to the HLA-A^*01:01 allele in AUC and achieves the best average AUC values on the six out of the seven k-mers (k=8,9,...,14, respectively represent the prediction task of a polypeptide consisting of k amino acids) except for the 9-mer prediction tasks. The source code and the fine-tuned feature extraction models are available at http://www.healthinformaticslab.org/supp/resources.php.

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).

Isolation of Neoantigen-Specific Human T Cell Receptors from Different Human and Murine Repertoires

Simple Summary

T cell-based immunotherapy has achieved remarkable clinical responses in patients with cancer. Neoepitope-specific T cells can specifically recognize mutated tumor cells and have led to tumor regression in mouse models and clinical studies. However, isolating neoepitope-specific T cell receptors (TCRs) from the patients’ own repertoire has shown limited success. Sourcing T cell repertoires, other than the patients’ own, has certain advantages: the availability of larger amounts of blood from healthy donors, circumventing tumor-related immunosuppression in patients, and including different donors to broaden the pool of specific T cells. Here, for the first time, a side-by-side comparison of three different TCR donor repertoires, including patients and HLA-matched allogenic healthy human repertoires, as well as repertoires of transgenic mice, is performed. Our results support recent studies that using not only healthy donor T cell repertoires, but also transgenic mice might be a viable strategy for isolating TCRs with known specificity directed against neoantigens for adoptive T cell therapy.

Abstract

(1) Background: Mutation-specific T cell receptor (TCR)-based adoptive T cell therapy represents a truly tumor-specific immunotherapeutic strategy. However, isolating neoepitope-specific TCRs remains a challenge. (2) Methods: We investigated, side by side, different TCR repertoires—patients’ peripheral lymphocytes (PBLs) and tumor-infiltrating lymphocytes (TILs), PBLs of healthy donors, and a humanized mouse model—to isolate neoepitope-specific TCRs against eight neoepitope candidates from a colon cancer and an ovarian cancer patient. Neoepitope candidates were used to stimulate T cells from different repertoires in vitro to generate neoepitope-specific T cells and isolate the specific TCRs. (3) Results: We isolated six TCRs from healthy donors, directed against four neoepitope candidates and one TCR from the murine T cell repertoire. Endogenous processing of one neoepitope, for which we isolated one TCR from both human and mouse-derived repertoires, could be shown. No neoepitope-specific TCR could be generated from the patients’ own repertoire. (4) Conclusion: Our data indicate that successful isolation of neoepitope-specific TCRs depends on various factors such as the heathy donor’s TCR repertoire or the presence of a tumor microenvironment allowing neoepitope-specific immune responses of the host. We show the advantage and feasibility of using healthy donor repertoires and humanized mouse TCR repertoires to generate mutation-specific TCRs with different specificities, especially in a setting when the availability of patient material is limited.