Improved structural method for T-cell cross-reactivity prediction

Unveiling cross-reactivity: implications for immune response modulation in cancer

Antigen recognition by CD8+ T-cell receptors (TCR) is crucial for immune responses to pathogens and tumors. TCRs are cross-reactive, a single TCR can recognize multiple peptide-Human Leukocyte Antigen (HLA) complexes. The study of cross-reactivity can support the development of therapies focusing on immune modulation, such as the expansion of pre-existing T-cell clones to fight pathogens and tumors. The peptide-HLA (pHLA) surface has previously been used to identify TCR cross-reactivities. In the present work, we sought to perform a comprehensive analysis of peptide-HLA by selecting thousands of human and viral epitopes. We profit from established docking models to identify features from different spatial perspectives of HLA-A*02:01, explore similarities between self and non-self epitopes, and list potential cross-reactive epitopes of therapeutic interest. A total of 2631 unique epitopes from representative viral proteins or human proteins were modeled. We were able to demonstrate that cross-reactive CDR3 sequences from public databases recognize epitopes with similar electrostatic potential, charge, and spatial location. Using data from published studies that measured T-cell reactivity to mutated epitopes, we observed a negative correlation between epitope dissimilarity and T-cell activation. Most analysed cancer epitopes were more similar to self epitopes, yet we identified features distinguishing those more similar to viral antigens. Finally, we enumerated potential cross-reactivities between tumoral and viral epitopes and highlighted some challenges in their identification for therapeutic use. Moreover, the thousands of peptide-HLA complexes generated in our work constitute a valuable resource to study T-cell cross-reactivity.

Mutations in glioblastoma proteins do not disrupt epitope presentation and recognition, maintaining a specific CD8 T cell immune response potential

Antigen-specific cytotoxic CD8 T cells are extremely effective in controlling tumor growth and have been the focus of immunotherapy approaches. We leverage in silico tools to investigate whether the occurrence of mutations in proteins previously described as immunogenic and highly expressed by glioblastoma multiforme (GBM), such as Epidermal Growth Factor Receptor (EGFR), Isocitrate Dehydrogenase 1 (IDH1), Phosphatase and Tensin homolog (PTEN) and Tumor Protein 53 (TP53), may be contributing to the differential presentation of immunogenic epitopes. We recovered Class I MHC binding information from wild-type and mutated proteins using the Immune Epitope Database (IEDB). After that, we built peptide-MHC (pMHC-I) models in HLA-arena, followed by hierarchical clustering analysis based on electrostatic surface features from each complex. We identified point mutations that are determinants for the presentation of a set of peptides from TP53 protein. We point to structural features in the pMHC-I complexes of wild-type and mutated peptides, which may play a role in the recognition of CD8 T cells. To further explore these features, we performed 100 ns molecular dynamics simulations for the peptide pairs (wt/mut) selected. In pursuit of novel therapeutic targets for GBM treatment, we selected peptides where our predictive results indicated that mutations would not disrupt epitope presentation, thereby maintaining a specific CD8 T cell immune response. These peptides hold potential for future GBM interventions, including peptide-based or mRNA vaccine development applications.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptides bound to the human MHC, HLA, has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within our curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer pHLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work may be applied towards predicting antigen immunogenicity, and receptor cross-reactivity.

CrossDome: an interactive R package to predict cross-reactivity risk using immunopeptidomics databases

