Membrane models of E. coli containing cyclic moieties in the aliphatic lipid chain

Investigation on the Synergy between Membrane Permeabilizing Amphiphilic α-Hydrazido Acids and Commonly Used Antibiotics against Drug-Resistant Bacteria

The growth of (multi)drug resistance in bacteria is among the most urgent global health issues. Monocationic amphiphilic α-hydrazido acid derivatives are structurally simple mimics of antimicrobial peptides (AMPs) with fewer drawbacks. Their mechanism of membrane permeabilization at subtoxic concentrations was found to begin with an initial electrostatic attraction of isolated amphiphile molecules to the phospholipid heads, followed by a rapid insertion of the apolar portions. As the accumulation into the bilayer proceeded, the membrane increased its fluidity and permeability without being subjected to major structural damage. After having ascertained that α-hydrazido acid amphiphiles do not interact with bacterial DNA, they were subjected to synergy evaluation for combinations with conventional antibiotics. Synergy was observed for combinations with tetracycline against sensitive S. aureus and E. coli, as well as with ciprofloxacin and colistin against resistant strains. Additivity with a remarkable recovery in activity of conventional antibiotics (from 2-fold to ≥32-fold) together with largely subtoxic concentrations of α-hydrazido acid derivatives was found for combinations with ciprofloxacin toward susceptible S. aureus and methicillin toward MRSa. However, no potentiation of conventional antibiotics was observed for combinations with linezolid and gentamicin against the corresponding resistant S. aureus and E. coli strains.

A Gram-negative-selective antibiotic that spares the gut microbiome

Infections caused by Gram-negative pathogens are increasingly prevalent and are typically treated with broad-spectrum antibiotics, resulting in disruption of the gut microbiome and susceptibility to secondary infections1-3. There is a critical need for antibiotics that are selective both for Gram-negative bacteria over Gram-positive bacteria, as well as for pathogenic bacteria over commensal bacteria. Here we report the design and discovery of lolamicin, a Gram-negative-specific antibiotic targeting the lipoprotein transport system. Lolamicin has activity against a panel of more than 130 multidrug-resistant clinical isolates, shows efficacy in multiple mouse models of acute pneumonia and septicaemia infection, and spares the gut microbiome in mice, preventing secondary infection with Clostridioides difficile. The selective killing of pathogenic Gram-negative bacteria by lolamicin is a consequence of low sequence homology for the target in pathogenic bacteria versus commensals; this doubly selective strategy can be a blueprint for the development of other microbiome-sparing antibiotics.

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.

Modulation of Phospholipase A2 Membrane Activity by Anti-inflammatory Drugs

The phospholipase A2 (PLA2) superfamily consists of lipolytic enzymes that hydrolyze specific cell membrane phospholipids and have long been considered a central hub of biosynthetic pathways, where their lipid metabolites exert a variety of physiological roles. A misregulated PLA2 activity is associated with mainly inflammatory-derived pathologies and thus has shown relevant therapeutic potential. Many natural and synthetic anti-inflammatory drugs (AIDs) have been proposed as direct modulators of PLA2 activity. However, despite the specific chemical properties that these drugs share in common, little is known about the indirect modulation able to finely tune membrane structural changes at the precise lipid-binding site. Here, we use a novel experimental strategy based on differential scanning calorimetry to systematically study the structural properties of lipid membrane systems during PLA2 cleavage and under the influence of several AIDs. For a better understanding of the AIDs-membrane interaction, we present a comprehensive and comparative set of molecular dynamics (MD) simulations. Our thermodynamic results clearly demonstrate that PLA2 cleavage is hindered by those AIDs that significantly reduce the lipid membrane cooperativity, while the rest of the AIDs oppositely tend to catalyze PLA2 activity to different extents. On the other hand, our MD simulations support experimental results by providing atomistic details on the binding, insertion, and dynamics of each AID on a pure lipid system; the drug efficacy to impact membrane cooperativity is related to the lipid order perturbation. This work suggests a membrane-based mechanism of action for diverse AIDs against PLA2 activity and provides relevant clues that must be considered in its modulation.

Effect of Glycone Diversity on the Interaction of Triterpenoid Saponins and Lipid Bilayers

Triterpenoid saponins are organic compounds widely available in the plant kingdom. These molecules have received extensive attention due to their antibacterial activity against both Gram-negative and Gram-positive bacteria. Recent studies identified the antibacterial activity of saponins closely relates to their interaction with bacterial membrane lipids; however, molecular details of this interaction remain unclear. Increased understanding of the mechanisms to disrupt bacterial lipid bilayers can help to mitigate development of antibiotic resistance. Here, we examined the effect of chemical structure and deprotonation states of saponin on its interaction with a bacterial membrane model using molecular dynamics simulations. We run multiple simulations with a ternary lipid mixture of POPE/POPG/DPPG (80/15/5 mol %) and different saponin molecules. While all saponin structures can permanently bind the membrane, their location and orientation inside the bilayer depend on the sugar chains attached to their backbone. Similarly, cluster formation and stability also depend on the chemical structure of the saponin molecule. Deprotonation site affects interactions with the bilayer by modulating hydrophilicity of the molecules. At the low concentrations simulated in this work, there is no statistically significant change in the membrane properties upon saponin(s) binding, but the molecules do preferentially partition to POPE lipid environment.

