TCR3d: The T cell receptor structural repertoire database

A naturally selected αβ T cell receptor binds HLA-DQ2 molecules without co-contacting the presented peptide

αβ T cell receptors (TCR) co-recognise peptide (p) antigens that are presented by major histocompatibility complex (MHC) molecules. While marked variations in TCR-p-MHC docking topologies have been observed from structural studies, the co-recognition paradigm has held fast. Using HLA-DQ2.5-peptide tetramers, here we identify a TRAV12-1+-TRBV5-1+ G9 TCR from human peripheral blood that binds HLA-DQ2.5 in a peptide-agnostic manner. The crystal structures of TCR-HLA-DQ2.5-peptide complexes show that the G9 TCR binds HLA-DQ2.5 in a reversed docking topology without contacting the peptide, with the TCR contacting the β1 region of HLA-DQ2.5 and distal from the peptide antigen binding cleft. High-throughput screening of HLA class I and II molecules finds the G9 TCR to be pan-HLA-DQ2 reactive, with leucine-55 of HLA-DQ2.5 being a key determinant underpinning G9 TCR specificity excluding other HLA-II allomorphs. Consistent with the functional assays, the interactions of the G9 TCR and HLA-DQ2.5 precludes CD4 binding, thereby impeding T cell activation. Collectively, we describe a naturally selected αβTCR from human peripheral blood that deviates from the TCR-p-MHC co-recognition paradigm.

TCR3d 2.0: expanding the T cell receptor structure database with new structures, tools and interactions

Recognition of antigens by T cell receptors (TCRs) is a key component of adaptive immunity. Understanding the structures of these TCR interactions provides major insights into immune protection and diseases, and enables design of therapeutics, vaccines and predictive modeling algorithms. Previously, we released TCR3d, a database and resource for structures of TCRs and their recognition. Due to the growth of available structures and categories of complexes, the content of TCR3d has expanded substantially in the past 5 years. This expansion includes new tables dedicated to TCR mimic antibody complex structures, TCR-CD3 complexes and annotated Class I and II peptide-MHC complexes. Additionally, tools are available for users to calculate docking geometries for input TCR and TCR mimic complex structures. The core tables of TCR-peptide-MHC complexes have grown by 50%, and include binding affinity data for experimentally determined structures. These major content and feature updates enhance TCR3d as a resource for immunology, therapeutics and structural biology research, and enable advanced approaches for predictive TCR modeling and design. TCR3d is available at: https://tcr3d.ibbr.umd.edu.

Cracking the code of adaptive immunity: The role of computational tools

In recent years, the advances in high-throughput and deep sequencing have generated a diverse amount of adaptive immune repertoire data. This surge in data has seen a proportional increase in computational methods aimed to characterize T cell receptor (TCR) repertoires. In this perspective, we will provide a brief commentary on the various domains of TCR repertoire analysis, their respective computational methods, and the ongoing challenges. Given the breadth of methods and applications of TCR analysis, we will focus our perspective on sequence-based computational methods.

Strengths and limitations of web servers for the modeling of TCRpMHC complexes

Cellular immunity relies on the ability of a T-cell receptor (TCR) to recognize a peptide (p) presented by a class I major histocompatibility complex (MHC) receptor on the surface of a cell. The TCR-peptide-MHC (TCRpMHC) interaction is a crucial step in activating T-cells, and the structural characteristics of these molecules play a significant role in determining the specificity and affinity of this interaction. Hence, obtaining 3D structures of TCRpMHC complexes offers valuable insights into various aspects of cellular immunity and can facilitate the development of T-cell-based immunotherapies. Here, we aimed to compare three popular web servers for modeling the structures of TCRpMHC complexes, namely ImmuneScape (IS), TCRpMHCmodels, and TCRmodel2, to examine their strengths and limitations. Each method employs a different modeling strategy, including docking, homology modeling, and deep learning. The accuracy of each method was evaluated by reproducing the 3D structures of a dataset of 87 TCRpMHC complexes with experimentally determined crystal structures available on the Protein Data Bank (PDB). All selected structures were limited to human MHC alleles, presenting a diverse set of peptide ligands. A detailed analysis of produced models was conducted using multiple metrics, including Root Mean Square Deviation (RMSD) and standardized assessments from CAPRI and DockQ. Special attention was given to the complementarity-determining region (CDR) loops of the TCRs and to the peptide ligands, which define most of the unique features and specificity of a given TCRpMHC interaction. Our study provides an optimistic view of the current state-of-the-art for TCRpMHC modeling but highlights some remaining challenges that must be addressed in order to support the future application of these tools for TCR engineering and computer-aided design of TCR-based immunotherapies.

Structural characterization and AlphaFold modeling of human T cell receptor recognition of NRAS cancer neoantigens

T cell receptors (TCRs) that recognize cancer neoantigens are important for anticancer immune responses and immunotherapy. Understanding the structural basis of TCR recognition of neoantigens provides insights into their exquisite specificity and can enable design of optimized TCRs. We determined crystal structures of a human TCR in complex with NRAS Q61K and Q61R neoantigen peptides and HLA-A1 major histocompatibility complex (MHC), revealing the molecular underpinnings for dual recognition and specificity versus wild-type NRAS peptide. We then used multiple versions of AlphaFold to model the corresponding complex structures, given the challenge of immune recognition for such methods. One implementation of AlphaFold2 (TCRmodel2) with additional sampling was able to generate accurate models of the complexes, while AlphaFold3 also showed strong performance, although success was lower for other complexes. This study provides insights into TCR recognition of a shared cancer neoantigen as well as the utility and practical considerations for using AlphaFold to model TCR-peptide-MHC complexes.

Predicting Antigen-Specificities of Orphan T Cell Receptors from Cancer Patients with TCRpcDist

Approaches to analyze and cluster T-cell receptor (TCR) repertoires to reflect antigen specificity are critical for the diagnosis and prognosis of immune-related diseases and the development of personalized therapies. Sequence-based approaches showed success but remain restrictive, especially when the amount of experimental data used for the training is scarce. Structure-based approaches which represent powerful alternatives, notably to optimize TCRs affinity toward specific epitopes, show limitations for large-scale predictions. To handle these challenges, TCRpcDist is presented, a 3D-based approach that calculates similarities between TCRs using a metric related to the physico-chemical properties of the loop residues predicted to interact with the epitope. By exploiting private and public datasets and comparing TCRpcDist with competing approaches, it is demonstrated that TCRpcDist can accurately identify groups of TCRs that are likely to bind the same epitopes. Importantly, the ability of TCRpcDist is experimentally validated to determine antigen specificities (neoantigens and tumor-associated antigens) of orphan tumor-infiltrating lymphocytes (TILs) in cancer patients. TCRpcDist is thus a promising approach to support TCR repertoire analysis and TCR deorphanization for individualized treatments including cancer immunotherapies.

Patient-derived xenografts and single-cell sequencing identifies three subtypes of tumor-reactive lymphocytes in uveal melanoma metastases

Uveal melanoma (UM) is a rare melanoma originating in the eye's uvea, with 50% of patients experiencing metastasis predominantly in the liver. In contrast to cutaneous melanoma, there is only a limited effectiveness of combined immune checkpoint therapies, and half of patients with uveal melanoma metastases succumb to disease within 2 years. This study aimed to provide a path toward enhancing immunotherapy efficacy by identifying and functionally validating tumor-reactive T cells in liver metastases of patients with UM. We employed single-cell RNA-seq of biopsies and tumor-infiltrating lymphocytes (TILs) to identify potential tumor-reactive T cells. Patient-derived xenograft (PDX) models of UM metastases were created from patients, and tumor sphere cultures were generated from these models for co-culture with autologous or MART1-specific HLA-matched allogenic TILs. Activated T cells were subjected to TCR-seq, and the TCRs were matched to those found in single-cell sequencing data from biopsies, expanded TILs, and in livers or spleens of PDX models injected with TILs. Our findings revealed that tumor-reactive T cells resided not only among activated and exhausted subsets of T cells, but also in a subset of cytotoxic effector cells. In conclusion, combining single-cell sequencing and functional analysis provides valuable insights into which T cells in UM may be useful for cell therapy amplification and marker selection.

