Molecular dynamics simulations of ion conductance in field-stabilized nanoscale lipid electropores

Structure and Dynamics of Aqueous 2-Aminothiazole/NaCl Electrolytes at Electrified Interfaces

A computational study was performed to investigate the dynamics of aqueous electrolytes containing organic corrosion inhibitors near electrified interfaces by using the constant-charge model in classical molecular dynamics simulations. The results showed that when inhibitors form films at the interface, the surface charge of the electrode causes displacement of the molecules, referred to as electroporation. The hydrophobicity of the inhibitor molecules affects both the stability of the films and their recovery time. This study highlights the value of computational investigations of the dynamics within inhibitor films as a complementary approach to the traditional focus on inhibitor-substrate interactions, leading to deeper insights into the mechanisms of corrosion inhibition mechanisms.

Oral TNF-α siRNA delivery via milk-derived exosomes for effective treatment of inflammatory bowel disease

Oral administration facilitates the direct delivery of drugs to lesions within the small intestine and colon, making it an ideal approach for treating patients with inflammatory bowel disease. However, multiple physical barriers impede the delivery of oral RNA drugs through the gastrointestinal tract. Herein, we developed a novel oral siRNA delivery system that protects nucleic acids in extreme environments by employing exosomes derived from milk to encapsulate tumor necrosis factor-alpha (TNF-α) siRNA completely. The remarkable structural stability of milk-derived exosomes (M-Exos), as opposed to those from HEK293T cells, makes them exceptional siRNA carriers. Results demonstrate that milk exosomes loaded with TNF-α siRNA (M-Exo/siR) can effectively inhibit the expression of TNF-α-related inflammatory cytokines. Moreover, given that milk exosomes are composed of unique lipids with high bioavailability, orally administered M-Exo/siR effectively reach colonic tissues, leading to decreased TNF-α expression and successful alleviation of colitis symptoms in a dextran sulfate sodium-induced inflammatory bowel disease murine model. Hence, milk-derived exosomes carrying TNF-α siRNA can be effectively employed to treat inflammatory bowel disease. Indeed, using exosomes naturally derived from milk may shift the current paradigm of oral gene delivery, including siRNA.

Prediction of ionic conductivity from adiabatic heating in non-equilibrium molecular dynamics on various test systems

CONTEXT

The evaluation of ionic conductivity through atomistic modeling typically involves calculating diffusion coefficients, which often necessitates simulations spanning several hundreds of nanoseconds. This study introduces a less computationally demanding approach based on non-equilibrium molecular dynamics applicable to a wide range of systems.

METHOD

Ionic conductivity is determined by evaluating the Joule heating effect recorded during non-equilibrium molecular dynamics (NEMD) simulations. These simulations which involve applying a uniform electric field using classical force fields in LAMMPS are conducted within the MedeA software environment. The conductivity value for a specific temperature can thus be obtained from a single simulation together with an estimation of the associated uncertainty. Guidelines for selecting NEMD parameters such as electric field intensity and initial temperature are proposed to satisfy linear irreversible transport.

RESULTS

The protocol presented in this study is applied to four different types of systems, namely, (i) molten NaCl, (ii) NaCl and LiCl aqueous solutions, (iii) solution of ionic liquid with two solvents, and (iv) NaX zeolites in the anhydrous and hydrated states. The main advantages of the proposed protocol are simplicity of implementation (eliminating the need to store individual ion trajectories), reliability (low electric field, linear response, no perturbation of the equations of motion by a thermostat), and a wide range of applications. The estimated contribution of field-induced drift motion of ions to kinetic energy appears very low, justifying the use of standard kinetic energy in the method. For each system, the reported influence of temperature, ion concentration, solvent nature, or hydration is correctly predicted.

Studying the roles of salt ions in the pore initiation and closure stages in the biomembrane electroporation

