Progress in Molecular Dynamics Simulations of Gram-Negative Bacterial Cell Envelopes

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Stability of Cytoplasmic Membrane of Escherichia coli Bacteria in Aqueous and Ethanolic Environment

The mechanism of action of any antibacterial agent or disinfectant depends largely on their interaction with the bacterial membrane. Herein, we use the SPICA (surface property fitting coarse graining) force-field and develop a coarse-grained (CG) model for the structure of the cytoplasmic membrane of Escherichia coli (E. coli) and its interaction with water and ethanol. We elucidate the impact of different concentrations of ethanol on the cytoplasmic membrane bilayers and vesicles of E. coli using the CG molecular dynamics (CG MD) simulations. Our modeling approach first focuses on the parametrization of the required force-field for POPG lipid and its interaction with water, ethanol, and POPE lipid. Subsequently, the structural stability of the E. coli bacterial membrane in the presence of high and low concentrations of ethanol is delineated. Both flat bilayers as well as vesicles of E. coli membrane were considered for the CG MD. Our results reveal that, at low ethanol concentrations (<30 mol %), the size of the E. coli vesicles increases with discernible deformations in their shapes. Because of ethanol-induced interdigitation, thinning of the E. coli vesicular membrane is also observed. However, at higher ethanol concentrations (>30 mol %), the integrity of the vesicles is lost because of deteriorating invasion of ethanol molecules into the vesicle bilayer and significant weakening of lipid-lipid interactions. At higher ethanol concentrations (40 and 70 mol %), both the multivesicle and single-vesicle bacterial membranes exhibit a similar rupturing pattern wherein the extraction of lipids from the membrane and formation of aggregates of the component lipids are observed. These aggregates consist of polar head groups of 3-5 POPE/POPG lipids with intertwined nonpolar tails.

Structural Model of a Porphyromonas gingivalis type IX Secretion System Shuttle Complex

Porphyromonas gingivalis is a gram-negative oral anaerobic pathogen and is one of the key causative agents of periodontitis. P. gingivalis utilises a range of virulence factors, including the cysteine protease RgpB, to drive pathogenesis and these are exported and attached to the cell surface via the type IX secretion system (T9SS). All cargo proteins possess a conserved C-terminal signal domain (CTD) which is recognised by the T9SS, and the outer membrane β-barrel protein PorV (PG0027/LptO) can interact with cargo proteins as they are exported to the bacterial surface. Using a combination of solution nuclear magnetic resonance (NMR) spectroscopy, biochemical analyses, machine-learning-based modelling and molecular dynamics (MD) simulations, we present a structural model of a PorV:RgpB-CTD complex from P. gingivalis. This is the first structural insight into CTD recognition by the T9SS and shows how the conserved motifs in the CTD are the primary sites that mediate binding. In PorV, interactions with extracellular surface loops are important for binding the CTD, and together these appear to cradle and lock RgpB-CTD in place. This work provides insight into cargo recognition by PorV but may also have important implications for understanding other aspects of type-IX dependent secretion.

Martini 3 Coarse-Grained Force Field for Carbohydrates

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.

Co-Exposure of Ambient Particulate Matter and Airborne Transmission Pathogens: The Impairment of the Upper Respiratory Systems

Recent evidence has pinpointed the positive relevance between air particulate matter (PM) pollution and epidemic spread. However, there are still significant knowledge gaps in understanding the transmission and infection of pathogens loaded on PMs, for example, the interactions between pathogens and pre-existing atmospheric PM and the health effects of co-exposure on the inhalation systems. Here, we unraveled the interactions between fine particulate matter (FPM) and Pseudomonas aeruginosa (P. aeruginosa) and evaluated the infection and detrimental effects of co-exposure on the upper respiratory systems in both in vitro and in vivo models. We uncovered the higher accessibility and invasive ability of pathogens to epithelial cells after loading on FPMs, compared with the single exposure. Furthermore, we designed a novel laboratory exposure model to simulate a real co-exposure scenario. Intriguingly, the co-exposure induced more serious functional damage and longer inflammatory reactions to the upper respiratory tract, including the nasal cavity and trachea. Collectively, our results provide a new point of view on the transmission and infection of pathogens loaded on FPMs and uncover the in vivo systematic impairments of the inhalation tract under co-exposure through a novel laboratory exposure model. Hence, this study sheds light on further investigations of the detrimental effects of air pollution and epidemic spread.