T-cell-based immunotherapies hold tremendous potential in the fight against cancer, thanks to their capacity to specifically targeting diseased cells. Nevertheless, this potential has been tempered with safety concerns regarding the possible recognition of unknown off-targets displayed by healthy cells. In a notorious example, engineered T-cells specific to MAGEA3 (EVDPIGHLY) also recognized a TITIN-derived peptide (ESDPIVAQY) expressed by cardiac cells, inducing lethal damage in melanoma patients. Such off-target toxicity has been related to T-cell cross-reactivity induced by molecular mimicry. In this context, there is growing interest in developing the means to avoid off-target toxicity, and to provide safer immunotherapy products. To this end, we present CrossDome, a multi-omics suite to predict the off-target toxicity risk of T-cell-based immunotherapies. Our suite provides two alternative protocols, i) a peptide-centered prediction, or ii) a TCR-centered prediction. As proof-of-principle, we evaluate our approach using 16 well-known cross-reactivity cases involving cancer-associated antigens. With CrossDome, the TITIN-derived peptide was predicted at the 99+ percentile rank among 36,000 scored candidates (p-value < 0.001). In addition, off-targets for all the 16 known cases were predicted within the top ranges of relatedness score on a Monte Carlo simulation with over 5 million putative peptide pairs, allowing us to determine a cut-off p-value for off-target toxicity risk. We also implemented a penalty system based on TCR hotspots, named contact map (CM). This TCR-centered approach improved upon the peptide-centered prediction on the MAGEA3-TITIN screening (e.g., from 27th to 6th, out of 36,000 ranked peptides). Next, we used an extended dataset of experimentally-determined cross-reactive peptides to evaluate alternative CrossDome protocols. The level of enrichment of validated cases among top 50 best-scored peptides was 63% for the peptide-centered protocol, and up to 82% for the TCR-centered protocol. Finally, we performed functional characterization of top ranking candidates, by integrating expression data, HLA binding, and immunogenicity predictions. CrossDome was designed as an R package for easy integration with antigen discovery pipelines, and an interactive web interface for users without coding experience. CrossDome is under active development, and it is available at https://github.com/AntunesLab/crossdome.

PepSim: T-cell cross-reactivity prediction via comparison of peptide sequence and peptide-HLA structure

Introduction

Peptide-HLA class I (pHLA) complexes on the surface of tumor cells can be targeted by cytotoxic T-cells to eliminate tumors, and this is one of the bases for T-cell-based immunotherapies. However, there exist cases where therapeutic T-cells directed towards tumor pHLA complexes may also recognize pHLAs from healthy normal cells. The process where the same T-cell clone recognizes more than one pHLA is referred to as T-cell cross-reactivity and this process is driven mainly by features that make pHLAs similar to each other. T-cell cross-reactivity prediction is critical for designing T-cell-based cancer immunotherapies that are both effective and safe.

Methods

Here we present PepSim, a novel score to predict T-cell cross-reactivity based on the structural and biochemical similarity of pHLAs.

Results and discussion

We show our method can accurately separate cross-reactive from non-crossreactive pHLAs in a diverse set of datasets including cancer, viral, and self-peptides. PepSim can be generalized to work on any dataset of class I peptide-HLAs and is freely available as a web server at pepsim.kavrakilab.org.

HLA3DB: comprehensive annotation of peptide/HLA complexes enables blind structure prediction of T cell epitopes

The class I proteins of the major histocompatibility complex (MHC-I) display epitopic peptides derived from endogenous proteins on the cell surface for immune surveillance. Accurate modeling of peptide/HLA (pHLA, the human MHC) structures has been mired by conformational diversity of the central peptide residues, which are critical for recognition by T cell receptors. Here, analysis of X-ray crystal structures within a curated database (HLA3DB) shows that pHLA complexes encompassing multiple HLA allotypes present a discrete set of peptide backbone conformations. Leveraging these representative backbones, we employ a regression model trained on terms of a physically relevant energy function to develop a comparative modeling approach for nonamer peptide/HLA structures named RepPred. Our method outperforms the top pHLA modeling approach by up to 19% in terms of structural accuracy, and consistently predicts blind targets not included in our training set. Insights from our work provide a framework for linking conformational diversity with antigen immunogenicity and receptor cross-reactivity.

MatchTope: A tool to predict the cross reactivity of peptides complexed with Major Histocompatibility Complex I

The therapeutic targeting of the immune system, for example in vaccinology and cancer treatment, is a challenging task and the subject of active research. Several in silico tools used for predicting immunogenicity are based on the analysis of peptide sequences binding to the Major Histocompatibility Complex (pMHC). However, few of these bioinformatics tools take into account the pMHC three-dimensional structure. Here, we describe a new bioinformatics tool, MatchTope, developed for predicting peptide similarity, which can trigger cross-reactivity events, by computing and analyzing the electrostatic potentials of pMHC complexes. We validated MatchTope by using previously published data from in vitro assays. We thereby demonstrate the strength of MatchTope for similarity prediction between targets derived from several pathogens as well as for indicating possible cross responses between self and tumor peptides. Our results suggest that MatchTope can enhance and speed up future studies in the fields of vaccinology and cancer immunotherapy.

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

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

SARS-Arena: Sequence and Structure-Guided Selection of Conserved Peptides from SARS-related Coronaviruses for Novel Vaccine Development

