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.     

Nanoscale, Electric Field-Driven Water Bridges in Vacuum Gaps and Lipid Bilayers

Formation of a water bridge across the lipid bilayer is the first stage of pore formation in molecular dynamic (MD) simulations of electroporation, suggesting that the intrusion of individual water molecules into the membrane interior is the initiation event in a sequence that leads to the formation of a conductive membrane pore. To delineate more clearly the role of water in membrane permeabilization, we conducted extensive MD simulations of water bridge formation, stabilization, and collapse in palmitoyloleoylphosphatidylcholine bilayers and in water–vacuum–water systems, in which two groups of water molecules are separated by a 2.8 nm vacuum gap, a simple analog of a phospholipid bilayer. Certain features, such as the exponential decrease in water bridge initiation time with increased external electric field, are similar in both systems. Other features, such as the relationship between water bridge lifetime and the diameter of the water bridge, are quite different between the two systems. Data such as these contribute to a better and more quantitative understanding of the relative roles of water and lipid in membrane electropore creation and annihilation, facilitating a mechanism-driven development of electroporation protocols. These methods can be extended to more complex, heterogeneous systems that include membrane proteins and intracellular and extracellular membrane attachments, leading to more accurate models of living cells in electric fields.     

On the Electroporation Thresholds of Lipid Bilayers: Molecular Dynamics Simulation Investigations

Electroporation relates to the cascade of events that follows the application of high electric fields and that leads to cell membrane permeabilization. Despite a wide range of applications, little is known about the electroporation threshold, which varies with membrane lipid composition. Here, using molecular dynamics simulations, we studied the response of dipalmitoyl-phosphatidylcholine, diphytanoyl-phosphocholine-ester and diphytanoyl-phosphocholine-ether lipid bilayers to an applied electric field. Comparing between lipids with acyl chains and methyl branched chains and between lipids with ether and ester linkages, which change drastically the membrane dipole potential, we found that in both cases the electroporation threshold differed substantially. We show, for the first time, that the electroporation threshold of a lipid bilayer depends not only on the “electrical” properties of the membrane, i.e., its dipole potential, but also on the properties of its component hydrophobic tails.     

Molecular-Level Characterization of Lipid Membrane Electroporation using Linearly Rising Current

We present experimental and theoretical results of electroporation of small patches of planar lipid bilayers by means of linearly rising current. The experiments were conducted on ~120-μm-diameter patches of planar phospholipid bilayers. The steadily increasing voltage across the bilayer imposed by linearly increasing current led to electroporation of the membrane for voltages above a few hundred millivolts. This method shows new molecular mechanisms of electroporation. We recorded small voltage drops preceding the breakdown of the bilayer due to irreversible electroporation. These voltage drops were often followed by a voltage re-rise within a fraction of a second. Modeling the observed phenomenon by equivalent electric circuits showed that these events relate to opening and closing of conducting pores through the bilayer. Molecular dynamics simulations performed under similar conditions indicate that each event is likely to correspond to the opening and closing of a single pore of about 5 nm in diameter, the conductance of which ranges in the 100-nS scale. This combined experimental and theoretical investigation provides a better quantitative characterization of the size, conductance and lifetime of pores created during lipid bilayer electroporation. Such a molecular insight should enable better control and tuning of electroporation parameters for a wide range of biomedical and biotechnological applications.     

Stochastic model for electric field-induced membrane pores. Electroporation

Electric impulses (1-20 kV cm-1, 1-5 microseconds) cause transient structural changes in biological membranes and lipid bilayers, leading to apparently reversible pore formation ( electroporation ) with cross-membrane material flow and, if two membranes are in contact, to irreversible membrane fusion ( electrofusion ). The fundamental process operative in electroporation and electrofusion is treated in terms of a periodic lipid block model, a block being a nearest-neighbour pair of lipid molecules in either of two states: (i) the polar head group in the bilayer plane or (ii) facing the centre of a pore (or defect site). The number of blocks in the pore wall is the stochastic variable of the model describing pore size and stability. The Helmholtz free energy function characterizing the transition probabilities of the various pore states contains the surface energies of the pore wall and the planar bilayer and, if an electric field is present, also a dielectric polarization term (dominated by the polarization of the water layer adjacent to the pore wall). Assuming a Poisson process the average number of blocks in a pore wall is given by the solution of a non-linear differential equation. At subcritical electric fields the average pore size is stationary and very small. At supercritical field strengths the pore radius increases and, reaching a critical pore size, the membrane ruptures (dielectric breakdown). If, however, the electric field is switched off, before the critical pore radius is reached, the pore apparently completely reseals to the closed bilayer configuration (reversible electroporation ).     

