Computer simulation of the rough lipopolysaccharide membrane of Pseudomonas aeruginosa

The Effect of Lipopolysaccharides on the Electrostatic Properties of Gram-Negative General Porins from Enterobacteriaceae

We investigated, by using all-atom molecular dynamics simulations, the effect of the outer membrane of Gram-negative bacteria, composed in the outer leaflet by polar/charged lipopolysaccharides (LPS), on the electrostatic properties of general porins from the Enterobacteriaceae family. General porins constitute the main path for the facilitated diffusion of polar antibiotics through the outer membrane. As model system we selected OmpK36 from Klebsiella pneumoniae, the ortholog of OmpC from Escherichia coli. This species presents high variability of amino acid composition of porins, with the effect to increase its resistance to the penetration of antibiotics. The various properties we analyzed seem to indicate that LPS acts as an independent layer without affecting the internal electrostatic properties of OmpK36. The only apparent effect on the microsecond time scale we sampled is the appearance of calcium ions, when present at moderate concentration in solution, inside the pore. However, we noticed increased fluctuations of the polarization density and only minor changes on its average value.

Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

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.

What have molecular simulations contributed to understanding of Gram-negative bacterial cell envelopes?

Bacterial cell envelopes are compositionally complex and crowded and while highly dynamic in some areas, their molecular motion is very limited, to the point of being almost static in others. Therefore, it is no real surprise that studying them at high resolution across a range of temporal and spatial scales requires a number of different techniques. Details at atomistic to molecular scales for up to tens of microseconds are now within range for molecular dynamics simulations. Here we review how such simulations have contributed to our current understanding of the cell envelopes of Gram-negative bacteria.

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.

Lipopolysaccharides at Solid and Liquid Interfaces: Models for Biophysical Studies of the Gram-negative Bacterial Outer Membrane

Lipopolysaccharides (LPSs) are a constitutive element of the cell envelope of Gram-negative bacteria, representing the main lipid in the external leaflet of their outer membrane (OM) lipid bilayer. These unique surface-exposed glycolipids play a central role in the interactions of Gram-negative organisms with their surrounding environment and represent a key element for protection against antimicrobials and the development of antibiotic resistance. The biophysical investigation of a wide range of different types of in vitro model membranes containing reconstituted LPS has revealed functional and structural properties of these peculiar membrane lipids, providing molecular-level details of their interaction with antimicrobial compounds. LPS assemblies reconstituted at interfaces represent a versatile tool to study the properties of the Gram-negative OM by exploiting several surface-sensitive techniques, in particular X-ray and neutron scattering, which can probe the structure of thin films with sub-nanometer resolution. This review provides an overview of different approaches employed to investigate structural and biophysical properties of LPS, focusing on studies on Langmuir monolayers of LPS at the air/liquid interface and a range of supported LPS-containing model membranes reconstituted at solid/liquid interfaces.

In silico investigation and surmounting of Lipopolysaccharide barrier in Gram-Negative Bacteria: How far has molecular dynamics Come?

Lipopolysaccharide (LPS), a main component of the outer membrane of Gram-negative bacteria, has crucial implications on both antibiotic resistance and the overstimulation of the host innate immune system. Fighting against these global concerns calls for the molecular understanding of the barrier function and immunostimulatory ability of LPS. Molecular dynamics (MD) simulations have become an invaluable tool for uncovering important findings in LPS research. While the reach of MD simulations for investigating the immunostimulatory ability of LPS has been already outlined, little attention has been paid to the role of MD simulations for exploring its barrier function and synthesis. Herein, we give an overview about the impact of MD simulations on gaining insight into the shield role and synthesis pathway of LPS, which have attracted considerable attention to discover molecules able to surmount antibiotic resistance, either circumventing LPS defenses or disrupting its synthesis. We specifically focus on the enhanced sampling and free energy calculation methods that have been combined with MD simulations to address such research. We also highlight the use of special-purpose MD supercomputers, the importance of appropriate LPS and ions parameterization to obtain reliable results, and the complementary views that MD and wet-lab experiments provide. Thereby, this work, which covers the last five years of research, apart from outlining the phenomena and strategies that are being explored, evidences the valuable insights that are gained by MD, which may be useful to advance antibiotic design, and what the prospects of this in silico method could be in LPS research.

