Giancarlo Croce

Interactions between T-Cell Receptors (TCRs) and antigenic peptides presented on Major Histocompatibility Complex (MHC) molecules are central to the immune recognition of infected and malignant cells. The complexity of the TCR sequence space, the flexibility of the TCR-epitope interface and the lack of a standardized framework to visualize TCR specificity have hindered a comprehensive understanding of the principles governing TCR-epitope recognition across a broad range of epitopes. Here, we introduce a fully interpretable probabilistic framework, termed TCR specificity profiles (TSPs), which captures fundamental properties distinguishing epitope-specific from baseline TCR repertoires. We demonstrate that TSPs unravel key determinants of TCR-epitope recognition specificity. By identifying and analyzing TCRs recognizing dozens of epitope variants, we show that TSPs accurately predict cross-reactivity and reveal how TCR specificity evolves with epitope sequence, binding mode and MHC restriction. TSPs further enable interpretation of machine learning tools and reveal how AlphaFold3 can be used to decipher key determinants of TCR-epitope recognition specificity.

T cells targeting epitopes in infectious diseases or cancer play a central role in spontaneous and therapy-induced immune responses. Epitope recognition is mediated by the binding of the T cell receptor (TCR), and TCRs recognizing clinically relevant epitopes are promising for T cell–based therapies. Starting from a TCR targeting the cancer-testis antigen NY-ESO-1_157–165 epitope, we built large phage display libraries of TCRs with randomized complementary determining region 3 of the β chain. The TCR libraries were panned against NY-ESO-1, which enabled us to collect thousands of epitope-specific TCR sequences. Leveraging these data, we trained a machine learning TCR-epitope interaction predictor and identified several epitope-specific TCRs from TCR repertoires. Cellular assays revealed that the predicted TCRs displayed activity toward NY-ESO-1 and no detectable cross-reactivity. Our work demonstrates how display technologies combined with TCR-epitope interaction predictors can effectively leverage large TCR repertoires for TCR discovery.

T cells have the ability to eliminate infected and cancer cells and play an essential role in cancer immunotherapy. T cell activation is elicited by the binding of the T cell receptor (TCR) to epitopes displayed on MHC molecules, and the TCR specificity is determined by the sequence of its α and β chains. Here, we collect and curate a dataset of 17,715 αβTCRs interacting with dozens of class I and class II epitopes. We use this curated data to develop MixTCRpred, an epitope-specific TCR-epitope interaction predictor. MixTCRpred accurately predicts TCRs recognizing several viral and cancer epitopes. MixTCRpred further provides a useful quality control tool for multiplexed single-cell TCR sequencing assays of epitope-specific T cells and pinpoints a substantial fraction of putative contaminants in public databases. Analysis of epitope-specific dual α T cells demonstrates that MixTCRpred can identify α chains mediating epitope recognition. Applying MixTCRpred to TCR repertoires from COVID-19 patients reveals enrichment of clonotypes predicted to bind an immunodominant SARS-CoV-2 epitope. Overall, MixTCRpred provides a robust tool to predict TCRs interacting with specific epitopes and interpret TCR-sequencing data from both bulk and epitope-specific T cells.

The emergence of new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major concern given their potential impact on the transmissibility and pathogenicity of the virus as well as the efficacy of therapeutic interventions. Here, we predict the mutability of all positions in SARS-CoV-2 protein domains to forecast the appearance of unseen variants. Using sequence data from other coronaviruses, preexisting to SARS-CoV-2, we build statistical models that not only capture amino acid conservation but also more complex patterns resulting from epistasis. We show that these models are notably superior to conservation profiles in estimating the already observable SARS-CoV-2 variability. In the receptor binding domain of the spike protein, we observe that the predicted mutability correlates well with experimental measures of protein stability and that both are reliable mutability predictors (receiver operating characteristic areas under the curve ∼0.8). Most interestingly, we observe an increasing agreement between our model and the observed variability as more data become available over time, proving the anticipatory capacity of our model. When combined with data concerning the immune response, our approach identifies positions where current variants of concern are highly overrepresented. These results could assist studies on viral evolution and future viral outbreaks and, in particular, guide the exploration and anticipation of potentially harmful future SARS-CoV-2 variants.