Molecular dynamics modeling of Pseudomonas aeruginosa outer membranes

Metabolic Changes in Pseudomonas oleovorans Isolated from Contaminated Construction Material Exposed to Varied Biocide Treatments

Biocide resistance poses a significant challenge in industrial processes, with bacteria like Pseudomonas oleovorans exhibiting intrinsic resistance to traditional antimicrobial agents. In this study, the impact of biocide exposure on the metabolome of two P. oleovorans strains, namely, P. oleovorans P4A, isolated from contaminated coating material, and P. oleovorans 1045 reference strain, were investigated. The strains were exposed to 2-Methylisothiazol-3(2H)-one (MI) MIT, 1,2-Benzisothiazol-3(2H)-one (BIT), and 5-chloro-2-methyl-isothiazol-3-one (CMIT) at two different sub-inhibitory concentrations and the lipids and polar and semipolar metabolites were analyzed by ultra-high performance liquid chromatography quadrupole time-of-flight mass spectrometry UPLC-Q-TOF/MS. Exposure to the BIT biocide induced significant metabolic modifications in P. oleovorans. Notable changes were observed in lipid and metabolite profiles, particularly in phospholipids, amino acid metabolism, and pathways related to stress response and adaptation. The 1045 strain showed more pronounced metabolic alterations than the P4A strain, suggesting potential implications for lipid, amino acid metabolism, energy metabolism, and stress adaptation. Improving our understanding of how different substances interact with bacteria is crucial for making antimicrobial chemicals more effective and addressing the challenges of resistance. We observed that different biocides trigged significantly different metabolic responses in these strains. Our study shows that metabolomics can be used as a tool for the investigation of metabolic mechanisms underlying biocide resistance, and thus in the development of targeted biocides. This in turn can have implications in combating biocide resistance in bacteria such as P. oleovorans.

Discovery of antibiotics that selectively kill metabolically dormant bacteria

There is a need to discover and develop non-toxic antibiotics that are effective against metabolically dormant bacteria, which underlie chronic infections and promote antibiotic resistance. Traditional antibiotic discovery has historically favored compounds effective against actively metabolizing cells, a property that is not predictive of efficacy in metabolically inactive contexts. Here, we combine a stationary-phase screening method with deep learning-powered virtual screens and toxicity filtering to discover compounds with lethality against metabolically dormant bacteria and favorable toxicity profiles. The most potent and structurally distinct compound without any obvious mechanistic liability was semapimod, an anti-inflammatory drug effective against stationary-phase E. coli and A. baumannii. Integrating microbiological assays, biochemical measurements, and single-cell microscopy, we show that semapimod selectively disrupts and permeabilizes the bacterial outer membrane by binding lipopolysaccharide. This work illustrates the value of harnessing non-traditional screening methods and deep learning models to identify non-toxic antibacterial compounds that are effective in infection-relevant contexts.

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.

The antibacterial properties of branched peptides based on poly(l-arginine): In vitro antibacterial evaluation and molecular dynamic simulations

The emergence of bacterial strains resistant to antibiotics is a major issue in the medical field. Antimicrobial peptides are widely studied as they do not generate as much resistant bacterial strains as conventional antibiotics and present a broad range of activity. Among them, the homopolypeptide poly(l-arginine) presents promising antibacterial properties, especially in the perspective of its use in biomaterials. Linear poly(l-arginine) has been extensively studied but the impact of its 3D structure remains unknown. In this study, the antibacterial properties of newly synthesized branched poly(l-arginine) peptides, belonging to the family of multiple antigenic peptides, are evaluated. First, in vitro activities of the peptides shows that branched poly(l-arginine) is more efficient than linear poly(l-arginine) containing the same number of arginine residues. Surprisingly, peptides with more arms and more residues are not the most effective. To better understand these unexpected results, interactions between these peptides and the membranes of Gram positive and Gram negative bacteria are simulated thanks to molecular dynamic. It is observed that the bacterial membrane is more distorted by the branched structure than by the linear one and by peptides containing smaller arms. This mechanism of action is in full agreement with in vitro results and suggest that our simulations form a robust model to evaluate peptide efficiency towards pathogenic bacteria.