Computational microbiology of bacteria: Advancements in molecular dynamics simulations

Microbiology is traditionally considered within the context of wet laboratory methodologies. Computational techniques have a great potential to contribute to microbiology. Here, we describe our loose definition of "computational microbiology" and provide a short survey focused on molecular dynamics simulations of bacterial systems that fall within this definition. It is our contention that increased compositional complexity and realistic levels of molecular crowding within simulated systems are key for bridging the divide between experimental and computational microbiology.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Computational Analysis Reveals Unique Binding Patterns of Oxygenated and Deoxygenated Myoglobin to the Outer Mitochondrial Membrane

Myoglobin (Mb) interaction with the outer mitochondrial membrane (OMM) promotes oxygen (O₂) release. However, comprehensive molecular details on specific contact regions of the OMM with oxygenated (oxy-) and deoxygenated (deoxy-)Mb are missing. We used molecular dynamics (MD) simulations to explore the interaction of oxy- and deoxy-Mb with the membrane lipids of the OMM in two lipid compositions: (a) a typical whole membrane on average, and (b) specifically the cardiolipin-enriched cristae region (contact site). Unrestrained relaxations showed that on average, both the oxy- and deoxy-Mb established more stable contacts with the lipids typical of the cristae contact site, then with those of the average OMM. However, in steered detachment simulations, deoxy-Mb clung more tightly to the average OMM, and oxy-Mb strongly preferred the contact sites of the OMM. The MD simulation analysis further indicated that a non-specific binding, mediated by local electrostatic interactions, existed between charged or polar groups of Mb and the membrane, for stable interaction. To the best of our knowledge, this is the first computational study providing the molecular details of the direct Mb-mitochondria interaction that assisted in distinguishing the preferred localization of oxy- and deoxy-Mb on the OMM. Our findings support the existing experimental evidence on Mb-mitochondrial association and shed more insights on Mb-mediated O₂ transport for cellular bioenergetics.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Cavitation bubble interaction with compliant structures on a microscale: A contribution to the understanding of bacterial cell lysis by cavitation treatment

Numerous studies have already shown that the process of cavitation can be successfully used for water treatment and eradication of bacteria. However, most of the relevant studies are being conducted on a macro scale, so the understanding of the processes at a fundamental level remains poor. In attempt to further elucidate the process of cavitation-assisted water treatment on a scale of a single bubble, the present paper numerically addresses interaction between a collapsing microbubble and a nearby compliant structure, that mechanically and structurally resembles a bacterial cell. A fluid-structure interaction methodology is employed, where compressible multiphase flow is considered and the bacterial cell wall is modeled as a multi-layered shell structure. Simulations are performed for two selected model structures, each resembling the main structural features of Gram-negative and Gram-positive bacterial cell envelopes. The contribution of two independent dimensionless geometric parameters is investigated, namely the bubble-cell distance δ and their size ratio ς. Three characteristic modes of bubble collapse dynamics and four modes of spatiotemporal occurrence of peak local stresses in the bacterial cell membrane are identified throughout the parameter space considered. The former range from the development of a weak and thin jet away from the cell to spherical bubble collapses. The results show that local stresses arising from bubble-induced loads can exceed poration thresholds of cell membranes and that bacterial cell damage could be explained solely by mechanical effects in absence of thermal and chemical ones. Based on this, the damage potential of a single microbubble for bacteria eradication is estimated, showing a higher resistance of the Gram-positive model organism to the nearby bubble collapse. Microstreaming is identified as the primary mechanical mechanism of bacterial cell damage, which in certain cases may be enhanced by the occurrence of shock waves during bubble collapse. The results are also discussed in the scope of bacteria eradication by cavitation treatment on a macro scale, where processes of hydrodynamic and ultrasonic cavitation are being employed.

Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale

Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.

Algorithm to catalogue topologies of dynamic lipid hydrogen-bond networks

Lipid membrane interfaces host reactions essential for the functioning of cells. The hydrogen-bonding environment at the membrane interface is particularly important for binding of proteins, drug molecules, and ions. We present here the implementation and applications of a depth-first search algorithm that analyzes dynamic lipid interaction networks. Lipid hydrogen-bond networks sampled transiently during simulations of lipid bilayers are clustered according to main types of topologies that characterize three-dimensional arrangements of lipids connected to each other via short water bridges. We characterize the dynamics of hydrogen-bonded lipid clusters in simulations of model POPE and POPE:POPG membranes that are often used for bacterial membrane proteins, in a model of the Escherichia coli membrane with six different lipid types, and in POPS membranes. We find that all lipids sample dynamic hydrogen-bonded networks with linear, star, or circular arrangements of the lipid headgroups, and larger networks with combinations of these three types of topologies. Overall, linear lipid-water bridges tend to be short. Water-mediated lipid clusters in all membranes with PE lipids tend to be somewhat small, with about four lipids in all membranes studied here. POPS membranes allow circular arrangements of three POPS lipids to be sampled frequently, and complex arrangements of linear, star, and circular paths may also be sampled. These findings suggest a molecular picture of the membrane interface whereby lipid molecules transiently connect in clusters with somewhat small spatial extension.

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

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

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

