Detergents with Scalable Properties Identify Noncanonical Lipopolysaccharide Binding to Bacterial Inner Membrane Proteins

Lipopolysaccharide (LPS) is vital for maintaining the outer membrane barrier in Gram-negative bacteria. LPS is also frequently obtained in complex with the inner membrane proteins after detergent purification. The question of whether or not LPS binding to inner membrane proteins not involved in outer membrane biogenesis reflects native lipid environments remains unclear. Here, we leverage the control of the hydrophilic-lipophilic balance and packing parameter concepts to chemically tune detergents that can be used to qualitatively differentiate the degree to which proteins copurify with phospholipids (PLs) and/or LPS. Given the scalable properties of these detergents, we demonstrate a detergent fine-tuning that enables the facile investigation of intact proteins and their complexes with lipids by native mass spectrometry (nMS). We conclude that LPS, a lipid that is believed to be important for outer membranes, can also affect the activity of membrane proteins that are currently not assigned to be involved in outer membrane biogenesis. Our results deliver a scalable detergent chemistry for a streamlined biophysical characterization of protein-lipid interactions, provide a rationale for the high affinity of LPS-protein binding, and identify noncanonical associations between LPS and inner membrane proteins with relevance for membrane biology and antibiotic research.

Peptidoglycan synthesis drives a single population of septal cell wall synthases during division in Bacillus subtilis

Bacterial cell division requires septal peptidoglycan (sPG) synthesis by the divisome complex. Treadmilling of the essential tubulin homologue FtsZ has been implicated in septal constriction, though its precise role remains unclear. Here we used live-cell single-molecule imaging of the divisome transpeptidase PBP2B to investigate sPG synthesis dynamics in Bacillus subtilis. In contrast to previous models, we observed a single population of processively moving PBP2B molecules whose motion is driven by peptidoglycan synthesis and is not associated with FtsZ treadmilling. However, despite the asynchronous motions of PBP2B and FtsZ, a partial dependence of PBP2B processivity on FtsZ treadmilling was observed. Additionally, through single-molecule counting experiments we provide evidence that the divisome synthesis complex is multimeric. Our results support a model for B. subtilis division where a multimeric synthesis complex follows a single track dependent on sPG synthesis whose activity and dynamics are asynchronous with FtsZ treadmilling.

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.

Structure of the native chemotaxis core signaling unit from phage E-protein lysed E. coli cells

IMPORTANCE

Bacterial chemotaxis is a ubiquitous behavior that enables cell movement toward or away from specific chemicals. It serves as an important model for understanding cell sensory signal transduction and motility. Characterization of the molecular mechanisms underlying chemotaxis is of fundamental interest and requires a high-resolution structural picture of the sensing machinery, the chemosensory array. In this study, we combine cryo-electron tomography and molecular simulation to present the complete structure of the core signaling unit, the basic building block of chemosensory arrays, from Escherichia coli. Our results provide new insight into previously poorly-resolved regions of the complex and offer a structural basis for designing new experiments to test mechanistic hypotheses.

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.

Structural basis of peptidoglycan synthesis by E. coli RodA-PBP2 complex

Peptidoglycan (PG) is an essential structural component of the bacterial cell wall that is synthetized during cell division and elongation. PG forms an extracellular polymer crucial for cellular viability, the synthesis of which is the target of many antibiotics. PG assembly requires a glycosyltransferase (GT) to generate a glycan polymer using a Lipid II substrate, which is then crosslinked to the existing PG via a transpeptidase (TP) reaction. A Shape, Elongation, Division and Sporulation (SEDS) GT enzyme and a Class B Penicillin Binding Protein (PBP) form the core of the multi-protein complex required for PG assembly. Here we used single particle cryo-electron microscopy to determine the structure of a cell elongation-specific E. coli RodA-PBP2 complex. We combine this information with biochemical, genetic, spectroscopic, and computational analyses to identify the Lipid II binding sites and propose a mechanism for Lipid II polymerization. Our data suggest a hypothesis for the movement of the glycan strand from the Lipid II polymerization site of RodA towards the TP site of PBP2, functionally linking these two central enzymatic activities required for cell wall peptidoglycan biosynthesis.

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.

PLSCR1 is a cell-autonomous defence factor against SARS-CoV-2 infection

