Molecular Dynamics Simulations of Lipid Membrane Electroporation

The permeability of cell membranes can be transiently increased following the application of external electric fields. Theoretical approaches such as molecular modeling provide a significant insight into the processes affecting, at the molecular level, the integrity of lipid cell membranes when these are subject to voltage gradients under similar conditions as those used in experiments. This article reports on the progress made so far using such simulations to model membrane—lipid bilayer—electroporation. We first describe the methods devised to perform in silico experiments of membranes subject to nanosecond, megavolt-per-meter pulsed electric fields and of membranes subject to charge imbalance, mimicking therefore the application of low-voltage, long-duration pulses. We show then that, at the molecular level, the two types of pulses produce similar effects: provided the TM voltage these pulses create are higher than a certain threshold, hydrophilic pores stabilized by the membrane lipid headgroups form within the nanosecond time scale across the lipid core. Similarly, when the pulses are switched off, the pores collapse (close) within similar time scales. It is shown that for similar TM voltages applied, both methods induce similar electric field distributions within the membrane core. The cascade of events following the application of the pulses, and taking place at the membrane, is a direct consequence of such an electric field distribution.     

Quantitative study of molecular transport due to electroporation: uptake of bovine serum albumin by erythrocyte ghosts.

Electroporation is believed to involve the creation of aqueous pathways in lipid bilayer membranes by transient elevation of the transmembrane voltage to approximately 1 V. Here, results are presented for a quantitative study of the number of bovine serum albumin (BSA) molecules transported into erythrocyte ghosts caused by electroportion. 1) Uptake of BSA was found to plateau at high field strength. However, this was not necessarily an absolute maximum in transport. Instead, it represented the maximum effect of increasing field strength for a particular pulse protocol. 2) Maximum uptake under any conditions used in this study corresponded to approximately one-fourth of apparent equilibrium with the external solution. 3) Multiple and longer pulses each increased uptake of BSA, where the total time integral of field strength correlated with uptake, independent of inter-pulse spacing. 4) Pre-pulse adsorption of BSA to ghost membranes appears to have increased transport. 5) Most transport of BSA probably occurred by electrically driven transport during pulses; post-pulse uptake occurred, but to a much lesser extent. Finally, approaches to increasing transport are discussed.     

Electroporation in dense cell suspension--theoretical and experimental analysis of ion diffusion and cell permeabilization

Electroporation is a process where increased permeability of cells exposed to an electric field is observed. It is used in many biomedical applications including electrogene transfection and electrochemotherapy. Although the increased permeability of the membrane is believed to be the result of pores due to an induced transmembrane voltage U(m), the exact molecular mechanisms are not fully explained. In this study we analyze transient conductivity changes during the electric pulses and increased membrane permeability for ions and molecules after the pulses in order to determine which parameters affect stabilization of pores, and to analyze the relation between transient pores and long-lived transport pores. By quantifying ion diffusion, fraction of transport pores f(per) was obtained. A simple model, which assumes a quadratic dependence of f(per) on E in the area where U(m)>U(c) very accurately describes experimental values, suggesting that f(per) increases with higher electric field due to larger permeabilized area and due to higher energy available for pore formation. The fraction of transport pores increases also with the number of pulses N, which suggest that each pulse contributes to formation of more and/or larger stable transport pores, whereas the number of transient pores does not depend on N.     

Life Cycle of an Electropore: Field-Dependent and Field-Independent Steps in Pore Creation and Annihilation

Electropermeabilization, an electric field-induced modification of the barrier functions of the cell membrane, is widely used in laboratories and increasingly in the clinic; but the mechanisms and physical structures associated with the electromanipulation of membrane permeability have not been definitively characterized. Indirect experimental observations of electrical conductance and small molecule transport as well as molecular dynamics simulations have led to models in which hydrophilic pores form in phospholipid bilayers with increased probability in the presence of an electric field. Presently available methods do not permit the direct, nanoscale examination of electroporated membranes that would confirm the existence of these structures. To facilitate the reconciliation of poration models with the observed properties of electropermeabilized lipid bilayers and cell membranes, we propose a scheme for characterizing the stages of electropore formation and resealing. This electropore life cycle, based on molecular dynamics simulations of phospholipid bilayers, defines a sequence of discrete steps in the electric field-driven restructuring of the membrane that leads to the formation of a head group-lined, aqueous pore and then, after the field is removed, to the dismantling of the pore and reassembly of the intact bilayer. Utilizing this scheme we can systematically analyze the interactions between the electric field and the bilayer components involved in pore initiation, construction and resealing. We find that the pore creation time depends strongly on the electric field gradient across the membrane interface and that the pore annihilation time is at least weakly dependent on the magnitude of the pore-initiating electric field and, in general, much longer than the pore creation time.     