Model architectures for bacterial membranes

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

Supplementary Information

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

All-Atom Modeling of Complex Cellular Membranes

Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.

Leaflet Asymmetry Modeling in the Lipid Composition of Escherichia coli Cytoplasmic Membranes

Lipid composition asymmetry between leaflets is important to cell function and plays a key role in the "positive inside" rule in transmembrane proteins. In this work, Escherichia coli inner plasma membrane models reflecting this asymmetry have been investigated at the early-log and stationary stages during the bacterial lifecycle using all-atom molecular dynamics simulations. The CHARMM36 lipid force field is used, and selected membrane properties are tested for variations between two leaflets and whole membranes. Our models include bacterial lipids with a cyclopropane moiety on the sn-2 acyl chain in the stationary membrane model. The PE/PG ratio for two leaflets reflects the "positive inside" rule of membrane proteins, set to 6.8 and 2.8 for the inner and outer leaflets of the two models, respectively. We are the first to model leaflet asymmetry in the lipid composition of E. coli cytoplasmic membranes and observe the effect on membrane properties in leaflets and whole membranes. Specifically, our results show that for the stationary phase bilayer, the surface area per lipid (SA/lipid) is larger, the thickness (2D_C and D_B) is smaller, the tilt angle is larger, the tilt modulus is smaller, and the deuterium order parameters (S_CD) of sn-1 and sn-2 tails are lower, compared to the early-log stage. Moreover, the stationary stage bilayer has a positive spontaneous curvature, while the early-log stage has a near flat spontaneous curvature. For leaflet asymmetry, the inner leaflet has a larger SA/lipid, a smaller thickness, a smaller elastic tilt modulus (a larger tilt angle), and lower S_CD, compared to the outer leaflet in both stages. Moreover, an asymmetric membrane involves a lipid tilt and a lateral extension, varying from a reference state of a pre-equilibrium membrane. This work encourages a more profound exploration of leaflet asymmetry in various other membrane models and how this might affect the structure and function of membrane-associated peptides and proteins.

Lpp positions peptidoglycan at the AcrA-TolC interface in the AcrAB-TolC multidrug efflux pump

The multidrug efflux pumps of Gram-negative bacteria are a class of complexes that span the periplasm, coupling both the inner and outer membranes to expel toxic molecules. The best-characterized example of these tripartite pumps is the AcrAB-TolC complex of Escherichia coli. However, how the complex interacts with the peptidoglycan (PG) cell wall, which is anchored to the outer membrane (OM) by Braun's lipoprotein (Lpp), is still largely unknown. In this work, we present molecular dynamics simulations of a complete, atomistic model of the AcrAB-TolC complex with the inner membrane, OM, and PG layers all present. We find that the PG localizes to the junction of AcrA and TolC, in agreement with recent cryo-tomography data. Free-energy calculations reveal that the positioning of PG is determined by the length and conformation of multiple Lpp copies anchoring it to the OM. The distance between the PG and OM measured in cryo-electron microscopy images of wild-type E. coli also agrees with the simulation-derived spacing. Sequence analysis of AcrA suggests a conserved role for interactions with PG in the assembly and stabilization of efflux pumps, one that may extend to other trans-envelope complexes as well.

Investigation of the Action of Peptides on Lipid Membranes

Calorimetric and incoherent neutron scattering methods were employed to investigate the action of magainin 2 and PGLa peptides on the phase behavior and molecular dynamics of lipids mimicking cytoplasmic membranes of Gram-negative bacteria. The impact of the peptides, tested individually and cooperatively by differential scanning calorimetry, presented a broadened peak, sometimes with a second shoulder, depicting the phase transition temperature around 21 °C. Neutron scattering revealed a small but significant variation of the membrane dynamics due to the peptides in both in-plane and out-of-plane directions. Although we did not find a clear hint for synergy in the interplay of the two peptides, the calorimetric and neutron data give compatible results in terms of a decrease of the enthalpy due to the presence of the peptides, which destabilize the membrane. The dynamics in the two directions was differentiated when the individual peptides were added to the membranes, but the impact was smaller when both peptides were added together.

Rational design of hyperstable antibacterial peptides for food preservation

We describe the design of peptides with properties like thermostability, pH stability, and antibacterial activity against a few bacterial food pathogens. Insights obtained from classical structure-function analysis of natural peptides and their mutants through antimicrobial and enzymatic assays are used to rationally develop a set of peptides. pH and thermostability assays were performed to demonstrate robust antimicrobial activity post-treatment with high temperatures and at wide pH ranges. We have also investigated the mode of action of these hyperstable peptides using membrane permeability assays, electron microscopy, and molecular dynamics simulations. Notably, through mutational studies, we show that these peptides elicit their antibacterial action via both membrane destabilization and inhibition of intracellular trypsin-the two functions attributable to separate peptide segments. Finally, toxicity studies and food preservation assays demonstrate the safety and efficacy of the designed peptides for food preservation. Overall, the study provides a general 'blueprint' for the development of stable antimicrobial peptides (AMPs). Insights obtained from this work may also be combined with combinatorial methods in high-throughput studies for future development of antimicrobials for various applications.

