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.

Protein-lipid charge interactions control the folding of outer membrane proteins into asymmetric membranes

Biological membranes consist of two leaflets of phospholipid molecules that form a bilayer, each leaflet comprising a distinct lipid composition. This asymmetry is created and maintained in vivo by dedicated biochemical pathways, but difficulties in creating stable asymmetric membranes in vitro have restricted our understanding of how bilayer asymmetry modulates the folding, stability and function of membrane proteins. In this study, we used cyclodextrin-mediated lipid exchange to generate liposomes with asymmetric bilayers and characterize the stability and folding kinetics of two bacterial outer membrane proteins (OMPs), OmpA and BamA. We found that excess negative charge in the outer leaflet of a liposome impedes their insertion and folding, while excess negative charge in the inner leaflet accelerates their folding relative to symmetric liposomes with the same membrane composition. Using molecular dynamics, mutational analysis and bioinformatics, we identified a positively charged patch critical for folding and stability. These results rationalize the well-known 'positive-outside' rule of OMPs and suggest insights into the mechanisms that drive OMP folding and assembly in vitro and in vivo.

Molecular dynamics simulations reveal membrane lipid interactions of the full-length lymphocyte specific kinase (Lck)

The membrane-bound lymphocyte-specific protein-tyrosine kinase (Lck) triggers T cell antigen receptor signalling to initiate adaptive immune responses. Despite many structure-function studies, the mode of action of Lck and the potential role of plasma membrane lipids in regulating Lck's activity remains elusive. Advances in molecular dynamics simulations of membrane proteins in complex lipid bilayers have opened a new perspective in gathering such information. Here, we have modelled the full-length Lck open and closed conformations  using data available from different crystalographic studies and simulated its interaction with the inner leaflet of the T cell plasma membrane. In both conformations, we found that the unstructured unique domain and the structured domains including the kinase interacted with the membrane with a preference for PIP lipids. Interestingly, our simulations suggest that the Lck-SH2 domain interacts with lipids differently in the open and closed Lck conformations, demonstrating that lipid interaction can potentially regulate Lck's conformation and in turn modulate T cell signalling. Additionally, the Lck-SH2 and kinase domain residues that significantly contacted PIP lipids are found to be conserved among the Src family of kinases, thereby potentially representing similar PIP interactions within the family.

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.

The guidance and adhesion protein FLRT2 dimerizes in cis via dual small-X3-small transmembrane motifs

Fibronectin Leucine-rich Repeat Transmembrane (FLRT 1-3) proteins are a family of broadly expressed single-spanning transmembrane receptors that play key roles in development. Their extracellular domains mediate homotypic cell-cell adhesion and heterotypic protein interactions with other receptors to regulate cell adhesion and guidance. These in trans FLRT interactions determine the formation of signaling complexes of varying complexity and function. Whether FLRTs also interact at the surface of the same cell, in cis, remains unknown. Here, molecular dynamics simulations reveal two dimerization motifs in the FLRT2 transmembrane helix. Single particle tracking experiments show that these Small-X₃-Small motifs synergize with a third dimerization motif encoded in the extracellular domain to permit the cis association and co-diffusion patterns of FLRT2 receptors on cells. These results may point to a competitive switching mechanism between in cis and in trans interactions, which suggests that homotypic FLRT interaction mirrors the functionalities of classic adhesion molecules.

Systematic simulation of the interactions of pleckstrin homology domains with membranes

Pleckstrin homology (PH) domains can recruit proteins to membranes by recognition of phosphatidylinositol phosphate (PIP) lipids. Several family members are linked to diseases including cancer. We report the systematic simulation of the interactions of 100 mammalian PH domains with PIP-containing membranes. The observed PIP interaction hotspots recapitulate crystallographic binding sites and reveal a number of insights: (i) The β1 and β2 strands and their connecting loop constitute the primary PIP interaction site but are typically supplemented by interactions at the β3-β4 and β5-β6 loops; (ii) we reveal exceptional cases such as the Exoc8 PH domain; (iii) PH domains adopt different membrane-bound orientations and induce clustering of anionic lipids; and (iv) beyond family-level insights, our dataset sheds new light on individual PH domains, e.g., by providing molecular detail of secondary PIP binding sites. This work provides a global view of PH domain/membrane association involving multivalent association with anionic lipids.

