Structure Based Prediction of Neoantigen Immunogenicity

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.

Cancer neoepitopes viewed through negative selection and peripheral tolerance: a new path to cancer vaccines

A proportion of somatic mutations in tumors create neoepitopes that can prime T cell responses that target the MHC I-neoepitope complexes on tumor cells, mediating tumor control or rejection. Despite the compelling centrality of neoepitopes to cancer immunity, we know remarkably little about what constitutes a neoepitope that can mediate tumor control in vivo and what distinguishes such a neoepitope from the vast majority of similar candidate neoepitopes that are inefficacious in vivo. Studies in mice as well as clinical trials have begun to reveal the unexpected paradoxes in this area. Because cancer neoepitopes straddle that ambiguous ground between self and non-self, some rules that are fundamental to immunology of frankly non-self antigens, such as viral or model antigens, do not appear to apply to neoepitopes. Because neoepitopes are so similar to self-epitopes, with only small changes that render them non-self, immune response to them is regulated at least partially the way immune response to self is regulated. Therefore, neoepitopes are viewed and understood here through the clarifying lens of negative thymic selection. Here, the emergent questions in the biology and clinical applications of neoepitopes are discussed critically and a mechanistic and testable framework that explains the complexity and translational potential of these wonderful antigens is proposed.

Accurate modeling of peptide-MHC structures with AlphaFold

Major histocompatibility complex (MHC) proteins present peptides on the cell surface for T cell surveillance. Reliable in silico prediction of which peptides would be presented and which T cell receptors would recognize them is an important problem in structural immunology. Here, we introduce an AlphaFold-based pipeline for predicting the three-dimensional structures of peptide-MHC complexes for class I and class II MHC molecules. Our method demonstrates high accuracy, outperforming existing tools in class I modeling accuracy and class II peptide register prediction. We validate its performance and utility with new experimental data on a recently described cancer neoantigen/wild-type peptide pair and explore applications toward improving peptide-MHC binding prediction.

Development and use of machine learning algorithms in vaccine target selection

Computer-aided discovery of vaccine targets has become a cornerstone of rational vaccine design. In this article, I discuss how Machine Learning (ML) can inform and guide key computational steps in rational vaccine design concerned with the identification of B and T cell epitopes and correlates of protection. I provide examples of ML models, as well as types of data and predictions for which they are built. I argue that interpretable ML has the potential to improve the identification of immunogens also as a tool for scientific discovery, by helping elucidate the molecular processes underlying vaccine-induced immune responses. I outline the limitations and challenges in terms of data availability and method development that need to be addressed to bridge the gap between advances in ML predictions and their translational application to vaccine design.

Peptonics: A new family of cell-protecting surfactants for the recombinant expression of therapeutic proteins in mammalian cell cultures

Polymer surfactants are key components of cell culture media as they prevent mechanical damage during fermentation in stirred bioreactors. Among cell-protecting surfactants, Pluronics are widely utilized in biomanufacturing to ensure high cell viability and productivity. Monodispersity of monomer sequence and length is critical for the effectiveness of Pluronics-since minor deviations can damage the cells-but is challenging to achieve due to the stochastic nature of polymerization. Responding to this challenge, this study introduces Peptonics, a novel family of peptide and peptoid surfactants whose monomer composition and sequence are designed to achieve high cell viability and productivity at a fraction of chain length and cost of Pluronics. A designed ensemble of Peptonics was initially characterized via light scattering and tensiometry to select sequences whose phase behavior and tensioactivity align with those of Pluronics. Selected sequences were evaluated as cell-protecting surfactants using Chinese hamster ovary (CHO) cells expressing therapeutic monoclonal antibodies (mAb). Peptonics IH-T1010, ih-T1010, and ih-T1020 afforded high cell density (up to 3 × 107 cells mL-1 ) and viability (up to 95% within 10 days of culture), while reducing the accumulation of ammonia (a toxic metabolite) by ≈10% compared to Pluronic F-68. Improved cell viability afforded high mAb titer (up to 5.5 mg mL-1 ) and extended the production window beyond 14 days; notably, Peptonic IH-T1020 decreased mAb fragmentation and aggregation ≈5%, and lowered the titer of host cell proteins by 16% compared to Pluronic F-68. These features can improve significantly the purification of mAbs, thus increasing their availability at a lower cost to patients.

CD8+ T Cell-Dependent Antitumor Activity In Vivo of a Mass Spectrometry-Identified Neoepitope despite Undetectable CD8+ Immunogenicity In Vitro

Identification of neoepitopes that can control tumor growth in vivo remains a challenge even 10 y after the first genomics-defined cancer neoepitopes were identified. In this study, we identify a neoepitope, resulting from a mutation in the junction plakoglobin (Jup) gene (chromosome 11), from the mouse colon cancer line MC38-FABF (C57BL/6). This neoepitope, Jup mutant (JupMUT), was detected during mass spectrometry of MHC class I-eluted peptides from the tumor. JupMUT has a predicted binding affinity of 564 nM for the Kb molecule and a higher predicted affinity of 82 nM for Db. However, whereas structural modeling of JupMUT and its unmutated counterpart Jup wild-type indicates that there are little conformational differences between the two epitopes bound to Db, large structural divergences are predicted between the two epitopes bound to Kb. Together with in vitro binding data with RMA-S cells, these data suggest that Kb rather than Db is the relevant MHC class I molecule of JupMUT. Immunization of naive C57BL/6 mice with JupMUT elicits CD8-dependent tumor control of a MC38-FABF challenge. Despite the CD8 dependence of JupMUT-mediated tumor control in vivo, CD8+ T cells from JupMUT-immunized mice do not produce higher levels of IFN-γ than do naive mice. The structural and immunological characteristics of JupMUT are substantially different from those of many other neoepitopes that have been shown to mediate tumor control.

Unsupervised and supervised AI on molecular dynamics simulations reveals complex characteristics of HLA-A2-peptide immunogenicity