Structural changes in the model of the outer cell membrane of Gram-negative bacteria interacting with melittin: an in situ spectroelectrochemical study

The cell membrane of Gram-negative bacteria interacting with an antimicrobial peptide presents a complex supramolecular assembly. Fabrication of models of bacterial cell membranes remains a large experimental challenge. Langmuir-Blodgett and Langmuir-Schaefer (LS-LB) transfer makes possible the deposition of multicomponent asymmetric lipid bilayers onto a gold surface. Two lipids: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and di[3-deoxy-D-manno-octulosonyl]-lipid A (KLA) were used to deposit a model of the outer membrane of Gram-negative bacteria on the Au(111) substrate. The use of gold as the solid substrate enables control of the membrane potential. Molecular scale changes in the model membrane exposed to physiological electric fields and interacting with melittin antimicrobial peptide are discussed in this paper. The interaction of the outer membrane with melittin leads to an increase in the membrane capacitance and permeability to ions and water. The stability of the outer membrane with bound melittin decreases at positive membrane potentials. In situ polarization modulation infrared reflection absorption spectroscopy is used to investigate membrane potential-dependent changes in the structure of the outer membrane interacting with melittin. The hydration of the ester carbonyl groups is not affected by the interaction with melittin. However, the orientation and hydrogen bond network with the carboxylate groups in KLA changes drastically after POPE-KLA bilayer interacts with melittin. We propose that the positively charged groups in the amino acids present at the C-terminus of the peptide interact directly with the polar head group of KLA. Simultaneously, the packing order in hydrocarbon chains in the membrane with bound melittin increases. A hydrophobic match between the chains in the lipids and the peptide, which spans the membrane, seems to be responsible for the ordering of the hydrocarbon chains region of the bilayer. The N-terminus enters into the hydrophobic region of the membrane and forms a channel to the hydrophilic head groups in POPE.

Osmotic stress and vesiculation as key mechanisms controlling bacterial sensitivity and resistance to TiO2 nanoparticles

Toxicity mechanisms of metal oxide nanoparticles towards bacteria and underlying roles of membrane composition are still debated. Herein, the response of lipopolysaccharide-truncated Escherichia coli K12 mutants to TiO2 nanoparticles (TiO2NPs, exposure in dark) is addressed at the molecular, single cell, and population levels by transcriptomics, fluorescence assays, cell nanomechanics and electrohydrodynamics. We show that outer core-free lipopolysaccharides featuring intact inner core increase cell sensitivity to TiO2NPs. TiO2NPs operate as membrane strippers, which induce osmotic stress, inactivate cell osmoregulation and initiate lipid peroxidation, which ultimately leads to genesis of membrane vesicles. In itself, truncation of lipopolysaccharide inner core triggers membrane permeabilization/depolarization, lipid peroxidation and hypervesiculation. In turn, it favors the regulation of TiO2NP-mediated changes in cell Turgor stress and leads to efficient vesicle-facilitated release of damaged membrane components. Remarkably, vesicles further act as electrostatic baits for TiO2NPs, thereby mitigating TiO2NPs toxicity. Altogether, we highlight antagonistic lipopolysaccharide-dependent bacterial responses to nanoparticles and we show that the destabilized membrane can generate unexpected resistance phenotype.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

Machine Learning Algorithm Identifies an Antibiotic Vocabulary for Permeating Gram-Negative Bacteria

