Outer Membrane Remodeling: The Structural Dynamics and Electrostatics of Rough Lipopolysaccharide Chemotypes

Martini-3 Coarse-Grained Models for the Bacterial Lipopolysaccharide Outer Membrane of Escherichia coli

The outer lipopolysaccharide (LPS) membrane of Gram-negative bacteria forms the main barrier for transport of antimicrobial molecules into the bacterial cell. In this study we develop coarse-grained models for the outer membrane of Escherichia coli in the Martini-3 framework. The coarse-grained model force field was parametrized and validated using all-atom simulations of symmetric membranes of lipid A and rough LPS as well as a complete asymmetric membrane of LPS with the O-antigen. The bonded parameters were obtained using an iterative refinement procedure with target bonded distributions obtained from all-atom simulations. The membrane thickness, area of the LPS, and density distributions for the different regions as well as the water and ion densities in Martini-3 simulations show excellent agreement with the all-atom data. Additionally the solvent accessible surface area for individual molecules in water was found to be in good agreement. The binding of calcium ions with phosphate and carboxylate moieties of LPS is accurately captured in the Martini-3 model, indicative of the integrity of the highly negatively charged LPS molecules in the outer membranes of Gram-negative bacteria. The melting transition of the coarse-grained lipid A membrane model was found to occur between 300 and 310 K, and the model captured variations in area per LPS, order parameter, and membrane thickness across the melting transition. Our study reveals that the proposed Martini-3 models for LPS are able to capture the physicochemical balance of the complex sugar architecture of the outer membrane of Escherichia coli. The coarse-grained models developed in this study would be useful for determining membrane protein interactions and permeation of potential antimicrobials through bacterial membranes at mesoscopic spatial and temporal scales.

Lipidome modulation by dietary omega-3 polyunsaturated fatty acid supplementation or selective soluble epoxide hydrolase inhibition suppresses rough LPS-accelerated glomerulonephritis in lupus-prone mice

Introduction

Lipopolysaccharide (LPS)-accelerated autoimmune glomerulonephritis (GN) in NZBWF1 mice is a preclinical model potentially applicable for investigating lipidome-modulating interventions against lupus. LPS can be expressed as one of two chemotypes: smooth LPS (S-LPS) or rough LPS (R-LPS) which is devoid of O-antigen polysaccharide sidechain. Since these chemotypes differentially affect toll-like receptor 4 (TLR4)-mediated immune cell responses, these differences may influence GN induction.

Methods

We initially compared the effects of subchronic intraperitoneal (i.p.) injection for 5 wk with 1) Salmonella S-LPS, 2) Salmonella R-LPS, or 3) saline vehicle (VEH) (Study 1) in female NZBWF1 mice. Based on the efficacy of R-LPS in inducing GN, we next used it to compare the impact of two lipidome-modulating interventions, ω-3 polyunsaturated fatty acid (PUFA) supplementation and soluble epoxide hydrolase (sEH) inhibition, on GN (Study 2). Specifically, effects of consuming ω-3 docosahexaenoic acid (DHA) (10 g/kg diet) and/or the sEH inhibitor 1-(4-trifluoro-methoxy-phenyl)-3-(1-propionylpiperidin-4-yl) urea (TPPU) (22.5 mg/kg diet ≈ 3 mg/kg/day) on R-LPS triggering were compared.

Results

In Study 1, R-LPS induced robust elevations in blood urea nitrogen, proteinuria, and hematuria that were not evident in VEH- or S-LPS-treated mice. R-LPS-treated mice further exhibited kidney histopathology including robust hypertrophy, hyperplasia, thickened membranes, lymphocytic accumulation containing B and T cells, and glomerular IgG deposition consistent with GN that was not evident in VEH- or SLPS-treated groups. R-LPS but not S-LPS induced spleen enlargement with lymphoid hyperplasia and inflammatory cell recruitment in the liver. In Study 2, resultant blood fatty acid profiles and epoxy fatty acid concentrations reflected the anticipated DHA- and TPPU-mediated lipidome changes, respectively. The relative rank order of R-LPS-induced GN severity among groups fed experimental diets based on proteinuria, hematuria, histopathologic scoring, and glomerular IgG deposition was: VEH/CON< R-LPS/DHA ≈ R-LPS/TPPU<<< R-LPS/TPPU+DHA ≈ R-LPS/CON. In contrast, these interventions had modest-to- negligible effects on R-LPS-induced splenomegaly, plasma antibody responses, liver inflammation, and inflammation-associated kidney gene expression.

