Ion transport through chemically induced pores in protein-free phospholipid membranes

Modulation of bacterial membranes and cellular macromolecules by dimethyl sulfoxide: A dose-dependent study providing novel insights

Using Escherichia coli as a model, this manuscript delves into the intricate interactions between dimethyl sulfoxide (DMSO) and membranes, cellular macromolecules, and the effects on various aspects of bacterial physiology. Given DMSO's wide-ranging use as a solvent in microbiology, we investigate the impacts of both non-growth inhibitory (1.0 % and 2.5 % v/v) and slightly growth-inhibitory (5.0 % v/v) concentrations of DMSO. The results demonstrate that DMSO causes alterations in bacterial membrane potential, influences the electrochemical characteristics of the cell surface, and exerts substantial effects on the composition and structure of cellular biomolecules. Genome-wide gene expression data from DMSO-treated E. coli was used to further investigate and bolster the results. The findings of this study provide valuable insights into the complex relationship between DMSO and biological systems, with potential implications in drug delivery and cellular manipulation. However, it is essential to exercise caution when utilizing DMSO to enhance the solubility and delivery of bioactive compounds, as even at low concentrations, DMSO exerts non-inert effects on cellular macromolecules and processes.

A thermodynamic description for physiological transmembrane transport

A general formulation for both passive and active transmembrane transport is derived from basic thermodynamical principles. The derivation takes into account the energy required for the motion of molecules across membranes and includes the possibility of modeling asymmetric flow. Transmembrane currents can then be described by the general model in the case of electrogenic flow. As it is desirable in new models, it is possible to derive other well-known expressions for transmembrane currents as particular cases of the general formulation. For instance, the conductance-based formulation for current turns out to be a linear approximation of the general formula for current. Also, under suitable assumptions, other formulas for current based on electrodiffusion, like the constant field approximation by Goldman, can be recovered from the general formulation. The applicability of the general formulations is illustrated first with fits to existing data, and after, with models of transmembrane potential dynamics for pacemaking cardiocytes and neurons. The general formulations presented here provide a common ground for the biophysical study of physiological phenomena that depend on transmembrane transport.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Molecular mechanism of the synergistic effects of vitrification solutions on the stability of phospholipid bilayers

The vitrification solutions used in the cryopreservation of biological samples aim to minimize the deleterious formation of ice by dehydrating cells and promoting the formation of the glassy state of water. They contain a mixture of different cryoprotective agents (CPAs) in water, typically polyhydroxylated alcohols and/or dimethyl sulfoxide (DMSO), which can damage cell membranes. Molecular dynamics simulations have been used to investigate the behavior of pure DPPC, pure DOPC, and mixed DOPC-β-sitosterol bilayers solvated in a vitrification solution containing glycerol, ethylene glycol, and DMSO at concentrations that approximate the widely used plant vitrification solution 2. As in the case of solutions containing a single CPA, the vitrification solution causes the bilayer to thin and become disordered, and pores form in the case of some bilayers. Importantly, the degree of thinning is, however, substantially reduced compared to solutions of DMSO containing the same total CPA concentration. The reduction in the damage done to the bilayers is a result of the ability of the polyhydroxylated species (especially glycerol) to form hydrogen bonds to the lipid and sterol molecules of the bilayer. A decrease in the amount of DMSO in the vitrification solution with a corresponding increase in the amount of glycerol or ethylene glycol diminishes further its damaging effect due to increased hydrogen bonding of the polyol species to the bilayer headgroups. These findings rationalize, to our knowledge for the first time, the synergistic effects of combining different CPAs, and form the basis for the optimization of vitrification solutions.

Plasma membrane and actin cytoskeleton as synergistic barriers to nanowire cell penetration

Nanowires are a rapidly emerging platform for manipulation of and material delivery directly into the cell cytosol. These high aspect ratio structures can breach the lipid membrane; however, the yield of penetrant structures is low, and the mechanism is largely unknown. In particular, some nanostructures appear to defeat the membrane transiently, while others can retain long-term access. Here, we examine if local dissolution of the lipid membrane, actin cytoskeleton, or both can enhance nanowire penetration. It is possible that, during cell contact, membrane rupture occurs; however, if the nanostructures do not penetrate the cytoskeleton, the membrane may reclose over a relatively short time frame. We show with quantitative analysis of the number of penetrating nanowires that the lipid bilayer and actin cytoskeleton are synergistic barriers to nanowire cell access, yet chemical poration through both is still insufficient to increase long-term access for adhered cells.

