Membrane-mediated regulation of the intrinsically disordered CD3ϵ cytoplasmic tail of the TCR

Quantitative Ensemble Interpretation of Membrane Paramagnetic Relaxation Enhancement (mPRE) for Studying Membrane-Associated Intrinsically Disordered Proteins

An understanding of the functional role played by a membrane-associated intrinsically disordered protein (IDP) requires characterization of its heterogeneous conformations as well as its poses relative to the membranes, which is of great interest but technically challenging. Here, we explore the membrane paramagnetic relaxation enhancement (mPRE) for constructing ensembles of IDPs that dynamically associate with membrane mimetics incorporating spin-labeled lipids. To accurately interpret the mPRE Γ2 rates, both the dynamics of IDPs and spin probe molecules are taken into account, with the latter described by a weighted three-dimensional (3D) grid model built based on all-atom simulations. The IDP internal conformations, orientations, and immersion depths in lipid bilayers are comprehensively optimized in the Γ2-based ensemble modeling. Our approach is tested and validated on the example of POPG bicelle-bound disordered cytoplasmic domain of CD3ε (CD3εCD), a component of the T-cell receptor (TCR) complex. The mPRE-derived CD3εCD ensemble provides new insights into the IDP-membrane fuzzy association, in particular for the tyrosine-based signaling motif that plays a critical role in TCR signaling. The comparative analysis of the ensembles for wild-type CD3εCD and mutants that mimic the mono- and dual-phosphorylation effects suggests a delicate membrane regulatory mechanism for activation and inhibition of the TCR activity.

A dynamic biomimetic model of the membrane-bound CD4-CD3-TCR complex during pMHC disengagement

The coordinated (dis)engagement of the membrane-bound T cell receptor (TCR)-CD3-CD4 complex from the peptide-major histocompatibility complex (pMHC) is fundamental to TCR signal transduction and T cell effector function. As such, an atomic-scale understanding would not only enhance our basic understanding of the adaptive immune response but would also accelerate the rational design of TCRs for immunotherapy. In this study, we explore the impact of the CD4 coreceptor on the TCR-pMHC (dis)engagement by constructing a molecular-level biomimetic model of the CD3-TCR-pMHC and CD4-CD3-TCR-pMHC complexes within a lipid bilayer. After allowing the system complexes to equilibrate (engage), we use steered molecular dynamics to dissociate (disengage) the pMHC. We find that 1) the CD4 confines the pMHC closer to the T cell by 1.8 nm at equilibrium; 2) CD4 confinement shifts the TCR along the MHC binding groove engaging a different set of amino acids and enhancing the TCR-pMHC bond lifetime; 3) the CD4 translocates under load increasing the interaction strength between the CD4-pMHC, CD4-TCR, and CD4-CD3; and 4) upon dissociation, the CD3-TCR complex undergoes structural oscillation and increased energetic fluctuation between the CD3-TCR and CD3-lipids. These atomic-level simulations provide mechanistic insight on how the CD4 coreceptor impacts TCR-pMHC (dis)engagement. More specifically, our results provide further support (enhanced bond lifetime) for a force-dependent kinetic proofreading model and identify an alternate set of amino acids in the TCR that dominate the TCR-pMHC interaction and could thus impact the design of TCRs for immunotherapy.

Spatial requirements for ITAM signaling in an intracellular natural killer cell model membrane