Membrane electroporation: a molecular dynamics simulation

We present results of molecular dynamics simulations of lipid bilayers under a high transverse electrical field aimed at investigating their electroporation. Several systems are studied, namely 1), a bare bilayer, 2), a bilayer containing a peptide nanotube channel, and 3), a system with a peripheral DNA double strand. In all systems, the applied transmembrane electric fields (0.5 V.nm(-1) and 1.0 V.nm(-1)) induce an electroporation of the lipid bilayer manifested by the formation of water wires and water channels across the membrane. The internal structures of the peptide nanotube assembly and that of the DNA strand are hardly modified under field. For system 2, no perturbation of the membrane is witnessed at the vicinity of the channel, which indicates that the interactions of the peptide with the nearby lipids stabilize the bilayer. For system 3, the DNA strand migrates to the interior of the membrane only after electroporation. Interestingly enough, switching of the external transmembrane potential in cases 1 and 2 for few nanoseconds is enough to allow for complete resealing and reconstitution of the bilayer. We provide evidence that the electric field induces a significant lateral stress on the bilayer, manifested by surface tensions of magnitudes in the order of 1 mN.m(-1). This study is believed to capture the essence of several dynamical phenomena observed experimentally and provides a framework for further developments and for new applications.     

Self-induced docking site of a deeply embedded peripheral membrane protein

As a first step toward understanding the principles of the targeting of C2 domains to membranes, we have carried out a molecular dynamics simulation of the C2 domain of cytosolic phospholipase A2 (cPLA2-C2) in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer at constant pressure and temperature (NPT, 300 K and 1 atm). Using the high-resolution crystal structure of cPLA2-C2 as a starting point, we embedded two copies of the C2 domain into a pre-equilibrated membrane at the depth and orientation previously defined by electron paramagnetic resonance (EPR). Noting that in the membrane-bound state the three calcium binding loops are complexed to two calcium ions, we initially restrained the calcium ions at the membrane depth determined by EPR. But the depth and orientation of the domains remained within EPR experimental errors when the restraints were later removed. We find that the thermally disordered, chemically heterogeneous interfacial zones of phosphatidylcholine bilayers allow local lipid remodeling to produce a nearly perfect match to the shape and polarity of the C2 domain, thereby enabling the C2 domain to assemble and optimize its own lipid docking site. The result is a cuplike docking site with a hydrophobic bottom and hydrophilic rim. Contrary to expectations, we did not find direct interactions between the protein-bound calcium ions and lipid headgroups, which were sterically excluded from the calcium binding cleft. Rather, the lipid phosphate groups provided outer-sphere calcium coordination through intervening water molecules. These results show that the combined use of high-resolution protein structures, EPR measurements, and molecular dynamics simulations provides a general approach for analyzing the molecular interactions between membrane-docked proteins and lipid bilayers.     

Electrical behavior and pore accumulation in a multicellular model for conventional and supra-electroporation

Extremely large but very short (20 kV/cm, 300 ns) electric field pulses were reported recently to non-thermally destroy melanoma tumors. The stated mechanism for field penetration into cells is pulse characteristic times faster than charge redistribution (displacement currents). Here we use a multicellular model with irregularly shaped, closely spaced cells to show that instead overwhelming pore creation (supra-electroporation) is dominant, with field penetration due to pores (ionic conduction currents) during most of the pulse. Moreover, the model's maximum membrane potential (about 1.2 V) is consistent with recent experimental observations on isolated cells. We also use the model to show that conventional electroporation resulting from 100 microsecond, 1 kV/cm pulses yields a spatially heterogeneous electroporation distribution. In contrast, the melanoma-destroying pulses cause nearly homogeneous electroporation of cells and their nuclear membranes. Electropores can persist for times much longer than the pulses, and are likely to be an important mechanism contributing to cell death.     

Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition

Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (T_m). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to T_m. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.     

The effect of calcium on the properties of charged phospholipid bilayers

