Coupling the role of lipids to the conformational dynamics of the ABC transporter P-glycoprotein.

The ATP-binding cassette transporter P-glycoprotein (P-gp) is a multidrug efflux pump that is overexpressed in a variety of cancers and associated with the drug-resistance phenomenon. P-gp structures were previously determined in detergent and in nanodiscs, in which different transmembrane helix conformations were found, "straight" and "kinked," respectively, indicating a possible role of the lipid environment on the P-gp structural ensemble. Here, we investigate the dynamic conformational ensembles and protein-lipid interactions of two human P-gp inward-open conformers, straight and kinked, employing all-atom molecular dynamics (MD) simulations in asymmetric multicomponent lipid bilayers that mimic the highly specialized hepatocyte membrane in which P-gp is expressed. The two conformers are found to differ in terms of the accessibility of the substrate cavity. The MD simulations show how cholesterol and different lipid species wedge, snorkel, and partially enter into the cavity of the straight P-gp conformer solved in detergent. However, access to the cavity of the kinked P-gp conformer solved in nanodiscs is restricted. Furthermore, the volume and dynamic fluctuations of the substrate cavity largely differ between the two P-gp conformers and are modulated by the presence (or absence) of cholesterol in the membrane and/or of ATP. From the mechanistic perspective, the findings indicate that the straight conformer likely precedes the kinked conformer in the functional working cycle of P-gp, with the latter conformation representing a post substrate-bound state. The inaccessibility of the main transmembrane cavity in the kinked conformer might be crucial in preventing substrate disengagement and transport withdrawal. Remarkably, in our unbiased MD simulations, one transmembrane helix (TM10) of the straight conformer underwent a spontaneous transition to a kinked conformation, underlining the relevance of both conformations in a native phospholipid environment and revealing structural descriptors defining the transition between the two P-gp conformers.

Anionic Phospholipids Stimulate the Proton Pumping Activity of the Plant Plasma Membrane P-Type H+-ATPase.

The activity of membrane proteins depends strongly on the surrounding lipid environment. Here, we characterize the lipid stimulation of the plant plasma membrane H+-ATPase Arabidopsis thaliana H+-ATPase isoform 2 (AHA2) upon purification and reconstitution into liposomes of defined lipid compositions. We show that the proton pumping activity of AHA2 is stimulated by anionic phospholipids, especially by phosphatidylserine. This activation was independent of the cytoplasmic C-terminal regulatory domain of the pump. Molecular dynamics simulations revealed several preferential contact sites for anionic phospholipids in the transmembrane domain of AHA2. These contact sites are partially conserved in functionally different P-type ATPases from different organisms, suggesting a general regulation mechanism by the membrane lipid environment. Our findings highlight the fact that anionic lipids play an important role in the control of H+-ATPase activity.

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.

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.

A Molecular Perspective on Mitochondrial Membrane Fusion: From the Key Players to Oligomerization and Tethering of Mitofusin.

Mitochondria are dynamic organelles characterized by an ultrastructural organization which is essential in maintaining their quality control and ensuring functional efficiency. The complex mitochondrial network is the result of the two ongoing forces of fusion and fission of inner and outer membranes. Understanding the functional details of mitochondrial dynamics is physiologically relevant as perturbations of this delicate equilibrium have critical consequences and involved in several neurological disorders. Molecular actors involved in this process are large GTPases from the dynamin-related protein family. They catalyze nucleotide-dependent membrane remodeling and are widely conserved from bacteria to higher eukaryotes. Although structural characterization of different family members has contributed in understanding molecular mechanisms of mitochondrial dynamics in more detail, the complete structure of some members as well as the precise assembly of functional oligomers remains largely unknown. As increasing structural data become available, the domain modularity across the dynamin superfamily emerged as a foundation for transfering the knowledge towards less characterized members. In this review, we will first provide an overview of the main actors involved in mitochondrial dynamics. We then discuss recent example of computational methodologies for the study of mitofusin oligomers, and present how the usage of integrative modeling in conjunction with biochemical data can be an asset in progressing the still challenging field of membrane dynamics.

Structural basis for triacylglyceride extraction from mycobacterial inner membrane by MFS transporter Rv1410.

Mycobacterium tuberculosis is protected from antibiotic therapy by a multi-layered hydrophobic cell envelope. Major facilitator superfamily (MFS) transporter Rv1410 and the periplasmic lipoprotein LprG are involved in transport of triacylglycerides (TAGs) that seal the mycomembrane. Here, we report a 2.7 Å structure of a mycobacterial Rv1410 homologue, which adopts an outward-facing conformation and exhibits unusual transmembrane helix 11 and 12 extensions that protrude ~20 Å into the periplasm. A small, very hydrophobic cavity suitable for lipid transport is constricted by a functionally important ion-lock likely involved in proton coupling. Combining mutational analyses and MD simulations, we propose that TAGs are extracted from the core of the inner membrane into the central cavity via lateral clefts present in the inward-facing conformation. The functional role of the periplasmic helix extensions is to channel the extracted TAG into the lipid binding pocket of LprG.