Amphipathic Barbiturates as Mimics of Antimicrobial Peptides and the Marine Natural Products Eusynstyelamides with Activity against Multi-resistant Clinical Isolates

We report a series of synthetic cationic amphipathic barbiturates inspired by the pharmacophore model of small antimicrobial peptides (AMPs) and the marine antimicrobials eusynstyelamides. These N,N'-dialkylated-5,5-disubstituted barbiturates consist of an achiral barbiturate scaffold with two cationic groups and two lipophilic side chains. Minimum inhibitory concentrations of 2-8 μg/mL were achieved against 30 multi-resistant clinical isolates of Gram-positive and Gram-negative bacteria, including isolates with extended spectrum β-lactamase-carbapenemase production. The guanidine barbiturate 7e (3,5-di-Br) demonstrated promising in vivo antibiotic efficacy in mice infected with clinical isolates of Escherichia coli and Klebsiella pneumoniae using a neutropenic peritonitis model. Mode of action studies showed a strong membrane disrupting effect and was supported by nuclear magnetic resonance and molecular dynamics simulations. The results express how the pharmacophore model of small AMPs and the structure of the marine eusynstyelamides can be used to design highly potent lead peptidomimetics against multi-resistant bacteria.

Algorithm to minimize MPI communications in the parallelized fast multipole method combined with molecular dynamics calculations

In the era of exascale supercomputers, large-scale, and long-time molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)-communication in the MPI-parallelized fast multipole method (FMM) combined with MD calculations under three-dimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimum-transferred data (MTD) method, should enable large-scale and long-time MD calculations to be calculated efficiently, under the condition of massive MPI-parallelization on exascale supercomputers.

Simulations of Diabetic and Non-Diabetic Peripheral Nerve Myelin Lipid Bilayers

The multilayered myelin sheath is a critical component of both central and peripheral nervous systems, forming a protective barrier against axonal damage and facilitating the movement of nervous impulses. It is primarily composed of cholesterol (CHL1), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), phosphatidylinositol (PI), sphingomyelin (SM), and galactosylceramide (GalCer) lipids. For rat sciatic nerve myelin (part of the peripheral nervous system, PNS), it has been found that cholesterol and unsaturated fatty acid contents are significantly lower in diabetic than in non-diabetic conditions. In this study, lipid compositions from experimental data are used to create four model rat sciatic nerve myelin lipid bilayers: PI-containing (non-diabetic and diabetic) and PS-containing (non-diabetic and diabetic), which were then simulated using the all-atom CHARMM36 force field. Simulation results of diabetic membranes indicate less rigid, more laterally expansive, and thinner bilayers as well as potentially reduced interactions between GalCer on opposing myelin leaflets, supporting a direct role of the cholesterol content decrease in instigating myelin deterioration and diabetic peripheral neuropathy. Compared to PI-lipids, PS-lipids were found to cause higher inter-lipid spacing and decreased order within membranes as a result of their smaller headgroup size and higher inter-lipid hydrogen bonding potential, which allow them to more frequently reside deeper in the membrane plane and produce pushing effects on other lipids. GalCer deuterium order parameters and non-diabetic headgroup-to-headgroup bilayer thicknesses were compared to experimental data, exhibiting close alignment, supporting the future usage of these models to study the PNS myelin sheath.

Evolution of the analytical scattering model of live Escherichia coli

A previously reported multi-scale model for (ultra-)small-angle X-ray (USAXS/SAXS) and (very) small-angle neutron scattering (VSANS/SANS) of live Escherichia coli was revised on the basis of compositional/metabolomic and ultrastructural constraints. The cellular body is modeled, as previously described, by an ellipsoid with multiple shells. However, scattering originating from flagella was replaced by a term accounting for the oligosaccharide cores of the lipopolysaccharide leaflet of the outer membrane including its cross-term with the cellular body. This was mainly motivated by (U)SAXS experiments showing indistinguishable scattering for bacteria in the presence and absence of flagella or fimbrae. The revised model succeeded in fitting USAXS/SAXS and differently contrasted VSANS/SANS data of E. coli ATCC 25922 over four orders of magnitude in length scale. Specifically, this approach provides detailed insight into structural features of the cellular envelope, including the distance of the inner and outer membranes, as well as the scattering length densities of all bacterial compartments. The model was also successfully applied to E. coli K12, used for the authors' original modeling, as well as for two other E. coli strains. Significant differences were detected between the different strains in terms of bacterial size, intermembrane distance and its positional fluctuations. These findings corroborate the general applicability of the approach outlined here to quantitatively study the effect of bactericidal compounds on ultrastructural features of Gram-negative bacteria without the need to resort to any invasive staining or labeling agents.

Vitamin B12-peptide nucleic acids use the BtuB receptor to pass through the Escherichia coli outer membrane

