Predicting HLA class I non-permissive amino acid residues substitutions

Comparative analysis of the variability of the human leukocyte antigen peptide-binding pockets in patients with acute leukaemia

The association between acute lymphoblastic leukaemia (ALL) and acute myeloid leukaemia (AML) and the human leukocyte antigens (HLA) has rarely been studied in terms of diversity of peptide-binding pockets. The objective of this study was to analyse whether motifs of HLA class I and class II peptide-binding pockets and/or their amino acid positions were differentially associated with ALL and AML. We included 849 patients from the Eurocord/European Blood and Marrow Transplant registry. The HLA peptide-binding pockets whose amino acid variability was analysed were B and F for HLA class I, P4, P6, and P9 for HLA-DRB1, and P4 and P9 for HLA-DQB1. The motif RFDRAY in P4 of HLA-DRB1*16:01/02/03/05 alleles and the motif YYVSY in P9 of HLA-DQB1*05:02/04/05 alleles, were statistically associated with ALL (corrected p value [p_c ] = 0.001 and p_c  = 0.035 respectively). The frequency of serine 57 in the P9 of HLA-DQB1 was higher in ALL (odds ratio 2.09, 95% confidence interval: 1.27-3.44; p_c  = 0.037). Our analysis suggests that specific motifs in terms of HLA class II pockets and amino acids might be unique to ALL. The associations identified in this study encourage further investigation oF the role of HLA peptide-binding pockets and their amino acids in immune processes underpinning acute leukaemia and ultimately in immunotherapy settings.

Molecular mimicry between tumor associated antigens and microbiota-derived epitopes

Background

The gut microbiota profile is unique for each individual and are composed by different bacteria species according to individual birth-to-infant transitions. In the last years, the local and systemic effects of microbiota on cancer onset, progression and response to treatments, such as immunotherapies, has been extensively described. Here we offer a new perspective, proposing a role for the microbiota based on the molecular mimicry of tumor associated antigens by microbiome-associated antigens.

Methods

In the present study we looked for homology between published TAAs and non-self microbiota-derived epitopes. Blast search for sequence homology was combined with extensive bioinformatics analyses.

Results

Several evidences for homology between TAAs and microbiota-derived antigens have been found. Strikingly, three cases of 100% homology between the paired sequences has been identified. The predicted average affinity to HLA molecules of microbiota-derived antigens is very high (< 100 nM). The structural conformation of the microbiota-derived epitopes is, in general, highly similar to the corresponding TAA. In some cases, it is identical and contact areas with both HLA and TCR chains are indistinguishable. Moreover, the spatial conformation of TCR-facing residues can be identical in paired TAA and microbiota-derived epitopes, with exactly the same values of planar as well as dihedral angles.

Conclusions

The data reported in the present study show for the first time the high homology in the linear sequence as well as in structure and conformation between TAAs and peptides derived from microbiota species of the Firmicutes and the Bacteroidetes phyla, which together account for 90% of gut microbiota. Cross-reacting CD8⁺ T cell responses are very likely induced. Therefore, the anti-microbiota T cell memory may turn out to be an anti-cancer T cell memory, able to control the growth of a cancer developed during the lifetime if the expressed TAA is similar to the microbiota epitope. This may ultimately represent a relevant selective advantage for cancer patients and may lead to a novel preventive anti-cancer vaccine strategy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03512-6.

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

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

Long-term memory T cells as preventive anticancer immunity elicited by TuA-derived heteroclitic peptides

