Multifaceted Activity of Fabimycin: Insights from Molecular Dynamics Studies on Bacterial Membrane Models

Membranes─cells' essential scaffolds─are valid molecular targets for substances with an antimicrobial effect. While certain substances, such as octenidine, have been developed to target membranes for antimicrobial purposes, the recently reported molecule, fabimycin (F2B)─a novel agent targeting drug-resistant Gram-negative bacteria─has not received adequate attention regarding its activity on membranes in the literature. The following study aims to investigate the effects of F2B on different bacterial membrane models, including simple planar bilayers and more complex bilayer systems that mimic the Escherichia coli shell equipped with double inner and outer bilayers. Our results show that F2B exhibited more pronounced interactions with bacterial membrane systems compared to the control PC system. Furthermore, we observed significant changes in local membrane property homeostasis in both the inner and outer membrane models, specifically in the case of lateral diffusion, membrane thickness, and membrane resilience (compressibility, tilt). Finally, our results showed that the effect of F2B differed in a complex system and a single membrane system. Our study provides new insights into the multifaceted activity of F2B, demonstrating its potential to disrupt bacterial membrane homeostasis, indicating that its activity extends the currently known mechanism of FabI enzyme inhibition. This disruption, coupled with the ability of F2B to penetrate the outer membrane layers, sheds new light on the behavior of this antimicrobial molecule. This highlights the importance of the interaction with the membrane, crucial in combating bacterial infections, particularly those caused by drug-resistant strains.

The molecular scale mechanism of deposition ice nucleation on silver iodide

Heterogeneous ice nucleation is a ubiquitous process in the natural and built environment. Deposition ice nucleation, i.e. heterogeneous ice nucleation that - according to the traditional view - occurs in a subsaturated water vapor environment and in the absence of supercooled water on the solid, ice-forming surface, is among the most important ice formation processes in high-altitude cirrus and mixed-phase clouds. Despite its importance, very little is known about the mechanism of deposition ice nucleation at the microscopic level. This study puts forward an adsorption-based mechanism for deposition ice nucleation through results from a combination of atomistic simulations, experiments and theoretical modelling. One of the most potent laboratory surrogates of ice nucleating particles, silver iodide, is used as a substrate for the simulations. We find that water initially adsorbs in clusters which merge and grow over time to form layers of supercooled water. Ice nucleation on silver iodide requires at minimum the adsorption of 4 molecular layers of water. Guided by the simulations we propose the following fundamental freezing steps: (1) Water molecules adsorb on the surface, forming nanodroplets. (2) The supercooled water nanodroplets merge into a continuous multilayer when they grow to about 3 molecular layers thick. (3) The layer continues to grow until the critical thickness for freezing is reached. (4) The critical ice cluster continues to grow.

Understanding the Molecular Mechanism of Anesthesia: Effect of General Anesthetics and Structurally Similar Non-Anesthetics on the Properties of Lipid Membranes

General anesthesia can be caused by various, chemically very different molecules, while several other molecules, many of which are structurally rather similar to them, do not exhibit anesthetic effects at all. To understand the origin of this difference and shed some light on the molecular mechanism of general anesthesia, we report here molecular dynamics simulations of the neat dipalmitoylphosphatidylcholine (DPPC) membrane as well as DPPC membranes containing the anesthetics diethyl ether and chloroform and the structurally similar non-anesthetics n-pentane and carbon tetrachloride, respectively. To also account for the pressure reversal of anesthesia, these simulations are performed both at 1 bar and at 600 bar. Our results indicate that all solutes considered prefer to stay both in the middle of the membrane and close to the boundary of the hydrocarbon domain, at the vicinity of the crowded region of the polar headgroups. However, this latter preference is considerably stronger for the (weakly polar) anesthetics than for the (apolar) non-anesthetics. Anesthetics staying in this outer preferred position increase the lateral separation between the lipid molecules, giving rise to a decrease of the lateral density. The lower lateral density leads to an increased mobility of the DPPC molecules, a decreased order of their tails, an increase of the free volume around this outer preferred position, and a decrease of the lateral pressure at the hydrocarbon side of the apolar/polar interface, a change that might well be in a causal relation with the occurrence of the anesthetic effect. All these changes are clearly reverted by the increase of pressure. Furthermore, non-anesthetics occur in this outer preferred position in a considerably smaller concentration and hence either induce such changes in a much weaker form or do not induce them at all.

Anesthetics significantly increase the amount of intramembrane water in lipid membranes

