A statistical model for multidimensional irreversible electroporation cell death in tissue

The additive effect of 17β-estradiol on the modulation of electrochemotherapy with calcium ions or cisplatin in human clear carcinoma cells

Calcium electroporation (CaEP) is an efficient approach for ovarian cancer treatment. It causes cell death by introducing elevated levels of calcium into cells. In this work, the research focused on two types of cell lines: CHO-K1, representing normal ovary cells, and OvBH-1, representing ovarian clear carcinoma cells. Those cell lines exhibited distinct reactions to calcium electroporation (CaEP). Also, we have evaluated the effects of 17β-estradiol following CaEP and electrochemotherapy (ECT) with cisplatin (CPP). The combination of ECT with CPP and CaEP with prior E2 preincubation resulted in approximately 23.55 % and 39 % decreases in cell survival compared to the control cells (exposed to CPP and CaCl2 alone) for ovarian cancer cells. The obtained results showed that ovarian cancer cells preincubated with 17β-estradiol after exposure to pulsed electric fields undergo primary necrosis. Additionally, preincubation of ovarian cancer cells with 17β-estradiol can significantly improve the effectiveness of both chemotherapy and electrochemotherapy involving cisplatin and calcium chloride.

In silico study about the influence of electroporation parameters on the cellular internalization, spatial uniformity, and cytotoxic effects of chemotherapeutic drugs using the Method of Fundamental Solutions

Reversible electroporation is a suitable technique to aid the internalization of medicaments in cancer tissues without inducing permanent cellular damage, allowing the enhancement of cytotoxic effects without incurring in electric-driven necrotic or apoptotic processes by the presence of non-reversible aqueous pores. An adequate selection of electroporation parameters acquires relevance to reach these goals and avoid opposite effects. This work applies the Method of Fundamental Solutions (MFS) for drug transport simulations in electroporated cancer tissues, using a continuum tumor cord approach and considering both electro-permeabilization and vasoconstriction effects. The MFS algorithm is validated with published results, obtaining satisfactory accuracy and convergence. Then, MFS simulations are executed to study the influence of electric field magnitude [Formula: see text], number of electroporation treatments [Formula: see text], and electroporation time [Formula: see text] on three assessment parameters of electrochemotherapy: the internationalization efficacy accounting for the ability of the therapy to introduce moles into viable cells, cell-kill capacity indicating the faculty to reduce the survival fraction of cancer cells, and distribution uniformity specifying the competence to supply drug homogeneously through the whole tissue domain. According to numerical results, when [Formula: see text] is the reversibility threshold, a positive influence on the first two parameters is only possible once specific values of [Formula: see text] and [Formula: see text] have been exceeded; when [Formula: see text] is just the irreversibility threshold, any combination of [Formula: see text] and [Formula: see text] is beneficial. On the other hand, the drug distribution uniformity is always adversely affected by the application of electric pulses, being this more noticeable as [Formula: see text], [Formula: see text], and [Formula: see text] increases.

Electroporation-Based Biopsy Treatment Planning with Numerical Models and Tissue Phantoms

Molecular sampling with vacuum-assisted tissue electroporation is a novel, minimally invasive method for molecular profiling of solid lesions. In this paper, we report on the design of the battery-powered pulsed electric field generator and electrode configuration for an electroporation-based molecular sampling device for skin cancer diagnostics. Using numerical models of skin electroporation corroborated by the potato tissue phantom model, we show that the electroporated tissue volume, which is the maximum volume for biomarker sampling, strongly depends on the electrode's geometry, needle electrode skin penetration depths, and the applied pulsed electric field protocol. In addition, using excised human basal cell carcinoma (BCC) tissues, we show that diffusion of proteins out of human BCC tissues into water strongly depends on the strength of the applied electric field and on the time after the field application. The developed numerical simulations, confirmed by experiments in potato tissue phantoms and excised human cancer lesions, provide essential tools for the development of electroporation-based molecular markers sampling devices for personalized skin cancer diagnostics.

Temperature Distribution on Classical Two Needles IRE Setup Versus a Single Needle Prototype

OBJECTIVES

Irreversible Electroporation (IRE) is a non-thermal minimally invasive cancer therapy used in the treatment of liver tumors. However, the therapy entails an electrical current flux which can be high enough to cause a noticeable temperature increase. Therefore, the analysis of the heat distribution is important: during any IRE treatment, the target area is intended to be treated with non-thermal effects, where existing thermal effects should not damage nearby sensitive structures. This article aims to compare the established two parallel needles electrode setup, used by FDA-approved electroporation delivering devices, to a single needle, multiple electrode prototype design.

METHODS

Levels and distributions of the temperature at different distances from the applicators during an IRE liver treatment were investigated. The prototype results were collated with already published in-vivo data. All electrode configurations were analyzed numerically in COMSOL Multiphysics for different pulse protocols.

RESULTS

The extension of coagulation necrosis predicted by the model matched available in-vivo data. While the maximum average temperature during pulsation was higher for the prototype (74 °C) than for the two-needle IRE setup (57 °C), the thickness of the coagulation necrosis around the conductive electrodes was in the same range for both configurations. However, the location differed completely: the necrosis engendered by the prototype was located inside the tumor, while the two-needle IRE setup created necrosis outside the tumor, potentially closer to sensitive structures.

CONCLUSION

The results highlighted the importance of heat distribution analysis for the design of new IRE needles as well as for IRE treatment planning. Proper analysis ensures that the non-thermal effects are maximized while minimizing any potential thermal damage to surrounding sensitive structures.

Influence of electric field, blood velocity, and pharmacokinetics on electrochemotherapy efficiency

The convective delivery of chemotherapeutic drugs in cancerous tissues is directly proportional to the blood perfusion rate, which in turns can be transiently reduced by the application of high-voltage and short-duration electric pulses due to vessel vasoconstriction. However, electric pulses can also increase vessel wall and cell membrane permeabilities, boosting the extravasation and cell internalization of drug. These opposite effects, as well as possible adverse impacts on the viability of tissues and endothelial cells, suggest the importance of conducting in silico studies about the influence of physical parameters involved in electric-mediated drug transport. In the present work, the global method of approximate particular solutions for axisymmetric domains, together with two solution schemes (Gauss-Seidel iterative and linearization+successive over-relaxation), is applied for the simulation of drug transport in electroporated cancer tissues, using a continuum tumor cord approach and considering both the electropermeabilization and vasoconstriction phenomena. The developed global method of approximate particular solutions algorithm is validated with numerical and experimental results previously published, obtaining a satisfactory accuracy and convergence. Then, a parametric study about the influence of electric field magnitude and inlet blood velocity on the internalization efficacy, drug distribution uniformity, and cell-kill capacity of the treatment, as expressed by the number of internalized moles into viable cells, homogeneity of exposure to bound intracellular drug, and cell survival fraction, respectively, is analyzed for three pharmacokinetic profiles, namely one-short tri-exponential, mono-exponential, and uniform. According to numerical results, the trade-off between vasoconstriction and electropermeabilization effects and, consequently, the influence of electric field magnitude and inlet blood velocity on the assessment parameters considered here (efficacy, uniformity, and cell-kill capacity) is different for each pharmacokinetic profile deemed.