Short modified oligonucleotides that bind in a sequence-specific way to messenger RNA essential for bacterial growth could be useful to fight bacterial infections. One such promising oligonucleotide is peptide nucleic acid (PNA), a synthetic DNA analog with a peptide-like backbone. However, the limitation precluding the use of oligonucleotides, including PNA, is that bacteria do not import them from the environment. We have shown that vitamin B₁₂, which most bacteria need to take up for growth, delivers PNAs to Escherichia coli cells when covalently linked with PNAs. Vitamin B₁₂ enters E. coli via a TonB-dependent transport system and is recognized by the outer-membrane vitamin B₁₂-specific BtuB receptor. We engineered the E. coli ΔbtuB mutant and found that transport of the vitamin B₁₂-PNA conjugate requires BtuB. Thus, the conjugate follows the same route through the outer membrane as taken by free vitamin B₁₂. From enhanced sampling all-atom molecular dynamics simulations, we determined the mechanism of conjugate permeation through BtuB. BtuB is a β-barrel occluded by its luminal domain. The potential of mean force shows that conjugate passage is unidirectional and its movement into the BtuB β-barrel is energetically favorable upon luminal domain unfolding. Inside BtuB, PNA extends making its permeation mechanically feasible. BtuB extracellular loops are actively involved in transport through an induced-fit mechanism. We prove that the vitamin B₁₂ transport system can be hijacked to enable PNA delivery to E. coli cells.

Native-like membrane models of E. coli polar lipid extract shed light on the importance of lipid composition complexity

BACKGROUND

Lipid-protein interactions stabilize protein oligomers, shape their structure, and modulate their function. Whereas in vitro experiments already account for the functional importance of lipids by using natural lipid extracts, in silico methods lack behind by embedding proteins in single component lipid bilayers. However, to accurately complement in vitro experiments with molecular details at very high spatio-temporal resolution, molecular dynamics simulations have to be performed in natural(-like) lipid environments.

RESULTS

To enable more accurate MD simulations, we have prepared four membrane models of E. coli polar lipid extract, a typical model organism, each at all-atom (CHARMM36) and coarse-grained (Martini3) representations. These models contain all main lipid headgroup types of the E. coli inner membrane, i.e., phosphatidylethanolamines, phosphatidylglycerols, and cardiolipins, symmetrically distributed between the membrane leaflets. The lipid tail (un)saturation and propanylation stereochemistry represent the bacterial lipid tail composition of E. coli grown at 37∘C until 3/4 of the log growth phase. The comparison of the Simple three lipid component models to the complex 14-lipid component model Avanti over a broad range of physiologically relevant temperatures revealed that the balance of lipid tail unsaturation and propanylation in different positions and inclusion of lipid tails of various length maintain realistic values for lipid mobility, membrane area compressibility, lipid ordering, lipid volume and area, and the bilayer thickness. The only Simple model that was able to satisfactory reproduce most of the structural properties of the complex Avanti model showed worse agreement of the activation energy of basal water permeation with the here performed measurements. The Martini3 models reflect extremely well both experimental and atomistic behavior of the E. coli polar lipid extract membranes. Aquaporin-1 embedded in our native(-like) membranes causes partial lipid ordering and membrane thinning in its vicinity. Moreover, aquaporin-1 attracts and temporarily binds negatively charged lipids, mainly cardiolipins, with a distinct cardiolipin binding site in the crevice at the contact site between two monomers, most probably stabilizing the tetrameric protein assembly.

CONCLUSIONS

The here prepared and validated membrane models of E. coli polar lipids extract revealed that lipid tail complexity, in terms of double bond and cyclopropane location and varying lipid tail length, is key to stabilize membrane properties over a broad temperature range. In addition, they build a solid basis for manifold future simulation studies on more realistic lipid membranes bridging the gap between simulations and experiments.

Preparing Membrane Proteins for Simulation Using CHARMM-GUI

Molecular dynamics simulations of membrane proteins have grown dramatically in the last 20 years. Running these simulations first requires embedding the protein's three-dimensional structure in a lipid bilayer of a suitable composition, one that resembles its native environment. This step is far from trivial, especially for modeling heterogeneous mixtures of lipids. CHARMM-GUI, a webserver for simulation system preparation greatly simplifies this step, allowing for the construction of complex heterogeneous and/or asymmetric membranes. Here, we demonstrate how to use CHARMM-GUI to build the membrane for the outer-membrane protein BamA.

Interactions of an Anionic Antimicrobial Peptide with Zinc(II): Application to Bacterial Mimetic Membranes

While the majority of known antimicrobial peptides are cationic, a small number consist of short Asp-rich sequences that are anionic. These require metal ions to become biologically active. Here, we report the study of the zinc complexes of the peptide GADDDDD (GAD5), an antimicrobial peptide. Using a combination of dynamic light scattering (DLS), ζ-potential, infrared, Raman, thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), and scanning electron microscopy (SEM), we find that adding zinc ions to GAD5 forces it into a compact structure. Higher amounts of zinc ions favor a larger structure, possibly a dimer. SEM images show that zinc ions reduce the size of the fibrillar structures of GAD5. TGA curves show that the addition of zinc ions increases the thermal stability of the structure of the peptide. TGA and DSC indicate that the association of GAD5 with a zwitterionic phospholipid in the presence of zinc ions is the most stable. The stability of that complex is due to the presence of a sharp endothermic peak in the 200-300 °C range, suggesting the presence of interlamellar water that is essential to the stabilization of the structure. These results indicate that the Zn-GAD5 complex prefers the bacteria-mimicking neutral (zwitterionic) membranes. In the presence of negatively charged phospholipids, the complex remains unordered and unstable. In terms of mechanism of action, the Zn-GAD5 complex promotes a possible endocytic uptake with respect to neutral (zwitterionic) membranes while promoting membrane disruption by forming pores with respect to negatively charged phospholipids.

