Membrane electroporation: The absolute rate equation and nanosecond time scale pore creation

Simulation study on electroporation of cancer cells in multicellular system

Electroporation (EP) of the normal cell and cancer cell both in single-cell and multicellular models was investigated by the meshed transport network method (MTNM) in this paper. The simulation results suggest that the cancer cell undergoes faster and more significant local EP than that of the corresponding normal cell induced by nanosecond pulsed electric fields (nsPEFs) both in single-cell and multicellular models. Furthermore, the results of the multicellular model indicate that there is a unidirectional neighboring effect in the multicellular model, meaning that cells at the center are affected and their pore formation is significantly reduced, but this effect is very weak for cells at the edges of the system. This means that the electric field selectively kills cells in different distribution locations. This work can provide guidance for the selection of parameters for the cancer cell EP process.

Numerical assessments of geometry, proximity and multi-electrode effects on electroporation in mitochondria and the endoplasmic reticulum to nanosecond electric pulses

Most simulations of electric field driven bioeffects have considered spherical cellular geometries or probed symmetrical structures for simplicity. This work assesses cellular transmembrane potential build-up and electroporation in a Jurkat cell that includes the endoplasmic reticulum (ER) and mitochondria, both of which have complex shapes, in response to external nanosecond electric pulses. The simulations are based on a time-domain nodal analysis that incorporates membrane poration utilizing the Smoluchowski model with angular-dependent changes in membrane conductivity. Consistent with prior experimental reports, the simulations show that the ER requires the largest electric field for electroporation, while the inner mitochondrial membrane (IMM) is the easiest membrane to porate. Our results suggest that the experimentally observed increase in intracellular calcium could be due to a calcium induced calcium release (CICR) process that is initiated by outer cell membrane breakdown. Repeated pulsing and/or using multiple electrodes are shown to create a stronger poration. The role of mutual coupling, screening, and proximity effects in bringing about electric field modifications is also probed. Finally, while including greater geometric details might refine predictions, the qualitative trends are expected to remain.

Electroporation in Cancer Therapy: A Simplified Model Derived from the Hodgkin-Huxley Model

Cancer remains a global health challenge, necessitating effective treatments with fewer side effects. Traditional methods such as chemotherapy and surgery often have complications. Pulsed electric fields and electroporation have emerged as promising approaches to mitigate these challenges. This study presents a comprehensive overview of electroporation as an innovative tool in cancer therapy, encompassing critical elements such as pulse generators and delivery devices. Furthermore, it introduces a simplified reversible electroporation model grounded in the Hodgkin-Huxley model. This model ensures resting potential stability by regulating ionic currents. When membrane charges reach the electroporation threshold, the model swiftly increases the fraction of open pores, resulting in a rapid rise in electroporation current. Conversely, as the transmembrane potential drops below the threshold, the model gradually reduces the fraction of open pores, leading to a gradual decline in electroporation current, indicating pore resealing. This model contributes to easier modeling and implementation of reversible electroporation dynamics, providing a valuable tool for further exploration of electroporation for cancer therapy.

Molecular Simulations Reveal the Free Energy Landscape and Transition State of Membrane Electroporation

The formation of pores over lipid membranes by the application of electric fields, termed membrane electroporation, is widely used in biotechnology and medicine to deliver drugs, vaccines, or genes into living cells. Continuum models for describing the free energy landscape of membrane electroporation were proposed decades ago, but they have never been tested against spatially detailed atomistic models. Using molecular dynamics (MD) simulations with a recently proposed reaction coordinate, we computed potentials of mean force of pore nucleation and pore expansion in lipid membranes at various transmembrane potentials. Whereas the free energies of pore expansion are compatible with previous continuum models, the experimentally important free energy barrier of pore nucleation is at variance with established models. The discrepancy originates from different geometries of the transition state; previous continuum models assumed the presence of a membrane-spanning defect throughout the process, whereas, according to the MD simulations, the transition state of pore nucleation is typically passed before a transmembrane defect has formed. A modified continuum model is presented that qualitatively agrees with the MD simulations. Using kinetics of pore opening together with transition state theory, our free energies of pore nucleation are in excellent agreement with previous experimental data.

Does the shape of the electric pulse matter in electroporation?

Electric pulses are widely used in biology, medicine, industry, and food processing. Numerous studies indicate that electroporation (EP) is a pulse-dependent process, and the electric pulse shape and duration strongly determine permeabilization efficacy. EP protocols are precisely planned in terms of the size and charge of the molecules, which will be delivered to the cell. In reversible and irreversible EP applications, rectangular or sine, polar or bipolar pulses are commonly used. The usage of pulses of the asymmetric shape is still limited to high voltage and low voltage (HV/LV) sequences in the context of gene delivery, while EP-based applications of ultra-short asymmetric pulses are just starting to emerge. This review emphasizes the importance and role of the pulse shape for membrane permeabilization by EP.

Nanosecond Pulsed Electric Field (nsPEF): Opening the Biotechnological Pandora’s Box

