Pre-T cell receptor self-MHC sampling restricts thymocyte dedifferentiation.

Programming T cells to distinguish self from non-self is a vital, multi-step process that occurs in the thymus1-4. Signalling through the pre-T cell receptor (preTCR), a CD3-associated heterodimer comprising an invariant pTα chain and a clone-specific ß chain, is a critical early checkpoint in thymocyte development within the αß T cell lineage5,6. PreTCRs arrayed on CD4-CD8- double-negative thymocytes ligate peptides bound to major histocompatibility complex molecules (pMHC) on thymic stroma, similar to αß T cell receptors that appear on CD4+CD8+ double-positive thymocytes, but via a different molecular docking strategy7-10. Here we show the consequences of these distinct interactions for thymocyte progression using synchronized fetal thymic progenitor cultures that differ in the presence or absence of pMHC on support stroma, and single-cell transcriptomes at key thymocyte developmental transitions. Although major histocompatibility complex (MHC)-negative stroma fosters αß T cell differentiation, the absence of preTCR-pMHC interactions leads to deviant thymocyte transcriptional programming associated with dedifferentiation. Highly proliferative double-negative and double-positive thymocyte subsets emerge, with antecedent characteristics of T cell lymphoblastic and myeloid malignancies. Compensatory upregulation of diverse MHC class Ib proteins in B2m/H2-Ab1 MHC-knockout mice partially safeguards in vivo thymocyte progression, although disseminated double-positive thymic tumours may develop with ageing. Thus, as well as promoting ß chain repertoire broadening for subsequent αß T cell receptor utilization, preTCR-pMHC interactions limit cellular plasticity to facilitate normal thymocyte differentiation and proliferation that, if absent, introduce developmental vulnerabilities.

Harnessing αß T cell receptor mechanobiology to achieve the promise of immuno-oncology.

T cell receptors (TCR) on cytolytic T lymphocytes (CTLs) recognize "foreign" antigens bound in the groove of major histocompatibility complex (MHC) molecules (H-2 in mouse and HLA in human) displayed on altered cells. These antigens are peptide fragments of proteins derived either from infectious pathogens or cellular transformations during cancer evolution. The conjoint ligand formed by the foreign peptide and MHC, termed pMHC, marks an aberrant cell as a target for CTL-mediated destruction. Recent data have provided compelling evidence that adaptive protection is achieved in a facile manner during immune surveillance when mechanical load consequent to cellular motion is applied to the bond formed between an αß TCR and its pMHC ligand arrayed on a disease-altered cell. Mechanobiology maximizes both TCR specificity and sensitivity in comparison to receptor ligation in the absence of force. While the field of immunotherapy has made advances to impact the survival of cancer patients, the latest information relevant to T cell targeting and mechanotransduction has yet to be applied for T cell monitoring and treatment of patients in the clinic. Here we review these data, and challenge scientists and physicians to apply critical biophysical parameters of TCR mechanobiology to the medical oncology field, broadening treatment success within and among various cancer types. We assert that TCRs with digital ligand-sensing performance capability directed at sparsely as well as luminously displayed tumor-specific neoantigens and certain tumor-associated antigens can improve effective cancer vaccine development and immunotherapy paradigms.

Molecular design of the γδT cell receptor ectodomain encodes biologically fit ligand recognition in the absence of mechanosensing.

High-acuity αßT cell receptor (TCR) recognition of peptides bound to major histocompatibility complex molecules (pMHCs) requires mechanosensing, a process whereby piconewton (pN) bioforces exert physical load on αßTCR-pMHC bonds to dynamically alter their lifetimes and foster digital sensitivity cellular signaling. While mechanotransduction is operative for both αßTCRs and pre-TCRs within the αßT lineage, its role in γδT cells is unknown. Here, we show that the human DP10.7 γδTCR specific for the sulfoglycolipid sulfatide bound to CD1d only sustains a significant load and undergoes force-induced structural transitions when the binding interface-distal γδ constant domain (C) module is replaced with that of αß. The chimeric γδ-αßTCR also signals more robustly than does the wild-type (WT) γδTCR, as revealed by RNA-sequencing (RNA-seq) analysis of TCR-transduced Rag2-/- thymocytes, consistent with structural, single-molecule, and molecular dynamics studies reflective of γδTCRs as mediating recognition via a more canonical immunoglobulin-like receptor interaction. Absence of robust, force-related catch bonds, as well as γδTCR structural transitions, implies that γδT cells do not use mechanosensing for ligand recognition. This distinction is consonant with the fact that their innate-type ligands, including markers of cellular stress, are expressed at a high copy number relative to the sparse pMHC ligands of αßT cells arrayed on activating target cells. We posit that mechanosensing emerged over ∼200 million years of vertebrate evolution to fulfill indispensable adaptive immune recognition requirements for pMHC in the αßT cell lineage that are unnecessary for the γδT cell lineage mechanism of non-pMHC ligand detection.

