Ionic conductivity of electroporated lipid bilayer membranes

Dielectric Dispersion Modulated Sensing of Yeast Suspension Electroporation

A specific pulsed electric field protocol can be used to induce electroporation. This is used in the food industry for yeast pasteurization, in laboratories for generic transfer and the medical field for cancer treatment. The sensing of electroporation can be done with simple ‘instantaneous’ voltage-current analysis. However, there are some intrinsic low-frequency phenomena superposing the electroporation current, such as electrode polarization. The biological media are non-homogeneous, giving them specific characterization in the broad frequency spectrum. For example, the cell barrier, i.e., cell membrane, causes so called β-dispersion in the frequency range of tens to thousands of kHz. Electroporation is a dynamic phenomenon characterized by altering the cell membrane permeability. In this work, we show that the impedance measurement at certain frequencies could be used to detect the occurrence of electroporation, i.e., dielectric dispersion modulated sensing. This approach may be used for the design and implementation of electroporation systems. Yeast suspension electroporation is simulated to show changes in the frequency spectrum. Moreover, the alteration depends on characteristics of the system. Three types of external buffers and their characteristics are evaluated.

On Modeling Diversity in Electrical Cellular Response: Data-Driven Approach

Electrical properties of biological cells and tissues possess valuable information that enabled numerous applications in biomedical engineering. The common foundation behind them is a numerical model that can predict electrical response of a single cell or a network of cells. We analyzed the past empirical observations to propose the first statistical model that accurately mimics biological diversity among animal cells, yeast cells, and bacteria. Based on membrane elasticity and cell migration mechanisms, we introduce a more realistic three-dimensional geometry generation procedure that captures membrane protrusions and retractions in adherent cells. Together, they form a model of diverse electrical response across multiple cell types. We experimentally verified the model with electrical impedance spectroscopy of a single human cervical carcinoma (HeLa) cell on a microelectrode array. The work is of particular relevance to medical diagnostic and therapeutic applications that involve exposure to electric and magnetic fields.

Membrane Electroporation and Electropermeabilization: Mechanisms and Models

Exposure of biological cells to high-voltage, short-duration electric pulses causes a transient increase in their plasma membrane permeability, allowing transmembrane transport of otherwise impermeant molecules. In recent years, large steps were made in the understanding of underlying events. Formation of aqueous pores in the lipid bilayer is now a widely recognized mechanism, but evidence is growing that changes to individual membrane lipids and proteins also contribute, substantiating the need for terminological distinction between electroporation and electropermeabilization. We first revisit experimental evidence for electrically induced membrane permeability, its correlation with transmembrane voltage, and continuum models of electropermeabilization that disregard the molecular-level structure and events. We then present insights from molecular-level modeling, particularly atomistic simulations that enhance understanding of pore formation, and evidence of chemical modifications of membrane lipids and functional modulation of membrane proteins affecting membrane permeability. Finally, we discuss the remaining challenges to our full understanding of electroporation and electropermeabilization.

Characterization of Cell Membrane Permeability In Vitro Part II: Computational Model of Electroporation-Mediated Membrane Transport

Electroporation is the process by which applied electric fields generate nanoscale defects in biological membranes to more efficiently deliver drugs and other small molecules into the cells. Due to the complexity of the process, computational models of cellular electroporation are difficult to validate against quantitative molecular uptake data. In part I of this two-part report, we describe a novel method for quantitatively determining cell membrane permeability and molecular membrane transport using fluorescence microscopy. Here, in part II, we use the data from part I to develop a two-stage ordinary differential equation model of cellular electroporation. We fit our model using experimental data from cells immersed in three buffer solutions and exposed to electric field strengths of 170 to 400 kV/m and pulse durations of 1 to 1000 μs. We report that a low-conductivity 4-(2-hydroxyethyl)-1 piperazineethanesulfonic acid buffer enables molecular transport into the cell to increase more rapidly than with phosphate-buffered saline or culture medium-based buffer. For multipulse schemes, our model suggests that the interpulse delay between two opposite polarity electric field pulses does not play an appreciable role in the resultant molecular uptake for delays up to 100 μs. Our model also predicts the per-pulse permeability enhancement decreases as a function of the pulse number. This is the first report of an ordinary differential equation model of electroporation to be validated with quantitative molecular uptake data and consider both membrane permeability and charging.