Self-Uptake Mechanism of Polymyxin-Based Lipopeptide against Gram-Negative Bacterial Membrane: Role of the First Adsorbed Lipopeptide

Research in the continuously increasing threat of polymyxin-resistant multidrug-resistant Pseudomonas aeruginosa, which causes severe infection in immunocompromised patients, has resulted in the development of several polymyxin-derived cyclic lipopeptides containing l-α-γ- diamino butyric acid-like FADDI-019 (F19). In this work, F19's insertion into a minimal model of the asymmetric outer membrane of the bacterium, which contained only penta-acylated lipid A (LipA) and lacked keto-d-octulosonic acid and O-antigens, in the top leaflet and phospholipids in the bottom leaflet, was studied. F19 exhibited all of the hallmarks of the self-uptake mechanism into the asymmetric bilayer. While a single monomer of the lipopeptide did not get partitioned into the inside of the bilayer, it competitively displaced Ca²⁺ from the membrane surface, observed as a decrease in Ca²⁺ coordination number with phosphate groups (1.89 vs 1.718), resulting in membrane destabilization. This resulted in an increment of the average defect size and the probability of interplay between lipid tails and hydrophobic residues of another F19. When more than one monomer was present in the system, the first monomer remained docked on the surface, while other monomers intercalated into the bilayer interior with their hydrophobic moieties "sleeved" by lipid acyl chains. The free energy barrier for partial insertion of the lipopeptide into a bilayer in the presence of surface-docked second F19 was recorded at ∼1.3 kcal/mol using two-dimensional (2D) well-tempered metadynamics, making it a low barrier process at 300 K. This study is an attempt to demonstrate the self-uptake mechanism of F19 during intercalation process into the bilayer interior, which may help in the design of better alternates for polymyxins to work against polymyxin resistance.

Simulation Study of the Effect of Antimicrobial Peptide Associations on the Mechanism of Action with Bacterial and Eukaryotic Membranes

Antimicrobial peptides (AMPs) can be directed to specific membranes based on differences in lipid composition. In this study, we performed atomistic and coarse-grained simulations of different numbers of the designed AMP adepantin-1 with a eukaryotic membrane, cytoplasmic Gram-positive and Gram-negative membranes, and an outer Gram-negative membrane. At the core of adepantin-1's behavior is its amphipathic α-helical structure, which was implemented in its design. The amphipathic structure promotes rapid self-association of peptide in water or upon binding to bacterial membranes. Aggregates initially make contact with the membrane via positively charged residues, but with insertion, the hydrophobic residues are exposed to the membrane's hydrophobic core. This adaptation alters the aggregate's stability, causing the peptides to diffuse in the polar region of the membrane, mostly remaining as a single peptide or pairing up to form an antiparallel dimer. Thus, the aggregate's proposed role is to aid in positioning the peptide into a favorable conformation for insertion. Simulations revealed the molecular basics of adepantin-1 binding to various membranes, and highlighted peptide aggregation as an important factor. These findings contribute to the development of novel anti-infective agents to combat the rapidly growing problem of bacterial resistance to antibiotics.

Interaction between Antimicrobial Peptide Magainin 2 and Nonlipid Components in the Bacterial Outer Envelope

Antimicrobial peptides (AMPs) offer advantages over conventional antibiotics; for example, bacteria develop more resistance to small-molecule antibiotics than to AMPs. The interaction of the AMPs with the lipopolysaccharide (LPS) layer of the Gram-negative bacteria cell envelope is not well understood. A MARTINI model was constructed of a Gram-negative bacterial outer membrane interacting with the AMP Magainin 2. In a 20 μs molecular dynamics (MD) simulation, the AMP diffused to the LPS layer of the cell envelope and remained there, suggesting interactions between the Magainin 2 and the LPS layer, causing the AMP to concentrate at that position. The free energy profile for the insertion of the Magainin 2 into the membrane was also calculated using umbrella sampling, which showed that the AMP positioned such that the cationic side chains of the AMP coordinated to the negatively charged phosphate groups of the LPS layer. These simulations indicate that the AMP Magainin 2 partition into the LPS layer of a bacterial membrane.

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.

Fighting Against Bacterial Lipopolysaccharide-Caused Infections through Molecular Dynamics Simulations: A Review