Exploring the potential of structure-based deep learning approaches for T cell receptor design

Deep learning methods, trained on the increasing set of available protein 3D structures and sequences, have substantially impacted the protein modeling and design field. These advancements have facilitated the creation of novel proteins, or the optimization of existing ones designed for specific functions, such as binding a target protein. Despite the demonstrated potential of such approaches in designing general protein binders, their application in designing immunotherapeutics remains relatively underexplored. A relevant application is the design of T cell receptors (TCRs). Given the crucial role of T cells in mediating immune responses, redirecting these cells to tumor or infected target cells through the engineering of TCRs has shown promising results in treating diseases, especially cancer. However, the computational design of TCR interactions presents challenges for current physics-based methods, particularly due to the unique natural characteristics of these interfaces, such as low affinity and cross-reactivity. For this reason, in this study, we explored the potential of two structure-based deep learning protein design methods, ProteinMPNN and ESM-IF1, in designing fixed-backbone TCRs for binding target antigenic peptides presented by the MHC through different design scenarios. To evaluate TCR designs, we employed a comprehensive set of sequence- and structure-based metrics, highlighting the benefits of these methods in comparison to classical physics-based design methods and identifying deficiencies for improvement.

B-cell-directed CAR T-cell therapy activates CD8+ cytotoxic CARneg bystander T cells in patients and nonhuman primates

ABSTRACT

Chimeric antigen receptor (CAR) T cells hold promise as a therapy for B-cell-derived malignancies, and despite their impressive initial response rates, a significant proportion of patients ultimately experience relapse. Although recent studies have explored the mechanisms of in vivo CAR T-cell function, little is understood about the activation of surrounding CARneg bystander T cells and their potential to enhance tumor responses. We performed single-cell RNA sequencing on nonhuman primate (NHP) and patient-derived T cells to identify the phenotypic and transcriptomic hallmarks of bystander activation of CARneg T cells following B-cell-targeted CAR T-cell therapy. Using a highly translatable CD20 CAR NHP model, we observed a distinct population of activated CD8+ CARneg T cells emerging during CAR T-cell expansion. These bystander CD8+ CARneg T cells exhibited a unique transcriptional signature with upregulation of natural killer-cell markers (KIR3DL2, CD160, and KLRD1), chemokines, and chemokine receptors (CCL5, XCL1, and CCR9), and downregulation of naïve T-cell-associated genes (SELL and CD28). A transcriptionally similar population was identified in patients after a tisagenlecleucel infusion. Mechanistic studies revealed that interleukin-2 (IL-2) and IL-15 exposure induced bystander-like CD8+ T cells in a dose-dependent manner. In vitro activated and patient-derived T cells with a bystander phenotype efficiently killed leukemic cells through a T-cell receptor-independent mechanism. Collectively, to our knowledge, these data provide the first comprehensive identification and profiling of CARneg bystander CD8+ T cells following B-cell-targeting CAR T-cell therapy and suggest a novel mechanism through which CAR T-cell infusion might trigger enhanced antileukemic responses. Patient samples were obtained from the trial #NCT03369353, registered at www.ClinicalTrials.gov.

A structure-guided approach to predict MHC-I restriction of T cell receptors for public antigens

Peptides presented by class I major histocompatibility complex (MHC-I) proteins provide biomarkers for therapeutic targeting using T cell receptors (TCRs), TCR-mimicking antibodies (TMAs), or other engineered protein binders. Despite the extreme sequence diversity of the Human Leucocyte Antigen (HLA, the human MHC), a given TCR or TMA is restricted to recognize epitopic peptides in the context of a limited set of different HLA allotypes. Here, guided by our analysis of 96 TCR:pHLA complex structures in the Protein Data Bank (PDB), we identify TCR contact residues and classify 148 common HLA allotypes into T-cell cross-reactivity groups (T-CREGs) on the basis of their interaction surface features. Insights from our work have actionable value for resolving MHC-I restriction of TCRs, guiding therapeutic expansion of existing therapies, and informing the selection of peptide targets for forthcoming immunotherapy modalities.

Structural characterization and AlphaFold modeling of human T cell receptor recognition of NRAS cancer neoantigens

T cell receptors (TCRs) that recognize cancer neoantigens are important for anti-cancer immune responses and immunotherapy. Understanding the structural basis of TCR recognition of neoantigens provides insights into their exquisite specificity and can enable design of optimized TCRs. We determined crystal structures of a human TCR in complex with NRAS Q61K and Q61R neoantigen peptides and HLA-A1 MHC, revealing the molecular underpinnings for dual recognition and specificity versus wild-type NRAS peptide. We then used multiple versions of AlphaFold to model the corresponding complex structures, given the challenge of immune recognition for such methods. Interestingly, one implementation of AlphaFold2 (TCRmodel2) was able to generate accurate models of the complexes, while AlphaFold3 also showed strong performance, although success was lower for other complexes. This study provides insights into TCR recognition of a shared cancer neoantigen, as well as the utility and practical considerations for using AlphaFold to model TCR-peptide-MHC complexes.

RACER-m leverages structural features for sparse T cell specificity prediction

Reliable prediction of T cell specificity against antigenic signatures is a formidable task, complicated by the immense diversity of T cell receptor and antigen sequence space and the resulting limited availability of training sets for inferential models. Recent modeling efforts have demonstrated the advantage of incorporating structural information to overcome the need for extensive training sequence data, yet disentangling the heterogeneous TCR-antigen interface to accurately predict MHC-allele-restricted TCR-peptide interactions has remained challenging. Here, we present RACER-m, a coarse-grained structural model leveraging key biophysical information from the diversity of publicly available TCR-antigen crystal structures. Explicit inclusion of structural content substantially reduces the required number of training examples and maintains reliable predictions of TCR-recognition specificity and sensitivity across diverse biological contexts. Our model capably identifies biophysically meaningful point-mutant peptides that affect binding affinity, distinguishing its ability in predicting TCR specificity of point-mutants from alternative sequence-based methods. Its application is broadly applicable to studies involving both closely related and structurally diverse TCR-peptide pairs.

Structure-Based Rational and General Strategy for Stabilizing Single-Chain T-Cell Receptors to Enhance Affinity

The T-cell receptor (TCR) is a crucial molecule in cellular immunity. The single-chain T-cell receptor (scTCR) is a potential format in TCR therapeutics because it eliminates the possibility of αβ-TCR mispairing. However, its poor stability and solubility impede the in vitro study and manufacturing of therapeutic applications. In this study, some conserved structural motifs are identified in variable domains regardless of germlines and species. Theoretical analysis helps to identify those unfavored factors and leads to a general strategy for stabilizing scTCRs by substituting residues at exact IMGT positions with beneficial propensities on the consensus sequence of germlines. Several representative scTCRs are displayed to achieve stability optimization and retain comparable binding affinities with the corresponding αβ-TCRs in the range of μM to pM. These results demonstrate that our strategies for scTCR engineering are capable of providing the affinity-enhanced and specificity-retained format, which are of great value in facilitating the development of TCR-related therapeutics.

DIRMC: a database of immunotherapy-related molecular characteristics