The host's immune system may be primed against antigens during the lifetime (e.g. microorganisms antigens-MoAs), and swiftly recalled upon growth of a tumor expressing antigens similar in sequence and structure. C57BL/6 mice were immunized in a preventive setting with tumor antigens (TuAs) or corresponding heteroclitic peptides specific for TC-1 and B16 cell lines. Immediately or 2-months after the end of the vaccination protocol, animals were implanted with cell lines. The specific anti-vaccine immune response as well as tumor growth were regularly evaluated for 2 months post-implantation. The preventive vaccination with TuA or their heteroclitic peptides (hPep) was able to delay (B16) or completely suppress (TC-1) tumor growth when cancer cells were implanted immediately after the end of the vaccination. More importantly, TC-1 tumor growth was significantly delayed, and suppressed in 6/8 animals, also when cells were implanted 2-months after the end of the vaccination. The vaccine-specific T cell response provided a strong immune correlate to the pattern of tumor growth. A preventive immunization with heteroclitic peptides resembling a TuA is able to strongly delay or even suppress tumor growth in a mouse model. More importantly, the same effect is observed also when tumor cells are implanted 2 months after the end of vaccination, which corresponds to 8 - 10 years in human life. The observed potent tumor control indicates that a memory T cell immunity elicited during the lifetime by a antigens similar to a TuA, i.e. viral antigens, may ultimately represent a great advantage for cancer patients and may lead to a novel preventive anti-cancer vaccine strategy.

MHC-Optimized Peptide Scaffold for Improved Antigen Presentation and Anti-Tumor Response

Tumor Associated Antigens (TAAs) may suffer from an immunological tolerance due to expression on normal cells. In order to potentiate their immunogenicity, heteroclitic peptides (htcPep) were designed according to prediction algorithms. In particular, specific modifications were introduced in peptide residues facing to TCR. Moreover, a MHC-optimized scaffold was designed for improved antigen presentation to TCR by H-2Db allele. The efficacy of such htcPep was assessed in C57BL/6 mice injected with syngeneic melanoma B16F10 or lung TC1 tumor cell lines, in combination with metronomic chemotherapy and immune checkpoint inhibitors. The immunogenicity of htcPep was significantly stronger than the corresponding wt peptide and the modification involving both MHC and TCR binding residues scored the strongest. In particular, the H-2Db-specific scaffold significantly potentiated the peptides' immunogenicity and control of tumor growth was comparable to wt peptide in a therapeutic setting. Overall, we demonstrated that modified TAAs show higher immunogenicity compared to wt peptide. In particular, the MHC-optimized scaffold can present different antigen sequences to TCR, retaining the conformational characteristics of the corresponding wt. Cross-reacting CD8+ T cells are elicited and efficiently kill tumor cells presenting the wild-type antigen. This novel approach can be of high clinical relevance in cancer vaccine development.

Identification and characterization of heteroclitic peptides in TCR-binding positions with improved HLA-binding efficacy

The antigenicity as well as the immunogenicity of tumor associated antigens (TAAs) may need to be potentiated in order to break the immunological tolerance. To this aim, heteroclitic peptides were designed introducing specific substitutions in the residue at position 4 (p4) binding to TCR. The effect of such modifications also on the affinity to the major histocompatibility class I (MHC-I) molecule was assessed. The Trp2 antigen, specific for the mouse melanoma B16F10 cells, as well as the HPV-E7 antigen, specific for the TC1 tumor cell lines, were used as models. Affinity of such heteroclitic peptides to HLA was predicted by bioinformatics tools and the most promising ones were validated by structural conformational and HLA binding analyses. Overall, we demonstrated that TAAs modified at the TCR-binding p4 residue are predicted to have higher affinity to MHC-I molecules. Experimental evaluation confirms the stronger binding, suggesting that this strategy may be very effective for designing new vaccines with improved antigenic efficacy.

HLA-A29 and Birdshot Uveitis: Further Down the Rabbit Hole

