Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not ...trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.
The outbreak of SARS-CoV-2 (2019-nCoV) virus has highlighted the need for fast and efficacious vaccine development. Stimulation of a proper immune response that leads to protection is highly ...dependent on presentation of epitopes to circulating T-cells via the HLA complex. SARS-CoV-2 is a large RNA virus and testing of all of its overlapping peptides in vitro to deconvolute an immune response is not feasible. Therefore HLA-binding prediction tools are often used to narrow down the number of peptides to test. We tested NetMHC suite tools' predictions by using an in vitro peptide-MHC stability assay. We assessed 777 peptides that were predicted to be good binders across 11 MHC alleles in a complex-stability assay and tested a selection of 19 epitope-HLA-binding prediction tools against the assay. In this investigation of potential SARS-CoV-2 epitopes we found that current prediction tools vary in performance when assessing binding stability, and they are highly dependent on the MHC allele in question. Designing a COVID-19 vaccine where only a few epitope targets are included is therefore a very challenging task. Here, we present 174 SARS-CoV-2 epitopes with high prediction binding scores, validated to bind stably to 11 HLA alleles. Our findings may contribute to the design of an efficacious vaccine against COVID-19.
Mutated isocitrate dehydrogenase 1 (IDH1) defines a molecularly distinct subtype of diffuse glioma
. The most common IDH1 mutation in gliomas affects codon 132 and encodes IDH1(R132H), which harbours ...a shared clonal neoepitope that is presented on major histocompatibility complex (MHC) class II
. An IDH1(R132H)-specific peptide vaccine (IDH1-vac) induces specific therapeutic T helper cell responses that are effective against IDH1(R132H)
tumours in syngeneic MHC-humanized mice
. Here we describe a multicentre, single-arm, open-label, first-in-humans phase I trial that we carried out in 33 patients with newly diagnosed World Health Organization grade 3 and 4 IDH1(R132H)
astrocytomas (Neurooncology Working Group of the German Cancer Society trial 16 (NOA16), ClinicalTrials.gov identifier NCT02454634). The trial met its primary safety endpoint, with vaccine-related adverse events restricted to grade 1. Vaccine-induced immune responses were observed in 93.3% of patients across multiple MHC alleles. Three-year progression-free and death-free rates were 0.63 and 0.84, respectively. Patients with immune responses showed a two-year progression-free rate of 0.82. Two patients without an immune response showed tumour progression within two years of first diagnosis. A mutation-specificity score that incorporates the duration and level of vaccine-induced IDH1(R132H)-specific T cell responses was associated with intratumoral presentation of the IDH1(R132H) neoantigen in pre-treatment tumour tissue. There was a high frequency of pseudoprogression, which indicates intratumoral inflammatory reactions. Pseudoprogression was associated with increased vaccine-induced peripheral T cell responses. Combined single-cell RNA and T cell receptor sequencing showed that tumour-infiltrating CD40LG
and CXCL13
T helper cell clusters in a patient with pseudoprogression were dominated by a single IDH1(R132H)-reactive T cell receptor.
Tumor-associated epitopes presented on MHC-I that can activate the immune system against cancer cells are typically identified from annotated protein-coding regions of the genome, but whether ...peptides originating from novel or unannotated open reading frames (nuORFs) can contribute to antitumor immune responses remains unclear. Here we show that peptides originating from nuORFs detected by ribosome profiling of malignant and healthy samples can be displayed on MHC-I of cancer cells, acting as additional sources of cancer antigens. We constructed a high-confidence database of translated nuORFs across tissues (nuORFdb) and used it to detect 3,555 translated nuORFs from MHC-I immunopeptidome mass spectrometry analysis, including peptides that result from somatic mutations in nuORFs of cancer samples as well as tumor-specific nuORFs translated in melanoma, chronic lymphocytic leukemia and glioblastoma. NuORFs are an unexplored pool of MHC-I-presented, tumor-specific peptides with potential as immunotherapy targets.
Binding of peptides to Major Histocompatibility Complex (MHC) molecules is the single most selective step in the recognition of pathogens by the cellular immune system. The human MHC class I system ...(HLA-I) is extremely polymorphic. The number of registered HLA-I molecules has now surpassed 1500. Characterizing the specificity of each separately would be a major undertaking.
Here, we have drawn on a large database of known peptide-HLA-I interactions to develop a bioinformatics method, which takes both peptide and HLA sequence information into account, and generates quantitative predictions of the affinity of any peptide-HLA-I interaction. Prospective experimental validation of peptides predicted to bind to previously untested HLA-I molecules, cross-validation, and retrospective prediction of known HIV immune epitopes and endogenous presented peptides, all successfully validate this method. We further demonstrate that the method can be applied to perform a clustering analysis of MHC specificities and suggest using this clustering to select particularly informative novel MHC molecules for future biochemical and functional analysis.
Encompassing all HLA molecules, this high-throughput computational method lends itself to epitope searches that are not only genome- and pathogen-wide, but also HLA-wide. Thus, it offers a truly global analysis of immune responses supporting rational development of vaccines and immunotherapy. It also promises to provide new basic insights into HLA structure-function relationships. The method is available at http://www.cbs.dtu.dk/services/NetMHCpan.
CD4 positive T helper cells control many aspects of specific immunity. These cells are specific for peptides derived from protein antigens and presented by molecules of the extremely polymorphic ...major histocompatibility complex (MHC) class II system. The identification of peptides that bind to MHC class II molecules is therefore of pivotal importance for rational discovery of immune epitopes. HLA-DR is a prominent example of a human MHC class II. Here, we present a method, NetMHCIIpan, that allows for pan-specific predictions of peptide binding to any HLA-DR molecule of known sequence. The method is derived from a large compilation of quantitative HLA-DR binding events covering 14 of the more than 500 known HLA-DR alleles. Taking both peptide and HLA sequence information into account, the method can generalize and predict peptide binding also for HLA-DR molecules where experimental data is absent. Validation of the method includes identification of endogenously derived HLA class II ligands, cross-validation, leave-one-molecule-out, and binding motif identification for hitherto uncharacterized HLA-DR molecules. The validation shows that the method can successfully predict binding for HLA-DR molecules-even in the absence of specific data for the particular molecule in question. Moreover, when compared to TEPITOPE, currently the only other publicly available prediction method aiming at providing broad HLA-DR allelic coverage, NetMHCIIpan performs equivalently for alleles included in the training of TEPITOPE while outperforming TEPITOPE on novel alleles. We propose that the method can be used to identify those hitherto uncharacterized alleles, which should be addressed experimentally in future updates of the method to cover the polymorphism of HLA-DR most efficiently. We thus conclude that the presented method meets the challenge of keeping up with the MHC polymorphism discovery rate and that it can be used to sample the MHC "space," enabling a highly efficient iterative process for improving MHC class II binding predictions.
Induction of a lasting protective immune response is dependent on presentation of epitopes to patrolling T cells through the HLA complex. While peptide:HLA (pHLA) complex affinity alone is widely ...exploited for epitope selection, we demonstrate that including the pHLA complex stability as a selection parameter can significantly reduce the high false discovery rate observed with predicted affinity. In this study, pHLA complex stability was measured on three common class I alleles and 1286 overlapping 9-mer peptides derived from the SARS-CoV-2 Spike protein. Peptides were pooled based on measured stability and predicted affinity. Strikingly, stability of the pHLA complex was shown to strongly select for immunogenic epitopes able to activate functional CD8
T cells. This result was observed across the three studied alleles and in both vaccinated and convalescent COVID-19 donors. Deconvolution of peptide pools showed that specific CD8
T cells recognized one or two dominant epitopes. Moreover, SARS-CoV-2 specific CD8
T cells were detected by tetramer-staining across multiple donors. In conclusion, we show that stability analysis of pHLA is a key factor for identifying immunogenic epitopes.
Binding of peptides to Major Histocompatibility class II (MHC-II) molecules play a central role in governing responses of the adaptive immune system. MHC-II molecules sample peptides from the ...extracellular space allowing the immune system to detect the presence of foreign microbes from this compartment. Predicting which peptides bind to an MHC-II molecule is therefore of pivotal importance for understanding the immune response and its effect on host-pathogen interactions. The experimental cost associated with characterizing the binding motif of an MHC-II molecule is significant and large efforts have therefore been placed in developing accurate computer methods capable of predicting this binding event. Prediction of peptide binding to MHC-II is complicated by the open binding cleft of the MHC-II molecule, allowing binding of peptides extending out of the binding groove. Moreover, the genes encoding the MHC molecules are immensely diverse leading to a large set of different MHC molecules each potentially binding a unique set of peptides. Characterizing each MHC-II molecule using peptide-screening binding assays is hence not a viable option.
Here, we present an MHC-II binding prediction algorithm aiming at dealing with these challenges. The method is a pan-specific version of the earlier published allele-specific NN-align algorithm and does not require any pre-alignment of the input data. This allows the method to benefit also from information from alleles covered by limited binding data. The method is evaluated on a large and diverse set of benchmark data, and is shown to significantly out-perform state-of-the-art MHC-II prediction methods. In particular, the method is found to boost the performance for alleles characterized by limited binding data where conventional allele-specific methods tend to achieve poor prediction accuracy.
The method thus shows great potential for efficient boosting the accuracy of MHC-II binding prediction, as accurate predictions can be obtained for novel alleles at highly reduced experimental costs. Pan-specific binding predictions can be obtained for all alleles with know protein sequence and the method can benefit by including data in the training from alleles even where only few binders are known. The method and benchmark data are available at http://www.cbs.dtu.dk/services/NetMHCIIpan-2.0.
The interaction of some human antibodies with heme results in posttranslational acquisition of binding to various self- and pathogen-derived antigens. The previous studies on this phenomenon were ...performed with oxidized heme (Fe
). In the present study, we elucidated the effect of other pathologically relevant species of heme, i.e., species that were formed after contact of heme with oxidizing agents such as hydrogen peroxide, situations in which heme's iron could acquire higher oxidation states. Our data reveal that hyperoxidized species of heme have a superior capacity to heme (Fe
) in triggering the autoreactivity of human IgG. Mechanistic studies demonstrated that oxidation status of iron was of critical importance for the heme's effect on antibodies. We also demonstrated that hyperoxidized heme species interacted at higher affinities with IgG and that this binding occurred through a different mechanism as compared to heme (Fe
). Regardless of their profound functional impact on the antigen-binding properties of antibodies, hyperoxidized species of heme did not affect Fc-mediated functions of IgG, such as binding to the neonatal Fc receptor. The obtained data contribute to a better understanding of the pathophysiological mechanism of hemolytic diseases and of the origin of elevated antibody autoreactivity in patients with some hemolytic disorders.
Glioblastoma (GBM) is resistant to multimodality therapeutic approaches. A high burden of tumor-specific mutant peptides (neoantigens) correlates with better survival and response to immunotherapies ...in selected solid tumors but how neoantigens impact clinical outcome in GBM remains unclear. Here, we exploit the similarity between tumor neoantigens and infectious disease-derived immune epitopes and apply a neoantigen fitness model for identifying high-quality neoantigens in a human pan-glioma dataset. We find that the neoantigen quality fitness model stratifies GBM patients with more favorable clinical outcome and, together with CD8
T lymphocytes tumor infiltration, identifies a GBM subgroup with the longest survival, which displays distinct genomic and transcriptomic features. Conversely, neither tumor neoantigen burden from a quantitative model nor the isolated enrichment of CD8
T lymphocytes were able to predict survival of GBM patients. This approach may guide optimal stratification of GBM patients for maximum response to immunotherapy.