Characterization of the membrane interactions of phospholipase Cγ reveals key features of the active enzyme

PLCγ enzymes are autoinhibited in resting cells and form key components of intracellular signaling that are also linked to disease development. Insights into physiological and aberrant activation of PLCγ require understanding of an active, membrane-bound form, which can hydrolyze inositol-lipid substrates. Here, we demonstrate that PLCγ1 cannot bind membranes unless the autoinhibition is disrupted. Through extensive molecular dynamics simulations and experimental evidence, we characterize membrane binding by the catalytic core domains and reveal previously unknown sites of lipid interaction. The identified sites act in synergy, overlap with autoinhibitory interfaces, and are shown to be critical for the phospholipase activity in cells. This work provides direct evidence that PLCγ1 is inhibited through obstruction of its membrane-binding surfaces by the regulatory region and that activation must shift PLCγ1 to a conformation competent for membrane binding. Knowledge of the critical sites of membrane interaction extends the mechanistic framework for activation, dysregulation, and therapeutic intervention.

Pocket delipidation induced by membrane tension or modification leads to a structurally analogous mechanosensitive channel state

The mechanosensitive ion channel of large conductance MscL gates in response to membrane tension changes. Lipid removal from transmembrane pockets leads to a concerted structural and functional MscL response, but it remains unknown whether there is a correlation between the tension-mediated state and the state derived by pocket delipidation in the absence of tension. Here, we combined pulsed electron paramagnetic resonance spectroscopy and hydrogen-deuterium exchange mass spectrometry, coupled with molecular dynamics simulations under membrane tension, to investigate the structural changes associated with the distinctively derived states. Whether it is tension- or modification-mediated pocket delipidation, we find that MscL samples a similar expanded subconducting state. This is the final step of the delipidation pathway, but only an intermediate stop on the tension-mediated path, with additional tension triggering further channel opening. Our findings hint at synergistic modes of regulation by lipid molecules in membrane tension-activated mechanosensitive channels.

Organization and Dynamics of the Red Blood Cell Band 3 Anion Exchanger SLC4A1: Insights From Molecular Dynamics Simulations

Molecular dynamics (MD) simulations have provided new insights into the organization and dynamics of the red blood cell Band 3 anion exchanger (AE1, SLC4A1). Band 3, like many solute carriers, works by an alternating access mode of transport where the protein rapidly (104/s) changes its conformation between outward and inward-facing states via a transient occluded anion-bound intermediate. While structural studies of membrane proteins usually reveal valuable structural information, these studies provide a static view often in the presence of detergents. Membrane transporters are embedded in a lipid bilayer and associated lipids play a role in their folding and function. In this review, we highlight MD simulations of Band 3 in realistic lipid bilayers that revealed specific lipid and protein interactions and were used to re-create a model of the Wright (Wr) blood group antigen complex of Band 3 and Glycophorin A. Current MD studies of Band 3 and related transporters are focused on describing the trajectory of substrate binding and translocation in real time. A structure of the intact Band 3 protein has yet to be achieved experimentally, but cryo-electron microscopy in combination with MD simulations holds promise to capture the conformational changes associated with anion transport in exquisite molecular detail.

Sphingomyelinase Disables Inactivation in Endogenous PIEZO1 Channels

Endogenous PIEZO1 channels of native endothelium lack the hallmark inactivation often seen when these channels are overexpressed in cell lines. Because prior work showed that the force of shear stress activates sphingomyelinase in endothelium, we considered if sphingomyelinase is relevant to endogenous PIEZO1. Patch clamping was used to quantify PIEZO1-mediated signals in freshly isolated murine endothelium exposed to the mechanical forces caused by shear stress and membrane stretch. Neutral sphingomyelinase inhibitors and genetic disruption of sphingomyelin phosphodiesterase 3 (SMPD3) cause PIEZO1 to switch to profoundly inactivating behavior. Ceramide (a key product of SMPD3) rescues non-inactivating channel behavior. Its co-product, phosphoryl choline, has no effect. In contrast to ceramide, sphingomyelin (the SMPD3 substrate) does not affect inactivation but alters channel force sensitivity. The data suggest that sphingomyelinase activity, ceramide, and sphingomyelin are determinants of native PIEZO gating that enable sustained activity.