Nanosecond Pulsed Electric Field (nsPEF) is an electrostimulation technique first developed in 1995; nsPEF requires the delivery of a series of pulses of high electric fields in the order of nanoseconds into biological tissues or cells. They primary effects in cells is the formation of membrane nanopores and the activation of ionic channels, leading to an incremental increase in cytoplasmic Ca2+ concentration, which triggers a signaling cascade producing a variety of effects: from apoptosis up to cell differentiation and proliferation. Further, nsPEF may affect organelles, making nsPEF a unique tool to manipulate and study cells. This technique is exploited in a broad spectrum of applications, such as: sterilization in the food industry, seed germination, anti-parasitic effects, wound healing, increased immune response, activation of neurons and myocites, cell proliferation, cellular phenotype manipulation, modulation of gene expression, and as a novel cancer treatment. This review thoroughly explores both nsPEF’s history and applications, with emphasis on the cellular effects from a biophysics perspective, highlighting the role of ionic channels as a mechanistic driver of the increase in cytoplasmic Ca2+ concentration.

Modeling the gene delivery process of the needle array-based tissue nanotransfection

Hollow needle array-based tissue nanotransfection (TNT) presents an in vivo transfection approach that directly translocate exogeneous genes to target tissues by using electric pulses. In this work, the gene delivery process of TNT was simulated and experimentally validated. We adopted the asymptotic method and cell-array-based model to investigate the electroporation behaviors of cells within the skin structure. The distribution of nonuniform electric field across the skin results in various electroporation behavior for each cell. Cells underneath the hollow microchannels of the needle exhibited the highest total pore numbers compared to others due to the stronger localized electric field. The percentage of electroporated cells within the skin structure, with pore radius over 10 nm, increases from 25% to 82% as the applied voltage increases from 100 to 150 V/mm. Furthermore, the gene delivery behavior across the skin tissue was investigated through the multilayer-stack-based model. The delivery distance increased nonlinearly as the applied voltage and pulse number increased, which mainly depends on the diffusion characteristics and electric conductivity of each layer. It was also found that the skin is required to be exfoliated prior to the TNT procedure to enhance the delivery depth. This work provides the foundation for transition from the study of murine skin to translation use in large animals and human settings.

Identification of electroporation sites in the complex lipid organization of the plasma membrane

The plasma membrane of a biological cell is a complex assembly of lipids and membrane proteins, which tightly regulate transmembrane transport. When a cell is exposed to strong electric field, the membrane integrity becomes transiently disrupted by formation of transmembrane pores. This phenomenon termed electroporation is already utilized in many rapidly developing applications in medicine including gene therapy, cancer treatment, and treatment of cardiac arrhythmias. However, the molecular mechanisms of electroporation are not yet sufficiently well understood; in particular, it is unclear where exactly pores form in the complex organization of the plasma membrane. In this study, we combine coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis to identify how formation of pores depends on the local lipid organization. We show that pores do not form homogeneously across the membrane, but colocalize with domains that have specific features, the most important being high density of polyunsaturated lipids. We further show that knowing the lipid organization is sufficient to reliably predict poration sites with machine learning. Additionally, by analysing poration kinetics with Bayesian survival analysis we show that poration does not depend solely on local lipid arrangement, but also on membrane mechanical properties and the polarity of the electric field. Finally, we discuss how the combination of atomistic and coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis can guide the design of future experiments and help us to develop an accurate description of plasma membrane electroporation on the whole-cell level. Achieving this will allow us to shift the optimization of electroporation applications from blind trial-and-error approaches to mechanistic-driven design.

A double bilayer to study the nonequilibrium environmental response of GIRK2 in complex states

G protein-gated inwardly rectifying potassium (GIRK) channels play essential roles in electrical signaling in neurons and muscle cells. Nonequilibrium environments provide crucial driving forces behind many cellular events. Here, we apply the antiparallel alignment double bilayer model to study GIRK2 in response to the time-dependent membrane potential. Using molecular dynamics and umbrella sampling, we examined the time-dependent environmental impact on the ion conduction, energy basis, and primary motions of GIRK2 in different complex states with phosphatidylinositol-4,5-bisphosphate (PIP2) and G-protein βγ subunits (Gβγ). The antiparallel alignment double bilayer model enables us to study the transport performance in inward and outward K+ and mixed K+ and Na+. We obtained the recoverable discharge process of GIRK2 complexed with both PIP2 and Gβγ, compared with occasional conduction under PIP2-only regulation. Calculations of potential of mean force suggest different regulation by the helix bundle crossing (HBC) gate and G-loop gate regarding different complex states and under a membrane potential. In a nonequilibrium environment, distinct functional rocking motions of GIRK2 were identified under strengthened correlations between the transmembrane helices and downstream cytoplasmic domains with binding of PIP2, cations, and Gβγ. The findings suggest the potential domain motions and dynamics associated with a nonequilibrium environment and highlight the application of the antiparallel alignment double bilayer model to investigate factors in an asymmetric environment.

Continuum analysis to assess field enhancements for tailoring electroporation driven by monopolar or bipolar pulsing based on nonuniformly distributed nanoparticles

Recent reports indicate that nanoparticle (NP) clusters near cell membranes could enhance local electric fields, leading to heightened electroporation. This aspect is quantitatively analyzed through numerical simulations whereby time dependent transmembrane potentials are first obtained on the basis of a distributed circuit mode, and the results then used to calculate pore distributions from continuum Smoluchowski theory. For completeness, both monopolar and bipolar nanosecond-range pulse responses are presented and discussed. Our results show strong increases in TMP with the presence of multiple NP clusters and demonstrate that enhanced poration could be possible even over sites far away from the poles at the short pulsing regime. Furthermore, our results demonstrate that nonuniform distributions would work to enable poration at regions far away from the poles. The NP clusters could thus act as distributed electrodes. Our results were roughly in line with recent experimental observations.

Nonlinear dispersive cell model for microdosimetry of nanosecond pulsed electric fields