Irreversible Electroporation Combined With Dendritic Cell-based Vaccines for the Treatment of Osteosarcoma

Osteosarcoma is the most common primary bone malignancy, and surgical resection combined with neoadjuvant chemotherapy is the gold-standard treatment for affected patients. Although the overall survival rates for patients with osteosarcoma currently range from 60% to 70%, outcomes remain disappointing for patients with recurrent, metastatic, or unresectable disease. Irreversible electroporation (IRE) is a novel ablation technique with the potential to elicit an immune response in solid tumors. Dendritic cell (DC)-based tumor vaccines have shown promising therapeutic efficacy in preclinical studies focused on osteosarcoma; however, only limited therapeutic efficacy has been observed in clinical trials. Thus, there is considerable potential therapeutic value in developing combination osteosarcoma treatments that involve IRE and DC-based tumor vaccines. In this review, we discuss recent advances in preclinical and clinical DC-based immunotherapies, as well as potential combinations of DC-based vaccines and IRE, that may improve therapeutic outcomes for patients with osteosarcoma.

Influence of electric pulse characteristics on the cellular internalization of chemotherapeutic drugs and cell survival fraction in electroporated and vasoconstricted cancer tissues using boundary element techniques

Electroporation has emerged as a suitable technique to induce the pore formation in the cell membrane of cancer tissues, facilitating the cellular internalization of chemotherapeutic drugs. An adequate selection of the electric pulse characteristics is crucial to guarantee the efficiency of this technique, minimizing the adverse effects. In the present work, the dual reciprocity boundary element method (DR-BEM) is applied for the simulation of drug transport in the extracellular and intracellular space of cancer tissues subjected to the application of controlled electric pulses, using a continuum tumour cord approach, and considering both the electro-permeabilization and vasoconstriction phenomena. The developed DR-BEM algorithm is validated with numerical and experimental results previously published, obtaining a satisfactory accuracy and convergence. Using the DR-BEM code, a study about the influence of the magnitude of electric field (E) and pulse spacing (dpulses) on the time behavior and spatial distribution of the internalized drug, as well as on the cell survival fraction, is carried out. In general, the change of drug concentration, drug exposure and cell survival fraction with the parameters E and dpulses is ruled by two important factors: the balance between the electro-permeabilization and vasoconstriction phenomena, and the relative importance of the sources of cell death (electric pulses and drug cytotoxicity); these two factors, in turn, significantly depend on the reversible and irreversible thresholds considered for the electric field.

Quantitative analysis of contribution of mild and moderate hyperthermia to thermal ablation and sensitization of irreversible electroporation of pancreatic cancer cells

INTRODUCTION

Irreversible electroporation (IRE) is an ablation modality that applies short, high-voltage electric pulses to unresectable cancers. Although considered a non-thermal technique, temperatures do increase during IRE. This temperature rise sensitizes tumor cells for electroporation as well as inducing partial direct thermal ablation.

AIM

To evaluate the extent to which mild and moderate hyperthermia enhance electroporation effects, and to establish and validate in a pilot study cell viability models (CVM) as function of both electroporation parameters and temperature in a relevant pancreatic cancer cell line.

METHODS

Several IRE-protocols were applied at different well-controlled temperature levels (37 °C ≤ T ≤ 46 °C) to evaluate temperature dependent cell viability at enhanced temperatures in comparison to cell viability at T = 37 °C. A realistic sigmoid CVM function was used based on thermal damage probability with Arrhenius Equation and cumulative equivalent minutes at 43 °C (CEM43°C) as arguments, and fitted to the experimental data using "Non-linear-least-squares"-analysis.

RESULTS

Mild (40 °C) and moderate (46 °C) hyperthermic temperatures boosted cell ablation with up to 30% and 95%, respectively, mainly around the IRE threshold E_th,50% electric-field strength that results in 50% cell viability. The CVM was successfully fitted to the experimental data.

CONCLUSION

Both mild- and moderate hyperthermia significantly boost the electroporation effect at electric-field strengths neighboring E_th,50%. Inclusion of temperature in the newly developed CVM correctly predicted both temperature-dependent cell viability and thermal ablation for pancreatic cancer cells exposed to a relevant range of electric-field strengths/pulse parameters and mild moderate hyperthermic temperatures.

Survival model database of human digestive system cells exposed to electroporation pulses: An in vitro and in silico study

Background and objectives

This study aimed to establish a mathematical survival model database containing cell-specific coefficients from human digestive system cells exposed to electroporation pulses (EPs).

Materials and methods

A total of 20 types of human digestive system cell lines were selected to investigate the effect of EPs on cell viability. Cell viability was measured after exposure to various pulse settings, and a cell survival model was established using the Peleg-Fermi model. Next, the cell-specific coefficients of each cell line were determined.

Results

Cell viability tended to decrease when exposed to stronger electric field strength (EFS), longer pulse duration, and more pulse number, but the decreasing tendency varied among different cell lines. When exposed to a lower EFS (<1,000 V/cm), only a slight decrease in cell viability occurred. All cell lines showed a similar tendency: the extent of electrical injury (EI) increased with the increase in pulse number and duration. However, there existed differences in heat sensitivity among organs.

Conclusions

This database can be used for the application of electroporation-based treatment (EBT) in the digestive system to predict cell survival and tissue injury distribution during the treatment.

Automated irreversible electroporated region prediction using deep neural network, a preliminary study for treatment planning

The primary purpose of cancer treatment with irreversible electroporation (IRE) is to maximize tumor damage and minimize surrounding healthy tissue damage. Finite element analysis is one of the popular ways to calculate electric field and cell kill probability in IRE. However, this method also has limitations. This paper will focus on using a deep neural network (DNN) in IRE to predict irreversible electroporated regions for treatment planning purposes. COMSOL Multiphysics was used to simulate the IRE. The electric conductivity change during IRE was considered to create accurate data sets of electric field distribution and cell kill probability distributions. We used eight pulses with a pulse width of 100 μs, frequency of 1 Hz, and different voltages. To create masks for DNN training, a 90% cell kill probability contour was used. After data set creation, U-Net architecture was trained to predict irreversible electroporated regions. In this study, the average U-Net DICE coefficient on test data was 0.96. Also, the average accuracy of U-Net for predicting irreversible electroporated regions was 0.97. As far as we are aware, this is the first time that U-Net was used to predict an irreversible electroporated region in IRE. The present study provides significant evidence for U-Net's use for predicting an irreversible electroporated region in treatment planning.