Nanopore-facilitated, voltage-driven phosphatidylserine translocation in lipid bilayers--in cells and in silico

Nanosecond, megavolt-per-meter pulses--higher power but lower total energy than the electroporative pulses used to introduce normally excluded material into biological cells--produce large intracellular electric fields without destructively charging the plasma membrane. Nanoelectropulse perturbation of mammalian cells causes translocation of phosphatidylserine (PS) to the outer face of the cell, intracellular calcium release, and in some cell types a subsequent progression to apoptosis. Experimental observations and molecular dynamics (MD) simulations of membranes in pulsed electric fields presented here support the hypothesis that nanoelectropulse-induced PS externalization is driven by the electric potential that appears across the lipid bilayer during a pulse and is facilitated by the poration of the membrane that occurs even during pulses as brief as 3 ns. MD simulations of phospholipid bilayers in supraphysiological electric fields show a tight association between PS externalization and membrane pore formation on a nanosecond time scale that is consistent with experimental evidence for electropermeabilization and anode-directed PS translocation after nanosecond electric pulse exposure, suggesting a molecular mechanism for nanoelectroporation and nanosecond PS externalization: electrophoretic migration of the negatively charged PS head group along the surface of nanometer-diameter electropores initiated by field-driven alignment of water dipoles at the membrane interface.     

Membrane electroporation--fast molecular exchange by electroosmosis.

Human and rabbit erythrocyte ghosts loaded with FITC-dextran (mol. mass = 10 kDa) and NBD-glucosamine (mol. mass = 342 Da) in buffers of different ionic strength and composition were subjected to electric pulses (intensity 0.7 kV/mm and decay half-time 1 ms) at 7-10 degrees C and 20-24 degrees C. The transfer of the fluorescent dyes from the interior of the ghosts through the electropores was observed by low light level video microscopy. The pulses caused the fluorescence to appear outside the membranes as a transient cylindrical cloud directed toward the negative electrode during the first video frame (17 ms). It was similar in both rabbit and human erythrocyte ghosts and at both temperatures but differs for the two dyes, the fluorescence cylinder is long and tall for the FITC-dextran and relatively short and thick for the NBD-glucosamine. The molecular exchange was 2-3 orders of magnitude faster within the first 17 ms after the pulse than the diffusional exchange. It decreased with increasing ionic strength. Formulae for the transfer of molecules by electroosmotic flow through the pores are in agreement with these observations. They allow estimation of the total area of pores with radii larger than that of the fluorescent dye during the pulse. The major conclusion is that electroosmosis is the dominating mechanism of molecular exchange in electroporation of erythrocyte ghosts.     

Synergistic Combination of Electrolysis and Electroporation for Tissue Ablation

Electrolysis, electrochemotherapy with reversible electroporation, nanosecond pulsed electric fields and irreversible electroporation are valuable non-thermal electricity based tissue ablation technologies. This paper reports results from the first large animal study of a new non-thermal tissue ablation technology that employs “Synergistic electrolysis and electroporation” (SEE). The goal of this pre-clinical study is to expand on earlier studies with small animals and use the pig liver to establish SEE treatment parameters of clinical utility. We examined two SEE methods. One of the methods employs multiple electrochemotherapy-type reversible electroporation magnitude pulses, designed in such a way that the charge delivered during the electroporation pulses generates the electrolytic products. The second SEE method combines the delivery of a small number of electrochemotherapy magnitude electroporation pulses with a low voltage electrolysis generating DC current in three different ways. We show that both methods can produce lesion with dimensions of clinical utility, without the need to inject drugs as in electrochemotherapy, faster than with conventional electrolysis and with lower electric fields than irreversible electroporation and nanosecond pulsed ablation.     

Modeling electroporation in a single cell