For applications based on nanosecond pulsed electric fields (nsPEFs), the underlying transmembrane potential (TMP) distribution on the plasma membrane is influenced by electroporation (EP) of the plasma membrane and dielectric dispersion (DP) of all cell compartments which is important for predicting the bioelectric effects. In this study, the temporal and spatial distribution of TMP on the plasma membrane induced by nsPEFs of various pulse durations (3 ns, 5 ns unipolar, 5 ns bipolar, and 10 ns) is investigated with the inclusion of both DP and EP. Based on the double-shelled dielectric spherical cell model, the Debye equation describing DP is transformed into the time-domain form with the introduction of polarization vector, and then we obtain the time course of TMP by solving the combination of Laplace equation and time-domain Debye equation. Next, the asymptotic version of the Smoluchowski equation is included to characterize the EP of plasma membrane in order to observe more profound electroporation effects with larger pore density and electroporated areas in consideration of both DP and EP. Through the simulation, it is clearer to understand the relationship between the applied nsPEFs and the induced bioelectric effects.

Nanofountain Probe Electroporation Enables Versatile Single-Cell Intracellular Delivery and Investigation of Postpulse Electropore Dynamics

Introducing exogenous molecules into cells with high efficiency and dosage control is a crucial step in basic research as well as clinical applications. Here, the capability of the nanofountain probe electroporation (NFP-E) system to deliver proteins and plasmids in a variety of continuous and primary cell types with appropriate dosage control is reported. It is shown that the NFP-E can achieve fine control over the relative expression of two cotransfected plasmids. Finally, the dynamics of electropore closure after the pulsing ends with the NFP-E is investigated. Localized electroporation has recently been utilized to demonstrate the converse process of delivery (sampling), in which a small volume of the cytosol is retrieved during electroporation without causing cell lysis. Single-cell temporal sampling confers the benefit of monitoring the same cell over time and can provide valuable insights into the mechanisms underlying processes such as stem cell differentiation and disease progression. NFP-E parameters that maximize the membrane resealing time, which is essential for increasing the sampled volume and in meeting the challenge of monitoring low copy number biomarkers, are identified. Its application in CRISPR/Cas9 gene editing, stem cell reprogramming, and single-cell sampling studies is envisioned.

Actin networks regulate the cell membrane permeability during electroporation

Transient physical disruption of cell membranes by electric pulses (or electroporation) has significance in biomedical and biological applications requiring the delivery of exogenous (bio)molecules to living cells. We demonstrate that actin networks regulate the cell membrane permeability during electroporation. Disruption of actin networks increases the uptake of membrane-impermeable molecules such as propidium iodide during electroporation. Our experiments at different temperatures ranging from 11 °C to 37 °C show that molecular uptake during electroporation increases with temperature. Furthermore, by examining the temperature-dependent kinetics of propidium iodide uptake, we infer that the activation energy barrier of electroporation is lowered when the actin networks are disrupted. Our numerical calculations of transmembrane voltage show that the reduced activation energy barrier for the cells with disrupted actin is not a consequence of the changes in transmembrane voltage associated with changes in the cell shape due to the disruption of actin, indicating that this could be due to changes in membrane mechanical properties. Our results suggest that the current theoretical models of electroporation should be advanced further by including the contributions of the cytoskeletal networks on the cell membrane permeability during the delivery of exogenous materials.

The interplay of excitation and electroporation in nanosecond pulse stimulation

Conventional electric stimuli of micro- and millisecond duration excite or activate cells at voltages 10-100 times below the electroporation threshold. This ratio is remarkably different for nanosecond electric pulses (nsEP), which caused excitation and activation only at or above the electroporation threshold in diverse cell lines, primary cardiomyocytes, neurons, and chromaffin cells. Depolarization to the excitation threshold often results from (or is assisted by) the loss of the resting membrane potential due to ion leaks across the membrane permeabilized by nsEP. Slow membrane resealing and the build-up of electroporation damages prevent repetitive excitation by nsEP. However, peripheral nerves and muscles are exempt from this rule and withstand multiple cycles of excitation by nsEP without the loss of function or signs of electroporation. We show that the damage-free excitation by nsEP may be enabled by the membrane charging time constant sufficiently large to (1) cap the peak transmembrane voltage during nsEP below the electroporation threshold, and (2) extend the post-nsEP depolarization long enough to activate voltage-gated ion channels. The low excitatory efficacy of nsEP compared to longer pulses makes them advantageous for medical applications where the neuromuscular excitation is an unwanted side effect, such as electroporation-based cancer and tissue ablation.

Irreversible Electroporation: Background, Theory, and Review of Recent Developments in Clinical Oncology

Irreversible electroporation (IRE) has established a clinical niche as an alternative to thermal ablation for the eradication of unresectable tumors, particularly those near critical vascular structures. IRE has been used in over 50 independent clinical trials and has shown clinical success when used as a standalone treatment and as a single component within combinatorial treatment paradigms. Recently, many studies evaluating IRE in larger patient cohorts and alongside other novel therapies have been reported. Here, we present the basic principles of reversible electroporation and IRE followed by a review of preclinical and clinical data with a focus on tumors in three organ systems in which IRE has shown great promise: the prostate, pancreas, and liver. Finally, we discuss alternative and future developments, which will likely further advance the use of IRE in the clinic.

The link between yeast cell wall porosity and plasma membrane permeability after PEF treatment