Modeling methods for treatment planning in overlapping electroporation treatments

OBJECTIVE

Irreversible electroporation (IRE) is a non thermal tissue ablation therapy which is induced by applying high voltage waveforms across electrode pairs. When multiple electrode pairs are sequentially used, the treatment volume (TV) is typically computed as the geometric union of the TVs of individual pairs. However, this method neglects that some regions are exposed to overlapping treatments. Recently, a model describing cell survival probability was introduced which effectively predicted TV with overlapping fields in vivo. However, treatment overlap has yet to be quantified. This study characterizes TV overlap in a controlled in vitro setup with the two existing methods which are compared to an adapted logistic model proposed here.

METHODS

CHO cells were immobilized in agarose gel. Initially, we characterized the electric field threshold and the cell survival probability for overlapping treatments. Subsequently, we created a 2D setup where we compared and validated the accuracy of the different methods in predicting the TV.

RESULTS

Overlap can reduce the electric field threshold required to induce cell death, particularly for treatments with low pulse number. However, it does not have a major impact on TV in the models assayed here, and all the studied methods predict TV with similar accuracy.

CONCLUSION

Treatment overlap has a minor influence in the TV for typical protocols found in IRE therapies.

SIGNIFICANCE

This study provides evidence that the modeling method used in most pre-clinical and clinical studies seems adequate.

Treatment Planning Optimization in Irreversible Electroporation for Complete Ablation of Variously Sized Cervical Tumors: A Numerical Study

Irreversible electroporation (IRE), a relatively new energy-based tumor ablation technology, has shown itself in the last decade to be able to safely ablate tumors with favorable clinical outcomes, yet little work has been done on optimizing the IRE protocol to variously sized tumors. Incomplete tumor ablation has been shown to be the main reason leading to the local recurrence and thus treatment failure. The goal of this study was to develop a general optimization approach to optimize the IRE protocol for cervical tumors in different sizes, while minimizing the damage to normal tissues. This kind of approach can lay a foundation for future personalized treatment of IRE. First, a statistical IRE cervical tumor death model was built using previous data in our group. Then, a multi-objective optimization problem model was built, in which the decision variables are five IRE-setting parameters, namely, the pulse strength (U), the length of active tip (H), the number of pulses delivered in one round between a pair of electrodes (A), the distance between electrodes (D), and the number of electrodes (N). The domains of the decision variables were determined based on the clinical experience. Finally, the problem model was solved by using nondominated sorting genetic algorithms II (NSGA-II) algorithm to give respective optimal protocol for three sizes of cervical tumors. Every protocol was assessed by the evaluation criterion established in the study to show the efficacy in a more straightforward way. The results of the study demonstrate this approach can theoretically provide the optimal IRE protocol for different sizes of tumors and may be generalizable to other types, sizes, and locations of tumors.

AC Pulsed Field Ablation Is Feasible and Safe in Atrial and Ventricular Settings: A Proof-of-Concept Chronic Animal Study

Introduction

Pulsed field ablation (PFA) exploits the delivery of short high-voltage shocks to induce cells death via irreversible electroporation. The therapy offers a potential paradigm shift for catheter ablation of cardiac arrhythmia. We designed an AC-burst generator and therapeutic strategy, based on the existing knowledge between efficacy and safety among different pulses. We performed a proof-of-concept chronic animal trial to test the feasibility and safety of our method and technology.

Methods

We employed 6 female swine - weight 53.75 ± 4.77 kg - in this study. With fluoroscopic and electroanatomical mapping assistance, we performed ECG-gated AC-PFA in the following settings: in the left atrium with a decapolar loop catheter with electrodes connected in bipolar fashion; across the interventricular septum applying energy between the distal electrodes of two tip catheters. After procedure and 4-week follow-up, the animals were euthanized, and the hearts were inspected for tissue changes and characterized. We perform finite element method simulation of our AC-PFA scenarios to corroborate our method and better interpret our findings.

Results

We applied square, 50% duty cycle, AC bursts of 100 μs duration, 100 kHz internal frequency, 900 V for 60 pulses in the atrium and 1500 V for 120 pulses in the septum. The inter-burst interval was determined by the native heart rhythm - 69 ± 9 bpm. Acute changes in the atrial and ventricular electrograms were immediately visible at the sites of AC-PFA - signals were elongated and reduced in amplitude (p < 0.0001) and tissue impedance dropped (p = 0.011). No adverse event (e.g., esophageal temperature rises or gas bubble streams) was observed - while twitching was avoided by addition of electrosurgical return electrodes. The implemented numerical simulations confirmed the non-thermal nature of our AC-PFA and provided specific information on the estimated treated area and need of pulse trains. The postmortem chest inspection showed no peripheral damage, but epicardial and endocardial discolorations at sites of ablation. T1-weighted scans revealed specific tissue changes in atria and ventricles, confirmed to be fibrotic scars via trichrome staining. We found isolated, transmural and continuous scars. A surviving cardiomyocyte core was visible in basal ventricular lesions.

Conclusion

We proved that our method and technology of AC-PFA is feasible and safe for atrial and ventricular myocardial ablation, supporting their systematic investigation into effectiveness evaluation for the treatment of cardiac arrhythmia. Further optimization, with energy titration or longer follow-up, is required for a robust atrial and ventricular AC-PFA.

The dose-response characteristics of four NTCP models: using a novel CT-based radiomic method to quantify radiation-induced lung density changes

Multiple competing normal tissue complication probability (NTCP) models have been proposed for predicting symptomatic radiation-induced lung injury in human. In this paper we tested the efficacy of four common NTCP models applied quantitatively to sub-clinical X-ray computed tomography (CT)-density changes in the lung following radiotherapy. Radiotherapy planning datasets and follow-up chest CTs were obtained in eight patients treated for targets within the lung or hilar region. Image pixel-wise radiation dose exposure versus change in observable CT Hounsfield units was recorded for early (2-5 months) and late (6-9 months) time-points. Four NTCP models, Lyman, Logistic, Weibull and Poisson, were fit to the population data. The quality of fits was assessed by five statistical criteria. All four models fit the data significantly (p < 0.05) well at early, late and cumulative time points. The Lyman model fitted best for early effects while the Weibull Model fitted best for late effects. No significant difference was found between the fits of the models and with respect to parameters D₅₀ and γ₅₀. The D₅₀ estimates were more robust than γ₅₀ to image registration error. For analyzing population-based sub-clinical CT pixel intensity-based dose response, all four models performed well.