Electroporation shows great potential in biology and biomedical applications. However, there is still a lack of reliable protocol for cell electroporation to achieve a high perforation efficiency due to the unclear influence mechanism of various factors, especially the salt ions in buffer solution. The tiny membrane structure of a cell and the electroporation scale make it difficult to monitor the electroporation process. In this study, we used both molecular dynamics (MD) simulation and experimental methods to explore the influence of salt ions on the electroporation process. Giant unilamellar vesicles (GUVs) were constructed as the model, and sodium chloride (NaCl) was selected as the representative salt ion in this study. The results show that the electroporation process follows lag-burst kinetics, where the lag period first appears after applying the electric field, followed by a rapid pore expansion. For the first time, we find that the salt ion plays opposite roles in different stages of the electroporation process. The accumulation of salt ions near the membrane surface provides an extra potential to promote the pore initiation, while the charge screening effect of the ions within the pore increases the line tension of the pore to induce the instability of the pore and lead to the closure. The GUV electroporation experiments obtain qualitatively consistent results with MD simulations. This work can provide guidance for the selection of parameters for cell electroporation process.

A molecular dynamics study of phospholipid membrane electroporation induced by bipolar pulses with different intervals

The mechanism of changes in cell electroporation (EP) during the intervals of bipolar pulses is still unclear, and few studies have investigated the effect of the intervals at the molecular level. In this study, EP induced by bipolar pulses (BP) with different intervals was investigated using all-atom molecular dynamics simulations. Firstly, EP was formed during the positive pulses of 2 ns and 0.5 V nm^-1, then the effects of various intervals of 0, 1, 5, and 10 ns on EP evolution were investigated, and the dynamic changes of different degrees of EP induced by the following negative pulses of 2 ns and 0.5 V nm^-1 were analyzed. The elimination effect of intervals was determined and it was related to the degrees of EP and the time of intervals. At the last moment of the intervals the phospholipid membrane was classified and quantitatively defined in three states according to the degrees of EP, namely, Resealing, Destabilizing and Retaining states. These states appeared due to the combined effect of both the positive pulse and the interval, and the states represent the degrees of EP which had different responses after applying the negative pulse. These results can improve our understanding of the fundamental mechanism of BP-induced EP.

Equivalence of Charge Imbalance and External Electric Fields during Free Energy Calculations of Membrane Electroporation

Electric fields across lipid membranes play important roles in physiology, medicine, and biotechnology, rationalizing the wide interest in modeling transmembrane potentials in molecular dynamics simulations. Transmembrane potentials have been implemented with external electric fields or by imposing charge imbalance between the two water compartments of a stacked double-membrane system. We compare the two methods in the context of membrane electroporation, which involves a large change of membrane structure and capacitance. We show that, given that Ewald electrostatics are defined with tinfoil boundary conditions, the two methods lead to (i) identical potentials of mean force (PMFs) of pore formation and expansion at various potentials, demonstrating that the two methods impose equivalent driving forces for large-scale transitions at membranes, and (ii) to identical polarization of water within thin water wires or open pores, suggesting that the two methods furthermore impose equivalent local electric fields. Without tinfoil boundary conditions, effects from external fields on pore formation are spuriously suppressed or even removed. Together, our study shows that both methods, external fields and charge imbalance, are well suitable for studying large-scale transitions of lipid membranes that involve changes of membrane capacitance. However, using charge imbalance is technically more challenging for maintaining a constant transmembrane potential since it requires updating of the charge imbalance as the membrane capacitance changes.

Assembly of transmembrane pores from mirror-image peptides

Tailored transmembrane alpha-helical pores with desired structural and functional versatility have promising applications in nanobiotechnology. Herein, we present a transmembrane pore DpPorA, based on the natural pore PorACj, built from D-amino acid α-helical peptides. Using single-channel current recordings, we show that DpPorA peptides self-assemble into uniform cation-selective pores in lipid membranes and exhibit properties distinct from their L-amino acid counterparts. DpPorA shows resistance to protease and acts as a functional nanopore sensor to detect cyclic sugars, polypeptides, and polymers. Fluorescence imaging reveals that DpPorA forms well-defined pores in giant unilamellar vesicles facilitating the transport of hydrophilic molecules. A second D-amino acid peptide based on the polysaccharide transporter Wza forms transient pores confirming sequence specificity in stable, functional pore formation. Finally, molecular dynamics simulations reveal the specific alpha-helical packing and surface charge conformation of the D-pores consistent with experimental observations. Our findings will aid the design of sophisticated pores for single-molecule sensing related technologies.