FcγRIIIa-FcεRIγ complexes, upon stimulation by antibodies, cluster to initiate intracellular signaling and activate natural killer (NK) cells. Intracellular signaling involves Lck phosphorylation of ITAMs of each monomer of a FcεRIγ homodimer in a FcγRIIIa-FcεRIγ complex and subsequent binding of two phosphotyrosines (pY) in tandem by a Syk family kinase. However, how FcR clustering triggers ITAM signaling is not resolved. Molecular modeling and dynamics (MD) simulations are applied to generate ensembles of structures of the FcγRIIIa and FcεRIγ homodimeric cytoplasmic tails of FcγRIIIa-FcεRIγ complexes based on the transmembrane helices and cytoplasmic tails spaced 120, 80, and 50 Å apart to model different extents of clustering. Site-identification by ligand competitive saturation method with Monte Carlo sampling (SILCS-MC) is used to model how Lck could phosphorylate a diversity of ITAM conformations. At 80 Å separation between FcγRIIIa-FcεRIγ complexes, Lck can perform multiple phosphorylations on individual and multiple ITAMs across complexes, including potential sequential phosphorylation events. Syk may then potentially bind the two pYs within a single ITAM in tandem in isolated FcγRIIIa-FcεRIγ complexes, as observed in CD3ε and ζ chains of T cell receptors by the Syk family kinase ZAP-70. In addition, at 50 Å separation between complexes, unique to natural killer cells over T cells, Syk could potentially bind in tandem to pYs in different ITAMs across FcγRIIIa-FcεRIγ complexes. Thus, we predict that an ensemble of spatial orientations of the ITAMS of FcγRIIIa-FcεRIγ complexes that occur upon clustering lead to ITAM phosphorylation by Lck and subsequent Syk activity thereby facilitating downstream signaling.

Intrinsic Disorder in the T Cell Receptor Creates Cooperativity and Controls ZAP70 Binding

Many immunoreceptors have cytoplasmic domains that are intrinsically disordered (i.e., have high configurational entropy), have multiple sites of posttranslational modification (e.g., tyrosine phosphorylation), and participate in nonlinear signaling pathways (e.g., exhibiting switch-like behavior). Several hypotheses to explain the origin of these nonlinearities fall under the broad hypothesis that modification at one site changes the immunoreceptor’s entropy, which in turn changes further modification dynamics. Here, we use coarse-grain simulation to study three scenarios, all related to the chains that constitute the T cell receptor (TCR). We find that first, if phosphorylation induces local changes in the flexibility of the TCR ζ-chain, this naturally leads to rate enhancements and cooperativity. Second, we find that TCR CD3ɛ can provide a switch by modulating its residence in the plasma membrane. By constraining our model to be consistent with the previous observation that both basic residues and phosphorylation control membrane residence, we find that there is only a moderate rate enhancement of 10% between first and subsequent phosphorylation events. Third, we find that volume constraints do not limit the number of ZAP70s that can bind the TCR but that entropic penalties lead to a 200-fold decrease in binding rate by the seventh ZAP70, potentially explaining the observation that each TCR has around six ZAP70 molecules bound after receptor triggering. In all three scenarios, our results demonstrate that phenomena that change an immunoreceptor chain’s entropy (stiffening, confinement to a membrane, and multiple simultaneous binding) can lead to nonlinearities (rate enhancement, switching, and negative cooperativity) in how the receptor participates in signaling. These polymer-entropy-driven nonlinearities may augment the nonlinearities that arise from, e.g., kinetic proofreading and cluster formation. They also suggest different design strategies for engineered receptors, e.g., whether or not to put signaling modules on one chain or multiple clustered chains.

Short Disordered Epitope of CRTAM Ig-Like V Domain as a Potential Target for Blocking Antibodies

Class-I Restricted T Cell-Associated Molecule (CRTAM) is a protein that is expressed after T cell activation. The interaction of CRTAM with its ligand, nectin-like 2 (Necl2), is required for the efficient production of IL-17, IL-22, and IFNγ by murine CD4 T cells, and it plays a role in optimal CD8 T and NK cell cytotoxicity. CRTAM promotes the pro-inflammatory cytokine profile; therefore, it may take part in the immunopathology of autoimmune diseases such as diabetes type 1 or colitis. Thus, antibodies that block the interaction between CRTAM and Necl2 would be useful for controlling the production of these inflammatory cytokines. In this work, using bioinformatics predictions, we identified three short disordered epitopes (sDE1-3) that are located in the Ig-like domains of murine CRTAM and are conserved in mammalian species. We performed a structural analysis by molecular dynamics simulations of sDE1 (QHPALKSSKY, Ig-like V), sDE2 (QRNGEKSVVK, Ig-like C1), and sDE3 (CSTERSKKPPPQI, Ig-like C1). sDE1, which is located within a loop of the contact interface of the heterotypic interaction with Nectl2, undergoes an order–disorder transition. On the contrary, even though sDE2 and sDE3 are flexible and also located within loops, they do not undergo order–disorder transitions. We evaluated the immunogenicity of sDE1 and sDE3 through the expression of these epitopes in chimeric L1 virus-like particles. We confirmed that sDE1 induces polyclonal antibodies that recognize the native folding of CRTAM expressed in activated murine CD4 T cells. In contrast, sDE3 induces polyclonal antibodies that recognize the recombinant protein hCRTAM-Fc, but not the native CRTAM. Thus, in this study, an exposed disordered epitope in the Ig-like V domain of CRTAM was identified as a potential site for therapeutic antibodies.