A numerical study on the effect of conductivity change in cell kill distribution in irreversible electroporation

Introduction: irreversible electroporation (IRE) is a tissue ablation technique and physical process used to kill the undesirable cells. In the IRE process by mathematical modelling we can calculate the cell kill probability and distribution inside the tissue. The purpose of the study is to determine the influence of electric conductivity change in the IRE process into the cell kill probability and distribution.

Methods: cell death probability and electric conductivity were calculated with COMSOL Multiphysics software package. 8 pulses with a frequency of 1 Hz, pulse width of 100 µs and electric field intensity from 1000 to 3000 V/Cm with steps of 500 V/Cm used as electric pulses.

Results: significantly, the electrical conductivity of tissue will increase during the time of pulse delivery. According to our results, electrical conductivity increased with an electric field intensity of pulses. By considering the effect of conductivity change on cell kill probability, the cell kill probability and distribution will change.

Conclusion: we believe that considering the impact of electric conductivity change on the cell kill probability will improve the accuracy of treatment outcome in the clinic for treatment with IRE.

Development of a Multi-Pulse Conductivity Model for Liver Tissue Treated With Pulsed Electric Fields

Pulsed electric field treatment modalities typically utilize multiple pulses to permeabilize biological tissue. This electroporation process induces conductivity changes in the tissue, which are indicative of the extent of electroporation. In this study, we characterized the electroporation-induced conductivity changes using all treatment pulses instead of solely the first pulse as in conventional conductivity models. Rabbit liver tissue was employed to study the tissue conductivity changes caused by multiple, 100 μs pulses delivered through flat plate electrodes. Voltage and current data were recorded during treatment and used to calculate the tissue conductivity during the entire pulsing process. Temperature data were also recorded to quantify the contribution of Joule heating to the conductivity according to the tissue temperature coefficient. By fitting all these data to a modified Heaviside function, where the two turning points (E₀, E₁) and the increase factor (A) are the main parameters, we calculated the conductivity as a function of the electric field (E), where the parameters of the Heaviside function (A and E₀) were functions of pulse number (N). With the resulting multi-factor conductivity model, a numerical electroporation simulation can predict the electrical current for multiple pulses more accurately than existing conductivity models. Moreover, the saturating behavior caused by electroporation can be explained by the saturation trends of the increase factor A in this model. The conductivity change induced by electroporation has a significant increase at about the first 30 pulses, then tends to saturate at 0.465 S/m. The proposed conductivity model can simulate the electroporation process more accurately than the conventional conductivity model. The electric field distribution computed using this model is essential for treatment planning in biomedical applications utilizing multiple pulsed electric fields, and the method proposed here, relating the pulse number to the conductivity through the variables in the Heaviside function, may be adapted to investigate the effect of other parameters, like pulse frequency and pulse width, on electroporation.

Simulation of the influence of voltage level and pulse spacing on the efficiency, aggressiveness and uniformity of the electroporation process in tissues using meshless techniques

Electroporation is a widely used method consisting of application of high-voltage, short-duration electric pulses to increase cell membrane permeability, allowing cellular internalization of medications. In this work, the influence of two primary parameters, voltage level (V) and pulse spacing (N), on electroporation efficiency, uniformity and aggressiveness, as quantified by the total mass transport to viable cells, intracellular concentration gradients and an aggressiveness factor introduced here, is studied by means of numerical simulations of drug transport in electroporated tissues. The global method of approximate particular solutions (Global MAPS) is used to solve the governing equations, together with domain scaling, singular value decomposition and smoothing algorithms, to address the ill-conditioning of the final system and suppress small scale oscillations. The accuracy of Global MAPS is evaluated by comparing the initial extracellular concentration, C_e , and final intracellular concentration, C_i , with previous finite volume method results, obtaining similar behavior of C_e and C_i along the tissue domain, with some differences for C_i in high-gradient zones. According to the Global MAPS results, the influence of V and N on C_i is only significant over a certain range, within which the largest drug transport to viable cells occurs. In general, both electroporation efficiency and aggressiveness change in nonuniform manner with V and decrease with N, whereas the electroporation uniformity decreases as V increases and N decreases. The contour plots obtained here can be considered useful tools to compare electroporation-based treatments in terms of their efficiency, aggressiveness and uniformity, assisting in the selection of a suitable treatment plan for cancer.

Genome mutation after the introduction of the gene editing by electroporation of Cas9 protein (GEEP) system into bovine putative zygotes

The present study was designed to investigate the effects of voltage strength on embryonic developmental rate and mutation efficiency in bovine putative zygotes during electroporation with the CRISPR/Cas9 system to target the MSTN gene at different time points after insemination. Results showed that there was no significant interaction between electroporation time and voltage strength on the embryonic cleavage and blastocyst formation rates. However, increasing the voltage strength to 20 V/mm to electroporate the zygotes at 10 h after the start of insemination yielded significantly lower blastocyst formation rates (P < 0.05) than those of the 10-V/mm electroporated zygotes. Mutation efficiency was then assessed in individual blastocysts by DNA sequence analysis of the target sites in the MSTN gene. A positive correlation between mutation rate and voltage strength was observed. The mutation efficiency in mutant blastocysts was significantly higher in the zygotes electroporated with 20 V/mm at 10 h after the start of insemination (P < 0.05) than in the zygotes electroporated at 15 h, irrespective of the voltage strength. We also noted that a certain number of blastocysts from zygotes that were electroporated with more than 15 V/mm at 10 h (4.8-16.7%) and 20 V/mm at 15 h (4.8%) were biallelic mutants. Our results suggest that the voltage strength during electroporation as well as electroporation time certainly have effects on the embryonic developmental rate and mutation efficiency in bovine putative zygotes.

Ascending Defibrillation Waveform Significantly Reduces Myocardial Morphological Damage and Injury Current

OBJECTIVES

This study tested the hypothesis that a biphasic defibrillation waveform with an ascending first phase (ASC) causes less myocardial damage by pathology and injury current than a standard biphasic truncated exponential (BTE) waveform in a swine model.

BACKGROUND

Although lifesaving, defibrillation shocks have significant iatrogenic effects that reduce their benefit for patient survival.

METHODS