Electroporation uses electric pulses to promote delivery of DNA and drugs into cells. This study presents a model of electroporation in a spherical cell exposed to an electric field. The model determines transmembrane potential, number of pores, and distribution of pore radii as functions of time and position on the cell surface. For a 1-ms, 40 kV/m pulse, electroporation consists of three stages: charging of the cell membrane (0-0.51 micros), creation of pores (0.51-1.43 micros), and evolution of pore radii (1.43 micros to 1 ms). This pulse creates approximately 341,000 pores, of which 97.8% are small ( approximately 1 nm radius) and 2.2% are large. The average radius of large pores is 22.8 +/- 18.7 nm, although some pores grow to 419 nm. The highest pore density occurs on the depolarized and hyperpolarized poles but the largest pores are on the border of the electroporated regions of the cell. Despite their much smaller number, large pores comprise 95.3% of the total pore area and contribute 66% to the increased cell conductance. For stronger pulses, pore area and cell conductance increase, but these increases are due to the creation of small pores; the number and size of large pores do not increase.     

The electrical breakdown of cell and lipid membranes: the similarity of phenomenologies

The current responses of human erythrocyte and L-cell membranes being subject to rectangular voltage pulses of 150-700 mV amplitude and 5 X 10(-3)-10 s duration were recorded by means of the patch-clamp method. The behaviour of planar lipid bilayer membranes of oxidized cholesterol and UO2(2+)-modified bilayers of azolectin in a high electric field was investigated for comparison. The gradual growth in the conductance (reversible electrical breakdown) was found for both the cell membranes and lipid bilayers of the compositions studied, with the application of voltage pulses of sufficient duration, to be completed by its drastic enhancement (irreversible breakdown). The time interval preceding the irreversible breakdown and the rate of increase in conductance during the reversible breakdown are determined by the amplitude of the voltage applied. The recovery of the initial properties of the membrane following the reversible breakdown consists of the two stages, the latter substantially differing by their characteristic times. The first very rapid stage (tau much less than 1 ms) reflects the lowering of the conductance of small pores with decreasing voltage across the membrane. The diminishing of the number and mean radii of the pores resulting in their complete disappearance occurs only at the second stage of membrane healing, which lasts several seconds or even minutes. The phenomenological similarity of the cell and lipid membrane breakdown indicates that pores developed during the electrical breakdown of biological membranes arise in their lipid matrices. The structure and the properties of the pores are discussed.     

Hydrofluoric and nitric acid transport through lipid bilayer membranes

Hydrofluoric and nitric acid transport through lipid bilayer membranes were studied by a combination of electrical conductance and pH electrode techniques. Transport occurs primarily by nonionic diffusion of molecular HF and HNO3. Membrane permeabilities to HF and HNO3 ranged from 10(-4) to 10(-3) cm . s-1, five to seven orders of magnitude higher than the permeabilities to NO-3, F- and H+. Our results are consistent with the hypothesis that F- transport through biological membranes occurs mainly by nonionic diffusion of HF. Our results also suggest that of the two principal components of 'acid rain', HNO3 may be more toxic than H2SO4.     

Mechanisms for the Intracellular Manipulation of Organelles by Conventional Electroporation

Conventional electroporation (EP) changes both the conductance and molecular permeability of the plasma membrane (PM) of cells and is a standard method for delivering both biologically active and probe molecules of a wide range of sizes into cells. However, the underlying mechanisms at the molecular and cellular levels remain controversial. Here we introduce a mathematical cell model that contains representative organelles (nucleus, endoplasmic reticulum, mitochondria) and includes a dynamic EP model, which describes formation, expansion, contraction, and destruction for the plasma and all organelle membranes. We show that conventional EP provides transient electrical pathways into the cell, sufficient to create significant intracellular fields. This emerging intracellular electrical field is a secondary effect due to EP and can cause transmembrane voltages at the organelles, which are large enough and long enough to gate organelle channels, and even sufficient, at some field strengths, for the poration of organelle membranes. This suggests an alternative to nanosecond pulsed electric fields for intracellular manipulations.     

Effect of actin cytoskeleton disruption on electric pulse‐induced apoptosis and electroporation in tumour cells

Electric pulses are known to affect the outer membrane and intracellular structures of tumour cells. By applying electrical pulses of 450 ns duration with electric field intensity of 8 kV/cm to HepG2 cells for 30 s, electric pulse‐induced changes in the integrity of the plasma membrane, apoptosis, viability and mitochondrial transmembrane potential were investigated. Results demonstrated that electric pulses induced cell apoptosis and necrosis accompanied with the decrease of mitochondrial transmembrane potential and the formation of pores in the membrane. The role of cytoskeleton in cellular response to electric pulses was investigated. We found that the apoptotic and necrosis percentages of cells in response to electric pulses decreased after cytoskeletal disruption. The electroporation of cell was not affected by cytoskeletal disruption. The results suggest that the disruption of actin skeleton is positive in protecting cells from killing by electric pulses, and the skeleton is not involved in the electroporation directly.     