Molecular dynamics study of lipid bilayers modeling outer and inner leaflets of plasma membranes of mouse hepatocytes. I. Differences in physicochemical properties between the two leaflets

Outer and inner leaflets of plasma cell membranes have different lipid compositions, and the membrane properties of each leaflet can differ from each other significantly due to these composition differences. However, because of the experimental difficulty in measuring the membrane properties for each leaflet separately, the differences are not well understood at a molecular level. In this study, we constructed two lipid bilayer systems, modeling outer and inner leaflets of plasma membranes of mouse hepatocytes based on experimental composition data. The ion concentration in the interlamellar water phase was also set to match the concentration in extra- and intracellular fluids. The differences in physical properties between the outer and inner leaflets of mouse hepatocyte cell membrane models were investigated by performing 1.2 μs-long all-atomistic molecular dynamics calculations under physiological temperature and pressure conditions (310.15 K and 1 atm). The calculated electron density profiles along the bilayer normal for each model bilayer system captured well the asymmetric feature of the experimental electron density profile across actual cell plasma membranes, indicating that our procedure of modeling the outer and inner leaflets of the cell plasma membranes was satisfactory. We found that compared to the outer leaflet model, the inner leaflet model had a very bulky and soft structure in the lateral direction. To confirm the differences, membrane fluidity was measured from the lateral diffusivity and relaxation times. The fluidity was significantly higher in the inner leaflet model than in the outer leaflet model. We also discuss two topics that are of wide interest in biology, i.e., the interdigitation of acyl tails of lipid molecules between two monolayers and the lateral concentration fluctuation of lipid molecules in the bilayers.

Extracellular loops of BtuB facilitate transport of vitamin B12 through the outer membrane of E. coli

Vitamin B12 (or cobalamin) is an enzymatic cofactor essential both for mammals and bacteria. However, cobalamin can be synthesized only by few microorganisms so most bacteria need to take it up from the environment through the TonB-dependent transport system. The first stage of cobalamin import to E. coli cells occurs through the outer-membrane receptor called BtuB. Vitamin B12 binds with high affinity to the extracellular side of the BtuB protein. BtuB forms a β-barrel with inner luminal domain and extracellular loops. To mechanically allow for cobalamin passage, the luminal domain needs to partially unfold with the help of the inner-membrane TonB protein. However, the mechanism of cobalamin permeation is unknown. Using all-atom molecular dynamics, we simulated the transport of cobalamin through the BtuB receptor embedded in an asymmetric and heterogeneous E. coli outer-membrane. To enhance conformational sampling of the BtuB loops, we developed the Gaussian force-simulated annealing method (GF-SA) and coupled it with umbrella sampling. We found that cobalamin needs to rotate in order to permeate through BtuB. We showed that the mobility of BtuB extracellular loops is crucial for cobalamin binding and transport and resembles an induced-fit mechanism. Loop mobility depends not only on the position of cobalamin but also on the extension of luminal domain. We provided atomistic details of cobalamin transport through the BtuB receptor showing the essential role of the mobility of BtuB extracellular loops. A similar TonB-dependent transport system is used also by many other compounds, such as haem and siderophores, and importantly, can be hijacked by natural antibiotics. Our work could have implications for future delivery of antibiotics to bacteria using this transport system.

A Minimal Membrane Metal Transport System: Dynamics and Energetics of mer Proteins

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+ -bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free-energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+ -bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse-grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Characteristics of a Folate Receptor-α Anchored into a Multilipid Bilayer Obtained from Atomistic Molecular Dynamics Simulations

Thorough computational description of the properties of membrane-anchored protein receptors, which are important for example in the context of active targeting drug delivery, may be achieved by models representing as close as possible the immediate environment of these macromolecules. An all-atom bilayer, including 35 different lipid types asymmetrically distributed among the two monolayers, is suggested as a model neoplastic cell membrane. One molecule of folate receptor-α (FRα) is anchored into its outer leaflet, and the behavior of the system is explored by atomistic molecular dynamics simulations. The total number of atoms in the model is ∼185 000. Three 1-μs-long simulations are carried out, where physiological conditions (310 K and 1 bar) are maintained with three different pressure scaling schemes. To evaluate the structure and the phase state of the membrane, the density profiles of the system, the average area per lipid, and the deuterium order parameter of the lipid tails are calculated. The bilayer is in liquid ordered state, and the specific arrangement varies between the three trajectories. The changes in the structure of FRα are investigated and are found time- and ensemble-dependent. The volume of the ligand binding pocket fluctuates with time, but this variation remains independent of the more global structural alterations. The latter are mostly "waving" motions of the protein, which periodically approaches and retreats from the membrane. The semi-isotropic pressure scaling perturbs the receptor most significantly, while the isotropic algorithm induces rather slow changes. Maintaining constant nonzero surface tension leads to behavior closest to the experimentally observed one.

Conformational Dynamics of AcrA Govern Multidrug Efflux Pump Assembly

