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

Calculating transmembrane voltage on the electric pulse-affected cancerous cell membrane: using molecular dynamics and finite element simulations

CONTEXT

Electroporation is a technique that creates electrically generated pores in the cell membrane by modifying transmembrane potential. In this work, the finite element method (FEM) was used to examine the induced transmembrane voltage (ITV) of a spherical-shaped MCF-7 cell, allowing researchers to determine the stationary ITV. A greater ITV than the critical value causes permeabilization of the membrane. Furthermore, the present study shows how a specific surface conductivity can act as a stand-in for the thin layer that constitutes a cell membrane as the barrier between extracellular and intracellular environments. Additionally, the distribution of ITV on the cell membrane and its maximum value were experimentally evaluated for a range of applied electric fields. Consequently, the entire cell surface area was electroporated 66% and 68% for molecular dynamics (MD) simulations and FEM, respectively, when the external electric field of 1500 V/cm was applied to the cell suspension using the previously indicated numerical methods. Furthermore, the lipid bilayers' molecular structure was changed, which led to the development of hydrophilic holes with a radius of 1.33 nm. Applying MD and FEM yielded threshold values for transmembrane voltage of 700 and 739 mV, respectively.

METHOD

Using MD simulations of palmitoyloleoyl-phosphatidylcholine (POPC), pores in cell membranes exposed to external electric fields were numerically investigated. The dependence on the electric field was estimated and developed, and the amount of the electroporated cell surface area matches the applied external electric field. To investigate more, a mathematical model based on an adaptive neuro-fuzzy inference system (ANFIS) is employed to predict the percent cell viability of cancerous cells after applying four pulses during electroporation. For MD simulations, ArgusLab, VMD, and GROMACS software packages were used. Moreover, for FEM analysis, COMSOL software package was used. Also, it is worth mentioning that for mathematical model, MATLAB software is used.

Determination of lethal electric field threshold for pulsed field ablation in ex vivo perfused porcine and human hearts

Introduction

Pulsed field ablation is an emerging modality for catheter-based cardiac ablation. The main mechanism of action is irreversible electroporation (IRE), a threshold-based phenomenon in which cells die after exposure to intense pulsed electric fields. Lethal electric field threshold for IRE is a tissue property that determines treatment feasibility and enables the development of new devices and therapeutic applications, but it is greatly dependent on the number of pulses and their duration.

Methods

In the study, lesions were generated by applying IRE in porcine and human left ventricles using a pair of parallel needle electrodes at different voltages (500-1500 V) and two different pulse waveforms: a proprietary biphasic waveform (Medtronic) and monophasic 48 × 100 μs pulses. The lethal electric field threshold, anisotropy ratio, and conductivity increase by electroporation were determined by numerical modeling, comparing the model outputs with segmented lesion images.

Results

