Counter-effects of Ethanol and Cholesterol on the Heterogeneous PSM–POPC Lipid Membrane: A Molecular Dynamics Simulation Study

Salting-out and Competitive Adsorption of Ethanol into Lipid Bilayer Membranes: Conflicting Effects of Salts on Ethanol-Membrane Interactions Studied by Molecular Dynamics Simulations

Small amphiphilic molecules, such as ethanol, disturb the structure of lipid bilayer membranes to increase the membrane permeability, which is important for applications such as drug delivery, disinfection, and fermentation. To investigate how and the extent to which coexisting salts affect membrane disturbance, we performed molecular dynamics (MD) simulations on lipid bilayer membranes composed of zwitterionic lipids in aqueous ethanol solutions containing 0-631 mM NaCl, KCl, and KI salts. The addition of salts at low concentrations induced cationic adsorption on the lipid membrane, which competes with ethanol adsorption, thereby reducing the hydrogen bonds between ethanol and lipid molecules. This competitive adsorption mitigated the membrane disturbance and decreased the permeation of ethanol molecules into the membrane. In contrast, higher salt concentrations enhanced the membrane disturbance and permeability, which was caused by the salting-out of ethanol from the aqueous phase to the lipid bilayer. These conflicting effects appearing at different concentrations were stronger with the chloride salts than with the iodide salt. Among the two chloride salts, NaCl and KCl, the latter showed a greater enhancement in ethanol permeation at high concentrations. This seeming anti-Hofmeister salting-out behavior resulted from greater Na+ adsorption, preventing the ethanol-lipid interactions.

Influence of cholesterol on hydrogen-bond dynamics of water molecules in lipid-bilayer systems at varying temperatures

Cholesterol (Chol) plays a crucial role in shaping the intricate physicochemical attributes of biomembranes, exerting a considerable influence on water molecules proximal to the membrane interface. In this study, we conducted molecular dynamics simulations on the bilayers of two lipid species, dipalmitoylphosphatidylcholine (DPPC) and palmitoyl sphingomyelin; they are distinct with respect to the structures of the hydrogen-bond (H-bond) acceptors. Our investigation focuses on the dynamic properties and H-bonds of water molecules in the lipid-membrane systems, with a particular emphasis on the influence of Chol at varying temperatures. Notably, in the gel phase at 303 K, the presence of Chol extends the lifetimes of H-bonds of the oxygen atoms acting as H-bond acceptors within DPPC with water molecules by a factor of 1.5-2.5. In the liquid-crystalline phase at 323 K, on the other hand, H-bonding dynamics with lipid membranes remain largely unaffected by Chol. This observed shift in H-bonding states serves as a crucial key to unraveling the subtle control mechanisms governing water dynamics in lipid-membrane systems.

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

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

Wrapping-Trapping versus Extraction Mechanism of Bactericidal Activity of MoS2 Nanosheets against Staphylococcus aureus Bacterial Membrane

The promising broad-spectrum antibacterial activity of two-dimensional molybdenum disulfide (2D MoS₂) has been widely recognized in the past decade. However, a comprehensive understanding of how the antibacterial pathways opted by the MoS₂ nanosheets varies with change in lipid compositions of different bacterial strains is imperative to harness their full antibacterial potential and remains unexplored thus far. Herein, we present an atomistic molecular dynamics (MD) study to investigate the distinct modes of antibacterial action of MoS₂ nanosheets against Staphylococcus aureus (S. aureus) under varying conditions. We observed that the freely dispersed nanosheets readily adhered to the bacterial membrane outer surface and opted for an unconventional surface directed "wrapping-trapping" mechanism at physiological temperature (i.e., 310 K). The adsorbed nanosheets mildly influenced the membrane structure by originating a compact packing of the lipid molecules present in its direct contact. Interestingly, these surface adsorbed nanosheets exhibited extensive phospholipid extraction to their surface, thereby inducing transmembrane water passage analogous to the cellular leakage, even at a slight increment of 20 K in the temperature. The strong van der Waals interactions between lipid fatty acyl tails and MoS₂ basal planes were primarily responsible for this destructive phospholipid extraction. In addition, the MoS₂ nanosheets bound to an imaginary substrate, controlling their vertical alignment, demonstrated a "nano-knives" action by spontaneously piercing inside the membrane core through their sharp corner, subsequently causing localized lipid ordering in their vicinity. The larger nanosheet produced a more profound deteriorating impact in all of the observed mechanisms. Keeping the existing knowledge about the bactericidal activity of 2D MoS₂ in view, our study concludes that their antibacterial activity is strongly governed by the lipid composition of the bacterial membrane and can be intensified either by controlling the nanosheet vertical alignment or by moderately warming up the systems.

