Search

Hyaluronan is an acidic heteropolysaccharide comprising alternating N-acetylglucosamine and glucuronic acid sugars that is ubiquitously expressed in the vertebrate extracellular matrix1. The high-molecular-mass polymer modulates essential physiological processes in health and disease, including cell differentiation, tissue homeostasis and angiogenesis2. Hyaluronan is synthesized by a membrane-embedded processive glycosyltransferase, hyaluronan synthase (HAS), which catalyses the synthesis and membrane translocation of hyaluronan from uridine diphosphate-activated precursors3,4. Here we describe five cryo-electron microscopy structures of a viral HAS homologue at different states during substrate binding and initiation of polymer synthesis. Combined with biochemical analyses and molecular dynamics simulations, our data reveal how HAS selects its substrates, hydrolyses the first substrate to prime the synthesis reaction, opens a hyaluronan-conducting transmembrane channel, ensures alternating substrate polymerization and coordinates hyaluronan inside its transmembrane pore. Our research suggests a detailed model for the formation of an acidic extracellular heteropolysaccharide and provides insights into the biosynthesis of one of the most abundant and essential glycosaminoglycans in the human body.

Patched1 (PTCH1) is the principal tumour suppressor protein of the mammalian Hedgehog (HH) signalling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the Class F G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating accessible cholesterol levels within ciliary membranes. Using extensive molecular dynamics (MD) simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15-20 kJ mol-1for cholesterol export. In simulations we identify cation binding sites within the PTCH1 transmembrane domain (TMD) which may provide the energetic impetus for cholesterol transport.In silicodata are coupled toin vivobiochemical assays of PTCH1 mutants to probe coupling between transmembrane motions and PTCH1 activity. Using complementary simulations of Dispatched1 (DISP1) we find that transition between ‘inward-open’ and solvent ‘occluded’ states is accompanied by Na+induced pinching of intracellular helical segments. Thus, our findings illuminate the energetics and ion-coupling stoichiometries of PTCH1 transport mechanisms, whereby 1-3 Na+or 2-3 K+couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Mycobacterium tuberculosis ( Mtb ) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

SummaryInoviruses are abundant filamentous phages infecting numerous prokaryotic phyla, where they can symbiotically promote host fitness and increase bacterial virulence. Due to their unique properties, inoviruses have also been utilised in biotechnology for phage display and as models for studying phase behaviour of colloidal rods. Inoviral phages secreted by bacteria can self-assemble into liquid crystalline droplets that protect bacterial cells in biofilms from antibiotics, however, factors governing the formation of such droplets and the mechanism of antibiotic protection are poorly understood. Here, we investigate the structural, biophysical, and protective properties of liquid crystalline droplets formed byPseudomonas aeruginosaandEscherichia coliinoviral phages. We report a cryo-EM structure of the capsid from the highly studiedE. colifd phage, revealing distinct biochemical properties of fd compared to Pf4 phage fromP. aeruginosa. We show that fd and Pf4 form liquid crystalline droplets with diverse morphologies governed by the underlying phage particle geometry and biophysics, rather than their surface biochemical properties. Finally, we show that these morphologically diverse droplets made of either phage can protect rod-shaped bacteria from antibiotic treatment, despite differing modes of association with cells. This study advances our understanding of phage assembly into liquid crystalline droplets, and provides insights into how filamentous molecules protect bacteria from extraneous molecules under crowding conditions, which are found in biofilms or on infected host tissues.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis, a disease that claims ~1.5 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. There is much to be understood about the Mtb cell envelope, a complicated barrier that antibiotics need to negotiate to enter the cell. Within this envelope, the plasma membrane is the ultimate obstacle and is proposed to be comprised of over 50% mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs), whose role in the membrane structure remains elusive. Here we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to drive biophysical properties of the Mycobacterial plasma membrane. To validate the model, we tested known anti-tubercular drugs and replicated previous experimental results. Our results shed new light into the organization of the Mycobacterial plasma membrane and provides a working model of this complex membrane to use for in silico studies. This opens the door for new methods to probe potential antibiotic targets and further understand membrane protein function.Abstract Figure

