T-Cell Receptor Binding Affects the Dynamics of the Peptide/MHC-I Complex

A Structural Voyage Toward the Landscape of Humoral and Cellular Immune Escapes of SARS-CoV-2

The genome-based surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the past nearly 5 years since its emergence has refreshed our understanding of virus evolution, especially on convergent co-evolution with the host. SARS-CoV-2 evolution has been characterized by the emergence of sets of mutations that affect the functional properties of the virus by altering its infectivity, virulence, transmissibility, and interactions with host immunity. This poses a huge challenge to global prevention and control measures based on drug treatment and vaccine application. As one of the key evasion strategies in response to the immune profile of the human population, there are overwhelming amounts of evidence for the reduced antibody neutralization of SARS-CoV-2 variants. Additionally, data also suggest that the levels of CD4+ and CD8+ T-cell responses against variants or sub-variants decrease in the populations, although non-negligible cross-T-cell responses are maintained. Herein, from the perspectives of structural immunology, we outline the characteristics and mechanisms of the T cell and antibody responses to SARS-CoV and its variants/sub-variants. The molecular bases for the impact of the immune escaping variants on the interaction of the epitopes with the key receptors in adaptive immunity, that is, major histocompatibility complex (MHC), T-cell receptor (TCR), and antibody are summarized and discussed, the knowledge of which will widen our understanding of this pandemic-threatening virus and assist the preparedness for Pathogen X in the future.

Quantifying conformational changes in the TCR:pMHC-I binding interface

Background

T cells form one of the key pillars of adaptive immunity. Using their surface bound T cell antigen receptors (TCRs), these cells screen millions of antigens presented by major histocompatibility complex (MHC) or MHC-like molecules. In other protein families, the dynamics of protein-protein interactions have important implications for protein function. Case studies of TCR:class I peptide-MHCs (pMHC-Is) structures have reported mixed results on whether the binding interfaces undergo conformational change during engagement and no robust statistical quantification has been done to generalise these results. Thus, it remains an open question of whether movement occurs in the binding interface that enables the recognition and activation of T cells.

Methods

In this work, we quantify the conformational changes in the TCR:pMHC-I binding interface by creating a dataset of 391 structures, comprising 22 TCRs, 19 MHC alleles, and 79 peptide structures in both unbound (apo) and bound (holo) conformations.

Results

In support of some case studies, we demonstrate that all complementarity determining region (CDR) loops move to a certain extent but only CDR3α and CDR3β loops modify their shape when binding pMHC-Is. We also map the contacts between TCRs and pMHC-Is, generating a novel fingerprint of TCRs on MHC molecules and show that the CDR3α tends to bind the N-terminus of the peptide and the CDR3β tends to bind the C-terminus of the peptide. Finally, we show that the presented peptides can undergo conformational changes when engaged by TCRs, as has been reported in past literature, but novelly show these changes depend on how the peptides are anchored in the MHC binding groove.

Conclusions

Our work has implications in understanding the behaviour of TCR:pMHC-I interactions and providing insights that can be used for modelling Tcell antigen specificity, an ongoing grand challenge in immunology.

Entropic analysis of antigen-specific CDR3 domains identifies essential binding motifs shared by CDR3s with different antigen specificities

Antigen-specific T cell receptor (TCR) sequences can have prognostic, predictive, and therapeutic value, but decoding the specificity of TCR recognition remains challenging. Unlike DNA strands that base pair, TCRs bind to their targets with different orientations and different lengths, which complicates comparisons. We present scanning parametrized by normalized TCR length (SPAN-TCR) to analyze antigen-specific TCR CDR3 sequences and identify patterns driving TCR-pMHC specificity. Using entropic analysis, SPAN-TCR identifies 2-mer motifs that decrease the diversity (entropy) of CDR3s. These motifs are the most common patterns that can predict CDR3 composition, and we identify "essential" motifs that decrease entropy in the same CDR3 α or β chain containing the 2-mer, and "super-essential" motifs that decrease entropy in both chains. Molecular dynamics analysis further suggests that these motifs may play important roles in binding. We then employ SPAN-TCR to resolve similarities in TCR repertoires against different antigens using public databases of TCR sequences.

Conformational flexibility of a free and TCR-bound pMHC-I protein investigated by long-term molecular dynamics simulations

BACKGROUND

Major histocompatibility complexes (MHCs) play a crucial role in the cell-mediated adaptive immune response as they present antigenic peptides (p) which are recognized by host T cells through a complex formation of the T cell receptor (TCR) with pMHC. In the present study, we report on changes in conformational flexibility within a pMHC molecule upon TCR binding by looking at molecular dynamics (MD) simulations of the free and the TCR-bound pMHC-I protein of the LC13-HLA-B*44:05-pEEYLQAFTY complex.