Understanding protective immunity to COVID-19 facilitates preparedness for future pandemics and combats new SARS-CoV-2 variants emerging in the human population. Neutralizing antibodies have been widely studied; however, on the basis of large-scale exome sequencing of protected versus severely ill patients with COVID-19, local cell-autonomous defence is also crucial1-4. Here we identify phospholipid scramblase 1 (PLSCR1) as a potent cell-autonomous restriction factor against live SARS-CoV-2 infection in parallel genome-wide CRISPR-Cas9 screens of human lung epithelia and hepatocytes before and after stimulation with interferon-γ (IFNγ). IFNγ-induced PLSCR1 not only restricted SARS-CoV-2 USA-WA1/2020, but was also effective against the Delta B.1.617.2 and Omicron BA.1 lineages. Its robust activity extended to other highly pathogenic coronaviruses, was functionally conserved in bats and mice, and interfered with the uptake of SARS-CoV-2 in both the endocytic and the TMPRSS2-dependent fusion routes. Whole-cell 4Pi single-molecule switching nanoscopy together with bipartite nano-reporter assays found that PLSCR1 directly targeted SARS-CoV-2-containing vesicles to prevent spike-mediated fusion and viral escape. A PLSCR1 C-terminal β-barrel domain-but not lipid scramblase activity-was essential for this fusogenic blockade. Our mechanistic studies, together with reports that COVID-associated PLSCR1 mutations are found in some susceptible people3,4, identify an anti-coronavirus protein that interferes at a late entry step before viral RNA is released into the host-cell cytosol.

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.

Activator-induced conformational changes regulate division-associated peptidoglycan amidases

AmiA and AmiB are peptidoglycan-hydrolyzing enzymes from Escherichia coli that are required to break the peptidoglycan layer during bacterial cell division and maintain integrity of the cell envelope. In vivo, the activity of AmiA and AmiB is tightly controlled through their interactions with the membrane-bound FtsEX-EnvC complex. Activation of AmiA and AmiB requires access to a groove in the amidase-activating LytM domain of EnvC which is gated by ATP-driven conformational changes in FtsEX-EnvC complex. Here, we present a high-resolution structure of the isolated AmiA protein, confirming that it is autoinhibited in the same manner as AmiB and AmiC, and a complex of the AmiB enzymatic domain bound to the activating EnvC LytM domain. In isolation, the active site of AmiA is blocked by an autoinhibitory helix that binds directly to the catalytic zinc and fills the volume expected to accommodate peptidoglycan binding. In the complex, binding of the EnvC LytM domain induces a conformational change that displaces the amidase autoinhibitory helix and reorganizes the active site for activity. Our structures, together with complementary mutagenesis work, defines the conformational changes required to activate AmiA and/or AmiB through their interaction with their cognate activator EnvC.

Structural basis for membrane attack complex inhibition by CD59

CD59 is an abundant immuno-regulatory receptor that protects human cells from damage during complement activation. Here we show how the receptor binds complement proteins C8 and C9 at the membrane to prevent insertion and polymerization of membrane attack complex (MAC) pores. We present cryo-electron microscopy structures of two inhibited MAC precursors known as C5b8 and C5b9. We discover that in both complexes, CD59 binds the pore-forming β-hairpins of C8 to form an intermolecular β-sheet that prevents membrane perforation. While bound to C8, CD59 deflects the cascading C9 β-hairpins, rerouting their trajectory into the membrane. Preventing insertion of C9 restricts structural transitions of subsequent monomers and indirectly halts MAC polymerization. We combine our structural data with cellular assays and molecular dynamics simulations to explain how the membrane environment impacts the dual roles of CD59 in controlling pore formation of MAC, and as a target of bacterial virulence factors which hijack CD59 to lyse human cells.

Characterization of a TatA/TatB binding site on the TatC component of the Escherichia coli twin arginine translocase