Alpha-helix nanopores have a range of potential applications and the inclusion of non-natural amino acids allows for modification. Here, the authors report on the creation of alpha-helix pores using D-amino acids and show the pores formed, have different properties to the L-counterparts and were resistant to proteases.

Pulsed Electric Fields Can Create Pores in the Voltage Sensors of Voltage-Gated Ion Channels

Pulsed electric fields are increasingly used in medicine to transiently increase the cell membrane permeability via electroporation to deliver therapeutic molecules into the cell. One type of event that contributes to this increase in membrane permeability is the formation of pores in the membrane lipid bilayer. However, electrophysiological measurements suggest that membrane proteins are affected as well, particularly voltage-gated ion channels (VGICs). The molecular mechanisms by which the electric field could affects these molecules remain unidentified. In this study, we used molecular dynamics simulations to unravel the molecular events that take place in different VGICs when exposing them to electric fields mimicking electroporation conditions. We show that electric fields can induce pores in the voltage-sensor domains (VSDs) of different VGICs and that these pores form more easily in some channels than in others. We demonstrate that poration is more likely in VSDs that are more hydrated and are electrostatically more favorable for the entry of ions. We further show that pores in VSDs can expand into so-called complex pores, which become stabilized by lipid headgroups. Our results suggest that such complex pores are considerably more stable than conventional lipid pores, and their formation can lead to severe unfolding of VSDs from the channel. We anticipate that such VSDs become dysfunctional and unable to respond to changes in transmembrane voltage, which is in agreement with previous electrophysiological measurements showing a decrease in the voltage-dependent transmembrane ionic currents after pulse treatment. Finally, we discuss the possibility of activation of VGICs by submicrosecond-duration pulses. Overall, our study reveals a new, to our knowledge, mechanism of electroporation through membranes containing VGICs.

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.

Controlling ionic conductivity through transprotein electropores in human aquaporin 4: a non-equilibrium molecular-dynamics study

Electroporation is a matter of intensive ongoing research interest, and a much-neglected topic in trans-membrane proteins, particularly in view of such promising potential applications in medicine and biotechnology. In particular, selected such novel and exciting applications are predicated on controlling ionic conductivity through electro-pores. Here, we scrutinise the mechanisms of ions' electric conductivity, by means of structural rearrangements, through quasi-stable electro-pores through human-AQP4 as a well-representative prototype of trans-membrane ionic conduction, achieving exquisite control over ionic permeability manipulated by the application of intense static electric fields.

Atomistic Simulations of Electroporation of Model Cell Membranes

Electroporation is a phenomenon that modifies the fundamental function of the cell since it perturbs transiently or permanently the integrity of its membrane. Today, this technique is applied in fields ranging from biology and biotechnology to medicine, e.g., for drug and gene delivery into cells, tumor therapy, etc., in which it made it to preclinical and clinical treatments. Experimentally, due to the complexity and heterogeneity of cell membranes, it is difficult to provide a description of the electroporation phenomenon in terms of atomically resolved structural and dynamical processes, a prerequisite to optimize its use. Atomistic modeling in general and molecular dynamics (MD) simulations in particular have proven to be an effective approach for providing such a level of detail. This chapter provides the reader with a comprehensive account of recent advances in using such a technique to complement conventional experimental approaches in characterizing several aspects of cell membranes electroporation.

Assessing the Structure and Stability of Transmembrane Oligomeric Intermediates of an α-Helical Toxin

