Footprints of antigen processing boost MHC class II natural ligand predictions

3D genome contributes to MHC-II neoantigen prediction

Reliable and ultra-fast DNA and RNA sequencing have been achieved with the emergence of high-throughput sequencing technology. When combining the results of DNA and RNA sequencing for tumor cells of cancer patients, neoantigens that potentially stimulate the immune response of either CD4+ or CD8+ T cells can be identified. However, due to the abundance of somatic mutations and the high polymorphic nature of human leukocyte antigen (HLA) it is challenging to accurately predict the neoantigens. Moreover, comparing to HLA-I presented peptides, the HLA-II presented peptides are more variable in length, making the prediction of HLA-II loaded neoantigens even harder. A number of computational approaches have been proposed to address this issue but none of them considers the DNA origin of the neoantigens from the perspective of 3D genome. Here we investigate the DNA origins of the immune-positive and non-negative HLA-II neoantigens in the context of 3D genome and discovered that the chromatin 3D architecture plays an important role in more effective HLA-II neoantigen prediction. We believe that the 3D genome information will help to increase the precision of HLA-II neoantigen discovery and eventually benefit precision and personalized medicine in cancer immunotherapy.

Comprehensive approaches to preclinical evaluation of monoclonal antibodies and their next-generation derivatives

Biotherapeutics hold great promise for the treatment of several diseases and offer innovative possibilities for new treatments that target previously unaddressed medical needs. Despite successful transitions from preclinical to clinical stages and regulatory approval, there are instances where adverse reactions arise, resulting in product withdrawals. As a result, it is essential to conduct thorough evaluations of safety and effectiveness on an individual basis. This article explores current practices, challenges, and future approaches in conducting comprehensive preclinical assessments to ensure the safety and efficacy of biotherapeutics including monoclonal antibodies, toxin-conjugates, bispecific antibodies, single-chain antibodies, Fc-engineered antibodies, antibody mimetics, and siRNA-antibody/peptide conjugates.

New light on the HLA-DR immunopeptidomic landscape

The set of peptides processed and presented by major histocompatibility complex class II molecules defines the immunopeptidome, and its characterization holds keys to understanding essential properties of the immune system. High-throughput mass spectrometry (MS) techniques enable interrogation of the diversity and complexity of the immunopeptidome at an unprecedented scale. Here, we analyzed a large set of MS immunopeptidomics data from 40 donors, 221 samples, covering 30 unique HLA-DR molecules. We identified likely co-immunoprecipitated HLA-DR irrelevant contaminants using state-of-the-art prediction methods and unveiled novel light on the properties of HLA antigen processing and presentation. The ligandome (HLA binders) was enriched in 15-mer peptides, and the contaminome (nonbinders) in longer peptides. Classification of singletons and nested sets showed that the first were enriched in contaminants. Investigating the source protein location of ligands revealed that only contaminants shared a positional bias. Regarding subcellular localization, nested peptides were found to be predominantly of endolysosomal origin, whereas singletons shared an equal distribution between the cytosolic and endolysosomal origin. According to antigen-processing signatures, no significant differences were observed between the cytosolic and endolysosomal ligands. Further, the sensitivity of MS immunopeptidomics was investigated by analyzing overlap and saturation between biological MS replicas, concluding that at least 5 replicas are needed to identify 80% of the immunopeptidome. Moreover, the overlap in immunopeptidome between donors was found to be very low both in terms of peptides and source proteins, the latter indicating a critical HLA bias in the antigen sampling in the HLA antigen presentation. Finally, the complementarity between MS and in silico approaches for comprehensively sampling the immunopeptidome was demonstrated.

Graph-pMHC: graph neural network approach to MHC class II peptide presentation and antibody immunogenicity

Antigen presentation on MHC class II (pMHCII presentation) plays an essential role in the adaptive immune response to extracellular pathogens and cancerous cells. But it can also reduce the efficacy of large-molecule drugs by triggering an anti-drug response. Significant progress has been made in pMHCII presentation modeling due to the collection of large-scale pMHC mass spectrometry datasets (ligandomes) and advances in machine learning. Here, we develop graph-pMHC, a graph neural network approach to predict pMHCII presentation. We derive adjacency matrices for pMHCII using Alphafold2-multimer and address the peptide-MHC binding groove alignment problem with a simple graph enumeration strategy. We demonstrate that graph-pMHC dramatically outperforms methods with suboptimal inductive biases, such as the multilayer-perceptron-based NetMHCIIpan-4.0 (+20.17% absolute average precision). Finally, we create an antibody drug immunogenicity dataset from clinical trial data and develop a method for measuring anti-antibody immunogenicity risk using pMHCII presentation models. Our model increases receiver operating characteristic curve (ROC)-area under the ROC curve (AUC) by 2.57% compared to just filtering peptides by hits in OASis alone for predicting antibody drug immunogenicity.

Structure-aware deep model for MHC-II peptide binding affinity prediction

The prediction of major histocompatibility complex (MHC)-peptide binding affinity is an important branch in immune bioinformatics, especially helpful in accelerating the design of disease vaccines and immunity therapy. Although deep learning-based solutions have yielded promising results on MHC-II molecules in recent years, these methods ignored structure knowledge from each peptide when employing the deep neural network models. Each peptide sequence has its specific combination order, so it is worth considering adding the structural information of the peptide sequence to the deep model training. In this work, we use positional encoding to represent the structural information of peptide sequences and validly combine the positional encoding with existing models by different strategies. Experiments on three datasets show that the introduction of position-coding information can further improve the performance built upon the existing model. The idea of introducing positional encoding to this field can provide important reference significance for the optimization of the deep network structure in the future.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY

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

In Silico Tools for Predicting Novel Epitopes

Identifying antigens within a pathogen is a critical task to develop effective vaccines and diagnostic methods, as well as understanding the evolution and adaptation to host immune responses. Historically, antigenicity was studied with experiments that evaluate the immune response against selected fragments of pathogens. Using this approach, the scientific community has gathered abundant information regarding which pathogenic fragments are immunogenic. The systematic collection of this data has enabled unraveling many of the fundamental rules underlying the properties defining epitopes and immunogenicity, and has resulted in the creation of a large panel of immunologically relevant predictive (in silico) tools. The development and application of such tools have proven to accelerate the identification of novel epitopes within biomedical applications reducing experimental costs. This chapter introduces some basic concepts about MHC presentation, T cell and B cell epitopes, the experimental efforts to determine those, and focuses on state-of-the-art methods for epitope prediction, highlighting their strengths and limitations, and catering instructions for their rational use.

An immunoinformatics study reveals a new BoLA-DR-restricted CD4+ T cell epitopes on the Gag protein of bovine leukemia virus