Lipopolysaccharide (LPS) is the primary component of the outer leaflet of Gram-negative bacterial outer membranes. LPS elicits an overwhelming immune response during infection, which can lead to life-threatening sepsis or septic shock for which no suitable treatment is available so far. As a result of the worldwide expanding multidrug-resistant bacteria, the occurrence and frequency of sepsis are expected to increase; thus, there is an urge to develop novel strategies for treating bacterial infections. In this regard, gaining an in-depth understanding about the ability of LPS to both stimulate the host immune system and interact with several molecules is crucial for fighting against LPS-caused infections and allowing for the rational design of novel antisepsis drugs, vaccines and LPS sequestration and detection methods. Molecular dynamics (MD) simulations, which are understood as being a computational microscope, have proven to be of significant value to understand LPS-related phenomena, driving and optimizing experimental research studies. In this work, a comprehensive review on the methods that can be combined with MD simulations, recently applied in LPS research, is provided. We focus especially on both enhanced sampling methods, which enable the exploration of more complex systems and access to larger time scales, and free energy calculation approaches. Thereby, apart from outlining several strategies for surmounting LPS-caused infections, this work reports the current state-of-the-art of the methods applied with MD simulations for moving a step forward in the development of such strategies.

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.

Integration of bio-responsive silver in 1D photonic crystals: towards the colorimetric detection of bacteria

The colour purity and versatility of fabrication of one-dimensional photonic crystals (1D PhCs) make them ideal candidates for colorimetric sensing of a variety of analytes. For instance, the detection of bacterial contaminants in food via colorimetric sensors can be highly appealing, as most of the existing detection techniques are in general time-consuming and the read-out requires specialised personnel. Here, we present a colorimetric sensor based on hybrid plasmonic/photonic 1D crystals. We demonstrate that the modification of the silver plasmon resonance brought about by the effective silver/bacterium interaction can be translated into the visible spectral region, producing a change in the structural colour. In addition, we observe a superior colorimetric sensitivity against the Gram negative Escherichia coli compared to the Gram positive Micrococcus luteus, a result that we attribute to the more efficient electrostatic interaction and cellular adhesion between the silver surface and the Gram-negative bacteria outer membrane. This approach demonstrates that in principle an easy colorimetric detection of bacterial contaminants can be achieved through the use of bio-responsive plasmonic materials, such as silver, whose selective electrostatic interaction with bacterial cell wall is well-known and occurs without the need of chemical functionalisation.

Novel probiotic lactic acid bacteria isolated from indigenous fermented foods from West Sumatera, Indonesia

Background and Aim:

Probiotics play an important role in maintaining a healthy gut and consequently promote good health. This study aimed to find novel probiotic lactic acid bacteria (LAB) from indigenous fermented foods of West Sumatera, Indonesia.

Materials and Methods:

This study utilized 10 LAB previously isolated from fermented buffalo milk (dadih), fermented fish (budu), and fermented cassava (tape) which have the ability to produce gamma-aminobutyric acid. The study commenced with the screening of LAB for certain properties, such as resistance to acid and bile salts, adhesion to mucosal surface, and antagonism against enteric pathogens (Escherichia coli, Salmonella Enteritidis, and Staphylococcus aureus). The promising isolates were identified through biochemical and gram staining methods.

Results:

All isolates in this study were potential novel probiotics. They survived at a pH level of 2.5 for 3 h (55.27-98.18%) and 6 h (50.98-84.91%). Survival in bile at a concentration of 0.3% was 39.90-58.61% and the survival rate was 28.38-52.11% at a concentration of 0.5%. The inhibitory diameter ranged from 8.75 to 11.54 mm for E. coli, 7.02 to 13.42 mm for S. aureus, and 12.49 to 19.00 mm for S. Enteritidis. All the isolates (84.5-92%) exhibited the ability to adhere to mucosal surfaces. This study revealed that all the isolates were potential probiotics but N16 proved to be superior because it was viable at a pH level of 2 (84.91%) and it had a good survival rate in bile salts assay (55.07%). This isolate was identified as Lactobacillus spp., Gram-positive bacilli bacteria, and tested negative in both the catalase and oxidase tests.

Conclusion:

All the isolates in this study may be used as probiotics, with isolate N16 (Lactobacillus spp.) as the most promising novel probiotic for poultry applications based on its ability to inhibit pathogenic bacteria.

Membrane targeting antimicrobial cyclic peptide nanotubes - an experimental and computational study