Cancer immunotherapy has brought about a revolutionary breakthrough in the field of cancer treatment. Immunotherapy has changed the treatment landscape for a variety of solid and hematologic malignancies. To assist researchers in efficiently uncovering valuable information related to cancer immunotherapy, we have presented a manually curated comprehensive database called DIRMC, which focuses on molecular features involved in cancer immunotherapy. All the content was collected manually from published literature, authoritative clinical trial data submitted by clinicians, some databases for drug target prediction such as DrugBank, and some experimentally confirmed high-throughput data sets for the characterization of immune-related molecular interactions in cancer, such as a curated database of T-cell receptor sequences with known antigen specificity (VDJdb), a pathology-associated TCR database (McPAS-TCR) et al. By constructing a fully connected functional network, ranging from cancer-related gene mutations to target genes to translated target proteins to protein regions or sites that may specifically affect protein function, we aim to comprehensively characterize molecular features related to cancer immunotherapy. We have developed the scoring criteria to assess the reliability of each MHC-peptide-T-cell receptor (TCR) interaction item to provide a reference for users. The database provides a user-friendly interface to browse and retrieve data by genes, target proteins, diseases and more. DIRMC also provides a download and submission page for researchers to access data of interest for further investigation or submit new interactions related to cancer immunotherapy targets. Furthermore, DIRMC provides a graphical interface to help users predict the binding affinity between their own peptide of interest and MHC or TCR. This database will provide researchers with a one-stop resource to understand cancer immunotherapy-related targets as well as data on MHC-peptide-TCR interactions. It aims to offer reliable molecular characteristics support for both the analysis of the current status of cancer immunotherapy and the development of new immunotherapy. DIRMC is available at http://www.dirmc.tech/. Database URL: http://www.dirmc.tech/.

Can AlphaFold's breakthrough in protein structure help decode the fundamental principles of adaptive cellular immunity?

T cells are essential immune cells responsible for identifying and eliminating pathogens. Through interactions between their T-cell antigen receptors (TCRs) and antigens presented by major histocompatibility complex molecules (MHCs) or MHC-like molecules, T cells discriminate foreign and self peptides. Determining the fundamental principles that govern these interactions has important implications in numerous medical contexts. However, reconstructing a map between T cells and their antagonist antigens remains an open challenge for the field of immunology, and success of in silico reconstructions of this relationship has remained incremental. In this Perspective, we discuss the role that new state-of-the-art deep-learning models for predicting protein structure may play in resolving some of the unanswered questions the field faces linking TCR and peptide-MHC properties to T-cell specificity. We provide a comprehensive overview of structural databases and the evolution of predictive models, and highlight the breakthrough AlphaFold provided the field.

Exploring the Potential of Structure-Based Deep Learning Approaches for T cell Receptor Design

Deep learning methods, trained on the increasing set of available protein 3D structures and sequences, have substantially impacted the protein modeling and design field. These advancements have facilitated the creation of novel proteins, or the optimization of existing ones designed for specific functions, such as binding a target protein. Despite the demonstrated potential of such approaches in designing general protein binders, their application in designing immunotherapeutics remains relatively unexplored. A relevant application is the design of T cell receptors (TCRs). Given the crucial role of T cells in mediating immune responses, redirecting these cells to tumor or infected target cells through the engineering of TCRs has shown promising results in treating diseases, especially cancer. However, the computational design of TCR interactions presents challenges for current physics-based methods, particularly due to the unique natural characteristics of these interfaces, such as low affinity and cross-reactivity. For this reason, in this study, we explored the potential of two structure-based deep learning protein design methods, ProteinMPNN and ESM-IF, in designing fixed-backbone TCRs for binding target antigenic peptides presented by the MHC through different design scenarios. To evaluate TCR designs, we employed a comprehensive set of sequence- and structure-based metrics, highlighting the benefits of these methods in comparison to classical physics-based design methods and identifying deficiencies for improvement.

Enhanced T cell receptor specificity through framework engineering

Development of T cell receptors (TCRs) as immunotherapeutics is hindered by inherent TCR cross-reactivity. Engineering more specific TCRs has proven challenging, as unlike antibodies, improving TCR affinity does not usually improve specificity. Although various protein design approaches have been explored to surmount this, mutations in TCR binding interfaces risk broadening specificity or introducing new reactivities. Here we explored if TCR specificity could alternatively be tuned through framework mutations distant from the interface. Studying the 868 TCR specific for the HIV SL9 epitope presented by HLA-A2, we used deep mutational scanning to identify a framework mutation above the mobile CDR3β loop. This glycine to proline mutation had no discernable impact on binding affinity or functional avidity towards the SL9 epitope but weakened recognition of SL9 escape variants and led to fewer responses in a SL9-derived positional scanning library. In contrast, an interfacial mutation near the tip of CDR3α that also did not impact affinity or functional avidity towards SL9 weakened specificity. Simulations indicated that the specificity-enhancing mutation functions by reducing the range of loop motions, limiting the ability of the TCR to adjust to different ligands. Although our results are likely to be TCR dependent, using framework engineering to control TCR loop motions may be a viable strategy for improving the specificity of TCR-based immunotherapies.

Structural basis for self-discrimination by neoantigen-specific TCRs

T cell receptors (TCR) are pivotal in mediating tumour cell cytolysis via recognition of mutation-derived tumour neoantigens (neoAgs) presented by major histocompatibility class-I (MHC-I). Understanding the factors governing the emergence of neoAg from somatic mutations is a major focus of current research. However, the structural and cellular determinants controlling TCR recognition of neoAgs remain poorly understood. This study describes the multi-level analysis of a model neoAg from the B16F10 murine melanoma, H2-Db/Hsf2 p.K72N68-76, as well as its cognate TCR 47BE7. Through cellular, molecular and structural studies we demonstrate that the p.K72N mutation enhances H2-Db binding, thereby improving cell surface presentation and stabilizing the TCR 47BE7 epitope. Furthermore, TCR 47BE7 exhibited high functional avidity and selectivity, attributable to a broad, stringent, binding interface enabling recognition of native B16F10 despite low antigen density. Our findings provide insight into the generation of anchor-residue modified neoAg, and emphasize the value of molecular and structural investigations of neoAg in diverse MHC-I contexts for advancing the understanding of neoAg immunogenicity.

A significance score for protein-protein interaction models through random docking

Comparing accuracies of structural protein-protein interaction (PPI) models for different complexes on an absolute scale is a challenge, requiring normalization of scores across structures of different sizes and shapes. To help address this challenge, we have developed a statistical significance metric for docking models, called random-docking (RD) p-value. This score evaluates a PPI model based on how likely a random docking process is to produce a model of better or equal accuracy. The binding partners are randomly docked against each other a large number of times, and the probability of sampling a model of equal or greater accuracy from this reference distribution is the RD p-value. Using a subset of top predicted models from CAPRI (Critical Assessment of PRediction of Interactions) rounds over 2017-2020, we find that the ease of achieving a given root mean squared deviation or DOCKQ score varies considerably by target; achieving the same relative metric can be thousands of times easier for one complex compared to another. In contrast, RD p-values inherently normalize scores for models of different complexes, making them globally comparable. Furthermore, one can calculate RD p-values after generating a reference distribution that accounts for prior information about the interface geometry, such as residues involved in binding, by giving the random-docking process access the same information. Thus, one can decouple improvements in prediction accuracy that arise solely from basic modeling constraints from those due to the rest of the method. We provide efficient code for computing RD p-values at https://github.com/Grigoryanlab/RDP.

Attention network for predicting T-cell receptor-peptide binding can associate attention with interpretable protein structural properties