An ASC waveform with an 8-ms linear ramp followed by an additional positive 0.5-ms decaying portion with amplitudes of 20 J (ASC 20J) and 25 J (ASC 25J) was used. The control was a 25-J BTE conventional waveform (BTE 25J) RESULTS: The ASC 20J and ASC 25J shocks were both successful in 6 of 6 pigs, but the BTE 25J was successful in only 6 of 14 pigs (p < 0.05). Post-shock ST-segment elevation (injury current) in the right ventricular electrode was significantly greater with BTE 25J than with ASC 20J and ASC 25J. With a blinded pathology reading, hemorrhage, inflammation, thrombi, and necrosis 24 h post-shock were significantly greater with BTE 25J than with ASC 20J and ASC 25J. Troponin levels were also markedly lower at 3, 4, 5, and 6 h post-shock.

CONCLUSIONS

Defibrillation shocks cause electrophysiological, histological, and biochemical signs of myocardial damage and necrosis. These signs of damage are markedly less for an ASC waveform than for a conventional BTE waveform.

Selective Inactivation of Pseudomonas aeruginosa and Staphylococcus epidermidis with Pulsed Electric Fields and Antibiotics

Objective: Increasing numbers of multidrug-resistant bacteria make many antibiotics ineffective; therefore, new approaches to combat microbial infections are needed. In addition, antibiotics are not selective-they kill pathogenic organisms as well as organisms that could positively contribute to wound healing (bio flora). Approach: Here we report on selective inactivation of Pseudomonas aeruginosa and Staphylococcus epidermidis, potential pathogens involved in wound infections with pulsed electric fields (PEFs) and antibiotics (mix of penicillin, streptomycin, and nystatin). Results: Using a Taguchi experimental design in vitro, we found that, under similar electric field strengths, the pulse duration is the most important parameter for P. aeruginosa inactivation, followed by the number of pulses and pulse frequency. P. aeruginosa, a potential severe pathogen, is more sensitive than the less pathogenic S. epidermidis to PEF (alone or in combination with antibiotics). Applying 200 pulses with a duration of 60 μs at 2.8 Hz, the minimum electric fields of 308.8 ± 28.3 and 378.4 ± 12.9 V/mm were required to inactive P. aeruginosa and S. epidermidis, respectively. Addition of antibiotics reduced the threshold for minimum electric fields required to inactivate the bacteria. Innovation: This study provides essential information, such as critical electric field parameters for bacteria inactivation, required for developing in vivo treatment and clinical protocols for using PEF for wound healing. Conclusion: A combination of PEFs with antibiotics reduces the electric field threshold required for bacteria disinfection. Such an approach simplifies devices required to disinfect large areas of infected wounds.

From Cell to Tissue Properties-Modeling Skin Electroporation With Pore and Local Transport Region Formation

Current models of tissue electroporation either describe tissue with its bulk properties or include cell level properties, but model only a few cells of simple shapes in low-volume fractions or are in two dimensions. We constructed a three-dimensional model of realistically shaped cells in realistic volume fractions. By using a 'unit cell' model, the equivalent dielectric properties of whole tissue could be calculated. We calculated the dielectric properties of electroporated skin. We modeled electroporation of single cells by pore formation on keratinocytes and on the papillary dermis which gave dielectric properties of the electroporated epidermis and papillary dermis. During skin electroporation, local transport regions are formed in the stratum corneum. We modeled local transport regions and increase in their radii or density which affected the dielectric properties of the stratum corneum. The final model of skin electroporation accurately describes measured electric current and voltage drop on the skin during electroporation with long low-voltage pulses. The model also accurately describes voltage drop on the skin during electroporation with short high-voltage pulses. However, our results indicate that during application of short high-voltage pulses additional processes may occur which increase the electric current. Our model connects the processes occurring at the level of cell membranes (pore formation), at the level of a skin layer (formation of local transport region in the stratum corneum) with the tissue (skin layers) and even level of organs (skin). Using a similar approach, electroporation of any tissue can be modeled, if the morphology of the tissue is known.

Development of a statistical model for cervical cancer cell death with irreversible electroporation in vitro

PURPOSE

The aim of this study was to develop a statistical model for cell death by irreversible electroporation (IRE) and to show that the statistic model is more accurate than the electric field threshold model in the literature using cervical cancer cells in vitro.

METHODS

HeLa cell line was cultured and treated with different IRE protocols in order to obtain data for modeling the statistical relationship between the cell death and pulse-setting parameters. In total, 340 in vitro experiments were performed with a commercial IRE pulse system, including a pulse generator and an electric cuvette. Trypan blue staining technique was used to evaluate cell death after 4 hours of incubation following IRE treatment. Peleg-Fermi model was used in the study to build the statistical relationship using the cell viability data obtained from the in vitro experiments. A finite element model of IRE for the electric field distribution was also built. Comparison of ablation zones between the statistical model and electric threshold model (drawn from the finite element model) was used to show the accuracy of the proposed statistical model in the description of the ablation zone and its applicability in different pulse-setting parameters.

RESULTS

The statistical models describing the relationships between HeLa cell death and pulse length and the number of pulses, respectively, were built. The values of the curve fitting parameters were obtained using the Peleg-Fermi model for the treatment of cervical cancer with IRE. The difference in the ablation zone between the statistical model and the electric threshold model was also illustrated to show the accuracy of the proposed statistical model in the representation of ablation zone in IRE.

CONCLUSIONS

This study concluded that: (1) the proposed statistical model accurately described the ablation zone of IRE with cervical cancer cells, and was more accurate compared with the electric field model; (2) the proposed statistical model was able to estimate the value of electric field threshold for the computer simulation of IRE in the treatment of cervical cancer; and (3) the proposed statistical model was able to express the change in ablation zone with the change in pulse-setting parameters.

Pulsed electric field-assisted sensitization of multidrug-resistant Candida albicans to antifungal drugs

AIM

Determine the influence of pH on the inactivation efficiency of Candida albicans in pulsed electric fields (PEF) and evaluate the possibilities for sensitization of a drug-resistant strain to antifungal drugs.

MATERIALS & METHODS

The effects of PEF (2.5-25 kVcm^-1) with fluconazole, terbinafine and naftifine were analyzed at a pH range of 3.0-9.0. Membrane permeabilization was determined by flow cytometry and propidium iodide.

RESULTS

PEF induced higher inactivation of C. albicans at low pH and increased sensitivity to terbinafine and naftifine to which the strain was initially resistant. Up to 5 log reduction in cell survival was achieved.

CONCLUSION

A proof of concept that electroporation can be used to sensitize drug-resistant microorganisms was presented, which is promising for treating biofilm-associated infections.

Anistropically varying conductivity in irreversible electroporation simulations

Background

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 propose a mathematical model that incorporates a tissue’s conductivity increasing more in the direction of the electrical field as this has been shown to occur in experiments.

Method