Protein membrane interactions play an important role in our understanding of diverse phenomena ranging from membrane-assisted protein aggregation to oligomerization and folding. Pore-forming toxins (PFTs) are the primary vehicle for infection by several strains of bacteria. These proteins which are expressed in a water-soluble form (monomers) bind to the target membrane and conformationally transform (protomers) and self-assemble to form a multimer transmembrane pore complex through a process of oligomerization. On the basis of the structure of the transmembrane domains, PFTs are broadly classified into β or α toxins. In contrast to β-PFTs, the paucity of available crystal structures coupled with the amphipathic nature of the transmembrane domains has hindered our understanding of α-PFT pore formation. In this article, we use molecular dynamics (MD) simulations to examine the process of pore formation of the bacterial α-PFT, cytolysin A from Escherichia coli (ClyA) in lipid bilayer membranes. Using atomistic MD simulations ranging from 50 to 500 ns, we show that transmembrane oligomeric intermediates or "arcs" form stable proteolipidic complexes consisting of protein arcs with toroidal lipids lining the free edges. By creating initial conditions where the lipids are contained within the arcs, we study the dynamics of spontaneous lipid evacuation and toroidal edge formation. This process occurs on the time scale of tens of nanoseconds, suggesting that once protomers oligomerize, transmembrane arcs are rapidly stabilized to form functional water channels capable of leakage. Using umbrella sampling with a coarse-grained molecular model, we obtain the free energy of insertion of a single protomer into the membrane. A single inserted protomer has a stabilization free energy of -52.9 ± 1.2 kJ/mol and forms a stable transmembrane water channel capable of leakage. Our simulations reveal that arcs are stable and viable intermediates that can occur during the pore-formation pathway for ClyA.

Damage-free peripheral nerve stimulation by 12-ns pulsed electric field

Modern technologies enable deep tissue focusing of nanosecond pulsed electric field (nsPEF) for non-invasive nerve and muscle stimulation. However, it is not known if PEF orders of magnitude shorter than the activation time of voltage-gated sodium channels (VGSC) would evoke action potentials (APs). One plausible scenario requires the loss of membrane integrity (electroporation) and resulting depolarization as an intermediate step. We report, for the first time, that the excitation of a peripheral nerve can be accomplished by 12-ns PEF without electroporation. 12-ns stimuli at 4.1-11 kV (3.3-8.8 kV/cm) evoked APs similarly to conventional stimuli (100-250 μs, 1-5 V, 103-515 V/m), except for having higher selectivity for the faster nerve fibers. Nerves sustained repeated tetanic stimulations (50 Hz or 100 Hz for 1 min) alternately by 12-ns PEF and by conventional pulses. Such tetani caused a modest AP decrease, to a similar extent for both types of stimuli. Nerve refractory properties were not affected. The lack of cumulative damages even from tens of thousands of 12-ns stimuli and the similarities with the conventional stimulation prove VGSC activation by nsPEF without nerve membrane damage.

Electropermeabilization by uni- or bipolar nanosecond electric pulses: The impact of extracellular conductivity

Cellular effects caused by nanosecond electric pulses (nsEP) can be reduced by an electric field reversal, a phenomenon known as bipolar cancellation. The reason for this cancellation effect remains unknown. We hypothesized that assisted membrane discharge is the mechanism for bipolar cancellation. CHO-K1 cells bathed in high (16.1mS/cm; HCS) or low (1.8mS/cm; LCS) conductivity solutions were exposed to either one unipolar (300-ns) or two opposite polarity (300+300-ns; bipolar) nsEP (4-40kV/cm) with increasing interpulse intervals (0.1-50μs). Time-lapse YO-PRO-1 (YP) uptake revealed enhanced membrane permeabilization in LCS compared to HCS at all tested voltages. The time-dependence of bipolar cancellation was similar in both solutions, using either identical (22kV/cm) or isoeffective nsEP treatments (12 and 32kV/cm for LCS and HCS, respectively). However, cancellation was significantly stronger in LCS when the bipolar nsEP had no, or very short (<1μs), interpulse intervals. Finally, bipolar cancellation was still present with interpulse intervals as long as 50μs, beyond the time expected for membrane discharge. Our findings do not support assisted membrane discharge as the mechanism for bipolar cancellation. Instead they exemplify the sustained action of nsEP that can be reversed long after the initial stimulus.