Multidrug efflux pumps of pathogenic, Gram-negative bacteria comprise an innate resistance mechanism and are key contributors to the emerging global pandemic of antibiotic resistance. Several increasingly detailed cryo-electron microscopy maps have been resolved of an entire efflux pump complex, AcrAB-TolC, resulting in atomistic structural models. Using a recent model, we have carried out nearly 40 μs of molecular dynamics simulations to study one of the key components of the protein complex AcrA, the membrane fusion protein that connects the inner-membrane-bound AcrB to the outer-membrane-bound TolC. We determined a three-dimensional potential of mean force (PMF) for AcrA, which displays two main conformational basins representing assembly competent and incompetent states. Corresponding experiments show that stabilizing mutations at an interdomain interface shift the dynamic equilibrium between these states to the incompetent one, disrupting pump assembly and function and resensitizing bacteria to existing antibiotics. The modulation of AcrA dynamics through pharmacological intervention therefore presents a promising route for the development of new antibiotics.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Response of microbial membranes to butanol: interdigitation vs. disorder

Biobutanol production by fermentation is potentially a sustainable alternative to butanol production from fossil fuels. However, the toxicity of butanol to fermentative bacteria, resulting largely from cell membrane fluidization, limits production titers and is a major factor limiting the uptake of the technology. Here, studies were undertaken, in vitro and in silico, on the butanol effects on a representative bacterial (i.e. Escherichia coli) inner cell membrane. A critical butanol : lipid ratio for stability of 2 : 1 was observed, computationally, consistent with complete interdigitation. However, at this ratio the bilayer was ∼20% thicker than for full interdigitation. Furthermore, butanol intercalation induced acyl chain bending and increased disorder, measured as a 27% lateral diffusivity increase experimentally in a supported lipid bilayer. There was also a monophasic Tm reduction in butanol-treated large unilamellar vesicles. Both behaviours are inconsistent with an interdigitated gel. Butanol thus causes only partial interdigitation at physiological temperatures, due to butanol accumulating at the phospholipid headgroups. Acyl tail disordering (i.e. splaying and bending) fills the subsequent voids. Finally, butanol short-circuits the bilayer and creates a coupled system where interdigitated and splayed phospholipids coexist. These findings will inform the design of strategies targeting bilayer stability for increasing biobutanol production titers.

Computational Modeling of Realistic Cell Membranes

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

Mechanisms by Which Lipids Influence Conformational Dynamics of the GlpG Intramembrane Protease

Rhomboid intramembrane proteases are bound to lipid membranes, where they dock and cleave other transmembrane substrates. How the lipid membrane surrounding the protease impacts the conformational dynamics of the protease is essential to understand because it informs on the reaction coordinate of substrate binding. Atomistic molecular dynamics simulations allow us to probe protein motions and characterize the coupling between protein and lipids. Simulations performed here on GlpG, the rhomboid protease from Escherichia coli, indicate that the thickness of the lipid membrane close to GlpG depends on both the composition of the lipid membrane and the conformation of GlpG. Transient binding of a lipid headgroup at the active site of the protease, as observed in some of the simulations reported here, suggests that a lipid headgroup might compete with the substrate for access to the GlpG active site. Interactions identified between lipid headgroups and the protein influence the dynamics of lipid interactions close to the substrate-binding site. These observations suggest that the lipid membrane environment shapes the energy profile of the substrate-docking region of the enzyme reaction coordinate.

Molecular Dynamics Investigation into the Effect of Zinc(II) on the Structure and Membrane Interactions of the Antimicrobial Peptide Clavanin A

Clavanin A (ClavA) is an antimicrobial peptide (AMP) whose antimicrobial activity is enhanced in the presence of Zn(II) ions. The antimicrobial activity of ClavA has been shown to increase 16-fold in the presence of Zn(II) ions. In this study, we investigate the potential sources of this enhancement, namely, the effect of Zn(II) binding on the helical conformation of ClavA and on the ClavA interaction with a model for gram-negative bacterial membranes. In addition, we investigate the effect of Zn(II) on the membrane mechanical properties. We employed all-atom equilibrium molecular dynamics simulations initiated from both fully helical and random coil structures of ClavA. We observe that Zn(II) can stabilize an existing helical conformation in the Zn(II)-binding region, but we do not observe induction of helical conformations in systems initiated in random coil configurations. Zn(II) binding to ClavA provides more favorable electrostatics for membrane association in the C-terminal region. This is evidenced by longer and stronger C-terminal-lipid interactions. Zn(II) is also capable of modulating the membrane properties in a manner which favors ClavA insertion and the potential for enhanced translocation into the cell. This work provides insights into the role of divalent metal cations in the antimicrobial activity of ClavA. This information can be used for the development of synthetic AMPs containing motifs that can bind metals (metalloAMPs) for therapeutic and medical purposes.

The absence of SigX results in impaired carbon metabolism and membrane fluidity in Pseudomonas aeruginosa