KdpFABC is a high-affinity prokaryotic K+ uptake system that forms a functional chimera between a channel-like subunit (KdpA) and a P-type ATPase (KdpB). At high K+ levels, KdpFABC needs to be inhibited to prevent excessive K+ accumulation to the point of toxicity. This is achieved by a phosphorylation of the serine residue in the TGES162 motif in the A domain of the pump subunit KdpB (KdpBS162-P). Here, we explore the structural basis of inhibition by KdpBS162 phosphorylation by determining the conformational landscape of KdpFABC under inhibiting and non-inhibiting conditions. Under turnover conditions, we identified a new inhibited KdpFABC state that we termed E1P tight, which is not part of the canonical Post-Albers transport cycle of P-type ATPases. It likely represents the biochemically described stalled E1P state adopted by KdpFABC upon KdpBS162 phosphorylation. The E1P tight state exhibits a compact fold of the three cytoplasmic domains and is likely adopted when the transition from high-energy E1P states to E2P states is unsuccessful. This study represents a structural characterization of a biologically relevant off-cycle state in the P-type ATPase family and supports the emerging discussion of P-type ATPase regulation by such states.

In this study we used coarse-grained molecular dynamics (CG MD) simulations to study protein-cholesterol interactions for different activation states of the A2A adenosine receptor (A2AR) and A1 adenosine receptor (A1R). Three different membranes were used, of which a plasma mimetic membrane best matched cholesterol binding sites previously detected in X-ray structures and computationally for the inactive state of A2AR. Here, cholesterol binds stably between TM5 and TM6 in the intracellular leaf and between the extracellular part of TM6 and TM7. The stability of the identified binding sites to A1R with CG MD simulations were further investigated using potential of mean force calculations combined with umbrella sampling.Cholesterol binding to certain cavities in active or inactive state can stabilize a conformation of active or inactive state respectively and is important to identify these cavities. For the active state the cholesterol binding sites have much longer residence times compared to the inactive state for both A2AR and A1R.We observed for the active A1R two cholesterol binding sites in the extracellular membrane leaf in a cavity between TM2/TM3 and along TM3 in the intracellular leaf. We found a partially overlapped binding area between TM6 and TM7 in the extracellular membrane leaf for active state and inactive A1R state which can be antagonized by cholesterol molecules although the site has a much longer residence time for active A1R. For the active A2AR we predicted a high residence time cholesterol position between TM1 and TM7 in the middle region of TMD which is different from the binding sites calculated in active or inactive A1R and for the inactive A2AR. For the inactive A2AR cholesterol binds to the same extracellular area between TM6 and TM7 as in the inactive A1R.The differences in stable, long residence time cholesterol binding sites between active and inactive states of A1R and for A2AR can be important for functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators which can bind through lipid pathways.

Lipid droplets (LDs) are essential for cellular lipid homeostasis by storing diverse neutral lipids (NLs), such as triacylglycerol (TAG), steryl esters (SE), and retinyl esters (RE). A proper assembly of TAG-containing LDs at the ER requires Seipin, a conserved protein often mutated in lipodystrophies. Here, we show that the yeast Seipin Sei1 and its partner Ldb16 also promote the storage of other NL in LDs. Importantly, this role of Sei1/Ldb16 is evolutionarily conserved as expression of human-Seipin restored normal SE-containing LDs in yeast Seipin mutants. As in the case of TAG, the formation of SE-containing LDs requires interactions between hydroxyl-residues in human Seipin or yeast Ldb16 with NL carboxyl esters. These findings provide a universal mechanism for Seipin-mediated LD formation and suggest a model for how Seipin distinguishes NLs from aliphatic phospholipid acyl chains in the center of the membrane bilayer.