The αßTCR mechanosensor exploits dynamic ectodomain allostery to optimize its ligand recognition site.

Each [Formula: see text]T cell receptor (TCR) functions as a mechanosensor. The TCR is comprised of a clonotypic TCR[Formula: see text] ligand-binding heterodimer and the noncovalently associated CD3 signaling subunits. When bound by ligand, an antigenic peptide arrayed by a major histocompatibility complex molecule (pMHC), the TCR[Formula: see text] has a longer bond lifetime under piconewton-level loads. The atomistic mechanism of this "catch bond" behavior is unknown. Here, we perform molecular dynamics simulation of a TCR[Formula: see text]-pMHC complex and its variants under physiologic loads to identify this mechanism and any attendant TCR[Formula: see text] domain allostery. The TCR[Formula: see text]-pMHC interface is dynamically maintained by contacts with a spectrum of occupancies, introducing a level of control via relative motion between Vα and Vß variable domains containing the pMHC-binding complementarity-determining region (CDR) loops. Without adequate load, the interfacial contacts are unstable, whereas applying sufficient load suppresses Vα-Vß motion, stabilizing the interface. A second level of control is exerted by Cα and Cß constant domains, especially Cß and its protruding FG-loop, that create mismatching interfaces among the four TCR[Formula: see text] domains and with a pMHC ligand. Applied load enhances fit through deformation of the TCR[Formula: see text] molecule. Thus, the catch bond involves the entire TCR[Formula: see text] conformation, interdomain motion, and interfacial contact dynamics, collectively. This multilayered architecture of the machinery fosters fine-tuning of cellular response to load and pMHC recognition. Since the germline-derived TCR[Formula: see text] ectodomain is structurally conserved, the proposed mechanism can be universally adopted to operate under load during immune surveillance by diverse [Formula: see text]TCRs constituting the T cell repertoire.

Molecular recognition of a host protein by NS1 of pandemic and seasonal influenza A viruses.

The 1918 influenza A virus (IAV) caused the most severe flu pandemic in recorded human history. Nonstructural protein 1 (NS1) is an important virulence factor of the 1918 IAV. NS1 antagonizes host defense mechanisms through interactions with multiple host factors. One pathway by which NS1 increases virulence is through the activation of phosphoinositide 3-kinase (PI3K) by binding to its p85ß subunit. Here we present the mechanism underlying the molecular recognition of the p85ß subunit by 1918 NS1. Using X-ray crystallography, we determine the structure of 1918 NS1 complexed with p85ß of human PI3K. We find that the 1918 NS1 effector domain (1918 NS1ED) undergoes a conformational change to bind p85ß. Using NMR relaxation dispersion and molecular dynamics simulation, we identify that free 1918 NS1ED exists in a dynamic equilibrium between p85ß-binding-competent and -incompetent conformations in the submillisecond timescale. Moreover, we discover that NS1ED proteins of 1918 (H1N1) and Udorn (H3N2) strains exhibit drastically different conformational dynamics and binding kinetics to p85ß. These results provide evidence of strain-dependent conformational dynamics of NS1. Using kinetic modeling based on the experimental data, we demonstrate that 1918 NS1ED can result in the faster hijacking of p85ß compared to Ud NS1ED, although the former has a lower affinity to p85ß than the latter. Our results suggest that the difference in binding kinetics may impact the competition with cellular antiviral responses for the activation of PI3K. We anticipate that our findings will increase the understanding of the strain-dependent behaviors of influenza NS1 proteins.

Structural basis for power stroke vs. Brownian ratchet mechanisms of motor proteins.