The twin arginine transport (Tat) pathway exports folded proteins across the cytoplasmic membranes of prokaryotes and the thylakoid membranes of chloroplasts. In Escherichia coli and other Gram-negative bacteria, the Tat machinery comprises TatA, TatB and TatC components. A Tat receptor complex, formed from all three proteins, binds Tat substrates, which triggers receptor organization and recruitment of further TatA molecules to form the active Tat translocon. The polytopic membrane protein TatC forms the core of the Tat receptor and harbours two binding sites for the sequence-related TatA and TatB proteins. A 'polar' cluster binding site, formed by TatC transmembrane helices (TMH) 5 and 6 is occupied by TatB in the resting receptor and exchanges for TatA during receptor activation. The second binding site, lying further along TMH6, is occupied by TatA in the resting state, but its functional relevance is unclear. Here we have probed the role of this second binding site through a programme of random and targeted mutagenesis. Characterization of three stably produced TatC variants, P221R, M222R and L225P, each of which is inactive for protein transport, demonstrated that the substitutions did not affect assembly of the Tat receptor. Moreover, the substitutions that we analysed did not abolish TatA or TatB binding to either binding site. Using targeted mutagenesis we introduced bulky substitutions into the TatA binding site. Molecular dynamics simulations and crosslinking analysis indicated that TatA binding at this site was substantially reduced by these amino acid changes, but TatC retained function. While it is not clear whether TatA binding at the TMH6 site is essential for Tat activity, the isolation of inactivating substitutions indicates that this region of the protein has a critical function.

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 TGES₁₆₂ motif in the A domain of the pump subunit KdpB (KdpB_S162-P). Here, we explore the structural basis of inhibition by KdpB_S162 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 KdpB_S162 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.

The dynamic interplay of PIP2 and ATP in the regulation of the KATP channel

ATP-sensitive potassium (KATP ) channels couple the intracellular ATP concentration to insulin secretion. KATP channel activity is inhibited by ATP binding to the Kir6.2 tetramer and activated by phosphatidylinositol 4,5-bisphosphate (PIP2 ). Here, we use molecular dynamics simulation, electrophysiology and fluorescence spectroscopy to show that ATP and PIP2 occupy different binding pockets that share a single amino acid residue, K39. When both ligands are present, simulations suggest that K39 shows a greater preference to co-ordinate with PIP2 than with ATP. They also predict that a neonatal diabetes mutation at K39 (K39R) increases the number of hydrogen bonds formed between K39 and PIP2 , potentially accounting for the reduced ATP inhibition observed in electrophysiological experiments. Our work suggests that PIP2 and ATP interact allosterically to regulate KATP channel activity. KEY POINTS: The KATP channel is activated by the binding of phosphatidylinositol 4,5-bisphosphate (PIP2 ) lipids and inactivated by the binding of ATP. K39 has the potential to bind to both PIP2 and ATP. A mutation to this residue (K39R) results in neonatal diabetes. This study uses patch-clamp fluorometry, electrophysiology and molecular dynamics simulation. We show that PIP2 competes with ATP for K39, and this reduces channel inhibition by ATP. We show that K39R increases channel affinity to PIP2 by increasing the number of hydrogen bonds with PIP2 , when compared with the wild-type K39. This therefore decreases KATP channel inhibition by ATP.

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.

Conformational dynamics of the membrane enzyme LspA upon antibiotic and substrate binding