Immunologic recognition of peptide antigens bound to class I major histocompatibility complex (MHC) molecules is essential to both novel immunotherapeutic development and human health at large. Current methods for predicting antigen peptide immunogenicity rely primarily on simple sequence representations, which allow for some understanding of immunogenic features but provide inadequate consideration of the full scale of molecular mechanisms tied to peptide recognition. We here characterize contributions that unsupervised and supervised artificial intelligence (AI) methods can make toward understanding and predicting MHC(HLA-A2)-peptide complex immunogenicity when applied to large ensembles of molecular dynamics simulations. We first show that an unsupervised AI method allows us to identify subtle features that drive immunogenicity differences between a cancer neoantigen and its wild-type peptide counterpart. Next, we demonstrate that a supervised AI method for class I MHC(HLA-A2)-peptide complex classification significantly outperforms a sequence model on small datasets corrected for trivial sequence correlations. Furthermore, we show that both unsupervised and supervised approaches reveal determinants of immunogenicity based on time-dependent molecular fluctuations and anchor position dynamics outside the MHC binding groove. We discuss implications of these structural and dynamic immunogenicity correlates for the induction of T cell responses and therapeutic T cell receptor design.

Structural basis for T cell recognition of cancer neoantigens and implications for predicting neoepitope immunogenicity

Adoptive cell therapy (ACT) with tumor-specific T cells has been shown to mediate durable cancer regression. Tumor-specific T cells are also the basis of other therapies, notably cancer vaccines. The main target of tumor-specific T cells are neoantigens resulting from mutations in self-antigens over the course of malignant transformation. The detection of neoantigens presents a major challenge to T cells because of their high structural similarity to self-antigens, and the need to avoid autoimmunity. How different a neoantigen must be from its wild-type parent for it to induce a T cell response is poorly understood. Here we review recent structural and biophysical studies of T cell receptor (TCR) recognition of shared cancer neoantigens derived from oncogenes, including p53R175H, KRASG12D, KRASG12V, HHATp8F, and PIK3CAH1047L. These studies have revealed that, in some cases, the oncogenic mutation improves antigen presentation by strengthening peptide-MHC binding. In other cases, the mutation is detected by direct interactions with TCR, or by energetically driven or other indirect strategies not requiring direct TCR contacts with the mutation. We also review antibodies designed to recognize peptide-MHC on cell surfaces (TCR-mimic antibodies) as an alternative to TCRs for targeting cancer neoantigens. Finally, we review recent computational advances in this area, including efforts to predict neoepitope immunogenicity and how these efforts may be advanced by structural information on peptide-MHC binding and peptide-MHC recognition by TCRs.

Personalised neoantigen-based therapy in colorectal cancer

Colorectal cancer (CRC) has become one of the most common tumours with high morbidity, mortality and distinctive evolution mechanism. The neoantigens arising from the somatic mutations have become considerable treatment targets in the management of CRC. As cancer-specific aberrant peptides, neoantigens can trigger the robust host immune response and exert anti-tumour effects while minimising the emergence of adverse events commonly associated with alternative therapeutic regimens. In this review, we summarised the mechanism, generation, identification and prognostic significance of neoantigens, as well as therapeutic strategies challenges of neoantigen-based therapy in CRC. The evidence suggests that the establishment of personalised neoantigen-based therapy holds great promise as an effective treatment approach for patients with CRC.

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.

Accelerating antibody discovery and design with artificial intelligence: Recent advances and prospects

Therapeutic antibodies are the largest class of biotherapeutics and have been successful in treating human diseases. However, the design and discovery of antibody drugs remains challenging and time-consuming. Recently, artificial intelligence technology has had an incredible impact on antibody design and discovery, resulting in significant advances in antibody discovery, optimization, and developability. This review summarizes major machine learning (ML) methods and their applications for computational predictors of antibody structure and antigen interface/interaction, as well as the evaluation of antibody developability. Additionally, this review addresses the current status of ML-based therapeutic antibodies under preclinical and clinical phases. While many challenges remain, ML may offer a new therapeutic option for the future direction of fully computational antibody design.

Neoantigen identification: Technological advances and challenges

Neoantigens have emerged as promising targets for cutting-edge immunotherapies, such as cancer vaccines and adoptive cell therapy. These neoantigens are unique to tumors and arise exclusively from somatic mutations or non-genomic aberrations in tumor proteins. They encompass a wide range of alterations, including genomic mutations, post-transcriptomic variants, and viral oncoproteins. With the advancements in technology, the identification of immunogenic neoantigens has seen rapid progress, raising new opportunities for enhancing their clinical significance. Prediction of neoantigens necessitates the acquisition of high-quality samples and sequencing data, followed by mutation calling. Subsequently, the pipeline involves integrating various tools that can predict the expression, processing, binding, and recognition potential of neoantigens. However, the continuous improvement of computational tools is constrained by the availability of datasets which contain validated immunogenic neoantigens. This review article aims to provide a comprehensive summary of the current knowledge as well as limitations in neoantigen prediction and validation. Additionally, it delves into the origin and biological role of neoantigens, offering a deeper understanding of their significance in the field of cancer immunotherapy. This article thus seeks to contribute to the ongoing efforts to harness neoantigens as powerful weapons in the fight against cancer.

A transfer-learning approach to predict antigen immunogenicity and T-cell receptor specificity

Antigen immunogenicity and the specificity of binding of T-cell receptors to antigens are key properties underlying effective immune responses. Here we propose diffRBM, an approach based on transfer learning and Restricted Boltzmann Machines, to build sequence-based predictive models of these properties. DiffRBM is designed to learn the distinctive patterns in amino-acid composition that, on the one hand, underlie the antigen's probability of triggering a response, and on the other hand the T-cell receptor's ability to bind to a given antigen. We show that the patterns learnt by diffRBM allow us to predict putative contact sites of the antigen-receptor complex. We also discriminate immunogenic and non-immunogenic antigens, antigen-specific and generic receptors, reaching performances that compare favorably to existing sequence-based predictors of antigen immunogenicity and T-cell receptor specificity.

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

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

Unraveling tumor specific neoantigen immunogenicity prediction: a comprehensive analysis

Introduction

Identification of tumor specific neoantigen (TSN) immunogenicity is crucial to develop peptide/mRNA based anti-tumoral vaccines and/or adoptive T-cell immunotherapies; thus, accurate in-silico classification/prioritization proves critical for cost-effective clinical applications. Several methods were proposed as TSNs immunogenicity predictors; however, comprehensive performance comparison is still lacking due to the absence of well documented and adequate TSN databases.

Methods

Here, by developing a new curated database having 199 TSNs with experimentally-validated MHC-I presentation and positive/negative immune response (ITSNdb), sixteen metrics were evaluated as immunogenicity predictors. In addition, by using a dataset emulating patient derived TSNs and immunotherapy cohorts containing predicted TSNs for tumor neoantigen burden (TNB) with outcome association, the metrics were evaluated as TSNs prioritizers and as immunotherapy response biomarkers.