The potency of anesthesia was directly linked to the partitioning of the drug molecules in cell membranes by Meyer and Overton. Many molecules interact with lipid bilayers and lead to structural and functional changes. It remains an open question which change in membrane properties is responsible for a potential anesthetic effect or if anesthetics act by binding to direct targets. We studied the effect of ethanol, diethyl ether and isoflurane on the water distribution in lipid bilayers by combining all-atom molecular dynamics simulations and neutron diffraction experiments. The simulations show strong membrane-drug interactions with partitioning coefficients of 38%, 92% and 100% for ethanol, diethyl ether and isoflurane, respectively, and provide evidence for an increased water partitioning in the membrane core. The amount of intramembrane water molecules was experimentally determined by selectively deuterium labeling lipids, anesthetic drug and water molecules in neutron diffraction experiments. Four additional water molecules per lipid were observed in the presence of ethanol. Diethyl ether and isoflurane were found to significantly increase the amount of intramembrane water by 25% (8 water molecules). This increase in intramembrane water may contribute to the non-specific interactions between anesthetics and lipid membranes.

Microcavity-Supported Lipid Bilayers; Evaluation of Drug-Lipid Membrane Interactions by Electrochemical Impedance and Fluorescence Correlation Spectroscopy

Many drugs have intracellular or membrane-associated targets, thus understanding their interaction with the cell membrane is of value in drug development. Cell-free tools used to predict membrane interactions should replicate the molecular organization of the membrane. Microcavity array-supported lipid bilayer (MSLB) platforms are versatile biophysical models of the cell membrane that combine liposome-like membrane fluidity with stability and addressability. We used an MSLB herein to interrogate drug-membrane interactions across seven drugs from different classes, including nonsteroidal anti-inflammatories: ibuprofen (Ibu) and diclofenac (Dic); antibiotics: rifampicin (Rif), levofloxacin (Levo), and pefloxacin (Pef); and bisphosphonates: alendronate (Ale) and clodronate (Clo). Fluorescence lifetime correlation spectroscopy (FLCS) and electrochemical impedance spectroscopy (EIS) were used to evaluate the impact of drug on 1,2-dioleyl- sn-glycerophosphocholine and binary bilayers over physiologically relevant drug concentrations. Although FLCS data revealed Ibu, Levo, Pef, Ale, and Clo had no impact on lipid lateral mobility, EIS, which is more sensitive to membrane structural change, indicated modest but significant decreases to membrane resistivity consistent with adsorption but weak penetration of drugs at the membrane. Ale and Clo, evaluated at pH 5.25, did not impact the impedance of the membrane except at concentrations exceeding 4 mM. Conversely, Dic and Rif dramatically altered bilayer fluidity, suggesting their translocation through the bilayer, and EIS data showed that resistivity of the membrane decreased substantially with increasing drug concentration. Capacitance changes to the bilayer in most cases were insignificant. Using a Langmuir-Freundlich model to fit the EIS data, we propose R_sat as an empirical value that reflects permeation. Overall, the data indicate that Ibu, Levo, and Pef adsorb at the interface of the lipid membrane but Dic and Rif interact strongly, permeating the membrane core modifying the water/ion permeability of the bilayer structure. These observations are discussed in the context of previously reported data on drug permeability and log P.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Effect of general anesthetics on the properties of lipid membranes of various compositions

Computer simulations of four lipid membranes of different compositions, namely neat DPPC and PSM, and equimolar DPPC-cholesterol and PSM-cholesterol mixtures, are performed in the presence and absence of the general anesthetics diethylether and sevoflurane both at 1 and 600 bar. The results are analyzed in order to identify membrane properties that are potentially related to the molecular mechanism of anesthesia, namely that change in the same way in any membrane with any anesthetics, and change oppositely with increasing pressure. We find that the lateral lipid density satisfies both criteria: it is decreased by anesthetics and increased by pressure. This anesthetic-induced swelling is attributed to only those anesthetic molecules that are located close to the boundary of the apolar phase. This lateral expansion is found to lead to increased lateral mobility of the lipids, an effect often thought to be related to general anesthesia; to an increased fraction of the free volume around the outer preferred position of anesthetics; and to the decrease of the lateral pressure in the nearby range of the ester and amide groups, a region into which anesthetic molecules already cannot penetrate. All these changes are reverted by the increase of pressure. Another important finding of this study is that cholesterol has an opposite effect on the membrane properties than anesthetics, and, correspondingly, these changes are less marked in the presence of cholesterol. Therefore, changes in the membrane that can lead to general anesthesia are expected to occur in the membrane domains of low cholesterol content.