Structure-Activity Relationship Study to Develop Peptide Amphiphiles as Species-Specific Antimicrobials

Antimicrobial peptide amphiphiles (PAs) are a promising class of molecules that can disrupt the bacterial membrane or act as drug nanocarriers. In this study, we prepared 33 PAs to establish supramolecular structure-activity relationships. We studied the morphology and activity of the nanostructures against different Gram-positive and Gram-negative bacterial strains (such as Staphylococcus aureus, Escherichia coli, Pseudomonas aeruginosa and Acinetobacter baumannii). Next, we used principal component analysis (PCA) to determine the key contributors to activity. We found that for S. aureus, the zeta potential was the major contributor to the activity while Gram-negative bacteria were more influenced by the partition coefficient (LogP) with the following order P. aeruginosa>E. coli>A. baumannii. We also performed a study of the mechanism of action of selected PAs on the bacterial membrane assessing the membrane permeability and depolarization, changes in zeta potential and overall integrity. We studied the toxicity of the nanostructures against mammalian cells. Finally, we performed an in vivo study using the wax moth larvae to determine the therapeutic efficacy of the active PAs. This study shows cationic PA nanostructures can be an intriguing platform for the development of nanoantibacterials.

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

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

Modeling the Effect on a Novel Fungal Peptaibol Placed in an All-Atom Bacterial Membrane Mimicking System via Accelerated Molecular Dynamics Simulations

We previously reported on a novel peptaibol, named Tripleurin XIIc (TPN), an 18-residue long sequence produced by the fungus Trichoderma pleuroti. We elucidated its 3D structure via classical and accelerated molecular dynamics simulation (aMD) methods and reported the folding dynamics of TPN in water and chloroform solvents. Peptaibols, in general, are insoluble in water, as they are amphipathic and may prefer hydrophobic environments like transmembrane regions. In this study, we attempted to use aMD simulations to model an all-atom bacterial membrane system while placing a TPN molecule in its vicinity. The results highlighted that TPN was able to introduce some disorder into the membrane and caused lipid clustering. It could also enter the transmembrane region from the water-bilayer interface. The structural dynamics of TPN in the transmembrane region revealed a single energetically stable conformation similar to the one obtained from water and chloroform solvent simulations reported by us previously. However, this linear structure was found to be at the local energy minimum (stable) in water but at a metastable intermediate state (higher energy) in chloroform. Therefore, it could be said that the water solvent can be successfully used for folding simulations of peptaibols.

A new nanobiotic: synthesis and characterization of an albumin nanoparticle with intrinsic antibiotic activity

Background and Objectives

With entering the "post-antibiotic era", antibiotic resistance is one of the most important problems in food security, health and medicine. Invention of nanoparticles with intrinsic antimicrobial activity has been provided a new tool to combat the problem, including some metal nanoparticles. But protein nanoparticles have been often used as nano-carrier for antibiotic drugs, not for their own antibiotic activity. In this article we have fabricated a very small BSA-NP without any chemical modification on BSA molecules showing antibacterial activity.

Materials and Methods

Bovine serum albumin nanoparticle (BSA-NP) was synthesized using botton-up approach, by dissolution of BSA in urea-containing Tris buffer for 60 min at 60°C. Then, the BSA solution was dialyzed against distilled water in order to nanoparticle formation. The resulted BSA-NP has been characterized by dynamic light scattering (DLS), field emission surface electron microscopy (FESEM), SDS-PAGE, Fourier transform infrared spectroscopy (FTIR) and UV-spectrophotometery. Minimum inhibitory concentration (MIC) method was used for evaluation of antibacterial activity of BSA-NP against Staphylococcus aureus and Pseudomonas aeruginosa.

Results

The results obtained by DLS technique indicated that BSA molecules were self-assembled into small aggregates with a hydrodynamic diameter of 23.23 ± 2.1 nm. With a small polydispersity index (PDI=0.522), the nanoparticles had good spherical uniformity. The nanoparticles made from a single type of protein molecule (BSA) and have a relatively transparent appearance. The BSA-NPs caused a decrease in cell growth of both P. aeruginosa and S. aureus. Moreover, they had a bacteriostatic effect on P. aeruginosa (MIC=112×10-5 μM).

