How peptide/MHC presence affects the dynamics of the LC13 T-cell receptor

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.

Engagement with the TCR induces plasticity in antigenic ligands bound to MHC class I and CD1 molecules

Complementarity-determining regions (CDRs) of αβ T-cell receptors (TCRs) sense peptide-bound MHC (pMHC) complexes via chemical interactions, thereby mediating antigen specificity and MHC restriction. Flexible finger-like movement of CDR loops contributes to the establishment of optimal interactions with pMHCs. In contrast, peptide ligands captured in MHC molecules are considered more static because of the rigid hydrogen-bond network that stabilizes peptide ligands in the antigen-binding groove of MHC molecules. An array of crystal structures delineating pMHC complexes in TCR-docked and TCR-undocked forms is now available, which enables us to assess TCR engagement-induced conformational changes in peptide ligands. In this short review, we overview conformational changes in MHC class I-bound peptide ligands upon TCR docking, followed by those for CD1-bound glycolipid ligands. Finally, we analyze the co-crystal structure of the TCR:lipopeptide-bound MHC class I complex that we recently reported. We argue that TCR engagement-induced conformational changes markedly occur in lipopeptide ligands, which are essential for exposure of a primary T-cell epitope to TCRs. These conformational changes are affected by amino acid residues, such as glycine, that do not interact directly with TCRs. Thus, ligand recognition by specific TCRs involves not only T-cell epitopes but also non-epitopic amino acid residues. In light of their critical function, we propose to refer to these residues as non-epitopic residues affecting ligand plasticity and antigenicity (NR-PA).

Can Glycosylation Mask the Detection of MHC Expressing p53 Peptides by T Cell Receptors?

Proteins of the major histocompatibility complex (MHC) class I, or human leukocyte antigen (HLA) in humans interact with endogenous peptides and present them to T cell receptors (TCR), which in turn tune the immune system to recognize and discriminate between self and foreign (non-self) peptides. Of especial importance are peptides derived from tumor-associated antigens. T cells recognizing these peptides are found in cancer patients, but not in cancer-free individuals. What stimulates this recognition, which is vital for the success of checkpoint based therapy? A peptide derived from the protein p53 (residues 161–169 or p161) was reported to show this behavior. T cells recognizing this unmodified peptide could be further stimulated in vitro to create effective cancer killing CTLs (cytotoxic T lymphocytes). We hypothesize that the underlying difference may arise from post-translational glycosylation of p161 in normal individuals, likely masking it against recognition by TCR. Defects in glycosylation in cancer cells may allow the presentation of the native peptide. We investigate the structural consequences of such peptide glycosylation by investigating the associated structural dynamics.

Unraveling the Allosteric Communication Mechanisms in T-Cell Receptor-Peptide-Loaded Major Histocompatibility Complex Dynamics Using Molecular Dynamics Simulations: An Approach Based on Dynamic Cross Correlation Maps and Residue Interaction Energy Calculations

Antigen presentation by major histocompatibility complex (MHC) proteins to T-cell receptors (TCRs) plays a crucial role in triggering the adaptive immune response. Most of our knowledge on TCR-peptide-loaded major histocompatibility complex (pMHC) interaction stemmed from experiments yielding static structures, yet the dynamic aspects of this molecular interaction are equally important to understand the underlying molecular mechanisms and to develop treatment strategies against diseases such as cancer and autoimmune diseases. To this end, computational biophysics studies including all-atom molecular dynamics simulations have provided useful insights; however, we still lack a basic understanding of an overall allosteric mechanism that results in conformational changes in the TCR and subsequent T-cell activation. Previous hydrogen-deuterium exchange and nuclear magnetic resonance studies provided clues regarding these molecular mechanisms, including global rigidification and allosteric effects on the constant domain of TCRs away from the pMHC interaction site. Here, we show that molecular dynamics simulations can be used to identify how this overall rigidification may be related to the allosteric communication within TCRs upon pMHC interaction via essential dynamics and nonbonded residue-residue interaction energy analyses. The residues taking part in the rigidification effect are highlighted with an intricate analysis on residue interaction changes, which lead to a detailed outline of the complex formation event. Our results indicate that residues of the Cβ domain of TCRs show significant differences in their nonbonded interactions upon complex formation. Moreover, the dynamic cross correlations between these residues are also increased, in line with their nonbonded interaction energy changes. Altogether, our approach may be valuable for elucidating intramolecular allosteric changes in the TCR structure upon pMHC interaction in molecular dynamics simulations.

Characterization of amino acid residues of T-cell receptors interacting with HLA-A*02-restricted antigen peptides

Background

The present study aimed to explore residues’ properties interacting with HLA-A*02-restricted peptides on T-cell receptors (TCRs) and their effects on bond types of interaction and binding free energy.

Methods

We searched the crystal structures of HLA-A*02-restricted peptide-TCR complexes from the Protein Data Bank (PDB) database and subsequently collected relevant parameters. We then employed Schrodinger to analyze the bond types of interaction and Gromacs 2019 to evaluate the TCR-antigen peptide complex’s molecular dynamics simulation. Finally, we compared the changes of bond types of interaction and binding free energy before and after residue substitution to ensure consistency of the conditions before and after residue substitution.

