Accurate modeling of peptide-MHC structures with AlphaFold.

Major histocompatibility complex (MHC) proteins present peptides on the cell surface for T cell surveillance. Reliable in silico prediction of which peptides would be presented and which T cell receptors would recognize them is an important problem in structural immunology. Here, we introduce an AlphaFold-based pipeline for predicting the three-dimensional structures of peptide-MHC complexes for class I and class II MHC molecules. Our method demonstrates high accuracy, outperforming existing tools in class I modeling accuracy and class II peptide register prediction. We validate its performance and utility with new experimental data on a recently described cancer neoantigen/wild-type peptide pair and explore applications toward improving peptide-MHC binding prediction.

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.

Physicochemical Heuristics for Identifying High Fidelity, Near-Native Structural Models of Peptide/MHC Complexes.

There is long-standing interest in accurately modeling the structural features of peptides bound and presented by class I MHC proteins. This interest has grown with the advent of rapid genome sequencing and the prospect of personalized, peptide-based cancer vaccines, as well as the development of molecular and cellular therapeutics based on T cell receptor recognition of peptide-MHC. However, while the speed and accessibility of peptide-MHC modeling has improved substantially over the years, improvements in accuracy have been modest. Accuracy is crucial in peptide-MHC modeling, as T cell receptors are highly sensitive to peptide conformation and capturing fine details is therefore necessary for useful models. Studying nonameric peptides presented by the common class I MHC protein HLA-A*02:01, here we addressed a key question common to modern modeling efforts: from a set of models (or decoys) generated through conformational sampling, which is best? We found that the common strategy of decoy selection by lowest energy can lead to substantial errors in predicted structures. We therefore adopted a data-driven approach and trained functions capable of predicting near native decoys with exceptionally high accuracy. Although our implementation is limited to nonamer/HLA-A*02:01 complexes, our results serve as an important proof of concept from which improvements can be made and, given the significance of HLA-A*02:01 and its preference for nonameric peptides, should have immediate utility in select immunotherapeutic and other efforts for which structural information would be advantageous.