ORIGINAL RESEARCH article
In Silico and Structural Analyses Demonstrate That Intrinsic Protein Motions Guide T Cell Receptor Complementarity Determining Region Loop Flexibility
- 1Division of Infection and Immunity and Systems Immunity Research Institute, Cardiff University School of Medicine, Cardiff, United Kingdom
- 2Immunocore, Abingdon, United Kingdom
- 3Department of Oncology, University Hospital of Lausanne, Lausanne, Switzerland
- 4Astbury Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom
- 5Department of Chemistry, University of Bath, Bath, United Kingdom
T-cell immunity is controlled by T cell receptor (TCR) binding to peptide major histocompatibility complexes (pMHCs). The nature of the interaction between these two proteins has been the subject of many investigations because of its central role in immunity against pathogens, cancer, in autoimmunity, and during organ transplant rejection. Crystal structures comparing unbound and pMHC-bound TCRs have revealed flexibility at the interaction interface, particularly from the perspective of the TCR. However, crystal structures represent only a snapshot of protein conformation that could be influenced through biologically irrelevant crystal lattice contacts and other factors. Here, we solved the structures of three unbound TCRs from multiple crystals. Superposition of identical TCR structures from different crystals revealed some conformation differences of up to 5 Å in individual complementarity determining region (CDR) loops that are similar to those that have previously been attributed to antigen engagement. We then used a combination of rigidity analysis and simulations of protein motion to reveal the theoretical potential of TCR CDR loop flexibility in unbound state. These simulations of protein motion support the notion that crystal structures may only offer an artifactual indication of TCR flexibility, influenced by crystallization conditions and crystal packing that is inconsistent with the theoretical potential of intrinsic TCR motions.
T-cells constitute our primary cellular defense against pathogenic challenge and play a major role in controlling neoplasms. The key molecular interface that enables T-cells to sense these threats is mediated by the clonally expressed T cell receptor (TCR) that classically distinguishes between self and foreign peptides. These peptides are derived from processed intra- and extra-cellular proteins, presented by highly diverse major histocompatibility complexes (pMHCs) on the surface of most nucleated cells (1). TCRs are required to respond to a vast number of potential foreign peptides that they have not encountered before, are unable to adapt to, and that can be presented by multiple MHCs (2, 3). More recently, it has been shown that the TCR can also recognize lipid antigens and metabolites presented by the invariant MHC-like cluster of differentiation 1d and MHC class I-related molecules, respectively (4, 5). Further evidence has implicated several other MHC-like molecules as antigenic targets for T-cells, exemplifying the extreme versatility of the TCR (6, 7).
In order to tackle this vast antigenic milieu, the gene rearrangement process that produces the TCR provides almost limitless possible TCR sequences through the recombination of TCR variable, joining and diversity genes (the germline encoded component), as well as addition and deletion of nucleotides (the somatic component). Additionally, through the combination of two chains (α and β) to form a heterodimeric αβ TCR, it is theoretically possible to generate ~1018 TCRs in humans. However, only a small fraction (<108) of these possibilities are ever expressed in any individual due to space limitations, suggesting that TCRs must be able to recognize multiple pMHCs to cover all potential pathogen encounters (8). Indeed, recent experimental evidence has confirmed this notion, demonstrating that TCRs can recognize millions of pMHCs with physiologically relevant sensitivity (9–13).
One key feature that is likely to facilitate this level of cross-reactivity is flexibility within the TCR complementarity-determining region loops (CDR loops) that form the antigen contact zone. Indeed, early thermodynamic evidence (14–16), combined with NMR spectroscopy (17–19) and fluorescence anisotropy (18, 20–22) has been used to indirectly and directly demonstrate TCR CDR loop motions. Alongside biophysical approaches, crystal structures comparing unbound and pMHC-bound TCRs have demonstrated that the TCR CDR loops can change shape, becoming stabilized upon binding (20, 23, 24). However, although X-ray crystallography provides unparalleled resolution of proteins too small for cryo-EM, the resulting snapshots represent only one conformational state that could be influenced by crystallographic artifacts. In order to extend the reach of atomic structures, computational modeling of protein motion has been used for almost 40 years to study a range of different systems ranging from enzymes, viral proteins, G protein coupled receptors and, more recently, immune receptors (21, 25–33). Application of this approach is beginning to shed light on the malleability of the TCR during pMHC binding, although questions still remain about the intrinsic flexibility of the TCR.
Here, we focused on two important unresolved questions. First, are conformational changes between unbound and pMHC-bound TCRs biologically relevant or crystal artifacts (especially considering that the CDR loops of unbound TCRs are more “exposed” for non-biologically relevant crystal contacts)? We addressed this issue by solving the structures of three different TCRs from multiple crystals. These included 12 datasets of the F11 TCR that recognizes a peptide from the influenza hemagglutinin protein (PKYVKQNTLKLAT) presented by HLA-DR*0101 (DR1-PKY), five data sets from the HA1.7 TCR that also recognizes DR1-PKY, and five datasets from the 003 TCR that recognizes an HIV p17 Gag-derived peptide (SLYNTVATL) presented by HLA-A*0201. These structures were compared to determine whether different loop conformations existed independent of ligand binding. Second, we investigated the intrinsic flexibility of the TCR CDR loops in unbound state by implementing geometric simulations of flexible TCR motion using a combination of rigidity analysis and coarse-grained elastic network normal mode analysis.
Overall, our data support the notion that crystal structures may represent an artifactual indication of TCR CDR loop flexibility that offers only a snapshot of the theoretical potential of TCR CDR loop motions. These data provide additional evidence contributing toward our understanding of the molecular mechanisms that mediate T-cell antigen discrimination and cross-reactivity.
Materials and Methods
Protein Expression, Refolding, and Purification
The F11, HA1.7, and 003 TCRs were generated as previously described (34), and the α and β chains were cloned into separate pGMT7 expression plasmids under the control of the T7 promoter. Each TCR was refolded and purified using methods that have been described previously (35).
Crystal Structure Determination
All protein crystals were grown at 18°C by vapor diffusion via the sitting drop technique. 200 nL of each TCR (10 mg/ml) in crystallization buffer (10 mM Tris pH 8.1 and 10 mM NaCl) was added to 200 nL of reservoir solution. The TCR crystals used in the structural investigations were grown in a variety of different conditions from PACT premier™ HT-96, JBScreen Classic HTS I, or TOPS (36) detailed in Table 1. Crystallization screens were conducted using an Art-Robbins Phoenix dispensing robot (Alpha Biotech Ltd., UK) and data were collected at 100 K at the diamond light source (DLS), Oxfordshire, UK using an ADSC Q315 CCD detector. Reflection intensities were estimated using XIA2 (37) and the data were analyzed with SCALA and the CCP4 package (38). Structures were solved with molecular replacement using PHASER (39). Sequences were adjusted with COOT (40) and the models were refined with REFMAC5. Graphical representations were prepared with PYMOL (41). Crystal contacts were determined using PYMOL and defined as intermolecular distances <4.0 Å. The reflection data and final model coordinates were deposited in the PDB database and are detailed in Tables 2–4.
Geometric Simulations of Flexible Motion
Amplitudes of motion in representative structures of the unbound TCRs solved here were simulated using a combination of rigidity analysis and coarse-grained elastic network normal mode analysis. Elnemo software (42) was used to obtain normal mode eigenvectors from coarse-grained elastic network modeling. FIRST/FRODA software (43, 44) was used to carry out rigidity analysis (FIRST) (45), which identified the noncovalent interaction network and labeling dihedral angles as locked or variable, and template-based geometric simulations of flexible motion (FRODA) (44) which project the all-atom structure over large amplitudes of motion, while maintaining local bonding and steric geometry.
Normal mode eigenvectors were generated in Elnemo in a one-site-per-residue coarse-graining using the Cα geometry of the input structure, placing springs of equal spring constant between all sites lying within an interaction distance cut-off of 12 Å. A rigidity analysis of the all-atom input structure was carried out in FIRST using the “pebble game” algorithm (43, 46), which matches degrees of freedom against bonding constraints in the molecular framework of the protein. Bonding constraints, include covalent, hydrophobic, and polar (hydrogen bond and salt bridge) interactions. As the strength of the polar interactions can be gauged from their geometry, the results of the analysis depend on an “energy cut-off” which selects the set of polar interactions to include in the constraint network (45). A cut-off of −3.0 kcal/mol was used in this study for simulations of flexible motion. We explored flexible motion biased along the 10 lowest-frequency nontrivial normal modes identified by Elnemo (modes 7–16; modes 1–6 are trivial rigid body motions).
Template-based geometric simulation of flexible motion, carried out using FRODA, explores the mobility of the all-atom structure by iterative perturbation and relaxation of atomic positions in parallel and antiparallel to the direction of normal mode eigenvectors. Several thousand iteration steps were carried out to generate large motion amplitudes. The simulation generates an initial phase of “easy” motion, where the bonding geometry is easily maintained, followed by the onset of “jamming” as the motion encounters steric and bonding constraints, which naturally limit its amplitude. The conformational changes of geometric simulations of TCRs projected using this method were observed, and compared with superpositions of crystallographic models in PyMOL.
TCR CDR Loop Flexibility Analysis Using X-Ray Crystallography
Several investigators, including ourselves, have previously used the structures of unbound and pMHC-bound TCRs to explore conformational changes that occur during pMHC ligand binding (20, 23, 24, 35). These studies have revealed that some TCRs undergo large conformational changes during binding, whereas others use a “lock-and-key”-type ligation strategy. However, questions remain over whether these changes accurately reflect how TCRs engage pMHC, or whether these observations are biased because they rely on a static image of a highly flexible protein interface that could be further affected by crystal lattice contacts, crystal packing, and/or crystallization conditions. To address this question, we solved multiple structures of three unbound TCRs at atomic resolution. We generated 12 structures of the F11 TCR, that recognizes a peptide from the influenza hemagglutinin protein (PKYVKQNTLKLAT) presented by HLA-DR*0101 (DR1-PKY), between 1.58 and 1.89 Å resolutions; 5 structures of the DR1-PKY specific HA1.7 TCR, between 2.31 and 2.98 Å resolution; and five structures of the 003 TCR, that recognizes an HIV-GAG-derived peptide (SLYNTVATL) presented by HLA-A*0201 (A2-SLY), between 1.26 and 1.37 Å resolution. All of the structures were solved with crystallographic Rwork/Rfree ratios within accepted limits as shown by the theoretically expected distribution (47). Statistical analysis and structure factors from two representative structures from each TCR are shown in Tables 2–4. The structures were refined by multiple individuals to avoid bias during refinements.
To accurately investigate TCR CDR loop movement during pMHC binding using crystal structures, the unbound TCR structures should, ideally, all be identical. We found that this was not the case for two of the three TCRs under investigation. Indeed, for the F11 TCR, we observed Cα backbone flexibility in the CDR2β loop, shifting by up to 3.6 Å in different structures of the same protein (Figure 1). For the HA1.7 TCR, we observed a larger shift in potential positions for the Cα backbone of the CDR3α loop, differing by up to 5 Å (Figure 2). In both cases, these shifts were comparable to loop movements reported in several other studies in which the unbound and pMHC-bound TCRs were compared (23). These findings demonstrate the potential for artifactual interpretation of the mechanism of TCR ligation using structures alone. However, this is not always the case. Indeed, the CDR loops of the third TCR included in our study, the 003 TCR, were superimposable in all the structures solved (Figure 3). B-factor analysis did not correlate with these loop movements (Figures 1–3).
Figure 1. Comparison of 12 unbound structures of the F11 T cell receptor (TCR). (A) Side-on view of the overall conformation of the F11 TCR (cartoon, colored by B-factor) including alignment of structures generated from 12 different crystals. (B) Top down view of an alignment of the complementarity determining region (CDR) loops (cartoon, colored by B-factor) from the 12 structures of unbound F11 TCR. The panels below demonstrate the two extreme conformations (conformation 1 and 2) of the CDR2β loop from the different structures, with the observed electron density map at 1σ. (C) Side-on view of the F11 TCR CDR2β loop (framework region in gray, apex of the loop in multiple colors) aligning all 12 structures with conformation 1 and conformation 2 labeled.
Figure 2. Comparison of five unbound structures of the HA1.7 T cell receptor (TCR). (A) Side-on view of the overall conformation of the HA1.7 TCR (cartoon, colored by B-factor) including alignment of structures generated from five different crystals. (B) Top down view of an alignment of the complementarity determining region (CDR) loops (cartoon, colored by B-factor) from the five structures of unbound HA1.7 TCR. The panels below demonstrate the two extreme conformations (conformation 1 and 2) of the CDR3α loop from the different structures, with the observed electron density map at 1σ. (C) Side-on view of the HA1.7 TCR CDR3α loop (framework region in gray, apex of the loop in multiple colors) aligning all five structures with conformation 1 and conformation 2 labeled.
Figure 3. Comparison of five unbound structures of the 003 T cell receptor (TCR). (A) Side-on view of the overall conformation of the 003 TCR (cartoon, colored by B-factor) including alignment of structures generated from five different crystals. (B) Top down view of an alignment of the complementarity determining region loops (cartoon, colored by B-factor) from the five structures of unbound 003 TCR.
Having observed positional differences in CDR loop positions in multiple, but not all, TCR datasets, we hypothesized that such differences might be explained by crystallographic artifacts. We therefore considered whether the resolution of the structures, or the crystal lattice contacts, had an impact on the nature of the loop movements observed in our structures. The 003 TCR structures, which exhibited practically identical CDR coordinates, were solved at an average resolution of 1.31 Å. In contrast, the HA1.7 TCR structures, in which a large loop movement was observed, were solved at 2.5 Å average resolution. This difference in resolution could partly explain why the interpretation of the position of the CDR3α loop of the HA.17 TCR varied between structures if the density was poor in this region. However, there was clear density supporting both conformations in the HA1.7 TCR structures (Figure 2B). Moreover, we also observed loop movement for the F11 TCR structures, which were solved at an average resolution of 1.55 Å, comparable to the higher resolution datasets for the 003 TCR. Last, we investigated whether stabilizing crystal lattice contacts might explain why some of the CDR loops appeared identical while others could shift between structures. For the F11 TCR, only the CDR2α and CDR2β loops were free from any lattice contacts (data not shown). Thus, the ability of the CDR2β loop to shift between structures could have been partly due to extra freedom imparted by individual crystal packing. Similarly, the HA1.7 TCR CDR3α loop was free from crystal lattice contacts and shifted between structures. However, several other loops in both the F11 and HA1.7 TCR structures were also free from crystal lattice contacts and did not shift between structures. Finally, none of the 003 TCR CDR loops, which were identical in each structure, made any potentially stabilizing crystal lattice contacts. Overall, neither the resolution of the structures, or the availability of stabilizing crystal lattice contacts were good predictors of TCR CDR loop shifts between structures. Thus, we conclude that CDR loop movements observed between unbound and pMHC-bound TCRs require further investigation in order to confirm whether they are artifactual, or a real part of the TCR-binding mechanism.
Investigating TCR CDR Loop Flexibility Using Rigidity Analysis and Simulations of Protein Motion
Direct measurements of protein flexibility at the single loop level is highly challenging because of the size (nm) and time (ms) of the movements. Indirect measurements using individually labeled amino acids are possible using NMR and other techniques, but these experiments are technically challenging, time consuming, and not universally available. As an alternative approach, computational modeling has developed rapidly over the past few years and has emerged as a useful technique to investigate protein motions (31). However, for T-cell recognition studies, most of these modeling approaches have focused on flexibility at the peptide-MHC, or TCR-pMHC interface rather than exploring the motional potential of the TCR in unbound state. In the one study that did use this method to investigate TCR flexibility, the authors found large differences in flexibility between two different TCRs that helped to explain the antigen-recognition mechanisms employed by each TCR, demonstrating the usefulness of this approach (21).
Here, we investigated the intrinsic rigidity and flexibility of the unbound TCR structure datasets using pebble-game rigidity analysis, elastic network modeling, and geometric simulations of flexible motion, using a combination of Elnemo and FIRST/FRODA software (44). FIRST software identifies the network of noncovalent constraints in the system, including both polar (hydrogen bond) and hydrophobic-tether interactions. Polar interactions are assigned strength in the range 0 to −10 kcal/mol based on their geometry. The set of polar interactions to include in the rigidity analysis is controlled by an energy cutoff parameter Ecut. We have demonstrated that biologically significant flexibility can be explored at cutoffs in the range −2 to −4 kcal/mol (48, 49). In order to determine the appropriate Ecut for TCRs, we performed rigidity analysis using the 003 TCR at cutoffs of −2 kcal/mol and −3 kcal/mol (data not shown). Analysis at a cutoff of −3 kcal/mol, but not −2 kcal/mol, demonstrated that the structure was largely flexible, with very few large rigid clusters. However, N and C terminal domains were still rich in noncovalent interactions maintaining the secondary, tertiary, and quaternary structure. We, therefore, explored flexible motion in all three TCR structures using the constraint network found at Ecut = −3 kcal/mol.
We used Elnemo software to identify the 10 lowest-frequency nontrivial normal modes in each TCR structure. We then used the FRODA module of FIRST to project the structure along each mode, while retaining the local covalent and noncovalent bonding geometry of the input structure. This represents intrinsic flexible motion of the structure which can easily be explored in solution. Recent work has shown that the character of motion identified using this method is consistent with conventional molecular dynamics (MD) simulations, while requiring minimal computational expense (a few CPU-hours) (49, 50). The lowest-frequency modes include substantial components of relative domain motions, in which the interdomain section of each chain (around residues 110–120) provides a flexible joint, as was also recently observed in the large dimeric enzyme Dcps (49). As a result, a structural overlay of conformations generated by FRODA (Figure 4) includes both domain-motion variations and local changes in loop geometry.
Figure 4. Flexibility analysis of the unbound F11, HA1.7, and 003 T cell receptor (TCRs). The TCRs were subjected to analysis using the FRODA/FIRST software packages to investigate protein flexibility. The top panels show the overall conformation of each TCR, using a single dataset, in which the TCR could flex over up to 2,000 frames, or until the run stopped because of steric clashes or overextended constraints. The bottom panel displays only the CDR loops of each TCR, with a distance scale to measure the motions, using the same analysis. (A) F11 TCR. (B) HA1.7 TCR. (C) 003 TCR. All structures are colored by protein domain, with the CDR1α in dark blue, CDR2α in blue, CDR3α in cyan, CDR1β in green, CDR2β in yellow, and CDR3β in orange.
To isolate the loop motions specifically, we carried out an alignment on the N-terminal domain for each structure and each chain (residues 6–110). The alignment was carried out in PyMOL on the non-loop residues of each domain (Table 5). This allowed visualization of CDR loop structural variations relative to a stable base of comparison, a set of 20 generated structural variants (Figure 5). For each normal mode, we selected the variants representing the natural limit of flexible motion parallel and antiparallel to the mode direction. This natural limit is the point at which covalent and noncovalent constraints (including steric contacts) start to limit the amplitude of the motion, such that further progress along the mode direction is “jammed”.
Figure 5. Flexible variations in the T cell receptor (TCR) complementarity determining region (CDR) loops using FRODA simulations. To isolate the loop motions, an alignment on the N-terminal domain for each TCR was performed (residues 6–110) to allow visualization of the loop structural variations relative to a stable base of comparison. Each panel shows a set of 20 structural variants aligned onto the initial crystal structure. Top panel TCRα CDR loops, bottom panel TCRβ CDR loops. For each normal mode, we selected the variants representing the natural limit of flexible motion parallel and antiparallel to the mode direction. The four loops are colored as follows: CDR1α in dark blue, CDR2α in blue, CDR3α in cyan, FWα in red, CDR1β in green, CDR2β in yellow, CDR3β in orange, and FWβ in purple. (A) F11 TCR. (B) HA1.7 TCR. (C) 003 TCR.
This analysis demonstrated that the scope for flexible variation in the loop geometries was substantial in all of the TCRs. Measurements of the maximal amplitude of the apex of each loop were conducted to provide estimations of the potential flexibility of each loop (Table 5). Although this analysis is an approximation and should be treated as such, the average maximal loop motion was slightly less (6.4 Å) for the 003 TCR compared to F11 (8.7 Å) and HA1.7 (9.0 Å), consistent with the structural data demonstrating greater rigidity in the 003 TCR. Further dissection of the data revealed that the longer somatically rearranged CDR3 loops had the most potential for loop motion (CDR3α: 10.2 Å, CDR3β: 12.2 Å) compared to the shorter germline encoded loops (CDR1α: 7.2 Å, CDR2α: 7.7 Å, Fwα: 8.7 Å, CDR1β:7.2 Å, CDR2β: 5.5 Å, Fwβ: 5.8 Å). These findings are consistent with the fact that the CDR3 loops generally make more interactions with the variable peptide component of the antigen, whereas the other loops are generally more focused toward the MHC surface. More generally, this analysis also indicated that all the structures, in solution, can explore large variations in loop geometry, providing an ensemble of flexible variations for conformational selection or induced-fit binding mechanisms. Overall, this analysis demonstrated, as expected, a large degree of potential motion, with the more disordered portions of the protein flexing more than those with secondary structure that was not apparent from the structural analysis (Figures 4 and 5).
The TCR governs T-cell specificity by discriminating between self and foreign peptides presented by MHC molecules. The finger-like CDR loops of the TCR are thought to meld around specific pMHCs, sampling the peptide cargo, and enabling T-cell triggering by ligands with sufficient affinity/dwell-time. This binding mode is also likely to facilitate TCR cross-reactivity by enabling the TCR to explore multiple conformations during ligand interrogation. However, the mechanism(s) that underpin the ability of T-cells to respond to millions of different pMHCs are still emerging. Several experiments have used atomic resolution structures to compare TCRs in unbound state and in complex with pMHC (23). These studies revealed conformational changes upon binding, supporting the idea that the CDR loops can flex to accommodate different peptide cargos using an induced fit mechanism. More recent data, using NMR, FRET, and MD support this view, but also demonstrate that the TCR-pMHC interface can be far more flexible than is apparent from the static image captured during X-ray crystallography (17–22).
First, we examined a very broad, but unanswered question: Is the conformation of a protein identical in every dataset collected and refined during X-ray crystallography experiments? This question is relevant to all structures solved by X-ray crystallography, but is particularly applicable when investigating ligand engagement by a receptor, as in the case of TCR-pMHC interaction. Several factors could affect the refined structure generated during this approach including; changes in the crystal packing between crystals, alterations in lattice contacts that could artificially stabilize protein regions, interpretation of the data during refinement, and differences in the protein preparation. In order to try to test some of these factors, we solved the structure of the same three TCRs from multiple crystals, grown in a range of conditions, from several different protein preparations, refined by different scientists.
Reassuringly, the overall conformation of each structure was virtually identical in all datasets tested. However, we observed several CDR loop re-organizations between structures for the HA1.7 and F11 TCRs, while all structures of the 003 TCR were identical at the level of the Cα backbone. These observations were not linked to the resolution of the structures, the crystal growing conditions, the availability of stabilizing crystal lattice contacts, or on who performed the refinement. Thus, we conclude that these loop movements represent real differences in the conformation of the CDR loops of the HA1.7 and F11 TCRs due to the intrinsic flexibility of these regions. We, therefore, recommend caution when using comparisons of unbound and pMHC-bound TCRs to describe binding mechanisms as these movements assume that the unbound structure of the TCR is a representative low energy state. Our observations suggest that the snapshot provided by X-ray crystallography may not be representative of CDR loop positions because of the highly dynamic nature of these disordered regions of the TCR.
To further investigate the flexible nature of the CDR loops, we measured large-amplitude protein motions in the three TCRs under investigation using FIRST/FRODA software. As expected, this analysis demonstrated a large degree of potential motion, with the more disordered portions of the protein flexing more than those with secondary structure. This analysis was far more revealing than the structural analysis alone, which only demonstrated structural mobility in some of the loops of just two of the TCRs (F11 and HA1.7). Rather, we observed large potential motions in all of the TCR CDR loops of all three TCRs, with more rigidity detected in the non-CDR loop portions of the TCR. This flexibility, which has been assumed, but not definitively proven for the TCR CDR loops, is consistent with the notion that the mechanism by which the TCR samples pMHC epitopes relies on flexibility at the interaction interface. Finer dissection of the motions of each individual CDR loop (including the Fw loop) demonstrated different maximal amplitudes at the apex of each loop. On average, the 003 TCR CDR loops moved slightly less compared to the F11 and HA1.7 TCR CDR loops, in line with the structural analysis. Furthermore, the somatically rearranged CDR3 loops in all of the TCRs studies were more mobile compared to the germline encoded CDR1, 2, and Fw loops. These findings are consistent with the observation that (1) the CDR3 loops are generally longer than the other loops, and (2) the CDR3 loops generally form the majority of the interactions with the variable peptide cargo, compared to the more MHC-centric CDR1, 2, and Fw loops. Thus, this extra level of flexibility in the CDR3 loops may represent an important mechanism enabling TCR cross-reactivity with multiple different peptides (2, 3, 10, 11).
These data have important implications for the general analysis of crystal structures, and more specifically for TCR antigen recognition. With the recent breakthroughs in MD and other modeling approaches, technologies that will rapidly develop in the near future, it seems rational to start pairing crystal structures with this type of analysis. Although theoretical, modeling approaches can provide another dimension of information to the complex and flexible amino acid network that governs the nature of protein–ligand dynamics. Findings from these analyses may reveal new areas of interest that can be tested experimentally. Our data, demonstrating the theoretical range of motion for unbound TCRs, are consistent with other modeling approaches that have focused on TCR-pMHC complexes (28–33), or TCRs alone (21). These studies have shown that the TCR-pMHC interface is highly flexible with some fixed interactions, but others that come and go as the TCR “rocks” on top of the pMHC (33). This binding mode is also congruous with recent experimental data demonstrating that TCRs are highly degenerate and can recognize many thousands, if not millions, of different peptide sequences (9–13). This enables T-cells to cross-react, thereby allowing a limited pool of TCR sequences within an individual to afford protection against the vast milieu of potential pathogenic peptide sequences that could be encountered (2, 3). Finally, our data reinforce the notion that some TCR CDR loops form a highly flexible and dynamic binding site, and that crystal structures alone may not be adequate to fully represent the complex mechanisms employed during pMHC ligation by the TCR. The fact that CDR loops can “move” between different free TCR structures of the same molecule suggests that caution is advisable when inferring binding mode based on the single snapshots provided by X-ray crystallography of ligated and unligated TCR.
PR, SW, and DC interpreted the data and designed experiments. CH, BM, VB, SH, RM, OV, AB, AF, PR, SW, and DC performed experiments. AG, AS, SW, and DC funded the study. AS, SW, and DC wrote the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
DC is a Wellcome Trust Career Development Fellow (WT095767). AS is a Wellcome Trust Senior Investigator (WT100327MA). BM is funded by The Wales Cancer Research Centre.
8. Arstila TP, Casrouge A, Baron V, Even J, Kanellopoulos J, Kourilsky P. A direct estimate of the human alphabeta T cell receptor diversity. Science (1999) 286:958–61. doi:10.1126/science.286.5441.958
9. Cole DK, Bulek AM, Dolton G, Schauenberg AJ, Szomolay B, Rittase W, et al. Hotspot autoimmune T cell receptor binding underlies pathogen and insulin peptide cross-reactivity. J Clin Invest (2016) 126:2191–204. doi:10.1172/JCI85679
10. Wooldridge L, Ekeruche-Makinde J, Van Den Berg HA, Skowera A, Miles JJ, Tan MP, et al. A single autoimmune T cell receptor recognizes more than a million different peptides. J Biol Chem (2012) 287:1168–77. doi:10.1074/jbc.M111.289488
12. Adams JJ, Narayanan S, Birnbaum ME, Sidhu SS, Blevins SJ, Gee MH, et al. Structural interplay between germline interactions and adaptive recognition determines the bandwidth of TCR-peptide-MHC cross-reactivity. Nat Immunol (2015) 17:87–94. doi:10.1038/ni.3310
13. Cole DK, van den Berg HA, Lloyd A, Crowther MD, Beck K, Ekeruche-Makinde J, et al. Structural mechanism underpinning cross-reactivity of a CD8+ T-cell clone that recognizes a peptide derived from human telomerase reverse transcriptase. J Biol Chem (2017) 292:802–13. doi:10.1074/jbc.M116.741603
14. Willcox BE, Gao GF, Wyer JR, Ladbury JE, Bell JI, Jakobsen BK, et al. TCR binding to peptide-MHC stabilizes a flexible recognition interface. Immunity (1999) 10:357–65. doi:10.1016/S1074-7613(00)80035-7
15. Boniface JJ, Reich Z, Lyons DS, Davis MM. Thermodynamics of T cell receptor binding to peptide-MHC: evidence for a general mechanism of molecular scanning. Proc Natl Acad Sci U S A (1999) 96:11446–51. doi:10.1073/pnas.96.20.11446
16. Krogsgaard M, Prado N, Adams EJ, He XL, Chow DC, Wilson DB, et al. Evidence that structural rearrangements and/or flexibility during TCR binding can contribute to T cell activation. Mol Cell (2003) 12:1367–78. doi:10.1016/S1097-2765(03)00474-X
17. Insaidoo FK, Zajicek J, Baker BM. A general and efficient approach for NMR studies of peptide dynamics in class I MHC peptide binding grooves. Biochemistry (2009) 48:9708–10. doi:10.1021/bi9008787
18. Hawse WF, De S, Greenwood AI, Nicholson LK, Zajicek J, Kovrigin EL, et al. TCR scanning of peptide/MHC through complementary matching of receptor and ligand molecular flexibility. J Immunol (2014) 192:2885–91. doi:10.4049/jimmunol.1302953
19. Varani L, Bankovich AJ, Liu CW, Colf LA, Jones LL, Kranz DM, et al. Solution mapping of T cell receptor docking footprints on peptide-MHC. Proc Natl Acad Sci U S A (2007) 104:13080–5. doi:10.1073/pnas.0703702104
20. Scott DR, Borbulevych OY, Piepenbrink KH, Corcelli SA, Baker BM. Disparate degrees of hypervariable loop flexibility control T-cell receptor cross-reactivity, specificity, and binding mechanism. J Mol Biol (2011) 414:385–400. doi:10.1016/j.jmb.2011.10.006
21. Ayres CM, Scott DR, Corcelli SA, Baker BM. Differential utilization of binding loop flexibility in T cell receptor ligand selection and cross-reactivity. Sci Rep (2016) 6:25070. doi:10.1038/srep25070
22. Scott DR, Vardeman CF, Corcelli SA, Baker BM. Limitations of time-resolved fluorescence suggested by molecular simulations: assessing the dynamics of T cell receptor binding loops. Biophys J (2012) 103:2532–40. doi:10.1016/j.bpj.2012.10.037
24. Holland CJ, Rizkallah PJ, Vollers S, Calvo-Calle JM, Madura F, Fuller A, et al. Minimal conformational plasticity enables TCR cross-reactivity to different MHC class II heterodimers. Sci Rep (2012) 2:629. doi:10.1038/srep00629
25. Ayres CM, Corcelli SA, Baker BM. Peptide and peptide-dependent motions in MHC proteins: immunological implications and biophysical underpinnings. Front Immunol (2017) 8:935. doi:10.3389/fimmu.2017.00935
27. Bermudez M, Mortier J, Rakers C, Sydow D, Wolber G. More than a look into a crystal ball: protein structure elucidation guided by molecular dynamics simulations. Drug Discov Today (2016) 21:1799–805. doi:10.1016/j.drudis.2016.07.001
28. Bello M, Correa-Basurto J. Energetic and flexibility properties captured by long molecular dynamics simulations of a membrane-embedded pMHCII–TCR complex. Mol Biosyst (2016) 12:1350–66. doi:10.1039/C6MB00058D
29. Knapp B, Dunbar J, Deane CM. Large scale characterization of the LC13 TCR and HLA-B8 structural landscape in reaction to 172 altered peptide ligands: a molecular dynamics simulation study. PLoS Comput Biol (2014) 10:e1003748. doi:10.1371/journal.pcbi.1003748
30. Knapp B, Dorffner G, Schreiner W. Early relaxation dynamics in the LC 13 T cell receptor in reaction to 172 altered peptide ligands: a molecular dynamics simulation study. PLoS One (2013) 8:e64464. doi:10.1371/journal.pone.0064464
33. Reboul CF, Meyer GR, Porebski BT, Borg NA, Buckle AM. Epitope flexibility and dynamic footprint revealed by molecular dynamics of a pMHC-TCR complex. PLoS Comput Biol (2012) 8:e1002404. doi:10.1371/journal.pcbi.1002404
35. Bulek AM, Cole DK, Skowera A, Dolton G, Gras S, Madura F, et al. Structural basis for the killing of human beta cells by CD8(+) T cells in type 1 diabetes. Nat Immunol (2012) 13:283–9. doi:10.1038/ni.2206
41. DeLano WL. The PyMOL Molecular Graphics System Version 1. Schrödinger LLC (2002). Available from: http://www.pymol.org
42. Suhre K, Sanejouand Y-H. ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement. Nucleic Acids Res (2004) 32:W610–4. doi:10.1093/nar/gkh368
45. Hespenheide BM, Rader AJ, Thorpe MF, Kuhn LA. Identifying protein folding cores from the evolution of flexible regions during unfolding. J Mol Graph Model (2002) 21:195–207. doi:10.1016/S1093-3263(02)00146-8
47. Tickle IJ, Laskowski RA, Moss DS. R(free) and the R(free) ratio. II. calculation of the expected values and variances of cross-validation statistics in macromolecular least-squares refinement. Acta Crystallogr D Biol Crystallogr (2000) 56:442–50. doi:10.1107/S0907444999016868
49. Wells SA, van der Kamp MW, McGeagh JD, Mulholland AJ. Structure and function in homodimeric enzymes: simulations of cooperative and independent functional motions. PLoS One (2015) 10:e0133372. doi:10.1371/journal.pone.0133372
Keywords: T-cells, T cell receptor, complementarity determining regions loops, protein flexibility, computational simulations, X-ray crystallography
Citation: Holland CJ, MacLachlan BJ, Bianchi V, Hesketh SJ, Morgan R, Vickery O, Bulek AM, Fuller A, Godkin A, Sewell AK, Rizkallah PJ, Wells S and Cole DK (2018) In Silico and Structural Analyses Demonstrate That Intrinsic Protein Motions Guide T Cell Receptor Complementarity Determining Region Loop Flexibility. Front. Immunol. 9:674. doi: 10.3389/fimmu.2018.00674
Received: 30 October 2017; Accepted: 19 March 2018;
Published: 11 April 2018
Edited by:Nick Gascoigne, National University of Singapore, Singapore
Reviewed by:Brian M. Baker, University of Notre Dame, United States
Dirk Michael Zajonc, La Jolla Institute for Allergy and Immunology (LJI), United States
Copyright: © 2018 Holland, MacLachlan, Bianchi, Hesketh, Morgan, Vickery, Bulek, Fuller, Godkin, Sewell, Rizkallah, Wells and Cole. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David K. Cole, firstname.lastname@example.org
†These authors have contributed equally to this work.