We have performed molecular dynamics simulations to investigate the structure and dynamics of charged bilayers as well as the distribution of counterions at the bilayer interface. For this, we have considered the negatively charged di-myristoyl-phosphatidyl-glycerol (DMPG) and di-myristoyl-phosphatidyl-serine (DMPS) bilayers as well as a protonated di-myristoyl-phosphatidyl-serine (DMPSH) bilayer. We were particularly interested in calcium ions due to their important role in biological systems. Simulations performed in the presence of calcium ions (DMPG, DMPS) or sodium ions (DMPS) were run for 45-60 ns. Simulation results for DMPG are compared with fluorescence measurements. The average areas per molecule were 47.4+/-0.5 A2 (DMPG with calcium), 47.3+/-0.5 A2 (DMPS with calcium), 51.3+/-1.0 A2 (DMPS with sodium) and 45.3+/-0.5 A2 (DMPSH). The structure of the negatively charged lipids is significantly affected by the counterions, where calcium ions have a more pronounced effect than sodium ions. Calcium ions were found to be tightly bound to the anionic groups of the lipid molecules and as such appear to constitute an integral part of the membrane interface on nanoseconds time scales. In contrast to sodium ions, calcium ions are localised in a narrow (approximately 10 A) band around the phosphate group. The interaction of calcium with the lipid molecules enhances the molecular packing of the PG and PS lipids. This observation is in good agreement with emission spectra of the membrane partitioning probe Laurdan in DMPG multilamellar vesicles that indicate an increase in the ordering of the DMPG bilayer due to the presence of calcium. Our results indicate that calcium ions, which often function as a second messengers in living cells have a pronounced effect on membrane structures, which may have implications during signal transduction events.     

Comparison of membrane electroporation and protein denature in response to pulsed electric field with different durations

In this paper, we compared the minimum potential differences in the electroporation of membrane lipid bilayers and the denaturation of membrane proteins in response to an intensive pulsed electric field with various pulse durations. Single skeletal muscle fibers were exposed to a pulsed external electric field. The field‐induced changes in the membrane integrity (leakage current) and the Na channel currents were monitored to identify the minimum electric field needed to damage the membrane lipid bilayer and the membrane proteins, respectively. We found that in response to a relatively long pulsed electric shock (longer than the membrane intrinsic time constant), a lower membrane potential was needed to electroporate the cell membrane than for denaturing the membrane proteins, while for a short pulse a higher membrane potential was needed. In other words, phospholipid bilayers are more sensitive to the electric field than the membrane proteins for a long pulsed shock, while for a short pulse the proteins become more vulnerable. We can predict that for a short or ultrashort pulsed electric shock, the minimum membrane potential required to start to denature the protein functions in the cell plasma membrane is lower than that which starts to reduce the membrane integrity. Bioelectromagnetics 34:253–263, 2013. © 2012 Wiley Periodicals, Inc.     

Electric Field-Driven Water Dipoles: Nanoscale Architecture of Electroporation

Electroporation is the formation of permeabilizing structures in the cell membrane under the influence of an externally imposed electric field. The resulting increased permeability of the membrane enables a wide range of biological applications, including the delivery of normally excluded substances into cells. While electroporation is used extensively in biology, biotechnology, and medicine, its molecular mechanism is not well understood. This lack of knowledge limits the ability to control and fine-tune the process. In this article we propose a novel molecular mechanism for the electroporation of a lipid bilayer based on energetics analysis. Using molecular dynamics simulations we demonstrate that pore formation is driven by the reorganization of the interfacial water molecules. Our energetics analysis and comparisons of simulations with and without the lipid bilayer show that the process of poration is driven by field-induced reorganization of water dipoles at the water-lipid or water-vacuum interfaces into more energetically favorable configurations, with their molecular dipoles oriented in the external field. Although the contributing role of water in electroporation has been noted previously, here we propose that interfacial water molecules are the main players in the process, its initiators and drivers. The role of the lipid layer, to a first-order approximation, is then reduced to a relatively passive barrier. This new view of electroporation simplifies the study of the problem, and opens up new opportunities in both theoretical modeling of the process and experimental research to better control or to use it in new, innovative ways.     

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.     

Self-electroporation as a model for fusion pore formation

The creation of a small opening called the fusion pore is a necessary prerequisite for neurotransmitter release from synaptic vesicles. It is known that high intensity electric fields can create pores in vesicles by a process called electroporation. Due to the presence of charged phosphatidylserine (PS) molecules on the inner leaflet of the cell membrane, an electric field that is strong enough to cause electroporation of a synaptic vesicle might be present. It was shown by K. Rosenheck [K. Rosenheck. Biophys J 75, 1237-1243 (1998)] that in a planar geometry, fields sufficient to cause electroporation can occur at intermembrane separations of less than approximately 3 nm. It is frequently found, however, that the cell membrane is not planar but caves inward at the locations where a vesicle is close to it. Indentation of the cell membrane in the fusion region was modelled as a hemisphere and a theoretical study of the electric field in the vicinity of the cell membrane taking into account the screening effect of dissolved ions in the cytoplasm was performed. It was discovered that fields crossing the electroporation threshold occurred at a distance of 2 nm or less, supporting the claim that electroporation could be a possible mechanism for fusion pore formation.     