The TCR is an allosterically regulated macromolecular machinery changing its conformation while working

The αβ T-cell receptor (TCR) is a multiprotein complex controlling the activation of T cells. Although the structure of the complete TCR is not known, cumulative evidence supports that the TCR cycles between different conformational states that are promoted either by thermal motion or by force. These structural transitions determine whether the TCR engages intracellular effectors or not, regulating TCR phosphorylation and signaling. As for other membrane receptors, ligand binding selects and stabilizes the TCR in active conformations, and/or switches the TCR to activating states that were not visited before ligand engagement. Here we review the main models of TCR allostery, that is, ligand binding at TCRαβ changes the structure at CD3 and ζ. (a) The ITAM and proline-rich sequence exposure model, in which the TCR's cytoplasmic tails shield each other and ligand binding exposes them for phosphorylation. (b) The membrane-ITAM model, in which the CD3ε and ζ tails are sequestered inside the membrane and again ligand binding exposes them. (c) The mechanosensor model in which ligand binding exerts force on the TCR, inducing structural changes that allow signaling. Since these models are complementary rather than competing, we propose a unified model that aims to incorporate all existing data.

Molecular Dynamics of the Recruitment of Immunoreceptor Signaling Module DAP12 Homodimer to Lipid Raft Boundary Regulated by PIP2

Lipid raft microdomain of the plasma membrane is implicated in various biological and pathological processes. The involvement of lipid raft in T cell receptor (TCR) signal transduction has been widely studied, whereas the role of these structures in immunoreceptor signaling by DAP12 in natural killing (NK) cells remains largely unknown. Here, we demonstrate that phosphatidylinositol 4,5-bisphosphate (PIP2) lipid localized to lipid raft boundary in our coarse-grained (CG) model raft-forming membrane, and this negatively charged lipid recruits DAP12 homodimer into lipid raft boundary through protein-lipid interaction between the basic-rich regions and signaling immunoreceptor tyrosine-based activation motifs (ITAMs) of DAP12 and PIP2. Furthermore, our results reveal that the protein-lipid interaction can be disrupted by Ca²⁺, which competitively binds to PIP2 instead of DAP12. As a result, the cytoplasmic region of DAP12 homodimer is dissociated from the membrane back to the nonraft domain, and the ITAMs are exposed to allow further downstream signaling. These findings provide fundamental insights to understand the mechanism of signal transduction in NK cells regulated by membrane microenvironment.

The structural basis for membrane assembly of immunoreceptor signalling complexes

Immunoreceptors are TM complexes that consist of separate ligand-binding and signal-transducing modules. Mounting evidence suggests that interactions with the local environment may influence the architecture of these TM domains, which assemble via crucial sets of conserved ionisable residues, and also control the peripheral association of immunoreceptor tyrosine-based activation motifs (ITAMs) whose phosphorylation triggers cytoplasmic signalling cascades. We now report a molecular dynamics (MD) simulation study of the archetypal T cell receptor (TCR) and its cluster of differentiation 3 (CD3) signalling partners, along with the analogous DNAX-activation protein of 12 kDa (DAP12)/natural killer group 2C (NKG2C) complex. Based on > 15 μs of explicitly solvated, atomic-resolution sampling, we explore molecular aspects of immunoreceptor complex stability in different functionally relevant states. A novel alchemical approach is used to simulate the cytoplasmic CD3ε tail at different depths within lipid bilayer models, revealing that the conformation and cytoplasmic exposure of ITAMs are highly sensitive to local enrichment by different lipid species and to phosphorylation. Furthermore, simulations of the TCR and DAP12 TM domains in various states of oligomerisation suggest that, during the early stages of assembly, stable membrane insertion is facilitated by the interfacial lipid/solvent environment and/or partial ionisation of charged residues. Collectively, our results indicate that the architecture and mechanisms of signal transduction in immunoreceptor complexes are tightly regulated by interactions with the microenvironment.