Bovine leukemia virus (BLV) is the causative agent of enzootic bovine leucosis (EBL), which has been reported worldwide. The expression of viral structural proteins: surface glycoprotein (gp51) and three core proteins - p15 (matrix), p24 (capsid), and p12 (nucleocapsid) induce a strong humoral and cellular immune response at first step of infection. CD4+ T-cell activation is generally induced by bovine leukocyte antigen (BoLA) region- positive antigen-presenting cells (APC) after processing of an exogenous viral antigen. Limited data are available on the BLV epitopes from the core proteins recognized by CD4+ T-cells. Thus, immunoinformatic analysis of Gag sequences obtained from 125 BLV isolates from Poland, Canada, Pakistan, Kazakhstan, Moldova and United States was performed to identify the presence of BoLA-DRB3 restricted CD4+ T-cell epitopes. The 379 15-mer overlapping peptides spanning the entire Gag sequence were run in BoLA-DRB3 allele-binding regions using a BoLA-DRB- peptide binding affinity prediction algorithm. The analysis identified 22 CD4+ T-cell peptide epitopes of variable length ranging from 17 to 22 amino acids. The predicted epitopes interacted with 73 different BoLA-DRB3 alleles found in BLV-infected cattle. Importantly, two epitopes were found to be linked with high proviral load in PBMC. A majority of dominant and subdominant epitopes showed high conservation across different viral strains, and therefore could be attractive targets for vaccine development.

Accurate prediction of HLA class II antigen presentation across all loci using tailored data acquisition and refined machine learning

Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules is crucial for rational development of immunotherapies and vaccines targeting CD4+ T cell activation. So far, most prediction methods for HLA class II antigen presentation have focused on HLA-DR because of limited availability of immunopeptidomics data for HLA-DQ and HLA-DP while not taking into account alternative peptide binding modes. We present an update to the NetMHCIIpan prediction method, which closes the performance gap between all three HLA class II loci. We accomplish this by first integrating large immunopeptidomics datasets describing the HLA class II specificity space across all loci using a refined machine learning framework that accommodates inverted peptide binders. Next, we apply targeted immunopeptidomics assays to generate data that covers additional HLA-DP specificities. The final method, NetMHCIIpan-4.3, achieves high accuracy and molecular coverage across all HLA class II allotypes.

Structural Framework for Analysis of CD4+ T-Cell Epitope Dominance in Viral Fusion Proteins

Antigen conformation shapes CD4+ T-cell specificity through mechanisms of antigen processing, and the consequences for immunity may rival those from conformational effects on antibody specificity. CD4+ T cells initiate and control immunity to pathogens and cancer and are at least partly responsible for immunopathology associated with infection, autoimmunity, and allergy. The primary trigger for CD4+ T-cell maturation is the presentation of an epitope peptide in the MHC class II antigen-presenting protein (MHCII), most commonly on an activated dendritic cell, and then the T-cell responses are recalled by subsequent presentations of the epitope peptide by the same or other antigen-presenting cells. Peptide presentation depends on the proteolytic fragmentation of the antigen in an endosomal/lysosomal compartment and concomitant loading of the fragments into the MHCII, a multistep mechanism called antigen processing and presentation. Although the role of peptide affinity for MHCII has been well studied, the role of proteolytic fragmentation has received less attention. In this Perspective, we will briefly summarize evidence that antigen resistance to unfolding and proteolytic fragmentation shapes the specificity of the CD4+ T-cell response to selected viral envelope proteins, identify several remarkable examples in which the immunodominant CD4+ epitopes most likely depend on the interaction of processing machinery with antigen conformation, and outline how knowledge of antigen conformation can inform future efforts to design vaccines.

Investigating the human and nonobese diabetic mouse MHC class II immunopeptidome using protein language modeling

MOTIVATION

Identifying peptides associated with the major histocompability complex class II (MHCII) is a central task in the evaluation of the immunoregulatory function of therapeutics and drug prototypes. MHCII-peptide presentation prediction has multiple biopharmaceutical applications, including the safety assessment of biologics and engineered derivatives in silico, or the fast progression of antigen-specific immunomodulatory drug discovery programs in immune disease and cancer. This has resulted in the collection of large-scale datasets on adaptive immune receptor antigenic responses and MHC-associated peptide proteomics. In parallel, recent deep learning algorithmic advances in protein language modeling have shown potential in leveraging large collections of sequence data and improve MHC presentation prediction.

RESULTS

Here, we train a compact transformer model (AEGIS) on human and mouse MHCII immunopeptidome data, including a preclinical murine model, and evaluate its performance on the peptide presentation prediction task. We show that the transformer performs on par with existing deep learning algorithms and that combining datasets from multiple organisms increases model performance. We trained variants of the model with and without MHCII information. In both alternatives, the inclusion of peptides presented by the I-Ag7 MHC class II molecule expressed by nonobese diabetic mice enabled for the first time the accurate in silico prediction of presented peptides in a preclinical type 1 diabetes model organism, which has promising therapeutic applications.

AVAILABILITY AND IMPLEMENTATION

The source code is available at https://github.com/Novartis/AEGIS.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

Machine learning reveals limited contribution of trans-only encoded variants to the HLA-DQ immunopeptidome

Human leukocyte antigen (HLA) class II antigen presentation is key for controlling and triggering T cell immune responses. HLA-DQ molecules, which are believed to play a major role in autoimmune diseases, are heterodimers that can be formed as both cis and trans variants depending on whether the α- and β-chains are encoded on the same (cis) or opposite (trans) chromosomes. So far, limited progress has been made for predicting HLA-DQ antigen presentation. In addition, the contribution of trans-only variants (i.e. variants not observed in the population as cis) in shaping the HLA-DQ immunopeptidome remains largely unresolved. Here, we seek to address these issues by integrating state-of-the-art immunoinformatics data mining models with large volumes of high-quality HLA-DQ specific mass spectrometry immunopeptidomics data. The analysis demonstrates highly improved predictive power and molecular coverage for models trained including these novel HLA-DQ data. More importantly, investigating the role of trans-only HLA-DQ variants reveals a limited to no contribution to the overall HLA-DQ immunopeptidome. In conclusion, this study furthers our understanding of HLA-DQ specificities and casts light on the relative role of cis versus trans-only HLA-DQ variants in the HLA class II antigen presentation space. The developed method, NetMHCIIpan-4.2, is available at https://services.healthtech.dtu.dk/services/NetMHCIIpan-4.2 .

Machine learning predictions of MHC-II specificities reveal alternative binding mode of class II epitopes