Results

Our results show high performance variability among methods, highlighting the need for substantial improvement. Deep learning predictors were top ranked on ITSNdb but show discrepancy on validation databases. In overall, current predicted TNB did not outperform existing biomarkers.

Conclusion

Recommendations for their clinical application and the ITSNdb are presented to promote development and comparison of computational TSNs immunogenicity predictors.

Usefulness of Docking and Molecular Dynamics in Selecting Tumor Neoantigens to Design Personalized Cancer Vaccines: A Proof of Concept

Personalized cancer vaccines based on neoantigens are a new and promising treatment for cancer; however, there are still multiple unresolved challenges to using this type of immunotherapy. Among these, the effective identification of immunogenic neoantigens stands out, since the in silico tools used generate a significant portion of false positives. Inclusion of molecular simulation techniques can refine the results these tools produce. In this work, we explored docking and molecular dynamics to study the association between the stability of peptide-HLA complexes and their immunogenicity, using as a proof of concept two HLA-A2-restricted neoantigens that were already evaluated in vitro. The results obtained were in accordance with the in vitro immunogenicity, since the immunogenic neoantigen ASTN1 remained bound at both ends to the HLA-A2 molecule. Additionally, molecular dynamic simulation suggests that position 1 of the peptide has a more relevant role in stabilizing the N-terminus than previously proposed. Likewise, the mutations may have a "delocalized" effect on the peptide-HLA interaction, which means that the mutated amino acid influences the intensity of the interactions of distant amino acids of the peptide with the HLA. These findings allow us to propose the inclusion of molecular simulation techniques to improve the identification of neoantigens for cancer vaccines.

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 peptide/HLA (pHLA, the human MHC) structures 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 a curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these representative backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer peptide/HLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in terms of structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work provide a framework for linking conformational diversity with antigen immunogenicity and receptor cross-reactivity.

Accurate modeling of peptide-MHC structures with AlphaFold

Major histocompatibility complex (MHC) proteins present peptides on the cell surface for T-cell surveillance. Reliable in silico prediction of which peptides would be presented and which T-cell receptors would recognize them is an important problem in structural immunology. Here, we introduce an AlphaFold-based pipeline for predicting the three-dimensional structures of peptide-MHC complexes for class I and class II MHC molecules. Our method demonstrates high accuracy, outperforming existing tools in class I modeling precision and class II peptide register prediction. We explore applications of this method towards improving peptide-MHC binding prediction.

Hidden hydrophobicity impacts polymer immunogenicity

Antibodies against poly(ethylene glycol) (PEG) have been found to be the culprit of side reactions and efficacy loss of a number of PEGylated drugs. Fundamental mechanisms of PEG immunogenicity and design principles for PEG alternatives still have not been fully explored. By using hydrophobic interaction chromatography (HIC) under varied salt conditions, we reveal the "hidden" hydrophobicity of those polymers which are generally considered as hydrophilic. A correlation between the hidden hydrophobicity of a polymer and its polymer immunogenicity is observed when this polymer is conjugated with an immunogenic protein. Such a correlation of hidden hydrophobicity vs. immunogenicity for a polymer also applies to corresponding polymer-protein conjugates. Atomistic molecular dynamics (MD) simulation results show a similar trend. Based on polyzwitterion modification and with this HIC technique, we are able to produce extremely low-immunogenic protein conjugates as their hydrophilicity is pushed to the limit and their hydrophobicity is eliminated, breaking the current barriers of eliminating anti-drug and anti-polymer antibodies.

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.

Mechanisms of resistance to immune checkpoint inhibitors in melanoma: What we have to overcome?

Marching into the second decade after the approval of ipilimumab, it is clear that immune checkpoint inhibitors (ICIs) have dramatically improved the prognosis of melanoma. Although the current edge is already high, with a 4-year OS% of 77.9% for adjuvant nivolumab and a 6.5-year OS% of 49% for nivolumab/ipilimumab combination in the metastatic setting, a high proportion of patients with advanced melanoma have no benefit from immunotherapy, or experience an early disease relapse/progression in the first few months of treatment, surviving much less. Reasonably, the primary and acquired resistance to ICIs has entered into the focus of clinical research with positive (e.g., nivolumab and relatlimab combination) and negative feedbacks (e.g., nivolumab with pegylated-IL2, pembrolizumab with T-VEC, nivolumab with epacadostat, and combinatorial triplets of BRAF/MEK inhibitors with immunotherapy). Many intrinsic (intracellular or intra-tumoral) but also extrinsic (systematic) events are considered to be involved in the development of this resistance to ICIs: i) melanoma cell immunogenicity (e.g., tumor mutational burden, antigen-processing machinery and immunogenic cell death, neoantigen affinity and heterogeneity, genomic instability, melanoma dedifferentiation and phenotypic plasticity), ii) immune cell trafficking, T-cell priming, and cell death evasion, iii) melanoma neovascularization, cellular TME components(e.g., T_regs, CAFs) and extracellular matrix modulation, iv) metabolic antagonism in the TME(highly glycolytic status, upregulated CD39/CD73/adenosine pathway, iDO-dependent tryptophan catabolism), v) T-cell exhaustion and negative immune checkpoints, and vi) gut microbiota. In the present overview, we discuss how these parameters compromise the efficacy of ICIs, with an emphasis on the lessons learned by the latest melanoma studies; and in parallel, we describe the main ongoing approaches to overcome the resistance to immunotherapy. Summarizing this information will improve the understanding of how these complicated dynamics contribute to immune escape and will help to develop more effective strategies on how anti-tumor immunity can surpass existing barriers of ICI-refractory melanoma.

InvitroSPI and a large database of proteasome-generated spliced and non-spliced peptides

Noncanonical epitopes presented by Human Leucocyte Antigen class I (HLA-I) complexes to CD8+ T cells attracted the spotlight in the research of novel immunotherapies against cancer, infection and autoimmunity. Proteasomes, which are the main producers of HLA-I-bound antigenic peptides, can catalyze both peptide hydrolysis and peptide splicing. The prediction of proteasome-generated spliced peptides is an objective that still requires a reliable (and large) database of non-spliced and spliced peptides produced by these proteases. Here, we present an extended database of proteasome-generated spliced and non-spliced peptides, which was obtained by analyzing in vitro digestions of 80 unique synthetic polypeptide substrates, measured by different mass spectrometers. Peptides were identified through invitroSPI method, which was validated through in silico and in vitro strategies. The peptide product database contains 16,631 unique peptide products (5,493 non-spliced, 6,453 cis-spliced and 4,685 trans-spliced peptide products), and a substrate sequence variety that is a valuable source for predictors of proteasome-catalyzed peptide hydrolysis and splicing. Potential artefacts and skewed results due to different identification and analysis strategies are discussed.

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.

