A molecular switch in RCK2 triggers sodium-dependent activation of KNa1.1 (KCNT1) potassium channels.

The Na+-activated K+ channel KNa1.1, encoded by the KCNT1 gene, is an important regulator of neuronal excitability. How intracellular Na+ ions bind and increase channel activity is not well understood. Analysis of KNa1.1 channel structures indicate that there is a large twisting of the ßN-αQ loop in the intracellular RCK2 domain between the inactive and Na+-activated conformations, with a lysine (K885, human subunit numbering) close enough to potentially form a salt bridge with an aspartate (D839) in ßL in the Na+-activated state. Concurrently, an aspartate (D884) adjacent in the same loop adopts a position within a pocket formed by the ßO strand. In carrying out mutagenesis and electrophysiology with human KNa1.1, we found that alanine substitution of selected residues in these regions resulted in almost negligible currents in the presence of up to 40 mM intracellular Na+. The exception was D884A, which resulted in constitutively active channels in both the presence and absence of intracellular Na+. Further mutagenesis of this site revealed an amino acid size-dependent effect. Substitutions at this site by an amino acid smaller than aspartate (D884V) also yielded constitutively active KNa1.1, and D884I had Na+ dependence similar to wild-type KNa1.1, while increasing the side-chain size larger than aspartate (D884E or D884F) yielded channels that could not be activated by up to 40 mM intracellular Na+. We conclude that Na+ binding results in a conformational change that accommodates D884 in the ßO pocket, which triggers further conformational changes in the RCK domains and channel activation.

mPPases create a conserved anionic membrane fingerprint as identified via multi-scale simulations.

Membrane-integral pyrophosphatases (mPPases) are membrane-bound enzymes responsible for hydrolysing inorganic pyrophosphate and translocating a cation across the membrane. Their function is essential for the infectivity of clinically relevant protozoan parasites and plant maturation. Recent developments have indicated that their mechanism is more complicated than previously thought and that the membrane environment may be important for their function. In this work, we use multiscale molecular dynamics simulations to demonstrate for the first time that mPPases form specific anionic lipid interactions at 4 sites at the distal and interfacial regions of the protein. These interactions are conserved in simulations of the mPPases from Thermotoga maritima, Vigna radiata and Clostridium leptum and characterised by interactions with positive residues on helices 1, 2, 3 and 4 for the distal site, or 9, 10, 13 and 14 for the interfacial site. Due to the importance of these helices in protein stability and function, these lipid interactions may play a crucial role in the mPPase mechanism and enable future structural and functional studies.

Multi-scale simulations of the T cell receptor reveal its lipid interactions, dynamics and the arrangement of its cytoplasmic region.

The T cell receptor (TCR-CD3) initiates T cell activation by binding to peptides of Major Histocompatibility Complexes (pMHC). The TCR-CD3 topology is well understood but the arrangement and dynamics of its cytoplasmic tails remains unknown, limiting our grasp of the signalling mechanism. Here, we use molecular dynamics simulations and modelling to investigate the entire TCR-CD3 embedded in a model membrane. Our study demonstrates conformational changes in the extracellular and transmembrane domains, and the arrangement of the TCR-CD3 cytoplasmic tails. The cytoplasmic tails formed highly interlaced structures while some tyrosines within the immunoreceptor tyrosine-based activation motifs (ITAMs) penetrated the hydrophobic core of the membrane. Interactions between the cytoplasmic tails and phosphatidylinositol phosphate lipids in the inner membrane leaflet led to the formation of a distinct anionic lipid fingerprint around the TCR-CD3. These results increase our understanding of the TCR-CD3 dynamics and the importance of membrane lipids in regulating T cell activation.

Molecular dynamics simulations of Piezo1 channel opening by increases in membrane tension.