Combinatorial diversity of Syk recruitment driven by its multivalent engagement with FcεRIγ

Syk/Zap70 family kinases are essential for signaling via multichain immune-recognition receptors such as tetrameric (αβγ2) FcεRI. Syk activation is generally attributed to cis binding of its tandem SH2 domains to dual phosphotyrosines within FcεRIγ-ITAMs (immunoreceptor tyrosine-based activation motifs). However, the mechanistic details of Syk docking on γ homodimers are unresolved. Here, we estimate that multivalent interactions for WT Syk improve cis-oriented binding by three orders of magnitude. We applied molecular dynamics (MD), hybrid MD/worm-like chain polymer modeling, and live cell imaging to evaluate relative binding and signaling output for all possible cis and trans Syk-FcεRIγ configurations. Syk binding is likely modulated during signaling by autophosphorylation on Y130 in interdomain A, since a Y130E phosphomimetic form of Syk is predicted to lead to reduced helicity of interdomain A and alter Syk's bias for cis binding. Experiments in reconstituted γ-KO cells, whose γ subunits are linked by disulfide bonds, as well as in cells expressing monomeric ITAM or hemITAM γ-chimeras, support model predictions that short distances between γ ITAM pairs are required for trans docking. We propose that the full range of docking configurations improves signaling efficiency by expanding the combinatorial possibilities for Syk recruitment, particularly under conditions of incomplete ITAM phosphorylation.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

The Neutral Sphingomyelinase 2 Is Required to Polarize and Sustain T Cell Receptor Signaling

By promoting ceramide release at the cytosolic membrane leaflet, the neutral sphingomyelinase 2 (NSM) is capable of organizing receptor and signalosome segregation. Its role in T cell receptor (TCR) signaling remained so far unknown. We now show that TCR-driven NSM activation is dispensable for TCR clustering and initial phosphorylation, but of crucial importance for further signal amplification. In particular, at low doses of TCR stimulatory antibodies, NSM is required for Ca²⁺ mobilization and T cell proliferation. NSM-deficient T cells lack sustained CD3ζ and ZAP-70 phosphorylation and are unable to polarize and stabilize their microtubular system. We identified PKCζ as the key NSM downstream effector in this second wave of TCR signaling supporting dynamics of microtubule-organizing center (MTOC). Ceramide supplementation rescued PKCζ membrane recruitment and MTOC translocation in NSM-deficient cells. These findings identify the NSM as essential in TCR signaling when dynamic cytoskeletal reorganization promotes continued lateral and vertical supply of TCR signaling components: CD3ζ, Zap70, and PKCζ, and functional immune synapses are organized and stabilized via MTOC polarization.

Membrane perturbing properties of toxin mycolactone from Mycobacterium ulcerans

Mycolactone is the exotoxin produced by Mycobacterium ulcerans and is the virulence factor behind the neglected tropical disease Buruli ulcer. The toxin has a broad spectrum of biological effects within the host organism, stemming from its interaction with at least two molecular targets and the inhibition of protein uptake into the endoplasmic reticulum. Although it has been shown that the toxin can passively permeate into host cells, it is clearly lipophilic. Association with lipid carriers would have substantial implications for the toxin’s distribution within a host organism, delivery to cellular targets, diagnostic susceptibility, and mechanisms of pathogenicity. Yet the toxin’s interactions with, and distribution in, lipids are unknown. Herein we have used coarse-grained molecular dynamics simulations, guided by all-atom simulations, to study the interaction of mycolactone with pure and mixed lipid membranes. Using established techniques, we calculated the toxin’s preferential localization, membrane translocation, and impact on membrane physical and dynamical properties. The computed water-octanol partition coefficient indicates that mycolactone prefers to be in an organic phase rather than in an aqueous environment. Our results show that in a solvated membrane environment the exotoxin mainly localizes in the water-membrane interface, with a preference for the glycerol moiety of lipids, consistent with the reported studies that found it in lipid extracts of the cell. The calculated association constant to the model membrane is similar to the reported association constant for Wiskott-Aldrich syndrome protein. Mycolactone is shown to modify the physical properties of membranes, lowering the transition temperature, compressibility modulus, and critical line tension at which pores can be stabilized. It also shows a tendency to behave as a linactant, a molecule that localizes at the boundary between different fluid lipid domains in membranes and promotes inter-mixing of domains. This property has implications for the toxin’s cellular access, T-cell immunosuppression, and therapeutic potential.

