The TCR binding site does move

Human Cytomegalovirus variant peptides adapt by decreasing their total coordination upon binding to a T cell receptor

The tertiary structure of the native Cytomegalovirus peptide (NLV) presented by HLA-A2 and bound to the RA14 T cell receptor was used as a reference for the calculation of atomic coordination differences of both the NLV as well as of a number of singly substituted NLV variants in the absence of TCR. Among the pMHC complexes, the native peptide was found to exhibit the highest total coordination difference in respect to the reference structure, suggesting that it experienced the widest structural adaptation upon recognition by the TCR. In addition, the peptide on the isolated NLV-MHC complex was over-coordinated as compared to the rest of the variants. Moreover, the trend was found to account for a set of measured dissociation constants and critical concentrations for target-cell lysis for all variants in complexation with RA14: functionally, all variant peptides were established to be either weak agonists or null peptides, while, at the same time, our current study established that they were also under-coordinated in respect to NLV. It could, thus, be argued that the most ‘efficient’ structural adaptation upon pMHC recognition by the TCR requires of the peptide to undergo the widest under-coordination possible. The main structural characteristic which differentiated the NLV in respect to the variants was a the presence of 16 oxygen atoms (waters) in the former׳s second coordination shell which accounted for over-coordination of roughly 100% and 30% in the O–O and C–O partials respectively. In fact, in the absence of second shell oxygens, the NLV peptide was decidedly under-coordinated in respect to all of the variants, as also suggested by the C–C partial.

Atomic Coordination Reflects Peptide Immunogenicity

We demonstrated that the immunological identity of variant peptides may be accurately predicted on the basis of atomic coordination of both unprotonated and protonated tertiary structures, provided that the structure of the native peptide (index) is known. The metric which was discovered to account for this discrimination is the coordination difference between the variant and the index; we also showed that increasing coordination difference in respect to the index was correlated to a correspondingly weakening immunological outcome of the variant. Additionally, we established that this metric quickly seizes to operate beyond the peptide scale, e.g., within a coordination shell inclusive of atoms up to a distance of 7 Å away from the peptide or over the entire pMHC-TCR complex. Analysis of molecular orbital interactions for a range of formal charges further revealed that the N-terminus of the agonists was always able to sustain a stable ammonium (NH3+) group which was consistently absent in antagonists. We deem that the presence of NH3+ constitutes a secondary observable with a biological consequence, signifying a change in T cell activation. While our analysis of protonated structures relied on the quantum chemical relaxation of the H species, the results were consistent across a wide range of peptide charge and spin polarization conditions.

The quantum chemical causality of pMHC-TCR biological avidity: Peptide atomic coordination data and the electronic state of agonist N termini

The quantum state of functional avidity of the synapse formed between a peptide-Major Histocompatibility Complex (pMHC) and a T cell receptor (TCR) is a subject not previously touched upon. Here we present atomic pair correlation meta-data based on crystalized tertiary structures of the Tax (HTLV-1) peptide along with three artificially altered variants, all of which were presented by the (Class I) HLA-A201 protein in complexation with the human (CD8(+)) A6TCR. The meta-data reveal the existence of a direct relationship between pMHC-TCR functional avidity (agonist/antagonist) and peptide pair distribution function (PDF). In this context, antagonist peptides are consistently under-coordinated in respect to Tax. Moreover, Density Functional Theory (DFT) datasets in the BLYP/TZ2P level of theory resulting from relaxation of the H species on peptide tertiary structures reveal that the coordination requirement of agonist peptides is also expressed as a physical observable of the protonation state of their N termini: agonistic peptides are always found to retain a stable ammonium (NH3 (+)) terminal group while antagonist peptides are not.

Quantum chemical calculations predict biological function: the case of T cell receptor interaction with a peptide/MHC class I

A combination of atomic correlation statistics and quantum chemical calculations are shown to predict biological function. In the present study, various antigenic peptide-Major Histocompatibility Complex (pMHC) ligands with near-identical stereochemistries, in complexation with the same T cell receptor (TCR), were found to consistently induce distinctly different quantum chemical behavior, directly dependent on the peptide's electron spin density and intrinsically expressed by the protonation state of the peptide's N-terminus. Furthermore, the cumulative coordination difference of any variant in respect to the native peptide was found to accurately reflect peptide biological function and immerges as the physical observable which is directly related to the immunological end-effect of pMHC-TCR interaction.

Early relaxation dynamics in the LC 13 T cell receptor in reaction to 172 altered peptide ligands: a molecular dynamics simulation study

The interaction between the T cell receptor and the major histocompatibility complex is one of the most important events in adaptive immunology. Although several different models for the activation process of the T cell via the T cell receptor have been proposed, it could not be shown that a structural mechanism, which discriminates between peptides of different immunogenicity levels, exists within the T cell receptor. In this study, we performed systematic molecular dynamics simulations of 172 closely related altered peptide ligands in the same T cell receptor/major histocompatibility complex system. Statistical evaluations yielded significant differences in the initial relaxation process between sets of peptides at four different immunogenicity levels.

T-cell receptors binding orientation over peptide/MHC class I is driven by long-range interactions

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.