Electric field-induced functional reductions in the K+ channels mainly resulted from supramembrane potential-mediated electroconformational changes.

The goal of this study is to distinguish the supramembrane potential difference-induced electroconformational changes from the huge transmembrane current-induced thermal damages in the delayed rectifier K+ channels. A double Vaseline-gap voltage clamp was used to deliver shock pulses and to monitor the channel currents. Three pairs of 4-ms shock pulses were used to mimic the electric shock by a power-line frequency electric field. Each pair consists of two pulses with the same magnitude, starting from 350 to 500 mV, but with opposite polarities. The shock pulse-generated transmembrane ion flux and the responding electric energy, the Joule heating, consumed in the cell membrane, as well as the effects on the K+ channel currents, were obtained. Results showed that huge transmembrane currents are not necessary to cause damages in the K+ channel proteins. In contrast, reductions in the K+ channel currents are directly related to the field-induced supramembrane potential differences. By a comparison with the shock field-induced Joule heating effects on cell membranes, the field-induced supramembrane potential difference plays a dominant role in damaging the K+ channels, resulting in electroconformational changes in the membrane proteins. In contrast, the shock field-induced huge transmembrane currents, therefore the thermal effects, play a secondary, trivial role.     

Model of creation and evolution of stable electropores for DNA delivery

Electroporation, in which electric pulses create transient pores in the cell membrane, is becoming an important technique for gene therapy. To enable entry of supercoiled DNA into cells, the pores should have sufficiently large radii (>10 nm), remain open long enough for the DNA chain to enter the cell (milliseconds), and should not cause membrane rupture. This study presents a model that can predict such macropores. The distinctive features of this model are the coupling of individual pores through membrane tension and the electrical force on the pores, which is applicable to pores of any size. The model is used to explore the process of pore creation and evolution and to determine the number and size of pores as a function of the pulse magnitude and duration. Next, our electroporation model is combined with a heuristic model of DNA uptake and used to predict the dependence of DNA uptake on pulsing parameters. Finally, the model is used to examine the mechanism of a two-pulse protocol, which was proposed specifically for gene delivery. The comparison between experimental results and the model suggests that this model is well-suited for the investigation of electroporation-mediated DNA delivery.     

A long-lived fusogenic state is induced in erythrocyte ghosts by electric pulses

Treatment of erythrocyte ghosts in random positions in a suspension with membrane fusion-inducing direct current electric field pulses causes the membranes to become fusogenic. Significant fusion yields are observed if the membranes are dielectrophoretically aligned into membrane-membrane contact with a weak alternating electric field as much as 5 min after the application of the pulses. This demonstrates that a long-lived membrane structural alteration is involved in this fusion mechanism. Other experiments indicate that the areas on the membrane which become fusogenic after treatment with the pulses may be very highly localized. The locations of these fusogenic areas coincide with where the trans-membrane electric field strength was greatest during the pulse. The fusogenic membrane alteration, or components thereof, in these areas laterally diffuses very slowly or not at all, or, to be fusogenic, must be present at concentrations in the membrane above a certain threshold. The loss of soluble 0.9-3-nm-diameter fluorescent probes from resealed cytoplasmic compartments of randomly positioned erythrocyte ghosts occurs through electric field pulse-induced pores only during a pulse but not between pulses or after a train of pulses if the probe diameter is 1.2 nm or greater. For a given pulse treatment of membranes in random positions in suspensions, an increase in ionic strength of the medium results in (a) a decrease in loss during the pulse, (b) no difference in loss between pulses, and (c) an increase in fusion yield when membrane-membrane contact is established. The latter two results (b and c) are incompatible with a fusion mechanism that proposes a simple relationship between electric field-induced pores and fusion.     

The effect of hexadecaprenol on molecular organisation and transport properties of model membranes