Spatiotemporal evolution of the clear cell renal cell carcinoma microenvironment links intra-tumoral heterogeneity to immune escape

BACKGROUND

Intratumoral heterogeneity (ITH) is a hallmark of clear cell renal cell carcinoma (ccRCC) that reflects the trajectory of evolution and influences clinical prognosis. Here, we seek to elucidate how ITH and tumor evolution during immune checkpoint inhibitor (ICI) treatment can lead to therapy resistance.

METHODS

Here, we completed a single-arm pilot study to examine the safety and feasibility of neoadjuvant nivolumab in patients with localized RCC. Primary endpoints were safety and feasibility of neoadjuvant nivolumab. Then, we spatiotemporally profiled the genomic and immunophenotypic characteristics of 29 ccRCC patients, including pre- and post-therapy samples from 17 ICI-treated patients. Deep multi-regional whole-exome and transcriptome sequencing were performed on 29 patients at different time points before and after ICI therapy. T cell repertoire was also monitored from tissue and peripheral blood collected from a subset of patients to study T cell clonal expansion during ICI therapy.

RESULTS

Angiogenesis, lymphocytic infiltration, and myeloid infiltration varied significantly across regions of the same patient, potentially confounding their utility as biomarkers of ICI response. Elevated ITH associated with a constellation of both genomic features (HLA LOH, CDKN2A/B loss) and microenvironmental features, including elevated myeloid expression, reduced peripheral T cell receptor (TCR) diversity, and putative neoantigen depletion. Hypothesizing that ITH may itself play a role in shaping ICI response, we derived a transcriptomic signature associated with neoantigen depletion that strongly associated with response to ICI and targeted therapy treatment in several independent clinical trial cohorts.

CONCLUSIONS

These results argue that genetic and immune heterogeneity jointly co-evolve and influence response to ICI in ccRCC. Our findings have implications for future biomarker development for ICI response across ccRCC and other solid tumors and highlight important features of tumor evolution under ICI treatment.

TRIAL REGISTRATION

The study was registered on ClinicalTrial.gov (NCT02595918) on November 4, 2015.

Therapeutic high affinity T cell receptor targeting a KRASG12D cancer neoantigen

Neoantigens derived from somatic mutations are specific to cancer cells and are ideal targets for cancer immunotherapy. KRAS is the most frequently mutated oncogene and drives the pathogenesis of several cancers. Here we show the identification and development of an affinity-enhanced T cell receptor (TCR) that recognizes a peptide derived from the most common KRAS mutant, KRAS^G12D, presented in the context of HLA-A*11:01. The affinity of the engineered TCR is increased by over one million-fold yet fully able to distinguish KRAS^G12D over KRAS^WT. While crystal structures reveal few discernible differences in TCR interactions with KRAS^WT versus KRAS^G12D, thermodynamic analysis and molecular dynamics simulations reveal that TCR specificity is driven by differences in indirect electrostatic interactions. The affinity enhanced TCR, fused to a humanized anti-CD3 scFv, enables selective killing of cancer cells expressing KRAS^G12D. Our work thus reveals a molecular mechanism that drives TCR selectivity and describes a soluble bispecific molecule with therapeutic potential against cancers harboring a common shared neoantigen.

Cancers often harbor mutations in genes encoding important regulatory proteins, but therapeutic targeting of these molecules proves difficult due to their high structural similarity to their non-mutated counterpart. Here authors show the engineering of T cell engaging bispecific protein able to selectively target cancer cells with a high-frequency mutation in the KRAS oncogene.

Recent Advances and Challenges in Cancer Immunotherapy

Simple Summary

Immunotherapy helps a person’s immune system to target tumor cells. Recent advances in cancer immunotherapy, including immune checkpoint inhibition, chimeric antigen receptor T-cell therapy and cancer vaccination, have changed the landscape of cancer treatment. These approaches have had profound success in certain cancer types but still fail in the majority of cases. This review will cover both successes and current challenges in cancer immunotherapy, as well as recent advances in the field of basic tumor immunology that will allow us to overcome resistance to existing treatments.

Abstract

Cancer immunotherapy has revolutionized the field of oncology in recent years. Harnessing the immune system to treat cancer has led to a large growth in the number of novel immunotherapeutic strategies, including immune checkpoint inhibition, chimeric antigen receptor T-cell therapy and cancer vaccination. In this review, we will discuss the current landscape of immuno-oncology research, with a focus on elements that influence immunotherapeutic outcomes. We will also highlight recent advances in basic aspects of tumor immunology, in particular, the role of the immunosuppressive cells within the tumor microenvironment in regulating antitumor immunity. Lastly, we will discuss how the understanding of basic tumor immunology can lead to the development of new immunotherapeutic strategies.

Functional landscapes of POLE and POLD1 mutations in checkpoint blockade-dependent antitumor immunity

Defects in pathways governing genomic fidelity have been linked to improved response to immune checkpoint blockade therapy (ICB). Pathogenic POLE/POLD1 mutations can cause hypermutation, yet how diverse mutations in POLE/POLD1 influence antitumor immunity following ICB is unclear. Here, we comprehensively determined the effect of POLE/POLD1 mutations in ICB and elucidated the mechanistic impact of these mutations on tumor immunity. Murine syngeneic tumors harboring Pole/Pold1 functional mutations displayed enhanced antitumor immunity and were sensitive to ICB. Patients with POLE/POLD1 mutated tumors harboring telltale mutational signatures respond better to ICB than patients harboring wild-type or signature-negative tumors. A mutant POLE/D1 function-associated signature-based model outperformed several traditional approaches for identifying POLE/POLD1 mutated patients that benefit from ICB. Strikingly, the spectrum of mutational signatures correlates with the biochemical features of neoantigens. Alterations that cause POLE/POLD1 function-associated signatures generate T cell receptor (TCR)-contact residues with increased hydrophobicity, potentially facilitating T cell recognition. Altogether, the functional landscapes of POLE/POLD1 mutations shape immunotherapy efficacy.

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

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