An investigation of the yeast cell resealing process was performed by studying the absorption of the tetraphenylphosphonium (TPP⁺) ion by the yeast Saccharomyces cerevisiae. It was shown that the main barrier for the uptake of such TPP⁺ ions is the cell wall. An increased rate of TPP⁺ absorption after treatment of such cells with a pulsed electric field (PEF) was observed only in intact cells, but not in spheroplasts. The investigation of the uptake of TPP⁺ in PEF treated cells exposed to TPP⁺ for different time intervals also showed the dependence of the absorption rate on the PEF strength. The modelling of the TPP⁺ uptake recovery has also shown that the characteristic decay time of the non-equilibrium (PEF induced) pores was approximately a few tens of seconds and this did not depend on the PEF strength. A further investigation of such cell membrane recovery process using a florescent SYTOX Green nucleic acid stain dye also showed that such membrane resealing takes place over a time that is like that occurring in the cell wall. It was thus concluded that the similar characteristic lifetimes of the non-equilibrium pores in the cell wall and membrane after exposure  to  PEF indicate a strong coupling between these parts of the cell.

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.

Visualization of Dynamic Sub-microsecond Changes in Membrane Potential

Direct observation of rapid membrane potential changes is critical to understand how complex neurological systems function. This knowledge is especially important when stimulation is achieved through an external stimulus meant to mimic a naturally occurring process. To enable exploration of this dynamic space, we developed an all-optical method for observing rapid changes in membrane potential at temporal resolutions of ∼25 ns. By applying a single 600-ns electric pulse, we observed sub-microsecond, continuous membrane charging and discharging dynamics. Close agreement between the acquired results and an analytical membrane-charging model validates the utility of this technique. This tool will deepen our understanding of the role of membrane potential dynamics in the regulation of many biological and chemical processes within living systems.

Electroporation of Skin Stratum Corneum Lipid Bilayer and Molecular Mechanism of Drug Transport: A Molecular Dynamics Study

The electroporation technique has been used significantly to increase drug permeation through the skin. This technique relies on the application of short-timed (microseconds to millisecond) electric fields (generally, order of 50--300 V) on the skin to create microscopic pores. However, the molecular mechanism of pore formation, resulting in an enhanced flux of active molecules through the skin, remains poorly understood. In this study, extensive atomistic molecular dynamics simulation of skin lipids [made up of ceramide (CER), cholesterol (CHOL), and free fatty acid (FFA)] has been performed at various external electric fields. We show for the first time the pore formation in the skin lipid bilayer during electroporation. We show the effect of the applied external electrical field (0.6-1.0 V/nm) on the pore formation dynamics in the lipid bilayer of different sizes (154, 616, and 2464 lipids) and compositions (CER/CHOL/FFA, 1:0:0, 1:0:1, 1:1:0, 1:1:1). The pore formation and resealing kinetics were different and were found to be highly dependent on the composition of the skin lipid bilayer. The pore formation time decreased with increase in the bilayer size. The pore sustaining electric field was found to be in the range of 0.20-0.25 V/nm for equimolar CER, CHOL, and FFA lipid bilayers. The skin lipid bilayer (1:1:1) sealed itself within 20 ns after the removal of the external electric field. We also present the molecular mechanism of enhancement of drug permeation in the presence of external field as compared to the passive diffusion. The molecular-level understanding obtained here could help in optimizing/designing the electroporation experiments for effective drug delivery. For a given skin composition and size of the drug molecule, the combination of pore formation time and pore growth model can be used to know a priori the desired electric field and time for the application of the electric field.

Core-shell magnetoelectric nanorobot - A remotely controlled probe for targeted cell manipulation

We have developed a remotely controlled dynamic process of manipulating targeted biological live cells using fabricated core-shell nanocomposites, which comprises of single crystalline ferromagnetic cores (CoFe₂O₄) coated with crystalline ferroelectric thin film shells (BaTiO₃). We demonstrate them as a unique family of inorganic magnetoelectric nanorobots (MENRs), controlled remotely by applied a.c. or d.c. magnetic fields, to perform cell targeting, permeation, and transport. Under a.c. magnetic field excitation (50 Oe, 60 Hz), the MENR acts as a localized electric periodic pulse generator and can permeate a series of misaligned cells, while aligning them to an equipotential mono-array by inducing inter-cellular signaling. Under a.c. magnetic field (40 Oe, 30 Hz) excitation, MENRs can be dynamically driven to a targeted cell, avoiding untargeted cells in the path, irrespective of cell density. D.C. magnetic field (−50 Oe) excitation causes the MENRs to act as thrust generator and exerts motion in a group of cells.

Lipid vesicles in pulsed electric fields: Fundamental principles of the membrane response and its biomedical applications

The present review focuses on the effects of pulsed electric fields on lipid vesicles ranging from giant unilamellar vesicles (GUVs) to small unilamellar vesicles (SUVs), from both fundamental and applicative perspectives. Lipid vesicles are the most popular model membrane systems for studying biophysical and biological processes in living cells. Furthermore, as vesicles are made from biocompatible and biodegradable materials, they provide a strategy to create safe and functionalized drug delivery systems in health-care applications. Exposure of lipid vesicles to pulsed electric fields is a common physical method to transiently increase the permeability of the lipid membrane. This method, termed electroporation, has shown many advantages for delivering exogenous molecules including drugs and genetic material into vesicles and living cells. In addition, electroporation can be applied to induce fusion between vesicles and/or cells. First, we discuss in detail how research on cell-size GUVs as model cell systems has provided novel insight into the basic mechanisms of cell electroporation and associated phenomena. Afterwards, we continue with a thorough overview how electroporation and electrofusion have been used as versatile methods to manipulate vesicles of all sizes in different biomedical applications. We conclude by summarizing the open questions in the field of electroporation and possible future directions for vesicles in the biomedical field.