Anionic Lipid Clustering-Mediated Bactericidal Activity and Selective Toxicity of Quaternary Ammonium-Substituted Polycationic Pullulan against the Staphylococcus aureus Bacterial Membrane

Non-amphiphilic polycations have recently been recognized to hold excellent antimicrobial potential with great mammalian cell compatibility. In a recent study, the excellent broad-spectrum bactericidal efficacy of a quaternary ammonium-substituted cationic pullulan (CP4) was demonstrated. Their selective toxicity and nominal probability to induce the acquisition of resistance among pathogens fulfill the fundamental requirements of new-generation antibacterials. However, there have been exiguous attempts in the literature to understand the antimicrobial activity of polycations against Gram-positive bacterial membranes. Here, for the first time, we have scrutinized the molecular level interactions of CP4 tetramers with a model Staphylococcus aureus membrane to understand their probable antibacterial function using molecular dynamics simulations. Our analysis reveals that the hydrophilic CP4 molecules are spontaneously adsorbed onto the membrane outer leaflet surface by virtue of strong electrostatic interactions and do not penetrate into the lipid tail hydrophobic region. This surface binding of CP4 is strengthened by the formation of anionic lipid-rich domains in their vicinity, causing lateral compositional heterogeneity. The major outcomes of the asymmetric accumulation of bulky polycationic CP4 on one leaflet are (i) anionic lipid segregation at the interaction site and (ii) a decrease in the cationic lipid acyl tail ordering and ease of water translocation across the lipid hydrophobic barrier. The membrane-CP4 interactions are strongly monitored by the ionic strength; a higher salt concentration weakens the binding of CP4 on the membrane surface. In addition, our study also substantiates the non-interacting behavior of CP4 oligomers with biomimetic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane, indicating their cell selectivity and specificity against pathogenic membranes.

Removal of tetracycline by aerobic granular sludge from marine aquaculture wastewater: A molecular dynamics investigation

Although biological treatment of marine aquaculture wastewater is promising, the fundamental principles driving the adsorption of tetracycline to microbial cell membrane are not well understood. Using a combination of experiments and molecular dynamics (MD) simulations, the mechanism underlying the biological removal of tetracycline from seawater was investigated. More than 90% tetracycline removal was achieved in an aerobic granular sludge system, with degradation accounting for 30% of total removal. A model of the tetracycline-dipalmitoylphosphatidylcholine lipid bilayers was established to elucidate the transport mechanism of tetracycline from bulk solution to microorganisms' cell membrane. 62% of the driving force for tetracycline adsorption on the cell membrane originated from electrostatic attraction. The electrophilic groups on tetracycline (amino and aromatic groups) were attracted to the phosphate groups in the cell membrane. Sodium ions, which are abundant in seawater, decreased the interaction energy between tetracycline and the cell membrane.

Impact of Surface Polarity on Lipid Assembly under Spatial Confinement