Mycolactone is a macrolide exotoxin secreted by Mycobacterium ulcerans, which causes a skin disease called Buruli ulcer, a neglected emerging disease. It is the third most common mycobacterial disease after tuberculosis and leprosy. Studies have shown how mycolactone plays a pivotal role in Buruli ulcer pathogenesis, and identified it as an attractive therapeutic target. This multifunctional cytotoxin exerts multiple local and global responses, including ulcerative, analgesic, and anti-inflammatory effects. Prompted by its lipid-like structure, we used extensive multi-resolution simulations to probe mycolactone’s interactions with model membranes. Our results suggest that mycolactone is sequestered in membranes where it alters several dynamical, physical, and mechanical properties. It also behaves as a linactant, localizing at the interface between lipid domains and decreasing the inter-domain line tension. Our results shed light on how mycolactone permeates host cell membranes and is distributed between lipid and aqueous environments. These findings have significant implications for the toxin’s distribution in the host environment and mechanisms of pathogenicity. Understanding the toxin’s distribution and mechanism of trafficking will have ramifications for targeted diagnostics, therapeutic approaches, and our understanding of Buruli ulcer pathogenesis.

Functions of intrinsic disorder in transmembrane proteins

Intrinsic disorder is common in integral membrane proteins, particularly in the intracellular domains. Despite this observation, these domains are not always recognized as being disordered. In this review, we will discuss the biological functions of intrinsically disordered regions of membrane proteins, and address why the flexibility afforded by disorder is mechanistically important. Intrinsically disordered regions are present in many common classes of membrane proteins including ion channels and transporters; G-protein coupled receptors (GPCRs), receptor tyrosine kinases and cytokine receptors. The functions of the disordered regions are many and varied. We will discuss selected examples including: (1) Organization of receptors, kinases, phosphatases and second messenger sources into signaling complexes. (2) Modulation of the membrane-embedded domain function by ball-and-chain like mechanisms. (3) Trafficking of membrane proteins. (4) Transient membrane associations. (5) Post-translational modifications most notably phosphorylation and (6) disorder-linked isoform dependent function. We finish the review by discussing the future challenges facing the membrane protein community regarding protein disorder.

Differential mast cell outcomes are sensitive to FcεRI-Syk binding kinetics

Cross-linking of immunoglobulin E-bound FcεRI triggers multiple cellular responses, including degranulation and cytokine production. Signaling is dependent on recruitment of Syk via docking of its dual SH2 domains to phosphorylated tyrosines within the FcεRI immunoreceptor tyrosine-based activation motifs. Using single-molecule imaging in live cells, we directly visualized and quantified the binding of individual mNeonGreen-tagged Syk molecules as they associated with the plasma membrane after FcεRI activation. We found that Syk colocalizes transiently to FcεRI and that Syk-FcεRI binding dynamics are independent of receptor aggregate size. Substitution of glutamic acid for tyrosine between the Syk SH2 domains (Syk-Y130E) led to an increased Syk-FcεRI off-rate, loss of site-specific Syk autophosphorylation, and impaired downstream signaling. Genome edited cells expressing only Syk-Y130E were deficient in antigen-stimulated calcium release, degranulation, and production of some cytokines (TNF-a, IL-3) but not others (MCP-1, IL-4). We propose that kinetic discrimination along the FcεRI signaling pathway occurs at the level of Syk-FcεRI interactions, with key outcomes dependent upon sufficiently long-lived Syk binding events.