Coarse-grained simulations uncover Gram-negative bacterial defense against polymyxins by the outer membrane

Polymyxin B1 in the Escherichia coli inner membrane: A complex story of protein and lipopolysaccharide-mediated insertion

The rise in multi-drug resistant Gram-negative bacterial infections has led to an increased need for "last-resort" antibiotics such as polymyxins. However, the emergence of polymyxin-resistant strains threatens to bring about a post-antibiotic era. Thus, there is a need to develop new polymyxin-based antibiotics, but a lack of knowledge of the mechanism of action of polymyxins hinders such efforts. It has recently been suggested that polymyxins induce cell lysis of the Gram-negative bacterial inner membrane (IM) by targeting trace amounts of lipopolysaccharide (LPS) localized there. We use multiscale molecular dynamics (MD), including long-timescale coarse-grained (CG) and all-atom (AA) simulations, to investigate the interactions of polymyxin B1 (PMB1) with bacterial IM models containing phospholipids (PLs), small quantities of LPS, and IM proteins. LPS was observed to (transiently) phase separate from PLs at multiple LPS concentrations, and associate with proteins in the IM. PMB1 spontaneously inserted into the IM and localized at the LPS-PL interface, where it cross-linked lipid headgroups via hydrogen bonds, sampling a wide range of interfacial environments. In the presence of membrane proteins, a small number of PMB1 molecules formed interactions with them, in a manner that was modulated by local LPS molecules. Electroporation-driven translocation of PMB1 via water-filled pores was favored at the protein-PL interface, supporting the 'destabilizing' role proteins may have within the IM. Overall, this in-depth characterization of PMB1 modes of interaction reveals how small amounts of mislocalized LPS may play a role in pre-lytic targeting and provides insights that may facilitate rational improvement of polymyxin-based antibiotics.

Seeing is believing: Illuminating the Gram-negative outer membrane with molecular dynamics simulations

Recent advances in molecular dynamics (MD) simulations have led to rapid improvement in our understanding of the molecular details of the outer membranes (OMs) of Gram-negative bacteria. In this review, we highlight the latest discoveries from MD simulations of OMs, shedding light on the dynamic nature of these bacteria's first line of defense. With the focus on cutting-edge approaches, we explore the OM's sensitivity to structural features, including divalent cations and membrane composition, which have emerged as crucial determinants of antimicrobial passage. Additionally, studies have provided novel insights into outer-membrane proteins (OMPs), revealing their intricate roles in substrate translocation and their distinct interactions with lipopolysaccharides (LPS) in the OM. Finally, we explore the challenging process of β-barrel membrane protein insertion, showcasing recent findings that have enhanced our grasp of this fundamental biological phenomenon.

Molecular dynamics simulation on regulation of liquid-liquid phase separation of repetitive peptides

Understanding the intricate interactions governing protein and peptide behavior in liquid-liquid phase separation (LLPS) is crucial for unraveling biological functions and dysfunctions. This study employs a residue-leveled coarse-grained molecular dynamics approach to simulate the phase separation of repetitive polyproline and polyarginine peptides (poly PR) with varying lengths and sequences in solution, considering different concentrations and temperatures. Our findings highlight the crucial role of sequence order in promoting LLPS in peptides with identical lengths of repetitive sequences. Interestingly, repetitive peptides containing fewer than 10 polyarginine repeats exhibit no LLPS, even at salt concentrations up to 3 M. Notably, our simulations align with experimental observations, pinpointing a salt concentration of 2.7 M for PR25-induced LLPS. Utilizing the same methodology, we predict the required salt concentrations for LLPS induction as 1.2 M, 1.5 M, and 2.7 M for PR12, PR15, and PR35, respectively. These predictions demonstrate good agreement with experimental results. Extending our investigation to include the peptide glutamine and arginine (GR15) in DNA solution, our simulations mirror experimental observations of phase separation. To unveil the molecular forces steering peptide phase separation, we introduce a dielectric constant modifier and hydrophobicity disruptor into poly PR systems. Our coarse-grained analysis includes an examination of temperature effects, leading to the inference that both hydrophobic and electrostatic interactions drive phase separation in peptide systems.

A molecular overview of the polymyxin-LPS interaction in the context of its mode of action and resistance development