Discussion

We show for the first time that absence of O-antigenic polysaccharide in R-LPS is critical to accelerated GN in lupus-prone mice. Furthermore, intervention by lipidome modulation through DHA feeding or sEH inhibition suppressed R-LPS-induced GN; however, these ameliorative effects were greatly diminished upon combining the treatments.

Genes mcr improve the intestinal fitness of pathogenic E. coli and balance their lifestyle to commensalism

BACKGROUND

The plasmid-mediated resistance gene mcr-1 confers colistin resistance in Escherichia coli and paves the way for the evolution to pan-drug resistance. We investigated the impact of mcr-1 in gut colonization in the absence of antibiotics using isogenic E. coli strains transformed with a plasmid encoding or devoid of mcr-1.

RESULTS

In gnotobiotic and conventional mice, mcr-1 significantly enhanced intestinal anchoring of E. coli but impaired their lethal effect. This improvement of intestinal fitness was associated with a downregulation of intestinal inflammatory markers and the preservation of intestinal microbiota composition. The mcr-1 gene mediated a cross-resistance to antimicrobial peptides secreted by the microbiota and intestinal epithelial cells (IECs), enhanced E. coli adhesion to IECs, and decreased the proinflammatory activity of both E. coli and its lipopolysaccharides.

CONCLUSION

Overall, mcr-1 changed multiple facets of bacterial behaviour and appeared as a factor enhancing commensal lifestyle and persistence in the gut even in the absence of antibiotics. Video Abstract.

New understanding of multidrug efflux and permeation in antibiotic resistance, persistence, and heteroresistance

Antibiotics effective against Gram-negative ESKAPE pathogens are a critical area of unmet need. Infections caused by these pathogens are not only difficult to treat but finding new therapies to overcome Gram-negative resistance is also a challenge. There are not enough antibiotics in development that target the most dangerous pathogens and there are not enough novel drugs in the pipeline. The major obstacle in the antibiotic discovery pipeline is the lack of understanding of how to breach antibiotic permeability barriers of Gram-negative pathogens. These barriers are created by active efflux pumps acting across both the inner and the outer membranes. Overproduction of efflux pumps alone or together with either modification of the outer membrane or antibiotic-inactivating enzymes and target mutations contribute to clinical levels of antibiotics resistance. Recent efforts have generated significant advances in the rationalization of compound efflux and permeation across the cell envelopes of Gram-negative pathogens. Combined with earlier studies and novel mathematical models, these efforts have led to a multilevel understanding of how antibiotics permeate these barriers and how multidrug efflux and permeation contribute to the development of antibiotic resistance and heteroresistance. Here, we discuss the new developments in this area.

Membrane-Specific Binding of 4 nm Lipid Nanoparticles Mediated by an Entropy-Driven Interaction Mechanism