Drug discovery faces a crisis. The industry has used up the "obvious" space in which to find novel drugs for biomedical applications, and productivity is declining. One strategy to combat this is rational approaches to expand the search space without relying on chemical intuition, to avoid rediscovery of similar spaces. In this work, we present proof of concept of an approach to rationally identify a "chemical vocabulary" related to a specific drug activity of interest without employing known rules. We focus on the pressing concern of multidrug resistance in Pseudomonas aeruginosa by searching for submolecules that promote compound entry into this bacterium. By synergizing theory, computation, and experiment, we validate our approach, explain the molecular mechanism behind identified fragments promoting compound entry, and select candidate compounds from an external library that display good permeation ability.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Molecular Simulations of Gram-Negative Bacterial Membranes Come of Age

Gram-negative bacteria are protected by a multicompartmental molecular architecture known as the cell envelope that contains two membranes and a thin cell wall. As the cell envelope controls influx and efflux of molecular species, in recent years both experimental and computational studies of such architectures have seen a resurgence due to the implications for antibiotic development. In this article we review recent progress in molecular simulations of bacterial membranes. We show that enormous progress has been made in terms of the lipidic and protein compositions of bacterial systems. The simulations have moved away from the traditional setup of one protein surrounded by a large patch of the same lipid type toward a more bio-logically representative viewpoint. Simulations with multiple cell envelope components are also emerging. We review some of the key method developments that have facilitated recent progress, discuss some current limitations, and offer a perspective on future directions.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

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.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Structure and Conformation of Wild-Type Bacterial Lipopolysaccharide Layers at Air-Water Interfaces

The outer membrane of Gram-negative bacteria is of great scientific interest because it mediates the action of antimicrobial agents. The membrane surface is composed of lipopolysaccharide (LPS) molecules with negatively charged oligosaccharide headgroups. To a certain fraction, LPSs additionally display linear polysaccharides termed O-side chains (OSCs). Structural studies on bacterial outer surfaces models, based on LPS monolayers at air-water interfaces, have so far dealt only with rough mutant LPSs lacking these OSCs. Here, we characterize monolayers of wild-type LPS from Escherichia coli O55:B5 featuring strain-specific OSCs in the presence of defined concentrations of monovalent and divalent ions. Pressure-area isotherms yield insight into in-plane molecular interactions and monolayer elastic moduli. Structural investigations by x-ray and neutron reflectometry reveal the saccharide conformation and allow quantifying the area per molecule and the fraction of LPS molecules carrying OSCs. The OSC conformation is satisfactorily described by the self-consistent field theory for end-grafted polymer brushes. The monolayers exhibit a significant structural response to divalent cations, which goes beyond generic electrostatic screening.

Atomistic and Coarse Grain Simulations of the Cell Envelope of Gram-Negative Bacteria: What Have We Learned?

Bacterial membranes, and those of Gram-negative bacteria in particular, are some of the most biochemically diverse membranes known. They incorporate a wide range of lipid types and proteins of varying sizes, architectures, and functions. While simpler biological membranes have been the focus of myriad simulation studies over the years that have yielded invaluable details to complement, and often to direct, ongoing experimental studies, simulations of complex bacterial membranes have been slower to emerge. However, the past few years have seen tremendous activity in this area, leading to advances such as the development of atomistic and coarse-grain models of the lipopolysaccharide (LPS) component of the outer membrane that are compatible with widely used simulation codes. In this Account, we review our contributions to the field of molecular simulations of the bacterial cell envelope, including the development of models of both membranes and the cell wall of Gram-negative bacteria, with a predominant focus on E. coli. At the atomistic level, simulations of chemically accurate models of both membranes have revealed the tightly cross-linked nature of the LPS headgroups and have shown that penetration of solutes through these regions is not as straightforward as the route through phospholipids. The energetic differences between the two routes have been calculated. Simulations of native outer membrane proteins in LPS-containing membranes have shown that the conformational dynamics of the proteins is not only slower in LPS but also different compared to in simpler models of phospholipid bilayers. These chemically more complex and consequently biologically more relevant models are leading to details of conformational dynamics that were previously inaccessible from simulations. Coarse-grain models have enabled simulations of multiprotein systems on time scales of microseconds, leading to insights not only into the rates of protein and lipid diffusion but also into the trends in their respective directions of flow. We find that the motions of LPS molecules are highly correlated with each other but also with outer membrane proteins embedded within the membrane. We have shown that the two leaflets of the outer membrane exhibit communication, whereby regions of low disorder in one leaflet correspond to regions of high disorder in the other. The cell wall remains a comparatively neglected component, although models of the E. coli peptidoglycan are now emerging, particularly at the atomistic level. Our simulations of Braun's lipoprotein have shown that bending and tilting of this protein afford a degree of variability in the gap between the cell wall and the OM. The noncovalent interactions with the cell wall of proteins such as OmpA can further influence the width of this gap by extension or contraction of their linker domains. Overall we have shown that the dynamics of proteins, lipids, and other molecular species within the outer membrane cannot be approximated using simpler phospholipid bilayers, if one is addressing questions regarding the in vivo behavior of Gram-negative bacteria. These membranes have their own unique chemical characteristics that cannot be decoupled from their biological functions.

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.