HLA class I alleles constitute established risk factors for non-infectious uveitis and preemptive genotyping of HLA class I alleles is standard practice in the diagnostic work-up. The HLA-A29 serotype is indispensable to Birdshot Uveitis (BU) and renders this enigmatic eye condition a unique model to better understand how the antigen processing and presentation machinery contributes to non-infectious uveitis or chronic inflammatory conditions in general. This review will discuss salient points regarding the protein structure of HLA-A29 and how key amino acid positions impact the peptide binding preference and interaction with T cells. We discuss to what extent the risk genes ERAP1 and ERAP2 uniquely affect HLA-A29 and how the discovery of a HLA-A29-specific submotif may impact autoantigen discovery. We further provide a compelling argument to solve the long-standing question why BU only affects HLA-A29-positive individuals from Western-European ancestry by exploiting data from the 1000 Genomes Project. We combine novel insights from structural and immunopeptidomic studies and discuss the functional implications of genetic associations across the HLA class I antigen presentation pathway to refine the etiological basis of Birdshot Uveitis.

Selecting Target Antigens for Cancer Vaccine Development

One of the principal goals of cancer immunotherapy is the development of efficient therapeutic cancer vaccines that are able to elicit an effector as well as memory T cell response specific to tumor antigens. In recent years, the attention has been focused on the personalization of cancer vaccines. However, the efficacy of therapeutic cancer vaccines is still disappointing despite the large number of vaccine strategies targeting different tumors that have been evaluated in recent years. While the preclinical data have frequently shown encouraging results, clinical trials have not provided satisfactory data to date. The main reason for such failures is the complexity of identifying specific target tumor antigens that should be unique or overexpressed only by the tumor cells compared to normal cells. Most of the tumor antigens included in cancer vaccines are non-mutated overexpressed self-antigens, eliciting mainly T cells with low-affinity T cell receptors (TCR) unable to mediate an effective anti-tumor response. In this review, the target tumor antigens employed in recent years in the development of therapeutic cancer vaccine strategies are described, along with potential new classes of tumor antigens such as the human endogenous retroviral elements (HERVs), unconventional antigens, and/or heteroclitic peptides.

Immunoinformatics and Vaccine Development: An Overview