Conclusion

In this study, using a green synthesis method, we succeeded in synthesizing a very small uniform BSA nanoparticles without any chemical modification on BSA molecules. It also has bacteriostatic properties against P. aeruginosa. Therefore, it is hypothesized that our BSA-NPs may be effective as a new approach to combat antibiotic resistance.

Modeling and Simulation of Bacterial Outer Membranes with Lipopolysaccharides and Capsular Polysaccharides

Capsule is one of the common virulence factors in Gram-negative bacteria protecting pathogens from host defenses and consists of long-chain capsular polysaccharides (CPS) anchored in the outer membrane (OM). Elucidating structural properties of CPS is important to understand its biological functions as well as the OM properties. However, the outer leaflet of the OM in current simulation studies is represented exclusively by LPS due to the complexity and diversity of CPS. In this work, representative Escherichia coli CPS, K_LPS (a lipid A-linked form) and K_PG (a phosphatidylglycerol-linked form), are modeled and incorporated into various symmetric bilayers with co-existing LPS in different ratios. All-atom molecular dynamics simulations of these systems have been conducted to characterize various bilayer properties. Incorporation of K_LPS makes the acyl chains of LPS more rigid and ordered, while incorporation of K_PG makes them less ordered and flexible. These results are consistent with the calculated area per lipid (APL) of LPS, in which the APL of LPS becomes smaller when K_LPS is incorporated, whereas it gets larger when K_PG is included. Torsional analysis reveals that the influence of the CPS presence on the conformational distributions of the glycosidic linkages of LPS is small, and minor differences are also detected for the inner and outer regions of the CPS. Combined with previously modeled enterobacterial common antigens (ECAs) in the form of mixed bilayers, this work provides more realistic OM models as well as the basis for characterization of interactions between the OM and OM proteins.

Effect of Local Stress on Accurate Modeling of Bacterial Outer Membranes Using All-Atom Molecular Dynamics

Biological membranes are fundamental components of living organisms that play an undeniable role in their survival. Molecular dynamics (MD) serves as an essential computational tool for studying biomembranes on molecular and atomistic scales. The status quo of MD simulations of biomembranes studies a nanometer-sized membrane patch periodically extended under periodic boundary conditions (PBCs). In nature, membranes are usually composed of different lipids in their two layers (referred to as leaflets). This compositional asymmetry imposes a fixed ratio of lipid numbers between the two leaflets in a periodically constrained membrane, which needs to be set appropriately. The widely adopted methods of defining a leaflet lipid ratio suffer from the lack of control over the mechanical tension of each leaflet, which could significantly influence research findings. In this study, we investigate the role of membrane-building protocol and the resulting initial stress state on the interaction between small molecules and asymmetric membranes. We model the outer membrane of Pseudomonas aeruginosa bacteria using two different building protocols and probe their interactions with the Pseudomonas quinolone signal (PQS). Our results show that differential stress could shift the position of free energy minimum for the PQS molecule between the two leaflets of the asymmetric membrane. This work provides critical insights into the relationship between the initial per-leaflet tension and the spontaneous intercalation of PQS.

Changes in fluidity of the E. coli outer membrane in response to temperature, divalent cations and polymyxin-B show two different mechanisms of membrane fluidity adaptation

The outer membrane (OM) is an essential component of the Gram-negative bacterial cell envelope. Restricted diffusion of integral OM proteins and lipopolysaccharide (LPS) that constitute the outer leaflet of the OM support a model in which the OM is in a semi-crystalline state. The low fluidity of the OM has been suggested to be an important property of this membrane that even contributes to cell rigidity. The LPS characteristics strongly determine the properties of the OM and the LPS layer fluidity has been measured using different techniques that require specific conditions or are technically challenging. Here, we characterize the Escherichia coli LPS fluidity by evaluating the lateral diffusion of the styryl dye FM4-64FX in fluorescence recovery after photobleaching experiments. This technique allowed us to determine the effect of different conditions and genetic backgrounds on the LPS fluidity. Our results show that a fraction of the LPS can slowly diffuse and that the fluidity of the LPS layer adapts by modifying the diffusion of the LPS and the fraction of mobile LPS molecules.

