Skip to content
BY-NC-ND 4.0 license Open Access Published by De Gruyter Open Access December 29, 2017

Electric field computation and measurements in the electroporation of inhomogeneous samples

  • Alessia Bernardis , Marco Bullo , Luca Giovanni Campana , Paolo Di Barba , Fabrizio Dughiero , Michele Forzan , Maria Evelina Mognaschi EMAIL logo , Paolo Sgarbossa and Elisabetta Sieni
From the journal Open Physics

Abstract

In clinical treatments of a class of tumors, e.g. skin tumors, the drug uptake of tumor tissue is helped by means of a pulsed electric field, which permeabilizes the cell membranes. This technique, which is called electroporation, exploits the conductivity of the tissues: however, the tumor tissue could be characterized by inhomogeneous areas, eventually causing a non-uniform distribution of current. In this paper, the authors propose a field model to predict the effect of tissue inhomogeneity, which can affect the current density distribution. In particular, finite-element simulations, considering non-linear conductivity against field relationship, are developed. Measurements on a set of samples subject to controlled inhomogeneity make it possible to assess the numerical model in view of identifying the equivalent resistance between pairs of electrodes.

1 Introduction

Electrochemotherapy (ECT) is a local cancer therapy that uses voltage pulses. Voltage pulses generate an electric field that reversibly permeabilizes cell membranes. This way the chemotherapeutic drugs uptake is improved [1, 2, 3]. This therapy is currently applied to skin tumors, such as melanomas, or chest wall recurrence of breast cancer [4, 5].

The electric field is applied to the patient according to standard procedures, by means of needle electrode, formed e.g. by 7 or 8 needles, that covers a surface of about 1-2 cm2 [2, 3]. ECT standard protocols uses 8 voltage pulses per electrode pair (rectangular pulses, 100 µs long at 5 kHz) in order to generate an electric field. The local strength of the electric field in the tissue depends on tissue inhomogeneity. In fact, in the treatment area the electric field intensity is homogeneous if electrical properties of the tissue are homogeneous. Conversely, if the electrical properties of the tissue are inhomogeneous, the electric field intensity depends on the local electrical resistivity. This behavior is not unusual in tumor tissues where both myxoid and fibrotic components can be found in the tumor area [6, 7]. Moreover, in clinical practice the two needles of a pair can be inserted in tissues with different electrical resistivity.

This paper systematically investigates the problem of the discontinuity of the tissue electrical properties by means of numerical simulations and experiments on suitable phantoms. Phantoms are used to design different electrical properties, varying the electrical resistivity. Finally, the electrical resistance is evaluated both experimentally, starting from voltage and current measurements, and using Finite Element simulations.

2 Measurement procedure

The inhomogeneous phantom used in experiments is obtained using materials with different electrical resistivity. One of these materials is a slice of potato. This tuber is used in order to show the electroporation occurrence. In fact, it is well known that potato becomes dark just few hours after cell electroporation [8, 9, 10, 11, 12] if the tuber is preserved at room temperature. Using voltage and current data recorded by means of the pulse generator, the sample resistance is computed [6, 7]. The results obtained with the experimental tests are compared with numerical results obtained using a 3D model.

2.1 Simulation model

The numerical model used in simulations is shown in Figure 1, and is composed of an internal parallelepiped (1×3.5×5 cm) divided into two sub-volumes immersed in a second parallelepiped (1×5.5×7 cm). The electric field is generated between two cylindrical electrodes 9 mm long, with diameter 0.6 mm, separated by d =7.3 mm. The two internal sub-volumes are characterized by an electrical resistivity ρi(E) dependent on the tissue type. The two needle electrodes are positioned in the potato sample and in the gel, respectively. A voltage of 730 V is applied between the needle pairs by means of Cliniporator EPS01-3.x version (Igea Clinical Biophysics S.p.A., Carpi, Modena Italy). The generator load can be designed as a resistive network schematized in Figure 2.