Lipid Flip-Flop and Pore Nucleation on Zwitterionic Bilayers are Asymmetric under Ionic Imbalance

Lipid flip-flop and its associated transient pore formation are key thermodynamic properties of living cell membranes. However, there is a lack of understanding of whether ionic imbalance that exists ubiquitously across cell membranes affects lipid flip-flop and its associated functions. Potential of mean force calculations show that the free-energy barrier of lipid flip-flop on the extracellular leaflet reduces with the presence of ionic imbalance, whereas the barrier on the intracellular leaflet is generally not affected. The linear decrease of the activation energy of lipid flip-flop on the extracellular leaflet is consistent with the experimentally measured conductance-voltage relationship of zwitterionic lipid bilayers. This suggests: 1) lipid flip-flop has a directionality under physiological conditions and phospholipids accumulate at a rate on the order of 105 µm^-2 h^-1 on the cytoplasmic side of cell membranes; 2) ion permeation across a lipid membrane is moderated by lipid flip-flop; 3) the energy barrier of pore formation is aligned with the weaker leaflet that has a lower energy of lipid flip-flop. The asymmetry of lipid flip-flop and pore nucleation may have substantial implications for protein translocation, signaling, enzymatic activities, vesicle fusion, and transportation of biomolecules on cell membranes.

Neuronal excitation and permeabilization by 200-ns pulsed electric field: An optical membrane potential study with FluoVolt dye

Electric field pulses of nano- and picosecond duration are a novel modality for neurostimulation, activation of Ca²⁺ signaling, and tissue ablation. However it is not known how such brief pulses activate voltage-gated ion channels. We studied excitation and electroporation of hippocampal neurons by 200-ns pulsed electric field (nsPEF), by means of time-lapse imaging of the optical membrane potential (OMP) with FluoVolt dye. Electroporation abruptly shifted OMP to a more depolarized level, which was reached within <1ms. The OMP recovery started rapidly (τ=8-12ms) but gradually slowed down (to τ>10s), so cells remained above the resting OMP level for at least 20-30s. Activation of voltage-gated sodium channels (VGSC) enhanced the depolarizing effect of electroporation, resulting in an additional tetrodotoxin-sensitive OMP peak in 4-5ms after nsPEF. Omitting Ca²⁺ in the extracellular solution did not reduce the depolarization, suggesting no contribution of voltage-gated calcium channels (VGCC). In 40% of neurons, nsPEF triggered a single action potential (AP), with the median threshold of 3kV/cm (range: 1.9-4kV/cm); no APs could be evoked by stimuli below the electroporation threshold (1.5-1.9kV/cm). VGSC opening could already be detected in 0.5ms after nsPEF, which is too fast to be mediated by the depolarizing effect of electroporation. The overlap of electroporation and AP thresholds does not necessarily reflect the causal relation, but suggests a low potency of nsPEF, as compared to conventional electrostimulation, for VGSC activation and AP induction.

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).

A statistical analytical model for hydrophilic electropore characterization: a comparison study

A new procedure for lipid pore characterization, which extracts the pore geometry from a statistical analysis of the tridimensional surface of the pore, is compared with other methods in the literature.

Geometrical Characterization of an Electropore from Water Positional Fluctuations

We present here a new method for calculating the radius of a transmembrane pore in a phospholipid bilayer. To compare size-related properties of pores in bilayers of various compositions, generated and maintained under different physical and chemical conditions, reference metrics are needed. Operational metrics can be associated with some observed behavior. For example, pore size can be defined by the largest object that will pass through the length of the pore. The novelty of the present approach resides in the characterization of electropore geometry via a statistical approach, based on essential dynamics rules. We define the pore size geometrically with an algorithm for determining the pore radius. In particular, we extract the radius from the tri-dimensional surface of a defined pore region. The method is applied to a pore formed in a phospholipid bilayer by application of an external electric field. Although the details described here are specific for lipid pores in molecular dynamics simulations, the method can be generalized for any kind of pores for which appropriate structural information is available.