The mitochondrial electron transport chain comprises a series of protein complexes embedded in the inner mitochondrial membrane that generate a proton motive force via oxidative phosphorylation, ultimately generating ATP. These protein complexes can oligomerize to form larger structures called supercomplexes. Cardiolipin (CL), a conical lipid, unique within eukaryotes to the inner mitochondrial membrane, has proven essential in maintaining the stability and function of supercomplexes. Monolysocardiolipin (MLCL) is a CL variant that accumulates in people with Barth syndrome (BTHS). BTHS is caused by defects in CL biosynthesis and characterised by abnormal mitochondrial bioenergetics and destabilised supercomplexes. However, the mechanisms by which MLCL causes pathogenesis remain unclear. Here, multiscale molecular dynamics characterise the interactions of CL and MLCL with yeast and mammalian mitochondrial supercomplexes containing Complex III (CIII) and Complex IV (CIV). Coarse-grained simulations reveal that both CL and MLCL bind to sites at the interface between CIII and CIV of the supercomplex. Free energy perturbation calculations show that MLCL interaction is weaker than that of CL and suggest that interaction with CIV drives this difference. Atomistic contact analyses show that, although interaction with CIII is similar for CL and MLCL, CIV makes more contacts with CL than MLCL, demonstrating that CL is a more successful “glue” between the two complexes. Simulations of the human CIII2CIV supercomplex show that this interface site is maintained between species. Our study suggests that MLCL accumulation in people with BTHS disrupts supercomplex stability by formation of relatively weak interactions at the interface lipid binding site.

The endoplasmic reticulum (ER) membrane protein complex (EMC) was identified over a decade ago in a genetic screen for ER protein homeostasis. The EMC inserts transmembrane domains (TMDs) with limited hydrophobicity. Two recent cryo‐EM structures, and a third model based on partial high‐ and low‐resolution structures, suggest how this is accomplished.

Cryo-electron microscopy (cryo-EM) enables the determination of membrane protein structures in native-like environments. Characterising how membrane proteins interact with the surrounding membrane lipid environment is assisted by resolution of lipid-like densities visible in cryo-EM maps. Nevertheless, establishing the molecular identity of putative lipid and/or detergent densities remains challenging. Here we present LipIDens, a pipeline for molecular dynamics (MD) simulation-assisted interpretation of lipid and lipid-like densities in cryo-EM structures. The pipeline integrates the implementation and analysis of multi-scale MD simulations for identification, ranking and refinement of lipid binding poses which superpose onto cryo-EM map densities. Thus, LipIDens enables direct integration of experimental and computational structural approaches to facilitate the interpretation of lipid-like cryo-EM densities and to reveal the molecular identities of protein-lipid interactions within a bilayer environment. The LipIDens code is open-source and embedded within a notebook format to assist automation and usability.

KdpFABC is a high-affinity prokaryotic K+ uptake system that forms a functional chimera between a channel-like subunit (KdpA) and a P-type ATPase (KdpB). At high K+ levels, KdpFABC needs to be inhibited to prevent excessive K+ accumulation to the point of toxicity. This is achieved by a phosphorylation of the serine residue in the TGES162 motif in the A domain of the pump subunit KdpB (KdpBS162-P). Here, we explore the structural basis of inhibition by KdpBS162 phosphorylation by determining the conformational landscape of KdpFABC under inhibiting and non-inhibiting conditions. Under turnover conditions, we identified a new inhibited KdpFABC conformation that we termed E1-P tight, which is not part of the canonical Post-Albers transport cycle of P-type ATPases. It likely represents the biochemically described stalled E1-P state adopted by KdpFABC upon KdpBS162 phosphorylation. The E1-P tight state exhibits a compact fold of the three cytoplasmic domains and is likely adopted when the transition from high-energy E1-P states to E2-P states is unsuccessful. This study represents a structural characterization of a biologically relevant off-cycle state in the P-type ATPase family and supports the emerging discussion of P-type ATPase regulation by such conformations.

Transport of proteins across and into membranes is a fundamental biological process with the vast majority being conducted by the ubiquitous Sec machinery. In bacteria, this is usually achieved when the SecY-complex engages the cytosolic ATPase SecA (secretion) or translating ribosomes (insertion). Great strides have been made towards understanding the mechanism of protein translocation. Yet, important questions remain – notably, the nature of the individual steps that constitute transport, and how the proton-motive force (PMF) across the plasma membrane contributes. Here, we apply a recently developed high-resolution protein transport assay to explore these questions. We find that pre-protein transport is limited primarily by the diffusion of arginine residues across the membrane, particularly in the context of bulky hydrophobic sequences. This specific effect of arginine, caused by its positive charge, is mitigated for lysine which can be deprotonated and transported across the membrane in its neutral form. These observations have interesting implications for the mechanism of protein secretion, suggesting a simple mechanism through which the PMF can aid transport by enabling a 'proton ratchet', wherein re-protonation of exiting lysine residues prevents channel re-entry, biasing transport in the outward direction.