Molecular dynamics (MD) simulations in the MARTINI model are used to study the assembly of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) molecules under spatial confinement, such as during solvent evaporation from ultrasmall (femtoliter quantity) droplets. The impact of surface polarity on molecular assembly is discussed in detail. To the best of our knowledge, this work represents the first of its kind. Our results reveal that solvent evaporation gives rise to the formation of well-defined stacks of lipid bilayers in a smectic alignment. These smectic mesophases form on both polar and nonpolar surfaces but with a notable distinction. On polar surfaces, the director of the stack is oriented perpendicular to the support surface. By contrast, the stacks orient at an angle on the nonpolar surfaces. The packing of head groups on surfaces and lipid molecular mobility exhibits significant differences as surface polarity changes. The role of glycerol in the assembly and stability is also revealed. The insights revealed from the simulation have a significant impact on additive manufacturing, biomaterials, model membranes, and engineering protocells. For example, POPC assemblies via evaporation of ultrasmall droplets were produced and characterized. The trends compare well with the bilayer stack models. The surface polarity influences the local morphology and structures at the interfaces, which could be rationalized via the molecule-surface interactions observed from simulations.

Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations

Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.

Mechanistic Insight on BioIL-Induced Structural Alterations in DMPC Lipid Bilayer

The emerging application risks of traditional ionic liquids (ILs) toward the ecosystem have changed the perception regarding their greenness. This resulted in the exploration of their more biocompatible alternatives known as biocompatible ILs (BioILs). Here, we have investigated the impact of two such biocompatible cholinium amino acid-based ILs on the structural behavior of model homogeneous DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) lipid bilayer using all-atom molecular dynamics simulation technique. Two classic cholinium-amino acid-based ILs, cholinium glycinate ([Ch][Gly]) and cholinium phenylalaninate ([Ch][Phe]), which differ only by the side chain lengths and hydrophobicity of the anions, have been utilized in the present work. Simultaneous analysis of the bilayer structural properties reveals that the existence of [Ch][Gly] BioIL above a particular concentration induces phase transition from fluid phase to gel phase in the DMPC lipid bilayer. Such a freezing of lipid bilayer upon the exposure to concentrated aqueous solution of [Ch][Gly] BioIL indicates the harmfulness of this BioIL toward the cell membranes majorly containing DMPC lipids, as the cell freezing can negatively affect its stability and functionality. Despite having a more hydrophobic amino acid side chain of [Phe]^- anion in [Ch][Phe], in the case of bilayer-[Ch][Phe] systems we observe the minimal impact of [Ch][Phe] BioIL on the DMPC bilayer properties up to 10 mol % concentration. In the presence of these BioIL, we observe the thickening of the bilayer and accumulation of the cations and anions of the BioILs at the interface of DMPC lipid heads and tails. The transfer free-energy profile of a [Phe]^- anion from aqueous phase to membrane center also indicates the anion partitioning at lipid head-tail interface and its inability to penetrate in the lipid membrane tail region. In contrast, the free-energy profile for a [Gly]^- anion offers a very high energy barrier to the insertion of [Gly]^- into the membrane interior, leading to accumulation of [Gly]^- anions at the lipid head-water region.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

The role of the envelope protein in the stability of a coronavirus model membrane against an ethanolic disinfectant

Ethanol is highly effective against various enveloped viruses and can disable the virus by disintegrating the protective envelope surrounding it. The interactions between the coronavirus envelope (E) protein and its membrane environment play key roles in the stability and function of the viral envelope. By using molecular dynamics simulation, we explore the underlying mechanism of ethanol-induced disruption of a model coronavirus membrane and, in detail, interactions of the E-protein and lipids. We model the membrane bilayer as N-palmitoyl-sphingomyelin and 1-palmitoyl-2-oleoylphosphatidylcholine lipids and the coronavirus E-protein. The study reveals that ethanol causes an increase in the lateral area of the bilayer along with thinning of the bilayer membrane and orientational disordering of lipid tails. Ethanol resides at the head-tail region of the membrane and enhances bilayer permeability. We found an envelope-protein-mediated increase in the ordering of lipid tails. Our simulations also provide important insights into the orientation of the envelope protein in a model membrane environment. At ∼25 mol. % of ethanol in the surrounding ethanol-water phase, we observe disintegration of the lipid bilayer and dislocation of the E-protein from the membrane environment.