With the rising incidences of antimicrobial resistance (AMR) and the diminishing options of novel antimicrobial agents, it is paramount to decipher the molecular mechanisms of action and the emergence of resistance to the existing drugs. Polymyxin, a cationic antimicrobial lipopeptide, is used to treat infections by Gram-negative bacterial pathogens as a last option. Though polymyxins were identified almost seventy years back, their use has been restricted owing to toxicity issues in humans. However, their clinical use has been increasing in recent times resulting in the rise of polymyxin resistance. Moreover, the detection of "mobile colistin resistance (mcr)" genes in the environment and their spread across the globe have complicated the scenario. The mechanism of polymyxin action and the development of resistance is not thoroughly understood. Specifically, the polymyxin-bacterial lipopolysaccharide (LPS) interaction is a challenging area of investigation. The use of advanced biophysical techniques and improvement in molecular dynamics simulation approaches have furthered our understanding of this interaction, which will help develop polymyxin analogs with better bactericidal effects and lesser toxicity in the future. In this review, we have delved deeper into the mechanisms of polymyxin-LPS interactions, highlighting several models proposed, and the mechanisms of polymyxin resistance development in some of the most critical Gram-negative pathogens.

Cell cycle dependent coordination of surface layer biogenesis in Caulobacter crescentus

Surface layers (S-layers) are proteinaceous, two-dimensional paracrystalline arrays that constitute a major component of the cell envelope in many prokaryotic species. In this study, we investigated S-layer biogenesis in the bacterial model organism Caulobacter crescentus. Fluorescence microscopy revealed localised incorporation of new S-layer at the poles and mid-cell, consistent with regions of cell growth in the cell cycle. Light microscopy and electron cryotomography investigations of drug-treated bacteria revealed that localised S-layer insertion is retained when cell division is inhibited, but is disrupted upon dysregulation of MreB or lipopolysaccharide. We further uncovered that S-layer biogenesis follows new peptidoglycan synthesis and localises to regions of high cell wall turnover. Finally, correlated cryo-light microscopy and electron cryotomographic analysis of regions of S-layer insertion showed the presence of discontinuities in the hexagonal S-layer lattice, contrasting with other S-layers completed by defined symmetric defects. Our findings present insights into how C. crescentus cells form an ordered S-layer on their surface in coordination with the biogenesis of other cell envelope components.

In-vitro antimicrobial activity of AF-DP protein and in-silico approach of cell membrane disruption

Microbial resistance against common antibiotics has become one of the most serious threats to human health. The increasing statistics on this problem show the necessity of finding a way to deal with it. In recent years, antimicrobial peptides with unique properties and the capability of targeting a wide range of pathogens, have been considered as a potential for replacing common antibiotics. A small chitin-binding protein with anticandidal activity was isolated from Moringa oleifera seeds by Neto and colleagues in 2017, which very much resembled antimicrobial peptides. In this study, the antimicrobial protein 'AF-DP' was identified and characterized. AF-DP was heterologously expressed, purified, and characterized, and its 3D structure was predicted. Six molecular dynamic simulations were performed to investigate how the protein interacts with Gram-negative inner and outer, Gram-positive, fungal, cancerous, and normal mammalian membranes. Also, its antimicrobial and anticancer activity was assessed in vitro via minimum inhibition concentration (MIC) and MTT assays, respectively. This protein with 111 amino acids and a total net charge (of 10.5) has been predicted to be mainly composed of alpha helix and random coils. Its MIC affecting the growth of Escherichia coli, Staphylococcus aureus, and Candida albicans was 30 µg/ml, 100 µg/ml, and 100 µg/ml, respectively; AF-DP showed anticancer activity against MCF-7 breast cancer cell line. Scanning electron microscopic analysis confirmed the creation of pores and scratches on the surface of the bacterial membrane. The results of this research show that AF-DP can be a candidate for the production of new drugs as an AMP with antimicrobial activity.Communicated by Ramaswamy H. Sarma.

Insights into colistin-mediated fluorescence labelling of bacterial LPS

Gram-negative bacterial infections are becoming untreatable due to their ability to mutate, and the gradual development of their resistance to the available antimicrobials. In recent times colistin, a drug of last resort, started losing its efficacy towards multidrug-resistant bacterial infections. Colistin targets bacterial endotoxin lipopolysaccharides (LPS) and destabilises the cytoplasmic membrane by disrupting the outer LPS membrane. In this study, we have tried to label the bacterial LPS, the main constituent of the cytoplasmic membrane of bacterial cells, to try to understand the interaction mechanism of LPS with colistin. The chemosensor, naphthaldehyde appended furfural (NAF) that selectively recognises colistin can label LPS, by showing its fluorescence signals. The computationally derived three-dimensional structure of LPS has been introduced to speculate on the possible binding mode of colistin with LPS, and this was also thoroughly studied with the help of quantum mechanics and molecular dynamics energy minimisation. Fluorescence microscopy and FE-SEM microscopic studies were also used to observe the change in the structural morphology of colistin-sensitive and resistant Salmonella typhi in different experimental conditions.