In Pseudomonas aeruginosa, SigX is an extra-cytoplasmic function σ factor that belongs to the cell wall stress response network. In previous studies, we made the puzzling observation that sigX mutant growth was severely affected in rich lysogeny broth (LB) but not in minimal medium. Here, through comparative transcriptomic and proteomic analysis, we show that the absence of SigX results in dysregulation of genes, whose products are mainly involved in transport, carbon and energy metabolisms. Production of most of these genes is controlled by carbon catabolite repression (CCR), a key regulatory system than ensures preferential carbon source uptake and utilization, substrate prioritization and metabolism. The strong CCR response elicited in LB was lowered in a sigX mutant, suggesting altered nutrient uptake. Since the absence of SigX affects membrane composition and fluidity, we suspected membrane changes to cause such phenotype. The detergent polysorbate 80 (PS80) can moderately destabilize the envelope resulting in non-specific increased nutrient intake. Remarkably, growth, membrane fluidity and expression of dysregulated genes in the sigX mutant strain were restored in LB supplemented with PS80. Altogether, these data suggest that SigX is indirectly involved in CCR regulation, possibly via its effects on membrane integrity and fluidity.

Distribution of mechanical stress in the Escherichia coli cell envelope

The cell envelope in Gram-negative bacteria comprises two distinct membranes with a cell wall between them. There has been a growing interest in understanding the mechanical adaptation of this cell envelope to the osmotic pressure (or turgor pressure), which is generated by the difference in the concentration of solutes between the cytoplasm and the external environment. However, it remains unexplored how the cell wall, the inner membrane (IM), and the outer membrane (OM) effectively protect the cell from this pressure by bearing the resulting surface tension, thus preventing the formation of inner membrane bulges, abnormal cell morphology, spheroplasts and cell lysis. In this study, we have used molecular dynamics (MD) simulations combined with experiments to resolve how and to what extent models of the IM, OM, and cell wall respond to changes in surface tension. We calculated the area compressibility modulus of all three components in simulations from tension-area isotherms. Experiments on monolayers mimicking individual leaflets of the IM and OM were also used to characterize their compressibility. While the membranes become softer as they expand, the cell wall exhibits significant strain stiffening at moderate to high tensions. We integrate these results into a model of the cell envelope in which the OM and cell wall share the tension at low turgor pressure (0.3 atm) but the tension in the cell wall dominates at high values (>1 atm).

Evolutionary specialization of MscCG, an MscS-like mechanosensitive channel, in amino acid transport in Corynebacterium glutamicum

MscCG, a mechanosensitive channel of Corynebacterium glutamicum provides a major export mechanism for glutamate in this Gram-positive bacterium, which has for many years been used for industrial production of glutamate and other amino acids. The functional characterization of MscCG is therefore, of great significance to understand its conductive properties for different amino acids. Here we report the first successful giant spheroplast preparation of C. glutamicum amenable to the patch clamp technique, which enabled us to investigate mechanosensitive channel activities of MscCG in the native membrane of this bacterium. Single channel recordings from these spheroplasts revealed the presence of three types of mechanosensitive channels, MscCG, MscCG2, and CgMscL, which differ largely from each other in their conductance and mechanosensitivity. MscCG has a relatively small conductance of ~340 pS followed by an intermediate MscCG2 conductance of ~1.0 nS and comparably very large conductance of 3.7 nS exhibited by CgMscL. By applying Laplace's law, we determined that very moderate membrane tension of ~5.5 mN/m was required for half activation of MscCG compared to ~12 mN/m required for half activation of both MscCG2 and CgMscL. Furthermore, by combining the micropipette aspiration technique with molecular dynamics simulations we measured mechanical properties of the C. glutamicum membrane, whose area elasticity module of K_A ≈ 15 mN/m is characteristic of a very soft membrane compared to the three times larger area expansion modulus of K_A ≈ 44 mN/m of the more elastic E. coli membrane. Moreover, we demonstrate that the "soft" properties of the C. glutamicum membrane have a significant impact on the MscCG gating characterized by a strong voltage-dependent hysteresis in the membrane of C. glutamicum compared to a complete absence of the hysteresis in the E. coli cell membrane. We thus propose that MscCG has evolved and adapted as an MscS-like channel to the mechanical properties of the C. glutamicum membrane enabling the channel to specialize in transport of amino acids such as glutamate, which are major osmolytes helping the bacterial cells survive extreme osmotic stress.

Are the short cationic lipopeptides bacterial membrane disruptors? Structure-Activity Relationship and molecular dynamic evaluation

Short cationic lipopeptides are amphiphilic molecules that exhibit antimicrobial activity mainly against Gram-positives. These compounds bind to bacterial membranes and disrupt their integrity. Here we examine the structure-activity relation (SAR) of lysine-based lipopeptides, with a prospect to rationally design more active compounds. The presented study aims to explain how antimicrobial activity of lipopeptides is affected by the charge of lipopeptide headgroup and the length of lipopeptide acyl chain. The obtained SAR models suggest that the lipophilicity of short synthetic cationic lipopeptides is the major factor that determines their antimicrobial activities. In order to link the differences in antimicrobial activity to the mechanism of action of lipopeptides containing one and two hydrophobic chains, we additionally performed molecular dynamic (MD) simulations. By using combined coarse-grained and all-atom simulations we also show that these compounds neither affect the organization of the membrane lipids nor aggregate to form separate phases. These results, along with the onset of antimicrobial activity of lipopeptides well below the critical micelle concentration (CMC), indicate that lipopeptides do not act in a simple detergent-like manner.