CD4⁺ T cells orchestrate the adaptive immune response against pathogens and cancer by recognizing epitopes presented on class II major histocompatibility complex (MHC-II) molecules. The high polymorphism of MHC-II genes represents an important hurdle toward accurate prediction and identification of CD4⁺ T cell epitopes. Here we collected and curated a dataset of 627,013 unique MHC-II ligands identified by mass spectrometry. This enabled us to precisely determine the binding motifs of 88 MHC-II alleles across humans, mice, cattle, and chickens. Analysis of these binding specificities combined with X-ray crystallography refined our understanding of the molecular determinants of MHC-II motifs and revealed a widespread reverse-binding mode in HLA-DP ligands. We then developed a machine-learning framework to accurately predict binding specificities and ligands of any MHC-II allele. This tool improves and expands predictions of CD4⁺ T cell epitopes and enables us to discover viral and bacterial epitopes following the aforementioned reverse-binding mode.

Contemplating immunopeptidomes to better predict them

The identification of T-cell epitopes is key for a complete molecular understanding of immune recognition mechanisms in infectious diseases, autoimmunity and cancer. T-cell epitopes further provide targets for personalized vaccines and T-cell therapy, with several therapeutic applications in cancer immunotherapy and elsewhere. T-cell epitopes consist of short peptides displayed on Major Histocompatibility Complex (MHC) molecules. The recent advances in mass spectrometry (MS) based technologies to profile the ensemble of peptides displayed on MHC molecules - the so-called immunopeptidome - had a major impact on our understanding of antigen presentation and MHC ligands. On the one hand, these techniques enabled researchers to directly identify hundreds of thousands of peptides presented on MHC molecules, including some that elicited T-cell recognition. On the other hand, the data collected in these experiments revealed fundamental properties of antigen presentation pathways and significantly improved our ability to predict naturally presented MHC ligands and T-cell epitopes across the wide spectrum of MHC alleles found in human and other organisms. Here we review recent computational developments to analyze experimentally determined immunopeptidomes and harness these data to improve our understanding of antigen presentation and MHC binding specificities, as well as our ability to predict MHC ligands. We further discuss the strengths and limitations of the latest approaches to move beyond predictions of antigen presentation and tackle the challenges of predicting TCR recognition and immunogenicity.

Assessment of commercial polymers with and without reactive groups using amino acid derivative reactivity assay based on both molar concentration approach and gravimetric approach

The amino acid derivative reactivity assay (ADRA), an alternative method for testing skin sensitization, has been established based on the molar concentration approach. However, the additional development of gravimetric concentration and fluorescence detection methods has expanded its range of application to mixtures, which cannot be evaluated using the conventional testing method, the direct peptide reactivity assay (DPRA). Although polymers are generally treated as mixtures, there have been no reports of actual polymer evaluations using alternative methods owing to their insolubility. Therefore, in this study, we evaluated skin sensitization potential of polymers, which is difficult to predict, using ADRA. As polymers have molecular weights ranging from several thousand to more than several tens of thousand Daltons, they are unlikely to cause skin sensitization due to their extremely low penetration into the skin, according to the 500-Da rule. However, if highly reactive functional groups remain at the ends or side chains of polymers, relatively low-molecular-weight polymer components may penetrate the skin to cause sensitization. Polymers can be roughly classified into three major types based on the features of their constituent monomers; we investigated the sensitization capacity of each type of polymer. Polymers with alert sensitization structures at their ends were classified as skin sensitizers, whereas those with no residual reactive groups were classified as nonsensitizers. Although polymers with a glycidyl group need to be evaluated carefully, we concluded that ADRA (0.5 mg/ml) is generally sufficient for polymer hazard assessment.

SphereFace Revived: Unifying Hyperspherical Face Recognition

This paper addresses the deep face recognition problem under an open-set protocol, where ideal face features are expected to have smaller maximal intra-class distance than minimal inter-class distance under a suitably chosen metric space. To this end, hyperspherical face recognition, as a promising line of research, has attracted increasing attention and gradually become a major focus in face recognition research. As one of the earliest works in hyperspherical face recognition, SphereFace explicitly proposed to learn face embeddings with large inter-class angular margin. However, SphereFace still suffers from severe training instability which limits its application in practice. In order to address this problem, we introduce a unified framework to understand large angular margin in hyperspherical face recognition. Under this framework, we extend the study of SphereFace and propose an improved variant with substantially better training stability - SphereFace-R. Specifically, we propose two novel ways to implement the multiplicative margin, and study SphereFace-R under three different feature normalization schemes (no feature normalization, hard feature normalization and soft feature normalization). We also propose an implementation strategy - "characteristic gradient detachment" - to stabilize training. Extensive experiments on SphereFace-R show that it is consistently better than or competitive with state-of-the-art methods.

The MHC Motif Atlas: a database of MHC binding specificities and ligands