Piezo1 is a mechanosensitive channel involved in many cellular functions and responsible for sensing shear stress and pressure forces in cells. Piezo1 has a unique trilobed topology with a curved membrane region in the closed state. It has been suggested that upon activation Piezo1 adopts a flattened conformation, but the molecular and structural changes underpinning the Piezo1 gating and opening mechanisms and how the channel senses forces in the membrane remain elusive. Here, we used molecular dynamics simulations to reveal the structural rearrangements that occur when Piezo1 moves from a closed to an open state in response to increased mechanical tension applied to a model membrane. We find that membrane stretching causes Piezo1 to flatten and expand its blade region, resulting in tilting and lateral movement of the pore-lining transmembrane helices 37 and 38. This is associated with the opening of the channel and movement of lipids out of the pore region. Our results reveal that because of the rather loose packing of Piezo1 pore region, movement of the lipids outside the pore region is critical for the opening of the pore. Our simulations also suggest synchronous flattening of the Piezo1 blades during Piezo1 activation. The flattened structure lifts the C-terminal extracellular domain up, exposing it more to the extracellular space. Our studies support the idea that it is the blade region of Piezo1 that senses tension in the membrane because pore opening failed in the absence of the blades. Additionally, our simulations reveal that upon opening, water molecules occupy lateral fenestrations in the cytosolic region of Piezo1, which might be likely paths for ion permeation. Our results provide a model for how mechanical force opens the Piezo1 channel and thus how it might couple mechanical force to biological response.

Modeling of full-length Piezo1 suggests importance of the proximal N-terminus for dome structure.

Piezo1 forms a mechanically activated calcium-permeable nonselective cation channel that is functionally important in many cell types. Structural data exist for C-terminal regions, but we lack information about N-terminal regions and how the entire channel interacts with the lipid bilayer. Here, we use computational approaches to predict the three-dimensional structure of the full-length Piezo1 and simulate it in an asymmetric membrane. A number of novel insights are suggested by the model: 1) Piezo1 creates a trilobed dome in the membrane that extends beyond the radius of the protein, 2) Piezo1 changes the lipid environment in its vicinity via preferential interactions with cholesterol and phosphatidylinositol 4,5-bisphosphate (PIP2) molecules, and 3) cholesterol changes the depth of the dome and PIP2 binding preference. In vitro alteration of cholesterol concentration inhibits Piezo1 activity in a manner complementing some of our computational findings. The data suggest the importance of N-terminal regions of Piezo1 for dome structure and membrane cholesterol and PIP2 interactions.

ß1 integrin is a sensor of blood flow direction.

Endothelial cell (EC) sensing of fluid shear stress direction is a critical determinant of vascular health and disease. Unidirectional flow induces EC alignment and vascular homeostasis, whereas bidirectional flow has pathophysiological effects. ECs express several mechanoreceptors that respond to flow, but the mechanism for sensing shear stress direction is poorly understood. We determined, by using in vitro flow systems and magnetic tweezers, that ß1 integrin is a key sensor of force direction because it is activated by unidirectional, but not bidirectional, shearing forces. ß1 integrin activation by unidirectional force was amplified in ECs that were pre-sheared in the same direction, indicating that alignment and ß1 integrin activity has a feedforward interaction, which is a hallmark of system stability. En face staining and EC-specific genetic deletion studies in the murine aorta revealed that ß1 integrin is activated and is essential for EC alignment at sites of unidirectional flow but is not activated at sites of bidirectional flow. In summary, ß1 integrin sensing of unidirectional force is a key mechanism for decoding blood flow mechanics to promote vascular homeostasis.This article has an associated First Person interview with the first author of the paper.

Force Sensing by Piezo Channels in Cardiovascular Health and Disease.

Mechanical forces are fundamental in cardiovascular biology, and deciphering the mechanisms by which they act remains a testing frontier in cardiovascular research. Here, we raise awareness of 2 recently discovered proteins, Piezo1 and Piezo2, which assemble as transmembrane triskelions to combine exquisite force sensing with regulated calcium influx. There is emerging evidence for their importance in endothelial shear stress sensing and secretion, NO generation, vascular tone, angiogenesis, atherosclerosis, vascular permeability and remodeling, blood pressure regulation, insulin sensitivity, exercise performance, and baroreceptor reflex, and there are early suggestions of relevance to cardiac fibroblasts and myocytes. Human genetic analysis points to significance in lymphatic disease, anemia, varicose veins, and potentially heart failure, hypertension, aneurysms, and stroke. These channels appear to be versatile force sensors, used creatively to inform various force-sensing situations. We discuss emergent concepts and controversies and suggest that the potential for new important understanding is substantial.

Molecular Simulations of Intact Anion Exchanger 1 Reveal Specific Domain and Lipid Interactions.