A permeability barrier surrounds taste buds in lingual epithelia

Epithelial tissues are characterized by specialized cell-cell junctions, typically localized to the apical regions of cells. These junctions are formed by interacting membrane proteins and by cytoskeletal and extracellular matrix components. Within the lingual epithelium, tight junctions join the apical tips of the gustatory sensory cells in taste buds. These junctions constitute a selective barrier that limits penetration of chemosensory stimuli into taste buds (Michlig et al. J Comp Neurol 502: 1003-1011, 2007). We tested the ability of chemical compounds to permeate into sensory end organs in the lingual epithelium. Our findings reveal a robust barrier that surrounds the entire body of taste buds, not limited to the apical tight junctions. This barrier prevents penetration of many, but not all, compounds, whether they are applied topically, injected into the parenchyma of the tongue, or circulating in the blood supply, into taste buds. Enzymatic treatments indicate that this barrier likely includes glycosaminoglycans, as it was disrupted by chondroitinase but, less effectively, by proteases. The barrier surrounding taste buds could also be disrupted by brief treatment of lingual tissue samples with DMSO. Brief exposure of lingual slices to DMSO did not affect the ability of taste buds within the slice to respond to chemical stimulation. The existence of a highly impermeable barrier surrounding taste buds and methods to break through this barrier may be relevant to basic research and to clinical treatments of taste.

Atomistic simulations of pore formation and closure in lipid bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

Luminal DMSO: effects on detrusor and urothelial/lamina propria function

DMSO is used as a treatment for interstitial cystitis and this study examined the effects of luminal DMSO treatment on bladder function and histology. Porcine bladder was incubated without (controls) or with DMSO (50%) applied to the luminal surface and the release of ATP, acetylcholine, and LDH assessed during incubation and in tissues strips after DMSO incubation. Luminally applied DMSO caused ATP, Ach, and LDH release from the urothelial surface during treatment, with loss of urothelial layers also evident histologically. In strips of urothelium/lamina propria from DMSO pretreated bladders the release of both ATP and Ach was depressed, while contractile responses to carbachol were enhanced. Detrusor muscle contractile responses to carbachol were not affected by DMSO pretreatment, but neurogenic responses to electrical field stimulation were enhanced. The presence of an intact urothelium/lamina propria inhibited detrusor contraction to carbachol by 53% and this inhibition was significantly reduced in DMSO pretreated tissues. Detection of LDH in the treatment medium suggests that DMSO permeabilised urothelial membranes causing leakage of cytosolic contents including ATP and Ach rather than enhancing release of these mediators. The increase in contractile response and high levels of ATP are consistent with initial flare up in IC/PBS symptoms after DMSO treatment.

The free energy of nanopores in tense membranes

Membrane nanopores are central players for a range of important cellular membrane remodeling processes as well as membrane rupture. Understanding pore formation in tense membranes requires comprehension of the molecular mechanism of pore formation and the associated free energy change as a function of the membrane tension. Here we propose a scheme to calculate the free energy change associated with the formation of a nanometer sized pore in molecular dynamics simulations as a function of membrane tension, which requires the calculation of only one computationally expensive potential of mean force. We show that membrane elastic theory can be used to estimate the pore formation free energy at different tension values from the free energy change in a relaxed membrane and the area expansion curves of the membranes. We have computed the pore formation free energy for a dipalmitoyl-phosphatidylcholine (DPPC) membrane at two different lateral pressure values, 1 bar and -40 bar, by calculating the potential of mean force acting on the head group of a single lipid molecule. Unrestrained simulations of the closing process confirm that the intermediate states along this reaction coordinate are reasonable and show that hydrophilic indentations spanning half the bilayer connected by a hydrophobic pore segment represent the corresponding high energy transition state. A comparison of the stability of simulated membranes to experiment at high loading rates show that, contrary to expectation, pores form too easily in small simulated membrane patches. This discrepancy originates from a combination of the absence of ions in the simulations and the small membrane size.