Figure 1 3D geometries: (a) three volumes with needles in potato and gel and (b) two volumes with both the needles in potato.
Figure 1

3D geometries: (a) three volumes with needles in potato and gel and (b) two volumes with both the needles in potato.

Figure 2 Electrical schema of the model.
Figure 2

Electrical schema of the model.

The electric field in the different volumes of the model is computed using Finite Element analysis solving a static conduction problem [13, 14, 15, 16, 17, 18]. Dirichlet conditions were imposed to each electrode’s surfaces [14, 16, 17, 18], i.e. constant electric potentials V1 = –V/2 and V2 = +V/2, respectivelly.

In order to obtain more accurate results [19, 20], the numerical model was solved considering potato resistivity dependent on electric field strength [10]:

ρ(E)=(σ0+σ1σ021+tanhkvEEth)1(1)

The typical range for conductivity is between 0.18 and 0.97 S/m, values measured on ‘electropored’ and ‘not electropored’ potatoes, respectively. Eth was evaluated experimentally.

The field problem was solved by means of a commercial software based on Finite Elements. A field model is of utmost importance when field quantities have to be evaluated and quantified precisely [21, 22].

2.2 Gel phantom preparation

The gel phantoms were produced according to a slightly modified procedure based on that proposed in [20] for the preparation of low-cost tissue-mimicking materials (spinal cord, SC and dura, D). Gelatin, water, agar, corn flour, potato flour, glycerin, sodium azide (NaN3) and sodium chloride (NaCl) were used as received. The list and amount of starting ingredients for the production of the phantoms is reported in Table 1 [23, 24].

Table 1

List of ingredients for the preparation of phantom materials.

D2DF1SC1
Glycerin2.00 g2.00 g0.90 g
Corn flour25.0 g-30.0 g
Potato flour-25.0g-
Gelatin--0.60 g
Agar1.00 g1.00 g-
Sodium azide0.80 g0.40 g0.30 g
Sodium chloride0.20g0.20 g-
Water(i) 20 ml + (ii) 50 ml

The procedure for the preparation of the tissue-mimic material follows four steps [23]:

  1. In a 50 ml beaker, the glycerin is added to 20 ml of deionized water at room temperature, then the corn or potato flour is slowly added while stirring vigorously;

  2. In an 800 ml beaker, the gelatin or agar is heated gradually to 90° C in a water bath with 50 ml of deionized water and NaCl. Heating is maintained until gelatin or agar is dissolved into a thick mixture;

  3. The viscous mixture from step (i) is added to the one obtained from step (ii). Heating and stirring are applied until the whole mixture turns semi-solid;

  4. After heating, the mixture is cooled to 50° C and NaN3 is added with mechanical stirring. The semisolid material is cast into boxes with suitable sizes for the experiments.

2.3 Material resistivity identification

The electrical resistivity, ρ, of the gel phantoms and potato tubers was evaluated experimentally using a plastic box with rectangular section (11 x× 9 mm, thickness 11.3 mm) and an electrode with two plate at L = 7.5 mm with rectangular geometry (width 10 mm) as shown in Figure 3. The voltage was applied to the electrode plates by means of a voltage pulse generator (manufactured by Igea, Carpi, Italy) and the corresponding current is measured. From voltage and current data the resistance was computed as in [7] and the resistivity was evaluated using Ohm’s law and the model for resistance, R, of a parallelepiped sample with length L (7.5 mm) and section A (10 × 11.3 mm):

ρ=RAL(2)
Figure 3 Experimental set-up for resistivity.
Figure 3

Experimental set-up for resistivity.

In order to evaluate gel resistivity the applied voltage pulses had an amplitude of 100 V (E = 130 V/cm) and 760 V (E = 1000 V/cm)).

2.4 Data analysis

Voltage and current data recorded by ECT equipment were used in order to evaluate the equivalent resistance at the cylindrical needle ends as in [7]. For the sake of a comparison, considering the computation model in Figure 1, the resistance at the electrode was evaluated as the ratio between imposed voltage V and the integral of the normal component of current density J evaluated on the transversal section, S, of the current density flux