The median threshold was 535 V/cm in porcine ((N = 51 lesions in n = 6 hearts) and 416 V/cm in the human donor hearts ((N = 21 lesions in n = 3 hearts) for the biphasic waveform. The median threshold value was 368 V/cm in porcine hearts ((N = 35 lesions in n = 9 hearts) cm for 48 × 100 μs pulses.

Discussion

The values obtained are compared with an extensive literature review of published lethal electric field thresholds in other tissues and were found to be lower than most other tissues, except for skeletal muscle. These findings, albeit preliminary, from a limited number of hearts suggest that treatments in humans with parameters optimized in pigs should result in equal or greater lesions.

Electroporation-Based Therapy for Brain Tumors: A Review

Electroporation-based therapy (EBT), as a high-voltage-pulse technology has been prevalent with favorable clinical outcomes in the treatment of various solid tumors. This review paper aims to promote the clinical translation of EBT for brain tumors. First, we briefly introduced the mechanism of pore formation in a cell membrane activated by external electric fields using a single cell model. Then, we summarized and discussed the current in vitro and in vivo preclinical studies, in terms of (1) the safety and effectiveness of EBT for brain tumors in animal models, and (2) the blood-brain barrier (BBB) disruption induced by EBT. Two therapeutic effects could be achieved in EBT for brain tumors simultaneously, i.e., the tumor ablation induced by irreversible electroporation (IRE) and transient BBB disruption induced by reversible electroporation (RE). The BBB disruption could potentially improve the uptake of antitumor drugs thereby enhancing brain tumor treatment. The challenges that hinder the application of EBT in the treatment of human brain tumors are discussed in the review paper as well.

Characterization of Nonlinearity and Dispersion in Tissue Impedance During High-Frequency Electroporation

OBJECTIVE

The use of high-voltage, high-frequency bipolar pulses (HFBPs) is an emerging electroporation-based therapy for the treatment of solid tumors. In this study, we quantify the extent of nonlinearity and dispersion during the HFBP treatment.

METHODS

We utilize flat-plate electrodes to capture the impedance of the porcine liver tissue during the delivery of a burst of HFBPs of widths 1 and 2  $\mu$s at different pulse amplitudes. Next, we fit the impedance data to a frequency-dependent parallel RC network to determine the conductivity and permittivity of the tissue as a function of frequency, for different applied electric fields. Finally, we present a simple model to approximate the field distribution in the tissue using the conductivity function at a frequency that could minimize the errors due to approximation with a nondispersive model.

RESULTS

The conductivity/permittivity of the tissue was plotted as a function of frequency for different electric fields. It was found that the extent of dispersion reduces with higher applied electric field magnitudes.

CONCLUSION

This is the first study to quantify dispersion and nonlinearity in the tissue during the HFBP treatment. The data have been used to predict the field distribution in a numerical model of the liver tissue utilizing two needle electrodes.

SIGNIFICANCE

The data and technique developed in this study to monitor the electrical properties of tissue during treatment can be used to generate treatment-planning models for future high-frequency electroporation therapies as well as provide insights regarding treatment effect.

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.

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.

Irreversible Electroporation (IRE) in Renal Tumors

Small renal masses (SRMs) have been traditionally managed with surgical resection. Minimally invasive nephron-sparing treatment methods are preferred to avoid harmful consequences of renal insufficiency, with partial nephrectomy (PN) considered the gold standard. With increase in the incidence of the SRMs and evolution of ablative technologies, percutaneous ablation is now considered a viable treatment alternative to surgical resection with comparable oncologic outcomes and better nephron-sparing property. Traditional thermal ablative techniques suffer from unique set of challenges in treating tumors near vessels or critical structures. Irreversible electroporation (IRE), with its non-thermal nature and connective tissue-sparing properties, has shown utility where traditional ablative techniques face challenges. This review presents the role of IRE in renal tumors based on the most relevant published literature on the IRE technology, animal studies, and human experience.

Irreversible electroporation: state of the art

The field of focal ablative therapy for the treatment of cancer is characterized by abundance of thermal ablative techniques that provide a minimally invasive treatment option in selected tumors. However, the unselective destruction inflicted by thermal ablation modalities can result in damage to vital structures in the vicinity of the tumor. Furthermore, the efficacy of thermal ablation intensity can be impaired due to thermal sink caused by large blood vessels in the proximity of the tumor. Irreversible electroporation (IRE) is a novel ablation modality based on the principle of electroporation or electropermeabilization, in which electric pulses are used to create nanoscale defects in the cell membrane. In theory, IRE has the potential of overcoming the aforementioned limitations of thermal ablation techniques. This review provides a description of the principle of IRE, combined with an overview of in vivo research performed to date in the liver, pancreas, kidney, and prostate.

Histological and Mathematical Analysis of the Irreversibly Electroporated Liver Tissue

Irreversible electroporation has clinically been used to treat various types of cancer. A plan on how to apply irreversible electroporation before practicing is very important to increase the ablation area and reduce the side effects. Several electrical models have been developed to predict the ablation area with applied electric energy. In this experiment, the static relationship between applied electric energy and ablated area was mathematically and experimentally investigated at 10 hours after applying irreversible electroporation. We performed the irreversible electroporation on the liver tissue of Sprague Dawley rats (male, 8 weeks, weighing 250-350 g). The ablated area was measured based on histological analysis and compared with the mathematical calculation from the electric energy, assuming that the tissue is homogeneous. The ablated area increased with the increase in applied electric energy. The numerically calculated contour lines of electric energy density overlapped well with the apoptotic area induced by the irreversible electroporation. The overlapped area clearly showed that the destructive threshold of apoptosis between electrodes is electric energy density level of 5.9 × 10⁵ J/m³. The results of the present study suggested that the clinical results of the irreversible electroporation on a liver tissue could be predicted through mathematical calculation.

Comparison of ablation defect on MR imaging with computer simulation estimated treatment zone following irreversible electroporation of patient prostate

To determine whether patient specific numerical simulations of irreversible electroporation (IRE) of the prostate correlates with the treatment effect seen on follow-up MR imaging. Computer models were created using intra-operative US images, post-treatment follow-up MR images and clinical data from six patients receiving IRE. Isoelectric contours drawn using simulation results were compared with MR imaging to estimate the energy threshold separating treated and untreated tissue. Simulation estimates of injury to the neurovascular bundle and rectum were compared with clinical follow-up and patient reported outcomes. At the electric field strength of 700 V/cm, simulation estimated electric field distribution was not different from the ablation defect seen on follow-up MR imaging (p = 0.43). Simulation predicted cross sectional area of treatment (mean 532.33 ± 142.32 mm(2)) corresponded well with the treatment zone seen on MR imaging (mean 540.16 ± 237.13 mm(2)). Simulation results did not suggest injury to the rectum or neurovascular bundle, matching clinical follow-up at 3 months. Computer simulation estimated zone of irreversible electroporation in the prostate at 700 V/cm was comparable to measurements made on follow-up MR imaging. Numerical simulation may aid treatment planning for irreversible electroporation of the prostate in patients.

Web-based tool for visualization of electric field distribution in deep-seated body structures and planning of electroporation-based treatments

BACKGROUND

Treatments based on electroporation are a new and promising approach to treating tumors, especially non-resectable ones. The success of the treatment is, however, heavily dependent on coverage of the entire tumor volume with a sufficiently high electric field. Ensuring complete coverage in the case of deep-seated tumors is not trivial and can in best way be ensured by patient-specific treatment planning. The basis of the treatment planning process consists of two complex tasks: medical image segmentation, and numerical modeling and optimization.

METHODS

In addition to previously developed segmentation algorithms for several tissues (human liver, hepatic vessels, bone tissue and canine brain) and the algorithms for numerical modeling and optimization of treatment parameters, we developed a web-based tool to facilitate the translation of the algorithms and their application in the clinic. The developed web-based tool automatically builds a 3D model of the target tissue from the medical images uploaded by the user and then uses this 3D model to optimize treatment parameters. The tool enables the user to validate the results of the automatic segmentation and make corrections if necessary before delivering the final treatment plan.

RESULTS

Evaluation of the tool was performed by five independent experts from four different institutions. During the evaluation, we gathered data concerning user experience and measured performance times for different components of the tool. Both user reports and performance times show significant reduction in treatment-planning complexity and time-consumption from 1-2 days to a few hours.

CONCLUSIONS

The presented web-based tool is intended to facilitate the treatment planning process and reduce the time needed for it. It is crucial for facilitating expansion of electroporation-based treatments in the clinic and ensuring reliable treatment for the patients. The additional value of the tool is the possibility of easy upgrade and integration of modules with new functionalities as they are developed.

A prospective Phase 2a pilot study investigating focal percutaneous irreversible electroporation (IRE) ablation by NanoKnife in patients with localised renal cell carcinoma (RCC) with delayed interval tumour resection (IRENE trial)

INTRODUCTION

Focal ablation therapy is playing an increasing role in oncology and may reduce the toxicity of current surgical treatments while achieving adequate oncological benefit. Irreversible electroporation (IRE) has been proposed to be tissue-selective with potential advantages compared with current thermal-ablation technologies or radiotherapy. The aim of this pilot trial is to determine the effectiveness and feasibility of focal percutaneous IRE in patients with localised renal cell cancer as a uro-oncological tumour model.

METHODS

Prospective, monocentric Phase 2a pilot study following current recommendations, including those of the International Working Group on Image-Guided Tumor Ablation. Twenty patients with kidney tumour (T1aN0M0) will be recruited. This sample permits an appropriate evaluation of the feasibility and effectiveness of image-guided percutaneous IRE ablation of locally confined kidney tumours as well as functional outcomes. Percutaneous biopsy for histopathology will be performed before IRE, with magnetic-resonance imaging one day before and 2, 7, 27 and 112 days after IRE; at 28 days after IRE the tumour region will be completely resected and analysed by ultra-thin-layer histology.

DISCUSSION

The IRENE study will investigate over a short-term observation period (by magnetic-resonance imaging, post-resection histology and assessment of technical feasibility) whether focal IRE, as a new ablation procedure for soft tissue, is feasible as a percutaneous, tissue-sparing method for complete ablation and cure of localised kidney tumours. Results from the kidney-tumour model can provide guidance for designing an effectiveness and feasibility trial to assess this new ablative technology, particularly in uro-oncology.

Feasibility of catheter-directed intraluminal irreversible electroporation of porcine ureter and acute outcomes in response to increasing energy delivery

PURPOSE

To evaluate the feasibility of focal intraluminal irreversible electroporation (IRE) in the ureter with a novel electrode catheter and to study the treatment effects in response to increasing pulse strength.

MATERIALS AND METHODS

Five IRE treatment settings were each evaluated twice for the ablation of normal ureter in 5 Yorkshire pigs (n = 1-4 ablations per animal; total of 10 ablations) with the use of a prototype device under ultrasound and fluoroscopic guidance. Animals received unilateral or bilateral treatment, limited to a maximum of 2 ablations in any 1 ureter. Treatment was delivered with increasing pulse strength (from 1,000 V to 3,000 V in increments of 500 V) while keeping the pulse duration (100 μs) and number of pulses (n = 90) constant. Ureter patency was assessed with antegrade ureteropyelography immediately following treatment. Animals were euthanized within 4 hours after treatment, and treated urinary tract was harvested for histopathologic analysis with hematoxylin and eosin and Masson trichrome stains.

RESULTS

IRE was successfully performed in all animals, without evidence of ureteral perforation. Hematoxylin and eosin analysis of IRE treatments demonstrated full-thickness ablation at higher field strengths (mucosa to the adventitia). Masson trichrome stains showed preservation of connective tissue at all field strengths.

CONCLUSIONS

Intraluminal catheter-directed IRE ablation is feasible and produces full-thickness ablation of normal ureters. There was no evidence of lumen perforation even at the maximum voltages evaluated.

In vivo irreversible electroporation kidney ablation: experimentally correlated numerical models

Irreversible electroporation (IRE) ablation uses brief electric pulses to kill a volume of tissue without damaging the structures contraindicated for surgical resection or thermal ablation, including blood vessels and ureters. IRE offers a targeted nephron-sparing approach for treating kidney tumors, but the relevant organ-specific electrical properties and cellular susceptibility to IRE electric pulses remain to be characterized. Here, a pulse protocol of 100 electric pulses, each 100 μs long, is delivered at 1 pulse/s to canine kidneys at three different voltage-to-distance ratios while measuring intrapulse current, completed 6 h before humane euthanasia. Numerical models were correlated with lesions and electrical measurements to determine electrical conductivity behavior and lethal electric field threshold. Three methods for modeling tissue response to the pulses were investigated (static, linear dynamic, and asymmetrical sigmoid dynamic), where the asymmetrical sigmoid dynamic conductivity function most accurately and precisely matched lesion dimensions, with a lethal electric field threshold of 575 ± 67 V/cm for the protocols used. The linear dynamic model also attains accurate predictions with a simpler function. These findings can aid renal IRE treatment planning under varying electrode geometries and pulse strengths. Histology showed a wholly necrotic core lesion at the highest electric fields, surrounded by a transitional perimeter of differential tissue viability dependent on renal structure.