Anion exchanger 1 (AE1) is responsible for the exchange of bicarbonate and chloride across the erythrocyte plasma membrane. Human AE1 consists of a cytoplasmic and a membrane domain joined by a 33-residue flexible linker. Crystal structures of the individual domains have been determined, but the intact AE1 structure remains elusive. In this study, we use molecular dynamics simulations and modeling to build intact AE1 structures in a complex lipid bilayer that resembles the native erythrocyte plasma membrane. AE1 models were evaluated using available experimental data to provide an atomistic view of the interaction and dynamics of the cytoplasmic domain, the membrane domain, and the connecting linker in a complete model of AE1 in a lipid bilayer. Anionic lipids were found to interact strongly with AE1 at specific amino acid residues that are linked to diseases and blood group antigens. Cholesterol was found in the dimeric interface of AE1, suggesting that it may regulate subunit interactions and anion transport.

Distinctive phosphoinositide- and Ca2+-binding properties of normal and cognitive performance-linked variant forms of KIBRA C2 domain.

Kidney- and brain-expressed protein (KIBRA), a multifunctional scaffold protein with around 20 known binding partners, is involved in memory and cognition, organ size control via the Hippo pathway, cell polarity, and membrane trafficking. KIBRA includes tandem N-terminal WW domains, a C2 domain, and motifs for binding atypical PKC and PDZ domains. A naturally occurring human KIBRA variant involving residue changes at positions 734 (Met-to-Ile) and 735 (Ser-to-Ala) within the C2 domain affects cognitive performance. We have elucidated 3D structures and calcium- and phosphoinositide-binding properties of human KIBRA C2 domain. Both WT and variant C2 adopt a canonical type I topology C2 domain fold. Neither Ca2+ nor any other metal ion was bound to WT or variant KIBRA C2 in crystal structures, and Ca2+ titration produced no significant reproducible changes in NMR spectra. NMR and X-ray diffraction data indicate that KIBRA C2 binds phosphoinositides via an atypical site involving ß-strands 5, 2, 1, and 8. Molecular dynamics simulations indicate that KIBRA C2 interacts with membranes via primary and secondary sites on the same domain face as the experimentally identified phosphoinositide-binding site. Our results indicate that KIBRA C2 domain association with membranes is calcium-independent and involves distinctive C2 domain-membrane relative orientations.

Interaction of the human erythrocyte Band 3 anion exchanger 1 (AE1, SLC4A1) with lipids and glycophorin A: Molecular organization of the Wright (Wr) blood group antigen.

The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5-bisphosphate (PIP2) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP2 molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers.

Structure and lipid-binding properties of the kindlin-3 pleckstrin homology domain.

Kindlins co-activate integrins alongside talin. They possess, like talin, a FERM domain (4.1-erythrin-radixin-moiesin domain) comprising F0-F3 subdomains, but with a pleckstrin homology (PH) domain inserted in the F2 subdomain that enables membrane association. We present the crystal structure of murine kindlin-3 PH domain determined at a resolution of 2.23 Šand characterise its lipid binding using biophysical and computational approaches. Molecular dynamics simulations suggest flexibility in the PH domain loops connecting ß-strands forming the putative phosphatidylinositol phosphate (PtdInsP)-binding site. Simulations with PtdInsP-containing bilayers reveal that the PH domain associates with PtdInsP molecules mainly via the positively charged surface presented by the ß1-ß2 loop and that it binds with somewhat higher affinity to PtdIns(3,4,5)P3 compared with PtdIns(4,5)P2 Surface plasmon resonance (SPR) with lipid headgroups immobilised and the PH domain as an analyte indicate affinities of 300 µM for PtdIns(3,4,5)P3 and 1 mM for PtdIns(4,5)P2 In contrast, SPR studies with an immobilised PH domain and lipid nanodiscs as the analyte show affinities of 0.40 µM for PtdIns(3,4,5)P3 and no affinity for PtdIns(4,5)P2 when the inositol phosphate constitutes 5% of the total lipids (∼5 molecules per nanodisc). Reducing the PtdIns(3,4,5)P3 composition to 1% abolishes nanodisc binding to the PH domain, as does site-directed mutagenesis of two lysines within the ß1-ß2 loop. Binding of PtdIns(3,4,5)P3 by a canonical PH domain, Grp1, is not similarly influenced by SPR experimental design. These data suggest a role for PtdIns(3,4,5)P3 clustering in the binding of some PH domains and not others, highlighting the importance of lipid mobility and clustering for the biophysical assessment of protein-membrane interactions.