R=V/SJ¯n¯dS(3)

Where J depends on the electric field and material conductivity, J = σE. This way the equivalent resistance at the needle ends is computed and depends on the electric field if the conductivity depends on electric field [10].

2.5 Evaluation of potato electroporation threshold

The electric field intensity used to obtain electroporation on the potato tuber was evaluated experimentally using the set-up in Figure 3. Potato cubes were positioned in the plastic support. For each sample, 8 voltage pulses at 5 kHz were applied. Different voltage amplitudes between 0 and 800 V (electric field between 0 to 1060 V/cm) were selected for each sample. Sample resistance and resistivity were computed from the voltage and current data. The potato color was observed 24 h after the pulse application. Samples were preserved at room temperature, by covering them with plastic film. The color of each potato sample was compared to its sample resistivity. Dark intensity was evaluated by means of ImageJ. For each image the median dark intensity was recorded. Experiments were repeated at least thrice.

2.6 Computational and experimental set-up

Using the models in Figure 1 three types of set-up of materials, potato and gels, were used in computational and experimental tests: (a) model in Figure 1(b) with volume Vp filled with potato, volume Vg with one gel from Table 2; (b) model in Figure 1(a) with volume V1 filled with potato, volume V2 with one gel from Table 2 and volume V3 empty; (c)model in Figure 1(a) with volume V1 filled with potato, volume V2 and V3 with one gel from Table 2.

Table 2

Resistivity and conductivity values for gels

(@130V/cm)(@1000V/cm)
Tissueρ[Ωm]σ[S/m]ρ[Ωm]σ[S/m]
D20.78 ± 0.091.30 ± 0.140.83 ± 0.031.20 ± 0.5
DF11.09 ± 0.020.92 ± 0.010.99 ± 0.061.01 ± 0.06
SC13.38 ± 0.060.3 ± 0.013.33 ± 0.20.3 ± 0.02

3 Results and discussion

Computational and experimental results are discussed below.

3.1 Gel electrical characteristics

Table 2 reports the gel resistivity and conductivity computed at two different electric field intensities using equation (3). The values evaluated from applying a uniform electric field at 130 V/cm and 1000 V/cm are comparable.

3.2 Potato characterization

Table 3 reports the resistivity for potato tuber evaluated by means of equation (3) varying the applied voltage. Each value is the average of at least 3 values. From these values of resistivity the conductivity values for electroporated (ρ1 = 1.05 Ωm, σ1 = 0.97 S/m) and not-electroporated (ρ0 = 6.6 Ωm, σ0 = 0.18 S/m) potato used in the equation (1) were evaluated. The average of dark intensity (DI) is also reported. In this case a scale between 0 and 255 was used, where a value close to 0 corresponds to a dark color.

Figure 4 Evaluation of electric field threshold for electroporation: (a) fit of experimental data and (b) potato color as a function of electric field.
Figure 4

Evaluation of electric field threshold for electroporation: (a) fit of experimental data and (b) potato color as a function of electric field.

Table 3

Potato resistivity, conductivity, and dark intensity (DI) at different voltage amplitude (electric field).

VE[V/cm]ρ[Ωm]σ[S/m]DI
00115.7 ±5.7
1001306.60 ± 2.790.18 ± 0.0798.1 ±11.4
1502004.12 ± 0.890.25 ± 0.0587.1 ±12.2
2002652.29 ± 0.500.46 ± 0.1249.9 ± 5.9
4005301.52 ± 0.210.67 ± 0.0950.6 ± 8.8
5006601.47 ± 0.630.78 ± 0.2654.7 ± 7.4
6008001.14 ± 0.260.92 ± 0.2054.9 ± 5.3
7009301.12 ± 0.230.93 ± 0.1940.7 ± 5.0
80010601.05 ± 0.170.97 ± 0.1335.0 ± 5.0

3.3 Inhomogeneous models