The highly polymorphic Major Histocompatibility Complex (MHC) genes are responsible for the binding and cell surface presentation of pathogen or cancer specific T-cell epitopes. This process is fundamental for eliciting T-cell recognition of infected or malignant cells. Epitopes displayed on MHC molecules further provide therapeutic targets for personalized cancer vaccines or adoptive T-cell therapy. To help visualizing, analyzing and comparing the different binding specificities of MHC molecules, we developed the MHC Motif Atlas (http://mhcmotifatlas.org/). This database contains information about thousands of class I and class II MHC molecules, including binding motifs, peptide length distributions, motifs of phosphorylated ligands, multiple specificities or links to X-ray crystallography structures. The database further enables users to download curated datasets of MHC ligands. By combining intuitive visualization of the main binding properties of MHC molecules together with access to more than a million ligands, the MHC Motif Atlas provides a central resource to analyze and interpret the binding specificities of MHC molecules.

The role of antigen expression in shaping the repertoire of HLA presented ligands

Human leukocyte antigen (HLA) presentation of peptides is a prerequisite of T cell immune activation. The understanding of the rules defining this event has large implications for our knowledge of basic immunology and for the rational design of immuno-therapeutics and vaccines. Historically, most of the available prediction methods have been solely focused on the information related to antigen processing and presentation. Recent work has, however, demonstrated that method performance can be boosted by integrating information related to antigen abundance. Here we expand on these later findings and develop an extended version of NetMHCpan, called NetMHCpanExp, integrating information on antigen abundance from RNA-Seq experiments. In line with earlier works, the model demonstrates improved performance for both HLA ligand and cancer neoantigen epitope prediction. Optimal results are obtained by use of sample-specific abundance information but also reference datasets can be applied with a limited performance drop. The developed tool is available at https://services.healthtech.dtu.dk/service.php?NetMHCpanExp-1.0.

CD4+ T-Cell Epitope Prediction by Combined Analysis of Antigen Conformational Flexibility and Peptide-MHCII Binding Affinity

Antigen processing in the class II MHC pathway depends on conventional proteolytic enzymes, potentially acting on antigens in native-like conformational states. CD4+ epitope dominance arises from a competition among antigen folding, proteolysis, and MHCII binding. Protease-sensitive sites, linear antibody epitopes, and CD4+ T-cell epitopes were mapped in plague vaccine candidate F1-V to evaluate the various contributions to CD4+ epitope dominance. Using X-ray crystal structures, antigen processing likelihood (APL) predicts CD4+ epitopes with significant accuracy for F1-V without considering peptide-MHCII binding affinity. We also show that APL achieves excellent performance over two benchmark antigen sets. The profiles of conformational flexibility derived from the X-ray crystal structures of the F1-V proteins, Caf1 and LcrV, were similar to the biochemical profiles of linear antibody epitope reactivity and protease sensitivity, suggesting that the role of structure in proteolysis was captured by the analysis of the crystal structures. The patterns of CD4+ T-cell epitope dominance in C57BL/6, CBA, and BALB/c mice were compared to epitope predictions based on APL, MHCII binding, or both. For a sample of 13 diverse antigens, the accuracy of epitope prediction by the combination of APL and I-A^b-MHCII-peptide affinity reached 36%. When MHCII allele specificity was also diverse, such as in human immunity, prediction of dominant epitopes by APL alone reached 42% when using a stringent scoring threshold. Because dominant CD4+ epitopes tend to occur in conformationally stable antigen domains, crystal structures typically are available for analysis by APL, and thus, the requirement for a crystal structure is not a severe limitation.

A high-throughput yeast display approach to profile pathogen proteomes for MHC-II binding

T cells play a critical role in the adaptive immune response, recognizing peptide antigens presented on the cell surface by major histocompatibility complex (MHC) proteins. While assessing peptides for MHC binding is an important component of probing these interactions, traditional assays for testing peptides of interest for MHC binding are limited in throughput. Here, we present a yeast display-based platform for assessing the binding of tens of thousands of user-defined peptides in a high-throughput manner. We apply this approach to assess a tiled library covering the SARS-CoV-2 proteome and four dengue virus serotypes for binding to human class II MHCs, including HLA-DR401, -DR402, and -DR404. While the peptide datasets show broad agreement with previously described MHC-binding motifs, they additionally reveal experimentally validated computational false positives and false negatives. We therefore present this approach as able to complement current experimental datasets and computational predictions. Further, our yeast display approach underlines design considerations for epitope identification experiments and serves as a framework for examining relationships between viral conservation and MHC binding, which can be used to identify potentially high-interest peptide binders from viral proteins. These results demonstrate the utility of our approach to determine peptide-MHC binding interactions in a manner that can supplement and potentially enhance current algorithm-based approaches.

The interdependence of machine learning and LC-MS approaches for an unbiased understanding of the cellular immunopeptidome

INTRODUCTION

The comprehensive collection of peptides presented by Major Histocompatibility Complex (MHC) molecules on the cell surface is collectively known as the immunopeptidome. The analysis and interpretation of such data sets holds great promise for furthering our understanding of basic immunology and adaptive immune activation and regulation, and for direct rational discovery of T cell antigens and the design of T-cell based therapeutics and vaccines. These applications are however challenged by the complex nature of immunopeptidome data.

AREAS COVERED

Here, we describe the benefits and shortcomings of applying liquid chromatography-tandem mass spectrometry (MS) to obtain large scale immunopeptidome data sets and illustrate how the accurate analysis and optimal interpretation of such data is reliant on the availability of refined and highly optimized machine learning approaches.

EXPERT OPINION

Further we demonstrate how the accuracy of immunoinformatics prediction methods within the field of MHC antigen presentation has benefited greatly from the availability of MS-immunopeptidomics data, and exemplify how optimal antigen discovery is best performed in a synergistic combination of MS experiments and such in silico models trained on large scale immunopeptidomics data.

Deciphering the landscape of phosphorylated HLA-II ligands

CD4⁺ T cell activation in infectious diseases and cancer is governed by the recognition of peptides presented on class II human leukocyte antigen (HLA-II) molecules. Therefore, HLA-II ligands represent promising targets for vaccine design and personalized cancer immunotherapy. Much work has been done to identify and predict unmodified peptides presented on HLA-II molecules. However, little is known about the presentation of phosphorylated HLA-II ligands. Here, we analyzed Mass Spectrometry HLA-II peptidomics data and identified 1,943 unique phosphorylated HLA-II ligands. This enabled us to precisely define phosphorylated binding motifs for more than 30 common HLA-II alleles and to explore various molecular properties of phosphorylated peptides. Our data were further used to develop the first predictor of phosphorylated peptide presentation on HLA-II molecules.

•

1,943 unique phosphorylated HLA-II ligands from MS HLA-II peptidomics data

•

Binding motifs of phosphorylated HLA-II ligands identified for more than 30 alleles

•

Predictor trained on phosphorylated peptides achieves higher accuracy

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

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

Adaptive selection signatures in river buffalo with emphasis on immune and major histocompatibility complex genes

River buffalo is an agriculturally important species with many traits, such as disease tolerance, which promote its use worldwide. Highly contiguous genome assemblies of the river buffalo, goat, pig, human and two cattle subspecies were aligned to study gene gains and losses and signs of positive selection. The gene families that have changed significantly in river buffalo since divergence from cattle play important roles in protein degradation, the olfactory receptor system, detoxification and the immune system. We used the branch site model in PAML to analyse single-copy orthologs to identify positively selected genes that may be involved in skin differentiation, mammary development and bone formation in the river buffalo branch. The high contiguity of the genomes enabled evaluation of differences among species in the major histocompatibility complex. We identified a Babesia-like L1 LINE insertion in the DRB1-like gene in the river buffalo and discuss the implication of this finding.

Targeted immune epitope prediction to HHLA2 and MAGEB5 protein variants as therapeutic approach to related viral diseases

Background

Targeted immunotherapy is mostly associated with cancer treatment wherein designed molecules engage signaling pathways and mutant proteins critical to the survival of the cell. One of several genetic approaches is the use of in silico methods to develop immune epitopes targeting specific antigenic regions on related mutant proteins. In a recent study we showed a functional association between the gamma retrovirus HERV-H Long Terminal Associating (HHLA1, HHLA2 and HHLA3) proteins and melanoma associated antigen of the B class proteins (MAGEB5), with a resultant decrease in expression of HLA class I and II immune variants. HLA-C and HLA-DRB5 were the main HLA class I and II Immune variants, respectively, that showed expression changes across viral samples of interest. Specific immune variants for HLA-C and HLA-DRB5 were filtered for the top ten based on their relative frequency of counts across the samples.

Results

Protein variants for HHLA1, HHLA2, HHLA3 and MAGEB5 were used to predict antigenic epitope peptides to immune peptide-MHC class I and II binding using artificial neural networks. For IC50 peptide scores (PS) ≥ 0.5 with a transformed binding ability between 0 and 1, the top 5 epitopes identified for all targeted genes HHLA1,2 & 3 and MAGEB5 were qualified as strong or weak binders according to the threshold. Domain analysis using NCBI Conserved Domain Database (CDD) identified HHLA2 with immunoglobulin-like domains (Ig_C1-set) and MAGEB5 with the MAGE Homology Domain (MHD). Linear regression showed a statistical correlation (P < 0.001) for HHLA2 and MAGEB5 predicted epitope peptides to HLA-C but not HLA-DRB5. The prediction model identified HLA-C variant 9 (HLA-C9, BAA08825.1 HLA-B*1511) at 1.1% as the most valuable immune target for clinical considerations. Identification of the 9-mer epitope peptide within the domain showed for HHLA2: YANRTSLFY (PS = 0.5837) and VLAYYLSSSQNTIIN (PS = 0.77) for HLA-C and HLA-DRB5, respectively and for MAGEB5, peptides: FVRLTYLEY (PS = 0.5293) and YPAHYQFLWGPRAYT (PS = 0.62) for HLA-C and HLA-DRB5, respectively.

Conclusion

Specific immune responses to targeted epitope peptides and their prediction models, suggested co-expression and co-evolution for HHLA2 and MAGEB5 in viral related diseases. HHLA2 and MAGEB5 could be considered markers for virus related tumors and targeted therapy for oncogenic diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12865-021-00440-w.

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action

BACKGROUND

HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR).

