The Hypervariable Loops of Free TCRs Sample Multiple Distinct Metastable Conformations in Solution

CD4+ and CD8+ αβ T cell antigen recognition is determined by the interaction between the TCR Complementarity Determining Region (CDR) loops and the peptide-presenting MHC complex. These T cells are known for their ability to recognize multiple pMHC complexes, and for a necessary promiscuity that is required for their selection and function in the periphery. Crystallographic studies have previously elucidated the role of structural interactions in TCR engagement, but our understanding of the dynamic process that occurs during TCR binding is limited. To better understand the dynamic states that exist for TCR CDR loops in solution, and how this relates to their states when in complex with pMHC, we simulated the 2C T cell receptor in solution using all-atom molecular dynamics in explicit water and constructed a Markov State Model for each of the CDR3α and CDR3β loops. These models reveal multiple metastable states for the CDR3 loops in solution. Simulation data and metastable states reproduce known CDR3β crystal conformations, and reveal several novel conformations suggesting that CDR3β bound states are the result of search processes from nearby pre-existing equilibrium conformational states. Similar simulations of the invariant, Type I Natural Killer T cell receptor NKT15, which engages the monomorphic, MHC-like CD1d ligand, demonstrate that iNKT TCRs also have distinct states, but comparatively restricted degrees of motion.

[1]  Andrew K. Sewell,et al.  Why must T cells be cross-reactive? , 2012, Nature Reviews Immunology.

[2]  Albert C. Pan,et al.  Building Markov state models along pathways to determine free energies and rates of transitions. , 2008, The Journal of chemical physics.

[3]  Cory M. Ayres,et al.  Differential utilization of binding loop flexibility in T cell receptor ligand selection and cross-reactivity , 2016, Scientific Reports.

[4]  D. W. Scott,et al.  Multivariate Density Estimation, Theory, Practice and Visualization , 1992 .

[5]  D. Mason,et al.  A very high level of crossreactivity is an essential feature of the T-cell receptor. , 1998, Immunology today.

[6]  K. M. Armstrong,et al.  Thermodynamics of T‐cell receptor–peptide/MHC interactions: progress and opportunities , 2008, Journal of molecular recognition : JMR.

[7]  Charlotte M. Deane,et al.  Large Scale Characterization of the LC13 TCR and HLA-B8 Structural Landscape in Reaction to 172 Altered Peptide Ligands: A Molecular Dynamics Simulation Study , 2014, PLoS Comput. Biol..

[8]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[9]  Cyril F. Reboul,et al.  Epitope Flexibility and Dynamic Footprint Revealed by Molecular Dynamics of a pMHC-TCR Complex , 2012, PLoS Comput. Biol..

[10]  J. Rossjohn,et al.  Recognition of CD1d-restricted antigens by natural killer T cells , 2012, Nature Reviews Immunology.

[11]  Bernard Malissen,et al.  A T cell receptor CDR3beta loop undergoes conformational changes of unprecedented magnitude upon binding to a peptide/MHC class I complex. , 2002, Immunity.

[12]  B. Baker,et al.  Limitations of time-resolved fluorescence suggested by molecular simulations: assessing the dynamics of T cell receptor binding loops. , 2012, Biophysical journal.

[13]  B. Baker,et al.  Disparate degrees of hypervariable loop flexibility control T-cell receptor cross-reactivity, specificity, and binding mechanism. , 2011, Journal of molecular biology.

[14]  Fiona Ross,et al.  NKT cell activation by diverse antigens and antigen-presenting cells , 2013 .

[15]  Marc Bonneville,et al.  Kinetic evidence for a ligand-binding-induced conformational transition in the T cell receptor , 2007, Proceedings of the National Academy of Sciences.

[16]  E. Adams,et al.  Phosphoantigen-induced conformational change of butyrophilin 3A1 (BTN3A1) and its implication on Vγ9Vδ2 T cell activation , 2017, Proceedings of the National Academy of Sciences.

[17]  Thomas J Lane,et al.  MDTraj: a modern, open library for the analysis of molecular dynamics trajectories , 2014, bioRxiv.

[18]  Vijay S Pande,et al.  Improvements in Markov State Model Construction Reveal Many Non-Native Interactions in the Folding of NTL9. , 2013, Journal of chemical theory and computation.

[19]  Bernhard Moser,et al.  On the T , 2006, Fuzzy Sets Syst..

[20]  Natalie A. Borg,et al.  A structural basis for selection and cross-species reactivity of the semi-invariant NKT cell receptor in CD1d/glycolipid recognition , 2006, The Journal of experimental medicine.

[21]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[22]  Robyn L Stanfield,et al.  How TCRs bind MHCs, peptides, and coreceptors. , 2006, Annual review of immunology.