Understanding how a T-cell receptor (TCR) recognizes its specific ligand peptide is crucial for gaining an insight into biological functions and disease mechanisms. Despite its importance, experimentally determining TCR-peptide-major histocompatibility complex (TCR-pMHC) interactions is expensive and time-consuming. To address this challenge, computational methods have been proposed, but they are typically evaluated by internal retrospective validation only, and few researchers have incorporated and tested an attention layer from language models into structural information. Therefore, in this study, we developed a machine learning model based on a modified version of Transformer, a source-target attention neural network, to predict the TCR-pMHC interaction solely from the amino acid sequences of the TCR complementarity-determining region (CDR) 3 and the peptide. This model achieved competitive performance on a benchmark dataset of the TCR-pMHC interaction, as well as on a truly new external dataset. Additionally, by analyzing the results of binding predictions, we associated the neural network weights with protein structural properties. By classifying the residues into large- and small-attention groups, we identified statistically significant properties associated with the largely attended residues such as hydrogen bonds within CDR3. The dataset that we created and the ability of our model to provide an interpretable prediction of TCR-peptide binding should increase our knowledge about molecular recognition and pave the way for designing new therapeutics.

pan-MHC and cross-Species Prediction of T Cell Receptor-Antigen Binding

Profiling the binding of T cell receptors (TCRs) of T cells to antigenic peptides presented by MHC proteins is one of the most important unsolved problems in modern immunology. Experimental methods to probe TCR-antigen interactions are slow, labor-intensive, costly, and yield moderate throughput. To address this problem, we developed pMTnet-omni, an Artificial Intelligence (AI) system based on hybrid protein sequence and structure information, to predict the pairing of TCRs of αβ T cells with peptide-MHC complexes (pMHCs). pMTnet-omni is capable of handling peptides presented by both class I and II pMHCs, and capable of handling both human and mouse TCR-pMHC pairs, through information sharing enabled this hybrid design. pMTnet-omni achieves a high overall Area Under the Curve of Receiver Operator Characteristics (AUROC) of 0.888, which surpasses competing tools by a large margin. We showed that pMTnet-omni can distinguish binding affinity of TCRs with similar sequences. Across a range of datasets from various biological contexts, pMTnet-omni characterized the longitudinal evolution and spatial heterogeneity of TCR-pMHC interactions and their functional impact. We successfully developed a biomarker based on pMTnet-omni for predicting immune-related adverse events of immune checkpoint inhibitor (ICI) treatment in a cohort of 57 ICI-treated patients. pMTnet-omni represents a major advance towards developing a clinically usable AI system for TCR-pMHC pairing prediction that can aid the design and implementation of TCR-based immunotherapeutics.

GGNpTCR: A Generative Graph Structure Neural Network for Predicting Immunogenic Peptides for T-cell Immune Response

Identifying the interactions between T-cell receptor (TCRs) and human antigens is a crucial step in developing new vaccines, diagnostics, and immunotherapy. Current methods primarily focus on learning binding patterns from known TCR binding repertoires by using sequence information alone without considering the binding specificity of new antigens or exogenous peptides that have not appeared in the training set. Furthermore, the spatial structure of antigens plays a critical role in immune studies and immunotherapy, which should be addressed properly in the identification of interacting TCR-antigen pairs. In this study, we introduced a novel deep learning framework based on generative graph structures, GGNpTCR, for predicting interactions between TCR and peptides from sequence information. Results of real data analysis indicate that our model achieved excellent prediction for new antigens unseen in the training data set, making significant improvements compared to existing methods. We also applied the model to a large COVID-19 data set with no antigens in the training data set, and the improvement was also significant. Furthermore, through incorporation of additional supervised mechanisms, GGNpTCR demonstrated the ability to precisely forecast the locations of peptide-TCR interactions within 3D configurations. This enhancement substantially improved the model's interpretability. In summary, based on the performance on multiple data sets, GGNpTCR has made significant progress in terms of performance, universality, and interpretability.

Structural principles of peptide-centric chimeric antigen receptor recognition guide therapeutic expansion

Peptide-centric chimeric antigen receptors (PC-CARs) recognize oncoprotein epitopes displayed by cell-surface human leukocyte antigens (HLAs) and offer a promising strategy for targeted cancer therapy. We have previously developed a PC-CAR targeting a neuroblastoma-associated PHOX2B peptide, leading to robust tumor cell lysis restricted by two common HLA allotypes. Here, we determine the 2.1-angstrom crystal structure of the PC-CAR-PHOX2B-HLA-A*24:02-β2m complex, which reveals the basis for antigen-specific recognition through interactions with CAR complementarity-determining regions (CDRs). This PC-CAR adopts a diagonal docking mode, where interactions with both conserved and polymorphic HLA framework residues permit recognition of multiple HLA allotypes from the A9 serological cross-reactive group, covering a combined global population frequency of up to 46.7%. Biochemical binding assays, molecular dynamics simulations, and structural and functional analyses demonstrate that high-affinity PC-CAR recognition of cross-reactive pHLAs necessitates the presentation of a specific peptide backbone, where subtle structural adaptations of the peptide are critical for high-affinity complex formation, and CAR T cell killing. Our results provide a molecular blueprint for engineering CARs with optimal recognition of tumor-associated antigens in the context of different HLAs, while minimizing cross-reactivity with self-epitopes.

Globally shared TCR repertoires within the tumor-infiltrating lymphocytes of patients with metastatic gynecologic cancer

Gynecologic cancer, including ovarian cancer and endometrial cancer, is characterized by morphological and molecular heterogeneity. Germline and somatic testing are available for patients to screen for pathogenic variants in genes such as BRCA1/2. Tissue expression levels of immunogenomic markers such as PD-L1 are also being used in clinical research. The basic therapeutic approach to gynecologic cancer combines surgery with chemotherapy. Immunotherapy, while not yet a mainstream treatment for gynecologic cancers, is advancing, with Dostarlimab recently receiving approval as a treatment for endometrial cancer. The goal remains to harness stimulated immune cells in the bloodstream to eradicate multiple metastases, a feat currently deemed challenging in a typical clinical setting. For the discovery of novel immunotherapy-based tumor targets, tumor-infiltrating lymphocytes (TILs) give a key insight on tumor-related immune activities by providing T cell receptor (TCR) sequences. Understanding the TCR repertoires of TILs in metastatic tissues and the circulation is important from an immunotherapy standpoint, as a subset of T cells in the blood have the potential to help kill tumor cells. To explore the relationship between distant tissue biopsy regions and blood circulation, we investigated the TCR beta chain (TCRβ) in bulk tumor and matched blood samples from 39 patients with gynecologic cancer. We found that the TCR clones of TILs at different tumor sites were globally shared within patients and had high overlap with the TCR clones in peripheral blood.

SARS-CoV-2 infection establishes a stable and age-independent CD8+ T cell response against a dominant nucleocapsid epitope using restricted T cell receptors

The resolution of SARS-CoV-2 replication hinges on cell-mediated immunity, wherein CD8+ T cells play a vital role. Nonetheless, the characterization of the specificity and TCR composition of CD8+ T cells targeting non-spike protein of SARS-CoV-2 before and after infection remains incomplete. Here, we analyzed CD8+ T cells recognizing six epitopes from the SARS-CoV-2 nucleocapsid (N) protein and found that SARS-CoV-2 infection slightly increased the frequencies of N-recognizing CD8+ T cells but significantly enhanced activation-induced proliferation compared to that of the uninfected donors. The frequencies of N-specific CD8+ T cells and their proliferative response to stimulation did not decrease over one year. We identified the N222-230 peptide (LLLDRLNQL, referred to as LLL thereafter) as a dominant epitope that elicited the greatest proliferative response from both convalescent and uninfected donors. Single-cell sequencing of T cell receptors (TCR) from LLL-specific CD8+ T cells revealed highly restricted Vα gene usage (TRAV12-2) with limited CDR3α motifs, supported by structural characterization of the TCR-LLL-HLA-A2 complex. Lastly, transcriptome analysis of LLL-specific CD8+ T cells from donors who had expansion (expanders) or no expansion (non-expanders) after in vitro stimulation identified increased chromatin modification and innate immune functions of CD8+ T cells in non-expanders. These results suggests that SARS-CoV-2 infection induces LLL-specific CD8+ T cell responses with a restricted TCR repertoire.