It was necessary to mathematically derive a valid form of the conductivity tensor such that it is dependent on the electrical field direction and can be easily implemented into numerical software. The derivation of a conductivity tensor that can take arbitrary functions for the conductivity in the directions tangent and normal to the electrical field is the main contribution of this paper. Numerical simulations were performed for isotropic-varying and anisotropic-varying conductivities to evaluate the importance of including the electrical field’s direction in the formulation for conductivity.

Results

By starting from previously published experimental results, this paper derived a general formulation for an anistropic-varying tensor for implementation into irreversible electroporation modeling software. The anistropic-varying tensor formulation allows the conductivity to take into consideration both electrical field direction and magnitude, as opposed to previous published works that only took into account electrical field magnitude.

The anisotropic formulation predicts roughly a five percent decrease in ablation size for the monopolar simulation and approximately a ten percent decrease in ablation size for the bipolar simulations. This is a positive result as previously reported results found the isotropic formulation to overpredict ablation size for both monopolar and bipolar simulations. Furthermore, it was also reported that the isotropic formulation overpredicts the ablation size more for the bipolar case than the monopolar case. Thus, our results are following the experimental trend by having a larger percentage change in volume for the bipolar case than the monopolar case.

Conclusions

The predicted volume of ablated cells decreased, and could be a possible explanation for the slight over-prediction seen by isotropic-varying formulations.

Predicting irreversible electroporation-induced tissue damage by means of magnetic resonance electrical impedance tomography

Irreversible electroporation (IRE) is gaining importance in routine clinical practice for nonthermal ablation of solid tumors. For its success, it is extremely important that the coverage and exposure time of the treated tumor to the electric field is within the specified range. Measurement of electric field distribution during the electroporation treatment can be achieved using magnetic resonance electrical impedance tomography (MREIT). Here, we show improved MREIT-enabled electroporation monitoring of IRE-treated tumors by predicting IRE-ablated tumor areas during IRE of mouse tumors in vivo. The in situ prediction is enabled by coupling MREIT with a corresponding Peleg-Fermi mathematical model to obtain more informative monitoring of IRE tissue ablation by providing cell death probability in the IRE-treated tumors. This technique can potentially be used in electroporation-based clinical applications, such as IRE tissue ablation and electrochemotherapy, to improve and assure the desired treatment outcome.

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.

APOBEC3A intratumoral DNA electroporation in mice

Human APOBEC3A (A3A) cytidine deaminase shows pro-apoptotic properties resulting from hypermutation of genomic DNA, induction of double-stranded DNA breaks (DSBs) and G1 cell cycle arrest. Given this, we evaluated the antitumor efficacy of A3A by intratumoral electroporation of an A3A expression plasmid. DNA was repeatedly electroporated into B16OVA, B16Luc tumors of C57BL/6J mice as well as the aggressive fibrosarcoma Sarc2 tumor of HLA-A*0201/DRB1*0101 transgenic mice using noninvasive plate electrodes. Intratumoral electroporation of A3A plasmid DNA resulted in regression of ~50% of small B16OVA melanoma tumors that did not rebound in the following 2 months without treatment. Larger or more aggressive tumors escaped regression when so treated. As APOBEC3A was much less efficient in provoking hypermutation and DSBs in B16OVA cells compared with human or quail cells, it is likely that APOBEC3A would be more efficient in a human setting than in a mouse model.

Electrolytic Effects During Tissue Ablation by Electroporation

Nonthermal irreversible electroporation is a new tissue ablation technique that consists of applying pulsed electric fields across cells to induce cell death by creating permanent defects in the cell membrane. Nonthermal irreversible electroporation is of interest because it allows treatment near sensitive tissue structures such as blood vessels and nerves. Two recent articles report that electrolytic reaction products at electrodes can be combined with electroporation pulses to augment and optimize tissue ablation. Those articles triggered a concern that the results of earlier studies on nonthermal irreversible electroporation may have been tainted by unaccounted for electrolytic effects. The goal of this study was to reexamine previous studies on nonthermal irreversible electroporation in the context of these articles. The study shows that the results from some of the earlier studies on nonthermal irreversible electroporation were affected by unaccounted for electrolysis, in particular the research with cells in cuvettes. It also shows that tissue ablation ascribed in the past to irreversible electroporation is actually caused by at least 3 different cytotoxic effects: irreversible electroporation without electrolysis, irreversible electroporation combined with electrolysis, and reversible electroporation combined with electrolysis. These different mechanisms may affect cell and tissue ablation in different ways, and the effects may depend on various clinical parameters such as the polarity of the electrodes, the charge delivered (voltage, number, and length of pulses), and the distance of the target tissue from the electrodes. Current clinical protocols employ ever-increasing numbers of electroporation pulses to values that are now an order of magnitude larger than those used in our first fundamental nonthermal irreversible electroporation studies in tissues. The different mechanisms of cell death, and the effect of the clinical parameters on the mechanisms may explain discrepancies between results of different clinical studies and should be taken into consideration in the design of optimal electroporation ablation protocols.

Dynamic finite-element model for efficient modelling of electric currents in electroporated tissue

In silico experiments (numerical simulations) are a valuable tool for non-invasive research of the influences of tissue properties, electrode placement and electric pulse delivery scenarios in the process of electroporation. The work described in this article was aimed at introducing time dependent effects into a finite element model developed specifically for electroporation. Reference measurements were made ex vivo on beef liver samples and experimental data were used both as an initial condition for simulation (applied pulse voltage) and as a reference value for numerical model calibration (measured pulse current). The developed numerical model is able to predict the time evolution of an electric pulse current within a 5% error over a broad range of applied pulse voltages, pulse durations and pulse repetition frequencies. Given the good agreement of the current flowing between the electrodes, we are confident that the results of our numerical model can be used both for detailed in silico research of electroporation mechanisms (giving researchers insight into time domain effects) and better treatment planning algorithms, which predict the outcome of treatment based on both spatial and temporal distributions of applied electric pulses.

The correlation between the electrode configuration and histopathology of irreversible electroporation ablations in prostate cancer patients

PURPOSE

Irreversible electroporation (IRE) is a novel minimally invasive therapy for prostate cancer using short electric pulses to ablate prostate tissue. The purpose of this study is to determine the IRE effects in prostate tissue and correlate electrode configuration with the histology of radical prostatectomy (RP) specimens. We hypothesize that the area within the electrode configuration is completely ablated and that the area within the electrode configuration is predictive for the ablated area after treatment.

METHODS

A prospective phase I/II study was conducted in 16 consecutive patients with histopathologically confirmed prostate cancer scheduled for RP. Focal or extended IRE treatment of the prostate was performed 4 weeks prior to RP. The locations of the electrodes were used to calculate the planned ablation zone. Following RP, the specimens were processed into whole-mount sections, histopathology (PA) was assessed and ablation zones were delineated. The area of the tissue alteration was determined by measuring the surface. The planned and the histological ablation zones were compared, analysed per individual patient and per protocol (focal vs. extended).