Lipid droplets (LDs) are universal lipid storage organelles with a core of neutral lipids, such as triacylglycerols, surrounded by a phospholipid monolayer. This unique architecture is generated during LD biogenesis at endoplasmic reticulum (ER) sites marked by Seipin, a conserved membrane protein mutated in lipodystrophy. Here structural, biochemical and molecular dynamics simulation approaches reveal the mechanism of LD formation by the yeast Seipin Sei1 and its membrane partner Ldb16. We show that Sei1 luminal domain assembles a homooligomeric ring, which, in contrast to other Seipins, is unable to concentrate triacylglycerol. Instead, Sei1 positions Ldb16, which concentrates triacylglycerol within the Sei1 ring through critical hydroxyl residues. Triacylglycerol recruitment to the complex is further promoted by Sei1 transmembrane segments, which also control Ldb16 stability. Thus, we propose that LD assembly by the Sei1/Ldb16 complex, and likely other Seipins, requires sequential triacylglycerol-concentrating steps via distinct elements in the ER membrane and lumen., Lipid droplet biogenesis is orchestrated by the conserved membrane protein Seipin via an unknown mechanism. Here, the authors use structural, biochemical and molecular dynamics simulation approaches to reveal the mechanism of lipid droplet formation by the yeast Seipin Sei1 and its partner Ldb16.

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the µ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

Lipids play important modulatory and structural roles for membrane proteins. Molecular dynamics simulations are frequently used to provide insights into the nature of these protein-lipid interactions. Systematic comparative analysis requires tools that provide algorithms for objective assessment of such interactions. We introduce PyLipID, a Python package for the identification and characterization of specific lipid interactions and binding sites on membrane proteins from molecular dynamics simulations. PyLipID uses a community analysis approach for binding site detection, calculating lipid residence times for both the individual protein residues and the detected binding sites. To assist structural analysis, PyLipID produces representative bound lipid poses from simulation data, using a density-based scoring function. To estimate residue contacts robustly, PyLipID uses a dual-cutoff scheme to differentiate between lipid conformational rearrangements while bound from full dissociation events. In addition to the characterization of protein-lipid interactions, PyLipID is applicable to analysis of the interactions of membrane proteins with other ligands. By combining automated analysis, efficient algorithms, and open-source distribution, PyLipID facilitates the systematic analysis of lipid interactions from large simulation data sets of multiple species of membrane proteins.

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-X3-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.FieldsStructural Biology and Biophysics / Computational Biology

Description, The structure of the hexameric fungal proton pump Pma1 reveals its mechanism of autoinhibition and interactions with lipids., The fungal plasma membrane H+-ATPase Pma1 is a vital enzyme, generating a proton-motive force that drives the import of essential nutrients. Autoinhibited Pma1 hexamers in the plasma membrane of starving fungi are activated by glucose signaling and subsequent phosphorylation of the autoinhibitory domain. As related P-type adenosine triphosphatases (ATPases) are not known to oligomerize, the physiological relevance of Pma1 hexamers remained unknown. We have determined the structure of hexameric Pma1 from Neurospora crassa by electron cryo-microscopy at 3.3-Å resolution, elucidating the molecular basis for hexamer formation and autoinhibition and providing a basis for structure-based drug development. Coarse-grained molecular dynamics simulations in a lipid bilayer suggest lipid-mediated contacts between monomers and a substantial protein-induced membrane deformation that could act as a proton-attracting funnel.