•

SMPD3 sphingomyelinase enables long-lasting PIEZO1 activity in response to force

•

Ceramide, a key lipid product of SMPD3, promotes long-lasting activity

•

Sphingomyelin, the SMPD3 substrate, does not affect the duration of activity

•

Sphingomyelin alters PIEZO1 force sensitivity

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.

Allosteric activation of T cell antigen receptor signaling by quaternary structure relaxation

The mechanism of T cell antigen receptor (TCR-CD3) signaling remains elusive. Here, we identify mutations in the transmembrane region of TCRβ or CD3ζ that augment peptide T cell antigen receptor (pMHC)-induced signaling not explicable by enhanced ligand binding, lateral diffusion, clustering, or co-receptor function. Using a biochemical assay and molecular dynamics simulation, we demonstrate that the gain-of-function mutations loosen the interaction between TCRαβ and CD3ζ. Similar to the activating mutations, pMHC binding reduces TCRαβ cohesion with CD3ζ. This event occurs prior to CD3ζ phosphorylation and at 0°C. Moreover, we demonstrate that soluble monovalent pMHC alone induces signaling and reduces TCRαβ cohesion with CD3ζ in membrane-bound or solubilised TCR-CD3. Our data provide compelling evidence that pMHC binding suffices to activate allosteric changes propagating from TCRαβ to the CD3 subunits, reconfiguring interchain transmembrane region interactions. These dynamic modifications could change the arrangement of TCR-CD3 boundary lipids to license CD3ζ phosphorylation and initiate signal propagation.

•

Mutations in TCRβ and CD3ζ TMRs that reduce their interaction augment signaling

•

pMHC and anti-CD3 binding to TCR-CD3 induce similar quaternary structure relaxation

•

Soluble monovalent pMHC alone signals and reduces TCRαβ cohesion with CD3ζ

•

Allosteric changes in TCR-CD3 dynamics instigate 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 (PIP₂) molecules, and 3) cholesterol changes the depth of the dome and PIP₂ 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 PIP₂ interactions.

Advances in understanding and in multi-disciplinary methodology used to assess lipid regulation of signalling cascades from the cancer cell plasma membrane

The lipid bilayer is a functional component of cells, forming a stable platform for the initiation of key biological processes, including cell signalling. There are distinct changes in the lipid composition of cell membranes during oncogenic transformation resulting in aberrant activation and inactivation of signalling transduction pathways. Studying the role of the cell membrane in cell signalling is challenging, since techniques are often limited to by timescale, resolution, sensitivity, and averaging. To overcome these limitations, combining 'computational', 'wet-lab' and 'semi-dry' approaches offers the best opportunity to resolving complex biological processes involved in membrane organisation. In this review, we highlight analytical tools that have been applied for the study of cell signalling initiation from the cancer cell membranes through computational microscopy, biological assays, and membrane biophysics. The cancer therapeutic potential of extracellular membrane-modulating agents, such as cholesterol-reducing agents is also discussed, as is the need for future collaborative inter-disciplinary research for studying the role of the cell membrane and its components in cancer therapy.

Distortion of the bilayer and dynamics of the BAM complex in lipid nanodiscs

The β-barrel assembly machinery (BAM) catalyses the folding and insertion of β-barrel outer membrane proteins (OMPs) into the outer membranes of Gram-negative bacteria by mechanisms that remain unclear. Here, we present an ensemble of cryoEM structures of the E. coli BamABCDE (BAM) complex in lipid nanodiscs, determined using multi-body refinement techniques. These structures, supported by single-molecule FRET measurements, describe a range of motions in the BAM complex, mostly localised within the periplasmic region of the major subunit BamA. The β-barrel domain of BamA is in a 'lateral open' conformation in all of the determined structures, suggesting that this is the most energetically favourable species in this bilayer. Strikingly, the BAM-containing lipid nanodisc is deformed, especially around BAM's lateral gate. This distortion is also captured in molecular dynamics simulations, and provides direct structural evidence for the lipid 'disruptase' activity of BAM, suggested to be an important part of its functional mechanism.

Inter-domain dynamics in the chaperone SurA and multi-site binding to its outer membrane protein clients