The Langmuir monolayer technique and voltammetric analysis were used to investigate the properties of model lipid membranes prepared from dioleoylphosphatidylcholine (DOPC), hexadecaprenol (C80), and their mixtures. Surface pressure-molecular area isotherms, current-voltage characteristics, and membrane conductance-temperature were measured. Molecular area isobars, specific molecular areas, excess free energy of mixing, collapse pressure and collapse area were determined for lipid monolayers. Membrane conductance, activation energy of ion migration across the membrane, and membrane permeability coefficient for chloride ions were determined for lipid bilayers. Hexadecaprenol decreases the activation energy and increases membrane conductance and membrane permeability coefficient. The results of monolayer and bilayer investigations show that some electrical, transport and packing properties of lipid membranes change under the influence of hexadecaprenol. The results indicate that hexadecaprenol modulates the molecular organisation of the membrane and that the specific molecular area of polyprenol molecules depends on the relative concentration of polyprenols in membranes. We suggest that hexadecaprenol modifies lipid membranes by the formation of fluid microdomains. The results also indicate that electrical transmembrane potential can accelerate the formation of pores in lipid bilayers modified by long chain polyprenols.     

Molecular dynamics simulations of hydrophilic pores in lipid bilayers

Hydrophilic pores are formed in peptide free lipid bilayers under mechanical stress. It has been proposed that the transport of ionic species across such membranes is largely determined by the existence of such meta-stable hydrophilic pores. To study the properties of these structures and understand the mechanism by which pore expansion leads to membrane rupture, a series of molecular dynamics simulations of a dipalmitoylphosphatidylcholine (DPPC) bilayer have been conducted. The system was simulated in two different states; first, as a bilayer containing a meta-stable pore and second, as an equilibrated bilayer without a pore. Surface tension in both cases was applied to study the formation and stability of hydrophilic pores inside the bilayers. It is observed that below a critical threshold tension of approximately 38 mN/m the pores are stabilized. The minimum radius at which a pore can be stabilized is 0.7 nm. Based on the critical threshold tension the line tension of the bilayer was estimated to be approximately 3 x 10(-11) N, in good agreement with experimental measurements. The flux of water molecules through these stabilized pores was analyzed, and the structure and size of the pores characterized. When the lateral pressure exceeds the threshold tension, the pores become unstable and start to expand causing the rupture of the membrane. In the simulations the mechanical threshold tension necessary to cause rupture of the membrane on a nanosecond timescale is much higher in the case of the equilibrated bilayers, as compared with membranes containing preexisting pores.     

A laser-temperature-jump method for the study of the rate of transfer of hydrophobic ions and carriers across the interface of thin lipid membranes

The first application of a laser-temperature-jump apparatus for the study of ion transport through planar (artificial) lipid membranes is described. The relaxation of the electric current is detected, either continuously at a constant applied voltage or discontinuously by a series of short voltage pulses. The second technique, a combined voltage- and temperature-jump method, is especially appropriate to investigate the kinetics of the adsorption/desorption process of hydrophobic ions and neutral carriers of cations at the membrane interface and to separate this phenomenon from the diffusion process through the unstirred aqueous layers adjacent to the membrane. The aim is to determine the rate-limiting step of transport. The permeation rate of the hydrophobic anion 2,4,6-trinitrophenolate is limited by the inner membrane barrier. For tetraphenylberate the rate constant of translocation across the inner barrier and that of desorption from the membrane into water are found to be of comparable magnitude. The membrane permeability of the neutral macrocyclic ion carrier enniatin B is strongly interface limited by its comparatively small rate of desorption into water. These results show that the frequently used a priori assumption of partition equilibrium at the membrane interfaces during transport is not justified.     

The effects of intense submicrosecond electrical pulses on cells

A simple electrical model for living cells predicts an increasing probability for electric field interactions with intracellular substructures of both prokaryotic and eukaryotic cells when the electric pulse duration is reduced into the sub-microsecond range. The validity of this hypothesis was verified experimentally by applying electrical pulses (durations 100 micros-60 ns, electric field intensities 3-150 kV/cm) to Jurkat cells suspended in physiologic buffer containing propidium iodide. Effects on Jurkat cells were assessed by means of temporally resolved fluorescence and light microscopy. For the longest applied pulses, immediate uptake of propidium iodide occurred consistent with electroporation as the cause of increased surface membrane permeability. For nanosecond pulses, more delayed propidium iodide uptake occurred with significantly later uptake of propidium iodide occurring after 60 ns pulses compared to 300 ns pulses. Cellular swelling occurred rapidly following 300 ns pulses, but was minimal following 60 ns pulses. These data indicate that submicrosecond pulses achieve temporally distinct effects on living cells compared to microsecond pulses. The longer pulses result in rapid permeability changes in the surface membrane that are relatively homogeneous across the cell population, consistent with electroporation, while shorter pulses cause surface membrane permeability changes that are temporally delayed and heterogeneous in their magnitude.