Potassium homeostasis is vital for all organisms, but is challenging in single-celled organisms like bacteria and yeast and immobile organisms like plants that constantly need to adapt to changing external conditions. KUP transporters facilitate potassium uptake by the co-transport of protons. Here, we uncover the molecular basis for transport in this widely distributed family. We identify the potassium importer KimA from Bacillus subtilis as a member of the KUP family, demonstrate that it functions as a K+/H+ symporter and report a 3.7 Å cryo-EM structure of the KimA homodimer in an inward-occluded, trans-inhibited conformation. By introducing point mutations, we identify key residues for potassium and proton binding, which are conserved among other KUP proteins., KUP transporters facilitate potassium uptake by the co-transport of protons and are key players in potassium homeostasis. Here authors identify the potassium importer KimA from Bacillus subtilis as a new member of the KUP transporter family and show the cryo-EM structure of KimA in an inward-occluded, trans-inhibited conformation.

Membrane proteins engage in a variety of contacts with their surrounding lipids, but distinguishing between specifically bound lipids, and non‐specific, annular interactions is a challenging problem. Applying native mass spectrometry to three membrane protein complexes with different lipid‐binding properties, we explore the ability of detergents to compete with lipids bound in different environments. We show that lipids in annular positions on the presenilin homologue protease are subject to constant exchange with detergent. By contrast, detergent‐resistant lipids bound at the dimer interface in the leucine transporter show decreased koff rates in molecular dynamics simulations. Turning to the lipid flippase MurJ, we find that addition of the natural substrate lipid‐II results in the formation of a 1:1 protein–lipid complex, where the lipid cannot be displaced by detergent from the highly protected active site. In summary, we distinguish annular from non‐annular lipids based on their exchange rates in solution., Separating the wheat from the chaff: Distinguishing specific and non‐specific lipid binding to membrane proteins is a major challenge. Non‐annular lipid interactions can be identified using a detergent competition assay and native mass spectrometry (MS). MD simulations reveal the location of specific lipid‐binding sites with slow k off rates, confirming the MS data.

The bacterial Sec translocon, SecYEG, associates with accessory proteins YidC and the SecDF-YajC subcomplex to form the bacterial holo-translocon (HTL). The HTL is a dynamic and flexible protein transport machine capable of coordinating protein secretion across the membrane and efficient lateral insertion of nascent membrane proteins. It has been hypothesized that a central lipid core facilitates the controlled passage of membrane proteins into the bilayer, ensuring the efficient formation of their native state. By performing small-angle neutron scattering on protein solubilized in “match-out” deuterated detergent, we have been able to interrogate a “naked” HTL complex, with the scattering contribution of the surrounding detergent micelle rendered invisible. Such an approach has allowed the confirmation of a lipid core within the HTL, which accommodates between 8 and 29 lipids. Coarse-grained molecular dynamics simulations of the HTL also demonstrate a dynamic, central pool of lipids. An opening at this lipid-rich region between YidC and the SecY lateral gate may provide an exit gateway for newly synthesized, correctly oriented, membrane protein helices, or even small bundles of helices, to emerge from the HTL.

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

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

Interactions with lipids can dramatically shape and define the activity of membrane proteins. Here, we describe tools that allow the identification of these interactions using molecular dynamics simulation. Additionally, we provide the details of how to use different methods to probe the affinity of these interactions.

Mitochondrial supercomplexes form around a conserved core of monomeric complex I and dimeric complex III; wherein a subunit of the former, NDUFA11, is conspicuously situated at the interface. We identified nduf-11 (B0491.5) as encoding the Caenorhabditis elegans homologue of NDUFA11. Animals homozygous for a CRISPR-Cas9-generated knockout allele of nduf-11 arrested at the second larval (L2) development stage. Reducing (but not eliminating) expression using RNAi allowed development to adulthood, enabling characterisation of the consequences: destabilisation of complex I and its supercomplexes and perturbation of respiratory function. The loss of NADH dehydrogenase activity was compensated by enhanced complex II activity, with the potential for detrimental reactive oxygen species (ROS) production. Cryo-electron tomography highlighted aberrant morphology of cristae and widening of both cristae junctions and the intermembrane space. The requirement of NDUF-11 for balanced respiration, mitochondrial morphology and development presumably arises due to its involvement in complex I and supercomplex maintenance. This highlights the importance of respiratory complex integrity for health and the potential for its perturbation to cause mitochondrial disease. This article has an associated First Person interview with Amber Knapp-Wilson, joint first author of the paper., Summary: Destabilisation of complex I and its supercomplexes by removal of the assembly factor NDUF-11 severely affects balanced respiration, mitochondrial morphology and, ultimately, whole animal physiology.

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-EM densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities by experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in e.g. cryo-EM derived membrane protein structures surrounded by a plethora of ancillary lipid densities.