METHODS

For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used.

RESULTS

All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found.

CONCLUSIONS

Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

MS-Based HLA-II Peptidomics Combined With Multiomics Will Aid the Development of Future Immunotherapies

Immunotherapies have emerged to treat diseases by selectively modulating a patient's immune response. Although the roles of T and B cells in adaptive immunity have been well studied, it remains difficult to select targets for immunotherapeutic strategies. Because human leukocyte antigen class II (HLA-II) peptides activate CD4+ T cells and regulate B cell activation, proliferation, and differentiation, these peptide antigens represent a class of potential immunotherapy targets and biomarkers. To better understand the molecular basis of how HLA-II antigen presentation is involved in disease progression and treatment, systematic HLA-II peptidomics combined with multiomic analyses of diverse cell types in healthy and diseased states is required. For this reason, MS-based innovations that facilitate investigations into the interplay between disease pathologies and the presentation of HLA-II peptides to CD4+ T cells will aid in the development of patient-focused immunotherapies.

Challenges targeting cancer neoantigens in 2021: a systematic literature review

Introduction: Cancer neoantigens represent important targets of cancer immunotherapy. The goal of cancer neoantigen vaccines is to induce neoantigen-specific immune responses and antitumor immunity while minimizing the potential for autoimmune toxicity. Advances in sequencing technologies, neoantigen prediction algorithms, and other technologies have dramatically improved the ability to identify and prioritize cancer neoantigens. Unfortunately, results from preclinical studies and early phase clinical trials highlight important challenges to the successful clinical translation of neoantigen cancer vaccines.Areas covered: In this review, we provide an overview of current strategies for the identification and prioritization of cancer neoantigens with a particular emphasis on the two most common strategies used for neoantigen identification: (1) direct identification of peptide ligands eluted from peptide-MHC complexes, and (2) next-generation sequencing combined with neoantigen prediction algorithms. We highlight the limitations of current neoantigen prediction pipelines, and discuss broader challenges associated with cancer neoantigen vaccines including tumor purity/heterogeneity and the immunosuppressive tumor microenvironment.Expert opinion: Despite current limitations, neoantigen prediction is likely to improve rapidly based on advances in sequencing, machine learning, and information sharing. The successful development of robust cancer neoantigen prediction strategies is likely to have a significant impact, with the potential to facilitate cancer neoantigen vaccine design.

HLA Class II Presentation Is Specifically Altered at Elevated Temperatures in the B-Lymphoblastic Cell Line JY

Human leukocyte antigen (HLA) molecules play critical roles in our adaptive immune system by signaling a cell's health status to the immune system, through presentation of small peptides. Understanding HLA biology is important because of its prominent role in autoimmune diseases and cancer immunotherapy. Although both the HLA class I and class II antigen processing and presentation pathways have been studied extensively, the fundamental rules in HLA class II antigen presentation still remain less understood. To clarify the mechanistic and adaptive differences between the HLA systems, we challenged a B lymphoblastic cell line (JY), widely used as model system in studying antigen presentation, with a high temperature treatment to mimic a "fever-like state", representing one of the most common physiological responses to infection. In the absence of real invading pathogenic peptides to present, we could focus on delineating the intrinsic HLA pathway adaptations in response to high temperature in this particular cell line. Following a three-pronged approach, we performed quantitative analyses of the proteome, the HLA class I ligandome, as well as the HLA class II ligandome. The data reveals that elevated temperature may already prepare these cells for an immune-like response through increased HLA class II presentation capacity and specific release of, from the invariant chain originating, CLIP peptides. Interestingly, at high temperature, prominent changes in the composition of the CLIP repertoire were observed, with enrichment of peptides containing C-terminal extensions beyond the CLIP-core region. Collectively, these illustrate intriguing temperature sensitive adaptations in this B cell line.

Applications of Machine and Deep Learning in Adaptive Immunity

Adaptive immunity is mediated by lymphocyte B and T cells, which respectively express a vast and diverse repertoire of B cell and T cell receptors and, in conjunction with peptide antigen presentation through major histocompatibility complexes (MHCs), can recognize and respond to pathogens and diseased cells. In recent years, advances in deep sequencing have led to a massive increase in the amount of adaptive immune receptor repertoire data; additionally, proteomics techniques have led to a wealth of data on peptide-MHC presentation. These large-scale data sets are now making it possible to train machine and deep learning models, which can be used to identify complex and high-dimensional patterns in immune repertoires. This article introduces adaptive immune repertoires and machine and deep learning related to biological sequence data and then summarizes the many applications in this field, which span from predicting the immunological status of a host to the antigen specificity of individual receptors and the engineering of immunotherapeutics.

Integral Use of Immunopeptidomics and Immunoinformatics for the Characterization of Antigen Presentation and Rational Identification of BoLA-DR-Presented Peptides and Epitopes

MHC peptide binding and presentation is the most selective event defining the landscape of T cell epitopes. Consequently, understanding the diversity of MHC alleles in a given population and the parameters that define the set of ligands that can be bound and presented by each of these alleles (the immunopeptidome) has an enormous impact on our capacity to predict and manipulate the potential of protein Ags to elicit functional T cell responses. Liquid chromatography-mass spectrometry analysis of MHC-eluted ligand data has proven to be a powerful technique for identifying such peptidomes, and methods integrating such data for prediction of Ag presentation have reached a high level of accuracy for both MHC class I and class II. In this study, we demonstrate how these techniques and prediction methods can be readily extended to the bovine leukocyte Ag class II DR locus (BoLA-DR). BoLA-DR binding motifs were characterized by eluted ligand data derived from bovine cell lines expressing a range of DRB3 alleles prevalent in Holstein-Friesian populations. The model generated (NetBoLAIIpan, available as a Web server at www.cbs.dtu.dk/services/NetBoLAIIpan) was shown to have unprecedented predictive power to identify known BoLA-DR-restricted CD4 epitopes. In summary, the results demonstrate the power of an integrated approach combining advanced mass spectrometry peptidomics with immunoinformatics for characterization of the BoLA-DR Ag presentation system and provide a prediction tool that can be used to assist in rational evaluation and selection of bovine CD4 T cell epitopes.