Conserved biophysical compatibility among the highly variable germline-encoded regions shapes TCR-MHC interactions

T cells are critically important components of the adaptive immune system primarily responsible for identifying and responding to pathogenic challenges. This recognition of pathogens is driven by the interaction between membrane-bound T cell receptors (TCRs) and antigenic peptides presented on major histocompatibility complex (MHC) molecules. The formation of the TCR-peptide-MHC complex (TCR-pMHC) involves interactions among germline-encoded and hypervariable amino acids. Germline-encoded and hypervariable regions can form contacts critical for complex formation, but only interactions between germline-encoded contacts are likely to be shared across many of all the possible productive TCR-pMHC complexes. Despite this, experimental investigation of these interactions have focused on only a small fraction of the possible interaction space. To address this, we analyzed every possible germline-encoded TCR-MHC contact in humans, thereby generating the first comprehensive characterization of these largely antigen-independent interactions. Our computational analysis suggests that germline-encoded TCR-MHC interactions that are conserved at the sequence level are rare due to the high amino acid diversity of the TCR CDR1 and CDR2 loops, and that such conservation is unlikely to dominate the dynamic protein-protein binding interface. Instead, we propose that binding properties such as the docking orientation are defined by regions of biophysical compatibility between these loops and the MHC surface.

TAPIR: a T-cell receptor language model for predicting rare and novel targets

T-cell receptors (TCRs) are involved in most human diseases, but linking their sequences with their targets remains an unsolved grand challenge in the field. In this study, we present TAPIR (T-cell receptor and Peptide Interaction Recognizer), a T-cell receptor (TCR) language model that predicts TCR-target interactions, with a focus on novel and rare targets. TAPIR employs deep convolutional neural network (CNN) encoders to process TCR and target sequences across flexible representations (e.g., beta-chain only, unknown MHC allele, etc.) and learns patterns of interactivity via several training tasks. This flexibility allows TAPIR to train on more than 50k either paired (alpha and beta chain) or unpaired TCRs (just alpha or beta chain) from public and proprietary databases against 1933 unique targets. TAPIR demonstrates state-of-the-art performance when predicting TCR interactivity against common benchmark targets and is the first method to demonstrate strong performance when predicting TCR interactivity against novel targets, where no examples are provided in training. TAPIR is also capable of predicting TCR interaction against MHC alleles in the absence of target information. Leveraging these capabilities, we apply TAPIR to cancer patient TCR repertoires and identify and validate a novel and potent anti-cancer T-cell receptor against a shared cancer neoantigen target (PIK3CA H1047L). We further show how TAPIR, when extended with a generative neural network, is capable of directly designing T-cell receptor sequences that interact with a target of interest.

Hydrophobic interactions dominate the recognition of a KRAS G12V neoantigen

Specificity remains a major challenge to current therapeutic strategies for cancer. Mutation associated neoantigens (MANAs) are products of genetic alterations, making them highly specific therapeutic targets. MANAs are HLA-presented (pHLA) peptides derived from intracellular mutant proteins that are otherwise inaccessible to antibody-based therapeutics. Here, we describe the cryo-EM structure of an antibody-MANA pHLA complex. Specifically, we determine a TCR mimic (TCRm) antibody bound to its MANA target, the KRASG12V peptide presented by HLA-A*03:01. Hydrophobic residues appear to account for the specificity of the mutant G12V residue. We also determine the structure of the wild-type G12 peptide bound to HLA-A*03:01, using X-ray crystallography. Based on these structures, we perform screens to validate the key residues required for peptide specificity. These experiments led us to a model for discrimination between the mutant and the wild-type peptides presented on HLA-A*03:01 based exclusively on hydrophobic interactions.

Large-scale template-based structural modeling of T-cell receptors with known antigen specificity reveals complementarity features

Introduction

T-cell receptor (TCR) recognition of foreign peptides presented by the major histocompatibility complex (MHC) initiates the adaptive immune response against pathogens. While a large number of TCR sequences specific to different antigenic peptides are known to date, the structural data describing the conformation and contacting residues for TCR-peptide-MHC complexes is relatively limited. In the present study we aim to extend and analyze the set of available structures by performing highly accurate template-based modeling of these complexes using TCR sequences with known specificity.

Methods

Identification of CDR3 sequences and their further clustering, based on available spatial structures, V- and J-genes of corresponding T-cell receptors, and epitopes, was performed using the VDJdb database. Modeling of the selected CDR3 loops was conducted using a stepwise introduction of single amino acid substitutions to the template PDB structures, followed by optimization of the TCR-peptide-MHC contacting interface using the Rosetta package applications. Statistical analysis and recursive feature elimination procedures were carried out on computed energy values and properties of contacting amino acid residues between CDR3 loops and peptides, using R.

Results

Using the set of 29 complex templates (including a template with SARS-CoV-2 antigen) and 732 specificity records, we built a database of 1585 model structures carrying substitutions in either TCRα or TCRβ chains with some models representing the result of different mutation pathways for the same final structure. This database allowed us to analyze features of amino acid contacts in TCR - peptide interfaces that govern antigen recognition preferences and interpret these interactions in terms of physicochemical properties of interacting residues.

Conclusion

Our results provide a methodology for creating high-quality TCR-peptide-MHC models for antigens of interest that can be utilized to predict TCR specificity.

Dynamics and survival associations of T cell receptor clusters in patients with pleural mesothelioma treated with immunotherapy

BACKGROUND

Immune checkpoint inhibitors (ICIs) are now a first-line treatment option for patients with pleural mesothelioma with the recent approval of ipilimumab and nivolumab. Mesothelioma has a low tumor mutation burden and no robust predictors of survival with ICI. Since ICIs enable adaptive antitumor immune responses, we investigated T-cell receptor (TCR) associations with survival in participants from two clinical trials treated with ICI.

METHODS

We included patients with pleural mesothelioma who were treated with nivolumab (NivoMes, NCT02497508) or nivolumab and ipilimumab (INITIATE, NCT03048474) after first-line therapy. TCR sequencing was performed with the ImmunoSEQ assay in 49 and 39 pretreatment and post-treatment patient peripheral blood mononuclear cell (PBMC) samples. These data were integrated with TCR sequences found in bulk RNAseq data by TRUST4 program in 45 and 35 pretreatment and post-treatment tumor biopsy samples and TCR sequences from over 600 healthy controls. The TCR sequences were clustered into groups of shared antigen specificity using GIANA. Associations of TCR clusters with overall survival were determined by cox proportional hazard analysis.

RESULTS

We identified 4.2 million and 12 thousand complementarity-determining region 3 (CDR3) sequences from PBMCs and tumors, respectively, in patients treated with ICI. These CDR3 sequences were integrated with 2.1 million publically available CDR3 sequences from healthy controls and clustered. ICI-enhanced T-cell infiltration and expanded T cell diversity in tumors. Cases with TCR clones in the top tertile in the pretreatment tissue or in circulation had significantly better survival than the bottom two tertiles (p<0.04). Furthermore, a high number of shared TCR clones between pretreatment tissue and in circulation was associated with improved survival (p=0.01). To potentially select antitumor clusters, we filtered for clusters that were (1) not found in healthy controls, (2) recurrent in multiple patients with mesothelioma, and (3) more prevalent in post-treatment than pretreatment samples. The detection of two-specific TCR clusters provided significant survival benefit compared with detection of 1 cluster (HR<0.001, p=0.026) or the detection of no TCR clusters (HR=0.10, p=0.002). These two clusters were not found in bulk tissue RNA-seq data and have not been reported in public CDR3 databases.