Lipid nanoparticles (LNPs) are a leading biomimetic drug delivery platform due to their distinctive advantages and highly tunable formulations. A mechanistic understanding of the interaction between LNPs and cell membranes is essential for developing the cell-targeted carriers for precision medicine. Here the interactions between sub 10 nm cationic LNPs (cLNPs; e.g., 4 nm in size) and varying model cell membranes are systematically investigated using molecular dynamics simulations. We find that the membrane-binding behavior of cLNPs is governed by a two-step mechanism that is initiated by direct contact followed by a more crucial lipid exchange (dissociation of cLNP's coating lipids and subsequent flip and intercalation into the membrane). Importantly, our simulations demonstrate that the membrane binding of cLNPs is an entropy-driven process, which thus enables cLNPs to differentiate between membranes having different lipid compositions (e.g., the outer and inner membranes of bacteria vs the red blood cell membranes). Accordingly, the possible strategies to drive the membrane-targeting behaviors of cLNPs, which mainly depend on the entropy change in the complicated entropy-enthalpy competition of the cLNP-membrane interaction process, are investigated. Our work unveils the molecular mechanism underlying the membrane selectivity of cLNPs and provides useful hints to develop cLNPs as membrane-targeting agents for precision medicine.

Surface Assessment via Grid Evaluation (SuAVE) for Every Surface Curvature and Cavity Shape

The surface assessment via grid evaluation (SuAVE) software was developed to account for the effect of curvature in the calculations of structural properties of chemical interfaces regardless of the chemical composition, asymmetry, and level of atom coarseness. It employs differential geometry techniques, enabling the representation of chemical surfaces as fully differentiable. In this article, we present novel developments of SuAVE to treat closed surfaces and complex cavity shapes. These developments expand the repertoire of curvature-dependent analyses already available in the previous version of SuAVE (e.g., area per lipid, density profiles, membrane thickness, deuterium-order parameters, volume per lipid, and surface curvature angle) to include new functionalities applicable to soft matter (e.g., sphericity, average radius, principal moment of inertia, and roundness) and crystalline porous materials (e.g., pore diameter, internal void volume, total area, and the total void volume of the unit cell structure). SuAVE can accurately handle chemical systems with high and low atom density as demonstrated for two distinct chemical systems: the lipid A vesicle and a set of selected metal–organic frameworks. The SuAVE software v2.0 is fully parallel and benefits from a compiler that supports OpenMP. SuAVE is freely available from https://github.com/SuAVE-Software/source and https://www.biomatsite.net/.

Model architectures for bacterial membranes

The complex composition of bacterial membranes has a significant impact on the understanding of pathogen function and their development towards antibiotic resistance. In addition to the inherent complexity and biosafety risks of studying biological pathogen membranes, the continual rise of antibiotic resistance and its significant economical and clinical consequences has motivated the development of numerous in vitro model membrane systems with tuneable compositions, geometries, and sizes. Approaches discussed in this review include liposomes, solid-supported bilayers, and computational simulations which have been used to explore various processes including drug-membrane interactions, lipid-protein interactions, host-pathogen interactions, and structure-induced bacterial pathogenesis. The advantages, limitations, and applicable analytical tools of all architectures are summarised with a perspective for future research efforts in architectural improvement and elucidation of resistance development strategies and membrane-targeting antibiotic mechanisms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-021-00913-7.

The Whole Is Bigger than the Sum of Its Parts: Drug Transport in the Context of Two Membranes with Active Efflux

Cell envelope plays a dual role in the life of bacteria by simultaneously protecting it from a hostile environment and facilitating access to beneficial molecules. At the heart of this ability lie the restrictive properties of the cellular membrane augmented by efflux transporters, which preclude intracellular penetration of most molecules except with the help of specialized uptake mediators. Recently, kinetic properties of the cell envelope came into focus driven on one hand by the urgent need in new antibiotics and, on the other hand, by experimental and theoretical advances in studies of transmembrane transport. A notable result from these studies is the development of a kinetic formalism that integrates the Michaelis-Menten behavior of individual transporters with transmembrane diffusion and offers a quantitative basis for the analysis of intracellular penetration of bioactive compounds. This review surveys key experimental and computational approaches to the investigation of transport by individual translocators and in whole cells, summarizes key findings from these studies and outlines implications for antibiotic discovery. Special emphasis is placed on Gram-negative bacteria, whose envelope contains two separate membranes. This feature sets these organisms apart from Gram-positive bacteria and eukaryotic cells by providing them with full benefits of the synergy between slow transmembrane diffusion and active efflux.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Atomistic Scale Effects of Lipopolysaccharide Modifications on Bacterial Outer Membrane Defenses

Lipopolysaccharides (LPS) are a main constituent of the outer membrane of Gram-negative bacteria. Salmonella enterica, like many other bacterial species, are able to chemically modify the structure of their LPS molecules through the PhoPQ pathway as a defense mechanism against the host immune response. These modifications make the outer membrane more resistant to antimicrobial peptides (AMPs), large lipophilic drugs, and cation depletion, and are crucial for survival within a host organism. It is believed that these LPS modifications prevent the penetration of large molecules and AMPs through a strengthening of lateral interactions between neighboring LPS molecules. Here, we performed a series of long-timescale molecular dynamics simulations to study how each of three key S. enterica lipid A modifications affect bilayer properties, with a focus on membrane structural characteristics, lateral interactions, and the divalent cation bridging network. Our results discern the unique impact each modification has on strengthening the bacterial outer membrane through effects such as increased hydrogen bonding and tighter lipid packing. Additionally, one of the modifications studied shifts Ca²⁺ from the lipid A region, replacing it as a major cross-linking agent between adjacent lipids and potentially making bacteria less susceptible to AMPs that competitively displace cations from the membrane surface. These results further improve our understanding of outer membrane chemical properties and help elucidate how outer membrane modification systems, such as PhoPQ in S. enterica, are able to alter bacterial virulence.

Molecular dynamics modeling of Pseudomonas aeruginosa outer membranes

Pseudomonas aeruginosa is a common Gram-negative bacterium and opportunistic human pathogen. The distinctive structure of its outer membrane (OM) and outer membrane vesicles (OMVs) plays a fundamental role in bacterial virulence, colonization ability, and antibiotic resistance. To provide critical insights into OM and OMV functionality, we conducted an all-atom molecular dynamics study of asymmetric membranes that are biologically relevant to P. aeruginosa. We hybridized a GLYCAM06-based lipopolysaccharides force field with the Stockholm lipids force field (Slipids) to model bilayer membranes with Lipid A molecules in one leaflet and physiologically relevant phospholipid molecules in the other, including 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG), and 1,2-dioleoyl-sn-glycero-3-phosphoglycerol (DOPG). In particular, a membrane with phospholipid composition representing the P. aeruginosa OM was constructed and modeled by mixing the physiologically dominant components. The detailed structure of membranes was characterized by area per lipid, transmembrane mass and charge densities, radial distribution function (RDF), deuterium order parameter (SCD) of acyl chains, and inclination angles of phosphates and disaccharide in Lipid A. The membrane fluidity in equilibrium and the hydration of functional groups were probed and characterized quantitatively. The consistent properties of the Lipid A leaflets in different membranes demonstrate its compatibility with various phospholipids present in the P. aeruginosa OM. The more ordered acyl chains of Lipid A compared to the cytoplasmic cell membrane contribute to the low permeability of bacterial outer membrane. The findings of this computational investigation of P. aeruginosa OM will further the understanding of microbial pathogenesis and enable future study of OMV biogenesis.