Figure 5 shows an example of the conduction field map obtained by attributing constant resistivity values to gel (data in Table 2) and potato. In this case, non-linear conductivity was considered for the potato samples. The field asymmetry is evident. Figure 6 shows the corresponding potato pieces obtained in experimental validation. The differences in the amplitudes of the dark areas at the boundaries with the gel (right side of each piece) are evident. They correspond to differences in the positions of the 390 V/cm equi-field level (orange line in Figure 5) in the potato domain in the simulation results. In particular, in the DF1 case, the 390 V/cm line at the boundary between the potato and the gel has a diameter of about 0.75 cm, while in the SC1 case it has a diameter of about 0.35 cm. Considering the electric field level at 265 V/cm, i.e. the level at which electroporation starts, the level line diameters are close to 1 cm and 1.5 cm, respectively, similar to the experimental values in Table 4.

Figure 5 Conduction field map considering a non-linear conductivity for potato, σ1, and constant for gels. σ2= σ3 equal to (a) 0.3 S/m and (b) 1 S/m.
Figure 5

Conduction field map considering a non-linear conductivity for potato, σ1, and constant for gels. σ2= σ3 equal to (a) 0.3 S/m and (b) 1 S/m.

Figure 6 Potato 24 h after pulse application In the model with three materials, as shown In Figure 1(a), using gels (a) DF1 and (b) SC1.
Figure 6

Potato 24 h after pulse application In the model with three materials, as shown In Figure 1(a), using gels (a) DF1 and (b) SC1.

Table 4

Resistance simulated and measured at electrode extremities considering different set-up, DI and distance XA

Needle positionRs [Ω]Rm [Ω]XA [cm]DI
(a) Two needles in potato
Big potato140180.0 ± 28.01.33 ± 0.1650
Medium potato142165.8 ± 25.11.11 ± 0.2454
Small potato145178.6 ± 6.21.10 ± 0.1758
(a1) Two needles in potato as in Figure 1(b)
Potato + D2136172.0 ± 14.41.55 ± 0.0845
Potato + DF1137170.6 ± 19.11.47 ± 0.2856
Potato + SC1141180.5 ± 10.21.37 ± 0.1264
(b) One needle in potato, one in gel
Potato and D2106145.1 ± 2.91.80 ± 0.0644
Potato and DF1120145.8 ± 9.81.72 ± 0.2649
Potato and SC1261283.1 ± 18.31.13 ± 0.3466
(c) One needle in potato, one in gel + surrounding gel
Potato and D2 + D2104133.9 ± 18.51.55 ± 0.0949
Potato and DF1 + DF1121143.3 ± 14.61.74 ± 0.3252
Potato and SC1 + SC1264276.7 ± 16.71.06 ± 0.2870

Table 4 includes findings of the tests performed using different combinations of gels and potato subject to a voltage of 730 V. It reports the resistance evaluated by means of (3) from computational models, and the experimental resistance. Moreover, the amplitude of the dark area at the boundary with the gel (right side of the potato pieces) for two- or three-material models, or in the middle between the two needles, is reported.

4 Conclusions

Tissue inhomogeneity is a relevant characteristic of tumors. From this background, the proposed Finite Elements model takes into account tissue inhomogeneities. The experiments were performed on a potato and the simulations are carried out accordingly. A good match between experimental results and simulations was found. In future works the proposed model will be extended to cancer tissues.

Acknowledgement

Authors are grateful to Igea S.p.A. (Carpi(MO) Italy), for pulse generator and electrodes. Project granted by CPDA138001 (Padua University) and partially made possible thanks to COST TD1104 action (www.electroporation.net).

References

[1] Mir L.M., Orlowski S., Mechanisms of electrochemotherapy, Advanced Drug Delivery Reviews, 1999, 35, 107–118.10.1016/S0169-409X(98)00066-0Search in Google Scholar PubMed