A Simplified Amino Acidic Alphabet to Unveil the T-Cells Receptors Antigens: A Computational Perspective

The exposure to pathogens triggers the activation of adaptive immune responses through antigens bound to surface receptors of antigen presenting cells (APCs). T cell receptors (TCR) are responsible for initiating the immune response through their physical direct interaction with antigen-bound receptors on the APCs surface. The study of T cell interactions with antigens is considered of crucial importance for the comprehension of the role of immune responses in cancer growth and for the subsequent design of immunomodulating anticancer drugs. RNA sequencing experiments performed on T cells represented a major breakthrough for this branch of experimental molecular biology. Apart from the gene expression levels, the hypervariable CDR3α/β sequences of the TCR loops can now be easily determined and modelled in the three dimensions, being the portions of TCR mainly responsible for the interaction with APC receptors. The most direct experimental method for the investigation of antigens would be based on peptide libraries, but their huge combinatorial nature, size, cost, and the difficulty of experimental fine tuning makes this approach complicated time consuming, and costly. We have implemented in silico methodology with the aim of moving from CDR3α/β sequences to a library of potentially antigenic peptides that can be used in immunologically oriented experiments to study T cells’ reactivity. To reduce the size of the library, we have verified the reproducibility of experimental benchmarks using the permutation of only six residues that can be considered representative of all ensembles of 20 natural amino acids. Such a simplified alphabet is able to correctly find the poses and chemical nature of original antigens within a small subset of ligands of potential interest. The newly generated library would have the advantage of leading to potentially antigenic ligands that would contribute to a better understanding of the chemical nature of TCR-antigen interactions. This step is crucial in the design of immunomodulators targeted towards T-cells response as well as in understanding the first principles of an immune response in several diseases, from cancer to autoimmune disorders.

Comprehensive analysis of T cell immunodominance and immunoprevalence of SARS-CoV-2 epitopes in COVID-19 cases

T cells are involved in control of SARS-CoV-2 infection. To establish the patterns of immunodominance of different SARS-CoV-2 antigens and precisely measure virus-specific CD4⁺ and CD8⁺ T cells, we study epitope-specific T cell responses of 99 convalescent coronavirus disease 2019 (COVID-19) cases. The SARS-CoV-2 proteome is probed using 1,925 peptides spanning the entire genome, ensuring an unbiased coverage of human leukocyte antigen (HLA) alleles for class II responses. For HLA class I, we study an additional 5,600 predicted binding epitopes for 28 prominent HLA class I alleles, accounting for wide global coverage. We identify several hundred HLA-restricted SARS-CoV-2-derived epitopes. Distinct patterns of immunodominance are observed, which differ for CD4⁺ T cells, CD8⁺ T cells, and antibodies. The class I and class II epitopes are combined into epitope megapools to facilitate identification and quantification of SARS-CoV-2-specific CD4⁺ and CD8⁺ T cells.

Competition-Based Cell Assay Employing Soluble T Cell Receptors to Assess MHC Class II Antigen Processing and Presentation

Accurate assessment of antigen-specific immune responses is critical in the development of safe and efficacious biotherapeutics and vaccines. Endosomal processing of a protein antigen followed by presentation on major histocompatibility complex (MHC) class II constitute necessary steps in the induction of CD4+ T cell immune responses. Current preclinical methods for assessing immunogenicity risk consist of in vitro cell-based assays and computational prediction tools. Cell-based assays are time and labor-intensive while in silico methodologies have limitations. Here, we propose a novel cell-based assay capable of investigating an antigen's endosomal processing and MHC class II presentation capabilities. This novel assay relies on competition between epitopes for MHC class II binding and employs labeled soluble T cell receptors (sTCRs) as detectors of epitope presentation.

Tumor microenvironment characterization in stage IV gastric cancer

Immunotherapy is remarkably affected by the immune environment of the principal tumor. Nonetheless, the immune environment's clinical relevance in stage IV gastric cancer (GC) is largely unknown. The gene expression profiles of 403 stage IV GC patients in the three cohorts: GEO (Gene Expression Omnibus, GSE84437 (n=292) and GSE62254 (n=77), and TCGA (The Cancer Genome Atlas, n=34) were used in the present study. Using four publicly available stage IV GC expression datasets, 29 immune signatures were expression profiled, and on this basis, we classified stage IV GC. The classification was conducted using the hierarchical clustering method. Three stage IV GC subtypes L, M, and H were identified representing low, medium, and high immunity, respectively. Immune correlation analysis of these three types revealed that Immune H exhibited a better prognostic outcome as well as a higher immune score compared with other subtypes. There was a noticeable difference in the three subgroups of HLA genes. Further, on comparing with other subtypes, CD86, CD80, CD274, CTLA4, PDCD1, and PDCD1LG2 had higher expression in the Immunity H subtype. In stage IV GC, potentially positive associations between immune and pathway activities were displayed, due to the enrichment of pathways including TNF signaling, Th-17 cell differentiation, and JAK-STAT signaling pathways in Immunity H vs Immunity L subtypes. External cohorts from TCGA cohort ratified these results. The identification of stage IV GC subtypes has potential clinical implications in stage IV GC treatment.

Comprehensive analysis of T cell immunodominance and immunoprevalence of SARS-CoV-2 epitopes in COVID-19 cases

T cells are involved in control of SARS-CoV-2 infection. To establish the patterns of immunodominance of different SARS-CoV-2 antigens, and precisely measure virus-specific CD4 ⁺ and CD8 ⁺ T cells, we studied epitope-specific T cell responses of approximately 100 convalescent COVID-19 cases. The SARS-CoV-2 proteome was probed using 1,925 peptides spanning the entire genome, ensuring an unbiased coverage of HLA alleles for class II responses. For HLA class I, we studied an additional 5,600 predicted binding epitopes for 28 prominent HLA class I alleles, accounting for wide global coverage. We identified several hundred HLA-restricted SARS-CoV-2-derived epitopes. Distinct patterns of immunodominance were observed, which differed for CD4 ⁺ T cells, CD8 ⁺ T cells, and antibodies. The class I and class II epitopes were combined into new epitope megapools to facilitate identification and quantification of SARS-CoV-2-specific CD4 ⁺ and CD8 ⁺ T cells.

Comparison of HLA ligand elution data and binding predictions reveals varying prediction performance for the multiple motifs recognized by HLA-DQ2.5