Ion-induced defect permeation of lipid membranes

We have explored the mechanisms of uncatalyzed membrane ion permeation using atomistic simulations and electrophysiological recordings. The solubility-diffusion mechanism of membrane charge transport has prevailed since the 1960s, despite inconsistencies in experimental observations and its lack of consideration for the flexible response of lipid bilayers. We show that direct lipid bilayer translocation of alkali metal cations, Cl(-), and a charged arginine side chain analog occurs via an ion-induced defect mechanism. Contrary to some previous suggestions, the arginine analog experiences a large free-energy barrier, very similar to those for Na(+), K(+), and Cl(-). Our simulations reveal that membrane perturbations, due to the movement of an ion, are central for explaining the permeation process, leading to both free-energy and diffusion-coefficient profiles that show little dependence on ion chemistry and charge, despite wide-ranging hydration energies and the membrane's dipole potential. The results yield membrane permeabilities that are in semiquantitative agreement with experiments in terms of both magnitude and selectivity. We conclude that ion-induced defect-mediated permeation may compete with transient pores as the dominant mechanism of uncatalyzed ion permeation, providing new understanding for the actions of a range of membrane-active peptides and proteins.

Molecular dynamics simulations of the interactions of DMSO, mono- and polyhydroxylated cryosolvents with a hydrated phospholipid bilayer

Molecular dynamics (MD) simulations have been used to investigate the interactions of a variety of hydroxylated cryosolvents (glycerol, propylene glycol and ethylene glycol), methanol and dimethyl sulfoxide (DMSO) in aqueous solution with a 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) bilayer in its fluid phase at 323K. Each cryosolvent induced lateral expansion of the membrane leading to thinning of the bilayer and resulting in disordering of the lipid hydrocarbon chains. Propylene glycol and DMSO were observed to exhibit a greater disordering effect on the structure of the membrane than the other three alcohols. Closer examination exposed a number of effects on the lipid bilayer as a function of the molecular size and hydrogen bonding capacity of the cryosolvents. Analyses of hydrogen bonds revealed that increased concentrations of the polyhydroxylated cryosolvents induced the formation of a cross-linked cryosolvent layer across the surface of the membrane bilayer. This effect was most pronounced for glycerol at sufficiently high concentrations, which displayed a comparatively enhanced capacity to induce cross-linking of lipid headgroups resulting in the formation of extensive hydrogen bonding bridges and the promotion of a dense cryosolvent layer across the phospholipid bilayer.

Lipid nanotechnology

Nanotechnology is a multidisciplinary field that covers a vast and diverse array of devices and machines derived from engineering, physics, materials science, chemistry and biology. These devices have found applications in biomedical sciences, such as targeted drug delivery, bio-imaging, sensing and diagnosis of pathologies at early stages. In these applications, nano-devices typically interface with the plasma membrane of cells. On the other hand, naturally occurring nanostructures in biology have been a source of inspiration for new nanotechnological designs and hybrid nanostructures made of biological and non-biological, organic and inorganic building blocks. Lipids, with their amphiphilicity, diversity of head and tail chemistry, and antifouling properties that block nonspecific binding to lipid-coated surfaces, provide a powerful toolbox for nanotechnology. This review discusses the progress in the emerging field of lipid nanotechnology.

Breaching the skin barrier--insights from molecular simulation of model membranes

Breaching the skin's barrier function by design is an important strategy for delivering drugs and vaccines to the body. However, while there are many proposed approaches for reversibly breaching the skin barrier, our understanding of the molecular processes involved is still rudimentary. Molecular simulation offers an unprecedented molecular-level resolution with an ability to reproduce molecular and bulk level properties. We review the basis of the molecular simulation methodology and give applications of relevance to the skin lipid barrier, focusing on permeation of molecules and chemical approaches for breaching the lipid barrier by design. The bulk kinetic model based on Fick's Law describing absorption of a drug through skin has been reconciled with statistical mechanical quantities such as the local excess chemical potential and local diffusion coefficient within the membrane structure. Applications of molecular simulation reviewed include investigations of the structure and dynamics of simple models of skin lipids, calculation of the permeability of molecules in simple model membranes, and mechanisms of action of the penetration enhancers, DMSO, ethanol and oleic acid. The studies reviewed illustrate the power and potential of molecular simulation to yield important physical insights, inform and rationalize experimental studies, and to predict structural changes, and kinetic and thermodynamic quantities.