RESULTS

All cells within the electrode configuration were completely ablated and consisted only of necrotic and fibrotic tissue without leaving any viable cells. The histological ablation zone was always larger than the electrodes configuration (2.9 times larger for the 3 electrodes configuration and 2.5 times larger for the ≥4 electrode configuration). These ablation effects extended beyond the prostatic capsule in the neurovascular bundle in 13 out of 15 cases.

CONCLUSIONS

IRE in prostate cancer results in completely ablated, sharply demarcated lesions with a histological ablation zone beyond the electrode configuration. No skip lesions were observed within the electrode configuration.

CLINICAL TRIALS

ClinicalTrials.gov Identifier: NCT01790451 https://clinicaltrials.gov/ct2/show/NCT01790451.

Effects of Cisplatin Electrochemotherapy on Human Neuroblastoma Cells

Electrochemotherapy is the usage of electroporation to introduce chemotherapeutic drugs through membrane pores into target cells for cancer treatment. The effectiveness of chemotherapeutic drugs would be increased dramatically when they are used in electrochemotherapy than standard chemotherapy. In the present study, we investigated the effects of cisplatin treatment with electroporation on human SH-SY5Y neuroblastoma cells. SH-SY5Y cells were treated with different concentrations (0.15-24 µg/mL) of cisplatin and then exposed to 1500 volts per centimeter (V/cm), 100 microseconds (µs) pulse duration, and 1 Hertz (Hz) electric pulses. Cisplatin alone showed a dose-dependent effect on cell viability. On the other hand, cisplatin + electroporation treatment was more effective than cisplatin treatment alone. Lower doses of cisplatin treatment with electroporation was as effective as higher doses of cisplatin treatment without electroporation. These results indicated that cisplatin cytotoxicity was potentiated after exposure of cells to high intensity electric pulses and low doses of cisplatin can be used with electroporation in the treatment of neuroblastoma.

A statistical model describing combined irreversible electroporation and electroporation-induced blood-brain barrier disruption

Background

Electroporation-based therapies such as electrochemotherapy (ECT) and irreversible electroporation (IRE) are emerging as promising tools for treatment of tumors. When applied to the brain, electroporation can also induce transient blood-brain-barrier (BBB) disruption in volumes extending beyond IRE, thus enabling efficient drug penetration. The main objective of this study was to develop a statistical model predicting cell death and BBB disruption induced by electroporation. This model can be used for individual treatment planning.

Material and methods

Cell death and BBB disruption models were developed based on the Peleg-Fermi model in combination with numerical models of the electric field. The model calculates the electric field thresholds for cell kill and BBB disruption and describes the dependence on the number of treatment pulses. The model was validated using in vivo experimental data consisting of rats brains MRIs post electroporation treatments.

Results

Linear regression analysis confirmed that the model described the IRE and BBB disruption volumes as a function of treatment pulses number (r² = 0.79; p < 0.008, r² = 0.91; p < 0.001). The results presented a strong plateau effect as the pulse number increased. The ratio between complete cell death and no cell death thresholds was relatively narrow (between 0.88-0.91) even for small numbers of pulses and depended weakly on the number of pulses. For BBB disruption, the ratio increased with the number of pulses. BBB disruption radii were on average 67% ± 11% larger than IRE volumes.

Conclusions

The statistical model can be used to describe the dependence of treatment-effects on the number of pulses independent of the experimental setup.

Towards electroporation based treatment planning considering electric field induced muscle contractions

The electric field threshold for muscle contraction is two orders of magnitudes lower than that for electroporation. Current electroporation treatment planning and electrode design studies focus on optimizing the delivery of electroporation electric fields to the targeted tissue. The goal of one part of this study was to investigate the relation between the volumes of tissue that experience electroporation electric fields in a targeted tissue volume and the volumes of tissue that experience muscle contraction inducing electric fields around the electroporated tissue volume, (V(MC)), during standard electroporation procedures and for various electroporation electrodes designs. The numerical analysis shows that conventional electroporation protocols and electrode design can generate muscle contraction inducing electric fields in surprisingly large volumes of non-target tissue, around the electroporation treated tissue. In studying various electrode configurations, we found that electrode placement in a structure we refer to as a "Current Cage" can substantially reduce the volume of non-target tissue exposed to electric fields above the muscle contraction threshold. In an experimental study on a tissue phantom we compare a commercial two parallel needle electroporation system with the Current Cage design. While tissue electroporated volumes were similar, V(MC) of tissue treated using the Current Cage design electrodes was an order of magnitude smaller than that using a commercially available system. An important aspect of the entire study is that it suggests the benefit of including the calculations of V(MC) for planning of electroporation based treatments such as DNA vaccination, electrochemotherapy and irreversible electroporation.

Mathematical Models Describing Chinese Hamster Ovary Cell Death Due to Electroporation In Vitro

Electroporation is a phenomenon used in the treatment of tumors by electrochemotherapy, non-thermal ablation with irreversible electroporation, and gene therapy. When treating patients, either predefined or variable electrode geometry is used. Optimal pulse parameters are predetermined for predefined electrode geometry, while they must be calculated for each specific case for variable electrode geometry. The position and number of electrodes are also determined for each patient. It is currently assumed that above a certain experimentally determined value of electric field, all cells are permeabilized/destroyed and under it they are unaffected. In this paper, mathematical models of survival in which the probability of cell death is continuously distributed from 0 to 100 % are proposed and evaluated. Experiments were performed on cell suspensions using electrical parameters similar to standard electrochemotherapy and irreversible electroporation parameters. The proportion of surviving cells was determined using clonogenic assay for assessing the ability of a cell to grow into a colony. Various mathematical models (first-order kinetics, Hülsheger, Peleg-Fermi, Weibull, logistic, adapted Gompertz, Geeraerd) were fitted to experimental data using a non-linear least-squares method. The fit was evaluated by calculating goodness of fit and by observing the trend of values of models' parameters. The most appropriate models of cell survival as a function of treatment time were the adapted Gompertz and the Geeraerd models and, as a function of the electric field, the logistic, adapted Gompertz and Peleg-Fermi models. The next steps to be performed are validation of the most appropriate models on tissues and determination of the models' predictive power.

Effect of blood vessel segmentation on the outcome of electroporation-based treatments of liver tumors