Binding prediction tools are commonly used to identify peptides presented on MHC class II molecules. Recently, a wealth of data in the form of naturally eluted ligands has become available and discrepancies between ligand elution data and binding predictions have been reported. Quantitative metrics for such comparisons are currently lacking. In this study, we assessed how efficiently MHC class II binding predictions can identify naturally eluted peptides, and investigated instances with discrepancies between the two methods in detail. We found that, in general, MHC class II eluted ligands are predicted to bind to their reported restriction element with high affinity. But, for several studies reporting an increased number of ligands that were not predicted to bind, we found that the reported MHC restriction was ambiguous. Additional analyses determined that most of the ligands predicted to not bind, are predicted to bind other co‐expressed MHC class II molecules. For selected alleles, we addressed discrepancies between elution data and binding predictions by experimental measurements and found that predicted and measured affinities correlate well. For DQA1*05:01/DQB1*02:01 (DQ2.5) however, binding predictions did miss several peptides that were determined experimentally to be binders. For these peptides and several known DQ2.5 binders, we determined key residues for conferring DQ2.5 binding capacity, which revealed that DQ2.5 utilizes two different binding motifs, of which only one is predicted effectively. These findings have important implications for the interpretation of ligand elution data and for the improvement of MHC class II binding predictions.

Improved prediction of HLA antigen presentation hotspots: Applications for immunogenicity risk assessment of therapeutic proteins

Immunogenicity risk assessment is a critical element in protein drug development. Currently, the risk assessment is most often performed using MHC-associated peptide proteomics (MAPPs) and/or T-cell activation assays. However, this is a highly costly procedure that encompasses limited sensitivity imposed by sample sizes, the MHC repertoire of the tested donor cohort and the experimental procedures applied. Recent work has suggested that these techniques could be complemented by accurate, high-throughput and cost-effective prediction of in silico models. However, this work covered a very limited set of therapeutic proteins and eluted ligand (EL) data. Here, we resolved these limitations by showcasing, in a broader setting, the versatility of in silico models for assessment of protein drug immunogenicity. A method for prediction of MHC class II antigen presentation was developed on the hereto largest available mass spectrometry (MS) HLA-DR EL data set. Using independent test sets, the performance of the method for prediction of HLA-DR antigen presentation hotspots was benchmarked. In particular, the method was showcased on a set of protein sequences including four therapeutic proteins and demonstrated to accurately predict the experimental MS hotspot regions at a significantly lower false-positive rate compared with other methods. This gain in performance was particularly pronounced when compared to the NetMHCIIpan-3.2 method trained on binding affinity data. These results suggest that in silico methods trained on MS HLA EL data can effectively and accurately be used to complement MAPPs assays for the risk assessment of protein drugs.

Repertoire-scale determination of class II MHC peptide binding via yeast display improves antigen prediction

CD4⁺ helper T cells contribute important functions to the immune response during pathogen infection and tumor formation by recognizing antigenic peptides presented by class II major histocompatibility complexes (MHC-II). While many computational algorithms for predicting peptide binding to MHC-II proteins have been reported, their performance varies greatly. Here we present a yeast-display-based platform that allows the identification of over an order of magnitude more unique MHC-II binders than comparable approaches. These peptides contain previously identified motifs, but also reveal new motifs that are validated by in vitro binding assays. Training of prediction algorithms with yeast-display library data improves the prediction of peptide-binding affinity and the identification of pathogen-associated and tumor-associated peptides. In summary, our yeast-display-based platform yields high-quality MHC-II-binding peptide datasets that can be used to improve the accuracy of MHC-II binding prediction algorithms, and potentially enhance our understanding of CD4⁺ T cell recognition.

Identifying peptides that can bind major histocompatibility complex II (MHC-II) is important for our understanding of T cell immunity and specificity. Here the authors present a yeast-display library screening approach that identifies more potential binders than various reported algorithms to help expand our understanding for antigen presentation.

A Systematic, Unbiased Mapping of CD8+ and CD4+ T Cell Epitopes in Yellow Fever Vaccinees

Examining CD8+ and CD4+ T cell responses after primary Yellow Fever vaccination in a cohort of 210 volunteers, we have identified and tetramer-validated 92 CD8+ and 50 CD4+ T cell epitopes, many inducing strong and prevalent (i.e., immunodominant) T cell responses. Restricted by 40 and 14 HLA-class I and II allotypes, respectively, these responses have wide population coverage and might be of considerable academic, diagnostic and therapeutic interest. The broad coverage of epitopes and HLA overcame the otherwise confounding effects of HLA diversity and non-HLA background providing the first evidence of T cell immunodomination in humans. Also, double-staining of CD4+ T cells with tetramers representing the same HLA-binding core, albeit with different flanking regions, demonstrated an extensive diversification of the specificities of many CD4+ T cell responses. We suggest that this could reduce the risk of pathogen escape, and that multi-tetramer staining is required to reveal the true magnitude and diversity of CD4+ T cell responses. Our T cell epitope discovery approach uses a combination of (1) overlapping peptides representing the entire Yellow Fever virus proteome to search for peptides containing CD4+ and/or CD8+ T cell epitopes, (2) predictors of peptide-HLA binding to suggest epitopes and their restricting HLA allotypes, (3) generation of peptide-HLA tetramers to identify T cell epitopes, and (4) analysis of ex vivo T cell responses to validate the same. This approach is systematic, exhaustive, and can be done in any individual of any HLA haplotype. It is all-inclusive in the sense that it includes all protein antigens and peptide epitopes, and encompasses both CD4+ and CD8+ T cell epitopes. It is efficient and, importantly, reduces the false discovery rate. The unbiased nature of the T cell epitope discovery approach presented here should support the refinement of future peptide-HLA class I and II predictors and tetramer technologies, which eventually should cover all HLA class I and II isotypes. We believe that future investigations of emerging pathogens (e.g., SARS-CoV-2) should include population-wide T cell epitope discovery using blood samples from patients, convalescents and/or long-term survivors, who might all hold important information on T cell epitopes and responses.

Enabling Routine MHC-II-Associated Peptide Proteomics for Risk Assessment of Drug-Induced Immunogenicity

Major histocompatibility complex-II (MHC-II)-Associated Peptide Proteomics (MAPPs) is a mass spectrometry-based approach to identify and relatively quantitate naturally processed and presented MHC-II-associated peptides that can potentially activate T cells and contribute to the immunogenicity of a drug. Acceptance of the MAPPs technology as an appropriate preclinical (and potentially clinical) immunogenicity risk assessment tool depends not only on its technical stability and robustness but also on the ability to compare results across experiments and donors. To this end, we developed a specialized MAPPs data processing pipeline, dataMAPPs, which presents complex mass spectrometric data sets in the form of heat maps (heatMAPPs), enabling rapid and convenient comparison between conditions and donors. A customized normalization procedure based on identified endogenous peptides standardizes signal intensities within and between donors and enables cross-experimental comparison. We evaluated the technical reproducibility of the MAPPs platform using tool compounds with respect to the most prominent experimental factors and found that the systematic biological differences across donors by far outweighed any technical source of variation. We illustrate the capability of the MAPPs platform to generate data that may be used for preclinical risk assessment of drug-induced immunogenicity and discuss its applicability in the clinics.