The periplasmic chaperone SurA plays a key role in outer membrane protein (OMP) biogenesis. E. coli SurA comprises a core domain and two peptidylprolyl isomerase domains (P1 and P2), but its mechanisms of client binding and chaperone function have remained unclear. Here, we use chemical cross-linking, hydrogen-deuterium exchange mass spectrometry, single-molecule FRET and molecular dynamics simulations to map the client binding site(s) on SurA and interrogate the role of conformational dynamics in OMP recognition. We demonstrate that SurA samples an array of conformations in solution in which P2 primarily lies closer to the core/P1 domains than suggested in the SurA crystal structure. OMP binding sites are located primarily in the core domain, and OMP binding results in conformational changes between the core/P1 domains. Together, the results suggest that unfolded OMP substrates bind in a cradle formed between the SurA domains, with structural flexibility between domains assisting OMP recognition, binding and release.

Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes

Association of peripheral proteins with lipid bilayers regulates membrane signaling and dynamics. Pleckstrin homology (PH) domains bind to phosphatidylinositol phosphate (PIP) molecules in membranes. The effects of local PIP enrichment on the interaction of PH domains with membranes is unclear. Molecular dynamics simulations allow estimation of the binding energy of GRP1 PH domain to PIP3-containing membranes. The free energy of interaction of the PH domain with more than two PIP3 molecules is comparable to experimental values, suggesting that PH domain binding involves local clustering of PIP molecules within membranes. We describe a mechanism of PH binding proceeding via an encounter state to two bound states which differ in the orientation of the protein relative to the membrane, these orientations depending on the local PIP concentration. These results suggest that nanoscale clustering of PIP molecules can control the strength and orientation of PH domain interaction in a concentration-dependent manner.

The Function of Membrane Integral Pyrophosphatases From Whole Organism to Single Molecule

Membrane integral pyrophosphatases (mPPases) are responsible for the hydrolysis of pyrophosphate. This enzymatic mechanism is coupled to the pumping of H+ or Na+ across membranes in a process that can be K+ dependent or independent. Understanding the movements and dynamics throughout the mPPase catalytic cycle is important, as this knowledge is essential for improving or impeding protein function. mPPases have been shown to play a crucial role in plant maturation and abiotic stress tolerance, and so have the potential to be engineered to improve plant survival, with implications for global food security. mPPases are also selectively toxic drug targets, which could be pharmacologically modulated to reduce the virulence of common human pathogens. The last few years have seen the publication of many new insights into the function and structure of mPPases. In particular, there is a new body of evidence that the catalytic cycle is more complex than originally proposed. There are structural and functional data supporting a mechanism involving half-of-the-sites reactivity, inter-subunit communication, and exit channel motions. A more advanced and in-depth understanding of mPPases has begun to be uncovered, leaving the field of research with multiple interesting avenues for further exploration and investigation.

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.

Membrane Recognition and Binding by the Phosphatidylinositol Phosphate Kinase PIP5K1A: A Multiscale Simulation Study

Phosphatidylinositol phosphates (PIPs) are lipid signaling molecules that play key roles in many cellular processes. PIP5K1A kinase catalyzes phosphorylation of PI4P to form PIP₂, which in turn interacts with membrane and membrane-associated proteins. We explore the mechanism of membrane binding by the PIP5K1A kinase using a multiscale molecular dynamics approach. Coarse-grained simulations show binding of monomeric PIP5K1A to a model cell membrane containing PI4P. PIP5K1A did not bind to zwitterionic or anionic membranes lacking PIP molecules. Initial encounter of kinase and bilayer was followed by reorientation to enable productive binding to the PI4P-containing membrane. The simulations suggest that unstructured regions may be important for the preferred orientation for membrane binding. Atomistic simulations indicated that the dimeric kinase could not bind to the membrane via both active sites at the same time, suggesting a conformational change in the protein and/or bilayer distortion may be needed for dual-site binding to occur.

•

PIP5K1A kinase interacts with PIP-containing membranes via its activation loop

•

PIP5K1A does not bind to zwitterionic or anionic membranes lacking PIP molecules

•

Initial encounter of protein and bilayer is followed by reorientation and binding

•

Dimeric PIP5K1A binds with membrane contacts via only one catalytic site at a time