The pandemic caused by the SARS-CoV-2 virus, the agent responsible for the COVID-19 disease, has affected millions of people worldwide. There is constant search for new therapies to either prevent or mitigate the disease. Fortunately, we have observed the successful development of multiple vaccines. Most of them are focused on one viral envelope protein, the spike protein. However, such focused approaches may contribute for the rise of new variants, fueled by the constant selection pressure on envelope proteins, and the widespread dispersion of coronaviruses in nature. Therefore, it is important to examine other proteins, preferentially those that are less susceptible to selection pressure, such as the nucleocapsid (N) protein. Even though the N protein is less accessible to humoral response, peptides from its conserved regions can be presented by class I Human Leukocyte Antigen (HLA) molecules, eliciting an immune response mediated by T-cells. Given the increased number of protein sequences deposited in biological databases daily and the N protein conservation among viral strains, computational methods can be leveraged to discover potential new targets for SARS-CoV-2 and SARS-CoV-related viruses. Here we developed SARS-Arena, a user-friendly computational pipeline that can be used by practitioners of different levels of expertise for novel vaccine development. SARS-Arena combines sequence-based methods and structure-based analyses to (i) perform multiple sequence alignment (MSA) of SARS-CoV-related N protein sequences, (ii) recover candidate peptides of different lengths from conserved protein regions, and (iii) model the 3D structure of the conserved peptides in the context of different HLAs. We present two main Jupyter Notebook workflows that can help in the identification of new T-cell targets against SARS-CoV viruses. In fact, in a cross-reactive case study, our workflows identified a conserved N protein peptide (SPRWYFYYL) recognized by CD8⁺ T-cells in the context of HLA-B7⁺. SARS-Arena is available at https://github.com/KavrakiLab/SARS-Arena.

Large-Scale Structure-Based Screening of Potential T Cell Cross-Reactivities Involving Peptide-Targets From BCG Vaccine and SARS-CoV-2

Although not being the first viral pandemic to affect humankind, we are now for the first time faced with a pandemic caused by a coronavirus. The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has been responsible for the COVID-19 pandemic, which caused more than 4.5 million deaths worldwide. Despite unprecedented efforts, with vaccines being developed in a record time, SARS-CoV-2 continues to spread worldwide with new variants arising in different countries. Such persistent spread is in part enabled by public resistance to vaccination in some countries, and limited access to vaccines in other countries. The limited vaccination coverage, the continued risk for resistant variants, and the existence of natural reservoirs for coronaviruses, highlight the importance of developing additional therapeutic strategies against SARS-CoV-2 and other coronaviruses. At the beginning of the pandemic it was suggested that countries with Bacillus Calmette-Guérin (BCG) vaccination programs could be associated with a reduced number and/or severity of COVID-19 cases. Preliminary studies have provided evidence for this relationship and further investigation is being conducted in ongoing clinical trials. The protection against SARS-CoV-2 induced by BCG vaccination may be mediated by cross-reactive T cell lymphocytes, which recognize peptides displayed by class I Human Leukocyte Antigens (HLA-I) on the surface of infected cells. In order to identify potential targets of T cell cross-reactivity, we implemented an in silico strategy combining sequence-based and structure-based methods to screen over 13,5 million possible cross-reactive peptide pairs from BCG and SARS-CoV-2. Our study produced (i) a list of immunogenic BCG-derived peptides that may prime T cell cross-reactivity against SARS-CoV-2, (ii) a large dataset of modeled peptide-HLA structures for the screened targets, and (iii) new computational methods for structure-based screenings that can be used by others in future studies. Our study expands the list of BCG peptides potentially involved in T cell cross-reactivity with SARS-CoV-2-derived peptides, and identifies multiple high-density "neighborhoods" of cross-reactive peptides which could be driving heterologous immunity induced by BCG vaccination, therefore providing insights for future vaccine development efforts.

Mechanistic diversity in MHC class I antigen recognition

Throughout its evolution, the human immune system has developed a plethora of strategies to diversify the antigenic peptide sequences that can be targeted by the CD8+ T cell response against pathogens and aberrations of self. Here we provide a general overview of the mechanisms that lead to the diversity of antigens presented by MHC class I complexes and their recognition by CD8+ T cells, together with a more detailed analysis of recent progress in two important areas that are highly controversial: the prevalence and immunological relevance of unconventional antigen peptides; and cross-recognition of antigenic peptides by the T cell receptors of CD8+ T cells.