Two mechanisms have been proposed for the function of motor proteins: The power stroke and the Brownian ratchet. The former refers to generation of a large downhill free energy gradient over which the motor protein moves nearly irreversibly in making a step, whereas the latter refers to biasing or rectifying the diffusive motion of the motor. Both mechanisms require input of free energy, which generally involves the processing of an ATP (adenosine 5'-triphosphate) molecule. Recent advances in experiments that reveal the details of the stepping motion of motor proteins, together with computer simulations of atomistic structures, have provided greater insights into the mechanisms. Here, we compare the various models of the power stroke and the Brownian ratchet that have been proposed. The 2 mechanisms are not mutually exclusive, and various motor proteins employ them to different extents to perform their biological function. As examples, we discuss linear motor proteins Kinesin-1 and myosin-V, and the rotary motor F1-ATPase, all of which involve a power stroke as the essential element of their stepping mechanism.

Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain.

Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their "fuzzy" nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy "hotspots" as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.

NMR: an essential structural tool for integrative studies of T cell development, pMHC ligand recognition and TCR mechanobiology.

Early studies of T cell structural biology using X-ray crystallography, surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) focused on a picture of the αßT cell receptor (αßTCR) component domains and their cognate ligands (peptides bound to MHC molecules, i.e. pMHCs) as static interaction partners. Moving forward requires integrating this corpus of data with dynamic technologies such as NMR, molecular dynamics (MD) simulations and real-time single molecule (SM) studies exemplified by optical tweezers (OT). NMR bridges relevant timescales and provides the potential for an all-atom dynamic description of αßTCR components prior to and during interactions with binding partners. SM techniques have opened up vistas in understanding the non-equilibrium nature of T cell signaling through the introduction of force-mediated binding measurements into the paradigm for T cell function. In this regard, bioforces consequent to T-lineage cell motility are now perceived as placing piconewton (pN)-level loads on single receptor-pMHC bonds to impact structural change and αßT-lineage biology, including peptide discrimination, cellular activation, and developmental progression. We discuss herein essential NMR technologies in illuminating the role of ligand binding in the preT cell receptor (preTCR), the αßTCR developmental precursor, and convergence of NMR, SM and MD data in advancing our comprehension of T cell development. More broadly we review the central hypothesis that the αßTCR is a mechanosensor, fostered by breakthrough NMR-based structural insights. Collectively, elucidating dynamic aspects through the integrative use of NMR, SM, and MD shall advance fundamental appreciation of the mechanism of T cell signaling as well as inform translational efforts in αßTCR and chimeric T cell (CAR-T) immunotherapies and T cell vaccinology.

Effect of Methylation on Local Mechanics and Hydration Structure of DNA.

Cytosine methylation affects mechanical properties of DNA and potentially alters the hydration fingerprint for recognition by proteins. The atomistic origin for these effects is not well understood, and we address this via all-atom molecular dynamics simulations. We find that the stiffness of the methylated dinucleotide step changes marginally, whereas the neighboring steps become stiffer. Stiffening is further enhanced for consecutively methylated steps, providing a mechanistic origin for the effect of hypermethylation. Steric interactions between the added methyl groups and the nonpolar groups of the neighboring nucleotides are responsible for the stiffening in most cases. By constructing hydration maps, we found that methylation also alters the surface hydration structure in distinct ways. Its resistance to deformation may contribute to the stiffening of DNA for deformational modes lacking steric interactions. These results highlight the sequence- and deformational-mode-dependent effects of cytosine methylation.

Molecular Mechanisms of Tight Binding through Fuzzy Interactions.

Many intrinsically disordered proteins (IDPs) form fuzzy complexes upon binding to their targets. Although many IDPs are weakly bound in fuzzy complexes, some IDPs form high-affinity complexes. One example is the nonstructural protein 1 (NS1) of the 1918 Spanish influenza A virus, which hijacks cellular CRKII through the strong binding affinity (Kd ∼10 nM) of its proline-rich motif (PRMNS1) to the N-terminal Src-homology 3 domain of CRKII. However, its molecular mechanism remains elusive. Here, we examine the interplay between structural disorder of a bound PRMNS1 and its long-range electrostatic interactions. Using x-ray crystallography and NMR spectroscopy, we found that PRMNS1 retains substantial conformational flexibility in the bound state. Moreover, molecular dynamics simulations showed that structural disorder of the bound PRMNS1 increases the number of electrostatic interactions and decreases the mean distances between the positively charged residues in PRMNS1 and the acidic residues in the N-terminal Src-homology 3 domain. These results are analyzed using a polyelectrostatic model. Our results provide an insight into the molecular recognition mechanism for a high-affinity fuzzy complex.