β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.

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 (PIP₂) 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 PIP₂ 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.

Human Band 3 (AE1, SLC4A1), an abundant 911 amino acid glycoprotein, catalyzes the exchange of bicarbonate and chloride across the red blood cell membrane, a process necessary for efficient respiration. Malfunction of Band 3 leads to inherited diseases such as Southeast Asian Ovalocytosis, hereditary spherocytosis and distal renal tubular acidosis. Despite much available structural and functional data about Band 3, key questions about the conformational changes associated with transport and the molecular details of its interaction with lipids and other proteins remain unanswered. In this study, we have used computer simulations to investigate the dynamics of Band 3 in lipid bilayers that resemble the red blood cell plasma membrane. Our results suggest that negatively charged phospholipids and cholesterol interact strongly with Band 3 forming an annulus around the protein. Glycophorin A (GPA) interacts with Band 3 to form the Wright (Wr) blood group antigen. We were able to recreate this complex and show that GPA promotes the clustering of Band 3 in red blood cell membranes. Understanding the molecular details of the interaction of Band 3 with GPA has provided new insights into the nature of the Wright blood group antigen.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

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 Ca²⁺ nor any other metal ion was bound to WT or variant KIBRA C2 in crystal structures, and Ca²⁺ 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.

Effects of Periplasmic Chaperones and Membrane Thickness on BamA-Catalyzed Outer-Membrane Protein Folding

The biogenesis of outer-membrane proteins (OMPs) in gram-negative bacteria involves delivery by periplasmic chaperones to the β-barrel assembly machinery (BAM), which catalyzes OMP insertion into the outer membrane. Here, we examine the effects of membrane thickness, the Escherichia coli periplasmic chaperones Skp and SurA, and BamA, the central subunit of the BAM complex, on the folding kinetics of a model OMP (tOmpA) using fluorescence spectroscopy, native mass spectrometry, and molecular dynamics simulations. We show that prefolded BamA promotes the release of tOmpA from Skp despite the nM affinity of the Skp:tOmpA complex. This activity is located in the BamA β-barrel domain, but is greater when full-length BamA is present, indicating that both the β-barrel and polypeptide transport-associated (POTRA) domains are required for maximal activity. By contrast, SurA is unable to release tOmpA from Skp, providing direct evidence against a sequential chaperone model. By varying lipid acyl chain length in synthetic liposomes we show that BamA has a greater catalytic effect on tOmpA folding in thicker bilayers, suggesting that BAM catalysis involves lowering of the kinetic barrier imposed by the hydrophobic thickness of the membrane. Consistent with this, molecular dynamics simulations reveal that increases in membrane thinning/disorder by the transmembrane domain of BamA is greatest in thicker bilayers. Finally, we demonstrate that cross-linking of the BamA barrel does not affect tOmpA folding kinetics in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) liposomes, suggesting that lateral gating of the BamA barrel and/or hybrid barrel formation is not required, at least for the assembly of a small 8-stranded OMP in vitro.

•

Mechanisms of OMP periplasmic transport and folding by BAM are poorly understood.

•

BamA catalyzes folding by reducing the kinetic barrier imposed by membrane thickness.

•

BamA proteoliposomes promote folding of Skp-bound tOmpA.

•

Lateral gating is not required for BamA-catalyzed folding of tOmpA in DMPC bilayers.

Dynamic interactions between a membrane binding protein and lipids induce fluctuating diffusivity

Pleckstrin homology (PH) domains are membrane-binding lipid recognition proteins that interact with phosphatidylinositol phosphate (PIP) molecules in eukaryotic cell membranes. Diffusion of PH domains plays a critical role in biological reactions on membrane surfaces. Although diffusivity can be estimated by long-time measurements, it lacks information on the short-time diffusive nature. We reveal two diffusive properties of a PH domain bound to the surface of a PIP-containing membrane using molecular dynamics simulations. One is fractional Brownian motion, attributed to the motion of the lipids with which the PH domain interacts. The other is temporally fluctuating diffusivity; that is, the short-time diffusivity of the bound protein changes substantially with time. Moreover, the diffusivity for short-time measurements is intrinsically different from that for long-time measurements. This fluctuating diffusivity results from dynamic changes in interactions between the PH domain and PIP molecules. Our results provide evidence that the complexity of protein-lipid interactions plays a crucial role in the diffusion of proteins on biological membrane surfaces. Changes in the diffusivity of PH domains and related membrane-bound proteins may in turn contribute to the formation/dissolution of protein complexes in membranes.