The ABC transporter MsbA adopts the wide inward-open conformation in E. coli cells.

Membrane proteins are currently investigated after detergent extraction from native cellular membranes and reconstitution into artificial liposomes or nanodiscs, thereby removing them from their physiological environment. However, to truly understand the biophysical properties of membrane proteins in a physiological environment, they must be investigated within living cells. Here, we used a spin-labeled nanobody to interrogate the conformational cycle of the ABC transporter MsbA by double electron-electron resonance. Unexpectedly, the wide inward-open conformation of MsbA, commonly considered a nonphysiological state, was found to be prominently populated in Escherichia coli cells. Molecular dynamics simulations revealed that extensive lateral portal opening is essential to provide access of its large natural substrate core lipid A to the binding cavity. Our work paves the way to investigate the conformational landscape of membrane proteins in cells.

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.

Structural dataset from microsecond-long simulations of yeast mitofusin Fzo1 in the context of membrane docking.

In this work we present a novel set of possible auto-oligomerisation states of yeast protein Fzo1 in the context of membrane docking. The dataset reports atomistic models and trajectories derived from a molecular dynamics study of the yeast mitofusin Fzo1, residues 101-855. The initial modelling was followed by coarse-grained molecular dynamics simulation to evaluate the stability and the dynamics of each structural model in a solvated membrane environment. Simulations were run for 1 µs and collected with GROMACS v5.0.4 using the martini v2.1 force field. For each structural model, the dataset comprises the production phase under semi-isotropic condition at 1 bar, 310 K and 150 mn NaCl. The integration step is 20 fs and coordinates have been saved every 1 ns. Each trajectory is associated with a ready-available visualization state for the VMD software. These structural detailed informations are a ready-available platform to plan integrative studies on the mitofusin Fzo1 and will aid the community to further elucidate the mitochondrial tethering process during membrane fusion. This dataset is based on the publication "Physics-based oligomeric models of the yeast mitofusin Fzo1 at the molecular scale in the context of membrane docking." (Brandner and De Vecchis et al., 2019)".

Physics-based oligomeric models of the yeast mitofusin Fzo1 at the molecular scale in the context of membrane docking.

Tethering and homotypic fusion of mitochondrial outer membranes is mediated by large GTPases of the dynamin-related proteins family called the mitofusins. The yeast mitofusin Fzo1 forms high molecular weight complexes and its assembly during membrane fusion likely involves the formation of high order complexes. Consistent with this possibility, mitofusins form oligomers in both cis (on the same lipid bilayer) and trans to mediate membrane attachment and fusion. Here, we utilize our recent Fzo1 model to investigate and discuss the formation of cis and trans mitofusin oligomers. We have built three distinct cis-assembly Fzo1 models that gave rise to three distinct trans-oligomeric models of mitofusin constructs. Each model involves two main components of mitofusin oligomerization: the GTPase and the trunk domains. The oligomeric models proposed in this study were further assessed for stability and dynamics in a membrane environment using a coarse-grained molecular dynamics (MD) simulation approach. A narrow opening 'head-to-head' cis-oligomerization (via the GTPase domain) followed by the antiparallel 'back-to-back' trans-associations (via the trunk domain) appears to be in agreement with all of the available experimental data. More broadly, this study opens new possibilities to start exploring cis and trans conformations for Fzo1 and mitofusins in general.

A membrane-inserted structural model of the yeast mitofusin Fzo1.

Mitofusins are large transmembrane GTPases of the dynamin-related protein family, and are required for the tethering and fusion of mitochondrial outer membranes. Their full-length structures remain unknown, which is a limiting factor in the study of outer membrane fusion. We investigated the structure and dynamics of the yeast mitofusin Fzo1 through a hybrid computational and experimental approach, combining molecular modelling and all-atom molecular dynamics simulations in a lipid bilayer with site-directed mutagenesis and in vivo functional assays. The predicted architecture of Fzo1 improves upon the current domain annotation, with a precise description of the helical spans linked by flexible hinges, which are likely of functional significance. In vivo site-directed mutagenesis validates salient aspects of this model, notably, the long-distance contacts and residues participating in hinges. GDP is predicted to interact with Fzo1 through the G1 and G4 motifs of the GTPase domain. The model reveals structural determinants critical for protein function, including regions that may be involved in GTPase domain-dependent rearrangements.