Characterization of folic acid-functionalized PLA-PEG nanomicelle to deliver Letrozole: A nanoinformatics study

Effective cancer treatment is currently the number one challenge to human health. To date, several treatment methods have been introduced for cancer cell targeting. Among the proposed new methods for attacking cancer cells, nanotechnology has attracted much attention. Hence, various nanocarriers have been developed for targeted delivery of available drugs and improve their effectiveness against malignant cells. The PLA-PEG functionalised with folic acid (PLA-PEG-FA) is one of the nanocarriers with a limited range of applications for targeting cancer cells. In this investigation, different types of in-silico methods such as molecular docking approach, molecular dynamics simulation and free energy calculations are employed to characterise the carriers studied. The effectiveness of PLA-PEG-FA and PLA-PEG in delivering Letrozole as an aromatase inhibitor in cancer cells is examined. It is found that in the presence of folic acid, the stability and cell membrane permeability of nanomicelle are increased. Therefore, PLA-PEG-FA can be considered as a versatile carrier that can increase the effectiveness of aromatase inhibitors (such as Letrozole) and reduce their side effects.

Microbiota-derived extracellular vesicles in interkingdom communication in the gut

The intestine is fundamental in controlling human health. Intestinal epithelial and immune cells are continuously exposed to millions of microbes that greatly impact on intestinal epithelial barrier and immune function. This microbial community, known as gut microbiota, is now recognized as an important partner of the human being that actively contribute to essential functions of the intestine but also of distal organs. In the gut ecosystem, bidirectional microbiota-host communication does not involve direct cell contacts. Both microbiota and host-derived extracellular vesicles (EVs) are key players of such interkingdom crosstalk. There is now accumulating body of evidence that bacterial secreted vesicles mediate microbiota functions by transporting and delivering into host cells effector molecules that modulate host signalling pathways and cell processes. Consequently, vesicles released by the gut microbiota may have great influence on health and disease. Here we review current knowledge on microbiota EVs and specifically highlight their role in controlling host metabolism, intestinal barrier integrity and immune training.

Iron coordination to pyochelin siderophore influences dynamics of FptA receptor from Pseudomonas aeruginosa: a molecular dynamics simulation study

FptA is a TonB-dependent transporter that permits the high affinity binding and transport of Fe(III)-pyochelin complex across the outer membrane of Pseudomonas aeruginosa. Molecular dynamics simulations were employed to FptA receptor and its complexes with pyochelin, and co-crystallized Fe(III)-pyochelin-ethanediol and Fe(III)-pyochelin-water embedded in dilauroyl phosphatidyl choline bilayer for the evaluation of their structural and dynamical properties. The evaluation of properties of the receptor bound to pyochelin molecule and Fe(III)-pyochelin complexes helped to figure out the iron coordination effect on the receptor properties. Moreover, comparison of these four simulation systems revealed further information on the dynamical changes occurred in extracellular loops, in particular loop-7 corresponding to the missing amino acid residues including the close-by loop-8 that was largely affected by the metal coordination to pyochelin. The binding of iron to pyochelin molecule affected the overall structure of the receptor therefore, evaluation fo the gyration radii and hydrogen bonding were evaluated as well as analysis of the pore size were also carried out to understand the effect of metal coordination on the dynamics of the helices which form a kind of translocation channel to transport the siderophore across the FptA protein into the periplasmic space. The properties of each component of the molecular systems were therefore observed to be perturbed by the incorporation of iron to the pyochelin molecule thus demonstrating that the bacteria use its receptor to abstract and transport iron from extracellular environment for its survival and that was made possible to understand at the molecular level through successful implementation of molecular dynamics simulations.

Gene Transfer Potential of Outer Membrane Vesicles of Gram-Negative Bacteria