Prediction of the Closed Conformation and Insights into the Mechanism of the Membrane Enzyme LpxR

Covalent modification of outer membrane lipids of Gram-negative bacteria can impact the ability of the bacterium to develop resistance to antibiotics as well as modulating the immune response of the host. The enzyme LpxR from Salmonella typhimurium is known to deacylate lipopolysaccharide molecules of the outer membrane; however, the mechanism of action is unknown. Here, we employ molecular dynamics and Monte Carlo simulations to study the conformational dynamics and substrate binding of LpxR in representative outer membrane models as well as detergent micelles. We examine the roles of conserved residues and provide an understanding of how LpxR binds its substrate. Our simulations predict that the catalytic H122 must be Nε-protonated for a single water molecule to occupy the space between it and the scissile bond, with a free binding energy of -8.5 kcal mol^-1. Furthermore, simulations of the protein within a micelle enable us to predict the structure of the putative "closed" protein. Our results highlight the need for including dynamics, a representative environment, and the consideration of multiple tautomeric and rotameric states of key residues in mechanistic studies; static structures alone do not tell the full story.

Molecular Mechanism of Uptake of Cationic Photoantimicrobial Phthalocyanine across Bacterial Membranes Revealed by Molecular Dynamics Simulations

Phthalocyanines are aromatic macrocyclic compounds, which are structurally related to porphyrins. In clinical practice, phthalocyanines are used in fluorescence imaging and photodynamic therapy of cancer and noncancer lesions. Certain forms of the substituted polycationic metallophthalocyanines have been previously shown to be active in photodynamic inactivation of both Gram-negative and Gram-positive bacteria; one of them is zinc octakis(cholinyl)phthalocyanine (ZnPcChol⁸⁺). However, the molecular details of how these compounds translocate across bacterial membranes still remain unclear. In the present work, we have developed a coarse-grained (CG) molecular model of ZnPcChol⁸⁺ within the framework of the popular MARTINI CG force field. The obtained model was used to probe the solvation behavior of phthalocyanine molecules, which agreed with experimental results. Subsequently, it was used to investigate the molecular details of interactions between phthalocyanines and membranes of various compositions. The results demonstrate that ZnPcChol⁸⁺ has high affinity to both the inner and the outer model membranes of Gram-negative bacteria, although this species does not show noticeable affinity to the 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphatidylcholine membrane. Furthermore, we found out that the process of ZnPcChol⁸⁺ penetration toward the center of the outer bacterial membrane is energetically favorable and leads to its overall disturbance and formation of the aqueous pore. Such intramembrane localization of ZnPcChol⁸⁺ suggests their twofold cytotoxic effect on bacterial cells: (1) via induction of lipid peroxidation by enhanced production of reactive oxygen species (i.e., photodynamic toxicity); (2) via rendering the bacterial membrane more permeable for additional Pc molecules as well as other compounds. We also found that the kinetics of penetration depends on the presence of phospholipid defects in the lipopolysaccharide leaflet of the outer membrane and the type of counterions, which stabilize it. Thus, the results of our simulations provide a detailed molecular view of ZnPcChol⁸⁺ "self-promoted uptake", the pathway previously proposed for some small molecules crossing the outer bacterial membrane.