How Alcoholic Disinfectants Affect Coronavirus Model Membranes: Membrane Fluidity, Permeability, and Disintegration

Atomistic molecular dynamics simulations have been carried out with a view to investigating the stability of the SARS-CoV-2 exterior membrane with respect to two common disinfectants, namely, aqueous solutions of ethanol and n-propanol. We used dipalmitoylphosphatidylcholine (DPPC) as a model membrane material and did simulations on both gel and liquid crystalline phases of membrane surrounded by aqueous solutions of varying alcohol concentrations (up to 17.5 mol %). While a moderate effect of alcohol on the gel phase of membrane is observed, its liquid crystalline phase is shown to be influenced dramatically by either alcohol. Our results show that aqueous solutions of only 5 and 10 mol % alcohol already have significant weakening effects on the membrane. The effects of n-propanol are always stronger than those of ethanol. The membrane changes its structure, when exposed to disinfectant solutions; uptake of alcohol causes it to swell laterally but to shrink vertically. At the same time, the orientational order of lipid tails decreases significantly. Metadynamics and grand-canonical ensemble simulations were done to calculate the free-energy profiles for permeation of alcohol and alcohol/water solubility in the DPPC. We found that the free-energy barrier to permeation of the DPPC liquid crystalline phase by all permeants is significantly lowered by alcohol uptake. At a disinfectant concentration of 10 mol %, it becomes insignificant enough to allow almost free passage of the disinfectant to the inside of the virus to cause damage there. It should be noted that the disinfectant also causes the barrier for water permeation to drop. Furthermore, the shrinking of the membrane thickness shortens the gap needed to be crossed by penetrants from outside the virus into its core. The lateral swelling also increases the average distance between head groups, which is a secondary barrier to membrane penetration, and hence further increases the penetration by disinfectants. At alcohol concentrations in the disinfectant solution above 15 mol %, we reliably observe disintegration of the DPPC membrane in its liquid crystalline phase.

Molecular Dynamics Evaluation of the Effect of Cholinium Phenylalaninate Biocompatible Ionic Liquid on Biomimetic Membranes

Toward the search of sustainable green solvents, choline amino acid ([Ch][AA]) ionic liquids (ILs), mainly derived from renewable feedstocks, have emerged as a promising atoxic alternative to the conventional solvents. Recent studies have shown the remarkably benign nature of cholinium-based ILs against biomimetic phospholipid membranes. However, few of the contemporaneous experimental studies have contradicted the aforesaid ecofriendly nature of these ILs with anions comprising longer alkyl or aromatic tails. Aiming to understand the influence of amino acid side-chain variation in a particular bio-IL on biomembranes, herein, we have evaluated the effect of cholinium phenylalaninate ([Ch][Phe]) IL on the structural stability of homogeneous biomimetic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) lipid bilayers using atomistic molecular dynamics simulations. Although we observe spontaneous intercalation of aromatic rings of [Phe]^- anions in the hydrophobic region of the bilayers, the polar backbone of the anion remains coordinated with the lipid polar part through strong electrostatic and H-bonding interactions. Besides, the [Ch]⁺ cations get accumulated at the lipid-water interface to counter the excess negative charge density. The intercalation of the anionic rings causes significant perturbations in the lipid structural arrangement while still maintaining the bilayer integrity. The quantitative evaluation to probe the deteriorating effect of this bio-IL application establishes anions as the principal component causing the observed structural perturbations. The analysis of the structural properties along with the free energy assessment reveals the higher efficacy of [Ch][Phe] bio-IL to perturb the POPE bilayer structure than the POPC bilayer.