Rational design of hybrid peptide with high antimicrobial property derived from Melittin and Lasioglossin

Hybridization of Antimicrobial peptides (AMPs) with unique abilities is now considered to improve AMPs' function as promising therapeutic candidates. In the current research, Lasioglossin (LL-III) with a high antimicrobial effect on Acinetobacter (A.) baumanni and Melittin with a high antimicrobial effect against Staphylococcus (S.) aureus were selected for designing a hybrid peptide with modified properties. In the present study, a hybrid peptide with modified properties was designed. Molecular dynamic (MD) and coarse-grained (CG) simulations were done to evaluate the stability and interaction of the hybrid peptide with related membrane models. In this study, a truncated Melittin peptide (11 amino acids) was fused to an LL-III peptide (15 amino acids) to raise the antimicrobial activity. A new hybrid peptide analog (LM1) was selected by replacing the arginine with isoleucine in the fifth position of truncated Melittin to raise the antimicrobial rate of the peptide. The potential for binding of the LM1 to lipid membrane (D factor) was increased from 2.02 related to Melittin to 3.62. Based on VMD results, the N-terminal of LM1 peptide related to LL-III was the alpha helix during 200 ns. However, the C-terminal region related to the Melittin peptide only at 50 ns was in alpha helix form. The RMSD of the LM1 peptide was in the range of 0.2 to 0.8, which, after 160 ns, reached stability. RMSD and RMSF results indicated no unwanted fluctuations during the 200 ns MD simulation. A significant movement of LM1 peptide inside the S. aureus membrane(4.76 nm) and A. baumanni membrane (3.2 nm) was observed by CG simulation. Our findings highlight the high stability of the designed hybrid peptide and its antimicrobial potential to halter A. baumanii and S. aureus bacteria.Communicated by Ramaswamy H. Sarma.

Computational investigation on lipid bilayer disruption induced by amphiphilic Janus nanoparticles: combined effect of Janus balance and charged lipid concentration

Janus nanoparticles (NPs) with charged/hydrophobic compartments have garnered attention for their potential antimicrobial activity. These NPs have been shown to disrupt lipid bilayers in experimental studies, yet the underlying mechanisms of this disruption at the particle-membrane interface remain unclear. To address this knowledge gap, the present study conducts a computational investigation to systematically examine the disruption of lipid bilayers induced by amphiphilic Janus NPs. The focus of this study is on the combined effects of the hydrophobicity of the Janus NP, referred to as the Janus balance, defined as the ratio of hydrophilic to hydrophobic surface coverage, and the concentration of charged phospholipids on the interactions between Janus NPs and lipid bilayers. Computational simulations were conducted using a coarse-grained molecular dynamics (MD) approach. The results of these MD simulations reveal that while the area change of the bilayer increases monotonically with the Janus balance, the effect of charged lipid concentration in the membrane is not easy to be predicted. Specifically, it was found that the concentration of negatively charged lipids is directly proportional to the intensity of membrane disruption. Conversely, positively charged lipids have a negligible effect on membrane defects. This study provides molecular insights into the significant role of Janus balance in the disruption of lipid bilayers by Janus NPs and supports the selectivity of Janus NPs for negatively charged lipid membranes. Furthermore, the anisotropic properties of Janus NPs were found to play a crucial role in their ability to disrupt the membrane via the combination of hydrophobic and electrostatic interactions. This finding is validated by testing the current Janus NP design on a bacterial membrane-mimicking model. This computational study may serve as a foundation for further studies aimed at optimizing the properties of Janus NPs for specific antimicrobial applications.

Polymyxin Stereochemistry and Its Role in Antibacterial Activity and Outer Membrane Disruption

With increasing rates of resistance toward commonly used antibiotics, especially among Gram-negative bacteria, there is renewed interested in polymyxins. Polymyxins are lipopeptide antibiotics with potent anti-Gram-negative activity and are generally believed to target lipid A, the lipopolysaccharide (LPS) anchor found in the outer membrane of Gram-negative bacteria. To characterize the stereochemical aspects of their mechanism(s) of action, we synthesized the full enantiomers of polymyxin B and the polymyxin B nonapeptide (PMBN). Both compounds were compared with the natural compounds in biological and biophysical assays, revealing strongly reduced antibacterial activity for the enantiomeric species. The enantiomeric compounds also exhibit reduced LPS binding, lower outer membrane (OM) permeabilization, and loss of synergetic potential. These findings provide new insights into the stereochemical requirements underlying the mechanisms of action of polymyxin B and PMBN.

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.