Molecular motor-powered shuttles along multi-walled carbon nanotube tracks.

As a complementary tool to nanofluidics, biomolecular-based transport is envisioned for nanotechnological devices. We report a new method for guiding microtubule shuttles on multi-walled carbon nanotube tracks, aligned by dielectrophoresis on a functionalized surface. In the absence of electric field and in fluid flow, alignment is maintained. The directed translocation of kinesin propelled microtubules has been investigated using fluorescence microscopy. To our knowledge, this is the first demonstration of microtubules gliding along carbon nanotubes.

Kinetic signature of fractal-like filament networks formed by orientational linear epitaxy.

We study a broad class of epitaxial assembly of filament networks on lattice surfaces. Over time, a scale-free behavior emerges with a 2.5-3 power-law exponent in filament length distribution. Partitioning between the power-law and exponential behaviors in a network can be used to find the stage and kinetic parameters of the assembly process. To analyze real-world networks, we develop a computer program that measures the network architecture in experimental images. Application to triaxial networks of collagen fibrils shows quantitative agreement with our model. Our unifying approach can be used for characterizing and controlling the network formation that is observed across biological and nonbiological systems.

Microtubule shuttles on kinesin-coated glass micro-wire tracks.

Gliding of microtubule filaments on surfaces coated with the motor protein kinesin has potential applications for nano-scale devices. The ability to guide the gliding direction in three dimensions allows the fabrication of tracks of arbitrary geometry in space. Here, we achieve this by using kinesin-coated glass wires of micrometer diameter range. Unlike previous methods in which the guiding tracks are fixed on flat two-dimensional surfaces, the flexibility of glass wires in shape and size facilitates building in-vitro devices that have deformable tracks.

Chain registry and load-dependent conformational dynamics of collagen.

Degradation of fibrillar collagen is critical for tissue maintenance. Yet, understanding collagen catabolism has been challenging partly due to a lack of atomistic picture for its load-dependent conformational dynamics, as both mechanical load and local unfolding of collagen affect its cleavage by matrix metalloproteinase (MMP). We use molecular dynamics simulation to find the most cleavage-prone arrangement of α chains in a collagen triple helix and find amino acids that modulate stability of the MMP cleavage domain depending on the chain registry within the molecule. The native-like state is mechanically inhomogeneous, where the cleavage site interfaces a stiff region and a locally unfolded and flexible region along the molecule. In contrast, a triple helix made of the stable glycine-proline-hydroxyproline motif is uniformly flexible and is dynamically stabilized by short-lived, low-occupancy hydrogen bonds. These results provide an atomistic basis for the mechanics, conformation, and stability of collagen that affect catabolism.

Evolutionary rewiring of the dynamic network underpinning allosteric epistasis in NS1 of influenza A virus.

Viral proteins frequently mutate to evade or antagonize host innate immune responses, yet the impact of these mutations on the molecular energy landscape remains unclear. Epistasis, the intramolecular communications between mutations, often renders the combined mutational effects unpredictable. Nonstructural protein 1 (NS1) is a major virulence factor of the influenza A virus (IAV) that activates host PI3K by binding to its p85ß subunit. Here, we present the deep analysis for the impact of evolutionary mutations in NS1 that emerged between the 1918 pandemic IAV strain and its descendant PR8 strain. Our analysis reveal how the mutations rewired inter-residue communications which underlies long-range allosteric and epistatic networks in NS1. Our findings show that PR8 NS1 binds to p85ß with approximately 10-fold greater affinity than 1918 NS1 due to allosteric mutational effects. Notably, these mutations also exhibited long-range epistatic effects. NMR chemical shift perturbation and methyl-axis order parameter analyses revealed that the mutations induced long-range structural and dynamic changes in PR8 NS1, enhancing its affinity to p85ß. Complementary MD simulations and graph-based network analysis uncover how these mutations rewire dynamic residue interaction networks, which underlies the long-range epistasis and allosteric effects on p85ß-binding affinity. Significantly, we find that conformational dynamics of residues with high betweenness centrality play a crucial role in communications between network communities and are highly conserved across influenza A virus evolution. These findings advance our mechanistic understanding of the allosteric and epistatic communications between distant residues and provides insight into their role in the molecular evolution of NS1.