A Microdosimetric Study of Electropulsation on Multiple Realistically Shaped Cells: Effect of Neighbours

Over the past decades, the effects of ultrashort-pulsed electric fields have been used to investigate their action in many medical applications (e.g. cancer, gene electrotransfer, drug delivery, electrofusion). Promising aspects of these pulses has led to several in vitro and in vivo experiments to clarify their action. Since the basic mechanisms of these pulses have not yet been fully clarified, scientific interest has focused on the development of numerical models at different levels of complexity: atomic (molecular dynamic simulations), microscopic (microdosimetry) and macroscopic (dosimetry). The aim of this work is to demonstrate that, in order to predict results at the cellular level, an accurate microdosimetry model is needed using a realistic cell shape, and with their position and packaging (cell density) characterised inside the medium.

Dependence of Electroporation Detection Threshold on Cell Radius: An Explanation to Observations Non Compatible with Schwan's Equation Model

It is widely accepted that electroporation occurs when the cell transmembrane voltage induced by an external applied electric field reaches a threshold. Under this assumption, in order to trigger electroporation in a spherical cell, Schwan's equation leads to an inversely proportional relationship between the cell radius and the minimum magnitude of the applied electric field. And, indeed, several publications report experimental evidences of an inverse relationship between the cell size and the field required to achieve electroporation. However, this dependence is not always observed or is not as steep as predicted by Schwan's equation. The present numerical study attempts to explain these observations that do not fit Schwan's equation on the basis of the interplay between cell membrane conductivity, permeability, and transmembrane voltage. For that, a single cell in suspension was modeled and the electric field necessary to achieve electroporation with a single pulse was determined according to two effectiveness criteria: a specific permeabilization level, understood as the relative area occupied by the pores during the pulse, and a final intracellular concentration of a molecule due to uptake by diffusion after the pulse, during membrane resealing. The results indicate that plausible model parameters can lead to divergent dependencies of the electric field threshold on the cell radius. These divergent dependencies were obtained through both criteria and using two different permeabilization models. This suggests that the interplay between cell membrane conductivity, permeability, and transmembrane voltage might be the cause of results which are noncompatible with the Schwan's equation model.

Imaging the dynamics of individual electropores

Electroporation is a widely used technique to permeabilize cell membranes. Despite its prevalence, our understanding of the mechanism of voltage-mediated pore formation is incomplete; methods capable of visualizing the time-dependent behavior of individual electropores would help improve our understanding of this process. Here, using optical single-channel recording, we track multiple isolated electropores in real time in planar droplet interface bilayers. We observe individual, mobile defects that fluctuate in size, exhibiting a range of dynamic behaviors. We observe fast (25 s(-1)) and slow (2 s(-1)) components in the gating of small electropores, with no apparent dependence on the applied potential. Furthermore, we find that electropores form preferentially in the liquid disordered phase. Our observations are in general supportive of the hydrophilic toroidal pore model of electroporation, but also reveal additional complexity in the interactions, dynamics, and energetics of electropores.

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.

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.

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (Lα) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the Lα phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the Lα phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Picosecond and Terahertz Perturbation of Interfacial Water and Electropermeabilization of Biological Membranes

Non-thermal probing and stimulation with subnanosecond electric pulses and terahertz electromagnetic radiation may lead to new, minimally invasive diagnostic and therapeutic procedures and to methods for remote monitoring and analysis of biological systems, including plants, animals, and humans. To effectively engineer these still-emerging tools, we need an understanding of the biophysical mechanisms underlying the responses that have been reported to these novel stimuli. We show here that subnanosecond (≤500 ps) electric pulses induce action potentials in neurons and cause calcium transients in neuroblastoma-glioma hybrid cells, and we report complementary molecular dynamics simulations of phospholipid bilayers in electric fields in which membrane permeabilization occurs in less than 1 ns. Water dipoles in the interior of these model membranes respond in less than 1 ps to permeabilizing electric potentials by aligning in the direction of the field, and they re-orient at terahertz frequencies to field reversals. The mechanism for subnanosecond lipid electropore formation is similar to that observed on longer time scales-energy-minimizing intrusions of interfacial water into the membrane interior and subsequent reorganization of the bilayer into hydrophilic, conductive structures.