CAD v1.0: Cancer Antigens Database Platform for Cancer Antigen Algorithm Development and Information Exploration

Cancer vaccines have gradually attracted attention for their tremendous preclinical and clinical performance. With the development of next-generation sequencing technologies and related algorithms, pipelines based on sequencing and machine learning methods have become mainstream in cancer antigen prediction; of particular focus are neoantigens, mutation peptides that only exist in tumor cells that lack central tolerance and have fewer side effects. The rapid prediction and filtering of neoantigen peptides are crucial to the development of neoantigen-based cancer vaccines. However, due to the lack of verified neoantigen datasets and insufficient research on the properties of neoantigens, neoantigen prediction algorithms still need to be improved. Here, we recruited verified cancer antigen peptides and collected as much relevant peptide information as possible. Then, we discussed the role of each dataset for algorithm improvement in cancer antigen research, especially neoantigen prediction. A platform, Cancer Antigens Database (CAD, http://cad.bio-it.cn/), was designed to facilitate users to perform a complete exploration of cancer antigens online.

Evaluating performance of existing computational models in predicting CD8+ T cell pathogenic epitopes and cancer neoantigens

T cell recognition of a cognate peptide-major histocompatibility complex (pMHC) presented on the surface of infected or malignant cells is of the utmost importance for mediating robust and long-term immune responses. Accurate predictions of cognate pMHC targets for T cell receptors would greatly facilitate identification of vaccine targets for both pathogenic diseases and personalized cancer immunotherapies. Predicting immunogenic peptides therefore has been at the center of intensive research for the past decades but has proven challenging. Although numerous models have been proposed, performance of these models has not been systematically evaluated and their success rate in predicting epitopes in the context of human pathology has not been measured and compared. In this study, we evaluated the performance of several publicly available models, in identifying immunogenic CD8+ T cell targets in the context of pathogens and cancers. We found that for predicting immunogenic peptides from an emerging virus such as severe acute respiratory syndrome coronavirus 2, none of the models perform substantially better than random or offer considerable improvement beyond HLA ligand prediction. We also observed suboptimal performance for predicting cancer neoantigens. Through investigation of potential factors associated with ill performance of models, we highlight several data- and model-associated issues. In particular, we observed that cross-HLA variation in the distribution of immunogenic and non-immunogenic peptides in the training data of the models seems to substantially confound the predictions. We additionally compared key parameters associated with immunogenicity between pathogenic peptides and cancer neoantigens and observed evidence for differences in the thresholds of binding affinity and stability, which suggested the need to modulate different features in identifying immunogenic pathogen versus cancer peptides. Overall, we demonstrate that accurate and reliable predictions of immunogenic CD8+ T cell targets remain unsolved; thus, we hope our work will guide users and model developers regarding potential pitfalls and unsettled questions in existing immunogenicity predictors.

The Importance of Being Presented: Target Validation by Immunopeptidomics for Epitope-Specific Immunotherapies

Presentation of tumor-specific or tumor-associated peptides by HLA class I molecules to CD8+ T cells is the foundation of epitope-centric cancer immunotherapies. While often in silico HLA binding predictions or in vitro immunogenicity assays are utilized to select candidates, mass spectrometry-based immunopeptidomics is currently the only method providing a direct proof of actual cell surface presentation. Despite much progress in the last decade, identification of such HLA-presented peptides remains challenging. Here we review typical workflows and current developments in the field of immunopeptidomics, highlight the challenges which remain to be solved and emphasize the importance of direct target validation for clinical immunotherapy development.

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.

Identification of neoantigens for individualized therapeutic cancer vaccines

Somatic mutations in cancer cells can generate tumour-specific neoepitopes, which are recognized by autologous T cells in the host. As neoepitopes are not subject to central immune tolerance and are not expressed in healthy tissues, they are attractive targets for therapeutic cancer vaccines. Because the vast majority of cancer mutations are unique to the individual patient, harnessing the full potential of this rich source of targets requires individualized treatment approaches. Many computational algorithms and machine-learning tools have been developed to identify mutations in sequence data, to prioritize those that are more likely to be recognized by T cells and to design tailored vaccines for every patient. In this Review, we fill the gaps between the understanding of basic mechanisms of T cell recognition of neoantigens and the computational approaches for discovery of somatic mutations and neoantigen prediction for cancer immunotherapy. We present a new classification of neoantigens, distinguishing between guarding, restrained and ignored neoantigens, based on how they confer proficient antitumour immunity in a given clinical context. Such context-based differentiation will contribute to a framework that connects neoantigen biology to the clinical setting and medical peculiarities of cancer, and will enable future neoantigen-based therapies to provide greater clinical benefit.

T Cell Recognition of Tumor Neoantigens and Insights Into T Cell Immunotherapy

In cancer, non-synonymous DNA base changes alter protein sequence and produce neoantigens that are detected by the immune system. For immune detection, neoantigens must first be presented on class I or II human leukocyte antigens (HLA) followed by recognition by peptide-specific receptors, exemplified by the T-cell receptor (TCR). Detection of neoantigens represents a unique challenge to the immune system due to their high similarity with endogenous 'self' proteins. Here, we review insights into how TCRs detect neoantigens from structural studies and delineate two broad mechanistic categories: 1) recognition of mutated 'self' peptides and 2) recognition of novel 'non-self' peptides generated through anchor residue modifications. While mutated 'self' peptides differ only by a single amino acid from an existing 'self' epitope, mutations that form anchor residues generate an entirely new epitope, hitherto unknown to the immune system. We review recent structural studies that highlight these structurally distinct mechanisms and discuss how they may lead to differential anti-tumor immune responses. We discuss how T cells specific for neoantigens derived from anchor mutations can be of high affinity and provide insights to their use in adoptive T cell transfer-based immunotherapy.

Beyond Sequencing: Prioritizing and Delivering Neoantigens for Cancer Vaccines