The increasing spread of multidrug-resistant pathogenic bacteria is one of the major threats to public health worldwide. Bacteria can acquire antibiotic resistance and virulence genes through horizontal gene transfer (HGT). A novel horizontal gene transfer mechanism mediated by outer membrane vesicles (OMVs) has been recently identified. OMVs are rounded nanostructures released during their growth by Gram-negative bacteria. Biologically active toxins and virulence factors are often entrapped within these vesicles that behave as molecular carriers. Recently, OMVs have been reported to contain DNA molecules, but little is known about the vesicle packaging, release, and transfer mechanisms. The present review highlights the role of OMVs in HGT processes in Gram-negative bacteria.

Amphiphilic Aminoglycosides as Medicinal Agents

The conjugation of hydrophobic group(s) to the polycationic hydrophilic core of the antibiotic drugs aminoglycosides (AGs), targeting ribosomal RNA, has led to the development of amphiphilic aminoglycosides (AAGs). These drugs exhibit numerous biological effects, including good antibacterial effects against susceptible and multidrug-resistant bacteria due to the targeting of bacterial membranes. In the first part of this review, we summarize our work in identifying and developing broad-spectrum antibacterial AAGs that constitute a new class of antibiotic agents acting on bacterial membranes. The target-shift strongly improves antibiotic activity against bacterial strains that are resistant to the parent AG drugs and to antibiotic drugs of other classes, and renders the emergence of resistant Pseudomonas aeruginosa strains highly difficult. Structure–activity and structure–eukaryotic cytotoxicity relationships, specificity and barriers that need to be crossed in their development as antibacterial agents are delineated, with a focus on their targets in membranes, lipopolysaccharides (LPS) and cardiolipin (CL), and the corresponding mode of action against Gram-negative bacteria. At the end of the first part, we summarize the other recent advances in the field of antibacterial AAGs, mainly published since 2016, with an emphasis on the emerging AAGs which are made of an AG core conjugated to an adjuvant or an antibiotic drug of another class (antibiotic hybrids). In the second part, we briefly illustrate other biological and biochemical effects of AAGs, i.e., their antifungal activity, their use as delivery vehicles of nucleic acids, of short peptide (polyamide) nucleic acids (PNAs) and of drugs, as well as their ability to cleave DNA at abasic sites and to inhibit the functioning of connexin hemichannels. Finally, we discuss some aspects of structure–activity relationships in order to explain and improve the target selectivity of AAGs.

Role of the lipid bilayer in outer membrane protein folding in Gram-negative bacteria

β-Barrel outer membrane proteins (OMPs) represent the major proteinaceous component of the outer membrane (OM) of Gram-negative bacteria. These proteins perform key roles in cell structure and morphology, nutrient acquisition, colonization and invasion, and protection against external toxic threats such as antibiotics. To become functional, OMPs must fold and insert into a crowded and asymmetric OM that lacks much freely accessible lipid. This feat is accomplished in the absence of an external energy source and is thought to be driven by the high thermodynamic stability of folded OMPs in the OM. With such a stable fold, the challenge that bacteria face in assembling OMPs into the OM is how to overcome the initial energy barrier of membrane insertion. In this review, we highlight the roles of the lipid environment and the OM in modulating the OMP-folding landscape and discuss the factors that guide folding in vitro and in vivo We particularly focus on the composition, architecture, and physical properties of the OM and how an understanding of the folding properties of OMPs in vitro can help explain the challenges they encounter during folding in vivo Current models of OMP biogenesis in the cellular environment are still in flux, but the stakes for improving the accuracy of these models are high. OMP folding is an essential process in all Gram-negative bacteria, and considering the looming crisis of widespread microbial drug resistance it is an attractive target. To bring down this vital OMP-supported barrier to antibiotics, we must first understand how bacterial cells build it.

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.

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.

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.

Accelerated Molecular Dynamics Applied to the Peptaibol Folding Problem