Band 3, the human red cell chloride/bicarbonate anion exchanger (AE1, SLC4A1), in a structural context.

The crystal structure of the dimeric membrane domain of human Band 3(1), the red cell chloride/bicarbonate anion exchanger 1 (AE1, SLC4A1), provides a structural context for over four decades of studies into this historic and important membrane glycoprotein. In this review, we highlight the key structural features responsible for anion binding and translocation and have integrated the following topological markers within the Band 3 structure: blood group antigens, N-glycosylation site, protease cleavage sites, inhibitor and chemical labeling sites, and the results of scanning cysteine and N-glycosylation mutagenesis. Locations of mutations linked to human disease, including those responsible for Southeast Asian ovalocytosis, hereditary stomatocytosis, hereditary spherocytosis, and distal renal tubular acidosis, provide molecular insights into their effect on Band 3 folding. Finally, molecular dynamics simulations of phosphatidylcholine self-assembled around Band 3 provide a view of this membrane protein within a lipid bilayer.

The Integrin Receptor in Biologically Relevant Bilayers: Insights from Molecular Dynamics Simulations.

Integrins are heterodimeric (αß) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2-F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2-F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes.

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes.

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the ß1/ß2, ß3/ß4, and ß5/ß6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.

Molecular dynamics simulations of the bacterial UraA H+-uracil symporter in lipid bilayers reveal a closed state and a selective interaction with cardiolipin.

The Escherichia coli UraA H+-uracil symporter is a member of the nucleobase/ascorbate transporter (NAT) family of proteins, and is responsible for the proton-driven uptake of uracil. Multiscale molecular dynamics simulations of the UraA symporter in phospholipid bilayers consisting of: 1) 1-palmitoyl 2-oleoyl-phosphatidylcholine (POPC); 2) 1-palmitoyl 2-oleoyl-phosphatidylethanolamine (POPE); and 3) a mixture of 75% POPE, 20% 1-palmitoyl 2-oleoyl-phosphatidylglycerol (POPG); and 5% 1-palmitoyl 2-oleoyl-diphosphatidylglycerol/cardiolipin (CL) to mimic the lipid composition of the bacterial inner membrane, were performed using the MARTINI coarse-grained force field to self-assemble lipids around the crystal structure of this membrane transport protein, followed by atomistic simulations. The overall fold of the protein in lipid bilayers remained similar to the crystal structure in detergent on the timescale of our simulations. Simulations were performed in the absence of uracil, and resulted in a closed state of the transporter, due to relative movement of the gate and core domains. Anionic lipids, including POPG and especially CL, were found to associate with UraA, involving interactions between specific basic residues in loop regions and phosphate oxygens of the CL head group. In particular, three CL binding sites were identified on UraA: two in the inner leaflet and a single site in the outer leaflet. Mutation of basic residues in the binding sites resulted in the loss of CL binding in the simulations. CL may play a role as a "proton trap" that channels protons to and from this transporter within CL-enriched areas of the inner bacterial membrane.

Interactions of phosphatase and tensin homologue (PTEN) proteins with phosphatidylinositol phosphates: insights from molecular dynamics simulations of PTEN and voltage sensitive phosphatase.