Properties of lipid electropores I: Molecular dynamics simulations of stabilized pores by constant charge imbalance

Molecular dynamics (MD) simulations have become a powerful tool to study electroporation (EP) in atomic detail. In the last decade, numerous MD studies have been conducted to model the effect of pulsed electric fields on membranes, providing molecular models of the EP process of lipid bilayers. Here we extend these investigations by modeling for the first time conditions comparable to experiments using long (μs-ms) low intensity (~kV/cm) pulses, by studying the characteristics of pores formed in lipid bilayers maintained at a constant surface tension and subject to constant charge imbalance. This enables the evaluation of structural (size) and electrical (conductance) properties of the pores formed, providing information hardly accessible directly by experiments. Extensive simulations of EP of simple phosphatidylcholine bilayers in 1M NaCl show that hydrophilic pores with stable radii (1-2.5 nm) form under transmembrane voltages between 420 and 630 mV, allowing for ionic conductance in the range of 6.4-29.5 nS. We discuss in particular these findings and characterize both convergence and size effects in the MD simulations. We further extend these studies in a follow-up paper (Rems et al., Bioelectrochemistry, Submitted), by proposing an improved continuum model of pore conductance consistent with the results from the MD simulations.

Mechanisms of transfer of bioactive molecules through the cell membrane by electroporation

A short review of biophysical mechanisms for electrotransfer of bioactive molecules through the cell membrane by using electroporation is presented. The concept of transient hydrophilic aqueous pores and membrane electroporation mechanisms of single cells and cells in suspension models are analyzed. Alongside the theoretical approach, some peculiarities of drug and gene electrotransfer into cells and applications in clinical trials are discussed.

Basic features of a cell electroporation model: illustrative behavior for two very different pulses

Science increasingly involves complex modeling. Here we describe a model for cell electroporation in which membrane properties are dynamically modified by poration. Spatial scales range from cell membrane thickness (5 nm) to a typical mammalian cell radius (10 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm), and can be used with idealized and experimental pulse waveforms. The model consists of traditional passive components and additional active components representing nonequilibrium processes. Model responses include measurable quantities: transmembrane voltage, membrane electrical conductance, and solute transport rates and amounts for the representative “long” and “short” pulses. The long pulse—1.5 kV/cm, 100 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μs—evolves two pore subpopulations with a valley at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sim}$$\end{document}∼5 nm, which separates the subpopulations that have peaks at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sim}$$\end{document}∼1.5 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sim}$$\end{document}∼12 nm radius. Such pulses are widely used in biological research, biotechnology, and medicine, including cancer therapy by drug delivery and nonthermal physical tumor ablation by causing necrosis. The short pulse—40 kV/cm, 10 ns—creates 80-fold more pores, all small (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<$$\end{document}<3 nm; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}∼1 nm peak). These nanosecond pulses ablate tumors by apoptosis. We demonstrate the model’s responses by illustrative electrical and poration behavior, and transport of calcein and propidium. We then identify extensions for expanding modeling capability. Structure-function results from MD can allow extrapolations that bring response specificity to cell membranes based on their lipid composition. After a pulse, changes in pore energy landscape can be included over seconds to minutes, by mechanisms such as cell swelling and pulse-induced chemical reactions that slowly alter pore behavior.

Transient oscillation of shape and membrane conductivity changes by field pulse-induced electroporation in nano-sized phospholipid vesicles

The results of electrooptical and conductometrical measurements on unilamellar lipid vesicles (of mean radius a = 90 nm), filled with 0.2 M NaCl solution, suspended in 0.33 M sucrose solution of 0.2 mM NaCl, and exposed to a stepwise decaying electric field (time constant τE = 154 μs) in the range 10 ≤ E0 (kV cm(-1)) ≤ 90, are analyzed in terms of cyclic changes in vesicle shape and vesicle membrane conductivity. The two peaks in the dichroitic turbidity relaxations reflect two cycles of rapid membrane electroporation and slower resealing of long-lived electropores. The field-induced changes reflect structural transitions between closed (C) and porated (P) membrane states, qualified by pores of type P1 and of type P2, respectively. The transient change in the membrane conductivity and the transient shape oscillation are based on changes in the pore density of the (larger) P2-pores along a hysteresis cycle. The P2-pore formation leads to transient net ion flows across the P2-pores and to transient changes in the membrane field. The kinetic data are numerically processed in terms of coupled structural relaxation modes. Using the torus-hole pore model, the mean inner pore radii are estimated to be r1 = 0.38 (±0.05) nm and r2 = 1.7 (±0.1) nm, respectively. The observation of a transient oscillation of membrane electroporation and of shape changes in a longer lasting external field pulse is suggestive of potential resonance enhancement, for instance, of electro-uptake by, and of electro-release of biogenic molecules from, biological cells in trains of long-lasting low-intensity voltage pulses.