Adrenal Chromaffin Cells Exposed to 5-ns Pulses Require Higher Electric Fields to Porate Intracellular Membranes than the Plasma Membrane: An Experimental and Modeling Study

Nanosecond-duration electric pulses (NEPs) can permeabilize the endoplasmic reticulum (ER), causing release of Ca²⁺ into the cytoplasm. This study used experimentation coupled with numerical modeling to understand the lack of Ca²⁺ mobilization from Ca²⁺-storing organelles in catecholamine-secreting adrenal chromaffin cells exposed to 5-ns pulses. Fluorescence imaging determined a threshold electric (E) field of 8 MV/m for mobilizing intracellular Ca²⁺ whereas whole-cell recordings of membrane conductance determined a threshold E-field of 3 MV/m for causing plasma membrane permeabilization. In contrast, a 2D numerical model of a chromaffin cell, which was constructed with internal structures representing a nucleus, mitochondrion, ER, and secretory granule, predicted that exposing the cell to the same 5-ns pulse electroporated the plasma and ER membranes at the same E-field amplitude, 3-4 MV/m. Agreement of the numerical simulations with the experimental results was obtained only when the ER interior conductivity was 30-fold lower than that of the cytoplasm and the ER membrane permittivity was twice that of the plasma membrane. A more realistic intracellular geometry for chromaffin cells in which structures representing multiple secretory granules and an ER showed slight differences in the thresholds necessary to porate the membranes of the secretory granules. We conclude that more sophisticated cell models together with knowledge of accurate dielectric properties are needed to understand the effects of NEPs on intracellular membranes in chromaffin cells, information that will be important for elucidating how NEPs porate organelle membranes in other cell types having a similarly complex cytoplasmic ultrastructure.

Uncertainty Quantification in Irreversible Electroporation Simulations

One recent area of cancer research is irreversible electroporation (IRE). Irreversible electroporation is a minimally invasive procedure where needle electrodes are inserted into the body to ablate tumor cells with electricity. The aim of this paper is to investigate how uncertainty in tissue and tumor conductivity propagate into final ablation predictions used for treatment planning. Two dimensional simulations were performed for a circular tumor surrounded by healthy tissue, and electroporated from two monopolar electrodes. The conductivity values were treated as random variables whose distributions were taken from published literature on the average and standard deviation of liver tissue and liver tumors. Three different Monte Carlo setups were simulated each at three different voltages. Average and standard deviation data was reported for a multitude of electrical field properties experienced by the tumor. Plots showing the variability in the electrical field distribution throughout the tumor are also presented.

Investigation of the morphological transition of a phospholipid bilayer membrane in an external electric field via molecular dynamics simulation

Elucidating the mechanisms for morphological transitions of the phospholipid bilayer membrane during cellular activity should lead to greater understanding of these membrane transitions and allow us to optimize biotechnologies such as drug delivery systems in organisms. To investigate the mechanism for and the dynamics of morphological changes in the phospholipid membrane, we performed molecular dynamics simulation of a phospholipid membrane with and without membrane protein under the influence of electric fields with different strengths. In the absence of membrane protein, it was possible to control the transition from one lamellar membrane morphology to another by applying a strong electric field. The strong electric field initially disordered the lipid molecules in the membrane, leading to the formation of a hydrophilic pore. The lipid molecules then spontaneously fused into a new lamellar membrane morphology. In the presence of membrane protein, a morphological transition from lamellar membrane to vesicle under the influence of a strong electric field was observed. Studying the complex transition dynamics associated with these changes in membrane morphology allowed us to gain deep insight into the electrofusion and electroporation that occur in the presence or absence of membrane protein, and the results obtained here should prove useful in work aimed at controlling membrane morphology. Graphical Abstract Memebrane morphological transition under the electric field of 0.6 V/nm with the membrane protein (down) and without membrane protein (up).

High-throughput Nuclear Delivery and Rapid Expression of DNA via Mechanical and Electrical Cell-Membrane Disruption

Nuclear transfection of DNA into mammalian cells is challenging yet critical for many biological and medical studies. Here, by combining cell squeezing and electric-field-driven transport in a device that integrates microfluidic channels with constrictions and microelectrodes, we demonstrate nuclear delivery of plasmid DNA within 1 hour after treatment, the most rapid DNA expression in a high-throughput setting (up to millions of cells per minute per device). Passing cells at high speed through microfluidic constrictions smaller than the cell diameter mechanically disrupts the cell membrane, allowing a subsequent electric field to further disrupt the nuclear envelope and drive DNA molecules into the cytoplasm and nucleus. By tracking the localization of the ESCRT-III (endosomal sorting complexes required for transport) protein CHMP4B, we show that the integrity of the nuclear envelope is recovered within 15 minutes of treatment. We also provide insight into subcellular delivery by comparing the performance of the disruption-and-field-enhanced method with those of conventional chemical, electroporation, and manual-injection systems.

Magneto-elasto-electroporation (MEEP): In-vitro visualization and numerical characteristics