RESULTS

We performed long-term MD simulations with a total simulation time of 8 µs, employing 10 independent 400 ns replicas for the free and the TCR-bound pMHC system. Upon TCR ligation, we observed a reduced dynamic flexibility in the central residues of the peptide and the MHC α1-helix, altered occurrences of hydrogen bonds between the peptide and the MHC, a reduced conformational entropy of the peptide-binding groove, as well as a decreased solvent accessible surface area.

CONCLUSIONS

In summary, our results from 8 µs MD simulations indicate a restricted conformational space of the MHC peptide-binding groove upon TCR ligation and suggest a minimum simulation time of approximately 100 ns for biomolecules of comparable complexity to draw meaningful conclusions. Given the relatively long total simulation time, our results contribute to a more detailed view on conformational flexibility properties of the investigated free and TCR-bound pMHC-I system.

Peptide-dependent tuning of major histocompatibility complex motional properties and the consequences for cellular immunity

T cell receptors (TCRs) and other receptors of the immune system recognize peptides presented by class I or class II major histocompatibility complex (MHC) proteins. Although we generally distinguish between the MHC protein and its peptide, at an atomic level the two form a structural composite, which allows peptides to influence MHC properties and vice versa. One consequence is the peptide-dependent tuning of MHC structural dynamics, which contributes to protein structural adaptability and influences how receptors identify and bind targets. Peptide-dependent tuning of MHC protein dynamics can impact processes such as antigenicity, TCR cross-reactivity, and T cell repertoire selection. Motional tuning extends beyond the binding groove, influencing peptide selection and exchange, as well as interactions with other immune receptors. Here, we review recent findings showing how peptides can affect the dynamic and adaptable nature of MHC proteins. We highlight consequences for immunity and demonstrate how MHC proteins have evolved to be highly sensitive dynamic reporters, with broad immunological consequences.

Immunoinformatics approach of epitope prediction for SARS-CoV-2

BACKGROUND

The novel coronavirus (SARS-CoV-2) caused lethal infections worldwide during an unprecedented pandemic. Identification of the candidate viral epitopes is the first step in the design of vaccines against the viral infection. Several immunoinformatic approaches were employed to identify the SARS-CoV-2 epitopes that bind specifically with the major histocompatibility molecules class I (MHC-I). We utilized immunoinformatic tools to analyze the whole viral protein sequences, to identify the SARS-CoV-2 epitopes responsible for binding to the most frequent human leukocyte antigen (HLA) alleles in the Egyptian population. These alleles were also found with high frequency in other populations worldwide.

RESULTS

Molecular docking approach showed that using the co-crystallized MHC-I and T cell receptor (TCR) instead of using MHC-I structure only, significantly enhanced docking scores and stabilized the conformation, as well as the binding affinity of the identified SARS-CoV-2 epitopes. Our approach directly predicts 7 potential vaccine subunits from the available SARS-CoV-2 spike and ORF1ab protein sequence. This prediction has been confirmed by published experimentally validated and in silico predicted spike epitope. On the other hand, we predicted novel epitopes (RDLPQGFSA and FCLEASFNY) showing high docking scores and antigenicity response with both MHC-I and TCR. Moreover, antigenicity, allergenicity, toxicity, and physicochemical properties of the predicted SARS-CoV-2 epitopes were evaluated via state-of-the-art bioinformatic approaches, showing high efficacy of the proposed epitopes as a vaccine candidate.

CONCLUSION

Our predicted SARS-CoV-2 epitopes can facilitate vaccine development to enhance the immunogenicity against SARS-CoV-2 and provide supportive data for further experimental validation. Our proposed molecular docking approach of exploiting both MHC and TCR structures can be used to identify potential epitopes for most microbial pathogens, provided the crystal structure of MHC co-crystallized with TCR.

In-silico design of peptides for inhibition of HLA-A*03-KLIETYFSK complex as a new drug design for treatment of multiples sclerosis disease