The search of new antibiotics, particularly with new mechanisms of action, is nowadays a very important public health issue, due to the worldwide increase of resistant pathogens. Within this effort, much research has been done on antimicrobial peptides, because having the membrane as a target, they represent a new antibiotic paradigm. Among these, cyclic peptides (CPs) made of sequences of D- and L-amino acids have emerged as a new class of potential antimicrobial peptides, due to their expected higher resistance to protease degradation. These CPs are planar structures that can form Self-assembled Cyclic Peptide Nanotubes (SCPNs), in particular in the presence of lipid membranes. Aiming at understanding their mechanism of action, we used biophysical experimental techniques (DSC and ATR-FTIR) together with Coarse-grained molecular dynamics (CG-MD) simulations, to characterize the interaction of these CPs with model membranes of different electrostatic charges' contents. DSC results revealed that the CPs show a strong interaction with negatively charged membranes, with differences in the strength of interactions depending on peptide and on membrane charge content, at odds with no or mild interactions with zwitterionic membranes. ATR-FTIR suggested that the peptides self-assemble at the membrane surface, adopting mainly a β-structure. The experiments with polarized light showed that in most cases they lie parallel to the membrane surface, but other forms and orientations are also apparent, depending on peptide structure and lipid:peptide ratio. The nanotube formation and orientation, as well as the dependence on membrane charge were also confirmed by the CG-MD simulations. These provide detail on the position and interactions, in agreement with the experimental results. Based on the findings reported here, we could proceed to the design and synthesis of a second-generation CPs, based on CP2 (soluble peptide), with increased activity and reduced toxicity.

Lipopolysaccharide Simulations Are Sensitive to Phosphate Charge and Ion Parameterization

The high proportion of lipopolysaccharide (LPS) molecules in the outer membrane of Gram-negative bacteria makes it a highly effective barrier to small molecules, antibiotic drugs, and other antimicrobial agents. Given this vital role in protecting bacteria from potentially hostile environments, simulations of LPS bilayers and outer membrane systems represent a critical tool for understanding the mechanisms of bacterial resistance and the development of new antibiotic compounds that circumvent these defenses. The basis of these simulations is parameterizations of LPS, which have been developed for all major molecular dynamics force fields. However, these parameterizations differ in both the protonation state of LPS and how LPS membranes behave in the presence of various ion species. To address these discrepancies and understand the effects of phosphate charge on bilayer properties, simulations were performed for multiple distinct LPS chemotypes with different ion parameterizations in both protonated or deprotonated lipid A states. These simulations show that bilayer properties, such as the area per lipid and inter-lipid hydrogen bonding, are highly influenced by the choice of phosphate group charges, cation type, and ion parameterization, with protonated LPS and monovalent cations with modified nonbonded parameters providing the best match to the experiments. Additionally, alchemical free energy simulations were performed to determine theoretical pK_a values for LPS and subsequently validated by ³¹P solid-state nuclear magnetic resonance experiments. Results from these complementary computational and experimental studies demonstrate that the protonated state dominates at physiological pH, contrary to the deprotonated form modeled by many LPS force fields. Overall, these results highlight the sensitivity of LPS simulations to phosphate charge and ion parameters while offering recommendations for how existing models should be updated for consistency between force fields as well as to best match experiments.

Self-Assembled Fluid Phase Floating Membranes with Tunable Water Interlayers

We present a reliable method for the fabrication of fluid phase, unsaturated lipid bilayers by self-assembly onto charged Self-Assembled Monolayer (SAM) surfaces with tunable membrane to surface aqueous interlayers. Initially, the formation of water interlayers between membranes and charged surfaces was characterized using a comparative series of bilayers deposited onto charged, self-assembled monolayers by sequential layer deposition. Using neutron reflectometry, a bilayer to surface water interlayer of ∼8 Å was found between the zwitterionic phospholipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) membrane and an anionic carboxyl terminated grafted SAM with the formation of this layer attributed to bilayer repulsion by hydration water on the SAM surface. Furthermore, we found we could significantly reduce the technical complexity of sample fabrication through self-assembly of planar membranes onto the SAM coated surfaces. Vesicle fusion onto carboxyl-terminated monolayers yielded high coverage (>95%) bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) which floated on a 7-11 Å solution interlayer between the membrane and the surface. The surface to membrane distance was then tuned via the addition of 200 mM NaCl to the bulk solution immersing a POPC floating membrane, which caused the water interlayer to swell reversibly to ∼33 Å. This study reveals that biomimetic membrane models can be readily self-assembled from solution onto functionalized surfaces without the use of polymer supports or tethers. Once assembled, surface to membrane distance can be tailored to the experimental requirements using physiological concentrations of electrolytes. These planar bilayers only very weakly interact with the substrate and are ideally suited for use as biomimetic models for accurate in vitro biochemical and biophysical studies, as well as for technological applications, such as biosensors.