A magnetically controlled elastically driven electroporation phenomenon, or magneto-elasto-electroporation (MEEP), is discovered while studying the interactions between core-shell magnetoelectric nanoparticles (CSMEN) and biological cells in the presence of an a.c. magnetic field. In this paper we report the effect of MEEP observed via a series of in-vitro experiments using core (CoFe2O4)-shell (BaTiO3) structured magnetoelectric nanoparticles and human epithelial cells (HEP2). The cell electroporation phenomenon and its correlation with the magnetic field modulated CSMEN are described in detail. The potential application of CSMEN in electroporation is confirmed by analyzing crystallographic phases, multiferroic properties of the fabricated CSMEN, influences of d.c. and a.c. magnetic fields on the CSMEN and cytotoxicity tests. The mathematical formalism to quantitatively describe the phenomena is also reported. The reported findings provide insights into the underlying MEEP mechanism and demonstrate the utility of CSMEN as an electric pulse-generating nano-probe in electroporation experiments with a potential application toward accurate and efficient targeted cell permeation.

Quantifying pulsed electric field-induced membrane nanoporation in single cells

Plasma membrane disruption can trigger a host of cellular activities. One commonly observed type of disruption is pore formation. Molecular dynamic (MD) simulations of simplified lipid membrane structures predict that controllably disrupting the membrane via nano-scale poration may be possible with nanosecond pulsed electric fields (nsPEF). Until recently, researchers hoping to verify this hypothesis experimentally have been limited to measuring the relatively slow process of fluorescent markers diffusing across the membrane, which is indirect evidence of nanoporation that could be channel-mediated. Leveraging recent advances in nonlinear optical microscopy, we elucidate the role of pulse parameters in nsPEF-induced membrane permeabilization in live cells. Unlike previous techniques, it is able to directly observe loss of membrane order at the onset of the pulse. We also develop a complementary theoretical model that relates increasing membrane permeabilization to membrane pore density. Due to the significantly improved spatial and temporal resolution possible with our imaging method, we are able to directly compare our experimental and theoretical results. Their agreement provides substantial evidence that nanoporation does occur and that its development is dictated by the electric field distribution.

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Mitigation of impedance changes due to electroporation therapy using bursts of high-frequency bipolar pulses

BACKGROUND

For electroporation-based therapies, accurate modeling of the electric field distribution within the target tissue is important for predicting the treatment volume. In response to conventional, unipolar pulses, the electrical impedance of a tissue varies as a function of the local electric field, leading to a redistribution of the field. These dynamic impedance changes, which depend on the tissue type and the applied electric field, need to be quantified a priori, making mathematical modeling complicated. Here, it is shown that the impedance changes during high-frequency, bipolar electroporation therapy are reduced, and the electric field distribution can be approximated using the analytical solution to Laplace's equation that is valid for a homogeneous medium of constant conductivity.

METHODS

Two methods were used to examine the agreement between the analytical solution to Laplace's equation and the electric fields generated by 100 µs unipolar pulses and bursts of 1 µs bipolar pulses. First, pulses were applied to potato tuber tissue while an infrared camera was used to monitor the temperature distribution in real-time as a corollary to the electric field distribution. The analytical solution was overlaid on the thermal images for a qualitative assessment of the electric fields. Second, potato ablations were performed and the lesion size was measured along the x- and y-axes. These values were compared to the analytical solution to quantify its ability to predict treatment outcomes. To analyze the dynamic impedance changes due to electroporation at different frequencies, electrical impedance measurements (1 Hz to 1 MHz) were made before and after the treatment of potato tissue.

RESULTS

For high-frequency bipolar burst treatment, the thermal images closely mirrored the constant electric field contours. The potato tissue lesions differed from the analytical solution by 39.7 ± 1.3 % (x-axis) and 6.87 ± 6.26 % (y-axis) for conventional unipolar pulses, and 15.46 ± 1.37 % (x-axis) and 3.63 ± 5.9 % (y-axis) for high- frequency bipolar pulses.

CONCLUSIONS

The electric field distributions due to high-frequency, bipolar electroporation pulses can be closely approximated with the homogeneous analytical solution. This paves way for modeling fields without prior characterization of non-linear tissue properties, and thereby simplifying electroporation procedures.

A transportable, inexpensive electroporator for in utero electroporation

Electroporation is a useful technique to study gene function during development but its broad application is hampered due to the expensive equipment needed. We describe the construction of a transportable, simple and inexpensive electroporator delivering square pulses with varying length and amplitude. The device was successfully used for in utero electroporation in mouse with a performance comparable to that of commercial products.

An engineered membrane to measure electroporation: effect of tethers and bioelectronic interface

This article reports on the construction and predictive models for a platform comprised of an engineered tethered membrane. The platform provides a controllable and physiologically relevant environment for the study of the electroporation process. The mixed self-assembled membrane is formed via a rapid solvent exchange technique. The membrane is tethered to the gold electrode and includes an ionic reservoir separating the membrane and gold surface. Above the membrane, there is an electrolyte solution, and a gold counterelectrode. A voltage is applied between the gold electrodes and the current measured. The current is dependent on the energy required to form aqueous pores and the conductance of each pore. A two-level predictive model, consisting of a macroscopic and a continuum model, is developed to relate the pore dynamics to the measured current. The macroscopic model consists of an equivalent circuit model of the tethered membrane, and asymptotic approximations to the Smoluchowski-Einstein equation of electroporation that is dependent on the pore conductance and the energy required to form aqueous pores. The continuum model is a generalized Poisson-Nernst-Planck (GPNP) system where an activity coefficient to account for steric effects of ions is added to the standard PNP system. The GPNP is used to evaluate the conductance of aqueous pores, and the electrical energy required to form the pores. As an outcome of the setup of the device and the two-level model, biologically important variables can be estimated from experimental measurements. To validate the accuracy of the two-level model, the predicted current is compared with experimentally measured current for different tethering densities.