KdpFABC, a high-affinity K+ pump, combines the ion channel KdpA and the P-type ATPase KdpB to secure survival at K+ limitation. Here, we apply a combination of cryo-EM, biochemical assays, and MD simulations to illuminate the mechanisms underlying transport and the coupling to ATP hydrolysis. We unambiguously show that ions are transported via an intersubunit tunnel through KdpA and KdpB. At the subunit interface, the tunnel is constricted by a phenylalanine, which, by polarized cation-π stacking, controls K+ entry into the canonical substrate binding site (CBS) of KdpB. Within the CBS, ATPase coupling is mediated by the charge distribution between an aspartate and a lysine. Interestingly, individual elements of the ion translocation mechanism of KdpFABC identified here are conserved among a wide variety of P-type ATPases from different families. This leads us to the hypothesis that KdpB might represent an early descendant of a common ancestor of cation pumps.

The fungal plasma membrane H+-ATPase Pma1 is a vital enzyme, generating a proton-motive force that drives the import of essential nutrients. Auto-inhibited Pma1 hexamers in starving fungi are activated by glucose signalling resulting in phosphorylation of the auto-inhibitory domain. As related P-type ATPases are not known to oligomerise, the physiological relevance of Pma1 hexamers remains unknown. We have determined the structure of hexameric Pma1 from Neurospora crassa by cryo-EM at 3.3 Å resolution, elucidating the molecular basis for hexamer formation and auto-inhibition, and providing a basis for structure-based drug development. Coarse-grained molecular dynamics simulations in a lipid bilayer suggest lipid-mediated contacts between monomers and a substantial protein-induced membrane deformation that could act as a proton-attracting funnel.

Integral membrane proteins are localised and/or regulated by lipids present in the surrounding bilayer. Whilst bacteria such as E. coli have relatively simple membranes when compared to eukaryotic cells, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different E. coli inner membrane proteins. Our data reveals a strong asymmetry between the membrane leaflets, with a marked increase of anionic lipid binding to the inner leaflet regions of membrane proteins, particularly for cardiolipin. From our simulations we identify over 700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high affinity cardiolipin binding site on (bacterial) membrane proteins, paving the way for a heuristic approach to defining more complex protein-lipid interactions.

Overdose deaths from synthetic opioids, such as fentanyl, have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist, that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by a novel binding mode to the μ-opioid receptor (MOPr).By employing coarse-grained molecular dynamics simulations and free energy calculations, we determined the routes by which fentanyl and morphine access the orthosteric pocket of MOPr.Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl takes a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices.This novel lipophilic route may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone.

Mitochondrial super-complexes form around a conserved core of monomeric complex I and dimeric complex III; wherein subunit NDUFA11, of the former, is conspicuously situated at the interface. We identifiedB0491.5(NDUF-11) as theC. eleganshomologue, of which animals homozygous for a CRISPR-Cas9 generated knockout allele arrested at the L2 development stage. Reducing expression by RNAi allowed development to the adult stage, enabling characterisation of the consequences: destabilisation of complex I and its super-complexes, and perturbation of respiratory function. The loss of NADH-dehydrogenase activity is compensated by enhanced complex II activity, resulting in excessive detrimental ROS production. Meanwhile, electron cryo-tomography highlight aberrant cristae morphology and widening of the inter-membrane space and cristae junctions. The requirement of NDUF-11 for balanced respiration, mitochondrial morphology and development highlights the importance of complex I/ super-complex maintenance. Their perturbation by this, or other means, is likely to be the cause of metabolic stress and disease.

Interactions with lipids can dramatically shape and define the activity of membrane proteins. Here, we describe tools that allow the identification of these interactions using molecular dynamics simulation. Additionally, we provide the details of how to use different methods to probe the affinity of these interactions.