CONCLUSIONS

We identified two unique TCR clusters that were associated with survival on treatment with ICI in patients with pleural mesothelioma. These clusters may enable approaches for antigen discovery and inform future targets for design of adoptive T cell therapies.

TCRmodel2: high-resolution modeling of T cell receptor recognition using deep learning

The cellular immune system, which is a critical component of human immunity, uses T cell receptors (TCRs) to recognize antigenic proteins in the form of peptides presented by major histocompatibility complex (MHC) proteins. Accurate definition of the structural basis of TCRs and their engagement of peptide-MHCs can provide major insights into normal and aberrant immunity, and can help guide the design of vaccines and immunotherapeutics. Given the limited amount of experimentally determined TCR-peptide-MHC structures and the vast amount of TCRs within each individual as well as antigenic targets, accurate computational modeling approaches are needed. Here, we report a major update to our web server, TCRmodel, which was originally developed to model unbound TCRs from sequence, to now model TCR-peptide-MHC complexes from sequence, utilizing several adaptations of AlphaFold. This method, named TCRmodel2, allows users to submit sequences through an easy-to-use interface and shows similar or greater accuracy than AlphaFold and other methods to model TCR-peptide-MHC complexes based on benchmarking. It can generate models of complexes in 15 minutes, and output models are provided with confidence scores and an integrated molecular viewer. TCRmodel2 is available at https://tcrmodel.ibbr.umd.edu.

Structural analysis of cancer-relevant TCR-CD3 and peptide-MHC complexes by cryoEM

The recognition of antigenic peptide-MHC (pMHC) molecules by T-cell receptors (TCR) initiates the T-cell mediated immune response. Structural characterization is key for understanding the specificity of TCR-pMHC interactions and informing the development of therapeutics. Despite the rapid rise of single particle cryoelectron microscopy (cryoEM), x-ray crystallography has remained the preferred method for structure determination of TCR-pMHC complexes. Here, we report cryoEM structures of two distinct full-length α/β TCR-CD3 complexes bound to their pMHC ligand, the cancer-testis antigen HLA-A2/MAGEA4 (230-239). We also determined cryoEM structures of pMHCs containing MAGEA4 (230-239) peptide and the closely related MAGEA8 (232-241) peptide in the absence of TCR, which provided a structural explanation for the MAGEA4 preference displayed by the TCRs. These findings provide insights into the TCR recognition of a clinically relevant cancer antigen and demonstrate the utility of cryoEM for high-resolution structural analysis of TCR-pMHC interactions.

Structural insights into protection against a SARS-CoV-2 spike variant by T cell receptor diversity

T cells play a crucial role in combatting SARS-CoV-2 and forming long-term memory responses to this coronavirus. The emergence of SARS-CoV-2 variants that can evade T cell immunity has raised concerns about vaccine efficacy and the risk of reinfection. Some SARS-CoV-2 T cell epitopes elicit clonally restricted CD8+ T cell responses characterized by T cell receptors (TCRs) that lack structural diversity. Mutations in such epitopes can lead to loss of recognition by most T cells specific for that epitope, facilitating viral escape. Here, we studied an HLA-A2-restricted spike protein epitope (RLQ) that elicits CD8+ T cell responses in COVID-19 convalescent patients characterized by highly diverse TCRs. We previously reported the structure of an RLQ-specific TCR (RLQ3) with greatly reduced recognition of the most common natural variant of the RLQ epitope (T1006I). Opposite to RLQ3, TCR RLQ7 recognizes T1006I with even higher functional avidity than the WT epitope. To explain the ability of RLQ7, but not RLQ3, to tolerate the T1006I mutation, we determined structures of RLQ7 bound to RLQ-HLA-A2 and T1006I-HLA-A2. These complexes show that there are multiple structural solutions to recognizing RLQ and thereby generating a clonally diverse T cell response to this epitope that assures protection against viral escape and T cell clonal loss.

UcTCRdb: An unconventional T cell receptor sequence database with online analysis functions

Unlike conventional major histocompatibility complex (MHC) class I and II molecules reactive T cells, the unconventional T cell subpopulations recognize various non-polymorphic antigen-presenting molecules and are typically characterized by simplified patterns of T cell receptors (TCRs), rapid effector responses and ‘public’ antigen specificities. Dissecting the recognition patterns of the non-MHC antigens by unconventional TCRs can help us further our understanding of the unconventional T cell immunity. The small size and irregularities of the released unconventional TCR sequences are far from high-quality to support systemic analysis of unconventional TCR repertoire. Here we present UcTCRdb, a database that contains 669,900 unconventional TCRs collected from 34 corresponding studies in humans, mice, and cattle. In UcTCRdb, users can interactively browse TCR features of different unconventional T cell subsets in different species, search and download sequences under different conditions. Additionally, basic and advanced online TCR analysis tools have been integrated into the database, which will facilitate the study of unconventional TCR patterns for users with different backgrounds. UcTCRdb is freely available at http://uctcrdb.cn/.

Immune repertoire and responses to neoadjuvant TCHP therapy in HER2-positive breast cancer

Background

Despite the introduction of trastuzumab, pathologic complete response (pCR) is not attained in approximately 30-40% of Human epithelial growth factor receptor-2-positive breast cancer. Tumor-infiltrating lymphocytes (TIL) have been suggested as a predictive marker of treatment response, albeit not always effective. We investigated the relationship between trastuzumab, docetaxel, carboplatin, and pertuzumab (TCHP) treatment and immune repertoire as a treatment response predictor.

Design

In all, 35 cases were divided into two experimental groups: 10 and 25 cases in the preliminary and main experiments, respectively. In the preliminary experiment, the biopsy tissues before TCHP treatment and the surgical tissues after TCHP treatment were compared. In the main experiment, the biopsy tissues before TCHP treatment were compared according to the TCHP treatment response.

Methods

The T-cell repertoire for TRA, TRB, TRG, and TRD, and B-cell repertoire for immunoglobulin heavy, immunoglobulin kappa, and immunoglobulin lambda were evaluated. Whole transcriptome sequencing was also performed.

Results

In the preliminary experiment, the density and richness of the T-cell receptor (TCR) and B-cell receptor (BCR) repertoires decreased after treatment, regardless of TCHP response. In the main experiment, the Shannon's entropy index, density, and length of CDR3 of the TCR and BCR repertoires did not differ significantly in patients who did and did not achieve pCR. The pCR and non-pCR subgroups according to the level of TILs revealed that the non-pCR/lowTIL group had a higher proportion of low-frequency clones than the pCR/lowTIL group in TRA (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 63% versus 45.3%; <0.01%, 32.9% versus 51.8%, p < 0.001) and TRB (non-pCR/lowTIL versus pCR/lowTIL, 0.01-0.1%, 26.5% versus 14.7%; <0.01%, 72.0% versus 84.1%, p < 0.001).

Conclusions

The role of the diversity, richness, and density of the TCR and BCR repertoires as predictive markers for TCHP response was not identified. Compositions of low-frequency clones could be candidates for predictive factors of TCHP response; however, validation studies and further research are necessary.

Immune repertoire profiling for disease pathobiology