Multiple sclerosis is recognized as a chronic inflammatory disease. Human leukocyte antigen (HLA) plays an important role in initiating adaptive immune responses. HLA class I is present in almost all nucleated cells and presents the cleaved endogenous peptide antigens to cytotoxic T cells. HLA-A*03 is one of the HLA class I alleles, which is reported as substantially related HLA to MS disease. In 2011, the structure of the HLA-A*03 in complex was identified with an immunodominant proteolipid protein (PLP) epitope (KLIETYFSK). This complex has been reported as an important autoantigen-presenting complex in MS pathogenesis. In this study, new peptides were designed to bind to this complex that may prevent specific pathogenic cytotoxic T cell binding to this autoantigen-presenting complex and CNS demyelination. Herein, 14 new helical peptides containing 19 amino acids were designed and their structures were predicted using the PEP-FOLD server. The binding of each designed peptide to the mentioned complex was then performed. A mutation approach was used by the BeAtMuSiC server to improve the binding affinity of the designed peptide. In each position, amino acid substitutions leading to an increase in the binding affinity of the peptide to the mentioned complex were determined. Finally, the resulting complexes were simulated for 40 ns using AMBER18 software. The results revealed that out of 14 designed peptides, "WRYWWKDWAKQFRQFYRWF" peptide exhibited the highest affinity for binding to the mentioned complex. This peptide can be considered as a potential drug to control multiple sclerosis disease in patients carrying the HLA-A*03 allele.

T Cell Epitope Prediction and Its Application to Immunotherapy

T cells play a crucial role in controlling and driving the immune response with their ability to discriminate peptides derived from healthy as well as pathogenic proteins. In this review, we focus on the currently available computational tools for epitope prediction, with a particular focus on tools aimed at identifying neoepitopes, i.e. cancer-specific peptides and their potential for use in immunotherapy for cancer treatment. This review will cover how these tools work, what kind of data they use, as well as pros and cons in their respective applications.

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.

Bioinformatics Prediction of SARS-CoV-2 Epitopes as Vaccine Candidates for the Colombian Population

Coronavirus disease (COVID-19) pandemic caused by the coronavirus SARS-CoV-2 represents an enormous challenge to global public health, with thousands of infections and deaths in over 200 countries worldwide. The purpose of this study was to identify SARS-CoV-2 epitopes with potential to interact in silico with the alleles of the human leukocyte antigen class I (HLA I) and class II (HLA II) commonly found in the Colombian population to promote both CD4 and CD8 immune responses against this virus. The generation and evaluation of the peptides in terms of HLA I and HLA II binding, immune response, toxicity and allergenicity were performed by using computer-aided tools, such as NetMHCpan 4.1, NetMHCIIpan 4.0, VaxiJem, ToxinPred and AllerTop. Furthermore, the interaction between the predicted epitopes with HLA I and HLA II proteins frequently found in the Colombian population was studied through molecular docking simulations in AutoDock Vina and interaction analysis in LigPlot+. One of the promising peptides proposed in this study is the HLA I epitope YQPYRVVVL, which displayed an estimated coverage of over 82% and 96% for the Colombian and worldwide population, respectively. These findings could be useful for the design of new epitope-vaccines that include Colombia among their population target.

Structure of Blood Coagulation Factor VIII in Complex With an Anti-C2 Domain Non-Classical, Pathogenic Antibody Inhibitor

Factor VIII (fVIII) is a procoagulant protein that binds to activated factor IX (fIXa) on platelet surfaces to form the intrinsic tenase complex. Due to the high immunogenicity of fVIII, generation of antibody inhibitors is a common occurrence in patients during hemophilia A treatment and spontaneously occurs in acquired hemophilia A patients. Non-classical antibody inhibitors, which block fVIII activation by thrombin and formation of the tenase complex, are the most common anti-C2 domain pathogenic inhibitors in hemophilia A murine models and have been identified in patient plasmas. In this study, we report on the X-ray crystal structure of a B domain-deleted bioengineered fVIII bound to the non-classical antibody inhibitor, G99. While binding to G99 does not disrupt the overall domain architecture of fVIII, the C2 domain undergoes an ~8 Å translocation that is concomitant with breaking multiple domain-domain interactions. Analysis of normalized B-factor values revealed several solvent-exposed loops in the C1 and C2 domains which experience a decrease in thermal motion in the presence of inhibitory antibodies. These results enhance our understanding on the structural nature of binding non-classical inhibitors and provide a structural dynamics-based rationale for cooperativity between anti-C1 and anti-C2 domain inhibitors.

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.

Exploring the different states of wild-type T-cell receptor and mutant conformational changes towards understanding the antigen recognition