Molecular Simulations of Gram-Negative Bacterial Membranes: A Vignette of Some Recent Successes

In the following review we use recent examples from the literature to discuss progress in the area of atomistic and coarse-grained molecular dynamics simulations of selected bacterial membranes and proteins, with a particular focus on Gram-negative bacteria. As structural biology continues to provide increasingly high-resolution data on the proteins that reside within these membranes, simulations have an important role to play in linking these data with the dynamical behavior and function of these proteins. In particular, in the last few years there has been significant progress in addressing the issue of biochemical complexity of bacterial membranes such that the heterogeneity of the lipid and protein components of these membranes are now being incorporated into molecular-level models. Thus, in future we can look forward to complementary data from structural biology and molecular simulations combining to provide key details of structure-dynamics-function relationships in bacterial membranes.

Computer simulations suggest direct and stable tip to tip interaction between the outer membrane channel TolC and the isolated docking domain of the multidrug RND efflux transporter AcrB

One way by which bacteria achieve antibiotics resistance is preventing drug access to its target molecule for example through an overproduction of multi-drug efflux pumps of the resistance nodulation division (RND) protein super family of which AcrAB-TolC in Escherichia coli is a prominent example. Although representing one of the best studied efflux systems, the question of how AcrB and TolC interact is still unclear as the available experimental data suggest that either both proteins interact in a tip to tip manner or do not interact at all but are instead connected by a hexamer of AcrA molecules. Addressing the question of TolC-AcrB interaction, we performed a series of 100 ns - 1 µs-molecular dynamics simulations of membrane-embedded TolC in presence of the isolated AcrB docking domain (AcrB(DD)). In 5/6 simulations we observe direct TolC-AcrB(DD) interaction that is only stable on the simulated time scale when both proteins engage in a tip to tip manner. At the same time we find TolC opening and closing freely on extracellular side while remaining closed at the inner periplasmic bottleneck region, suggesting that either the simulated time is too short or additional components are required to unlock TolC.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

The membranes of Gram-negative bacteria: progress in molecular modelling and simulation

Molecular modelling and simulations have been employed to study the membranes of Gram-negative bacteria for over 20 years. Proteins native to these membranes, as well as antimicrobial peptides and drug molecules have been studied using molecular dynamics simulations in simple models of membranes, usually only comprising one lipid species. Thus, traditionally, the simulations have reflected the majority of in vitro membrane experimental setups, enabling observations from the latter to be rationalized at the molecular level. In the last few years, the sophistication and complexity of membrane models have improved considerably, such that the heterogeneity of the lipid and protein composition of the membranes can now be considered both at the atomistic and coarse-grain levels of granularity. Importantly this means relevant biology is now being retained in the models, thereby linking the in silico and in vivo scenarios. We discuss recent progress in simulations of proteins in simple lipid bilayers, more complex membrane models and finally describe some efforts to overcome timescale limitations of atomistic molecular dynamics simulations of bacterial membranes.