Effects of dimethyl sulfoxide in cholesterol-containing lipid membranes: a comparative study of experiments in silico and with cells

Dimethyl sulfoxide (DMSO) has been known to enhance cell membrane permeability of drugs or DNA. Molecular dynamics (MD) simulations with single-component lipid bilayers predicted the existence of three regimes of action of DMSO: membrane loosening, pore formation and bilayer collapse. We show here that these modes of action are also reproduced in the presence of cholesterol in the bilayer, and we provide a description at the atomic detail of the DMSO-mediated process of pore formation in cholesterol-containing lipid membranes. We also successfully explore the applicability of DMSO to promote plasma membrane permeability to water, calcium ions (Ca(2+)) and Yo-Pro-1 iodide (Yo-Pro-1) in living cell membranes. The experimental results on cells in culture can be easily explained according to the three expected regimes: in the presence of low doses of DMSO, the membrane of the cells exhibits undulations but no permeability increase can be detected, while at intermediate DMSO concentrations cells are permeabilized to water and calcium but not to larger molecules as Yo-Pro-1. These two behaviors can be associated to the MD-predicted consequences of the effects of the DMSO at low and intermediate DMSO concentrations. At larger DMSO concentrations, permeabilization is larger, as even Yo-Pro-1 can enter the cells as predicted by the DMSO-induced membrane-destructuring effects described in the MD simulations.

Ion transport through dimethyl sulfoxide (DMSO) induced transient water pores in cell membranes

It is well known that dimethyl sulphoxide (DMSO) increases membrane permeability, which makes it widely used as a vehicle to facilitate drug delivery across biological membranes. However, the mechanism of how DMSO increases membrane permeability has not been well understood. Recently, molecular dynamics simulations have demonstrated that DMSO can induce water pores in biological membranes, but no direct experimental evidence is so far available to prove the simulation result. Using FluxOR Tl⁺ influx assay and intracellular Ca²⁺ imaging technique, we studied the effect of DMSO on Tl⁺ and Ca²⁺ permeation across cell membranes. Upon application of DMSO on CHO-K1 cell line, Tl⁺ influx was transiently increased in a dose-dependent manner. The increase in Tl⁺ permeability induced by DMSO was not changed in the presence of blockers for K⁺ channel and Na⁺-K⁺ ATPase, suggesting that Tl⁺ permeates through transient water pores induced by DMSO to enter into the cell. In addition, Ca²⁺ permeability was significantly increased upon application of DMSO, indicating that the transient water pores induced by DMSO were non-selective pores. Furthermore, similar results could be obtained from RAW264.7 macrophage cell line. Therefore, this study provided experimental evidence to support the prediction that DMSO can induce transient water pores in cell membranes, which in turn facilitates the transport of active substances across membranes.

Interactions of a hydrophobically modified polycation with zwitterionic lipid membranes

The interactions between synthetic polycations and phospholipid bilayers play an important role in some biophysical applications such as gene delivery or antibacterial usage. Despite extensive investigation into the nature of these interactions, their physical and molecular bases remain poorly understood. In this Article, we present the results of our studies on the impact of a hydrophobically modified strong polycation on the properties of a zwitterionic bilayer used as a model of the mammalian cellular membrane. The study was carried out using a set of complementary experimental methods and molecular dynamic (MD) simulations. A new polycation, poly(allyl-N,N-dimethyl-N-hexylammonium chloride) (polymer 3), was synthesized, and its interactions with liposomes composed of 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) were examined using dynamic light scattering (DLS), zeta potential measurements, and cryo-transmission electron microscopy (cryo-TEM). Our results have shown that polymer 3 can efficiently associate with and insert into the POPC membrane. However, it does not change its lamellar structure, as was demonstrated by cryo-TEM. The influence of polymer 3 on the membrane functionality was studied by leakage experiments applying a fluorescence dye (calcein) encapsulated in the phospholipid vesicles. The MD simulations of model systems reveal that polymer 3 promotes formation of hydrophilic pores in the membrane, thus increasing considerably its permeability.