Electroporation-mediated gene delivery

Electroporation has been used extensively to transfer DNA to bacteria, yeast, and mammalian cells in culture for the past 30 years. Over this time, numerous advances have been made, from using fields to facilitate cell fusion, delivery of chemotherapeutic drugs to cells and tissues, and most importantly, gene and drug delivery in living tissues from rodents to man. Electroporation uses electrical fields to transiently destabilize the membrane allowing the entry of normally impermeable macromolecules into the cytoplasm. Surprisingly, at the appropriate field strengths, the application of these fields to tissues results in little, if any, damage or trauma. Indeed, electroporation has even been used successfully in human trials for gene delivery for the treatment of tumors and for vaccine development. Electroporation can lead to between 100 and 1000-fold increases in gene delivery and expression and can also increase both the distribution of cells taking up and expressing the DNA as well as the absolute amount of gene product per cell (likely due to increased delivery of plasmids into each cell). Effective electroporation depends on electric field parameters, electrode design, the tissues and cells being targeted, and the plasmids that are being transferred themselves. Most importantly, there is no single combination of these variables that leads to greatest efficacy in every situation; optimization is required in every new setting. Electroporation-mediated in vivo gene delivery has proven highly effective in vaccine production, transgene expression, enzyme replacement, and control of a variety of cancers. Almost any tissue can be targeted with electroporation, including muscle, skin, heart, liver, lung, and vasculature. This chapter will provide an overview of the theory of electroporation for the delivery of DNA both in individual cells and in tissues and its application for in vivo gene delivery in a number of animal models.

Electroporation in food processing and biorefinery

Electroporation is a method of treatment of plant tissue that due to its nonthermal nature enables preservation of the natural quality, colour and vitamin composition of food products. The range of processes where electroporation was shown to preserve quality, increase extract yield or optimize energy input into the process is overwhelming, though not exhausted; e.g. extraction of valuable compounds and juices, dehydration, cryopreservation, etc. Electroporation is--due to its antimicrobial action--a subject of research as one stage of the pasteurization or sterilization process, as well as a method of plant metabolism stimulation. This paper provides an overview of electroporation as applied to plant materials and electroporation applications in food processing, a quick summary of the basic technical aspects on the topic, and a brief discussion on perspectives for future research and development in the field. The paper is a review in the very broadest sense of the word, written with the purpose of orienting the interested newcomer to the field of electroporation applications in food technology towards the pertinent, highly relevant and more in-depth literature from the respective subdomains of electroporation research.

The importance of membrane defects-lessons from simulations