Anti-infectives in Drug Delivery-Overcoming the Gram-Negative Bacterial Cell Envelope

Infectious diseases are becoming a major menace to the state of health worldwide, with difficulties in effective treatment especially of nosocomial infections caused by Gram-negative bacteria being increasingly reported. Inadequate permeation of anti-infectives into or across the Gram-negative bacterial cell envelope, due to its intrinsic barrier function as well as barrier enhancement mediated by resistance mechanisms, can be identified as one of the major reasons for insufficient therapeutic effects. Several in vitro, in silico, and in cellulo models are currently employed to increase the knowledge of anti-infective transport processes into or across the bacterial cell envelope; however, all such models exhibit drawbacks or have limitations with respect to the information they are able to provide. Thus, new approaches which allow for more comprehensive characterization of anti-infective permeation processes (and as such, would be usable as screening methods in early drug discovery and development) are desperately needed. Furthermore, delivery methods or technologies capable of enhancing anti-infective permeation into or across the bacterial cell envelope are required. In this respect, particle-based carrier systems have already been shown to provide the opportunity to overcome compound-related difficulties and allow for targeted delivery. In addition, formulations combining efflux pump inhibitors or antimicrobial peptides with anti-infectives show promise in the restoration of antibiotic activity in resistant bacterial strains. Despite considerable progress in this field however, the design of carriers to specifically enhance transport across the bacterial envelope or to target difficult-to-treat (e.g., intracellular) infections remains an urgently needed area of improvement. What follows is a summary and evaluation of the state of the art of both bacterial permeation models and advanced anti-infective formulation strategies, together with an outlook for future directions in these fields.

Energetics of Endotoxin Recognition in the Toll-Like Receptor 4 Innate Immune Response

Bacterial outer membrane lipopolysaccharide (LPS) potently stimulates the mammalian innate immune system, and can lead to sepsis, the primary cause of death from infections. LPS is sensed by Toll-like receptor 4 (TLR4) in complex with its lipid-binding coreceptor MD-2, but subtle structural variations in LPS can profoundly modulate the response. To better understand the mechanism of LPS-induced stimulation and bacterial evasion, we have calculated the binding affinity to MD-2 of agonistic and antagonistic LPS variants including lipid A, lipid IVa, and synthetic antagonist Eritoran, and provide evidence that the coreceptor is a molecular switch that undergoes ligand-induced conformational changes to appropriately activate or inhibit the receptor complex. The plasticity of the coreceptor binding cavity is shown to be essential for distinguishing between ligands, whilst similar calculations for a model bacterial LPS bilayer reveal the "membrane-like" nature of the protein cavity. The ability to predict the activity of LPS variants should facilitate the rational design of TLR4 therapeutics.

Membrane Protein Structure, Function, and Dynamics: a Perspective from Experiments and Theory

Membrane proteins mediate processes that are fundamental for the flourishing of biological cells. Membrane-embedded transporters move ions and larger solutes across membranes; receptors mediate communication between the cell and its environment and membrane-embedded enzymes catalyze chemical reactions. Understanding these mechanisms of action requires knowledge of how the proteins couple to their fluid, hydrated lipid membrane environment. We present here current studies in computational and experimental membrane protein biophysics, and show how they address outstanding challenges in understanding the complex environmental effects on the structure, function, and dynamics of membrane proteins.

Interaction of the antimicrobial peptide polymyxin B1 with both membranes of E. coli: a molecular dynamics study