Empirical and Rational Design of T Cell Receptor-Based Immunotherapies

The use of T cells reactive with intracellular tumor-associated or tumor-specific antigens has been a promising strategy for cancer immunotherapies in the past three decades, but the approach has been constrained by a limited understanding of the T cell receptor's (TCR) complex functions and specificities. Newer TCR and T cell-based approaches are in development, including engineered adoptive T cells with enhanced TCR affinities, TCR mimic antibodies, and T cell-redirecting bispecific agents. These new therapeutic modalities are exciting opportunities by which TCR recognition can be further exploited for therapeutic benefit. In this review we summarize the development of TCR-based therapeutic strategies and focus on balancing efficacy and potency versus specificity, and hence, possible toxicity, of these powerful therapeutic modalities.

Predicting Cross-Reactivity and Antigen Specificity of T Cell Receptors

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

HLA-Arena: A Customizable Environment for the Structural Modeling and Analysis of Peptide-HLA Complexes for Cancer Immunotherapy

PURPOSE

HLA protein receptors play a key role in cellular immunity. They bind intracellular peptides and display them for recognition by T-cell lymphocytes. Because T-cell activation is partially driven by structural features of these peptide-HLA complexes, their structural modeling and analysis are becoming central components of cancer immunotherapy projects. Unfortunately, this kind of analysis is limited by the small number of experimentally determined structures of peptide-HLA complexes. Overcoming this limitation requires developing novel computational methods to model and analyze peptide-HLA structures.

METHODS

Here we describe a new platform for the structural modeling and analysis of peptide-HLA complexes, called HLA-Arena, which we have implemented using Jupyter Notebook and Docker. It is a customizable environment that facilitates the use of computational tools, such as APE-Gen and DINC, which we have previously applied to peptide-HLA complexes. By integrating other commonly used tools, such as MODELLER and MHCflurry, this environment includes support for diverse tasks in structural modeling, analysis, and visualization.

RESULTS

To illustrate the capabilities of HLA-Arena, we describe 3 example workflows applied to peptide-HLA complexes. Leveraging the strengths of our tools, DINC and APE-Gen, the first 2 workflows show how to perform geometry prediction for peptide-HLA complexes and structure-based binding prediction, respectively. The third workflow presents an example of large-scale virtual screening of peptides for multiple HLA alleles.

CONCLUSION

These workflows illustrate the potential benefits of HLA-Arena for the structural modeling and analysis of peptide-HLA complexes. Because HLA-Arena can easily be integrated within larger computational pipelines, we expect its potential impact to vastly increase. For instance, it could be used to conduct structural analyses for personalized cancer immunotherapy, neoantigen discovery, or vaccine development.

pMHC Structural Comparisons as a Pivotal Element to Detect and Validate T-Cell Targets for Vaccine Development and Immunotherapy-A New Methodological Proposal

The search for epitopes that will effectively trigger an immune response remains the "El Dorado" for immunologists. The development of promising immunotherapeutic approaches requires the appropriate targets to elicit a proper immune response. Considering the high degree of HLA/TCR diversity, as well as the heterogeneity of viral and tumor proteins, this number will invariably be higher than ideal to test. It is known that the recognition of a peptide-MHC (pMHC) by the T-cell receptor is performed entirely in a structural fashion, where the atomic interactions of both structures, pMHC and TCR, dictate the fate of the process. However, epitopes with a similar composition of amino acids can produce dissimilar surfaces. Conversely, sequences with no conspicuous similarities can exhibit similar TCR interaction surfaces. In the last decade, our group developed a database and in silico structural methods to extract molecular fingerprints that trigger T-cell immune responses, mainly referring to physicochemical similarities, which could explain the immunogenic differences presented by different pMHC-I complexes. Here, we propose an immunoinformatic approach that considers a structural level of information, combined with an experimental technology that simulates the presentation of epitopes for a T cell, to improve vaccine production and immunotherapy efficacy.

Machine Learning for Cancer Immunotherapies Based on Epitope Recognition by T Cell Receptors