Neoantigens are tumor-specific proteins and peptides that can be highly immunogenic. Immune-mediated tumor rejection is strongly associated with cytotoxic responses to neoantigen-derived peptides in noncovalent association with self-HLA molecules. Neoantigen-based therapies, such as adoptive T cell transfer, have shown the potential to induce remission of treatment-resistant metastatic disease in select patients. Cancer vaccines are similarly designed to elicit or amplify antigen-specific T cell populations and stimulate directed antitumor immunity, but the selection and prioritization of the neoantigens remains a challenge. Bioinformatic algorithms can predict tumor neoantigens from somatic mutations, insertion-deletions, and other aberrant peptide products, but this often leads to hundreds of potential neoepitopes, all unique for that tumor. Selecting neoantigens for cancer vaccines is complicated by the technical challenges of neoepitope discovery, the diversity of HLA molecules, and intratumoral heterogeneity of passenger mutations leading to immune escape. Despite strong preclinical evidence, few neoantigen cancer vaccines tested in vivo have generated epitope-specific T cell populations, suggesting suboptimal immune system activation. In this chapter, we review factors affecting the prioritization and delivery of candidate neoantigens in the design of therapeutic and preventive cancer vaccines and consider synergism with standard chemotherapies.

Self-mediated positive selection of T cells sets an obstacle to the recognition of nonself

Adaptive immune recognition is mediated by the binding of peptide-human leukocyte antigen complexes by T cells. Positive selection of T cells in the thymus is a fundamental step in the generation of a responding T cell repertoire: only those T cells survive that recognize human peptides presented on the surface of cortical thymic epithelial cells. We propose that while this step is essential for optimal immune function, the process results in a defective T cell repertoire because it is mediated by self-peptides. To test our hypothesis, we focused on amino acid motifs of peptides in contact with T cell receptors. We found that motifs rarely or not found in the human proteome are unlikely to be recognized by the immune system just like the ones that are not expressed in cortical thymic epithelial cells or not presented on their surface. Peptides carrying such motifs were especially dissimilar to human proteins. Importantly, we present our main findings on two independent T cell activation datasets and directly demonstrate the absence of naïve T cells in the repertoire of healthy individuals. We also show that T cell cross-reactivity is unable to compensate for the absence of positively selected T cells. Additionally, we show that the proposed mechanism could influence the risk for different infectious diseases. In sum, our results suggest a side effect of T cell positive selection, which could explain the nonresponsiveness to many nonself peptides and could improve the understanding of adaptive immune recognition.

Toward Guided Mutagenesis: Gaussian Process Regression Predicts MHC Class II Antigen Mutant Binding

Antigen-specific immunotherapies (ASI) require successful loading and presentation of antigen peptides into the major histocompatibility complex (MHC) binding cleft. One route of ASI design is to mutate native antigens for either stronger or weaker binding interaction to MHC. Exploring all possible mutations is costly both experimentally and computationally. To reduce experimental and computational expense, here we investigate the minimal amount of prior data required to accurately predict the relative binding affinity of point mutations for peptide-MHC class II (pMHCII) binding. Using data from different residue subsets, we interpolate pMHCII mutant binding affinities by Gaussian process (GP) regression of residue volume and hydrophobicity. We apply GP regression to an experimental data set from the Immune Epitope Database, and theoretical data sets from NetMHCIIpan and Free Energy Perturbation calculations. We find that GP regression can predict binding affinities of nine neutral residues from a six-residue subset with an average R² coefficient of determination value of 0.62 ± 0.04 (±95% CI), average error of 0.09 ± 0.01 kcal/mol (±95% CI), and with an receiver operating characteristic (ROC) AUC value of 0.92 for binary classification of enhanced or diminished binding affinity. Similarly, metrics increase to an R² value of 0.69 ± 0.04, average error of 0.07 ± 0.01 kcal/mol, and an ROC AUC value of 0.94 for predicting seven neutral residues from an eight-residue subset. Our work finds that prediction is most accurate for neutral residues at anchor residue sites without register shift. This work holds relevance to predicting pMHCII binding and accelerating ASI design.

Leveraging Immunopeptidomics To Study and Combat Infectious Disease

T cells must recognize pathogen-derived peptides bound to major histocompatibility complexes (MHCs) in order to initiate a cell-mediated immune response against an infection, or to support the development of high-affinity antibody responses. Identifying antigens presented on MHCs by infected cells and professional antigen-presenting cells (APCs) during infection may therefore provide a route toward developing new vaccines. Peptides bound to MHCs can be identified at whole-proteome scale using mass spectrometry-a technique referred to as "immunopeptidomics." This technique has emerged as a powerful tool for identifying potential vaccine targets in the context of many infectious diseases. In this review, we discuss the contributions immunopeptidomic studies have made to understanding antigen presentation and T cell priming in the context of infection and the potential for immunopeptidomics to inform the development of vaccines to address pressing global health problems in infectious disease.

neoDL: a novel neoantigen intrinsic feature-based deep learning model identifies IDH wild-type glioblastomas with the longest survival

Background

Neoantigen based personalized immune therapies achieve promising results in melanoma and lung cancer, but few neoantigen based models perform well in IDH wild-type GBM, and the association between neoantigen intrinsic features and prognosis remain unclear in IDH wild-type GBM. We presented a novel neoantigen intrinsic feature-based deep learning model (neoDL) to stratify IDH wild-type GBMs into subgroups with different survivals.

Results

We first derived intrinsic features for each neoantigen associated with survival, followed by applying neoDL in TCGA data cohort(AUC = 0.988, p value < 0.0001). Leave one out cross validation (LOOCV) in TCGA demonstrated that neoDL successfully classified IDH wild-type GBMs into different prognostic subgroups, which was further validated in an independent data cohort from Asian population. Long-term survival IDH wild-type GBMs identified by neoDL were found characterized by 12 protective neoantigen intrinsic features and enriched in development and cell cycle.

Conclusions

The model can be therapeutically exploited to identify IDH wild-type GBM with good prognosis who will most likely benefit from neoantigen based personalized immunetherapy. Furthermore, the prognostic intrinsic features of the neoantigens inferred from this study can be used for identifying neoantigens with high potentials of immunogenicity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04301-6.

DeepAntigen: a novel method for neoantigen prioritization via 3D genome and deep sparse learning

MOTIVATION

The mutations of cancers can encode the seeds of their own destruction, in the form of T-cell recognizable immunogenic peptides, also known as neoantigens. It is computationally challenging, however, to accurately prioritize the potential neoantigen candidates according to their ability of activating the T-cell immunoresponse, especially when the somatic mutations are abundant. Although a few neoantigen prioritization methods have been proposed to address this issue, advanced machine learning model that is specifically designed to tackle this problem is still lacking. Moreover, none of the existing methods considers the original DNA loci of the neoantigens in the perspective of 3D genome which may provide key information for inferring neoantigens' immunogenicity.