Electroporation-based treatments rely on increasing the permeability of the cell membrane by high voltage electric pulses applied to tissue via electrodes. To ensure that the whole tumor is covered with sufficiently high electric field, accurate numerical models are built based on individual patient anatomy. Extraction of patient's anatomy through segmentation of medical images inevitably produces some errors. In order to ensure the robustness of treatment planning, it is necessary to evaluate the potential effect of such errors on the electric field distribution. In this work we focus on determining the effect of errors in automatic segmentation of hepatic vessels on the electric field distribution in electroporation-based treatments in the liver. First, a numerical analysis was performed on a simple 'sphere and cylinder' model for tumors and vessels of different sizes and relative positions. Second, an analysis of two models extracted from medical images of real patients in which we introduced variations of an error of the automatic vessel segmentation method was performed. The results obtained from a simple model indicate that ignoring the vessels when calculating the electric field distribution can cause insufficient coverage of the tumor with electric fields. Results of this study indicate that this effect happens for small (10 mm) and medium-sized (30 mm) tumors, especially in the absence of a central electrode inserted in the tumor. The results obtained from the real-case models also show higher negative impact of automatic vessel segmentation errors on the electric field distribution when the central electrode is absent. However, the average error of the automatic vessel segmentation did not have an impact on the electric field distribution if the central electrode was present. This suggests the algorithm is robust enough to be used in creating a model for treatment parameter optimization, but with a central electrode.

Tissue heterogeneity in structure and conductivity contribute to cell survival during irreversible electroporation ablation by "electric field sinks"

Irreversible electroporation (IRE) is an emerging, minimally invasive technique for solid tumors ablation, under clinical investigation for cancer therapy. IRE affects only the cell membrane, killing cells while preserving the extracellular matrix structure. Current reports indicate tumors recurrence rate after IRE averaging 31% of the cases, of which 10% are local recurrences. The mechanisms for these recurrences are not known and new explanations for incomplete cell death are needed. Using finite elements method for electric field distribution, we show that presence of vascular structures with blood leads to the redistribution of electric fields leading to the areas with more than 60% reduced electric field strength in proximity to large blood vessels and clustered vessel structures. In an in vivo rat model of liver IRE ablation, we show that cells located in the proximity of larger vessel structures and in proximity of clustered vessel structures appear less affected by IRE ablation than cells in the tissue parenchyma or in the proximity of small, more isolated vessels. These findings suggest a role for “electric field sinks” in local tumors recurrences after IRE and emphasize the importance of the precise mapping of the targeted organ structure and conductivity for planning of electroporation procedures.

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.

A numerical investigation of the electric and thermal cell kill distributions in electroporation-based therapies in tissue

Electroporation-based therapies are powerful biotechnological tools for enhancing the delivery of exogeneous agents or killing tissue with pulsed electric fields (PEFs). Electrochemotherapy (ECT) and gene therapy based on gene electrotransfer (EGT) both use reversible electroporation to deliver chemotherapeutics or plasmid DNA into cells, respectively. In both ECT and EGT, the goal is to permeabilize the cell membrane while maintaining high cell viability in order to facilitate drug or gene transport into the cell cytoplasm and induce a therapeutic response. Irreversible electroporation (IRE) results in cell kill due to exposure to PEFs without drugs and is under clinical evaluation for treating otherwise unresectable tumors. These PEF therapies rely mainly on the electric field distributions and do not require changes in tissue temperature for their effectiveness. However, in immediate vicinity of the electrodes the treatment may results in cell kill due to thermal damage because of the inhomogeneous electric field distribution and high current density during the electroporation-based therapies. Therefore, the main objective of this numerical study is to evaluate the influence of pulse number and electrical conductivity in the predicted cell kill zone due to irreversible electroporation and thermal damage. Specifically, we simulated a typical IRE protocol that employs ninety 100-µs PEFs. Our results confirm that it is possible to achieve predominant cell kill due to electroporation if the PEF parameters are chosen carefully. However, if either the pulse number and/or the tissue conductivity are too high, there is also potential to achieve cell kill due to thermal damage in the immediate vicinity of the electrodes. Therefore, it is critical for physicians to be mindful of placement of electrodes with respect to critical tissue structures and treatment parameters in order to maintain the non-thermal benefits of electroporation and prevent unnecessary damage to surrounding healthy tissue, critical vascular structures, and/or adjacent organs.

Planning irreversible electroporation in the porcine kidney: are numerical simulations reliable for predicting empiric ablation outcomes?

PURPOSE

Numerical simulations are used for treatment planning in clinical applications of irreversible electroporation (IRE) to determine ablation size and shape. To assess the reliability of simulations for treatment planning, we compared simulation results with empiric outcomes of renal IRE using computed tomography (CT) and histology in an animal model.

METHODS

The ablation size and shape for six different IRE parameter sets (70-90 pulses, 2,000-2,700 V, 70-100 µs) for monopolar and bipolar electrodes was simulated using a numerical model. Employing these treatment parameters, 35 CT-guided IRE ablations were created in both kidneys of six pigs and followed up with CT immediately and after 24 h. Histopathology was analyzed from postablation day 1.

RESULTS

Ablation zones on CT measured 81 ± 18 % (day 0, p ≤ 0.05) and 115 ± 18 % (day 1, p ≤ 0.09) of the simulated size for monopolar electrodes, and 190 ± 33 % (day 0, p ≤ 0.001) and 234 ± 12 % (day 1, p ≤ 0.0001) for bipolar electrodes. Histopathology indicated smaller ablation zones than simulated (71 ± 41 %, p ≤ 0.047) and measured on CT (47 ± 16 %, p ≤ 0.005) with complete ablation of kidney parenchyma within the central zone and incomplete ablation in the periphery.

CONCLUSION

Both numerical simulations for planning renal IRE and CT measurements may overestimate the size of ablation compared to histology, and ablation effects may be incomplete in the periphery.

Rapid dramatic alterations to the tumor microstructure in pancreatic cancer following irreversible electroporation ablation

AIM

NanoKnife(®) (Angiodynamics, Inc., NY, USA) or irreversible electroporation (IRE) is a newly available ablation technique to induce the formation of nanoscale pores within the cell membrane in targeted tissues. The purpose of this study was to elucidate morphological alterations following 30 min of IRE ablation in a mouse model of pancreatic cancer.

MATERIALS & METHODS

Immunohistochemistry markers were compared with diffusion-weighted MRI apparent diffusion coefficient measurements before and after IRE ablation.

RESULTS

Immunohistochemistry apoptosis index measurements were significantly higher in IRE-treated tumors than in controls. Rapid tissue alterations after 30 min of IRE ablation procedures (structural and morphological alterations along with significantly elevated apoptosis markers) were consistently observed and well correlated to apparent diffusion coefficient measurements.

DISCUSSION

This imaging assay offers the potential to serve as an in vivo biomarker for noninvasive detection of tumor response following IRE ablation.