Molecular dynamics simulation study of the effect of halothane on mixed DPPC/DPPE phospholipid membranes

We report results of a molecular dynamics simulation study of the effect of one general anesthetic, halothane, on some properties of mixed DPPC/DPPE phospholipid membranes. This is a suitable model for the study of simple, two-phospholipid membrane systems. From the simulation runs, we determined several membrane properties for five different molecular proportions of DPPC/DPPE. The effect of halothane on the studied membrane properties (area per lipid molecule, density of membrane, order parameter, etc.) was rather small. The distribution of halothane is not uniform through the bilayer thickness. Instead, there is a maximum of anesthetic concentration around 1.2 nm from the center of the membrane. The anesthetic molecule is located close to the phospholipid headgroups. The position of the halothane density maximum depends slightly on the DPPC/DPPE molar proportion. Snapshots taken over the plane of the membrane, as well as calculated two-dimensional radial distribution functions show that the anesthetic has no preference for either phospholipid (DPPC or DPPE). Our results indicate that this anesthetic molecule has only small effects on DPPC/DPPE mixed membranes. In addition, halothane displays no preferential location around DPPC or DPPE. This is probably due to the hydrophobic nature of halothane and to the fact that the chosen phospholipids have the same hydrophobic tails.

Interactions between model cell membranes and the neuroactive drug propofol

Vibrational sum frequency spectroscopy (VSFS) complemented by surface pressure isotherm and neutron reflectometry (NR) experiments were employed to investigate the interactions between propofol, a small amphiphilic molecule that currently is the most common general anaesthetic drug, and phospholipid monolayers. A series of biologically relevant saturated phospholipids of varying chain length from C₁₈ to C₁₄ were spread on either pure water or propofol (2,6-bis(1-methylethyl)phenol) solution in a Langmuir trough, and the change in the molecular structure of the film, induced by the interaction with propofol, was studied with respect to the surface pressure. The results from the surface pressure isotherm experiments revealed that propofol, as long as it remains at the interface, enhances the fluidity of the phospholipid monolayer. The VSF spectra demonstrate that for each phospholipid the amount of propofol in the monolayer region decreases with increasing surface pressure. Such squeeze out is in contrast to the enhanced interactions that can be exhibited by more complex amphiphilic molecules such as peptides. At surface pressures of 22-25 mN m^-1, which are relevant for biological cell membranes, most of the propofol has been expelled from the monolayer, especially in the case of the C₁₆ and C₁₈ phospholipids that adopt a liquid condensed phase packing of its alkyl tails. At lower surface pressures of 5 mN m^-1, the effect of propofol on the structure of the alkyl tails is enhanced when the phospholipids are present in a liquid expanded phase. Specifically, for the C₁₆ phospholipid, NR data reveal that propofol is located exclusively in the head group region, which is rationalized in the context of previous studies. The results imply a non-homogeneous distribution of propofol in the plane of real cell membranes, which is an inference that requires urgent testing and may help to explain why such low concentration of the drug are required to induce general anaesthesia.

Propagation of a thermo-mechanical perturbation on a lipid membrane

The propagation of sound waves on lipid monolayers supported on water has been previously studied during the melting transition. Since changes in volume, area, and compressibility in lipid membranes have biological relevance, the observed sound propagation is of paramount importance. However, it is unknown what would occur on a lipid bilayer, which is a more approximate model of a cell membrane. With the aim to answer this relevant question, we built an experimental setup to assemble long artificial lipid membranes. We found that if these membranes are heated in order to force local melting, a thermo-mechanical perturbation propagates a long distance. Our findings may support the existence of solitary waves, postulated to explain the propagation of isentropic signals together with the action potential in neurons.

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Exploring the Effects on Lipid Bilayer Induced by Noble Gases via Molecular Dynamics Simulations

Noble gases seem to have no significant effect on the anesthetic targets due to their simple, spherical shape. However, xenon has strong narcotic efficacy and can be used clinically, while other noble gases cannot. The mechanism remains unclear. Here, we performed molecular dynamics simulations on phospholipid bilayers with four kinds of noble gases to elucidate the difference of their effects on the membrane. Our results showed that the sequence of effects on membrane exerted by noble gases from weak to strong was Ne, Ar, Kr and Xe, the same order as their relative narcotic potencies as well as their lipid/water partition percentages. Compared with the other three kinds of noble gases, more xenon molecules were distributed between the lipid tails and headgroups, resulting in membrane’s lateral expansion and lipid tail disorder. It may contribute to xenon’s strong anesthetic potency. The results are well consistent with the membrane mediated mechanism of general anesthesia.