RESULTS

In this study, we discovered that DNA loci of the immunopositive and immunonegative MHC-I neoantigens have distinct spatial distribution patterns across the genome. We therefore used the 3D genome information along with an ensemble pMHC-I coding strategy, and developed a group feature selection-based deep sparse neural network model (DNN-GFS) that is optimized for neoantigen prioritization. DNN-GFS demonstrated increased neoantigen prioritization power comparing to existing sequence-based approaches. We also developed a webserver named deepAntigen (http://yishi.sjtu.edu.cn/deepAntigen) that implements the DNN-GFS as well as other machine learning methods. We believe that this work provides a new perspective toward more accurate neoantigen prediction which eventually contribute to personalized cancer immunotherapy.

AVAILABILITY AND IMPLEMENTATION

Data and implementation are available on webserver: http://yishi.sjtu.edu.cn/deepAntigen.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

ERAP2 Increases the Abundance of a Peptide Submotif Highly Selective for the Birdshot Uveitis-Associated HLA-A29

Birdshot Uveitis (BU) is a blinding inflammatory eye condition that only affects HLA-A29-positive individuals. Genetic association studies linked ERAP2 with BU, an aminopeptidase which trims peptides before their presentation by HLA class I at the cell surface, which suggests that ERAP2-dependent peptide presentation by HLA-A29 drives the pathogenesis of BU. However, it remains poorly understood whether the effects of ERAP2 on the HLA-A29 peptidome are distinct from its effect on other HLA allotypes. To address this, we focused on the effects of ERAP2 on the immunopeptidome in patient-derived antigen presenting cells. Using complementary HLA-A29-based and pan-class I immunopurifications, isotope-labeled naturally processed and presented HLA-bound peptides were sequenced by mass spectrometry. We show that the effects of ERAP2 on the N-terminus of ligands of HLA-A29 are shared across endogenous HLA allotypes, but discover and replicate that one peptide motif generated in the presence of ERAP2 is specifically bound by HLA-A29. This motif can be found in the amino acid sequence of putative autoantigens. We further show evidence for internal sequence specificity for ERAP2 imprinted in the immunopeptidome. These results reveal that ERAP2 can generate an HLA-A29-specific antigen repertoire, which supports that antigen presentation is a key disease pathway in BU.

Characterization of CD8+ T-cell responses to non-anchor-type HLA class I neoantigens with single amino-acid substitutions

CD8⁺ T cells are capable of recognizing mutation-derived neoantigens displayed by HLA class I molecules, thereby exhibiting the ability to distinguish between cancer and normal cells. However, accumulating evidence has shown that only a small fraction of nonsynonymous somatic mutations give rise to clinically relevant neoantigens. The properties of such neoantigens, which must be presented by HLA and immunogenic to induce a T-cell response, remain elusive. In this study, we explored the HLA class I ligandome of a human cancer cell line with microsatellite instability using a proteogenomic approach. The results demonstrated that neoantigens accounted for only 0.34% of the HLA class I ligandome, and most neoantigens were encoded by genes with abundant expression. Thereafter, T-cell responses were prioritized, and immunodominant neoantigens were defined using naive CD8⁺ T cells derived from healthy donors. AKF9, an immunogenic neoantigen with a mutation at a non-anchor position, formed a stable peptide-HLA complex. T-cell responses were analyzed against a panel of AKF9 variants with single amino-acid substitutions, in which mutations did not alter the high HLA-binding affinity and stability. The responses varied across individuals, demonstrating the impact of heterogeneous T-cell repertoires in this human cancer model. Moreover, responses were biased toward a variant group with large structural changes compared to the wild-type peptide. Thus, naive T-cell induction can be attributed to multiple determinants. Combining structural dissimilarity with gene-expression levels, HLA-binding affinity, and stability may further help prioritize the immunogenicity of non-anchor-type neoantigens.

Qualification of tumour mutational burden by targeted next-generation sequencing as a biomarker in hepatocellular carcinoma

BACKGROUND & AIMS

Tumour mutational burden (TMB) predicts improved response and survival to immunotherapy. In this pilot study, we optimized targeted next-generation sequencing (tNGS) to estimate TMB in hepatocellular carcinoma (HCC).

METHODS

We sequenced 48 non-paired samples (21 fresh-frozen [FF] and 27 paraffin-embedded [FFPE]), among which 11 FFPE samples were pretreated with uracil-DNA glycosylase (UDG). Thirty samples satisfied post-sequencing quality control. High/low TMB was defined by median number of mutations/Mb (Mut/Mb), across different minimum allele frequency (MAF) thresholds (≥0.05, ≥0.1 and ≥0.2).

RESULTS

Eligible patients (n = 29) were cirrhotic (84%) with TNM stage I-II HCC (75%). FFPE samples had higher TMB (median 958.39 vs 2.51 Mut/Mb, P < .0001), estimated deamination counts (median 1335.50 vs 0, P < .0001) and C > T transitions at CpG sites (median 60.3% vs 9.1%, P = .002) compared to FF. UDG-treated samples had lower TMB (median 4019.92 vs 353 Mut/Mb, P = .041) and deamination counts (median 6393.5 vs 328.5, P = .041) vs untreated FFPE. At 0.2 MAF threshold with UDG treatment, median TMB was 5.48 (range 1.68-16.07) and did not correlate with salient pathologic features of HCC, including survival.

CONCLUSION

While tNGS on fresh HCC samples appears to be the optimal source of tumour DNA, the low median TMB values observed may limit the role of TMB as a predictor of response to immunotherapy in HCC.

High-throughput modeling and scoring of TCR-pMHC complexes to predict cross-reactive peptides

MOTIVATION

The binding of T cell receptors (TCRs) to their target peptide MHC (pMHC) ligands initializes the cell-mediated immune response. In autoimmune diseases such as multiple sclerosis, the TCR erroneously recognizes self-peptides as foreign and activates an immune response against healthy cells. Such responses can be triggered by cross-recognition of the autoreactive TCR with foreign peptides. Hence, it would be desirable to identify such foreign-antigen triggers to provide a mechanistic understanding of autoimmune diseases. However, the large sequence space of foreign antigens presents an obstacle in the identification of cross-reactive peptides.

RESULTS

Here, we present an in silico modeling and scoring method which exploits the structural properties of TCR-pMHC complexes to predict the binding of cross-reactive peptides. We analyzed three mouse TCRs and one human TCR isolated from a patient with multiple sclerosis. Cross-reactive peptides for these TCRs were previously identified via yeast display coupled with deep sequencing, providing a robust dataset for evaluating our method. Modeling query peptides in their associated TCR-pMHC crystal structures, our method accurately selected the top binding peptides from sets containing more than a hundred thousand unique peptides.

AVAILABILITY AND IMPLEMENTATION

Analyses were performed using custom Python and R scripts available at https://github.com/tborrman/antigen-predict.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Structure-Based Modeling of SARS-CoV-2 Peptide/HLA-A02 Antigens

SARS-CoV-2-specific CD4 and CD8 T cells have been shown to be present in individuals with acute, mild, and asymptomatic Coronavirus disease (COVID-19). Toward the development of diagnostic and therapeutic tools to fight COVID-19, it is important to predict and characterize T cell epitopes expressed by SARS-CoV-2. Here, we use RosettaMHC, a comparative modeling approach which leverages existing structures of peptide/MHC complexes available in the Protein Data Bank, to derive accurate 3D models for putative SARS-CoV-2 CD8 epitopes. We outline an application of our method to model 8-10 residue epitopic peptides predicted to bind to the common allele HLA-A*02:01, and we make our models publicly available through an online database (https://rosettamhc.chemistry.ucsc.edu). We further compare electrostatic surfaces with models of homologous peptide/HLA-A*02:01 complexes from human common cold coronavirus strains to identify epitopes which may be recognized by a shared pool of cross-reactive TCRs. As more detailed studies on antigen-specific T cell recognition become available, RosettaMHC models can be used to understand the link between peptide/HLA complex structure and surface chemistry with immunogenicity, in the context of SARS-CoV-2 infection.

Predicting Cross-Reactivity and Antigen Specificity of T Cell Receptors

Adaptive immune recognition is mediated by specific interactions between heterodimeric T cell receptors (TCRs) and their cognate peptide-MHC (pMHC) ligands, and the methods to accurately predict TCR:pMHC interaction would have profound clinical, therapeutic and pharmaceutical applications. Herein, we review recent developments in predicting cross-reactivity and antigen specificity of TCR recognition. We discuss current experimental and computational approaches to investigate cross-reactivity and antigen-specificity of TCRs and highlight how integrating kinetic, biophysical and structural features may offer valuable insights in modeling immunogenicity. We further underscore the close inter-relationship of these two interconnected notions and the need to investigate each in the light of the other for a better understanding of T cell responsiveness for the effective clinical applications.

Main Strategies for the Identification of Neoantigens

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

Point mutation screening of tumor neoantigens and peptide-induced specific cytotoxic T lymphocytes using The Cancer Genome Atlas database

The aim of the present study was to use The Cancer Genome Atlas (TCGA) database to identify tumor neoantigens, combined with a bioinformatics analysis to design and analyze antigen epitope peptides. Epitopes were screened using immunogenicity tests to identify the ideal epitope peptides to target tumor neoantigens, which can specifically activate the immune response of T cells. The high-frequency mutation loci (top 10) of colorectal, lung and liver cancer genes were screened using TCGA database. The antigenic epitope peptides with high affinity for major histocompatibility complex molecules were selected and synthesized using computer prediction algorithms, and were subsequently detected using flow cytometry. The cytotoxicity of specific cytotoxic T lymphocytes (CTLs) on peptide-loaded T2 cells was initially verified using interferon IFN-γ detection and a calcein-acetoxymethyl (Cal-AM) release assay. Tumor cell lines expressing point mutations in KRAS, TP53 and CTNNB1 genes were constructed respectively, and the cytotoxicity of peptide-induced specific CTLs on wild-type and mutant tumor cells was verified using a Cal-AM release assay and carboxyfluorescein succinimidyl ester-propidium iodide staining. The high-frequency gene mutation loci of KRAS proto-oncogene (KRAS) G12V, tumor protein p53 (TP53) R158L and catenin β1 (CTNNB1) K335I were identified in TCGA database. A total of 3 groups of wild-type and mutant peptides were screened using a peptide prediction algorithm. The CTNNB1 group had a strong affinity for the human leukocyte antigen-A2 molecule, as determined using flow cytometry. The IFN-γ secretion of specific CTLs in the CTNNB1 group was the highest, followed by the TP53 and the KRAS groups. The killing rate of mutant peptide-induced specific CTLs on peptide-loaded T2 cells in the CTNNB1 group was higher compared with that observed in the other groups. The killing rate of specific CTLs induced by mutant peptides present on tumor cells was higher compared with that induced by wild-type peptides. However, when compared with the TP53 and KRAS groups, specific CTLs induced by mutant peptides in the CTNNB1 group had more potent cytotoxicity towards mutant and wild-type tumor cells. In conclusion, point mutant tumor neoantigens screened in the three groups improved the cytotoxicity of specific T cells, and the mutant peptides in the CTNNB1 group were more prominent, indicating that they may activate the cellular immune response more readily.

Role of in silico structural modeling in predicting immunogenic neoepitopes for cancer vaccine development

In prior studies, we delineated the landscape of neoantigens arising from nonsynonymous point mutations in a murine pancreatic cancer model, Panc02. We developed a peptide vaccine by targeting neoantigens predicted using a pipeline that incorporates the MHC binding algorithm NetMHC. The vaccine, when combined with immune checkpoint modulators, elicited a robust neoepitope-specific antitumor immune response and led to tumor clearance. However, only a small fraction of the predicted neoepitopes induced T cell immunity, similarly to that reported for neoantigen vaccines tested in clinical studies. While these studies have used binding affinities to MHC I as surrogates for T cell immunity, this approach does not include spatial information on the mutated residue that is crucial for TCR activation. Here, we investigate conformational alterations in and around the MHC binding groove induced by selected minimal neoepitopes, and we examine the influence of a given mutated residue as a function of its spatial position. We found that structural parameters, including the solvent-accessible surface area (SASA) of the neoepitope and the position and spatial configuration of the mutated residue within the sequence, can be used to improve the prediction of immunogenic neoepitopes for inclusion in cancer vaccines.

Structural parameters, including the solvent exposed surface area of the neoepitope and the position and spatial configuration of the mutated residue can be used to improve the prediction of immunogenic neoepitopes for inclusion in cancer vaccines.