[23]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[24]  Sotaro Fuchigami,et al.  Slow dynamics in protein fluctuations revealed by time-structure based independent component analysis: the case of domain motions. , 2011, The Journal of chemical physics.

[25]  Ariel Fernández,et al.  WRAPPA : A Screening Tool for Candidate Dehydron Identification , 2011 .

[26]  Robyn L. Stanfield,et al.  An αβ T Cell Receptor Structure at 2.5 Å and Its Orientation in the TCR-MHC Complex , 1996, Science.

[27]  Y. Shoenfeld,et al.  [The T cell receptor]. , 1985, Harefuah.

[28]  Kurt H Piepenbrink,et al.  T cell receptor cross-reactivity directed by antigen-dependent tuning of peptide-MHC molecular flexibility. , 2009, Immunity.

[29]  S. Wells,et al.  In Silico and Structural Analyses Demonstrate That Intrinsic Protein Motions Guide T Cell Receptor Complementarity Determining Region Loop Flexibility , 2018, Front. Immunol..

[30]  Z Reich,et al.  Thermodynamics of T cell receptor binding to peptide-MHC: evidence for a general mechanism of molecular scanning. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Frank Noé,et al.  Markov models of molecular kinetics: generation and validation. , 2011, The Journal of chemical physics.

[32]  Yi Li,et al.  T-cell Receptor Specificity Maintained by Altered Thermodynamics* , 2013, The Journal of Biological Chemistry.

[33]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[34]  D. Kranz,et al.  T cell receptors: affinities, cross-reactivities, and a conformer model. , 2004, Molecular immunology.

[35]  Vijay S. Pande,et al.  Everything you wanted to know about Markov State Models but were afraid to ask. , 2010, Methods.

[36]  K. Garcia,et al.  A functional hot spot for antigen recognition in a superagonist TCR/MHC complex. , 2000, Immunity.

[37]  E Yvonne Jones,et al.  The structural dynamics and energetics of an immunodominant T cell receptor are programmed by its Vbeta domain. , 2008, Immunity.

[38]  Kathryn M Hart,et al.  Discovery of multiple hidden allosteric sites by combining Markov state models and experiments , 2015, Proceedings of the National Academy of Sciences.

[39]  K. Garcia,et al.  How a Single T Cell Receptor Recognizes Both Self and Foreign MHC , 2007, Cell.

[40]  P. Deuflhard,et al.  Robust Perron cluster analysis in conformation dynamics , 2005 .

[41]  L R Pease,et al.  Structural basis of plasticity in T cell receptor recognition of a self peptide-MHC antigen. , 1998, Science.

[42]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[43]  Jennifer Maynard,et al.  Structural evidence for a germline-encoded T cell receptor–major histocompatibility complex interaction 'codon' , 2007, Nature Immunology.

[44]  K. Christopher Garcia,et al.  How the T Cell Receptor Sees Antigen—A Structural View , 2005, Cell.

[45]  Arup K. Chakraborty,et al.  Fast on-rates allow short dwell time ligands to activate T cells , 2010, Proceedings of the National Academy of Sciences.

[46]  Diwakar Shukla,et al.  Transition path theory analysis of c-Src kinase activation , 2016, Proceedings of the National Academy of Sciences.

[47]  Mark M. Davis,et al.  Two-step binding mechanism for T-cell receptor recognition of peptide–MHC , 2002, Nature.

[48]  B. Baker,et al.  Structural and dynamic control of T‐cell receptor specificity, cross‐reactivity, and binding mechanism , 2012, Immunological reviews.

[49]  G. Besra,et al.  Human and Mouse Type I Natural Killer T Cell Antigen Receptors Exhibit Different Fine Specificities for CD1d-Antigen Complex* , 2012, The Journal of Biological Chemistry.

[50]  A. Sewell,et al.  Minimal conformational plasticity enables TCR cross-reactivity to different MHC class II heterodimers , 2012, Scientific Reports.

[51]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[52]  Ariel Fernández,et al.  Dehydron analysis: quantifying the effect of hydrophobic groups on the strength and stability of hydrogen bonds. , 2010, Advances in experimental medicine and biology.

[53]  B K Jakobsen,et al.  TCR binding to peptide-MHC stabilizes a flexible recognition interface. , 1999, Immunity.

[54]  Thomas J Lane,et al.  MSMBuilder2: Modeling Conformational Dynamics at the Picosecond to Millisecond Scale. , 2011, Journal of chemical theory and computation.

[55]  James McCluskey,et al.  Differential recognition of CD1d-alpha-galactosyl ceramide by the V beta 8.2 and V beta 7 semi-invariant NKT T cell receptors. , 2009, Immunity.

[56]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[57]  Clemencia Pinilla,et al.  Specificity and degeneracy of T cells. , 2004, Molecular immunology.