Antimicrobial peptides are small, cationic proteins that can induce lysis of bacterial cells through interaction with their membranes. Different mechanisms for cell lysis have been proposed, but these models tend to neglect the role of the chemical composition of the membrane, which differs between bacterial species and can be heterogeneous even within a single cell. Moreover, the cell envelope of Gram-negative bacteria such as E. coli contains two membranes with differing compositions. To this end, we report the first molecular dynamics simulation study of the interaction of the antimicrobial peptide, polymyxin B1 with complex models of both the inner and outer membranes of E. coli. The results of >16 microseconds of simulation predict that polymyxin B1 is likely to interact with the membranes via distinct mechanisms. The lipopeptides aggregate in the lipopolysaccharide headgroup region of the outer membrane with limited tendency for insertion within the lipid A tails. In contrast, the lipopeptides readily insert into the inner membrane core, and the concomitant increased hydration may be responsible for bilayer destabilization and antimicrobial function. Given the urgent need to develop novel, potent antibiotics, the results presented here reveal key mechanistic details that may be exploited for future rational drug development.

Insights into the species-specific TLR4 signaling mechanism in response to Rhodobacter sphaeroides lipid A detection

TLR4 in complex with MD2 senses the presence of lipid A (LA) and initiates a signaling cascade that curb the infection. This complex is evolutionarily conserved and can initiate the immune system in response to a variety of LAs. In this study, molecular dynamics simulation (25 ns) was performed to elucidate the differential behavior of TLR4/MD2 complex in response to Rhodobacter sphaeroides lipid A (RsLA). Penta-acyl chain-containing RsLA is at the verge of agonist (6 acyl-chains) and antagonist (4 acyl-chains) structure, and activates the TLR4 pathway in horses and hamsters, while inhibiting in humans and murine. In the time-evolved coordinates, the promising factors that dictated the differential response included the local and global mobility pattern of complexes, solvent-accessible surface area of ligand, and surface charge distributions of TLR4 and MD2. We showed that the GlcN1-GlcN2 backbone acquires agonist (3FXI)-like configurations in horses and hamsters, while acquiring antagonist (2E59)-like configurations in humans and murine systems. Moreover, analysis of F126 behavior in the MD2 F126 loop (amino acids 123-129) and loop EF (81-89) suggested that certain sequence variations also contribute to species-specific response. This study underlines the TLR4 signaling mechanism and provides new therapeutic opportunities.

Lipopolysaccharide membrane building and simulation

While membrane simulations are widely employed to study the structure and dynamics of various lipid bilayers and membrane proteins in the bilayers, simulations of lipopolysaccharides (LPS) in membrane environments have been limited due to their structural complexity, difficulties in building LPS-membrane systems, and lack of the appropriate molecular force fields. In this work, as a first step to extend CHARMM-GUI Membrane Builder to incorporate LPS molecules and to explore their structures and dynamics in membrane environments using molecular dynamics simulations, we describe step-by-step procedures to build LPS bilayer systems using CHARMM and the recently developed CHARMM carbohydrate and lipid force fields. Such procedures are illustrated by building various bilayers of Escherichia coli R1.O6 LPS and the presentation of preliminary simulation results in terms of per-LPS area and density distributions of various components along the membrane normal.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

Modeling the electrostatic potential of asymmetric lipopolysaccharide membranes: the MEMPOT algorithm implemented in DelPhi

Four chemotypes of the rough lipopolysaccharides (LPS) membrane from Pseudomonas aeruginosa were investigated by a combined approach of explicit water molecular dynamics (MD) simulations and Poisson-Boltzmann continuum electrostatics with the goal to deliver the distribution of the electrostatic potential across the membrane. For the purpose of this investigation, a new tool for modeling the electrostatic potential profile along the axis normal to the membrane, MEMbrane POTential (MEMPOT), was developed and implemented in DelPhi. Applying MEMPOT on the snapshots obtained by MD simulations, two observations were made: (a) the average electrostatic potential has a complex profile but is mostly positive inside the membrane due to the presence of Ca(2+) ions, which overcompensate for the negative potential created by lipid phosphate groups; and (b) correct modeling of the electrostatic potential profile across the membrane requires taking into account the water phase, while neglecting it (vacuum calculations) results in dramatic changes including a reversal of the sign of the potential inside the membrane. Furthermore, using DelPhi to assign different dielectric constants for different regions of the LPS membranes, it was investigated whether a single frame structure before MD simulations with appropriate dielectric constants for the lipid tails, inner, and the external leaflet regions, can deliver the same average electrostatic potential distribution as obtained from the MD-generated ensemble of structures. Indeed, this can be attained by using smaller dielectric constant for the tail and inner leaflet regions (mostly hydrophobic) than for the external leaflet region (hydrophilic) and the optimal dielectric constant values are chemotype-specific.

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.

