Synthesis of lipid-linked precursors of the bacterial cell wall is governed by a feedback control mechanism in Pseudomonas aeruginosa.

Many bacterial surface glycans such as the peptidoglycan (PG) cell wall are built from monomeric units linked to a polyprenyl lipid carrier. How this limiting carrier is distributed among competing pathways has remained unclear. Here we describe the isolation of hyperactive variants of Pseudomonas aeruginosa MraY, the enzyme that forms the first lipid-linked PG precursor. These variants result in the elevated production of the final PG precursor lipid II in cells and are hyperactive in vitro. The activated MraY variants have substitutions that map to a cavity on the extracellular side of the dimer interface, far from the active site. Our structural and molecular dynamics results suggest that this cavity is a binding site for externalized lipid II. Overall, our results support a model in which excess externalized lipid II allosterically inhibits MraY, providing a feedback mechanism that prevents the sequestration of lipid carrier in the PG biogenesis pathway.

Biophysical basis of filamentous phage tactoid-mediated antibiotic tolerance in P. aeruginosa.

Inoviruses are filamentous phages infecting numerous prokaryotic phyla. Inoviruses can self-assemble into mesoscale structures with liquid-crystalline order, termed tactoids, which protect bacterial cells in Pseudomonas aeruginosa biofilms from antibiotics. Here, we investigate the structural, biophysical, and protective properties of tactoids formed by the P. aeruginosa phage Pf4 and Escherichia coli phage fd. A cryo-EM structure of the capsid from fd revealed distinct biochemical properties compared to Pf4. Fd and Pf4 formed tactoids with different morphologies that arise from differing phage geometries and packing densities, which in turn gave rise to different tactoid emergent properties. Finally, we showed that tactoids formed by either phage protect rod-shaped bacteria from antibiotic treatment, and that direct association with a tactoid is required for protection, demonstrating the formation of a diffusion barrier by the tactoid. This study provides insights into how filamentous molecules protect bacteria from extraneous substances in biofilms and in host-associated infections.

LipIDens: simulation assisted interpretation of lipid densities in cryo-EM structures of membrane proteins.

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. We demonstrate this by application of our open-source LipIDens code to ten diverse membrane protein structures which exhibit lipid-like densities.

LptM promotes oxidative maturation of the lipopolysaccharide translocon by substrate binding mimicry.

Insertion of lipopolysaccharide (LPS) into the bacterial outer membrane (OM) is mediated by a druggable OM translocon consisting of a ß-barrel membrane protein, LptD, and a lipoprotein, LptE. The ß-barrel assembly machinery (BAM) assembles LptD together with LptE at the OM. In the enterobacterium Escherichia coli, formation of two native disulfide bonds in LptD controls translocon activation. Here we report the discovery of LptM (formerly YifL), a lipoprotein conserved in Enterobacteriaceae, that assembles together with LptD and LptE at the BAM complex. LptM stabilizes a conformation of LptD that can efficiently acquire native disulfide bonds, whereas its inactivation makes disulfide bond isomerization by DsbC become essential for viability. Our structural prediction and biochemical analyses indicate that LptM binds to sites in both LptD and LptE that are proposed to coordinate LPS insertion into the OM. These results suggest that, by mimicking LPS binding, LptM facilitates oxidative maturation of LptD, thereby activating the LPS translocon.

A feedback control mechanism governs the synthesis of lipid-linked precursors of the bacterial cell wall.

Many bacterial surface glycans such as the peptidoglycan (PG) cell wall, O-antigens, and capsules are built from monomeric units linked to a polyprenyl lipid carrier. How this limiting lipid carrier is effectively distributed among competing pathways has remained unclear for some time. Here, we describe the isolation and characterization of hyperactive variants of Pseudomonas aeruginosa MraY, the essential and conserved enzyme catalyzing the formation of the first lipid-linked PG precursor called lipid I. These variants result in the elevated production of the final PG precursor lipid II in cells and are hyperactive in a purified system. Amino acid substitutions within the activated MraY variants unexpectedly map to a cavity on the extracellular side of the dimer interface, far from the active site. Our structural evidence and molecular dynamics simulations suggest that the cavity is a binding site for lipid II molecules that have been transported to the outer leaflet of the membrane. Overall, our results support a model in which excess externalized lipid II allosterically inhibits MraY, providing a feedback mechanism to prevent the sequestration of lipid carrier in the PG biogenesis pathway. MraY belongs to the broadly distributed polyprenyl-phosphate N-acetylhexosamine 1-phosphate transferase (PNPT) superfamily of enzymes. We therefore propose that similar feedback mechanisms may be widely employed to coordinate precursor supply with demand by polymerases, thereby optimizing the partitioning of lipid carriers between competing glycan biogenesis pathways.

The energetics and ion coupling of cholesterol transport through Patched1.

