Modeling T cell receptor recognition of CD1-lipid and MR1-metabolite complexes
© Pierce et al.; licensee BioMed Central Ltd. 2014
Received: 15 June 2014
Accepted: 22 September 2014
Published: 26 September 2014
T cell receptors (TCRs) can recognize diverse lipid and metabolite antigens presented by MHC-like molecules CD1 and MR1, and the molecular basis of many of these interactions has not been determined. Here we applied our protein docking algorithm TCRFlexDock, previously developed to perform docking of TCRs to peptide-MHC (pMHC) molecules, to predict the binding of αβ and γδ TCRs to CD1 and MR1, starting with the structures of the unbound molecules.
Evaluating against TCR-CD1d complexes with crystal structures, we achieved near-native structures in the top 20 models for two out of four cases, and an acceptable-rated prediction for a third case. We also predicted the structure of an interaction between a MAIT TCR and MR1-antigen that has not been structurally characterized, yielding a top-ranked model that agreed remarkably with a characterized TCR-MR1-antigen structure that has a nearly identical TCR α chain but a different β chain, highlighting the likely dominance of the conserved α chain in MR1-antigen recognition. Docking performance was improved by re-training our scoring function with a set of TCR-pMHC complexes, and for a case with an outlier binding mode, we found that alternative docking start positions improved predictive accuracy. We then performed unbound docking with two mycolyl-lipid specific TCRs that recognize lipid-bound CD1b, which represent a class of interactions that is not structurally characterized. Highly-ranked models of these complexes showed remarkable agreement between their binding topologies, as expected based on their shared germline sequences, while differences in residue-level interactions with their respective antigens point to possible mechanisms underlying their distinct specificities.
Together these results indicate that flexible docking simulations can provide accurate models and atomic-level insights into TCR recognition of MHC-like molecules presenting lipid and other small molecule antigens.
T cell receptors (TCRs) display remarkable versatility in their ability to specifically recognize a wide array of structurally and chemically diverse antigens. This is highlighted by a number of studies showing that, in addition to well-characterized recognition of peptide-major histocompatibility complex (pMHC) , many TCRs engage MHC-like molecules CD1, which present a variety of lipid antigens , and MR1, which presents vitamin B metabolites . At the sequence level, these TCRs can be restricted or diverse in their germ-line chain usage, depending in part on the T cell type, antigen, and antigen-presenting molecule .
The understanding of the molecular basis of TCR interactions with CD1 and MR1 has been greatly advanced by a number of crystallographic studies that have elucidated the interface sites, key contacts, and binding modes of several of these complexes. This includes type I and type II natural killer T cell (NKT) TCR interactions with CD1d-lipid [4–6], mucosal-associated invariant T cell (MAIT) TCR interfaces with unliganded  and ligand-bound MR1 , and recent studies of two γδ TCRs bound to CD1d-lipid [9, 10]. While these have provided a clear view of the antigen recognition underlying several of the invariant and diverse T cell subsets, owing to their variety and the effort required in experimental structure determination, the crystal structures of a number of key complexes have not yet been solved.
Here we describe the adaptation of the TCR docking algorithm, TCRFlexDock, previously shown to produce highly accurate TCR-pMHC structural predictions , to predict complexes of TCRs with CD1 and MR1 with antigen (Ag). We assessed unbound docking performance using four known structures of TCRs with CD1d, showing that for most cases accurate models were produced. We also predicted the complex between a MAIT TCR and MR1-Ag, an interaction that has not yet been structurally described, and assessing models against a related MAIT TCR bound to MR1-Ag, found strong agreement between predicted complex and the likely conserved binding mode. Applying TCRFlexDock to predict the structure of two germline-encoded, mycolyl lipid-reactive (GEM) TCRs bound to CD1b, we produced models that provide insight into both their shared sequence features and distinct antigen specificities.
Prediction of TCR-CD1d and TCR-MR1 complexes
The TCR-CD1d-Ag and TCR-MR1-Ag test cases
Binding RMSD1, Å
Type I NKT
Type II NKT
The predictive docking performance starting from the “start1” site
Start RMSD, Å
T20 RMSD3, Å
To investigate whether the TCR CDR loop backbone structure had a significant relationship with overall docking accuracy for these cases, we compared RMSD for modeled versus bound TCR residue backbone atoms (for the subset of TCR residues within 10 Å of CD1-Ag) with rigid-body ligand RMSD for the 1000 docking models of two test cases (Additional file 1: Figure S1). Surprisingly, there was little relationship between these values, indicating that although a bound-like TCR backbone will improve the hit rate of the docking simulation as we previously noted , the sampling of CDR loop conformations during docking includes bound-like conformations for both hits and non-hits, and a near-native docking position (i.e. low ligand RMSD) does not necessarily lead to accurately-modeled CDR loops.
Retraining the TCRFlexDock scoring function
We re-trained the ZRANK scoring function using the previously reported set of TCR-pMHC test cases and docking results  to determine whether such an optimized function would improve success on the TCR-CD1 and TCR-MR1 cases in Table 1. The re-weighted function we derived (named ZRT, for ZRANK TCR) indeed led to improved predictive performance in Table 2; though the first hit for case 4MNG was still ranked relatively low, its rank improved (rank 107, versus 235 for ZR2), while the near-hit was ranked 8 by ZRT (versus 10 by ZR2). The new TCR-pMHC derived weights had a lower van der Waals attractive weight relative to the other terms compared with the previous ZRANK function , likely due to the lower shape complementarity of TCR-pMHC interactions versus protein-protein interfaces in general . The top-ranked ZRT prediction for test case 4EI5 had a 2.21 Å interface RMSD from native; its structure is shown in Additional file 2: Figure S2. Despite some deviations of the sulfatide antigen head group and CDR3β with respect to the complex crystal structure, the XV19 TCR variable domains as a whole, as well as several key interface side chains, are positioned similarly to those of the known complex.
Utilizing alternate docking start sites
Prediction of TCR complexes with CD1b-GMM and CD1b-MA
Recently the unbound crystal structures of two clones of germline mycolyl-lipid reactive (GEM) αβ TCRs were reported : clone 18 (PDB code 4G8E), which binds CD1b-MA, and clone 42 (PDB code 4G8F) which binds CD1b-GMM. As the complex structures have not been reported, and the structure of CD1b-GMM was reported in a previous study , we used TCRFlexDock to predict these two structures representing this class of interactions. As the α chains are nearly identical between the two GEM TCR clones and likely dominate the interactions with CD1b-Ag , we computed the distances between the models ranked in the top 20 by ZRT for the two TCRs based on the RMSDs between their germline CDRα loop positions (with CD1b superposed; Additional file 6: Figure S5). We found that the 4G8E and 4G8F models (for brevity, we refer to unbound TCR PDB codes to represent these two complexes) were highly similar to each other, with many of the top 20 models less than 4 Å apart. 4G8E model 3 and 4G8F model 12 exhibited the highest similarity in germ line CDRα positioning (1.37 Å RMSD) over CD1b-Ag, and fell into the largest cluster of models in both sets of top 20 predictions (the bottom left cluster in Additional file 7: Figure S6). On this basis, as well as notable contacts with antigen as described below, we selected these models for further analysis.
We have demonstrated that our TCR docking algorithm, TCRFlexDock, is capable of producing accurate models of TCRs recognizing MHC-like molecules presenting a variety of non-peptide ligands. By retraining our scoring function, we were able to improve the re-ranking performance for test cases representing this class of complexes. Further improvements, for instance to improve the ranking of hits for the 4MNG or 3HUJ cases and performance using pooled results from alternate docking start sites, could potentially be achieved through inclusion of other terms (e.g. long-range partial charge rather than full charge electrostatics), and further optimization of the pair potential parameters and antigen atom typing. Additionally, though antigen bond torsion angles were minimized during Rosetta docking simulations, explicit inclusion of lipid head group rotamers (based on pre-generated conformers) could represent another route for improving docking performance.
We report models of CD1b-bound GEM TCRs from TCRFlexDock that account for a likely shared docking orientation as well fine differences in antigen recognition. The docking searches for these cases were focused on the CD1b site analogous to the pMHC interface engaged by most TCRs; while outlier TCR binding modes have been observed for some non-MHC ligands  and would reduce change of docking success (as seen for test case 3HUJ), mutagenesis evidence strongly suggests that TCRs engage CD1b-Ag with a footprint similar to TCR-pMHC interactions . Two previously reported TCR-CD1b-Ag modeling attempts (using other TCR clones) were based on homology with TCR-pMHC complex structures [25, 28], rather than a flexible docking search; notably, in one of these studies the TCR loops exhibited considerable clash with antigen in the modeled complexes , a challenge that is explicitly addressed during the conformational search of TCRFlexDock. Within our set of TCR-CD1b-Ag models, the dominant cluster was similar to the Class I MHC-bound Yae62 TCR (Figure 5 and Additional file 7: Figure S6). However, we did identify alternative models, including several that were similar to the TRBV6 MAIT TCR-MR1-Ag and ELS4 TCR-HLA-peptide complexes. As noted by others , the MAIT, ELS4 and GEM TCRs share TRAV1-2 genes and the former two exhibit similar docking orientations over their respective ligands , yet given the structural and residue-level differences between CD1b, MR1 and MHC, there is a significant possibility that TCR-CD1b-Ag interactions exhibit a distinct binding conformation, as we identified using TCRFlexDock. Experimental structural characterization of GEM TCR-CD1b-Ag complexes will allow determination of which models are most accurate.
While this manuscript was under review, the x-ray structure of a TRBV20-containing MAIT TCR (corresponding to the unbound TCR we used for test case 4L4T) in complex with MR1 and antigen was reported . Though the antigens of our docking simulation and the solved structure are not identical (thus preventing an analysis of native contacts and hits), we evaluated interface backbone RMSD between our models and the solved structure (PDB id 4PJ8), and the top-ranked model by ZRT score indeed had a low interface RMSD with respect to the experimentally determined complex structure (1.68 Å; Additional file 9: Figure S8).
These updates to the TCRFlexDock algorithm, as well as the docking test cases we assembled, represent an advance in modeling TCR-ligand recognition, extending capabilities beyond TCR-pMHC modeling methods reported by ourselves  and others [30, 31] to modeling complexes with lipid and metabolite antigens presented by MHC-like molecules. Our models of MAIT TCR-MR1-Ag and GEM TCR-CD1b-Ag complexes, which have not been described crystallographically, provide likely mechanisms of ligand engagement by these TCRs, and highlight the ability of advanced protein docking algorithms to complement experimental techniques in probing the structural basis of molecular recognition. These methods can be applied to model other TCRs bound to MHC-like ligands, as more of these interactions and TCRs become characterized, such as a γδ TCR clone that was recently found to engage the endothelial protein C receptor .
Input structures and TCR docking
TCR docking test cases were identified based on identical sequences between TCRs and ligands from unbound and bound PDB entries. Lipid and metabolite antigens were also evaluated for unbound and bound structural matches. For test case 4LHU, the lipid molecules were not identical (JLS for the bound CD1d-α-GalCer, AGH for the unbound CD1d-α-GalCer) but this was due to differences in buried hydrocarbon chains, while the exposed head groups were the same. Prior to docking, missing unbound TCR CDR loop residues were added to PDB files 4LFH (test case 4LHU), 4MNH (test case 4MNG), 4G8E, and 4G8F using Modeller , keeping the remainder of the TCR fixed. The unbound CD1b-MA structure was generated by truncation of the GMM antigen in the unbound CD1b-GMM structure (PDB code 1UQS).
Flexible TCR docking was performed in the same manner as presented previously , using a modified version of RosettaDock to generate 1000 docking models starting with unbound structures in a predetermined starting orientation (described in more detail below), followed by re-scoring docking models using ZRANK. RosettaDock simulations were performed using an iterative procedure (written with RosettaScripts) where the CDR3 loops undergo a perturbation using the kinetic closure (KIC) loop modeling protocol in Rosetta, followed by Monte Carlo-based sampling of local rigid-body docking orientation, side chain rotamer optimization, and two rounds of KIC refinement of all CDR loops. Parameters for lipid and antigens were generated using Open Babel  to convert from PDB to Mol2 format, followed by conversion to Rosetta parameter files using the molfile_to_params.py script . This process was also used to determine antigen torsion angles that were minimized by Rosetta during the docking process. Docking was performed on a Linux cluster with 2.8 GHz AMD Opteron cores; producing 50 docking models on a single core (20 such jobs were run in parallel for each test case) took approximately 10 hours on average (12 minutes per model). For ZRANK, antigen partial charge parameters were obtained from polar hydrogen Mol2 files from Open Babel, while IFACE and ACE atom types were assigned based on congruence with amino acid atoms.
Initial and alternate docking input positioning
As with our previous TCR docking study, for initial docking start position (corresponding to the “start1” site), TCRs were aligned with pseudo-symmetric axis perpendicular to the plane of the MHC-like helices, positioned over the MHC-like helix centroid at a distance of 25 Å, with a 45° crossing angle. We also tested an alternate “start2” docking start site; for this site, with respect to “start1”, the TCR was shifted by 20 Å along the MHC-like helix axis toward the α1 helix C-term, rotated to a 0° crossing angle, and tilted by 25° toward the MHC-like helix axis C-term. Additional “start3” and “start4” sites included 10 Å and 25 Å TCR shifts, crossing angles of 22.5° and -15°, and tilts of 12.5° and 40°, respectively.
Evaluation of docking models
Docking models were evaluated using CAPRI criteria to assess predictions as “incorrect”, “acceptable”, “medium accuracy” and “high accuracy” . As with our previous TCR docking study, models assessed as medium or high accuracy are referred to as hits. Lipid and metabolite antigen atoms were included in calculation of bound contacts with TCR residues, but were omitted from RMSD calculations (only protein backbone atoms were used).
Retraining ZRANK scoring function
As for the original ZRANK implementation , we employed a downhill simplex to select weights for the terms of the ZRANK scoring function, in this case maximizing the average area under the receiver operating characteristic curve (AUC) of 20 TCR-pMHC docking test cases , distinguishing hit from non-hit predictions among the top 30 ranked models for each TCR-pMHC case. To avoid local minima, the downhill simplex was run from 12,000 randomly generated starting points. For three-fold cross-validation within the TCR-pMHC set, training and testing sets were selected such that cases with the same TCR were not used simultaneously in both sets. The retrained scoring function weights (and ZR2 weight values  in parentheses, for comparison) used to select TCR docking models are:
van der Waals attractive: 0.027 (1.0)
van der Waals repulsive: 1.33 (0.23)
Electrostatics short-range attractive: 0.35 (0.57)
Electrostatics short-range repulsive: 0.29 (0.56)
Electrostatics long-range attractive: 0.30 (1.09)
Electrostatics long-range repulsive: 0.21 (0.29)
ACE: 0.64 (0.7)
IFACE: 0.27 (0.38)
Plots were generated using gnuplot (http://www.gnuplot.info), and molecular visualizations were produced using PyMOL (http://www.pymol.org). Clustering and generation of the unrooted tree were performed using the APE package  in R (http://www.r-project.org).
We thank Mary Ellen Gipson-Fitzpatrick (Boston University) for computing support, and Lawrence Stern for valuable discussions.
This work was funded by NIH grant GM103773.
- Rudolph MG, Stanfield RL, Wilson IA: How TCRs bind MHCs, peptides, and coreceptors. Annu Rev Immunol. 2006, 24: 419-466. 10.1146/annurev.immunol.23.021704.115658.View ArticlePubMedGoogle Scholar
- Adams EJ: Lipid presentation by human CD1 molecules and the diverse T cell populations that respond to them. Curr Opin Immunol. 2014, 26: 1-6.View ArticlePubMedGoogle Scholar
- Kjer-Nielsen L, Patel O, Corbett AJ, Le Nours J, Meehan B, Liu L, Bhati M, Chen Z, Kostenko L, Reantragoon R, Williamson NA, Purcell AW, Dudek NL, McConville MJ, O'Hair RA, Khairallah GN, Godfrey DI, Fairlie DP, Rossjohn J, McCluskey J: MR1 presents microbial vitamin B metabolites to MAIT cells. Nature. 2012, 491 (7426): 717-723.PubMedGoogle Scholar
- Borg NA, Wun KS, Kjer-Nielsen L, Wilce MC, Pellicci DG, Koh R, Besra GS, Bharadwaj M, Godfrey DI, McCluskey J, Rossjohn J: CD1d-lipid-antigen recognition by the semi-invariant NKT T-cell receptor. Nature. 2007, 448 (7149): 44-49. 10.1038/nature05907.View ArticlePubMedGoogle Scholar
- Patel O, Pellicci DG, Gras S, Sandoval-Romero ML, Uldrich AP, Mallevaey T, Clarke AJ, Le Nours J, Theodossis A, Cardell SL, Gapin L, Godfrey DI, Rossjohn J: Recognition of CD1d-sulfatide mediated by a type II natural killer T cell antigen receptor. Nat Immunol. 2012, 13 (9): 857-863. 10.1038/ni.2372.View ArticlePubMedGoogle Scholar
- Girardi E, Maricic I, Wang J, Mac TT, Iyer P, Kumar V, Zajonc DM: Type II natural killer T cells use features of both innate-like and conventional T cells to recognize sulfatide self antigens. Nat Immunol. 2012, 13 (9): 851-856. 10.1038/ni.2371.View ArticlePubMed CentralPubMedGoogle Scholar
- Lopez-Sagaseta J, Dulberger CL, Crooks JE, Parks CD, Luoma AM, McFedries A, Van Rhijn I, Saghatelian A, Adams EJ: The molecular basis for Mucosal-Associated Invariant T cell recognition of MR1 proteins. Proc Natl Acad Sci U S A. 2013, 110 (19): E1771-E1778. 10.1073/pnas.1222678110.View ArticlePubMed CentralPubMedGoogle Scholar
- Patel O, Kjer-Nielsen L, Le Nours J, Eckle SB, Birkinshaw R, Beddoe T, Corbett AJ, Liu L, Miles JJ, Meehan B, Reantragoon R, Sandoval-Romero ML, Sullivan LC, Brooks AG, Chen Z, Fairlie DP, McCluskey J, Rossjohn J: Recognition of vitamin B metabolites by mucosal-associated invariant T cells. Nat Commun. 2013, 4: 2142-PubMedGoogle Scholar
- Uldrich AP, Le Nours J, Pellicci DG, Gherardin NA, McPherson KG, Lim RT, Patel O, Beddoe T, Gras S, Rossjohn J, Godfrey DI: CD1d-lipid antigen recognition by the gammadelta TCR. Nat Immunol. 2013, 14 (11): 1137-1145. 10.1038/ni.2713.View ArticlePubMedGoogle Scholar
- Luoma AM, Castro CD, Mayassi T, Bembinster LA, Bai L, Picard D, Anderson B, Scharf L, Kung JE, Sibener LV, Savage PB, Jabri B, Bendelac A, Adams EJ: Crystal structure of Vdelta1 T cell receptor in complex with CD1d-sulfatide shows MHC-like recognition of a self-lipid by human gammadelta T cells. Immunity. 2013, 39 (6): 1032-1042. 10.1016/j.immuni.2013.11.001.View ArticlePubMedGoogle Scholar
- Pierce BG, Weng Z: A flexible docking approach for prediction of T cell receptor-peptide-MHC complexes. Protein Sci. 2013, 22 (1): 35-46. 10.1002/pro.2181.View ArticlePubMed CentralPubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.View ArticlePubMed CentralPubMedGoogle Scholar
- Eckle SB, Turner SJ, Rossjohn J, McCluskey J: Predisposed alphabeta T cell antigen receptor recognition of MHC and MHC-I like molecules?. Curr Opin Immunol. 2013, 25 (5): 653-659. 10.1016/j.coi.2013.07.010.View ArticlePubMedGoogle Scholar
- Lopez-Sagaseta J, Dulberger CL, McFedries A, Cushman M, Saghatelian A, Adams EJ: MAIT recognition of a stimulatory bacterial antigen bound to MR1. J Immunol. 2013, 191 (10): 5268-5277. 10.4049/jimmunol.1301958.View ArticlePubMedGoogle Scholar
- Pellicci DG, Patel O, Kjer-Nielsen L, Pang SS, Sullivan LC, Kyparissoudis K, Brooks AG, Reid HH, Gras S, Lucet IS, Koh R, Smyth MJ, Mallevaey T, Matsuda JL, Gapin L, McCluskey J, Godfrey DI, Rossjohn J: Differential recognition of CD1d-alpha-galactosyl ceramide by the V beta 8.2 and V beta 7 semi-invariant NKT T cell receptors. Immunity. 2009, 31 (1): 47-59. 10.1016/j.immuni.2009.04.018.View ArticlePubMed CentralPubMedGoogle Scholar
- Kjer-Nielsen L, Borg NA, Pellicci DG, Beddoe T, Kostenko L, Clements CS, Williamson NA, Smyth MJ, Besra GS, Reid HH, Bharadwaj M, Godfrey DI, Rossjohn J, McCluskey J: A structural basis for selection and cross-species reactivity of the semi-invariant NKT cell receptor in CD1d/glycolipid recognition. J Exp Med. 2006, 203 (3): 661-673. 10.1084/jem.20051777.View ArticlePubMed CentralPubMedGoogle Scholar
- Koch M, Stronge VS, Shepherd D, Gadola SD, Mathew B, Ritter G, Fersht AR, Besra GS, Schmidt RR, Jones EY, Cerundolo V: The crystal structure of human CD1d with and without alpha-galactosylceramide. Nat Immunol. 2005, 6 (8): 819-826. 10.1038/ni1225.View ArticlePubMedGoogle Scholar
- Zajonc DM, Maricic I, Wu D, Halder R, Roy K, Wong CH, Kumar V, Wilson IA: Structural basis for CD1d presentation of a sulfatide derived from myelin and its implications for autoimmunity. J Exp Med. 2005, 202 (11): 1517-1526. 10.1084/jem.20051625.View ArticlePubMed CentralPubMedGoogle Scholar
- Reantragoon R, Kjer-Nielsen L, Patel O, Chen Z, Illing PT, Bhati M, Kostenko L, Bharadwaj M, Meehan B, Hansen TH, Godfrey DI, Rossjohn J, McCluskey J: Structural insight into MR1-mediated recognition of the mucosal associated invariant T cell receptor. J Exp Med. 2012, 209 (4): 761-774. 10.1084/jem.20112095.View ArticlePubMed CentralPubMedGoogle Scholar
- Hwang H, Pierce B, Mintseris J, Janin J, Weng Z: Protein-protein docking benchmark version 3.0. Proteins. 2008, 73 (3): 705-709. 10.1002/prot.22106.View ArticlePubMed CentralPubMedGoogle Scholar
- Mendez R, Leplae R, De Maria L, Wodak SJ: Assessment of blind predictions of protein-protein interactions: current status of docking methods. Proteins. 2003, 52 (1): 51-67. 10.1002/prot.10393.View ArticlePubMedGoogle Scholar
- Pierce B, Weng Z: A combination of rescoring and refinement significantly improves protein docking performance. Proteins. 2008, 72 (1): 270-279. 10.1002/prot.21920.View ArticlePubMed CentralPubMedGoogle Scholar
- Garcia KC, Degano M, Pease LR, Huang M, Peterson PA, Teyton L, Wilson IA: Structural basis of plasticity in T cell receptor recognition of a self peptide-MHC antigen. Science. 1998, 279 (5354): 1166-1172. 10.1126/science.279.5354.1166.View ArticlePubMedGoogle Scholar
- Van Rhijn I, Kasmar A, de Jong A, Gras S, Bhati M, Doorenspleet ME, de Vries N, Godfrey DI, Altman JD, de Jager W, Rossjohn J, Moody DB: A conserved human T cell population targets mycobacterial antigens presented by CD1b. Nat Immunol. 2013, 14 (7): 706-713. 10.1038/ni.2630.View ArticlePubMedGoogle Scholar
- Batuwangala T, Shepherd D, Gadola SD, Gibson KJ, Zaccai NR, Fersht AR, Besra GS, Cerundolo V, Jones EY: The crystal structure of human CD1b with a bound bacterial glycolipid. J Immunol. 2004, 172 (4): 2382-2388. 10.4049/jimmunol.172.4.2382.View ArticlePubMedGoogle Scholar
- Adams EJ, Chien YH, Garcia KC: Structure of a gammadelta T cell receptor in complex with the nonclassical MHC T22. Science. 2005, 308 (5719): 227-231. 10.1126/science.1106885.View ArticlePubMedGoogle Scholar
- Melian A, Watts GF, Shamshiev A, De Libero G, Clatworthy A, Vincent M, Brenner MB, Behar S, Niazi K, Modlin RL, Almo S, Ostrov D, Nathenson SG, Porcelli SA: Molecular recognition of human CD1b antigen complexes: evidence for a common pattern of interaction with alpha beta TCRs. J Immunol. 2000, 165 (8): 4494-4504. 10.4049/jimmunol.165.8.4494.View ArticlePubMedGoogle Scholar
- Grant EP, Degano M, Rosat JP, Stenger S, Modlin RL, Wilson IA, Porcelli SA, Brenner MB: Molecular recognition of lipid antigens by T cell receptors. J Exp Med. 1999, 189 (1): 195-205. 10.1084/jem.189.1.195.View ArticlePubMed CentralPubMedGoogle Scholar
- Eckle SB, Birkinshaw RW, Kostenko L, Corbett AJ, McWilliam HE, Reantragoon R, Chen Z, Gherardin NA, Beddoe T, Liu L, Patel O, Meehan B, Fairlie DP, Villadangos JA, Godfrey DI, Kjer-Nielsen L, McCluskey J, Rossjohn J: A molecular basis underpinning the T cell receptor heterogeneity of mucosal-associated invariant T cells. J Exp Med. 2014, 211 (8): 1585-1600. 10.1084/jem.20140484.View ArticlePubMed CentralPubMedGoogle Scholar
- Kato Z, Stern JN, Nakamura HK, Kuwata K, Kondo N, Strominger JL: Positioning of autoimmune TCR-Ob.2 F3 and TCR-Ob.3D1 on the MBP85-99/HLA-DR2 complex. Proc Natl Acad Sci U S A. 2008, 105 (40): 15523-15528. 10.1073/pnas.0807338105.View ArticlePubMed CentralPubMedGoogle Scholar
- Leimgruber A, Ferber M, Irving M, Hussain-Kahn H, Wieckowski S, Derre L, Rufer N, Zoete V, Michielin O: TCRep 3D: an automated in silico approach to study the structural properties of TCR repertoires. PLoS One. 2011, 6 (10): e26301-10.1371/journal.pone.0026301.View ArticlePubMed CentralPubMedGoogle Scholar
- Willcox CR, Pitard V, Netzer S, Couzi L, Salim M, Silberzahn T, Moreau JF, Hayday AC, Willcox BE, Dechanet-Merville J: Cytomegalovirus and tumor stress surveillance by binding of a human gammadelta T cell antigen receptor to endothelial protein C receptor. Nat Immunol. 2012, 13 (9): 872-879. 10.1038/ni.2394.View ArticlePubMedGoogle Scholar
- Eswar N, Webb B, Marti-Renom MA, Madhusudhan MS, Eramian D, Shen MY, Pieper U, Sali A: Comparative protein structure modeling using Modeller. Curr Protoc Bioinform. 2006, Chapter 5: Unit 5-Unit 6.Google Scholar
- O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR: Open Babel: an open chemical toolbox. J Cheminform. 2011, 3: 33-10.1186/1758-2946-3-33.View ArticlePubMed CentralPubMedGoogle Scholar
- Lemmon G, Meiler J: Rosetta Ligand docking with flexible XML protocols. Methods Mol Biol. 2012, 819: 143-155. 10.1007/978-1-61779-465-0_10.View ArticlePubMed CentralPubMedGoogle Scholar
- Pierce B, Weng Z: ZRANK: reranking protein docking predictions with an optimized energy function. Proteins. 2007, 67 (4): 1078-1086. 10.1002/prot.21373.View ArticlePubMedGoogle Scholar
- Paradis E, Claude J, Strimmer K: APE: analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004, 20 (2): 289-290. 10.1093/bioinformatics/btg412.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.