A Parametric Study Delineating Irreversible Electroporation from Thermal Damage Based on a Minimally Invasive Intracranial Procedure
© Garcia et al; licensee BioMed Central Ltd. 2011
Received: 18 January 2011
Accepted: 30 April 2011
Published: 30 April 2011
Irreversible electroporation (IRE) is a new minimally invasive technique to kill undesirable tissue in a non-thermal manner. In order to maximize the benefits from an IRE procedure, the pulse parameters and electrode configuration must be optimized to achieve complete coverage of the targeted tissue while preventing thermal damage due to excessive Joule heating.
We developed numerical simulations of typical protocols based on a previously published computed tomographic (CT) guided in vivo procedure. These models were adapted to assess the effects of temperature, electroporation, pulse duration, and repetition rate on the volumes of tissue undergoing IRE alone or in superposition with thermal damage.
Nine different combinations of voltage and pulse frequency were investigated, five of which resulted in IRE alone while four produced IRE in superposition with thermal damage.
The parametric study evaluated the influence of pulse frequency and applied voltage on treatment volumes, and refined a proposed method to delineate IRE from thermal damage. We confirm that determining an IRE treatment protocol requires incorporating all the physical effects of electroporation, and that these effects may have significant implications in treatment planning and outcome assessment. The goal of the manuscript is to provide the reader with the numerical methods to assess multiple-pulse electroporation treatment protocols in order to isolate IRE from thermal damage and capitalize on the benefits of a non-thermal mode of tissue ablation.
Irreversible electroporation is a new technique for the focal ablation of undesirable tissue using high voltage, low energy electric pulses [1, 2]. An IRE treatment involves placing electrodes within the region of interest and delivering a series of electric pulses that are microseconds in duration . The pulses create an electric field that induces an increase in the resting transmembrane potential (TMP) of the cells in the tissue . The induced increase in the TMP is dependent on the electric pulse (e.g. strength, duration, repetition rate, shape, and number) and physical configuration of the electrodes used to deliver the pulses. Depending on the magnitude of the induced TMP, as well as its duration and repetition rate for induction, the electric pulses can have no effect, transiently increase membrane permeability, or cause cell death . Spatially, for a given set of conditions, the TMP and therefore the degree of electroporation is dependent on the local electric field to which the cells are exposed. Because the transitions in cellular response to the electric pulses are sudden, the treated regions are sharply delineated. Consequently, numerical models that simulate the electric field distributions in tissue are needed to predict the treated region [6–8].
There have been several studies evaluating the efficacy and safety of IRE in treating both experimental and spontaneous tumors. Al-Sakere et al. subcutaneously implanted sarcoma tumors in mice and achieved a complete response in 12 of 13 tumors with IRE treatment . Guo et al. achieved regression of hepatocellular carcinoma tumors implanted in liver in 9 out of 10 rats treated with IRE . Neal et al. implanted human mammary tumors orthotopically in mice and produced a complete response in 5 of 7 tumors with IRE which verified that IRE can be used in a heterogeneous environment . In a clinical series of IRE based therapies, our group has long-term follow-up on canine patients with spontaneous tumors. One canine patient was treated with IRE and radiation therapy for a non-resectable, high-grade glioma, resulting in complete remission of the tumor at four months . Another canine patient with a focal histiocytic sarcoma has been in complete remission for 8 months since completion of the last IRE treatment .
One of the main advantages of IRE over other focal ablation techniques is that the therapy does not use thermal damage from Joule heating to kill the cells. As a result, major blood vessels, extracellular matrix and other critical structures are spared [1, 2]. Because electroporation based therapies require high-voltage pulses to be administered to the tissue, thermistors and thermocouples may become damaged during treatment. Therefore, previous investigations into the thermal aspects of electroporation based therapies have relied on numerical modeling, typically using a modified Pennes' Bioheat equation with an added Joule heating term to predict the thermal effects. There have been several theoretical attempts in the literature to investigate the thermal response of tissues to electroporation-based treatments and assess the degree, if any, of thermal damage. In some studies, the authors calculate the pulse time required to reach a maximum temperature of 50°C, which they assume is when instantaneous thermal damage will occur [13, 14]. Others calculate the equivalent thermal dose or thermal damage associated with one or multiple pulses to determine the amount, if any, of tissue damage due to exposure of the tissue to elevated temperatures [4, 15–19]. Finally, other papers show the equivalent thermal dose for an 80-pulse IRE treatment . Pliquett et al. performed a qualitative assessment of thermal effects induced by electroporation by using temperature-sensitive liquid crystals that change colors at 40°C, 45°C, and 50°C . Although these theoretical and qualitative analyses are very powerful and well-grounded, to the best of our knowledge there is no experimental data for actual temperature changes during IRE pulse administration to prove that cell death occurs independent of classical thermal-induced mechanisms. This data is vital in order to validate the numerical models and better predict the temperature changes during a procedure for thermal damage assessment. In addition, the numerical models assessing thermal damage in the literature do not simultaneously incorporate the significant changes in the electrical conductivity of the tissue due to temperature changes as well as electroporation. Therefore, models of electroporation-based protocols that include the electrical conductivity changes and do not assume that the heat will dissipate by the beginning of the following pulse are needed in order capture the entire thermal effects of a procedure. It is then possible to maximize the sparing of critical structures in the brain and other organs and to determine the upper limit of the IRE treatment, above which thermal damage ensues. Accurate prediction of all treatment associated effects is vital to the development and implementation of optimized treatment protocols.
Our group has confirmed the safety of intracranial IRE procedures in three experimental canines . These procedures were performed through craniectomy defects to expose the cortex (grey matter) and allow for the insertion of the electrodes in the brain. We have also correlated numerical models with 3D lesion reconstructions in order to establish electric field intensities needed to kill grey matter . These studies have shown that IRE has the potential to treat intracranial disorders in canine and human patients. In the present study, we use a previously reported treatment performed through 1.2 mm diameter burr holes with CT guidance placement within a subcortical neuroanatomic target as the basis for a parametric study .
The parametric study in brain tissue evaluates the effects that the change in tissue electrical conductivity due to electroporation and the thermal effects have on the electric field distribution. It further simulates treatment volumes for similar procedures performed at three frequencies that have been used clinically in other tissues including prostate, kidney and lungs [24–27]. This study demonstrates how one can use an Arrhenius analysis to relate temperature and length of exposure during electroporation-based procedures.
A neuromuscular blocker was administered to suppress patient motion prior to the IRE treatment . A focal ablative IRE lesion was created in the white matter of the brain using the NanoKnife® generator (Angiodynamics, Queensbury, NY USA) . Two 1-mm diameter blunt tip electrodes with 5 mm exposure were inserted into the brain through the burr holes with a center-to-center separation distance of 5 mm . After insertion of the electrodes, four sets of twenty, 50 μs long, electric pulses were delivered with an applied voltage of 500 V . The polarity of the electrodes was alternated between the sets to minimize charge build-up on the electrode surface . These parameters were determined from our previous in vivo intracranial IRE procedures which showed that they were sufficient to ablate grey matter [21, 22, 28]. The NanoKnife® pulses were synchronized with the canine's heart rate in order to prevent ventricular fibrillation or other cardiac arrhythmias and were delivered in trains of ten . Due to recharging demands of the capacitors, each train of ten pulses was delivered 3.5 seconds after completion of the previous train.
Immediately after completion of the pulses, the subject was imaged with CT and a 0.2 T magnetic resonance imaging (MRI) system. The animal was then humanely euthanized 2 hours post-IRE by intravenous barbiturates and the brain was harvested and fixed in 10% buffered formalin. The fixed ex vivo brain was later imaged on a 7.0 T MRI for a more detailed analysis of the lesion produced.
Numerical models can be used for treatment planning to ensure that only the targeted regions are ablated . In order for this to be accurate, one has to know both the physical properties of the tissue, the electric field distribution, and the electric field threshold needed for IRE. This study examined two sets of models. The first was developed to replicate the experimental procedure and used the temperature and current data to calibrate the properties and behavior of the tissue in response to the electric pulses. After calibrating the model with properties based on the experimental procedure, the model was adjusted to simulate treatments at three pulse repetition rates (0.5, 1, and 4 Hz) and three voltages (500, 1000, and 1500 V) for up to 80 pulses. The computations were performed with a commercial finite element package (Comsol Multiphysics, v.3.5a, Stockholm, Sweden).
Electric Field Distribution
Initially, the 3D simulation was solved for a negligible fraction of the total treatment duration under homogenous tissue conditions in order to establish a baseline electric field distribution. The homogeneous electric field map provides the starting values for the dynamic conductivity function . In our function, we assumed that the conductivity would increase by a factor of 3.0 due to electroporation, since this is similar to the reported factor in other organs during electroporation [30, 33–37]. Additionally, this factor matches the experimental current (data not shown) that was measured by the NanoKnife® after the transient membrane charging effects had settled during the delivery of the pulses .
Physical properties used in the numerical simulations
α, temperature coefficient
k, thermal conductivity
W m-1 K-1
Cp , heat capacity
J kg-1 K-1
q''', metabolic heat generation
Cb , heat capacity
J kg-1 K-1
ρb , density
wb , perfusion rate
σ, electrical conductivity
k, thermal conductivity
W m-1 K-1
Cp , heat capacity
J kg-1 K-1
σ, electrical conductivity
k, thermal conductivity
W m -1 K -1
Cp , heat capacity
J kg-1 K-1
Instead of modeling eighty individual pulses, we modified the approach to have a continuous delivery of the electric field since we assume that once the conductivity increased due to electroporation it would not revert back. Using this approach eliminates the need to manipulate the time steps in order to ensure that the microsecond pulses are captured by the solver. This helps the simulation run faster and smoother since there are no abrupt changes due to the pulses. In order to deliver the same amount of energy as in the pulsed approach, we multiplied the Joule heating by the duty cycle (duration/period) of the pulse in the tissue and insulation domains. This ensures that at the onset of each pulse, equal amounts of thermal energy have been deposited in the tissue using either approach.
where k is the thermal conductivity of the tissue, T is the temperature above the arterial temperature (T a = 37°C), w b is the blood perfusion rate, C b is the heat capacity of the blood, ρ b the blood density, q''' is the metabolic heat generation, ρ is the tissue density, and C p is the heat capacity of the tissue. Several thermal boundary conditions can be employed to study the heat exchange between the electrodes and the tissue [13, 17, 38]. In our models, the electrodes were considered as heat sinks, , which dissipate heat from the tissue through the electrodes to the environment [19, 22].
Thermal Damage Distribution
where ζ is the frequency factor, E a the activation energy, R the universal gas constant, T(t) is the temperature distribution and τ is the heating time [4, 15, 44, 45]. It has been shown that Ω = 0.53 is the threshold for burn injuries in blood-perfused skin [45–47]. We have adapted the Arrhenius equation, which traditionally has been used to study burn injuries in skin and transdermal drug delivery using electroporation, to investigate therapeutic IRE.
Experimental Temperature Distribution
Figure 5 also includes the calculated temperature (dashed) distribution from the calibration model at the two locations where the thermal probes were positioned experimentally. This numerical simulation replicated all aspects of the experimental procedure including the four sets of twenty pulses and the 3.5 seconds delay after the first ten pulses in each set due to the recharging demands of the capacitors. Even though the starting temperature was set to 33°C, we scaled the resulting initial temperatures to match the experimental values in order to provide a more objective comparison. From this figure, it is clear that the temperatures calculated with the numerical model were marginally higher than the measured ones. This calibration model was used as the basis for the parametric study since we were able to closely match the experimental and calculated temperature and electrical current.
Parametric Study Model
Volumes (cm3) of tissue treated with IRE for static σ(T) and dynamic σ(E, T) conductivities
σ(T) t= 0
σ(E, T) t= 0
We report the volumes of tissue treated with IRE as those exposed to a minimum electric field of 500 V/cm, which was found to be the IRE threshold in grey matter for similar pulse parameters to those used in this parametric study . Although we used a 500 V/cm threshold for our calculations, other researchers could adapt the numerical model to investigate the results of IRE in other tissues where the threshold could be different. In order to provide insight to the reader, we modeled the electrical conductivity of the tissue with static and dynamic functions. In the static function, σ(T), the electrical conductivity of the tissue was assumed to be homogeneous and dependent only on the temperature. The dynamic function, σ(E , T), incorporated the dependency of the electrical conductivity on temperature and electroporation. Applying 500 V at 0.5, 1, and 4 Hz resulted in IRE treated volumes between 0.179 - 0.182 cm3 for the static function and 0.293 cm3 for the dynamic function. The IRE treated volumes ranged between 0.424 - 0.460 cm3 for the static function and between 0.706 - 0.732 cm3 for the dynamic function when the applied voltage was 1000 V. Finally, applying 1500 V generated IRE volumes between 0.683 - 0.835 cm3 for the static models and between 1.134 - 1.296 cm3 for the dynamic ones. The IRE treatment volume increased 55% - 69% when the dynamic conductivity function was incorporated as compared to the static conductivity function. The results show the importance of using a conductivity function that takes into account all the relevant physical phenomena that occurs during electroporation in order to provide accurate treatment planning [22, 30–32]. Researchers and physicians should be aware of the increase in treatment volumes due to electroporation and temperature based conductivity changes when performing treatment planning as has been described by several groups in the field [21, 30–33]. Other groups have developed algorithms that are capable of determining optimum electrode configuration and optimum amplitude of the electric pulses for treatment planning of electroporation-based therapies [52, 53].
Electrical Current Distribution
Irreversible Electroporation v. Thermal Effects
Traditional Thermal Assessment
Thermal Damage Assessment
An IRE treatment using applied voltages of 500 V and 1000 V did not result in any thermal damage when delivered at 0.5 Hz. For these cases there was sufficient time for the heat to dissipate through conduction and blood perfusion prior to the onset of the following pulse. However, 1500 V pulses delivered at 0.5 Hz resulted in thermal damage in 0.052 cm3 (4.6%) of the IRE treated tissue. In Figure 9B, there are virtually identical IRE treatment volumes for the 1 Hz repetition rate as the 0.5 Hz of Figure 9A, but when applying 1500 V, there is some thermal damage generated within 20 seconds that affects about 0.183 cm3 of tissue, approximately 16% of the IRE volume. Finally, Figure 9C displays the IRE and thermal damage volumes for the 4 Hz treatment. In this case, thermal damage occurs in 0.094 cm3 (12.8%) of the tissue when applying 1000 V, and approximately 29% of the IRE volume is thermally damaged (0.376 cm3) by increasing the voltage to 1500 V. Additionally, if one focuses on the first seconds of the 1500 V, there is also an increase in the IRE lesion volume due to the increase in the temperature, and thus the electric conductivity.
We previously reported on the first experience applying IRE to the deep subcortical white matter of canine brain . In that procedure electrodes were placed under CT-guidance through minimally invasive 1.2 mm diameter burr holes in order to produce a lesion. Temperatures were measured during this procedure, including the location in close proximity to where the lesion was produced. The low temperatures measured by our system confirmed the unique, non-thermal mode of IRE cell death. The ex vivo lesion volume was smaller than that observed from the in vivo MRIs due to elimination of edema as well as brain shrinkage during the fixation process. There was also limited time for the lesion to evolve relative to our previous work since the experimental aim of this study was to perform the procedure deep in the brain and evaluate the thermal effects, and therefore did not include the 3-day survival .
The ability of IRE to focally ablate small volumes of brain tissue in a minimally invasive fashion has significant potential clinical implications for the treatment of brain diseases in which destruction of focal neuroanatomic target is desired, such as some forms of epilepsy or central neuropathic pain syndromes . We have shown in previous studies the ability to safely produce lesions in the grey matter of the brain cortex [21, 22]. However, many of the potential central nervous system targets may reside deep within the brain, including the white matter . Therefore, it is important to show the ability of IRE to produce a lesion deep within the white matter of the brain. To the best of our knowledge, we performed the first report of a CT-guided intracranial IRE treatment, as well as the first showing IRE pulses may be delivered within the deep white matter of the brain without causing significant edema . We believe that the rapid implementation, minimally invasive nature, and precision offered by image guided IRE will be the preferred treatment delivery platform for future applications of this technology in the brain.
It should be noted that the thermal effects are most prevalent closest to the electrodes, where the electric field magnitude is also highest. Therefore, any thermal damage induced by an IRE procedure should occur within the targeted ablation volume and will not eliminate the effectiveness of the treatment. However, IREs unique non-thermal mechanisms are the key to its ability to be implemented in the vicinity of sensitive structures such as blood vessels and major nerves, a major limitation to resection and thermal therapies. Therefore, a comprehensive and accurate understanding of the potential thermal effects and/or damage is essential to ensure maintaining of these advantages and mitigating the challenges associated with thermal therapies.
Based on the electrodes configuration, measured electrical current and temperature in one canine, we developed a parametric study to investigate the effect of pulse frequency on three different applied voltages of 500, 1000, and 1500 V. The parametric study provides a reliable method to develop treatment protocols to ensure the IRE protocol achieves localized cell death independently from thermal damage. The study was based on pulse frequency, and confirms that if pulses are delivered too rapidly, thermal damage ensues and many of the benefits from this technology will not be optimized for the patient treatment.
The described method of this study takes pulse parameters (frequency, magnitude, and number of pulses) into account in addition to the dynamic changes in tissue electrical conductivity due to temperature increase as well as electroporation. Furthermore, the model accounts for the biological processes of the Pennes' Bioheat Equation, including metabolic heat generation and blood perfusion. Several researchers have demonstrated that blood perfusion is compromised after electroporation in organs outside the central nervous system, thus the heat dissipation from blood convection will be reduced, and it becomes even more important to decrease the frequency of pulse delivery in order to allow for thorough heat dissipation through conduction [55–57]. Even though the effects of pulse duration, electrode exposure, and separation distance were not explicitly investigated in this manuscript, the method can be readily adapted in order to select protocols that do not generate thermal damage in superposition with IRE. Thus, it is necessary that models are explored for each particular application in order to optimize the treatment protocol and better predict the treatment outcome.
Several values have been reported in the literature describing the amount that the electric conductivity of brain tissue increases per degree Celsius [39, 58, 59]. To be conservative, we selected 3.2%°C-1 as the temperature coefficient in this study as reported by Duck et al. . This value is higher than other reported values in the literature that range between 1.4 - 2.0%°C-1 [58, 59], resulting in higher calculated temperatures for our models. In order to assess the effect of a lower temperature coefficient, we simulated treatments with a 1.6%°C-1 value since it is half of the magnitude used in the parametric study and it is still within the range reported in the literature. The volumes of tissue treated with IRE at the onset of the pulses were identical to the ones reported in Table 3 with the dynamic function. At the completion of the pulses, the 1000 and 1500 V applied voltages resulted in smaller volumes of tissue treated with IRE compared to the values reported in Table 3. As with the 3.2%°C-1 temperature coefficient, there were no significant increases in the predicted IRE treatment volume for the 500 V trials. Applying 1000 and 1500 V resulted in IRE treatment volume increases of 1.26 - 3.27% and 1.91 - 8.85% versus those calculated at the onset of the pulses, respectively (before thermal effects begin). The increase in electric conductivity and thus IRE treatment volumes due to the thermal effects are moderate compared to the effects of electroporation. Nevertheless, the electric conductivity dependence on temperature must be incorporated from a thermal perspective in order to optimize IRE protocols, while minimizing any potential thermal damage. IRE is an emerging focal ablation technique and it is vital for researchers and physicians to work together in developing the numerical models for predictable treatment planning. The models presented here provide insight to the role of electroporation and temperature in the resulting volumes of tissue ablated with IRE alone or in superposition with thermal damage. The aim of this work was to provide the reader with numerical methods capable of evaluating pulse parameters used clinically to maximize the benefits of a non-thermal mode of tissue ablation. The numerical methods presented are capable of delineating volumes of tissue undergoing IRE from volumes undergoing thermal damage as a function of time. In this manner, the time point at which different treatment protocols achieve IRE while preventing thermal damage can be determined. It is important for researchers and physicians to be aware of the upper limit of IRE in order to maximize the benefits of a non-thermal model of tissue ablation. Future work should correlate the electric field distribution from these numerical models with reconstructed IRE lesions as seen in MRI and histopathology in order to generate an electric field threshold for brain tissue for clinical use. Future investigations should also determine the electrical magnitudes at which the increase in electrical conductivity occurs for grey matter, white matter, and pathologic brain tissue.
We present the results of a parametric study in brain tissue that investigates 3 voltages delivered at 3 different frequencies which have been used clinically in other tissues such as prostate, kidney, and lungs [24–27]. These numerical simulations were based on an in vivo experimental procedure where a lesion was produced in the white matter of brain . The procedure was performed in a minimally invasive fashion through 1.2 mm diameter burr holes and electrode placement was confirmed with CT imaging . For the first time, the current and temperature were measured together in real-time during the delivery of the pulses and were used together as the basis for the numerical models. The models included all relevant pulse parameters and dynamic changes during treatment, and were capable of determining whether the lesions occurred due to IRE alone or in superposition with thermal damage. IRE alone allows preservation of the major vasculature, extracellular matrix and other critical structures, while achieving cell death in a target location. We hope our results provide physicians and researchers a way to assess individual protocols in order to capitalize on the benefits of this non-thermal mode of tissue ablation.
This work was supported in part by the Coulter Foundation. PAG and RVD acknowledge support from NSF CBET-0933335 towards the modeling effort presented in the manuscript. The authors thank Gregory Daniel, Carolina Ricco, John Robertson, Dana Calicott, Barbara Kafka and Stephanie Milburn for their assistance in surgery. We also thank Barbara Wheeler for her help in processing the brain tissue and John Pearce for his guidance in numerical modeling.
- Al-Sakere B, Andre F, Bernat C, Connault E, Opolon P, Davalos RV, Rubinsky B, Mir LM: Tumor ablation with irreversible electroporation. PLoS ONE 2007, 2(11):e1135. 10.1371/journal.pone.0001135View ArticleGoogle Scholar
- Edd JF, Horowitz L, Davalos RV, Mir LM, Rubinsky B: In vivo results of a new focal tissue ablation technique: irreversible electroporation. IEEE Trans Biomed Eng 2006, 53(7):1409–15. 10.1109/TBME.2006.873745View ArticleGoogle Scholar
- Lee EW, Loh CT, Kee ST: Imaging guided percutaneous irreversible electroporation: ultrasound and immunohistological correlation. Technol Cancer Res Treat 2007, 6(4):287–94.View ArticleGoogle Scholar
- Davalos RV, Mir LM, Rubinsky B: Tissue ablation with irreversible electroporation. Ann Biomed Eng 2005, 33(2):223–31. 10.1007/s10439-005-8981-8View ArticleGoogle Scholar
- Rubinsky B: Irreversible Electroporation in Medicine. Technology in Cancer Research and Treatement 2007, 6(4):255–260.View ArticleGoogle Scholar
- Miklavcic D, Semrov D, Mekid H, Mir LM: A validated model of in vivo electric field distribution in tissues for electrochemotherapy and for DNA electrotransfer for gene therapy. Biochimica et Biophysica Acta 2000, 1523(1):73–83.View ArticleGoogle Scholar
- Miklavcic D, Beravs K, Semrov D, Cemazar M, Demsar F, Sersa G: The importance of electric field distribution for effective in vivo electroporation of tissues. Biophysical Journal 1998, 74(5):2152–2158. 10.1016/S0006-3495(98)77924-XView ArticleGoogle Scholar
- Edd JF, Davalos RV: Mathematical modeling of irreversible electroporation for treatment planning. Technology in Cancer Research and Treatment 2007, 6: 275–286.View ArticleGoogle Scholar
- Guo Y, Zhang Y, Klein R, Nijm GM, Sahakian AV, Omary RA, Yang GY, Larson AC: Irreversible electroporation therapy in the liver: longitudinal efficacy studies in a rat model of hepatocellular carcinoma. Cancer Res 2010, 70(4):1555–63. 10.1158/0008-5472.CAN-09-3067View ArticleGoogle Scholar
- Neal RE II, Singh R, Hatcher HC, Kock ND, Torti SV, Davalos RV: Treatment of breast cancer through the application of irreversible electroporation using a novel minimally invasive single needle electrode. Breast Cancer Res Treat 2010, 123(1):295–301. 10.1007/s10549-010-0803-5View ArticleGoogle Scholar
- Garcia PA, Pancotto T, Rossmeisl JH, Henao-Guerrero N, Gustafson NR, Daniel GB, Robertson JL, Ellis TL, Davalos RV: Non-thermal irreversible electroporation (NTIRE) and adjuvant fractionated radiotherapeutic multimodal therapy for intracranial malignant glioma in a canine patient. Technol Cancer Res Treat 2011, 10(1):73–83.Google Scholar
- Neal RE, Rossmeisl JH, Garcia PA, Lanz OI, Henao-Guerrero N, Davalos RV: Successful Treatment of a Large Soft Tissue Sarcoma With Irreversible Electroporation. J Clin Oncol 2011, 14: 14.Google Scholar
- Davalos RV, Rubinsky B, Mir LM: Theoretical analysis of the thermal effects during in vivo tissue electroporation. Bioelectrochemistry 2003, 61(1–2):99–107. 10.1016/j.bioelechem.2003.07.001View ArticleGoogle Scholar
- Davalos RV, Rubinsky B: Temperature considerations during irreversible electroporation. International Journal of Heat and Mass Transfer 2008, 51(23–24):5617–5622. 10.1016/j.ijheatmasstransfer.2008.04.046View ArticleMATHGoogle Scholar
- Tropea BI, Lee RC: Thermal injury kinetics in electrical trauma. J Biomech Eng 1992, 114(2):241–50. 10.1115/1.2891378View ArticleGoogle Scholar
- Lee RC, Zhang D, Hannig J: Biophysical injury mechanisms in electrical shock trauma. Annual Review of Biomedical Engineering 2000, 2(1):477–509. 10.1146/annurev.bioeng.2.1.477View ArticleGoogle Scholar
- Becker SM, Kuznetsov AV: Numerical modeling of in vivo plate electroporation thermal dose assessment. Journal of Biomechanical Engineering 2006, 128(1):76–84. 10.1115/1.2132375View ArticleGoogle Scholar
- Shafiee H, Garcia PA, Davalos RV: A Preliminary Study to Delineate Irreversible Electroporation From Thermal Damage Using the Arrhenius Equation. Journal of Biomechanical Engineering 2009, 131(7):074509. 10.1115/1.3143027View ArticleGoogle Scholar
- Lackovic I, Magjarevic R, Miklavcic D: Three-dimensional finite-element analysis of joule heating in electrochemotherapy and in vivo gene electrotransfer. IEEE Transactions on Dielectrics and Electrical Insulation 2009, 16(5):1338–1347.View ArticleGoogle Scholar
- Pliquett U: Joule heating during solid tissue electroporation. Med Biol Eng Comput 2003, 41(2):215–9. 10.1007/BF02344892View ArticleGoogle Scholar
- Ellis TL, Garcia PA, Rossmeisl JH, Henao-Guerrero N, Robertson J, Davalos RV: Nonthermal irreversible electroporation for intracranial surgical applications. J Neurosurg 2011, 114: 681–688. 10.3171/2010.5.JNS091448View ArticleGoogle Scholar
- Garcia PA, Rossmeisl JH, Neal RE II, Ellis TL, Olson J, Henao-Guerrero N, Robertson J, Davalos RV: Intracranial nonthermal irreversible electroporation: In vivo analysis. J Membr Biol 2010, 236(1):127–136. 10.1007/s00232-010-9284-zView ArticleGoogle Scholar
- Garcia PA, Neal RE II, Rossmeisl JH, Davalos RV: Non-thermal irreversible electroporation for deep intracranial disorders. Conf Proc IEEE Eng Med Biol Soc 2010, 1: 2743–6.Google Scholar
- Lee EW, Thai S, Kee ST: Irreversible Electroporation: A Novel Image-Guided Cancer Therapy. Gut Liver 2010, 4(Suppl 1):S99-S104.View ArticleGoogle Scholar
- Ball C, Thomson KR, Kavnoudias H: Irreversible electroporation: a new challenge in "out of operating theater" anesthesia. Anesth Analg 2010, 110(5):1305–9. 10.1213/ANE.0b013e3181d27b30View ArticleGoogle Scholar
- Pech M, Janitzky A, Wendler JJ, Strang C, Blaschke S, Dudeck O, Ricke J, Liehr UB: Irreversible Electroporation of Renal Cell Carcinoma: A First-in-Man Phase I Clinical Study. Cardiovasc Intervent Radiol 2010.Google Scholar
- Onik G, Rubinsky B: Irreversible Electroporation. In Irreversible Electroporation: First Patient Experience Focal Therapy of Prostate Cancer. Edited by: Onik G, Rubinsky B. Springer Berlin Heidelberg; 2010:235–247.View ArticleGoogle Scholar
- Garcia PA, Rossmeisl JH, Robertson J, Ellis TL, Davalos RV: Pilot study of irreversible electroporation for intracranial surgery. Conf Proc IEEE Eng Med Biol Soc 2009, 1: 6513–6.Google Scholar
- Latikka J, Kuurne T, Eskola H: Conductivity of living intracranial tissues. Phys Med Biol 2001, 46(6):1611–6. 10.1088/0031-9155/46/6/302View ArticleGoogle Scholar
- Davalos RV, Otten DM, Mir LM, Rubinsky B: Electrical impedance tomography for imaging tissue electroporation. IEEE Transactions on Biomedical Engineering 2004, 51(5):761–767. 10.1109/TBME.2004.824148View ArticleGoogle Scholar
- Ivorra A, Rubinsky B: In vivo electrical impedance measurements during and after electroporation of rat liver. Bioelectrochemistry 2007, 70(2):287–95. 10.1016/j.bioelechem.2006.10.005View ArticleGoogle Scholar
- Sel D, Cukjati D, Batiuskaite D, Slivnik T, Mir LM, Miklavcic D: Sequential finite element model of tissue electropermeabilization. IEEE Trans Biomed Eng 2005, 52(5):816–27. 10.1109/TBME.2005.845212View ArticleGoogle Scholar
- Ivorra A, Al-Sakere B, Rubinsky B, Mir LM: In vivo electrical conductivity measurements during and after tumor electroporation: conductivity changes reflect the treatment outcome. Phys Med Biol 2009, 54(19):5949–63. 10.1088/0031-9155/54/19/019View ArticleGoogle Scholar
- Miklavcic D, Snoj M, Zupanic A, Kos B, Cemazar M, Kropivnik M, Bracko M, Pecnik T, Gadzijev E, Sersa G: Towards treatment planning and treatment of deep-seated solid tumors by electrochemotherapy. BioMedical Engineering OnLine 2010, 9(1):10. 10.1186/1475-925X-9-10View ArticleGoogle Scholar
- Pavselj N, Bregar Z, Cukjati D, Batiuskaite D, Mir LM, Miklavcic D: The course of tissue permeabilization studied on a mathematical model of a subcutaneous tumor in small animals. IEEE Trans Biomed Eng 2005, 52(8):1373–81. 10.1109/TBME.2005.851524View ArticleGoogle Scholar
- Sel D, Lebar AM, Miklavcic D: Feasibility of employing model-based optimization of pulse amplitude and electrode distance for effective tumor electropermeabilization. IEEE Trans Biomed Eng 2007, 54(5):773–81.View ArticleGoogle Scholar
- Cukjati D, Batiuskaite D, Andre F, Miklavcic D, Mir LM: Real time electroporation control for accurate and safe in vivo non-viral gene therapy. Bioelectrochemistry 2007, 70(2):501–7. 10.1016/j.bioelechem.2006.11.001View ArticleGoogle Scholar
- Becker SM, Kuznetsov AV: Thermal damage reduction associated with in vivo skin electroporation: A numerical investigation justifying aggressive pre-cooling. International Journal of Heat and Mass Transfer 2007, 50: 105–116. 10.1016/j.ijheatmasstransfer.2006.06.030View ArticleMATHGoogle Scholar
- Duck FA: Physical Properties of Tissues: A Comprehensive Reference Book. San Diego: Academic Press; 1990.Google Scholar
- Werner J, Buse M: Temperature profiles with respect to inhomogeneity and geometry of the human body. J Appl Physiol 1988, 65(3):1110–8.Google Scholar
- Uzuka T, Tanaka R, Takahashi H, Kakinuma K, Matsuda J, Kato K: Planning of hyperthermic treatment for malignant glioma using computer simulation. Int J Hyperthermia 2001, 17(2):114–22. 10.1080/02656730010010676View ArticleGoogle Scholar
- Cosman ER Jr, Cosman ER Sr: Electric and thermal field effects in tissue around radiofrequency electrodes. Pain Med 2005, 6(6):405–24. 10.1111/j.1526-4637.2005.00076.xView ArticleGoogle Scholar
- Diller KR: Modeling of bioheat transfer processes at high and low temperatures. In Bioengineering heat transfer. Volume 32. Edited by: Choi YI. Boston: Academic Press, Inc; 1992:157–357.View ArticleGoogle Scholar
- Sherar MD, Moriarty JA, Kolios MC, Chen JC, Peters RD, Ang LC, Hinks RS, Henkelman RM, Bronskill MJ, Kucharcyk W: Comparison of thermal damage calculated using magnetic resonance thermometry, with magnetic resonance imaging post-treatment and histology, after interstitial microwave thermal therapy of rabbit brain. Phys Med Biol 2000, 45(12):3563–76. 10.1088/0031-9155/45/12/304View ArticleGoogle Scholar
- Henriques FC, Moritz AR: Studies in thermal injuries: the predictability and the significance of thermally induced rate processes leading to irreversible epidermal damage. Arch Pathol 1947, 43: 489–502.Google Scholar
- Diller KR, Hayes LJ: A finite element model of burn injury in blood-perfused skin. J Biomech Eng 1983, 105(3):300–7. 10.1115/1.3138423View ArticleGoogle Scholar
- Jiang SC, Ma N, Li HJ, Zhang XX: Effects of thermal properties and geometrical dimensions on skin burn injuries. Burns 2002, 28(8):713–7. 10.1016/S0305-4179(02)00104-3View ArticleGoogle Scholar
- Brown SL, Hunt JW, Hill RP: Differential thermal sensitivity of tumor and normal tissue microvascular response during hyperthermia. International Journal of Hyperthermia 1992, 8(4):501–514. 10.3109/02656739209037988View ArticleGoogle Scholar
- Borrelli MJ, Thompson LL, Cain CA, Dewey WC: Time-temperature analysis of cell killing of BHK cells heated at temperatures in the range of 43.5 degrees C to 57.0 degrees C. International Journal of Radiation Oncology Biology Physics 1990, 19(2):389–399. 10.1016/0360-3016(90)90548-XView ArticleGoogle Scholar
- Jacques SL, Newman C, He XY: Thermal coagulation of tissues: liver studies indicate a distribution of rate parameters, not a single rate parameter, describes the coagulation process. Advances in Biological Heat and Mass Transfer by New York: ASME: J McGrath 1991, 71–73.Google Scholar
- Neal RE II, Davalos RV: The feasibility of irreversible electroporation for the treatment of breast cancer and other heterogeneous systems. Ann Biomed Eng 2009, 37(12):2615–25. 10.1007/s10439-009-9796-9View ArticleGoogle Scholar
- Kos B, Zupanic A, Kotnik T, Snoj M, Sersa G, Miklavcic D: Robustness of treatment planning for electrochemotherapy of deep-seated tumors. J Membr Biol 2010, 236(1):147–53. 10.1007/s00232-010-9274-1View ArticleGoogle Scholar
- Zupanic A, Corovic S, Miklavcic D: Optimization of electrode position and electric pulse amplitude in electrochemotherapy. Radiol Oncol 2008, 42(2):93–101. 10.2478/v10019-008-0005-5View ArticleGoogle Scholar
- Nicolaidis S: Neurosurgical treatments of intractable pain. Metabolism 2010, 59(Suppl 1):S27–31.View ArticleGoogle Scholar
- Jarm T, Cemazar M, Miklavcic D, Sersa G: Antivascular effects of electrochemotherapy: implications in treatment of bleeding metastases. Expert Review of Anticancer Therapy 2010, 10(5):729–746. 10.1586/era.10.43View ArticleGoogle Scholar
- Sersa G, Cemazar M, Miklavcic D, Chaplin DJ: Tumor blood flow modifying effect of electrochemotherapy with bleomycin. Anticancer Research 1999, 19(5B):4017–4022.Google Scholar
- Sersa G, Jarm T, Kotnik T, Coer A, Podkrajsek M, Sentjurc M, Miklavcic D, Kadivec M, Kranjc S, Secerov A, Cemazar M: Vascular disrupting action of electroporation and electrochemotherapy with bleomycin in murine sarcoma. British Journal of Cancer 2008, 98(2):388–398. 10.1038/sj.bjc.6604168View ArticleGoogle Scholar
- Pethig R: Dielectric Properties of Biological Materials: Biophysical and Medical Applications. IEEE Transactions on Electrical Insulation 1984, EI-19(5):453–474.View ArticleGoogle Scholar
- Foster KR, Schwan HP: Dielectric properties of tissues and biological materials: a critical review. Crit Rev Biomed Eng 1989, 17(1):25–104.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.