Recognition of proteolytic peptide fragments presented by major histocompatibility complex (MHC) on target cells by T-cell receptor (TCR) is among the most important interactions in the adaptive immune system. Several computational studies have been performed to investigate conformational and dynamical properties of TCRs for enhanced immunogenicity. Here, we present the large-scale molecular dynamics (MD) simulation studies of the two comprehensive systems consisting of the wild-type and mutant IG4 TCR in complex with the tumor epitope NY-ESO peptide (SLLMWITQC) and analyzed for mapping conformational changes of TCR in the states prior to antigen binding, upon antigen binding and after the antigen was released. All of the simulations were performed with different states of TCRs for each 1000 ns of simulation time, providing six simulations for time duration of 6000 ns (6µs). We show that rather than undergoing most critical conformational changes upon antigen binding, the high proportion of complementarity-determining region (CDR) loops change by comparatively small amount. The hypervariable CDRα3 and CDRβ3 loops showed significant structural changes. Interestingly, the TCR β chain loops showed the least changes, which is reliable with recent implications that β domain of TCR may propel antigen interaction. The mutant shows higher rigidity than wild-type even in released state; expose an induced fit mechanism occurring from the re-structuring of CDRα3 loop and can allow enhanced binding affinity of the peptide antigen. Additionally, we show that CDRα3 loop and peptide contacts are an adaptive feature of affinity enhanced mutant TCR.Communicated by Ramaswamy H. Sarma.

Molecular mechanism of Be2+-ion binding to HLA-DP2: tetrahedral coordination, conformational changes and multi-ion binding

Be small and positive: the smaller size and higher charge of the Be²⁺-ion results in strong binding between the M2 peptide and the β-chain of HLA-DP2, which induces conformational changes at the periphery suitable for TCR binding.

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.

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

The interaction between T-cell receptors (TCRs) of T-cells and potentially immunogenic peptides presented by MHCs of antigen presenting cells is one of the most important mechanisms of the adaptive human immune system. A large number of structural simulations of the TCR/peptide/MHC system have been carried out. However, to date no study has investigated the differences of the dynamics between free TCRs and pMHC bound TCRs on a large scale. Here we present a study totalling 37 100 ns investigating the LC13 TCR in its free form as well as in complex with HLA-B*08:01 and different peptides. Our results show that the dynamics of the bound and unbound LC13 TCR differ significantly. This is reflected in (a) expected results such as an increased flexibility and increased solvent accessible surface of the CDRs of unbound TCR simulations but also in (b) less expected results such as lower CDR distances and compactness as well as alteration in the hydrogen bond network around CDR3α of unbound TCR simulations. Our study further emphasises the structural flexibility of TCRs and confirms the importance of the CDR3 loops for the adoption to MHC.

Avoiding False Positive Conclusions in Molecular Simulation: The Importance of Replicas

Molecular simulations are a computational technique used to investigate the dynamics of proteins and other molecules. The free energy landscape of these simulations is often rugged, and minor differences in the initial velocities, floating-point precision, or underlying hardware can cause identical simulations (replicas) to take different paths in the landscape. In this study we investigated the magnitude of these effects based on 310 000 ns of simulation time. We performed 100 identically parametrized replicas of 3000 ns each for a small 10 amino acid system as well as 100 identically parametrized replicas of 100 ns each for an 827 residue T-cell receptor/MHC system. Comparing randomly chosen subgroups within these replica sets, we estimated the reproducibility and reliability that can be achieved by a given number of replicas at a given simulation time. These results demonstrate that conclusions drawn from single simulations are often not reproducible and that conclusions drawn from multiple shorter replicas are more reliable than those from a single longer simulation. The actual number of replicas needed will always depend on the question asked and the level of reliability sought. On the basis of our data, it appears that a good rule of thumb is to perform a minimum of five to 10 replicas.

The partial dissociation of MHC class I-bound peptides exposes their N terminus to trimming by endoplasmic reticulum aminopeptidase 1

Endoplasmic reticulum aminopeptidase 1 (ERAP1) and ERAP2 process N-terminally extended antigenic precursors for optimal loading onto major histocompatibility complex class I (MHC I) molecules. We and others have demonstrated that ERAP1 processes peptides bound to MHC I, but the underlying mechanism is unknown. To this end, we utilized single-chain trimers (SCT) of the ovalbumin-derived epitope SIINFEKL (SL8) tethered to the H2-Kb MHC I determinant from mouse and introduced three substitutions, E63A, K66A, and W167A, at the A-pocket of the peptide-binding groove in the MHC I heavy chain, which interact with the N termini of peptides. These variants significantly decreased SL8-presenting SCT at the cell surface in the presence of ERAP1, but did not affect overall SCT expression, indicating that ERAP1 trims the SL8 N terminus. Comparison of the X-ray crystal structures of WT and three variant SCTs revealed only minor perturbations of the peptide-binding domain in the variants. However, molecular dynamics simulations suggested that SL8 can dissociate partially within a sub-microsecond timescale, exposing its N terminus to the solvent. We also found that the C terminus of MHC I-bound SL8 remains deeply buried in the F-pocket of MHC I. Furthermore, free-energy calculations revealed that the three SCT variants exhibit lower free-energy barriers of N terminus dissociation than the WT Kb Taken together, our results are consistent with a previously observed model in which the partial dissociation of bound peptides from MHC I exposes their N terminus to trimming by ERAP1, whereas their C terminus is anchored at the F-pocket.