Load-based divergence in the dynamic allostery of two TCRs recognizing the same pMHC.

Increasing evidence suggests that mechanical load on the αß T cell receptor (TCR) is crucial for recognizing the antigenic peptide-loaded major histocompatibility complex (pMHC) molecule. Our recent all-atom molecular dynamics (MD) simulations revealed that the inter-domain motion of the TCR is responsible for the load-induced catch bond behavior of the TCR-pMHC complex and peptide discrimination. To further examine the generality of the mechanism, we perform all-atom MD simulations of the B7 TCR under different conditions for comparison with our previous simulations of the A6 TCR. The two TCRs recognize the same pMHC and have similar interfaces with pMHC in crystal structures. We find that the B7 TCR-pMHC interface stabilizes under ~15-pN load using a conserved dynamic allostery mechanism that involves the asymmetric motion of the TCR chassis. However, despite forming comparable contacts with pMHC as A6 in the crystal structure, B7 has fewer high-occupancy contacts with pMHC during the simulation. These results suggest that the dynamic allostery common to the TCRαß chassis can amplify slight differences in interfacial contacts into distinctive mechanical responses and potentially nuanced biological outcomes.

Asymmetric framework motion of TCRαß controls load-dependent peptide discrimination.

Mechanical force is critical for the interaction between an αß T cell receptor (TCR) and a peptide-bound major histocompatibility complex (pMHC) molecule to initiate productive T-cell activation. However, the underlying mechanism remains unclear. We use all-atom molecular dynamics simulations to examine the A6 TCR bound to HLA-A*02:01 presenting agonist or antagonist peptides under different extensions to simulate the effects of applied load on the complex, elucidating their divergent biological responses. We found that TCR α and ß chains move asymmetrically, which impacts the interface with pMHC, in particular the peptide-sensing CDR3 loops. For the wild-type agonist, the complex stabilizes in a load-dependent manner while antagonists destabilize it. Simulations of the Cß FG-loop deletion, which reduces the catch bond response, and simulations with in silico mutant peptides further support the observed behaviors. The present results highlight the combined role of interdomain motion, fluctuating forces, and interfacial contacts in determining the mechanical response and fine peptide discrimination by a TCR, thereby resolving the conundrum of nearly identical crystal structures of TCRαß-pMHC agonist and antagonist complexes.

Parsing digital or analog TCR performance through piconewton forces.

αß T cell receptors (TCRs) principally recognize aberrant peptides bound to major histocompatibility complex molecules (pMHCs) on unhealthy cells, amplifying specificity and sensitivity through physical load placed on the TCR-pMHC bond during immunosurveillance. To understand this mechanobiology, TCRs stimulated by abundantly and sparsely arrayed epitopes (NP366-374/Db and PA224-233/Db, respectively) following in vivo influenza A virus infection were studied with optical tweezers. While certain NP repertoire CD8 T lymphocytes require many ligands for activation, others are digital, needing just few. Conversely, all PA TCRs perform digitally, exhibiting pronounced bond lifetime increases through sustained, energizing volleys of structural transitioning. Optimal digital performance is superior in vivo, correlating with ERK phosphorylation, CD3 loss, and activation marker up-regulation in vitro. Given neoantigen array paucity, digital TCRs are likely critical for immunotherapies.

CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

Modular aspects of kinesin force generation machinery.

The motor head of kinesin carries out microtubule binding, ATP hydrolysis, and force generation. Despite a high level of sequence and structural conservation, subtle variations in subdomains of the motor head determine family-specific properties. In particular, both Kinesin-1 (Kin-1) and Kinesin-5 (Kin-5) walk processively to the microtubule plus-end, yet show distinct motility characteristics suitable for their functions. We studied chimeric Kin-1/Kin-5 constructs with a combination of single molecule motility assays and molecular dynamics simulations to demonstrate that Kin-5 possesses a force-generating element similar to Kin-1, i.e., the cover-neck bundle. Furthermore, the Kin-5 neck linker makes additional contacts with the core of the motor head via loop L13, which putatively compensates for the shorter cover-neck bundle of Kin-5. Our results indicate that Kin-1 is mechanically optimized for individual cargo transport, whereas Kin-5 does not necessarily maximize its mechanical performance. Its biochemical rates and enhanced force sensitivity may instead be beneficial for operation in a group of motors. Such variations in subdomains would be a strategy for achieving diversity in motility with the conserved motor head.