The use of enhanced sampling molecular dynamics simulations to facilitate the folding of proteins is a relatively new approach which has quickly gained momentum in recent years. Accelerated molecular dynamics (aMD) can elucidate the dynamic path from the unfolded state to the near-native state, “flattened” by introducing a non-negative boost to the potential. Alamethicin F30/3 (Alm F30/3), chosen in this study, belongs to the class of peptaibols that are 7–20 residue long, non-ribosomally synthesized, amphipathic molecules that show interesting membrane perturbing activity. The recent studies undertaken on the Alm molecules and their transmembrane channels have been reviewed. Three consecutive simulations of ~900 ns each were carried out where N-terminal folding could be observed within the first 100 ns, while C-terminal folding could only be achieved almost after 800 ns. It took ~1 μs to attain the near-native conformation with stronger potential boost which may take several μs worth of classical MD to produce the same results. The Alm F30/3 hexamer channel was also simulated in an E. coli mimicking membrane under an external electric field that correlates with previous experiments. It can be concluded that aMD simulation techniques are suited to elucidate peptaibol structures and to understand their folding dynamics.

Dewetting-induced formation and mechanical properties of synthetic bacterial outer membrane models (GUVs) with controlled inner-leaflet lipid composition

The double-membrane cellular envelope of Gram-negative bacteria enables them to endure harsh environments and represents a barrier to many clinically available antibiotics. The outer membrane (OM) is exposed to the environment and is the first point of contact involved in bacterial processes such as signaling, pathogenesis, and motility. As in the cytoplasmic membrane, the OM in Gram-negative bacteria has a phospholipid-rich inner leaflet and an outer leaflet that is predominantly composed of lipopolysaccharide (LPS). We report on a microfluidic technique for fabricating monodisperse asymmetric giant unilamellar vesicles (GUVs) possessing the Gram-negative bacterial OM lipid composition. Our continuous microfluidic fabrication technique generates 50-150 μm diameter water-in-oil-in-water double emulsions at high-throughput. The water-oil and oil-water interfaces facilitate the self-assembly of phospholipid and LPS molecules to create the inner and outer leaflets of the lipid bilayer, respectively. The double emulsions have ultrathin oil shells, which minimizes the amount of residual organic solvent that remains trapped between the leaflets of the GUV membrane. An extraction process by ethanol and micropipette aspiration of the ultrathin oil shells triggers an adhesive interaction between the two lipid monolayers assembled on the water-oil and oil-water interfaces (i.e., dewetting transition), forcing them to contact and form a lipid bilayer membrane. The effect of different inner-leaflet lipid compositions on the emulsion/vesicle stability and the dewetting transition is investigated. We also report on the values for bending and area expansion moduli of synthetic asymmetric model membranes with lipid composition/architecture that is physiologically relevant to the OM in Pseudomonas aeruginosa bacteria.

Molecular conformation affects the interaction of the Pseudomonas quinolone signal with the bacterial outer membrane

Gram-negative bacteria produce outer-membrane vesicles (OMVs) that package genetic elements, virulence factors, and cell-to-cell communication signaling compounds. Despite their importance in many disease-related processes, how these versatile structures are formed is incompletely understood. A self-produced secreted small molecule, the Pseudomonas quinolone signal (PQS), has been shown to initiate OMV formation in Pseudomonas aeruginosa by interacting with the outer membrane (OM) and inducing its curvature. Other bacterial species have also been shown to respond to PQS, supporting a common biophysical mechanism. Here, we conducted molecular dynamics simulations to elucidate the specific interactions between PQS and a model P. aeruginosa OM at the atomistic scale. We discovered two characteristic states of PQS interacting with the biologically relevant membrane, namely attachment to the membrane surface and insertion into the lipid A leaflet. The hydrogen bonds between PQS and the lipid A phosphates drove the PQS-membrane association. An analysis of PQS trajectory and molecular conformation revealed sequential events critical for spontaneous insertion, including probing, docking, folding, and insertion. Remarkably, PQS bent its hydrophobic side chain into a closed conformation to lower the energy barrier for penetration through the hydrophilic headgroup zone of the lipid A leaflet, which was confirmed by the potential of mean force (PMF) measurements. Attachment and insertion were simultaneously observed in the simulation with multiple PQS molecules. Our findings uncover a sequence of molecular interactions that drive PQS insertion into the bacterial OM and provide important insight into the biophysical mechanism of small molecule-induced OMV biogenesis.