An Intelligent Strategy with All-Atom Molecular Dynamics Simulations for the Design of Lipopeptides against Multidrug-Resistant Pseudomonas aeruginosa

Multidrug-resistant Gram-negative bacteria seriously threaten modern medicine due to the lack of efficacious therapeutic options. Their outer membrane (OM) is an essential protective fortress to exclude many antibiotics. Unfortunately, current structural biology methods are not able to resolve the membrane structure and it is difficult to examine the specific interaction between the OM and small molecules. These limitations hinder mechanistic understanding of antibiotic penetration through the OM and antibiotic discovery. Here, we developed biologically relevant OM models by quantitatively determining membrane lipidomics of Pseudomonas aeruginosa and elucidated how lipopolysaccharide modifications and OM vesicles mediated resistance to polymyxins. Supported by chemical biology and pharmacological assays, our multiscale molecular dynamics simulations provide an intelligent platform to quantify the membrane-penetrating thermodynamics of peptides and predict their antimicrobial activity. Through experimental validations with our in-house polymyxin analogue library, our computational strategy may have significant potential in accelerating the discovery of lipopeptides against bacterial "superbugs".

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.

Interactions between polymyxin B and various bacterial membrane mimics: A molecular dynamics study

Polymyxin B (PMB) is clinically used as a last-line therapy against life-threatening Gram-negative "superbugs". However, thorough understanding of the membrane actions of PMB at a molecular level is still lacking. In this work, a variety of bacterial membrane mimics with varying lipid compositions were built, and their interactions with PMB were systematically investigated using coarse-grained molecular dynamics simulation. PMB demonstrated characteristic preference to specific lipid species during its interaction with different membrane systems, such as the rough mutant lipipolysacchrides (Re LPS) preference in an outer membrane (OM) or the cardiolipin and POPG affinity in an inner membrane (IM). As a result of the lipid-specific actions, complicated membrane interaction states of PMB were observed, including adsorption on the OM surface. In contrast, for the IM or a mutative OM containing "impurity lipids" like POPE, POPG or lipid A, it could insert into the membrane via its acyl chain. Such actions of PMB influence the structure and lipid mobility of the membrane. In particular, the OM-bound PMB breaks the synchronous movement of Re LPS molecules in the outer leaflet and makes them diffuse more randomly, while its insertion into IM blocks the phospholipid diffusion and makes the membrane more homogeneous in the trajectory space. Our results provide insight into the action mechanism of PMB at a membrane level and a foundation for developing novel and safer polymyxin strategies for better clinical use.

Polymyxin B1 within the E. coli cell envelope: insights from molecular dynamics simulations

Polymyxins are used as last-resort antibiotics, where other treatments have been ineffectual due to antibiotic resistance. However, resistance to polymyxins has also been now reported, therefore it is instructive to characterise at the molecular level, the mechanisms of action of polymyxins. Here we review insights into these mechanisms from molecular dynamics simulations and discuss the utility of simulations as a complementary technique to  experimental methodologies.

A generic force field for simulating native protein structures using dissipative particle dynamics

A coarse-grained force field for molecular dynamics simulations of the native structures of proteins in a dissipative particle dynamics (DPD) framework is developed. The parameters for bonded interactions are derived by mapping the bonds and angles for 20 amino acids onto target distributions obtained from fully atomistic simulations in explicit solvent. A dual-basin potential is introduced for stabilizing backbone angles, to cover a wide spectrum of protein secondary structures. The backbone dihedral potential enables folding of the protein from an unfolded initial state to the folded native structure. The proposed force field is validated by evaluating the structural properties of several model peptides and proteins including the SARS-CoV-2 fusion peptide, consisting of α-helices, β-sheets, loops and turns. Detailed comparisons with fully atomistic simulations are carried out to assess the ability of the proposed force field to stabilize the different secondary structures present in proteins. The compact conformations of the native states were evident from the radius of gyration and the high intensity peaks of the root mean square deviation histograms, which were found to be within 0.4 nm. The Ramachandran-like energy landscape on the phase space of backbone angles (θ) and dihedrals (ϕ) effectively captured the conformational phase space of α-helices at ∼(ϕ = 50°,θ = 90°) and β-strands at ∼(ϕ = ±180°,θ = 90-120°). Furthermore, the residue-residue native contacts were also well reproduced by the proposed DPD model. The applicability of the model to multidomain complexes was assessed using lysozyme and a large α-helical bacterial pore-forming toxin, cytolysin A. Our study illustrates that the proposed force field is generic, and can potentially be extended for efficient in silico investigations of membrane bound polypeptides and proteins using DPD simulations.