Lymphocytes consist of highly heterogeneous populations, each expressing a specific cell surface receptor corresponding to a particular antigen. Lymphocytes are both the cause and regulator of various diseases, including autoimmune/allergic diseases, lifestyle diseases, neurodegenerative diseases, and cancers. Recently, immune repertoire sequencing has attracted much attention because it helps obtain global profiles of the immune receptor sequences of infiltrating T and B cells in specimens. Immune repertoire sequencing not only helps deepen our understanding of the molecular mechanisms of immune-related pathology but also assists in discovering novel therapeutic modalities for diseases, thereby shedding colorful light on otherwise tiny monotonous cells when observed under a microscope. In this review article, we introduce and detail the background and methodology of immune repertoire sequencing and summarize recent scientific achievements in association with human diseases. Future perspectives on this genetic technique in the field of histopathological research will also be discussed.

Prediction of Peptide and TCR CDR3 Loops in Formation of Class I MHC-Peptide-TCR Complexes Using Molecular Models with Solvation

Formation of major histocompatibility (MHC)-peptide-T cell receptor (TCR) complexes is central to initiation of an adaptive immune response. These complexes form through initial stabilization of the MHC fold via binding of a short peptide, and subsequent interaction of the TCR to form a ternary complex, with contacts made predominantly through the complementarity-determining region (CDR) loops of the TCR. Stimulation of an immune response is central to cancer immunotherapy. This approach depends on identification of the appropriate combinations of MHC molecules, peptides, and TCRs to elicit an antitumor immune response. This prediction is a current challenge in computational biochemistry. In this chapter, we introduce a predictive method that involves generation of multiple peptides and TCR CDR 3 loop conformations, solvation of these conformers in the context of the MHC-peptide-TCR ternary complex, extraction of parameters from the generated complexes, and use of an AI model to evaluate the potential for the assembled ternary complex to support an immune response.

Deaza-modification of MR1 ligands modulates recognition by MR1-restricted T cells

MR1-restricted T (MR1T) cells recognize microbial small molecule metabolites presented on the MHC Class I-like molecule MR1 and have been implicated in early effector responses to microbial infection. As a result, there is considerable interest in identifying chemical properties of metabolite ligands that permit recognition by MR1T cells, for consideration in therapeutic or vaccine applications. Here, we made chemical modifications to known MR1 ligands to evaluate the effect on MR1T cell activation. Specifically, we modified 6,7-dimethyl-8-D-ribityllumazine (DMRL) to generate 6,7-dimethyl-8-D-ribityldeazalumazine (DZ), and then further derivatized DZ to determine the requirements for retaining MR1 surface stabilization and agonistic properties. Interestingly, the IFN-γ response toward DZ varied widely across a panel of T cell receptor (TCR)-diverse MR1T cell clones; while one clone was agnostic toward the modification, most displayed either an enhancement or depletion of IFN-γ production when compared with its response to DMRL. To gain insight into a putative mechanism behind this phenomenon, we used in silico molecular docking techniques for DMRL and its derivatives and performed molecular dynamics simulations of the complexes. In assessing the dynamics of each ligand in the MR1 pocket, we found that DMRL and DZ exhibit differential dynamics of both the ribityl moiety and the aromatic backbone, which may contribute to ligand recognition. Together, our results support an emerging hypothesis for flexibility in MR1:ligand-MR1T TCR interactions and enable further exploration of the relationship between MR1:ligand structures and MR1T cell recognition for downstream applications targeting MR1T cells.

A class-mismatched TCR bypasses MHC restriction via an unorthodox but fully functional binding geometry

MHC restriction, which describes the binding of TCRs from CD4⁺ T cells to class II MHC proteins and TCRs from CD8⁺ T cells to class I MHC proteins, is a hallmark of immunology. Seemingly rare TCRs that break this paradigm exist, but mechanistic insight into their behavior is lacking. TIL1383I is a prototypical class-mismatched TCR, cloned from a CD4⁺ T cell but recognizing the tyrosinase tumor antigen presented by the class I MHC HLA-A2 in a fully functional manner. Here we find that TIL1383I binds this class I target with a highly atypical geometry. Despite unorthodox binding, TCR signaling, antigen specificity, and the ability to use CD8 are maintained. Structurally, a key feature of TIL1383I is an exceptionally long CDR3β loop that mediates functions that are traditionally performed separately by hypervariable and germline loops in canonical TCR structures. Our findings thus expand the range of known TCR binding geometries compatible with normal function and specificity, provide insight into the determinants of MHC restriction, and may help guide TCR selection and engineering for immunotherapy.

Machine-Learning-Assisted Analysis of TCR Profiling Data Unveils Cross-Reactivity between SARS-CoV-2 and a Wide Spectrum of Pathogens and Other Diseases

During the last two years, the emergence of SARS-CoV-2 has led to millions of deaths worldwide, with a devastating socio-economic impact on a global scale. The scientific community's focus has recently shifted towards the association of the T cell immunological repertoire with COVID-19 progression and severity, by utilising T cell receptor sequencing (TCR-Seq) assays. The Multiplexed Identification of T cell Receptor Antigen (MIRA) dataset, which is a subset of the immunoACCESS study, provides thousands of TCRs that can specifically recognise SARS-CoV-2 epitopes. Our study proposes a novel Machine Learning (ML)-assisted approach for analysing TCR-Seq data from the antigens' point of view, with the ability to unveil key antigens that can accurately distinguish between MIRA COVID-19-convalescent and healthy individuals based on differences in the triggered immune response. Some SARS-CoV-2 antigens were found to exhibit equal levels of recognition by MIRA TCRs in both convalescent and healthy cohorts, leading to the assumption of putative cross-reactivity between SARS-CoV-2 and other infectious agents. This hypothesis was tested by combining MIRA with other public TCR profiling repositories that host assays and sequencing data concerning a plethora of pathogens. Our study provides evidence regarding putative cross-reactivity between SARS-CoV-2 and a wide spectrum of pathogens and diseases, with M. tuberculosis and Influenza virus exhibiting the highest levels of cross-reactivity. These results can potentially shift the emphasis of immunological studies towards an increased application of TCR profiling assays that have the potential to uncover key mechanisms of cell-mediated immune response against pathogens and diseases.

Benchmarking AlphaFold for protein complex modeling reveals accuracy determinants

High-resolution experimental structural determination of protein-protein interactions has led to valuable mechanistic insights, yet due to the massive number of interactions and experimental limitations there is a need for computational methods that can accurately model their structures. Here we explore the use of the recently developed deep learning method, AlphaFold, to predict structures of protein complexes from sequence. With a benchmark of 152 diverse heterodimeric protein complexes, multiple implementations and parameters of AlphaFold were tested for accuracy. Remarkably, many cases (43%) had near-native models (medium or high critical assessment of predicted interactions accuracy) generated as top-ranked predictions by AlphaFold, greatly surpassing the performance of unbound protein-protein docking (9% success rate for near-native top-ranked models), however AlphaFold modeling of antibody-antigen complexes within our set was unsuccessful. We identified sequence and structural features associated with lack of AlphaFold success, and we also investigated the impact of multiple sequence alignment input. Benchmarking of a multimer-optimized version of AlphaFold (AlphaFold-Multimer) with a set of recently released antibody-antigen structures confirmed a low rate of success for antibody-antigen complexes (11% success), and we found that T cell receptor-antigen complexes are likewise not accurately modeled by that algorithm, showing that adaptive immune recognition poses a challenge for the current AlphaFold algorithm and model. Overall, our study demonstrates that end-to-end deep learning can accurately model many transient protein complexes, and highlights areas of improvement for future developments to reliably model any protein-protein interaction of interest.

Proteome-Scale Screening to Identify High-Expression Signal Peptides with Minimal N-Terminus Biases via Yeast Display