Effects of nanosecond pulse electric fields on cellular elasticity

We investigated the effects of a single 60 nanosecond pulsed electric field (nsPEF) of low (15 kV/cm) and high (60 kV/cm) field strengths on cellular morphology and membrane elasticity in Jurkat cells using fluorescent microscopy and atomic force microscopy (AFM). We performed force displacement measurements on cells using AFM and calculated the Young's modulus for membrane elasticity. Differential effects were observed depending upon pulsing conditions. We found that a single nsPEF of low field strength did not induce any apparent cytoskeletal breakdown and had minor morphological changes. Interestingly, force measurements and calculation of Young's modulus showed a significant decrease in membrane elasticity. A single nsPEF of high field strength induced stark morphological changes due to disruption of the actin cytoskeleton and a marked decrease in elasticity likely caused by irreversible membrane damage. We suggest that the cellular morphology is mainly dependent on stabilization by the actin cytoskeleton, while the elasticity changes are partially dependent on the cytoskeletal integrity.

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 μ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 μs--evolves two pore subpopulations with a valley at ~5 nm, which separates the subpopulations that have peaks at ~1.5 and ~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 (<3 nm; ~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.

Ion transport into cells exposed to monopolar and bipolar nanosecond pulses

Experiments with CHO cells exposed to 60 and 300 ns pulsed electric fields with amplitudes in the range from several kV/cm to tens of kV/cm showed a decrease of the uptake of calcium ions by more than an order of magnitude when, immediately after a first pulse, a second one of opposite polarity was applied. This effect is assumed to be due to the reversal of the electrophoretic transport of ions through the electroporated membrane during the second phase of the bipolar pulse. This assumption, however, is only valid if electrophoresis is the dominant transport mechanism, rather than diffusion. Comparison of calculated calcium ion currents with experimental results showed that for nanosecond pulses, electrophoresis is at least as important as diffusion. By delaying the second pulse with respect to the first one, the effect of reverse electrophoresis is reduced. Consequently, separating nanosecond pulses of opposite polarity by up to approximately hundred microseconds allows us to vary the uptake of ions from very small values to those obtained with two pulses of the same polarity. The measured calcium ion uptake obtained with bipolar pulses also allowed us to determine the membrane pore recovery time. The calculated recovery time constants are on the order of 10 μs.

Cell electrofusion using nanosecond electric pulses

Electrofusion is an efficient method for fusing cells using short-duration high-voltage electric pulses. However, electrofusion yields are very low when fusion partner cells differ considerably in their size, since the extent of electroporation (consequently membrane fusogenic state) with conventionally used microsecond pulses depends proportionally on the cell radius. We here propose a new and innovative approach to fuse cells with shorter, nanosecond (ns) pulses. Using numerical calculations we demonstrate that ns pulses can induce selective electroporation of the contact areas between cells (i.e. the target areas), regardless of the cell size. We then confirm experimentally on B16-F1 and CHO cell lines that electrofusion of cells with either equal or different size by using ns pulses is indeed feasible. Based on our results we expect that ns pulses can improve fusion yields in electrofusion of cells with different size, such as myeloma cells and B lymphocytes in hybridoma technology.

Numerical study of lipid translocation driven by nanoporation due to multiple high-intensity, ultrashort electrical pulses

The dynamical translocation of lipids from one leaflet to another due to membrane permeabilization driven by nanosecond, high-intensity (>100kV/cm) electrical pulses has been probed. Our simulations show that lipid molecules can translocate by diffusion through water-filled nanopores which form following high voltage application. Our focus is on multiple pulsing, and such simulations are relevant to gauge the time duration over which nanopores might remain open, and facilitate continued lipid translocations and membrane transport. Our results are indicative of a N(½) scaling with pulse number for the pore radius. These results bode well for the use of pulse trains in biomedical applications, not only due to cumulative behaviors and in reducing electric intensities and pulsing hardware, but also due to the possibility of long-lived thermo-electric physics near the membrane, and the possibility for pore coalescence.

Conducting nanosponge electroporation for affordable and high-efficiency disinfection of bacteria and viruses in water

High-efficiency, affordable, and low energy water disinfection methods are in great need to prevent diarrheal illness, which is one of the top five leading causes of death over the world. Traditional water disinfection methods have drawbacks including carcinogenic disinfection byproducts formation, energy and time intensiveness, and pathogen recovery. Here, we report an innovative method that achieves high-efficiency water disinfection by introducing nanomaterial-assisted electroporation implemented by a conducting nanosponge filtration device. The use of one-dimensional (1D) nanomaterials allows electroporation to occur at only several volts, which is 2 to 3 orders of magnitude lower than that in traditional electroporation applications. The disinfection mechanism of electroporation prevents harmful byproduct formation and ensures a fast treatment speed of 15,000 L/(h·m(2)), which is equal to a contact time of 1 s. The conducting nanosponge made from low-cost polyurethane sponge coated with carbon nanotubes and silver nanowires ensures the device's affordability. This method achieves more than 6 log (99.9999%) removal of four model bacteria, including Escherichia coli, Salmonella enterica Typhimirium, Enterococcus faecalis, and Bacillus subtilis, and more than 2 log (99%) removal of one model virus, bacteriophage MS2, with a low energy consumption of only 100 J/L.

Nucleofection of expression vectors induces a robust interferon response and inhibition of cell proliferation

The interferon (IFN) response, induced as a side effect after transfection of nucleic acids into mammalian cells, is known but inadequately described. We followed the IFN response, the fate of cells, and the possible mechanisms leading to this response in NIH3T3 mouse fibroblasts after DNA nucleofection. The gateway destination vector, phGf, and its derivatives encoding toxic and non-toxic variants of the minor structural proteins of polyomaviruses, VP2 and VP3, were used. DNA vector sequences induced in cells the production of high levels of IFN and the upregulation of the IFN-inducible genes, Mx-1, STAT1, IRF1, and IRF7. The IFN response was not restricted to phGf-derived plasmids. In nucleofected cells, upregulation of the modified γ-histone 2A.X indicating DNA damage and inhibition of cell proliferation were also observed. Although 3T3 cells expressed the Toll-like receptor-9 (TLR9) and vectors used for nucleofection contained unmethylated CpGs, signaling leading to IFN induction was found to be TLR9 independent. However, the early activation of nuclear factor-kappa B suggested the participation of this transcription factor in IFN induction. Surprisingly, in contrast to nucleofection, transfection using a cationic polymer induced only a poor IFN response. Together, the results point to a strong side effect of nucleofection.

Inactivation of Pseudomonas putida by pulsed electric field treatment: a study on the correlation of treatment parameters and inactivation efficiency in the short-pulse range

An important issue for an economic application of the pulsed electric field treatment for bacterial decontamination of wastewater is the specific treatment energy needed for effective reduction of bacterial populations. The present experimental study performed in a field amplitude range of 40 > E > 200 kV/cm and for a suspension conductivity of 0.01 = κ(e) > 0.2 S/m focusses on the application of short pulses, 25 ns > T > 10 μs, of rectangular, bipolar and exponential shape and was made on Pseudomonas putida, which is a typical and widespread wastewater microorganism. The comparison of inactivation results with calculations of the temporal and azimuthal membrane charging dynamics using the model of Pauly and Schwan revealed that for efficient inactivation, membrane segments at the cell equator have to be charged quickly and to a sufficiently high value, on the order of 0.5 V. After fulfilling this basic condition by an appropriate choice of pulse field strength and duration, the log rate of inactivation for a given suspension conductivity of 0.2 S/m was found to be independent of the duration of individual pulses for constant treatment energy expenditure. Moreover, experimental results suggest that even pulse shape plays a minor role in inactivation efficiency. The variation of the suspension conductivity resulted in comparable inactivation performance of identical pulse parameters if the product of pulse duration and number of pulses was the same, i.e., required treatment energy can be linearly downscaled for lower conductivities, provided that pulse amplitude and duration are selected for entire membrane surface permeabilization.

Electroporation of intracellular liposomes using nanosecond electric pulses--a theoretical study

Nanosecond (ns) electric pulses of sufficient amplitude can provoke electroporation of intracellular organelles. This paper investigates whether such pulses could provide a method for controlled intracellular release of a content of small internalized artificial lipid vesicles (liposomes). To estimate the pulse parameters needed to selectively electroporate liposomes while keeping the plasma and nuclear membranes intact, we constructed a numerical model of a biological cell containing a nucleus and liposomes of different sizes (with radii from 50 to 500 nm), which were placed in various sites in the cytoplasm. Our results show that under physiological conditions selective electroporation is only possible for the largest liposomes and when using very short pulses (few ns). By increasing the liposome interior conductivity and/or decreasing the cytoplasmic conductivity, selective electroporation of even smaller liposomes could be achieved. The location of the liposomes inside the cell does not play a significant role, meaning that liposomes of similar size could all be electroporated simultaneously. Our results indicate the possibility of using ns pulse treatment for liposomal drug release.

Thresholds for phosphatidylserine externalization in Chinese hamster ovarian cells following exposure to nanosecond pulsed electrical fields (nsPEF)

High-amplitude, MV/m, nanosecond pulsed electric fields (nsPEF) have been hypothesized to cause nanoporation of the plasma membrane. Phosphatidylserine (PS) externalization has been observed on the outer leaflet of the membrane shortly after nsPEF exposure, suggesting local structural changes in the membrane. In this study, we utilized fluorescently-tagged Annexin V to observe the externalization of PS on the plasma membrane of isolated Chinese Hamster Ovary (CHO) cells following exposure to nsPEF. A series of experiments were performed to determine the dosimetric trends of PS expression caused by nsPEF as a function of pulse duration, τ, delivered field strength, E_D, and pulse number, n. To accurately estimate dose thresholds for cellular response, data were reduced to a set of binary responses and ED50s were estimated using Probit analysis. Probit analysis results revealed that PS externalization followed the non-linear trend of (τ*E_D²)⁻¹ for high amplitudes, but failed to predict low amplitude responses. A second set of experiments was performed to determine the nsPEF parameters necessary to cause observable calcium uptake, using cells preloaded with calcium green (CaGr), and membrane permeability, using FM1-43 dye. Calcium influx and FM1-43 uptake were found to always be observed at lower nsPEF exposure parameters compared to PS externalization. These findings suggest that multiple, higher amplitude and longer pulse exposures may generate pores of larger diameter enabling lateral diffusion of PS; whereas, smaller pores induced by fewer, lower amplitude and short pulse width exposures may only allow extracellular calcium and FM1-43 uptake.