Variable Regions of Antibodies and T-Cell Receptors May Not Be Sufficient in Molecular Simulations Investigating Binding

Antibodies and T-cell receptors are important proteins of the immune system that share similar structures. Both contain variable and constant regions. Insight into the dynamics of their binding can be provided by computational simulations. For these simulations the constant regions are often removed to save runtime as binding occurs in the variable regions. Here we present the first study to investigate the effect of removing the constant regions from antibodies and T-cell receptors on such simulations. We performed simulations of an antibody/antigen and T-cell receptor/MHC system with and without constant regions using 10 replicas of 100 ns of each of the four setups. We found that simulations without constant regions show significantly different behavior compared to simulations with constant regions. If the constant regions are not included in the simulations alterations in the binding interface hydrogen bonds and even partial unbinding can occur. These results indicate that constant regions should be included in antibody and T-cell receptor simulations for reliable conclusions to be drawn.

Partial Dissociation of Truncated Peptides Influences the Structural Dynamics of the MHCI Binding Groove

Antigen processing on MHCI involves the exchange of low-affinity peptides by high-affinity, immunodominant ones. This peptide editing process is mediated by tapasin and ERAAP at the peptide C- and N-terminus, respectively. Since tapasin does not contact the peptide directly, a sensing mechanism involving conformational changes likely allows tapasin to distinguish antigen-loaded MHCI molecules from those occupied by weakly bound, non-specific peptides. To understand this mechanism at the atomic level, we performed molecular dynamics simulations of MHCI allele B*44:02 loaded with peptides truncated or modified at the C- or N-terminus. We show that the deletion of peptide anchor residues leads to reversible, partial dissociation of the peptide from MHCI on the microsecond timescale. Fluctuations in the MHCI α2-1 helix segment, bordering the binding groove and cradled by tapasin in the PLC, are influenced by the peptide C-terminus occupying the nearby F-pocket. Simulations of tapasin complexed with MHCI bound to a low-affinity peptide show that tapasin widens the MHCI binding groove near the peptide C-terminus and weakens the attractive forces between MHCI and the peptide. Our simulations thus provide a detailed, spatially resolved picture of MHCI plasticity, revealing how peptide loading status can affect key structural regions contacting tapasin.

DynaDom: structure-based prediction of T cell receptor inter-domain and T cell receptor-peptide-MHC (class I) association angles

Background

T cell receptor (TCR) molecules are involved in the adaptive immune response as they distinguish between self- and foreign-peptides, presented in major histocompatibility complex molecules (pMHC). Former studies showed that the association angles of the TCR variable domains (Vα/Vβ) can differ significantly and change upon binding to the pMHC complex. These changes can be described as a rotation of the domains around a general Center of Rotation, characterized by the interaction of two highly conserved glutamine residues.

Methods

We developed a computational method, DynaDom, for the prediction of TCR Vα/Vβ inter-domain and TCR/pMHC orientations in TCRpMHC complexes, which allows predicting the orientation of multiple protein-domains. In addition, we implemented a new approach to predict the correct orientation of the carboxamide endgroups in glutamine and asparagine residues, which can also be used as an external, independent tool.

Results

The approach was evaluated for the remodeling of 75 and 53 experimental structures of TCR and TCRpMHC (class I) complexes, respectively. We show that the DynaDom method predicts the correct orientation of the TCR Vα/Vβ angles in 96 and 89% of the cases, for the poses with the best RMSD and best interaction energy, respectively. For the concurrent prediction of the TCR Vα/Vβ and pMHC orientations, the respective rates reached 74 and 72%. Through an exhaustive analysis, we could show that the pMHC placement can be further improved by a straightforward, yet very time intensive extension of the current approach.

Conclusions

The results obtained in the present remodeling study prove the suitability of our approach for interdomain-angle optimization. In addition, the high prediction rate obtained specifically for the energetically highest ranked poses further demonstrates that our method is a powerful candidate for blind prediction. Therefore it should be well suited as part of any accurate atomistic modeling pipeline for TCRpMHC complexes and potentially other large molecular assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-016-0071-7) contains supplementary material, which is available to authorized users.