Results

The main sites on the antigen peptides that formed the intermolecular interaction [hydrogen bond (HB) and pi stack] with TCRs were P4, P8, P2, and P6. The hydrophobicity of the amino acids inside or outside the disulfide bond of TCRs may be related to the intermolecular interaction and binding free energy between TCRs and peptides. Residues located outside the disulfide bond of TCR α or β chains and forming pi stack force played favorable roles in the complex intermolecular interaction and binding free energy. The residues of the TCR α or β chains that interacted with peptides were replaced by alanine (Ala) or glycine (Gly), and their intermolecular binding free energy of the complex had been improved. However, it had nothing to do with the formation of HB.

Conclusions

The findings of this study suggest that the hydrophobic nature of the amino acids inside or outside the disulfide bonds on the TCR may be associated with the intermolecular interaction and binding between the TCR and polypeptide. The residues located outside the TCR α or β single-chain disulfide bond and forming the pi-stack force showed a beneficial effect on the intermolecular interaction and binding of the complex. In addition, the part of the residues on the TCR α or β single chain that produced bond types of interaction with the polypeptide after being replaced by Ala or Gly, the intermolecular binding free energy of the complex was increased, regardless of whether HB was formed.

Molecular dynamics simulations provide molecular insights into the role of HLA-B51 in Behçet's disease pathogenesis

Behçet's disease is an inflammatory disorder of unknown etiology. Genetic tendency has an important role in its pathogenesis, and HLA-B51, a class I MHC antigen, has been recognized as the strongest susceptibility factor for Behçet's disease. Despite the confirmation of the association of HLA-B51 with Behçet's disease in different populations, its pathogenic mechanisms remain elusive. HLA-B51 differs in only two amino acids from HLA-B52, other split antigen of HLA-B5, which is not associated with Behçet's disease. These two amino acids are located in the B pocket of the antigen-binding groove, which occupies the second amino acids of the bound peptides. To understand the nature of the HLA-peptide interactions, differences in structure and dynamics of two HLA alleles were investigated by molecular dynamics simulations using YAYDGKDYI, LPRSTVINI, and IPYQDLPHL peptides. For HLA-B51, all bound peptides fluctuated to larger extent than HLA-B52. Free energy profiles of unbinding process for YAYDGKDYI by steered molecular dynamics simulations showed that unbinding from HLA-B52 results in greater free energy differences than HLA-B51. These results suggest the possibility of an instability of HLA-B51 associated with the repertoire of peptides, and this finding may provide significant insight to its pathogenic role in Behçet's disease.

Methods for sequence and structural analysis of B and T cell receptor repertoires

B cell receptors (BCRs) and T cell receptors (TCRs) make up an essential network of defense molecules that, collectively, can distinguish self from non-self and facilitate destruction of antigen-bearing cells such as pathogens or tumors. The analysis of BCR and TCR repertoires plays an important role in both basic immunology as well as in biotechnology. Because the repertoires are highly diverse, specialized software methods are needed to extract meaningful information from BCR and TCR sequence data. Here, we review recent developments in bioinformatics tools for analysis of BCR and TCR repertoires, with an emphasis on those that incorporate structural features. After describing the recent sequencing technologies for immune receptor repertoires, we survey structural modeling methods for BCR and TCRs, along with methods for clustering such models. We review downstream analyses, including BCR and TCR epitope prediction, antibody-antigen docking and TCR-peptide-MHC Modeling. We also briefly discuss molecular dynamics in this context.

MHC binding affects the dynamics of different T-cell receptors in different ways

T cells use their T-cell receptors (TCRs) to scan other cells for antigenic peptides presented by MHC molecules (pMHC). If a TCR encounters a pMHC, it can trigger a signalling pathway that could lead to the activation of the T cell and the initiation of an immune response. It is currently not clear how the binding of pMHC to the TCR initiates signalling within the T cell. One hypothesis is that conformational changes in the TCR lead to further downstream signalling. Here we investigate four different TCRs in their free state as well as in their pMHC bound state using large scale molecular simulations totalling 26 000 ns. We find that the dynamical features within TCRs differ significantly between unbound TCR and TCR/pMHC simulations. However, apart from expected results such as reduced solvent accessibility and flexibility of the interface residues, these features are not conserved among different TCR types. The presence of a pMHC alone is not sufficient to cause cross-TCR-conserved dynamical features within a TCR. Our results argue against models of TCR triggering involving conserved allosteric conformational changes.

The interaction between T-cells and other cells is one of the most important interactions in the human immune system. If T-cells are not triggered major parts of the immune system cannot be activated or are not working effectively. Despite many years of research the exact mechanism of how a T-cell is initially triggered is not clear. One hypothesis is that conformational changes within the T-cell receptor (TCR) can cause further downstream signalling within the T-cell. In this study we computationally investigate the dynamics of four different TCRs in their free and bound configuration. Our large scale simulations show that all four TCRs react to binding in different ways. In some TCRs mainly the areas close to the binding region are affected while in other TCRs areas further apart from the binding region are also affected. Our results argue against a conserved structural activation mechanism across different types of TCRs.