Trans-envelope multidrug efflux pumps of Gram-negative bacteria and their synergism with the outer membrane barrier

Antibiotic resistance is a serious threat to public health. Significant efforts are currently directed toward containment of the spread of resistance, finding new therapeutic options concerning resistant human and animal pathogens, and addressing the gaps in the fundamental understanding of mechanisms of resistance. Experimental data and kinetic modeling revealed a major factor in resistance, the synergy between active efflux and the low permeability barrier of the outer membrane, which dramatically reduces the intracellular accumulation of many antibiotics. The structural and mechanistic particularities of trans-envelope efflux pumps amplify the effectiveness of cell envelopes as permeability barriers. An important feature of this synergism is that efflux pumps and the outer membrane barriers are mechanistically independent and select antibiotics based on different physicochemical properties. The synergism amplifies even weak polyspecificity of multidrug efflux pumps and creates a major hurdle in the discovery and development of new therapeutics against Gram-negative pathogens.

Synergy between Active Efflux and Outer Membrane Diffusion Defines Rules of Antibiotic Permeation into Gram-Negative Bacteria

Gram-negative bacteria are notoriously resistant to antibiotics, but the extent of the resistance varies broadly between species. We report that in significant human pathogens Acinetobacter baumannii, Pseudomonas aeruginosa, and Burkholderia spp., the differences in antibiotic resistance are largely defined by their penetration into the cell. For all tested antibiotics, the intracellular penetration was determined by a synergistic relationship between active efflux and the permeability barrier. We found that the outer membrane (OM) and efflux pumps select compounds on the basis of distinct properties and together universally protect bacteria from structurally diverse antibiotics. On the basis of their interactions with the permeability barriers, antibiotics can be divided into four clusters that occupy defined physicochemical spaces. Our results suggest that rules of intracellular penetration are intrinsic to these clusters. The identified specificities in the permeability barriers should help in the designing of successful therapeutic strategies against antibiotic-resistant pathogens.IMPORTANCE Multidrug-resistant strains of Gram-negative pathogens rapidly spread in clinics. Significant efforts worldwide are currently directed to finding the rules of permeation of antibiotics across two membrane envelopes of these bacteria. This study created the tools for analysis of and identified the major differences in antibacterial activities that distinguish the permeability barriers of P. aeruginosa, A. baumannii, Burkholderia thailandensis, and B. cepacia We conclude that synergy between active efflux and the outer membrane barrier universally protects Gram-negative bacteria from antibiotics. We also found that the diversity of antibiotics affected by active efflux and outer membrane barriers is broader than previously thought and that antibiotics cluster according to specific biological determinants such as the requirement of specific porins in the OM, targeting of the OM, or specific recognition by efflux pumps. No universal rules of antibiotic permeation into Gram-negative bacteria apparently exist. Our results suggest that antibiotic clusters are defined by specific rules of permeation and that further studies could lead to their discovery.