Theoretical and experimental analysis of electroporated membrane conductance in cell suspension

An intense electric field can be applied to increase the membrane conductance G(m) and consequently, the conductivity of cell suspension. This phenomenon is called electroporation. This mechanism is used in a wide range of medical applications, genetic engineering, and therapies. Conductivity measurements of cell suspensions were carried out during application of electric fields from 40 to 165 kV/m. Experimental results were analyzed with two electroporation models: the asymptotic electroporation model was used to estimate G(m) at the beginning and at the end of electric field pulse, and the extended Kinosita electroporation model to increase G(m) linearly in time. The maximum G(m) was 1-7 × 10(4) S/m(2), and the critical angle (when the G(m) is insignificant) was 50°-65°. In addition, the sensitivity of electroporated membrane conductance to extracellular and cytoplasmatic conductivity and cell radius has been studied. This study showed that external conductivity and cell radius are important parameters affecting the pore-opening phenomenon. However, if the cell radius is larger than 7 μm in low conductivity medium, the cell dimensions are not so important.

Numerical study of the electroporation pulse shape effect on molecular uptake of biological cells

Background

In order to reduce the side-effects of chemotherapy, combined chemotherapy-electroporation (electrochemotherapy) has been suggested. Electroporation, application of appropriate electric pulses to biological cells, can significantly enhance molecular uptake of cells due to formation of transient pores in the cell membrane. It was experimentally demonstrated that the efficiency of electroporation is under the control of electric pulse parameters. However, the theoretical basis for these experimental results is not fully explained. In order to predict the outcome of experiments and optimize the efficiency of electroporation before each treatment, we developed a model to investigate the effect of pulse shape on efficiency of electroporation.

Results

Our model is based on a developed chemical-kinetics scheme and trapezium barrier model, while self-consistency was taken into account. This model is further supplemented with a molecular transport model to acquire the molecular uptake of cells. The investigated pulse shapes in this study were unipolar rectangular pulses with different rise and fall times, triangular, sinusoidal and bipolar rectangular pulses and also sinusoidal modulated unipolar pulses with different percentages of modulation. The obtained results from our modelling and simulations are in good agreement with previously published experimental results.

Conclusions

We therefore conclude that this model can be used to predict the effects of arbitrarily shaped electroporation pulses on cell membrane conductivity and molecular transport across the cell membrane.

The current-voltage relation for electropores with conductivity gradients

In electroporation, an electric field transiently permeabilizes the cell membrane to gain access to the cytoplasm, and to deliver active agents such as DNA, proteins, and drug molecules. Past work suggests that the permeabilization is caused by the formation of aqueous, conducting pores on the lipid membrane, which are also known as electropores. The current-voltage relation across the membrane-bound pores is critical for understanding and predicting electroporation. In this work, we solve the Nernst-Planck equations in a geometry encompassing an isolated electropore to investigate this relation. In particular, we study cases where the intra- and extracellular electrical conductivities differ. We first derive an analytical solution, which is subsequently validated with a direct numerical simulation using a finite volume method. The main result of the current work is a formula for the effective pore resistance as a function of the pore radius, the membrane thickness, and the intra- and extracellular conductivities. This formula can be incorporated into whole-cell or planar-membrane electroporation models for system-level prediction and understanding.

Gene electrotransfer: from biophysical mechanisms to in vivo applications : Part 1- Biophysical mechanisms

Electropulsation is one of the nonviral methods successfully used to deliver genes into living cells in vitro and in vivo. This approach shows promise in the field of gene and cellular therapies. The present review focuses on the processes supporting gene electrotransfer in vitro. In the first part, we will report the events occurring before, during, and after pulse application in the specific field of plasmid DNA electrotransfer at the cell level. A critical discussion of the present theoretical considerations about membrane electropermeabilization and the transient structures involved in the plasmid uptake follows in a second part.

Efficiency of the delivery of small charged molecules into cells in vitro