[2] Mir L.M., Gehl J., Sersa G., Collins C.G., et al., Standard operating procedures of the electrochemotherapy: Instructions for the use of bleomycin or cisplatin administered either systemically or locally and electric pulses delivered by the CliniporatorTM by means of invasive or non-invasive electrodes, EJC Supplements, 2006, 4, 14–25.10.1016/j.ejcsup.2006.08.003Search in Google Scholar

[3] Marty M., Sersa G., Garbay J.R., Gehl J., et al., Electrochemotherapy – An easy, highly effective and safe treatment of cutaneous and subcutaneous metastases: Results of ESOPE (European Standard Operating Procedures of Electrochemotherapy) study, European Journal of Cancer Supplements, 2006, 4, 3–13.10.1016/j.ejcsup.2006.08.002Search in Google Scholar

[4] Campana L., Mocellin S., Basso M., Puccetti O., et al., Bleomycin-Based Electrochemotherapy: Clinical Outcome from a Single Institution’s Experience with 52 Patients, Annals of Surgical Oncology, 2009, 16, 191–199.10.1245/s10434-008-0204-8Search in Google Scholar PubMed

[5] Sersa G., Cufer T., Paulin S.M., Cemazar M., et al., Electrochemotherapy of chest wall breast cancer recurrence. Cancer Treat. Rev, 2012, 38, 379–386.10.1016/j.ctrv.2011.07.006Search in Google Scholar PubMed

[6] Tosi A.L., Campana L.G., Dughiero F., Forzan M., et al., Microscopic histological characteristics of soft tissue sarcomas: analysis of tissue features and electrical resistance, Medical & Biological Engineering & Computing 2017, 55, 1097–1108.10.1007/s11517-016-1573-ySearch in Google Scholar PubMed

[7] Campana L.G., Cesari M., Dughiero F., Forzan M., et al., Electrical resistance of human soft tissue sarcomas: an ex vivo study on surgical specimens, Medical & Biological Engineering & Computing, 2016, 54, 773–787.10.1007/s11517-015-1368-6Search in Google Scholar PubMed

[8] Hjouj M., Rubinsky B., Magnetic Resonance Imaging Characteristics of Nonthermal Irreversible Electroporation in Vegetable Tissue, J Membrane Biol, 2010, 236, 137–146.10.1007/s00232-010-9281-2Search in Google Scholar PubMed

[9] Breton M., Mir L.M., Microsecond and nanosecond electric pulses in cancer treatments, Bioelectromagnetics 2011, n/a.10.1002/bem.20692Search in Google Scholar PubMed

[10] Breton M., Buret F., Krahenbuhl L., Leguebe M., et al., Non-Linear Steady-State Electrical Current Modeling for the Electropermeabilization of Biological Tissue, IEEE Transactions on Magnetics, 2015, 51, 1–4.10.1109/TMAG.2014.2351836Search in Google Scholar

[11] Ivorra A., Mir L.M., Rubinsky B., Electric Field Redistribution due to Conductivity Changes during Tissue Electroporation: Experiments with a Simple Vegetal Model, In: Dössel, O., Schlegel, W. (Eds.), World Congress on Medical Physics and Biomedical Engineering (7-12 September 2009, Munich, Germany), Springer Berlin Heidelberg, 2010, 59–62.10.1007/978-3-642-03895-2_18Search in Google Scholar

[12] Castellví Q., Banús J., Ivorra A., 3D Assessment of Irreversible Electroporation Treatments in Vegetal Models, In: Jarm, T., Kramar, P. (Eds.), 1st World Congress on Electroporation and Pulsed Electric Fields in Biology, Medicine and Food & Environmental Technologies, Springer Singapore, Singapore, 2016, 294–297.10.1007/978-981-287-817-5_65Search in Google Scholar

[13] Ongaro A., Campana L.G., De Mattei M., Dughiero F., et al., Effect of Electrode Distance in Electrochemotherapy: From Numerical Model to in Vitro Tests, In: Jarm, T., Kramar, P. (Eds.), 1st World Congress on Electroporation and Pulsed Electric Fields in Biology, Medicine and Food & Environmental Technologies, Springer Singapore, 2016, 167–170.10.1007/978-981-287-817-5_37Search in Google Scholar