Calcium-binding to non-ionic lipids: crystal structure of a calcium chloride complex of geraniol.

X-ray diffraction data were used to determine the crystal structure of a calcium chloride complex of geraniol. The geraniol molecules assume a bilayer arrangement, with channels of calcium and chloride ions separating the bilayers. Each calcium ion is coordinated to the hydroxyl groups of two symmetry-related geraniol molecules and to four chloride ions. Our results demonstrate that hydrophobic interactions within a lipid bilayer can lead to an arrangement of hydroxyl groups suitable for binding calcium ions. Similar interactions may be involved in the calcium-binding sites on membrane surfaces.     

Ion solvation and structural stability in a sodium channel investigated by molecular dynamics calculations

The stability and ion binding properties of the homo-tetrameric pore domain of a prokaryotic, voltage-gated sodium channel are studied by extensive all-atom molecular dynamics simulations, with the channel protein being embedded in a fully hydrated lipid bilayer. It is found that Na(+) ion presents in a mostly hydrated state inside the wide pore of the selectivity filter of the sodium channel, in sharp contrast to the nearly fully dehydrated state for K(+) ions in potassium channels. Our results also indicate that Na(+) ions make contact with only one or two out of the four polypeptide chains forming the selectivity filter, and surprisingly, the selectivity filter exhibits robust stability for various initial ion configurations even in the absence of ions. These findings are quite different from those in potassium channels. Furthermore, an electric field above 0.5V/nm is suggested to be able to induce Na(+) permeation through the selectivity filter.     

External TEA block of shaker K+ channels is coupled to the movement of K+ ions within the selectivity filter

Recent molecular dynamic simulations and electrostatic calculations suggested that the external TEA binding site in K+ channels is outside the membrane electric field. However, it has been known for some time that external TEA block of Shaker K+ channels is voltage dependent. To reconcile these two results, we reexamined the voltage dependence of block of Shaker K+ channels by external TEA. We found that the voltage dependence of TEA block all but disappeared in solutions in which K+ ions were replaced by Rb+. These and other results with various concentrations of internal K+ and Rb+ ions suggest that the external TEA binding site is not within the membrane electric field and that the voltage dependence of TEA block in K+ solutions arises through a coupling with the movement of K+ ions through part of the membrane electric field. Our results suggest that external TEA block is coupled to two opposing voltage-dependent movements of K+ ions in the pore: (a) an inward shift of the average position of ions in the selectivity filter equivalent to a single ion moving approximately 37% into the pore from the external surface; and (b) a movement of internal K+ ions into a vestibule binding site located approximately 13% into the membrane electric field measured from the internal surface. The minimal voltage dependence of external TEA block in Rb+ solutions results from a minimal occupancy of the vestibule site by Rb+ ions and because the energy profile of the selectivity filter favors a more inward distribution of Rb+ occupancy.     

Role of cholesterol flip-flop in oxidized lipid bilayers

We performed a series of molecular dynamics simulations of cholesterol (Chol) in nonoxidized 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and in binary mixtures of PLPC-oxidized-lipid-bilayers with 0–50% Chol concentration and oxidized lipids with hydroperoxide and aldehyde oxidized functional groups. From the 60 unbiased molecular dynamics simulations (total of 161 μs), we found that Chol inhibited pore formation in the aldehyde-containing oxidized lipid bilayers at concentrations greater than 11%. For both pure PLPC bilayer and bilayers with hydroperoxide lipids, no pores were observed at any Chol concentration. Furthermore, increasing cholesterol concentration led to a change of phase state from the liquid-disordered to the liquid-ordered phase. This condensing effect of Chol was observed in all systems. Data analysis shows that the addition of Chol results in an increase in bilayer thickness. Interestingly, we observed Chol flip-flop only in the aldehyde-containing lipid bilayer but neither in the PLPC nor the hydroperoxide bilayers. Umbrella-sampling simulations were performed to calculate the translocation free energies and the Chol flip-flop rates. The results show that Chol’s flip-flop rate depends on the lipid bilayer type, and the highest rate are found in aldehyde bilayers. As the main finding, we shown that Chol stabilizes the oxidized lipid bilayer by confining the distribution of the oxidized functional groups.     

Synergistic effect of electric field and lipid oxidation on the permeability of cell membranes

BackgroundStrong electric fields are known to affect cell membrane permeability, which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on how this affects the cell membrane permeability.MethodWe here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.Results and conclusionsWe show how oxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.General significanceThis study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.