The phosphatase and tensin homologue (PTEN) and the Ciona intestinalis voltage sensitive phosphatase (Ci-VSP) are both phosphatidylinositol phosphate (PIP) phosphatases that contain a C2 domain. PTEN is a tumor suppressor protein that acts as a phosphatase on PIP3 in mammalian cell membranes. It contains two principal domains: a phosphatase domain (PD) and a C2 domain. Despite detailed structural and functional characterization, less is known about its mechanism of interaction with PIP-containing lipid bilayers. Ci-VSP consists of an N-terminal transmembrane voltage sensor domain and a C-terminal PTEN domain, which in turn contains a PD and a C2 domain. The nature of the interaction of the PTEN domain of Ci-VSP with membranes has not been well established. We have used multiscale molecular dynamics simulations to define the interaction mechanisms of PTEN and of the Ci-VSP PTEN domains with PIP-containing lipid bilayers. Our results suggest a novel mechanism of association of the PTEN with such bilayers, in which an initial electrostatics-driven encounter of the protein and bilayer is followed by reorientation of the protein to optimize its interactions with PIP molecules in the membrane. Although a PIP3 molecule binds close to the active site of PTEN, our simulations suggest a further conformational change of the protein may be required for catalytically productive binding to occur. Ci-VSP interacted with membranes in an orientation comparable to that of PTEN but bound directly to PIP-containing membranes without a subsequent reorientation step. Again, PIP3 bound close to the active site of the Ci-VSP PD, but not in a catalytically productive manner. Interactions of Ci-VSP with the bilayer induced clustering of PIP molecules around the protein.

Interactions of peripheral proteins with model membranes as viewed by molecular dynamics simulations.

Many cellular signalling and related events are triggered by the association of peripheral proteins with anionic lipids in the cell membrane (e.g. phosphatidylinositol phosphates or PIPs). This association frequently occurs via lipid-binding modules, e.g. pleckstrin homology (PH), C2 and four-point-one, ezrin, radixin, moesin (FERM) domains, present in peripheral and cytosolic proteins. Multiscale simulation approaches that combine coarse-grained and atomistic MD simulations may now be applied with confidence to investigate the molecular mechanisms of the association of peripheral proteins with model bilayers. Comparisons with experimental data indicate that such simulations can predict specific peripheral protein-lipid interactions. We discuss the application of multiscale MD simulation and related approaches to investigate the association of peripheral proteins which contain PH, C2 or FERM-binding modules with lipid bilayers of differing phospholipid composition, including bilayers containing multiple PIP molecules.

Conformational changes in talin on binding to anionic phospholipid membranes facilitate signaling by integrin transmembrane helices.

Integrins are heterodimeric (αß) cell surface receptors that are activated to a high affinity state by the formation of a complex involving the α/ß integrin transmembrane helix dimer, the head domain of talin (a cytoplasmic protein that links integrins to actin), and the membrane. The talin head domain contains four sub-domains (F0, F1, F2 and F3) with a long cationic loop inserted in the F1 domain. Here, we model the binding and interactions of the complete talin head domain with a phospholipid bilayer, using multiscale molecular dynamics simulations. The role of the inserted F1 loop, which is missing from the crystal structure of the talin head, PDB:3IVF, is explored. The results show that the talin head domain binds to the membrane predominantly via cationic regions on the F2 and F3 subdomains and the F1 loop. Upon binding, the intact talin head adopts a novel V-shaped conformation which optimizes its interactions with the membrane. Simulations of the complex of talin with the integrin α/ß TM helix dimer in a membrane, show how this complex promotes a rearrangement, and eventual dissociation of, the integrin α and ß transmembrane helices. A model for the talin-mediated integrin activation is proposed which describes how the mutual interplay of interactions between transmembrane helices, the cytoplasmic talin protein, and the lipid bilayer promotes integrin inside-out activation.

Multiscale simulations suggest a mechanism for integrin inside-out activation.

Integrins are large cell-surface adhesion receptors that can be activated to a high affinity state by the formation of an intracellular complex between the integrin ß-subunit tail, the membrane, and talin. The F2 and F3 subdomains of the talin head play a key role in formation of this complex. Here, activation of the integrin αIIb/ß3 dimer by the talin head domain was probed using multiscale molecular dynamics simulations. A number of novel insights emerge from these studies, including (i) the importance of the integrin αIIb subunit F992 and F993 residues in stabilizing the "off" state of the αIIb/ß3 dimer, (ii) a crucial role for negatively charged groups in the F2-F3/membrane interaction, (iii) binding of the talin F2-F3 domain to negatively charged lipid headgroups in the membrane induces a reorientation of the ß transmembrane (TM) domain, (iv) an increase in the tilt angle of the ß TM domain relative to the bilayer normal helps to destabilize the α/ß TM interaction and promote a scissor-like movement of the integrin TM helices. These results, combined with various published experimental observations, suggest a model for the mechanism of inside-out activation of integrins by talin.