Patched1 (PTCH1) is a tumor suppressor protein of the mammalian Hedgehog (HH) signaling pathway, implicated in embryogenesis and tissue homeostasis. PTCH1 inhibits the G protein-coupled receptor Smoothened (SMO) via a debated mechanism involving modulating ciliary cholesterol accessibility. Using extensive molecular dynamics simulations and free energy calculations to evaluate cholesterol transport through PTCH1, we find an energetic barrier of ~15 to 20 kilojoule per mole for cholesterol export. In silico data are coupled to in vivo biochemical assays of PTCH1 mutants to probe coupling between cation binding sites, transmembrane motions, and PTCH1 activity. Using complementary simulations of Dispatched1, 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 one to three Na+ or two to three K+ couple to cholesterol export, and provide the first molecular description of transitions between distinct transport states.

Cyclic di-AMP traps proton-coupled K+ transporters of the KUP family in an inward-occluded conformation.

Cyclic di-AMP is the only known essential second messenger in bacteria and archaea, regulating different proteins indispensable for numerous physiological processes. In particular, it controls various potassium and osmolyte transporters involved in osmoregulation. In Bacillus subtilis, the K+/H+ symporter KimA of the KUP family is inactivated by c-di-AMP. KimA sustains survival at potassium limitation at low external pH by mediating potassium ion uptake. However, at elevated intracellular K+ concentrations, further K+ accumulation would be toxic. In this study, we reveal the molecular basis of how c-di-AMP binding inhibits KimA. We report cryo-EM structures of KimA with bound c-di-AMP in detergent solution and reconstituted in amphipols. By combining structural data with functional assays and molecular dynamics simulations we reveal how c-di-AMP modulates transport. We show that an intracellular loop in the transmembrane domain interacts with c-di-AMP bound to the adjacent cytosolic domain. This reduces the mobility of transmembrane helices at the cytosolic side of the K+ binding site and therefore traps KimA in an inward-occluded conformation.

The Energetics and Ion Coupling of Cholesterol Transport Through Patched1.

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 -1 for 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 silico data are coupled to in vivo biochemical 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.

Comparative Study of Receptor-, Receptor State-, and Membrane-Dependent Cholesterol Binding Sites in A2A and A1 Adenosine Receptors Using Coarse-Grained Molecular Dynamics Simulations.

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 the A1 adenosine receptor (A1R) and predict new cholesterol binding sites indicating amino acid residues with a high residence time in three biologically relevant membranes. Compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)-cholesterol and POPC-phosphatidylinositol-bisphosphate (PIP2)-cholesterol, the plasma mimetic membrane best described the cholesterol binding sites previously detected for the inactive state of A2AR and revealed the binding sites with long-lasting amino acid residues. We observed that using the plasma mimetic membrane and plotting residues with cholesterol residence time ≥2 µs, our CG MD simulations captured most obviously the cholesterol-protein interactions. For the inactive A2AR, we identified one more binding site in which cholesterol is bound to residues with a long residence time compared to the previously detected, for the active A1R, three binding sites, and for the inactive A1R, two binding sites. We calculated that for the active states, cholesterol binds to residues with a much longer residence time compared to the inactive state for both A2AR and A1R. The stability of the identified binding sites to A1R or A2AR with CG MD simulations was additionally investigated with potential of mean force calculations using umbrella sampling. We observed that the binding sites with residues to which cholesterol has a long residence time in A2AR have shallow binding free energy minima compared to the related binding sites in A1R, suggesting a stronger binding for cholesterol to A1R. The differences in binding sites in which cholesterol is stabilized and interacts with residues with a long residence time between active and inactive states of A1R and A2AR can be important for differences in functional activity and orthosteric agonist or antagonist affinity and can be used for the design of allosteric modulators, which can bind through lipid pathways. We observed a stronger binding for cholesterol to A1R (i.e., generally higher association rates) compared to A2AR, which remains to be demonstrated. For the active states, cholesterol binds to residues with much longer residence times compared to the inactive state for both A2AR and A1R. Taken together, binding sites of active A1R may be considered as promising allosteric targets.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane.

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.

Cardiolipin, and not monolysocardiolipin, preferentially binds to the interface of complexes III and IV.

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.

Inhibited KdpFABC transitions into an E1 off-cycle state.

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.

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology.

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.

Seipin concentrates distinct neutral lipids via interactions with their acyl chain carboxyl esters.

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

Rate-limiting transport of positively charged arginine residues through the Sec-machinery is integral to the mechanism of protein secretion.

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.

Peptidoglycan biosynthesis is driven by lipid transfer along enzyme-substrate affinity gradients.

Maintenance of bacterial cell shape and resistance to osmotic stress by the peptidoglycan (PG) renders PG biosynthetic enzymes and precursors attractive targets for combating bacterial infections. Here, by applying native mass spectrometry, we elucidate the effects of lipid substrates on the PG membrane enzymes MraY, MurG, and MurJ. We show that dimerization of MraY is coupled with binding of the carrier lipid substrate undecaprenyl phosphate (C55-P). Further, we demonstrate the use of native MS for biosynthetic reaction monitoring and find that the passage of substrates and products is controlled by the relative binding affinities of the different membrane enzymes. Overall, we provide a molecular view of how PG membrane enzymes convey lipid precursors through favourable binding events and highlight possible opportunities for intervention.

Structure, substrate recognition and initiation of hyaluronan synthase.

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.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

PyLipID: A Python Package for Analysis of Protein-Lipid Interactions from Molecular Dynamics Simulations.

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.