[14] Ongaro A., Campana L.G., De Mattei M., Dughiero F., et al., Evaluation of the Electroporation Efficiency of a Grid Electrode for Electrochemotherapy: From Numerical Model to In Vitro Tests, Technology in Cancer Research & Treatment, 2016, 15, 296–307.10.1177/1533034615582350Search in Google Scholar PubMed

[15] Corovic S., Lackovic I., Sustaric P., Sustar T., et al., Modeling of electric field distribution in tissues during electroporation, BioMedical Engineering OnLine, 2013, 12, 16.10.1186/1475-925X-12-16Search in Google Scholar PubMed PubMed Central

[16] Pavselj N., Miklavcic D., Numerical Models of Skin Electropermeabilization Taking into Account Conductivity Changes and the Presence of Local Transport Regions, IEEE Transactions on Plasma Science, 2008, 36, 1650–1658.10.1109/TPS.2008.928715Search in Google Scholar

[17] Campana L.G., Di Barba P., Dughiero F., Rossi C.R., et al., Optimal Needle Positioning for Electrochemotherapy: A Constrained Multiobjective Strategy, IEEE Transactions on Magnetics, 2013, 49, 2141–2144.10.1109/TMAG.2013.2241031Search in Google Scholar

[18] Corovic S., Zupanic A., Miklavcic D., Numerical Modeling and Optimization of Electric Field Distribution in Subcutaneous Tumor Treated with Electrochemotherapy Using Needle Electrodes, IEEE Transactions on Plasma Science, 2008, 36, 1665–1672.10.1109/TPS.2008.2000996Search in Google Scholar

[19] Di Barba P., Dolezel I., Mognaschi M.E., Savini A., et al., Non-Linear Multi-Physics Analysis and Multi-Objective Optimization in Electroheating Applications, IEEE Transactions on Magnetics, 2014, 50, 673–676.10.1109/TMAG.2013.2286491Search in Google Scholar

[20] Di Barba P., Dolezel I., Karban P., Kus P., et al., Multiphysics field analysis and multiobjective design optimization: a benchmark problem, Inverse Problems in Science and Engineering, 2014, 22, 1214–1225.10.1080/17415977.2013.860590Search in Google Scholar

[21] Mognaschi M., Field models in low-frequency bioelectromagnetics, PRZEGLĄD ELEKTROTECHNICZNY 2016, 1, 3–6.10.15199/48.2016.12.01Search in Google Scholar

[22] Mognaschi M., Di Barba P., Magenes G., Lenzi A., et al., Field models and numerical dosimetry inside an extremely-low-frequency electromagnetic bioreactor: the theoretical link between the electromagnetically induced mechanical forces and the biological mechanisms of the cell tensegrity, SpringerPlus, 2014, 3, 473.10.1186/2193-1801-3-473Search in Google Scholar PubMed PubMed Central

[23] Campana L.G., Di Barba P., Mognaschi M.E., Bullo M., Dughiero F., Forzan M., et al., Electrical resistance in inhomogeneous samples during electroporation, in: Proceedings of the 14th International Conference on Synthesis, Modeling, Analysis and Simulation Methods and Applications to Circuit Design-SMACD (12-15 June 2017, Taormina, Italy), 2017, 1-4.10.1109/SMACD.2017.7981596Search in Google Scholar

[24] Mobashsher A.T., Abbosh A.M., Three-Dimensional Human Head Phantom With Realistic Electrical Properties and Anatomy, IEEE Antennas and Wireless Propagation Letters, 2014, 13, 1401–1404.10.1109/LAWP.2014.2340409Search in Google Scholar

Received: 2017-11-2
Accepted: 2017-11-12
Published Online: 2017-12-29

© 2017 A. Bernardis et al.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Downloaded on 30.3.2024 from https://www.degruyter.com/document/doi/10.1515/phys-2017-0092/html
Scroll to top button