New amphiphilic neamine derivatives active against resistant Pseudomonas aeruginosa and their interactions with lipopolysaccharides

The development of novel antimicrobial agents is urgently required to curb the widespread emergence of multidrug-resistant bacteria like colistin-resistant Pseudomonas aeruginosa. We previously synthesized a series of amphiphilic neamine derivatives active against bacterial membranes, among which 3',6-di-O-[(2"-naphthyl)propyl]neamine (3',6-di2NP), 3',6-di-O-[(2"-naphthyl)butyl]neamine (3',6-di2NB), and 3',6-di-O-nonylneamine (3',6-diNn) showed high levels of activity and low levels of cytotoxicity (L. Zimmermann et al., J. Med. Chem. 56:7691-7705, 2013). We have now further characterized the activity of these derivatives against colistin-resistant P. aeruginosa and studied their mode of action; specifically, we characterized their ability to interact with lipopolysaccharide (LPS) and to alter the bacterial outer membrane (OM). The three amphiphilic neamine derivatives were active against clinical colistin-resistant strains (MICs, about 2 to 8 μg/ml), The most active one (3',6-diNn) was bactericidal at its MIC and inhibited biofilm formation at 2-fold its MIC. They cooperatively bound to LPSs, increasing the outer membrane permeability. Grafting long and linear alkyl chains (nonyl) optimized binding to LPS and outer membrane permeabilization. The effects of amphiphilic neamine derivatives on LPS micelles suggest changes in the cross-bridging of lipopolysaccharides and disordering in the hydrophobic core of the micelles. The molecular shape of the 3',6-dialkyl neamine derivatives induced by the nature of the grafted hydrophobic moieties (naphthylalkyl instead of alkyl) and the flexibility of the hydrophobic moiety are critical for their fluidifying effect and their ability to displace cations bridging LPS. Results from this work could be exploited for the development of new amphiphilic neamine derivatives active against colistin-resistant P. aeruginosa.

Molecular dynamics and NMR spectroscopy studies of E. coli lipopolysaccharide structure and dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From (1)H,(1)H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.

Extracellular polymeric substances from copper-tolerance Sinorhizobium meliloti immobilize Cu²⁺

The copper tolerance gene of wild-type heavy metal-tolerance Sinorhizobium meliloti CCNWSX0020 was mutated by transposon Tn5-a. The mutant was sensitive up to 1.4mM Cu(2+). Production, components, surface morphology, and functional groups of extracellular polymeric substances (EPS) of the wild-type strains were compared with sensitive mutant in immobilization of Cu(2+). EPS produced by S. meliloti CCNWSX0020 restricts uptake of Cu(2+). The cell wall EPS were categorized based on the compactness and fastness: soluble EPS (S-EPS), loosely bound EPS (LB-EPS), and tightly bound EPS (TB-EPS). LB-EPS played a more important role than S-EPS and TB-EPS in Cu(2+) immobilization. Scanning electron microscopy (SEM) analysis LB-EPS had rough surface and many honeycomb pores, making them conducive to copper entry; therefore, they may play a role as a microbial protective barrier. Fourier transform-infrared (FT-IR) analysis further confirm that proteins and carbohydrates were the main extracellular compounds which had functional groups such as carboxyl (COOH), hydroxyl (OH), and amide (NH), primarily involved in metal ion binding.

Physical interactions of fish protamine and antisepsis peptide drugs with bacterial membranes revealed by combination of specular x-ray reflectivity and grazing-incidence x-ray fluorescence

As a defined model of outer membranes of gram negative bacteria, we investigated the interaction of monolayers of lipopolysacchrides from Salmonella enterica rough strains R90 (LPS Ra) with natural and synthetic peptides. The fine structures perpendicular to the membrane plane and the ion distribution near the interface were determined by specular x-ray reflectivity (XRR) and grazing-incidence x-ray fluorescence (GIXF) in the presence and absence of divalent cations. The unique combination of XRR and GIXF allows for the quantitative identification of different modes of interactions in a high spatial resolution, which cannot be assessed by other experimental methods. Natural fish protamine disrupts the stratified membrane structures in the absence of Ca(2+) ions, while staying away from the membrane surface in the presence of Ca(2+) ions. In contrast, synthetic antisepsis peptide Pep 19-2.5 weakly adsorbs to the membrane and stays near the uncharged sugar units even in the absence of Ca(2+). In the presence of Ca(2+), Pep 19-2.5 can reach the negatively charged inner core without destroying the barrier capability against ions.

Lipid A from lipopolysaccharide recognition: structure, dynamics and cooperativity by molecular dynamics simulations

Molecular dynamics simulations of Lipid A and its natural precursor Lipid IVA from E.coli have been carried out free in solution, bound to the myeliod differentiation protein 2 (MD2) and in the complex of MD2 with the toll like receptor 4 (TLR4). In addition, simulations of the ligand free MD2 and MD2-TLR4 complex were performed. A structural and energetic characterization of the bound and unbound states of Lipid A/IVA was generated. As the crystal structures depict, the main driving force for MD2-Lipid A/IVA are the hydrophobic interactions between the aliphatic tails and the MD2 cavity. The charged phosphate groups do strongly interact with positively charged residues, located at the surface of MD2. However, they are not essential for keeping the lipids in the cavity, indicating a more prominent role in binding recognition and ionic interactions with TLR4 at the MD2/TLR4 interface. Interestingly, in the absence of any ligand MD2 rapidly closes, blocking the binding cavity. The presence of TLR4, though changing the dynamics, was not able to impede the aforementioned closing event. We hypothesize that fluctuations of the H1 region are essential for this phenomenon, and it is plausible that an equilibrium between the open and closed states exists, although the lengths of our simulations are not sufficient to encompass the reversible process. The MD2/Lipid A-TLR4 complex simulations show that the presence of the ligand energetically stabilizes the complex relative to the ligand-free structures, indicating cooperativity in the binding process.

Supported lipopolysaccharide bilayers

In this report, the formation of supported lipopolysaccharide bilayers (LPS-SLBs) is studied with extracted native and glycoengineered LPS from Escherichia coli ( E. coli ) and Salmonella enterica sv typhimurium ( S. typhimurium ) to assemble a platform that allows measurement of LPS membrane structure and the detection of membrane tethered saccharide-protein interactions. We present quartz crystal microbalance with dissipation monitoring (QCM-D) and fluorescence recovery after photobleaching (FRAP) characterization of LPS-SLBs with different LPS species, having, for example, different molecular weights, that show successful formation of SLBs through vesicle fusion on SiO(2) surfaces with LPS fractions up to 50 wt %. The thickness of the LPS bilayers were investigated with AFM force-distance measurements which showed only a slight thickness increase compared to pure POPC SLBs. The E. coli LPS were chosen to study the saccharide-protein interaction between the Htype II glycan epitope and the Ralstonia solanacearum lectin (RSL). RSL specifically recognizes fucose sugars, which are present in the used Htype II glycan epitope and absent in the epitopes LPS1 and EY2. We show via fluorescence microscopy that the specific, but weak and multivalent interaction can be detected and discriminated on the LPS-SLB platform.