On the interpretation of reflectivity data from lipid bilayers in terms of molecular-dynamics models

Neutron and X-ray reflectivity of model membranes is increasingly used as a tool for the study of membrane structures and dynamics. As the systems under study become more complex, and as long, all-atom molecular-dynamics (MD) simulations of membranes become more available, there is increasing interest in the use of MD simulations in the analysis of reflectometry data from membranes. In order to perform this, it is necessary to produce a model of the complete interface, including not only the MD-derived structure of the membrane, but also the supporting substrate and any other interfacial layers that may be present. Here, it is shown that this is best performed by first producing a model of the occupied volume across the entire interface, and then converting this into a scattering length density (SLD) profile, rather than by splicing together the separate SLD profiles from the substrate layers and the membrane, since the latter approach can lead to discontinuities in the SLD profile and subsequent artefacts in the reflectivity calculation. It is also shown how the MD-derived membrane structure should be corrected to account for lower than optimal coverage and out-of-plane membrane fluctuations. Finally, the method of including the entire membrane structure in the reflectivity calculation is compared with an alternative approach in which the membrane components are approximated by functional forms, with only the component volumes being extracted from the simulation. It is shown that using only the fragment volumes is insufficient for a typical neutron data set of a single deuteration measured at several water contrasts, and that either weighting the model by including more structural information from the fit, or a larger data set involving a range of deuterations, are required to satisfactorily define the problem.

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.

Neuromuscular junction formation and maintenance is an essential biological process, the disruption of which leads to congenital myasthenic syndromes and premature death. Dok7 is a key member in formation, maintenance and signaling in neuromuscular junctions. Dok7 is a peripheral membrane protein that is necessary for full activation of the receptor tyrosine kinase MuSK, a receptor tyrosine kinase residing in the postsynaptic membrane. The structure of Dok7 consists of both a PH domain and a PTB domain. The interaction of Dok7 with cell membranes is not well understood. Here, we use molecular simulations to show that the Dok7 PH domain preferentially binds to PIP lipid molecules when associating with a membrane. Dok7 interacts with the bilayer in both a canonical binding mode and an alternative binding mode. Our simulations also demonstrate the presence of both a canonical and an atypical binding site for PIPs in agreement with recent crystallographic studies of the ASAP1 PH domain. Analysis of density landscapes for the interaction of the Dok7 PH domain with PIP-containing lipid bilayers is also able to identify both canonical and alternative binding modes. Our improved understanding of how Dok7 interacts with a membrane is key to examining Dok7/MuSK signaling.

Interactions of Pleckstrin Homology Domains with Membranes: Adding Back the Bilayer via High-Throughput Molecular Dynamics

A molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

Association of Peripheral Membrane Proteins with Membranes: Free Energy of Binding of GRP1 PH Domain with Phosphatidylinositol Phosphate-Containing Model Bilayers

Understanding the energetics of peripheral protein-membrane interactions is important to many areas of biophysical chemistry and cell biology. Estimating free-energy landscapes by molecular dynamics (MD) simulation is challenging for such systems, especially when membrane recognition involves complex lipids, e.g., phosphatidylinositol phosphates (PIPs). We combined coarse-grained MD simulations with umbrella sampling to quantify the binding of the well-explored GRP1 pleckstrin homology (PH) domain to model membranes containing PIP molecules. The experimentally observed preference of GRP1-PH for PIP3 over PIP2 was reproduced. Mutation of a key residue (K273A) within the canonical PIP-binding site significantly reduced the free energy of PIP binding. The presence of a noncanonical PIP-interaction site, observed experimentally in other PH domains but not previously in GRP1-PH, was also revealed. These studies demonstrate how combining coarse-grained simulations and umbrella sampling can unmask the molecular basis of the energetics of interactions between peripheral membrane proteins and complex cellular membranes.

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.

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.