The effectiveness of the delivery of small charged molecules, including anticancer drugs into MH22 hepatoma cells in vitro was investigated. It was shown that for each kind of small molecules one can find a specific set of pulse strength-duration combinations that define electrotransfer of chosen compounds into the same amount of electroporated cells. Analysis of experimental data from the point of theory of hydrophilic aqueous pores and the estimation of the contribution of the electrostatic Born's energy to the change in free energy suggests that the main factors defining small molecules transfer through the membrane are: the charge and size of molecules, the permittivities of external medium, membrane material, and the electropores respectively as well as the size of electropores. The joint impact of all mentioned factors on transfer efficiency is essential.

Kinetics of transmembrane transport of small molecules into electropermeabilized cells

The transport of propidium iodide into electropermeabilized Chinese hamster ovary cells was monitored with a photomultiplier tube during and after the electric pulse. The influence of pulse amplitude and duration on the transport kinetics was investigated with time resolutions from 200 ns to 4 ms in intervals from 400 micros to 8 s. The transport became detectable as early as 60 micros after the start of the pulse, continued for tens of seconds after the pulse, and was faster and larger for higher pulse amplitudes and/or longer pulse durations. With fixed pulse parameters, transport into confluent monolayers of cells was slower than transport into suspended cells. Different time courses of fluorescence increase were observed during and at various times after the pulse, reflecting different transport mechanisms and ongoing membrane resealing. The data were compared to theoretical predictions of the Nernst-Planck equation. After a delay of 60 micros, the time course of fluorescence during the pulse was approximately linear, supporting a mainly electrophoretic solution of the Nernst-Planck equation. The time course after the pulse agreed with diffusional solution of the Nernst-Planck equation if the membrane resealing was assumed to consist of three distinct components, with time constants in the range of tens of microseconds, hundreds of microseconds, and tens of seconds, respectively.

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.

Electrical energy required to form large conducting pores

This study computes the contribution of the externally induced transmembrane potential to the energy of large, highly conductive pores. This work was undertaken because the pore energy formulas existing in the literature predict qualitatively different behavior of large pores: the original formula proposed by Abidor et al. in 1979 implies that the electrical force expanding the pore increases linearly with pore radius, while later extensions of this formula imply that this force decreases to zero for large pores. Starting from the Maxwell stress tensors, our study derives the formula for the mechanical work required to deform a dielectric body in an ionic solution with steady-state electric current. This formula is related to a boundary value problem (BVP) governing electric potentials and fields in a proximity of a pore. Computer simulations yield estimates of the electrical energy for pores of two different shapes: cylindrical and toroidal. In both cases, the energy increases linearly for pore radii above approximately 20 nm, implying that the electrical force expanding the pore asymptotes to a constant value for large pores. This result is different from either of the two energy formulas mentioned above. Our study traces the source of this disagreement to approximations made by previous studies, which are suitable only for small pores. Therefore, this study provides a better understanding of the energy of large pores, which is needed for designing pulsing protocols for DNA delivery.

Digression on membrane electroporation for drug and gene delivery

Membrane electroporation (ME) defines an electrical technique to render lipid membranes porous and permeable, transiently and reversibly, by external voltage pulses. Although there are numerous applications of ME to manipulate cells, organelles and tissues in cell biology, biotechnology and medicine, yet the molecular mechanism of ME is only slowly being understood. A general chemical- thermodynamical approach for the quantitative description of cell membrane electroporation has been developed to provide the framework to quantitatively rationalize electroporative cell transformation and electroporative uptake of drug-like dyes into cells, as well as electrolyte efflux from salt-filled electroporated vesicles. Mechanistically, the electroporative transfer of gene and drug-like dyes involves the coupling between an interactive contact formation of the permeates with the cell surface membrane and the structural electroporation-resealing cycle C <--> (P) where C is the closed and (P) represents a number of different porated membrane states, respectively. The experimentally accessible concentration fraction f(p) = [(P)] / ([C] + [(P)]) of porous states is related to thermodynamic and electro-mechanic parameters such as temperature and the electric field strength, membrane rigidity or curvature. The results of the theoretical approach, mainly based on electrooptical data of lipid vesicles, have been successfully used to analyze single cells and to specify conditions for the practical purpose of direct electroporative gene transfer and drug delivery, in particular in the new medical disciplines of electroporative chemotherapy and electroporative gene vaccination.