Unbiased Characterization of Peptide-HLA Class II Interactions Based on Large-Scale Peptide Microarrays; Assessment of the Impact on HLA Class II Ligand and Epitope Prediction

Human Leukocyte Antigen class II (HLA-II) molecules present peptides to T lymphocytes and play an important role in adaptive immune responses. Characterizing the binding specificity of single HLA-II molecules has profound impacts for understanding cellular immunity, identifying the cause of autoimmune diseases, for immunotherapeutics, and vaccine development. Here, novel high-density peptide microarray technology combined with machine learning techniques were used to address this task at an unprecedented level of high-throughput. Microarrays with over 200,000 defined peptides were assayed with four exemplary HLA-II molecules. Machine learning was applied to mine the signals. The comparison of identified binding motifs, and power for predicting eluted ligands and CD4+ epitope datasets to that obtained using NetMHCIIpan-3.2, confirmed a high quality of the chip readout. These results suggest that the proposed microarray technology offers a novel and unique platform for large-scale unbiased interrogation of peptide binding preferences of HLA-II molecules.

Vaccine candidate designed against carcinoembryonic antigen-related cell adhesion molecules using immunoinformatics tools

Carcinoembryonic antigen-related cell adhesion (CEACAM) molecules belong to a family of membrane glycoproteins that mediate intercellular interactions influencing cellular growth, immune cell activation, apoptosis, and tumor suppression. Several family members (CEACAM1, CEACAM5, and CEACAM6) are highly expressed in cancers, and they share a conserved N-terminal domain that serves as an attractive target for cancer immunotherapy. A multi-epitope vaccine candidate against this conserved domain has been developed using immunoinformatics tools. Specifically, several epitopes predicted to interact with MHC class I and II molecules were linked together with appropriate linkers. The tertiary structure of the vaccine is generated by homology and ab initio modeling. Molecular docking of epitopes to MHC structures have revealed that the lowest energy conformations are the epitopes bound to the antigen-binding groove of the MHC molecules. Subsequent molecular dynamics simulation has confirmed the stability of the binding conformations in solution. The predicted vaccine has relatively high antigenicity and low allergenicity, suggesting that it is an ideal candidate for further refinement and development.Communicated by Ramaswamy H. Sarma.

NetMHCpan-4.1 and NetMHCIIpan-4.0: improved predictions of MHC antigen presentation by concurrent motif deconvolution and integration of MS MHC eluted ligand data

Major histocompatibility complex (MHC) molecules are expressed on the cell surface, where they present peptides to T cells, which gives them a key role in the development of T-cell immune responses. MHC molecules come in two main variants: MHC Class I (MHC-I) and MHC Class II (MHC-II). MHC-I predominantly present peptides derived from intracellular proteins, whereas MHC-II predominantly presents peptides from extracellular proteins. In both cases, the binding between MHC and antigenic peptides is the most selective step in the antigen presentation pathway. Therefore, the prediction of peptide binding to MHC is a powerful utility to predict the possible specificity of a T-cell immune response. Commonly MHC binding prediction tools are trained on binding affinity or mass spectrometry-eluted ligands. Recent studies have however demonstrated how the integration of both data types can boost predictive performances. Inspired by this, we here present NetMHCpan-4.1 and NetMHCIIpan-4.0, two web servers created to predict binding between peptides and MHC-I and MHC-II, respectively. Both methods exploit tailored machine learning strategies to integrate different training data types, resulting in state-of-the-art performance and outperforming their competitors. The servers are available at http://www.cbs.dtu.dk/services/NetMHCpan-4.1/ and http://www.cbs.dtu.dk/services/NetMHCIIpan-4.0/.

Immunoinformatics: Predicting Peptide-MHC Binding

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

Immunopeptidomic Data Integration to Artificial Neural Networks Enhances Protein-Drug Immunogenicity Prediction

Recombinant DNA technology has, in the last decades, contributed to a vast expansion of the use of protein drugs as pharmaceutical agents. However, such biological drugs can lead to the formation of anti-drug antibodies (ADAs) that may result in adverse effects, including allergic reactions and compromised therapeutic efficacy. Production of ADAs is most often associated with activation of CD4 T cell responses resulting from proteolysis of the biotherapeutic and loading of drug-specific peptides into major histocompatibility complex (MHC) class II on professional antigen-presenting cells. Recently, readouts from MHC-associated peptide proteomics (MAPPs) assays have been shown to correlate with the presence of CD4 T cell epitopes. However, the limited sensitivity of MAPPs challenges its use as an immunogenicity biomarker. In this work, MAPPs data was used to construct an artificial neural network (ANN) model for MHC class II antigen presentation. Using Infliximab and Rituximab as showcase stories, the model demonstrated an unprecedented performance for predicting MAPPs and CD4 T cell epitopes in the context of protein-drug immunogenicity, complementing results from MAPPs assays and outperforming conventional prediction models trained on binding affinity data.

Improved Prediction of MHC II Antigen Presentation through Integration and Motif Deconvolution of Mass Spectrometry MHC Eluted Ligand Data

Major histocompatibility complex II (MHC II) molecules play a vital role in the onset and control of cellular immunity. In a highly selective process, MHC II presents peptides derived from exogenous antigens on the surface of antigen-presenting cells for T cell scrutiny. Understanding the rules defining this presentation holds critical insights into the regulation and potential manipulation of the cellular immune system. Here, we apply the NNAlign_MA machine learning framework to analyze and integrate large-scale eluted MHC II ligand mass spectrometry (MS) data sets to advance prediction of CD4+ epitopes. NNAlign_MA allows integration of mixed data types, handling ligands with multiple potential allele annotations, encoding of ligand context, leveraging information between data sets, and has pan-specific power allowing accurate predictions outside the set of molecules included in the training data. Applying this framework, we identified accurate binding motifs of more than 50 MHC class II molecules described by MS data, particularly expanding coverage for DP and DQ beyond that obtained using current MS motif deconvolution techniques. Furthermore, in large-scale benchmarking, the final model termed NetMHCIIpan-4.0 demonstrated improved performance beyond current state-of-the-art predictors for ligand and CD4+ T cell epitope prediction. These results suggest that NNAlign_MA and NetMHCIIpan-4.0 are powerful tools for analysis of immunopeptidome MS data, prediction of T cell epitopes, and development of personalized immunotherapies.

T Cell Epitope Predictions

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