Symporters are proteins that are responsible for the co-transport of ions and small molecule solutes across cell membranes. UraA is an example of a symporter, and is responsible for the proton-driven uptake of uracil in bacteria like E. coli. Despite its importance as a member of a large family of nucleobase/ascorbate transporters (NAT) and the existence of structural and functional data, the mechanism by which UraA transports uracil across the bacterial membrane, and in particular the role of its diverse and complex lipid environment in the transport mechanism, remains elusive. In this study, we have used a multiscale computational methodology to examine the dynamics of UraA and to elucidate its interactions with lipids that resemble its native environment in the bacterial inner membrane. Our results demonstrate that negatively-charged lipids in the membrane (phosphatidylglycerol and cardiolipin) associate preferentially with UraA and may play a role in its function. Additionally, our simulations resulted in a closed state of UraA, a likely intermediate in the transport mechanism that may not be readily accessible by experimental methods.

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 PIP₃ 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 PIP₃ 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, PIP₃ 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 the auxilin-1 PTEN-like domain with model membranes result in nanoclustering of phosphatidyl inositol phosphates

Auxilin-1 is a neuron-specific membrane-binding protein involved in a late stage of clathrin-mediated endocytosis. It recruits Hsc70, thus initiating uncoating of the clathrin-coated vesicles. Interactions of auxilin-1 with the vesicle membrane are crucial for this function and are mediated via an N-terminal PTEN-like domain. We have used multiscale molecular dynamics simulations to probe the interactions of the auxilin-1 PTEN-like domain with lipid bilayers containing differing phospholipid composition, including bilayers containing phosphatidyl inositol phosphates. Our results suggest a novel, to our knowledge, model for the auxilin/membrane encounter and subsequent interactions. Negatively charged lipids (especially PIP2) enhance binding of auxilin to lipid bilayers and facilitate its correct orientation relative to the membrane. Mutations in three basic residues (R301E/R307E/K311E) of the C2 subdomain of the PTEN-like domain perturbed its interaction with the bilayer, changing its orientation. The interaction of membrane-bound auxilin-1 PTEN-like domain with negatively charged lipid headgroups results in nanoclustering of PIP2 molecules in the adjacent bilayer leaflet.

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.

Transmission of signals across the cell membrane is an essential process for all living organisms. Integrins are one example of cell surface receptors (αβ) which, uniquely, form a bidirectional signalling pathway across the membrane. Integrins are crucial for many cellular processes and play key roles in pathological defects such as cardiovascular diseases and cancer. They are activated to a high affinity state by the intracellular protein talin in a process known as ‘inside-out activation’. Despite their importance and the existence of functional and structural data, the mechanism by which talin activates integrin remains elusive. In this study we use a multi-scale computational approach, which combines coarse-grained and atomistic molecular dynamics simulations, to suggest how the formation of the complex between the talin head domain, the cell membrane and the integrin moves the integrin equilibrium towards an active state. Our results show that conformational changes within the talin head domains optimize its interactions with the cell membrane. Upon binding to the integrin, talin facilitates rearrangement of the integrin TM region thus promoting integrin activation. This study also provides a demonstration of the strengths of a computational multi-scale approach in studies of membrane interactions and receptor conformational changes and associated proteins that enable transmembrane signaling.

A helix heterodimer in a lipid bilayer: prediction of the structure of an integrin transmembrane domain via multiscale simulations

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling.

The structure of the talin/integrin complex at a lipid bilayer: an NMR and MD simulation study

Integrins are cell surface receptors crucial for cell migration and adhesion. They are activated by interactions of the talin head domain with the membrane surface and the integrin β cytoplasmic tail. Here, we use coarse-grained molecular dynamic simulations and nuclear magnetic resonance spectroscopy to elucidate the membrane-binding surfaces of the talin head (F2-F3) domain. In particular, we show that mutations in the four basic residues (K258E, K274E, R276E, and K280E) in the F2 binding surface reduce the affinity of the F2-F3 for the membrane and modify its orientation relative to the bilayer. Our results highlight the key role of anionic lipids in talin/membrane interactions. Simulation of the F2-F3 in complex with the α/β transmembrane dimer reveals information for its orientation relative to the membrane. Our studies suggest that the perturbed orientation of talin relative to the membrane in the F2 mutant would be expected to in turn perturb talin/integrin interactions.