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

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

Distinctive phosphoinositide- and Ca²⁺-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.

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

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

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

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

Modelling 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

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

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