In silico specificity determination of neoantigen-reactive T-lymphocytes

Effective personalized immunotherapies of the future will need to capture not only the peculiarities of the patient’s tumor but also of his immune response to it. In this study, using results of in vitro high-throughput specificity assays, and combining comparative models of pMHCs and TCRs using molecular docking, we have constructed all-atom models for the putative complexes of all their possible pairwise TCR-pMHC combinations. For the models obtained we have calculated a dataset of physics-based scores and have trained binary classifiers that perform better compared to their solely sequence-based counterparts. These structure-based classifiers pinpoint the most prominent energetic terms and structural features characterizing the type of protein-protein interactions that underlies the immune recognition of tumors by T cells.