Lipoprotein signal peptidase (LspA) is an aspartyl protease that cleaves the transmembrane helix signal peptide of lipoproteins as part of the lipoprotein-processing pathway. Members of this pathway are excellent targets for the development of antibiotic therapeutics because they are essential in Gram-negative bacteria, are important for virulence in Gram-positive bacteria, and may not develop antibiotic resistance. Here, we report the conformational dynamics of LspA in the apo state and bound to the antibiotic globomycin determined using molecular dynamics simulations and electron paramagnetic resonance. The periplasmic helix fluctuates on the nanosecond timescale and samples unique conformations in the different states. In the apo state, the dominant conformation is the most closed and occludes the charged active site from the lipid bilayer. With antibiotic bound there are multiple binding modes with the dominant conformation of the periplasmic helix in a more open conformation. The different conformations observed in both bound and apo states indicate a flexible and adaptable active site, which explains how LspA accommodates and processes such a variety of substrates.

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 (C₅₅-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.

Bacterial cell wall enzymes and their precursors are critical targets for antibiotic development. Here, the authors investigate several biosynthetic enzymes with their substrates and show that the passage of substrates and products in the pathway is controlled by their relative binding affinities.

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.

Structural basis of lipopolysaccharide maturation by the O-antigen ligase

The outer membrane of Gram-negative bacteria has an external leaflet that is largely composed of lipopolysaccharide, which provides a selective permeation barrier, particularly against antimicrobials1. The final and crucial step in the biosynthesis of lipopolysaccharide is the addition of a species-dependent O-antigen to the lipid A core oligosaccharide, which is catalysed by the O-antigen ligase WaaL2. Here we present structures of WaaL from Cupriavidus metallidurans, both in the apo state and in complex with its lipid carrier undecaprenyl pyrophosphate, determined by single-particle cryo-electron microscopy. The structures reveal that WaaL comprises 12 transmembrane helices and a predominantly α-helical periplasmic region, which we show contains many of the conserved residues that are required for catalysis. We observe a conserved fold within the GT-C family of glycosyltransferases and hypothesize that they have a common mechanism for shuttling the undecaprenyl-based carrier to and from the active site. The structures, combined with genetic, biochemical, bioinformatics and molecular dynamics simulation experiments, offer molecular details on how the ligands come in apposition, and allows us to propose a mechanistic model for catalysis. Together, our work provides a structural basis for lipopolysaccharide maturation in a member of the GT-C superfamily of glycosyltransferases.

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.

Interrogation of the Pathogen Box reveals small molecule ligands against the mycobacterial trehalose transporter LpqY-SugABC

Tuberculosis, caused by Mycobacterium tuberculosis, claims ∼1.5 million lives annually. Effective chemotherapy is essential to control TB, however the emergence of drug-resistant strains of TB have seriously threatened global attempts to control and eradicate this deadly pathogen. Trehalose recycling via the LpqY-SugABC importer is essential for the virulence and survival of Mtb and inhibiting or hijacking this transport system is an attractive approach for the development of novel anti-tubercular and diagnostic agents. Therefore, we interrogated the drug-like compounds in the open-source Medicines for Malaria Pathogen Box and successfully identified seven compounds from the TB, kinetoplastids and reference compound disease sets that recognise LpqY. The molecules have diverse chemical scaffolds, are not specific trehalose analogues, and may be used as novel templates to facilitate the development of therapeutics that kill Mtb with a novel mechanism of action via the mycobacterial trehalose LpqY-SugABC transport system.

Interrogation of the Pathogen Box identified diverse chemical scaffolds against the mycobacterial trehalose transporter.

NaViA: a program for the visual analysis of complex mass spectra

MOTIVATION

Native mass spectrometry is now a well-established method for the investigation of protein complexes, specifically their subunit stoichiometry and ligand binding properties. Recent advances allowing the analysis of complex mixtures lead to an increasing diversity and complexity in the spectra obtained. These spectra can be time-consuming to tackle through manual assignment and challenging for automated approaches.

RESULTS

Native Mass Spectrometry Visual Analyser is a web-based tool to augment the manual process of peak assignment. In addition to matching masses to the stoichiometry of its component subunits, it allows raw data processing, assignment and annotation and permits mass spectra to be shared with their respective interpretation.

AVAILABILITY AND IMPLEMENTATION

NaViA is open-source and can be accessed online under https://navia.ms. The source code and documentation can be accessed at https://github.com/d-que/navia, under the BSD 2-Clause licence.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High-resolution mapping of metal ions reveals principles of surface layer assembly in Caulobacter crescentus cells

Surface layers (S-layers) are proteinaceous crystalline coats that constitute the outermost component of most prokaryotic cell envelopes. In this study, we have investigated the role of metal ions in the formation of the Caulobacter crescentus S-layer using high-resolution structural and cell biology techniques, as well as molecular simulations. Utilizing optical microscopy of fluorescently tagged S-layers, we show that calcium ions facilitate S-layer lattice formation and cell-surface binding. We report all-atom molecular dynamics simulations of the S-layer lattice, revealing the importance of bound metal ions. Finally, using electron cryomicroscopy and long-wavelength X-ray diffraction experiments, we mapped the positions of metal ions in the S-layer at near-atomic resolution, supporting our insights from the cellular and simulations data. Our findings contribute to the understanding of how C. crescentus cells form a regularly arranged S-layer on their surface, with implications on fundamental S-layer biology and the synthetic biology of self-assembling biomaterials.

•

Live imaging shows Ca²⁺-dependent expansion of the C. crescentus S-layer

•

Molecular simulations reveal Ca²⁺-binding properties of the S-layer

•

Ca²⁺ ion mapping in three-dimensional crystals using in-vacuum X-ray anomalous diffraction

•

Ca²⁺ replacement by Ho³⁺ allows cryo-EM mapping of heavy metals

CG2AT2: an Enhanced Fragment-Based Approach for Serial Multi-scale Molecular Dynamics Simulations

Coarse-grained molecular dynamics provides a means for simulating the assembly and interactions of macromolecular complexes at a reduced level of representation, thereby allowing both longer timescale and larger sized simulations. Here, we describe an enhanced fragment-based protocol for converting macromolecular complexes from coarse-grained to atomistic resolution, for further refinement and analysis. While the focus is upon systems that comprise an integral membrane protein embedded in a phospholipid bilayer, the technique is also suitable for membrane-anchored and soluble protein/nucleotide complexes. Overall, this provides a method for generating an accurate and well-equilibrated atomic-level description of a macromolecular complex. The approach is evaluated using a diverse test set of 11 system configurations of varying size and complexity. Simulations are assessed in terms of protein stereochemistry, conformational drift, lipid/protein interactions, and lipid dynamics.

Relative Affinities of Protein-Cholesterol Interactions from Equilibrium Molecular Dynamics Simulations

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-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using 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 an array of ancillary lipid densities.

Mechanism of lipid droplet formation by the yeast Sei1/Ldb16 Seipin complex

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.

From Bench to Biomolecular Simulation: Phospholipid Modulation of Potassium Channels

Potassium (K+) ion channels are crucial in numerous cellular processes as they hyperpolarise a cell through K+ conductance, returning a cell to its resting potential. K+ channel mutations result in multiple clinical complications such as arrhythmia, neonatal diabetes and migraines. Since 1995, the regulation of K+ channels by phospholipids has been heavily studied using a range of interdisciplinary methods such as cellular electrophysiology, structural biology and computational modelling. As a result, K+ channels are model proteins for the analysis of protein-lipid interactions. In this review, we will focus on the roles of lipids in the regulation of K+ channels, and how atomic-level structures, along with experimental techniques and molecular simulations, have helped guide our understanding of the importance of phospholipid interactions.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

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.

Switching Cytolytic Nanopores into Antimicrobial Fractal Ruptures by a Single Side Chain Mutation

Disruption of cell membranes is a fundamental host defense response found in virtually all forms of life. The molecular mechanisms vary but generally lead to energetically favored circular nanopores. Here, we report an elaborate fractal rupture pattern induced by a single side-chain mutation in ultrashort (8-11-mers) helical peptides, which otherwise form transmembrane pores. In contrast to known mechanisms, this mode of membrane disruption is restricted to the upper leaflet of the bilayer where it exhibits propagating fronts of peptide-lipid interfaces that are strikingly similar to viscous instabilities in fluid flow. The two distinct disruption modes, pores and fractal patterns, are both strongly antimicrobial, but only the fractal rupture is nonhemolytic. The results offer wide implications for elucidating differential membrane targeting phenomena defined at the nanoscale.

Alternative Architecture of the E. coli Chemosensory Array

Chemotactic responses in motile bacteria are the result of sophisticated signal transduction by large, highly organized arrays of sensory proteins. Despite tremendous progress in the understanding of chemosensory array structure and function, a structural basis for the heightened sensitivity of networked chemoreceptors is not yet complete. Here, we present cryo-electron tomography visualisations of native-state chemosensory arrays in E. coli minicells. Strikingly, these arrays appear to exhibit a p2-symmetric array architecture that differs markedly from the p6-symmetric architecture previously described in E. coli. Based on this data, we propose molecular models of this alternative architecture and the canonical p6-symmetric assembly. We evaluate our observations and each model in the context of previously published data, assessing the functional implications of an alternative architecture and effects for future studies.

Membrane Binding of Antimicrobial Peptides Is Modulated by Lipid Charge Modification

Peptide interactions with lipid bilayers play a key role in a range of biological processes and depend on electrostatic interactions between charged amino acids and lipid headgroups. Antimicrobial peptides (AMPs) initiate the killing of bacteria by binding to and destabilizing their membranes. The multiple peptide resistance factor (MprF) provides a defense mechanism for bacteria against a broad range of AMPs. MprF reduces the negative charge of bacterial membranes through enzymatic conversion of the anionic lipid phosphatidyl glycerol (PG) to either zwitterionic alanyl-phosphatidyl glycerol (Ala-PG) or cationic lysyl-phosphatidyl glycerol (Lys-PG). The resulting change in the membrane charge is suggested to reduce the binding of AMPs to membranes, thus impeding downstream AMP activity. Using coarse-grained molecular dynamics to investigate the effects of these modified lipids on AMP binding to model membranes, we show that AMPs have substantially reduced affinity for model membranes containing Ala-PG or Lys-PG. More than 5000 simulations in total are used to define the relationship between lipid bilayer composition, peptide sequence (using five different membrane-active peptides), and peptide binding to membranes. The degree of interaction of a peptide with a membrane correlates with the membrane surface charge density. Free energy profile (potential of mean force) calculations reveal that the lipid modifications due to MprF alter the energy barrier to peptide helix penetration of the bilayer. These results will offer a guide to the design of novel peptides, which addresses the issue of resistance via MprF-mediated membrane modification.

Structural basis of trehalose recognition by the mycobacterial LpqY-SugABC transporter

The Mycobacterium tuberculosis (Mtb) LpqY-SugABC ATP-binding cassette transporter is a recycling system that imports trehalose released during remodeling of the Mtb cell-envelope. As this process is essential for the virulence of the Mtb pathogen, it may represent an important target for tuberculosis drug and diagnostic development, but the transporter specificity and molecular determinants of substrate recognition are unknown. To address this, we have determined the structural and biochemical basis of how mycobacteria transport trehalose using a combination of crystallography, saturation transfer difference NMR, molecular dynamics, site-directed mutagenesis, biochemical/biophysical assays, and the synthesis of trehalose analogs. This analysis pinpoints key residues of the LpqY substrate binding lipoprotein that dictate substrate-specific recognition and has revealed which disaccharide modifications are tolerated. These findings provide critical insights into how the essential Mtb LpqY-SugABC transporter reuses trehalose and modified analogs and specifies a framework that can be exploited for the design of new antitubercular agents and/or diagnostic tools.

Structures of the stator complex that drives rotation of the bacterial flagellum

The bacterial flagellum is the prototypical protein nanomachine and comprises a rotating helical propeller attached to a membrane-embedded motor complex. The motor consists of a central rotor surrounded by stator units that couple ion flow across the cytoplasmic membrane to generate torque. Here, we present the structures of the stator complexes from Clostridium sporogenes, Bacillus subtilis and Vibrio mimicus, allowing interpretation of the extensive body of data on stator mechanism. The structures reveal an unexpected asymmetric A5B2 subunit assembly where the five A subunits enclose the two B subunits. Comparison to structures of other ion-driven motors indicates that this A5B2 architecture is fundamental to bacterial systems that couple energy from ion flow to generate mechanical work at a distance and suggests that such events involve rotation in the motor structures.

The SARS-COV-2 Spike Protein Binds Sialic Acids and Enables Rapid Detection in a Lateral Flow Point of Care Diagnostic Device

There is an urgent need to understand the behavior of the novel coronavirus (SARS-COV-2), which is the causative agent of COVID-19, and to develop point-of-care diagnostics. Here, a glyconanoparticle platform is used to discover that N-acetyl neuraminic acid has affinity toward the SARS-COV-2 spike glycoprotein, demonstrating its glycan-binding function. Optimization of the particle size and coating enabled detection of the spike glycoprotein in lateral flow and showed selectivity over the SARS-COV-1 spike protein. Using a virus-like particle and a pseudotyped lentivirus model, paper-based lateral flow detection was demonstrated in under 30 min, showing the potential of this system as a low-cost detection platform.

Evaluating inositol phospholipid interactions with inward rectifier potassium channels and characterising their role in disease

Membrane proteins are frequently modulated by specific protein-lipid interactions. The activation of human inward rectifying potassium (hKir) channels by phosphoinositides (PI) has been well characterised. Here, we apply a coarse-grained molecular dynamics free-energy perturbation (CG-FEP) protocol to capture the energetics of binding of PI lipids to hKir channels. By using either a single- or multi-step approach, we establish a consistent value for the binding of PIP2 to hKir channels, relative to the binding of the bulk phosphatidylcholine phospholipid. Furthermore, by perturbing amino acid side chains on hKir6.2, we show that the neonatal diabetes mutation E179K increases PIP2 affinity, while the congenital hyperinsulinism mutation K67N results in a reduced affinity. We show good agreement with electrophysiological data where E179K exhibits a reduction in neomycin sensitivity, implying that PIP2 binds more tightly E179K channels. This illustrates the application of CG-FEP to compare affinities between lipid species, and for annotating amino acid residues.

Conformational dynamics of a G protein-coupled receptor helix 8 in lipid membranes

G protein-coupled receptors (GPCRs) are the largest and pharmaceutically most important class of membrane proteins encoded in the human genome, characterized by a seven-transmembrane helix architecture and a C-terminal amphipathic helix 8 (H8). In a minority of GPCR structures solved to date, H8 either is absent or adopts an unusual conformation. The controversial existence of H8 of the class A GPCR neurotensin receptor 1 (NTS1) has been examined here for the nonthermostabilized receptor in a functionally supporting membrane environment using electron paramagnetic resonance, molecular dynamics simulations, and circular dichroism. Lipid-protein interactions with phosphatidylserine and phosphatidylethanolamine lipids, in particular, stabilize the residues 374 to 390 of NTS1 into forming a helix. Furthermore, introduction of a helix-breaking proline residue in H8 elicited an increase in ß-arrestin-NTS1 interactions observed in pull-down assays, suggesting that the structure and/or dynamics of H8 might play an important role in GPCR signaling.

New insights into KATP channel gene mutations and neonatal diabetes mellitus

The ATP-sensitive potassium channel (K_ATP channel) couples blood levels of glucose to insulin secretion from pancreatic β-cells. K_ATP channel closure triggers a cascade of events that results in insulin release. Metabolically generated changes in the intracellular concentrations of adenosine nucleotides are integral to this regulation, with ATP and ADP closing the channel and MgATP and MgADP increasing channel activity. Activating mutations in the genes encoding either of the two types of K_ATP channel subunit (Kir6.2 and SUR1) result in neonatal diabetes mellitus, whereas loss-of-function mutations cause hyperinsulinaemic hypoglycaemia of infancy. Sulfonylurea and glinide drugs, which bind to SUR1, close the channel through a pathway independent of ATP and are now the primary therapy for neonatal diabetes mellitus caused by mutations in the genes encoding K_ATP channel subunits. Insight into the molecular details of drug and nucleotide regulation of channel activity has been illuminated by cryo-electron microscopy structures that reveal the atomic-level organization of the K_ATP channel complex. Here we review how these structures aid our understanding of how the various mutations in the genes encoding Kir6.2 (KCNJ11) and SUR1 (ABCC8) lead to a reduction in ATP inhibition and thereby neonatal diabetes mellitus. We also provide an update on known mutations and sulfonylurea therapy in neonatal diabetes mellitus.

Atomistic mechanism of transmembrane helix association

Transmembrane helix association is a fundamental step in the folding of helical membrane proteins. The prototypical example of this association is formation of the glycophorin dimer. While its structure and stability have been well-characterized experimentally, the detailed assembly mechanism is harder to obtain. Here, we use all-atom simulations within phospholipid membrane to study glycophorin association. We find that initial association results in the formation of a non-native intermediate, separated by a significant free energy barrier from the dimer with a native binding interface. We have used transition-path sampling to determine the association mechanism. We find that the mechanism of the initial bimolecular association to form the intermediate state can be mediated by many possible contacts, but seems to be particularly favoured by formation of non-native contacts between the C-termini of the two helices. On the other hand, the contacts which are key to determining progression from the intermediate to the native state are those which define the native binding interface, reminiscent of the role played by native contacts in determining folding of globular proteins. As a check on the simulations, we have computed association and dissociation rates from the transition-path sampling. We obtain results in reasonable accord with available experimental data, after correcting for differences in native state stability. Our results yield an atomistic description of the mechanism for a simple prototype of helical membrane protein folding.

Many important cellular functions are performed by membrane proteins, and in particular by association of proteins via transmembrane helices. However, the mechanism of how the helices associate has been challenging to study, by either experiment or simulation. Here, we use advanced molecular simulation methods to overcome the slow time scales involved in helix association and dissociation and obtain a view of the association mechanism in atomic detail. We show that association occurs via an initially non-native dimer, before proceeding to the native state, and we validate our results by comparison to available experimental kinetic data. Our methods will also aid in the study of the assembly mechanism of larger transmembrane proteins via molecular simulation.