Hybrid One-Dimensional Plasmonic-Photonic Crystals for Optical Detection of Bacterial Contaminants

Photonic crystal-based biosensors hold great promise as low-cost devices for real-time monitoring of a variety of biotargets, for example, bacterial contaminants in food. Here, we report the proof-of-concept for a new colorimetric sensor of bacterial contamination, which is based on a novel hybrid plasmonic-photonic device. Our system consists of a layer of silver, a plasmonic metal exhibiting a well-known bioactivity, on top of a one-dimensional photonic crystal. We attribute the bioresponsivity to the formation of polarization charges at the Ag/bacterium interface within a sort of "bio-doping" mechanism. Interestingly, this triggers a blue shift in the photonic response. As an example, we assessed the validity of our approach by detecting one of the most hazardous contaminants, Escherichia coli. This work demonstrates that our device can be a low-cost and portable platform for the detection of common bacterial contaminants.

Multiscale modeling of innate immune receptors: Endotoxin recognition and regulation by host defense peptides

The innate immune system provides a first line of defense against foreign microorganisms, and is typified by the Toll-like receptor (TLR) family. TLR4 is of particular interest, since over-stimulation of its pathway by excess lipopolysaccharide (LPS) molecules from the outer membranes of Gram-negative bacteria can result in sepsis, which causes millions of deaths each year. In this review, we outline our use of molecular simulation approaches to gain a better understanding of the determinants of LPS recognition, towards the search for novel immunotherapeutics. We first describe how atomic-resolution simulations have enabled us to elucidate the regulatory conformational changes in TLR4 associated with different LPS analogues, and hence a means to rationalize experimental structure-activity data. Furthermore, multiscale modelling strategies have provided a detailed description of the thermodynamics and intermediate structures associated with the entire TLR4 relay - which consists of a number of transient receptor/coreceptor complexes - allowing us trace the pathway of LPS transfer from bacterial membranes to the terminal receptor complex at the plasma membrane surface. Finally, we describe our efforts to leverage these computational models, in order to elucidate previously undisclosed anti-inflammatory mechanisms of endogenous host-defense peptides found in wounds. Collectively, this work represents a promising avenue for the development of novel anti-septic treatments, inspired by nature's innate defense strategies.

Outer Membrane Proteins OmpA, FhuA, OmpF, EstA, BtuB, and OmpX Have Unique Lipopolysaccharide Fingerprints

The outer membrane of Gram-negative bacteria has a highly complex asymmetrical architecture, containing a mixture of phospholipids in the inner leaflet and almost exclusively lipopolysaccharide (LPS) molecules in the outer leaflet. In E. coli, the outer membrane contains a wide range of proteins with a β barrel architecture, that vary in size from the smallest having eight strands to larger barrels composed of 22 strands. Here we report coarse-grained molecular dynamics simulations of six proteins from the E. coli outer membrane OmpA, OmpX, BtuB, FhuA, OmpF, and EstA in a range of membrane environments, which are representative of the in vivo conditions for different strains of E. coli. We show that each protein has a unique pattern of interaction with the surrounding membrane, which is influenced by the composition of the protein, the level of LPS in the outer leaflet, and the differing mobilities of the lipids in the two leaflets of the membrane. Overall we present analyses from over 200 μs of simulation for each protein.

Physical Properties of Bacterial Outer Membrane Models: Neutron Reflectometry & Molecular Simulation