Signal peptides are critical for the efficient expression and routing of extracellular and secreted proteins. Most protein production and screening technologies rely upon a relatively small set of signal peptides. Despite their central role in biotechnology, there are limited studies comprehensively examining the interplay between signal peptides and expressed protein sequences. Here, we describe a high-throughput method to screen novel signal peptides that maintain a high degree of surface expression across a range of protein scaffolds with highly variable N-termini. We find that the canonical signal peptide used in yeast surface display, derived from Aga2p, fails to achieve high surface expression for 42.5% of constructs containing diverse N-termini. To circumvent this, we have identified two novel signal peptides derived from endogenous yeast proteins, SRL1 and KISH, which are highly tolerant to diverse N-terminal sequences. This pipeline can be used to expand our understanding of signal peptide function, identify improved signal peptides for protein expression, and refine the computational tools used for signal peptide prediction.

Immunogenicity and therapeutic targeting of a public neoantigen derived from mutated PIK3CA

Public neoantigens (NeoAgs) represent an elite class of shared cancer-specific epitopes derived from recurrently mutated driver genes. Here we describe a high-throughput platform combining single-cell transcriptomic and T cell receptor (TCR) sequencing to establish whether mutant PIK3CA, among the most frequently genomically altered driver oncogenes, generates an immunogenic public NeoAg. Using this strategy, we developed a panel of TCRs that recognize an endogenously processed neopeptide encompassing a common PIK3CA hotspot mutation restricted by the prevalent human leukocyte antigen (HLA)-A*03:01 allele. Mechanistically, immunogenicity to this public NeoAg arises from enhanced neopeptide/HLA complex stability caused by a preferred HLA anchor substitution. Structural studies indicated that the HLA-bound neopeptide presents a comparatively 'featureless' surface dominated by the peptide's backbone. To bind this epitope with high specificity and affinity, we discovered that a lead TCR clinical candidate engages the neopeptide through an extended interface facilitated by an unusually long CDR3β loop. In patients with diverse malignancies, we observed NeoAg clonal conservation and spontaneous immunogenicity to the neoepitope. Finally, adoptive transfer of TCR-engineered T cells led to tumor regression in vivo in mice bearing PIK3CA-mutant tumors but not wild-type PIK3CA tumors. Together, these findings establish the immunogenicity and therapeutic potential of a mutant PIK3CA-derived public NeoAg.

T cell receptors (TCRs) employ diverse strategies to target a p53 cancer neoantigen

Adoptive cell therapy with tumor-specific T cells can mediate durable cancer regression. The prime target of tumor-specific T cells are neoantigens arising from mutations in self-proteins during malignant transformation. To understand T cell recognition of cancer neoantigens at the atomic level, we studied oligoclonal T cell receptors (TCRs) that recognize a neoepitope arising from a driver mutation in the p53 oncogene (p53R175H) presented by the major histocompatibility complex class I molecule HLA-A2. We previously reported the structures of three p53R175H-specific TCRs (38-10, 12-6, and 1a2) bound to p53R175H and HLA-A2. The structures showed that these TCRs discriminate between WT and mutant p53 by forming extensive interactions with the R175H mutation. Here, we report the structure of a fourth p53R175H-specific TCR (6-11) in complex with p53R175H and HLA-A2. In contrast to 38-10, 12-6, and 1a2, TCR 6-11 makes no direct contacts with the R175H mutation, yet is still able to distinguish mutant from WT p53. Structure-based in silico mutagenesis revealed that the 60-fold loss in 6-11 binding affinity for WT p53 compared to p53R175H is mainly due to the higher energetic cost of desolvating R175 in the WT p53 peptide during complex formation than H175 in the mutant. This indirect strategy for preferential neoantigen recognition by 6-11 is fundamentally different from the direct strategies employed by other TCRs and highlights the multiplicity of solutions to recognizing p53R175H with sufficient selectivity to mediate T cell killing of tumor but not normal cells.

Structural Prediction of Peptide-MHC Binding Modes

The immune system is constantly protecting its host from the invasion of pathogens and the development of cancer cells. The specific CD8⁺ T-cell immune response against virus-infected cells and tumor cells is based on the T-cell receptor recognition of antigenic peptides bound to class I major histocompatibility complexes (MHC) at the surface of antigen presenting cells. Consequently, the peptide binding specificities of the highly polymorphic MHC have important implications for the design of vaccines, for the treatment of autoimmune diseases, and for personalized cancer immunotherapy. Evidence-based machine-learning approaches have been successfully used for the prediction of peptide binders and are currently being developed for the prediction of peptide immunogenicity. However, understanding and modeling the structural details of peptide/MHC binding is crucial for a better understanding of the molecular mechanisms triggering the immunological processes, estimating peptide/MHC affinity using universal physics-based approaches, and driving the design of novel peptide ligands. Unfortunately, due to the large diversity of MHC allotypes and possible peptides, the growing number of 3D structures of peptide/MHC (pMHC) complexes in the Protein Data Bank only covers a small fraction of the possibilities. Consequently, there is a growing need for rapid and efficient approaches to predict 3D structures of pMHC complexes. Here, we review the key characteristics of the 3D structure of pMHC complexes before listing databases and other sources of information on pMHC structures and MHC specificities. Finally, we discuss some of the most prominent pMHC docking software.

OUP accepted manuscript

The study and manipulation of T cell receptors (TCRs) is central to multiple fields across basic and translational immunology research. Produced by V(D)J recombination, TCRs are often only recorded in the literature and data repositories as a combination of their V and J gene symbols, plus their hypervariable CDR3 amino acid sequence. However, numerous applications require full-length coding nucleotide sequences. Here we present Stitchr, a software tool developed to specifically address this limitation. Given minimal V/J/CDR3 information, Stitchr produces complete coding sequences representing a fully spliced TCR cDNA. Due to its modular design, Stitchr can be used for TCR engineering using either published germline or novel/modified variable and constant region sequences. Sequences produced by Stitchr were validated by synthesizing and transducing TCR sequences into Jurkat cells, recapitulating the expected antigen specificity of the parental TCR. Using a companion script, Thimble, we demonstrate that Stitchr can process a million TCRs in under ten minutes using a standard desktop personal computer. By systematizing the production and modification of TCR sequences, we propose that Stitchr will increase the speed, repeatability, and reproducibility of TCR research. Stitchr is available on GitHub.

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.

In-silico design of peptides for inhibition of HLA-A*03-KLIETYFSK complex as a new drug design for treatment of multiples sclerosis disease

Multiple sclerosis is recognized as a chronic inflammatory disease. Human leukocyte antigen (HLA) plays an important role in initiating adaptive immune responses. HLA class I is present in almost all nucleated cells and presents the cleaved endogenous peptide antigens to cytotoxic T cells. HLA-A*03 is one of the HLA class I alleles, which is reported as substantially related HLA to MS disease. In 2011, the structure of the HLA-A*03 in complex was identified with an immunodominant proteolipid protein (PLP) epitope (KLIETYFSK). This complex has been reported as an important autoantigen-presenting complex in MS pathogenesis. In this study, new peptides were designed to bind to this complex that may prevent specific pathogenic cytotoxic T cell binding to this autoantigen-presenting complex and CNS demyelination. Herein, 14 new helical peptides containing 19 amino acids were designed and their structures were predicted using the PEP-FOLD server. The binding of each designed peptide to the mentioned complex was then performed. A mutation approach was used by the BeAtMuSiC server to improve the binding affinity of the designed peptide. In each position, amino acid substitutions leading to an increase in the binding affinity of the peptide to the mentioned complex were determined. Finally, the resulting complexes were simulated for 40 ns using AMBER18 software. The results revealed that out of 14 designed peptides, "WRYWWKDWAKQFRQFYRWF" peptide exhibited the highest affinity for binding to the mentioned complex. This peptide can be considered as a potential drug to control multiple sclerosis disease in patients carrying the HLA-A*03 allele.