In the last years, immunotherapies have shown tremendous success as treatments for multiple types of cancer. However, there are still many obstacles to overcome in order to increase response rates and identify effective therapies for every individual patient. Since there are many possibilities to boost a patient's immune response against a tumor and not all can be covered, this review is focused on T cell receptor-mediated therapies. CD8+ T cells can detect and destroy malignant cells by binding to peptides presented on cell surfaces by MHC (major histocompatibility complex) class I molecules. CD4+ T cells can also mediate powerful immune responses but their peptide recognition by MHC class II molecules is more complex, which is why the attention has been focused on CD8+ T cells. Therapies based on the power of T cells can, on the one hand, enhance T cell recognition by introducing TCRs that preferentially direct T cells to tumor sites (so called TCR-T therapy) or through vaccination to induce T cells in vivo. On the other hand, T cell activity can be improved by immune checkpoint inhibition or other means that help create a microenvironment favorable for cytotoxic T cell activity. The manifold ways in which the immune system and cancer interact with each other require not only the use of large omics datasets from gene, to transcript, to protein, and to peptide but also make the application of machine learning methods inevitable. Currently, discovering and selecting suitable TCRs is a very costly and work intensive in vitro process. To facilitate this process and to additionally allow for highly personalized therapies that can simultaneously target multiple patient-specific antigens, especially neoepitopes, breakthrough computational methods for predicting antigen presentation and TCR binding are urgently required. Particularly, potential cross-reactivity is a major consideration since off-target toxicity can pose a major threat to patient safety. The current speed at which not only datasets grow and are made available to the public, but also at which new machine learning methods evolve, is assuring that computational approaches will be able to help to solve problems that immunotherapies are still facing.

T-cell-receptor cross-recognition and strategies to select safe T-cell receptors for clinical translation

Adoptive transfer of T-cell-receptor (TCR)-transduced T cells has shown promising results for cancer treatment, but has also produced severe immunotoxicities caused by on-target as well as off-target TCR recognition. Off-target toxicities are related to the ability of a single T cell to cross-recognize and respond to several different peptide–major histocompatibility complex (pMHC) antigens; a property that is essential for providing broad antigenic coverage despite a confined number of unique TCRs in the human body. However, this degeneracy makes it incredibly difficult to account for the range of targets that any TCR might recognize, which represents a major challenge for the clinical development of therapeutic TCRs. The prospect of using affinity-optimized TCRs has been impeded due to observations that affinity enhancement might alter the specificity of a TCR, thereby increasing the risk that it will cross-recognize endogenous tissue. Strategies for selecting safe TCRs for the clinic have included functional assessment after individual incubations with tissue-derived primary cells or with peptides substituted with single amino acids. However, these strategies have not been able to predict cross-recognition sufficiently, leading to fatal cross-reactivity in clinical trials. Novel technologies have emerged that enable extensive characterization of the exact interaction points of a TCR with pMHC, which provides a foundation from which to make predictions of the cross-recognition potential of individual TCRs. This review describes current advances in strategies for dissecting the molecular interaction points of TCRs, focusing on their potential as tools for predicting cross-recognition of TCRs in clinical development.

•

T-cell-receptor (TCR) degeneracy plays a fundamental role in the capacity of our immune systems to recognize foreign antigens.

•

TCR cross-reactivity provides an inherent risk in TCR–gene transfer cell therapies.

•

Advances in description of TCR cross-recognition can guide the selection process for TCRs into clinical use.

T-Cell Receptor Cognate Target Prediction Based on Paired α and β Chain Sequence and Structural CDR Loop Similarities

T-cell receptors (TCR) mediate immune responses recognizing peptides in complex with major histocompatibility complexes (pMHC) displayed on the surface of cells. Resolving the challenge of predicting the cognate pMHC target of a TCR would benefit many applications in the field of immunology, including vaccine design/discovery and the development of immunotherapies. Here, we developed a model for prediction of TCR targets based on similarity to a database of TCRs with known targets. Benchmarking the model on a large set of TCRs with known target, we demonstrated how the predictive performance is increased (i) by focusing on CDRs rather than the full length TCR protein sequences, (ii) by incorporating information from paired α and β chains, and (iii) integrating information for all 6 CDR loops rather than just CDR3. Finally, we show how integration of the structure of CDR loops, as obtained through homology modeling, boosts the predictive power of the model, in particular in situations where no high-similarity TCRs are available for the query. These findings demonstrate that TCRs that bind to the same target also share, to a very high degree, sequence, and structural features. This observation has profound impact for future development of prediction models for TCR-pMHC interactions and for the use of such models for the rational design of T cell based therapies.

A new in silico approach to investigate molecular aspects of factor IX missense causative mutations and their impact on the hemophilia B severity