The defects and pores within lipid membranes are scientifically interesting and have a number of biological applications. Although lipid bilayers are extremely thin hydrophobic barriers, just ∼3 nm thick, they include diverse chemistry and have complex structures. Bilayers are soft and dynamic, and as a result, they can bend and deform in response to different stimuli by means of structural changes in their component lipids. Though defects occur within these structures, their transience and small size have made it difficult to characterize them. However, with recent advances in computer power and computational modeling techniques, researchers can now use simulations as a powerful tool to probe the mechanism and energies of defect and pore formation in a number of situations. In this Account, we present results from our detailed molecular dynamics computer simulations of hydrophilic pores and related defects in lipid bilayers at an atomistic level. Electroporation can be used to increase the permeability of cellular membranes, with potential therapeutic applications. Atomistic simulations of electroporation have illustrated the molecular details of this process, including the importance of water dipole interactions at the water-membrane interface. Characterization of the lipid-protein interactions provides an important tool for understanding transmembrane protein structure and thermodynamic stability. Atomistic simulations give a detailed picture of the free energies of model peptides and side chains in lipid membranes; the energetic cost of defect formation strongly influences the energies of interactions between lipids and polar and charged residues. Many antimicrobial peptides form hydrophilic pores in lipid membranes, killing bacteria or cancer cells. On the basis of simulation data, at least some of these peptides form defects and pores near the center of the bilayer, with a common disordered structure where hydrated headgroups form an approximately toroidal shape. The localization and trafficking of lipids supports general membrane structure and a number of important signaling cascades, such as those involving ceramide, diacylglycerol, and cholesterol. Atomistic simulations have determined the rates and free energies of lipid flip-flop. During the flip-flop of most phosphatidylcholine lipids, a hydrophilic pore forms when the headgroup moves near the center of the bilayer. Simulations have provided novel insight into many features of defects and pores in lipid membranes. Simulation data from very different systems and models show how water penetration and defect formation can determine the free energies of many membrane processes. Bilayers can deform and allow transient defects and pores when exposed to a diverse range of stimuli. Future work will explore many aspects of membrane defects with increased resolution and scope, including the study of more complex lipid mixtures, membrane domains, and large-scale membrane remodeling. Such studies will examine processes including vesicle budding and fusion, non-bilayer lipid phases, and interactions between lipid bilayers and other biomolecules. Simulations provide information that complements experimental studies, allowing microscopic insight into experimental observations and suggesting novel hypotheses and experiments. These studies should enable a deeper understanding of the role of lipid bilayers in cellular biology and support the development of future lipid-based biotechnology.

Effects of dimethyl sulfoxide on lipid membrane electroporation

Pores can be generated in lipid membranes by the application of an external electric field or by the addition of particular chemicals such as dimethyl sulfoxide (DMSO). Molecular dynamics (MD) has been shown to be a useful tool for unveiling many aspects of pore formation in lipid membranes in both situations. By means of MD simulations, we address the formation of electropores in cholesterol-containing lipid bilayers under the influence of DMSO. We show how a combination of physical and chemical mechanisms leads to more favorable conditions for generating membrane pores and, in particular, how the addition of DMSO to the medium significantly reduces the minimum electric field required to electroporate a lipid membrane. The strong alteration of membrane transversal properties and the energetic stabilization of the hydrophobic pore stage by DMSO provide the physicochemical mechanisms that explain this effect.

Cancellation of cellular responses to nanoelectroporation by reversing the stimulus polarity

Nanoelectroporation of biomembranes is an effect of high-voltage, nanosecond-duration electric pulses (nsEP). It occurs both in the plasma membrane and inside the cell, and nanoporated membranes are distinguished by ion-selective and potential-sensitive permeability. Here we report a novel phenomenon of bioeffects cancellation that puts nsEP cardinally apart from the conventional electroporation and electrostimulation by milli- and microsecond pulses. We compared the effects of 60- and 300-ns monopolar, nearly rectangular nsEP on intracellular Ca(2+) mobilization and cell survival with those of bipolar 60 + 60 and 300 + 300 ns pulses. For diverse endpoints, exposure conditions, pulse numbers (1-60), and amplitudes (15-60 kV/cm), the addition of the second phase cancelled the effects of the first phase. The overall effect of bipolar pulses was profoundly reduced, despite delivering twofold more energy. Cancellation also took place when two phases were separated into two independent nsEP of opposite polarities; it gradually tapered out as the interval between two nsEP increased, but was still present even at a 10-µs interval. The phenomenon of cancellation is unique for nsEP and has not been predicted by the equivalent circuit, transport lattice, and molecular dynamics models of electroporation. The existing paradigms of membrane permeabilization by nsEP will need to be modified. Here we discuss the possible involvement of the assisted membrane discharge, two-step oxidation of membrane phospholipids, and reverse transmembrane ion transport mechanisms. Cancellation impacts nsEP applications in cancer therapy, electrostimulation, and biotechnology, and provides new insights into effects of more complex waveforms, including pulsed electromagnetic emissions.