Modelling of noble anaesthetic gases and high hydrostatic pressure effects in lipid bilayers

Our objective was to study molecular processes that might be responsible for inert gas narcosis and high-pressure nervous syndrome. The classical molecular dynamics trajectories (200 ns) of dioleoylphosphatidylcholine (DOPC) bilayers simulated by the Berger force field were evaluated for water and the atomic distribution of noble gases around DOPC molecules in the pressure range of 1-1000 bar and at a temperature of 310 K. Xenon and argon have been tested as model gases for general anaesthetics, and neon has been investigated for distortions that are potentially responsible for neurological tremors in hyperbaric conditions. The analysis of stacked radial pair distribution functions of DOPC headgroup atoms revealed the explicit solvation potential of the gas molecules, which correlates with their dimensions. The orientational dynamics of water molecules at the biomolecular interface should be considered as an influential factor, while excessive solvation effects appearing in the lumen of membrane-embedded ion channels could be a possible cause of inert gas narcosis. All the noble gases tested exhibit similar order parameter patterns for both DOPC acyl chains, which are opposite of the patterns found for the order parameter curve at high hydrostatic pressures in intact bilayers. This finding supports the 'critical volume' hypothesis of anaesthesia pressure reversal. The irregular lipid headgroup-water boundary observed in DOPC bilayers saturated with neon in the pressure range of 1-100 bar could be associated with the possible manifestation of neurological tremors at the atomic scale. The non-immobiliser neon also demonstrated the highest momentum impact on the normal component of the DOPC diffusion coefficient representing the monolayer undulation rate, which indicates that enhanced diffusivity rather than atomic size is the key factor.

The effect of anaesthetics on the properties of a lipid membrane in the biologically relevant phase: a computer simulation study

Phospholipid membranes containing four different general anaesthetic molecules are simulated in the biologically relevant L_α phase at atmospheric and high pressures.

Plasma membranes as heat stress sensors: from lipid-controlled molecular switches to therapeutic applications

The classic heat shock (stress) response (HSR) was originally attributed to protein denaturation. However, heat shock protein (Hsp) induction occurs in many circumstances where no protein denaturation is observed. Recently considerable evidence has been accumulated to the favor of the "Membrane Sensor Hypothesis" which predicts that the level of Hsps can be changed as a result of alterations to the plasma membrane. This is especially pertinent to mild heat shock, such as occurs in fever. In this condition the sensitivity of many transient receptor potential (TRP) channels is particularly notable. Small temperature stresses can modulate TRP gating significantly and this is influenced by lipids. In addition, stress hormones often modify plasma membrane structure and function and thus initiate a cascade of events, which may affect HSR. The major transactivator heat shock factor-1 integrates the signals originating from the plasma membrane and orchestrates the expression of individual heat shock genes. We describe how these observations can be tested at the molecular level, for example, with the use of membrane perturbers and through computational calculations. An important fact which now starts to be addressed is that membranes are not homogeneous nor do all cells react identically. Lipidomics and cell profiling are beginning to address the above two points. Finally, we observe that a deregulated HSR is found in a large number of important diseases where more detailed knowledge of the molecular mechanisms involved may offer timely opportunities for clinical interventions and new, innovative drug treatments. This article is part of a Special Issue entitled: Membrane Structure and Function: Relevance in the Cell's Physiology, Pathology and Therapy.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

Biophysical changes induced by xenon on phospholipid bilayers

Structural and dynamic changes in cell membrane properties induced by xenon, a volatile anesthetic molecule, may affect the function of membrane-mediated proteins, providing a hypothesis for the mechanism of general anesthetic action. Here, we use molecular dynamics simulation and differential scanning calorimetry to examine the biophysical and thermodynamic effects of xenon on model lipid membranes. Our results indicate that xenon atoms preferentially localize in the hydrophobic core of the lipid bilayer, inducing substantial increases in the area per lipid and bilayer thickness. Xenon depresses the membrane gel-liquid crystalline phase transition temperature, increasing membrane fluidity and lipid head group spacing, while inducing net local ordering effects in a small region of the lipid carbon tails and modulating the bilayer lateral pressure profile. Our results are consistent with a role for nonspecific, lipid bilayer-mediated mechanisms in producing xenon's general anesthetic action.