The outer membrane (OM) of Gram-negative bacteria is an asymmetric bilayer having phospholipids in the inner leaflet and lipopolysaccharides in the outer leaflet. This unique asymmetry and the complex carbohydrates in lipopolysaccharides make it a daunting task to study the asymmetrical OM structure and dynamics, its interactions with OM proteins, and its roles in translocation of substrates, including antibiotics. In this study, we combine neutron reflectometry and molecular simulation to explore the physical properties of OM mimetics. There is excellent agreement between experiment and simulation, allowing experimental testing of the conclusions from simulations studies and also atomistic interpretation of the behavior of experimental model systems, such as the degree of lipid asymmetry, the lipid component (tail, head, and sugar) profiles along the bilayer normal, and lateral packing (i.e., average surface area per lipid). Therefore, the combination of both approaches provides a powerful new means to explore the biological and biophysical behavior of the bacterial OM.

Fosfomycin Permeation through the Outer Membrane Porin OmpF

Fosfomycin is a frequently prescribed drug in the treatment of acute urinary tract infections. It enters the bacterial cytoplasm and inhibits the biosynthesis of peptidoglycans by targeting the MurA enzyme. Despite extensive pharmacological studies and clinical use, the permeability of fosfomycin across the bacterial outer membrane is largely unexplored. Here, we investigate the fosfomycin permeability across the outer membrane of Gram-negative bacteria by electrophysiology experiments as well as by all-atom molecular dynamics simulations including free-energy and applied-field techniques. Notably, in an electrophysiological zero-current assay as well as in the molecular simulations, we found that fosfomycin can rapidly permeate the abundant Escherichia coli porin OmpF. Furthermore, two triple mutants in the constriction region of the porin have been investigated. The permeation rates through these mutants are slightly lower than that of the wild type but fosfomycin can still permeate. Altogether, this work unravels molecular details of fosfomycin permeation through the outer membrane porin OmpF of E. coli and moreover provides hints for understanding the translocation of phosphonic acid antibiotics through other outer membrane pores.

Molecular Dynamics Simulations Are Redefining Our View of Peptides Interacting with Biological Membranes

Ever since the first molecular mechanics computer simulations of biological molecules became possible, there has been the dream to study all complex biological phenomena in silico, simply bypassing the enormous experimental challenges and their associated costs. For this, two inherent requirements need to be met: First, the time scales achievable in simulations must reach up to the millisecond range and even longer. Second, the computational model must accurately reproduce what is measured experimentally. Despite some recent successes, the general consensus in the field to date has been that neither of these conditions have yet been met and that the dream will be realized, if at all, only in the distant future. In this Account, we show that this view is wrong; instead, we are actually in the middle of the in silico molecular dynamics (MD) revolution, which is reshaping how we think about protein function. The example explored in this Account is a recent advance in the field of membrane-active peptides (MAPs). MD simulations have succeeded in accurately capturing the process of peptide binding, folding, and partitioning into lipid bilayers as well as revealing how channels form spontaneously from polypeptide fragments and conduct ionic and other cargo across membranes, all at atomic resolution. These game-changing advances have been made possible by a combination of steadily advancing computational power, more efficient algorithms and techniques, clever accelerated sampling schemes, and thorough experimental verifications. The great advantage of MD is the spatial and temporal resolution, directly providing a molecular movie of a protein undergoing folding and cycling through a functional process. This is especially important for proteins with transitory functional states, such as pore-forming MAPs. Recent successes are demonstrated here for the large class of antimicrobial peptides (AMPs). These short peptides are an essential part of the nonadaptive immune system for many organisms, ubiquitous in nature, and of particular interest to the pharmaceutical industry in the age of rising bacterial resistance to conventional antibiotic treatments. Unlike integral membrane proteins, AMPs are sufficiently small to allow converged sampling with the unbiased high-temperature sampling methodology outlined here and are relatively easy to handle experimentally. At the same time, AMPs exhibit a wealth of complex and poorly understood interactions with lipid bilayers, which allow not only tuning and validation of the simulation methodology but also advancement of our knowledge of protein-lipid interactions at a fundamental level. Space constraints limit our discussion to AMPs, but the MD methodologies outlined here can be applied to all phenomena involving peptides in membranes, including cell-penetrating peptides, signaling peptides, viral channel forming peptides, and fusion peptides, as well as ab initio membrane protein folding and assembly. For these systems, the promise of MD simulations to predict the structure of channels and to provide complete-atomic-detail trajectories of the mechanistic processes underlying their biological functions appears to rapidly become a reality. The current challenge is to design joint experimental and computational benchmarks to verify and tune MD force fields. With this, MD will finally fulfill its promise to become an inexpensive, powerful, and easy-to-use tool providing atomic-detail insights to researchers as part of their investigations into membrane biophysics and beyond.