The use of vaccines have resulted in a remarkable improvement in global health. It has saved several lives, reduced treatment costs and raised the quality of animal and human lives. Current traditional vaccines came empirically with either vague or completely no knowledge of how they modulate our immune system. Even at the face of potential vaccine design advance, immune-related concerns (as seen with specific vulnerable populations, cases of emerging/re-emerging infectious disease, pathogens with complex lifecycle and antigenic variability, need for personalized vaccinations, and concerns for vaccines' immunological safety -specifically vaccine likelihood to trigger non-antigen-specific responses that may cause autoimmunity and vaccine allergy) are being raised. And these concerns have driven immunologists toward research for a better approach to vaccine design that will consider these challenges. Currently, immunoinformatics has paved the way for a better understanding of some infectious disease pathogenesis, diagnosis, immune system response and computational vaccinology. The importance of this immunoinformatics in the study of infectious diseases is diverse in terms of computational approaches used, but is united by common qualities related to host–pathogen relationship. Bioinformatics methods are also used to assign functions to uncharacterized genes which can be targeted as a candidate in vaccine design and can be a better approach toward the inclusion of women that are pregnant into vaccine trials and programs. The essence of this review is to give insight into the need to focus on novel computational, experimental and computation-driven experimental approaches for studying of host–pathogen interactions and thus making a case for its use in vaccine development.

GradDock: rapid simulation and tailored ranking functions for peptide-MHC Class I docking

Motivation

The identification of T-cell epitopes has many profound translational applications in the areas of transplantation, disease diagnosis, vaccine/therapeutic protein development and personalized immunotherapy. While data-driven methods have been widely used for the prediction of peptide binders with notable successes, the structural modeling of peptide binding to MHC molecules is crucial for understanding the underlying molecular mechanism of the immunological processes.

Results

We developed GradDock, a structure-based method for the rapid and accurate modeling of peptide binding to MHC Class I (pMHC-I). GradDock explicitly models diverse unbound peptides in vacuo and inserts them into the MHC-I groove through a steered gradient descent with a topological correction process. The simulation process yields diverse structural conformations including native-like peptides. We completely revised the Rosetta score terms and developed a new ranking function specifically for pMHC-I. Using the diverse peptides, a linear programming approach is applied to find the optimal weights for the individual Rosetta score terms. Our examination revealed that a refinement of the dihedral angles and a modification of the repulsion can dramatically improve the modeling quality. GradDock is five-times faster than a Rosetta-based docking approach for pMHC-I. We also demonstrate that the predictive capability of GradDock with the re-weighted Rosetta ranking function is consistently more accurate than the Rosetta-based method with the standard Rosetta score (approximately three-times better for a cross-docking set).

Availability and implementation

GradDock is freely available for academic purposes. The program and the ranking score weights for Rosetta are available at http://bel.kaist.ac.kr/research/GradDock.

Contact

hskim76@kaist.ac.kr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Unique true predicted neoantigens (TPNAs) correlates with anti-tumor immune control in HCC patients

Background

A novel prediction algorithm is needed for the identification of effective tumor associated mutated neoantigens. Only those with no homology to self wild type antigens are true predicted neoantigens (TPNAs) and can elicit an antitumor T cell response, not attenuated by central tolerance. To this aim, the mutational landscape was evaluated in HCV-associated hepatocellular carcinoma.

Methods

Liver tumor biopsies and adjacent non-tumor liver tissues were obtained from 9 HCV-chronically infected subjects and subjected to RNA-Seq analysis. Mutant peptides were derived from single nucleotide variations and TPNAs were predicted using two prediction servers (e.g. NetTepi and NetMHCstabpan) by comparison with corresponding wild-type sequences, non-related self and pathogen-related antigens. Immunological confirmation was obtained in preclinical as well as clinical setting.

Results

The development of such an improved algorithm resulted in a handful of TPNAs despite the large number of predicted neoantigens. Furthermore, TPNAs may share homology to pathogen’s antigens and be targeted by a pre-existing T cell immunity. Cross-reactivity between such antigens was confirmed in an experimental pre-clinical setting. Finally, TPNAs homologous to pathogen’s antigens were found in the only HCC long-term survival patient, suggesting a correlation between the pre-existing T cell immunity specific for these TPNAs and the favourable clinical outcome.

Conclusions

The new algorithm allowed the identification of the very few TPNAs in cancer cells, and those targeted by a pre-existing immunity strongly correlated with long-term survival. Only such TPNAs represent the optimal candidates for immunotherapy strategies.

Electronic supplementary material

The online version of this article (10.1186/s12967-018-1662-9) contains supplementary material, which is available to authorized users.

Identification of high-risk amino-acid substitutions in hematopoietic cell transplantation: a challenging task

Allogeneic hematopoietic cell transplantation (HCT) offers the potential to cure hematologic malignancies. In the absence of an HLA-matched donor, HLA mismatched unrelated donors may be used, although risks of GvHD and treatment-related mortality (TRM) are higher. Identification and avoidance of amino-acid substitution and position types (AASPT) conferring higher risks of TRM and GvHD would potentially improve the success of transplantation from single HLA mismatched unrelated donors. Using random forest and logistic regression analyses, we identified 19 AASPT associated with greater risks for at least one adverse transplant outcome: grade III-IV acute GvHD, TRM, lower disease-free survival or worse overall survival relative to HLA-matched unrelated donors and to other AASPT. When tested in an independent validation cohort of 3530 patients, none of the AASPT from the training set were validated as high risk, however. Review of the literature shows that failure to validate original observations is the rule and not the exception in immunobiology and emphasizes the importance of independent validation before clinical application. Our current data do not support avoiding any specific class I AASPT for unrelated donors. Additional studies should be performed to fully understand the role of AASPT in HCT outcomes.

Functional distance between recipient and donor HLA-DPB1 determines nonpermissive mismatches in unrelated HCT

The role of HLA amino acid (AA) polymorphism for the outcome of hematopoietic cell transplantation (HCT) is controversial, in particular for HLA class II. Here, we investigated this question in nonpermissive HLA-DPB1 T-cell epitope (TCE) mismatches reflected by numerical functional distance (FD) scores, assignable to all HLA-DPB1 alleles based on the combined impact of 12 polymorphic AAs. We calculated the difference in FD scores (ΔFD) of mismatched HLA-DPB1 alleles in patients and their 10/10 HLA-matched unrelated donors of 379 HCTs performed at our center for acute leukemia or myelodysplastic syndrome. Receiver-operator curve-based stratification into 2 ΔFD subgroups showed a significantly higher percentage of nonpermissive TCE mismatches for ΔFD >2.665, compared with ΔFD ≤2.665 (88% vs 25%, P < .0001). In multivariate analysis, ΔFD >2.665 was significantly associated with overall survival (hazard ratio [HR], 1.40; 95% confidence interval [CI], 1.05-1.87; P < .021) and event-free survival (HR, 1.39; 95% CI, 1.05-1.82; P < .021), compared with ΔFD ≤2.665. These associations were stronger than those observed for TCE mismatches. There was a marked but not statistically significant increase in the hazards of relapse and nonrelapse mortality in the high ΔFD subgroup, whereas no differences were observed for acute and chronic graft-versus-host disease. Seven nonconservative AA substitutions in peptide-binding positions had a significantly stronger impact on ΔFD compared with 5 others (P = .0025), demonstrating qualitative differences in the relative impact of AA polymorphism in HLA-DPB1. The novel concept of ΔFD sheds new light onto nonpermissive HLA-DPB1 mismatches in unrelated HCT.

Key role of human leukocyte antigen in modulating human immunodeficiency virus progression: An overview of the possible applications

Host and viral factors deeply influence the human immunodeficiency virus (HIV) disease progression. Among them human leukocyte antigen (HLA) locus plays a key role at different levels. In fact, genes of the HLA locus have shown the peculiar capability to modulate both innate and adaptive immune responses. In particular, HLA class I molecules are recognized by CD8⁺ T-cells and natural killers (NK) cells towards the interaction with T cell receptor (TCR) and Killer Immunoglobulin Receptor (KIR) 3DL1 respectively. Polymorphisms within the different HLA alleles generate structural changes in HLA class I peptide-binding pockets. Amino acid changes in the peptide-binding pocket lead to the presentation of a different set of peptides to T and NK cells. This review summarizes the role of HLA in HIV progression toward acquired immunodeficiency disease syndrome and its receptors. Recently, many studies have been focused on determining the HLA binding-peptides. The novel use of immune-informatics tools, from the prediction of the HLA-bound peptides to the modification of the HLA-receptor complexes, is considered. A better knowledge of HLA peptide presentation and recognition are allowing new strategies for immune response manipulation to be applied against HIV virus.

HLA-C Incompatibilities in Allogeneic Unrelated Hematopoietic Stem Cell Transplantation

An increasingly larger fraction of patients with hematological diseases are treated by hematopoietic stem cells transplantation (HSCT) from HLA matched unrelated donors. Polymorphisms of HLA genes represent a major barrier to HSCT because HLA-A, -B, -C and DRB1 incompatibilities confer a higher risk of acute graft-versus-host disease (aGVHD) and mortality. Although >22 million volunteer HLA-typed donors are available worldwide, still a significant number of patients do not find a highly matched HSC donor. Because of the large haplotypic diversity in HLA-B–C associations, incompatibilities occur most frequently at HLA-C, so that unrelated donors with a single HLA-C mismatch often represent the only possible choice. The ratio of HLA-C-mismatched HSCT over the total number of transplants varies from 15 to 30%, as determined in 12 multicenter studies. Six multicenter studies involving >1800 patients have reported a 21–43% increase in mortality risk. By using in vitro cellular assays, a large heterogeneity in T-cell allorecognition has been observed. Yet the permissiveness of individual HLA-C mismatches remains poorly defined. It could be linked to the position and nature of the mismatched residues on HLA-C molecules, but also to variability in the expression levels of the mismatched alleles. The permissive C*03:03–03:04 mismatch is characterized by full compatibility at residues 9, 97, 99, 116, 152, 156, and 163 reported to be key positions influencing T-cell allorecognition. With a single difference among these seven key residues the C*07:01–07:02 mismatch might also be considered by analogy as permissive. High variability of HLA-C expression as determined by quantitative RT-PCR has been observed within individual allotypes and shows some correlation with A–B–C–DRB1 haplotypes. Thus in addition to the position of mismatched amino acid residues, expression level of patient’s mismatched HLA-C allotype might influence T-cell allorecognition, with patients low expression-C alleles representing possible permissive mismatches.

MHC universal cells survive in an allogeneic environment after incompatible transplantation

Cell, tissue, and organ transplants are commonly performed for the treatment of different diseases. However, major histocompatibility complex (MHC) diversity often prevents complete donor-recipient matching, resulting in graft rejection. This study evaluates in a preclinical model the capacity of MHC class I-silenced cells to engraft and grow upon allogeneic transplantation. Short hairpin RNA targeting β2-microglobulin (RN_shβ2m) was delivered into fibroblasts derived from LEW/Ztm (RT1^l) (RT1-A^l) rats using a lentiviral-based vector. MHC class I (RT1-A-) expressing and -silenced cells were injected subcutaneously in LEW rats (RT1^l) and MHC-congenic LEW.1W rats (RT1^u), respectively. Cell engraftment and the status of the immune response were monitored for eight weeks after transplantation. In contrast to RT1-A-expressing cells, RT1-A-silenced fibroblasts became engrafted and were still detectable eight weeks after allogeneic transplantation. Plasma levels of proinflammatory cytokines IL-1α, IL-1β, IL-6, TNF-α, and IFN-γ were significantly higher in animals transplanted with RT1-A-expressing cells than in those receiving RT1-A-silenced cells. Furthermore, alloantigen-specific T-cell proliferation rates derived from rats receiving RT1-A-expressing cells were higher than those in rats transplanted with RT1-A-silenced cells. These data suggest that silencing MHC class I expression might overcome the histocompatibility barrier, potentially opening up new avenues in the field of cell transplantation and regenerative medicine.

Immunoinformatic docking approach for the analysis of KIR3DL1/HLA-B interaction

KIR3DL1 is among the most interesting receptors studied, within the killer immunoglobulin receptor (KIR) family. Human leukocyte antigen (HLA) class I Bw4 epitope inhibits strongly Natural Killer (NK) cell's activity through interaction with KIR3DL1 receptor, while Bw6 generally does not. This interaction has been indicated to play an important role in the immune control of different viral infectious diseases. However, the structural interaction between the KIR3DL1 receptor and different HLA-B alleles has been scarcely studied. To understand the complexity of KIR3DL1-HLA-B interaction, HLA-B alleles carrying Bw4/Bw6 epitope and KIR3DL1∗001 allele in presence of different peptides has been evaluated by using a structural immunoinformatic approach. Different energy minimization force fields (ff) have been tested and NOVA ff enables the successful prediction of ligand-receptor interaction. HLA-B alleles carrying Bw4 epitope present the highest capability of interaction with KIR3DL1∗001 compared to the HLA-B alleles presenting Bw6. The presence of the epitope Bw4 determines a conformational change which leads to a stronger interaction between nonpolymorphic arginine at position 79 of HLA-B and KIR3DL1∗001 136–142 loop. The data shed new light on the modalities of KIR3DL1 interaction with HLA-B alleles essential for the modulation of NK immune-mediated response.