Electrochemical impedance spectroscopy of tethered bilayer membranes

The electrochemical impedance spectra (EIS) of tethered bilayer membranes (tBLMs) were analyzed, and the analytical solution for the spectral response of membranes containing natural or artificially introduced defects was derived. The analysis carried out in this work shows that the EIS features of an individual membrane defect cannot be modeled by conventional electrical elements. The primary reason for this is the complex nature of impedance of the submembrane ionic reservoir separating the phospholipid layer and the solid support. We demonstrate that its EIS response, in the case of radially symmetric defects, is described by the Hankel functions of a complex variable. Therefore, neither the impedance of the submembrane reservoir nor the total impedance of tBLMs can be modeled using the conventional elements of the equivalent electrical circuits of interfaces. There are, however, some limiting cases in which the complexity of the EIS response of the submembrane space reduces. In the high frequency limit, the EIS response of a submembrane space that surrounds the defect transforms into a response of a constant phase element (CPE) with the exponent (α) value of 0.5. The onset of this transformation is, beside other parameters, dependent on the defect size. Large-sized defects push the frequency limit lower, therefore, the EIS spectra exhibiting CPE behavior with α ≈ 0.5, can serve as a diagnostic criterion for the presence of such defects. In the low frequency limit, the response is dependent on the density of the defects, and it transforms into the capacitive impedance if the area occupied by a defect is finite. The higher the defect density, the higher the frequency edge at which the onset of the capacitive behavior is observed. Consequently, the presented analysis provides practical tools to evaluate the defect density in tBLMs, which could be utilized in tBLM-based biosensor applications. Alternatively, if the parameters of the defects, e.g., ion channels, such as the diameter and the conductance are known, the EIS data analysis provides a possibility to estimate other physical parameters of the system, such as thickness of the submembrane reservoir and its conductance. Finally, current analysis demonstrates a possibility to discriminate between the situations, in which the membrane defects are evenly distributed or clustered on the surface of tBLMs. Such sensitivity of EIS could be used for elucidation of the mechanisms of interaction between the proteins and the membranes.

Nanoscopic description of biomembrane electrostatics: results of molecular dynamics simulations and fluorescence probing

Electrostatic fields generated on and inside biological membranes are recognized to play a fundamental role in key processes of cell functioning. Their understanding requires an adequate description on the level of elementary charges and the reconstruction of electrostatic potentials by integration over all elementary interactions. Out of all the available research tools, only molecular dynamics simulations are capable of this, extending from the atomic to the mesoscopic level of description on the required time and space scale. A complementary approach is that offered by molecular probe methods, with the application of electrochromic dyes. Highly sensitive to intermolecular interactions, they generate integrated signals arising from electric fields produced by elementary charges at the sites of their location. This review is an attempt to provide a critical analysis of these two approaches and their present and potential applications. The results obtained by both methods are consistent in that they both show an extremely complex profile of the electric field in the membrane. The nanoscopic view, with two-dimensional averaging over the bilayer plane and formal separation of the electrostatic potential into surface (Psi(s)), dipole (Psi(d)) and transmembrane (Psi(t)) potentials, is constructive in the analysis of different functional properties of membranes.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.

Chemically induced phospholipid translocation across biological membranes

Chemical means of manipulating the distribution of lipids across biological membranes is of considerable interest for many biomedical applications as a characteristic lipid distribution is vital for numerous cellular functions. Here we employ atomic-scale molecular simulations to shed light on the ability of certain amphiphilic compounds to promote lipid translocation (flip-flops) across membranes. We show that chemically induced lipid flip-flops are most likely pore-mediated: the actual flip-flop event is a very fast process (time scales of tens of nanoseconds) once a transient water defect has been induced by the amphiphilic chemical (dimethylsulfoxide in this instance). Our findings are consistent with available experimental observations and further emphasize the importance of transient membrane defects for chemical control of lipid distribution across cell membranes.