Factor IX (encoded by F9) is a protein in the coagulation process, where its lack or deficiency leads to hemophilia B. This condition has been much less studied than hemophilia A, especially in Latin America. We analyzed the structural and functional impact of 54 missense mutations (18 reported by us previously, and 36 other mutations from the Factor IX database) through molecular modeling approaches. To accomplish this task, we examine the electrostatic patterns, hydrophobicity/hydrophilicity, disulfide, and H-bond differences of the Factor IX structures harboring the missense mutations found, correlating them with their clinical effects. The 54 mutated sequences were modeled and their physicochemical features were determined and used as input in clusterization tools. The electrostatic pattern seems to influence in disease severity, especially for mutations investigated in epidermal growth factors 1 and 2 (EGF1/2) domains. The combined use of all physicochemical information improved the clustering of structures associated to similar phenotypes, especially for mutations from GLA and EGF1-2 domains. The effect of mutations in the disease phenotype severity seems to be a complex interplay of molecular features, each one contributing to different impacts. This highlights that previous studies and tools analyzing individually single features for single mutations are missing elements that fulfill the whole picture.

General Prediction of Peptide-MHC Binding Modes Using Incremental Docking: A Proof of Concept

The class I major histocompatibility complex (MHC) is capable of binding peptides derived from intracellular proteins and displaying them at the cell surface. The recognition of these peptide-MHC (pMHC) complexes by T-cells is the cornerstone of cellular immunity, enabling the elimination of infected or tumoral cells. T-cell-based immunotherapies against cancer, which leverage this mechanism, can greatly benefit from structural analyses of pMHC complexes. Several attempts have been made to use molecular docking for such analyses, but pMHC structure remains too challenging for even state-of-the-art docking tools. To overcome these limitations, we describe the use of an incremental meta-docking approach for structural prediction of pMHC complexes. Previous methods applied in this context used specific constraints to reduce the complexity of this prediction problem, at the expense of generality. Our strategy makes no assumption and can potentially be used to predict binding modes for any pMHC complex. Our method has been tested in a re-docking experiment, reproducing the binding modes of 25 pMHC complexes whose crystal structures are available. This study is a proof of concept that incremental docking strategies can lead to general geometry prediction of pMHC complexes, with potential applications for immunotherapy against cancer or infectious diseases.

Interpreting T-Cell Cross-reactivity through Structure: Implications for TCR-Based Cancer Immunotherapy

Immunotherapy has become one of the most promising avenues for cancer treatment, making use of the patient’s own immune system to eliminate cancer cells. Clinical trials with T-cell-based immunotherapies have shown dramatic tumor regressions, being effective in multiple cancer types and for many different patients. Unfortunately, this progress was tempered by reports of serious (even fatal) side effects. Such therapies rely on the use of cytotoxic T-cell lymphocytes, an essential part of the adaptive immune system. Cytotoxic T-cells are regularly involved in surveillance and are capable of both eliminating diseased cells and generating protective immunological memory. The specificity of a given T-cell is determined through the structural interaction between the T-cell receptor (TCR) and a peptide-loaded major histocompatibility complex (MHC); i.e., an intracellular peptide–ligand displayed at the cell surface by an MHC molecule. However, a given TCR can recognize different peptide–MHC (pMHC) complexes, which can sometimes trigger an unwanted response that is referred to as T-cell cross-reactivity. This has become a major safety issue in TCR-based immunotherapies, following reports of melanoma-specific T-cells causing cytotoxic damage to healthy tissues (e.g., heart and nervous system). T-cell cross-reactivity has been extensively studied in the context of viral immunology and tissue transplantation. Growing evidence suggests that it is largely driven by structural similarities of seemingly unrelated pMHC complexes. Here, we review recent reports about the existence of pMHC “hot-spots” for cross-reactivity and propose the existence of a TCR interaction profile (i.e., a refinement of a more general TCR footprint in which some amino acid residues are more important than others in triggering T-cell cross-reactivity). We also make use of available structural data and pMHC models to interpret previously reported cross-reactivity patterns among virus-derived peptides. Our study provides further evidence that structural analyses of pMHC complexes can be used to assess the intrinsic likelihood of cross-reactivity among peptide-targets. Furthermore, we hypothesize that some apparent inconsistencies in reported cross-reactivities, such as a preferential directionality, might also be driven by particular structural features of the targeted pMHC complex. Finally, we explain why TCR-based immunotherapy provides a special context in which meaningful T-cell cross-reactivity predictions can be made.