International audience; The sarco(endo)plasmic reticulum Ca 2+ -ATPase (SERCA) is a P-type ATPase that transports Ca 2+ from the cytosol into the sarco(endo)plasmic reticulum (SR/ER) lumen, driven by ATP. This primary transport activity depends on tight coupling between movements of the transmembrane helices forming the two Ca 2+ -binding sites and the cytosolic headpiece mediating ATP hydrolysis. We have addressed the molecular basis for this intramolecular communication by analyzing the structure and functional properties of the SERCA mutant E340A. The mutated Glu340 residue is strictly conserved among the P-type ATPase family of membrane transporters and is located at a seemingly strategic position at the interface between the phosphorylation domain and the cytosolic ends of 5 of SERCA’s 10 transmembrane helices. The mutant displays a marked slowing of the Ca 2+ -binding kinetics, and its crystal structure in the presence of Ca 2+ and ATP analog reveals a rotated headpiece, altered connectivity between the cytosolic domains, and an altered hydrogen bonding pattern around residue 340. Supported by molecular dynamics simulations, we conclude that the E340A mutation causes a stabilization of the Ca 2+ sites in a more occluded state, hence displaying slowed dynamics. This finding underpins a crucial role of Glu340 in interdomain communication between the headpiece and the Ca 2+ -binding transmembrane region.

Significance We present a crystal structure, functional data, and molecular dynamics (MD) simulations of the sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) mutant E340A. The mutation slows Ca2+-binding kinetics, and the structural differences between wild type and E340A indicate that the mutation disrupts a central interdomain “communication hub” governing Ca2+ binding/dissociation. MD simulations reveal altered dynamics in regions mediating Ca2+ occlusion, a critical step in SERCA’s alternating access mechanism. The mutation stabilizes a more occluded state of the Ca2+ sites. The strict conservation of Glu340 among P-type ATPases is the result of its critical role in interdomain communication between the cytosolic headpiece and the transmembrane domain, ensuring a delicate balance between dynamics of ion binding, occlusion, and release—key steps in the transport process., The sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) is a P-type ATPase that transports Ca2+ from the cytosol into the sarco(endo)plasmic reticulum (SR/ER) lumen, driven by ATP. This primary transport activity depends on tight coupling between movements of the transmembrane helices forming the two Ca2+-binding sites and the cytosolic headpiece mediating ATP hydrolysis. We have addressed the molecular basis for this intramolecular communication by analyzing the structure and functional properties of the SERCA mutant E340A. The mutated Glu340 residue is strictly conserved among the P-type ATPase family of membrane transporters and is located at a seemingly strategic position at the interface between the phosphorylation domain and the cytosolic ends of 5 of SERCA’s 10 transmembrane helices. The mutant displays a marked slowing of the Ca2+-binding kinetics, and its crystal structure in the presence of Ca2+ and ATP analog reveals a rotated headpiece, altered connectivity between the cytosolic domains, and an altered hydrogen bonding pattern around residue 340. Supported by molecular dynamics simulations, we conclude that the E340A mutation causes a stabilization of the Ca2+ sites in a more occluded state, hence displaying slowed dynamics. This finding underpins a crucial role of Glu340 in interdomain communication between the headpiece and the Ca2+-binding transmembrane region.

The sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) is a P-type ATPase that transports Ca2+from the cytosol into the SR/ER lumen, driven by ATP. This primary transport activity depends on tight coupling between movements of the transmembrane helices forming the two Ca2+binding sites and of the cytosolic headpiece mediating ATP hydrolysis. We have addressed the molecular basis for this intramolecular communication by analyzing the structure and functional properties of the SERCA mutant E340A. The mutated Glu340 residue is strictly conserved among the P-type ATPase family of membrane transporters and is located at a seemingly strategic position at the interface between the phosphorylation domain and the cytosolic ends of five out of SERCA’s ten transmembrane helices. The mutant displays a marked slowing of the Ca2+-binding kinetics, and its crystal structure in the presence of Ca2+and ATP analogue reveals a rotated headpiece, altered connectivity between the cytosolic domains and altered hydrogen bonding pattern around residue 340. Supported by molecular dynamics simulations, we conclude that the E340A mutation causes a stabilization of the Ca2+sites in a more occluded state, hence displaying slowed dynamics. This finding underpins a crucial role of Glu340 in inter-domain communication between the headpiece and the Ca2+-binding transmembrane region.