Submit manuscript...
Journal of
eISSN: 2373-633X

Cancer Prevention & Current Research

Research Article Volume 6 Issue 3

Simulation study of electric field distribution in brain tumors

Sri Harsha Venuturumilli,1 Raji Sundararajan2

1IIT-Madras, India
2Purdue University, USA

Correspondence: Raji Sundararajan, 401 N. Grant St., West Lafayette IN 47907, USA

Received: August 29, 2016 | Published: December 9, 2016

Citation: Venuturumilli SH, Sundararajan R. Simulation study of electric field distribution in brain tumors. J Cancer Prev Curr Res. 2016;6(3):452-454. DOI: 10.15406/jcpcr.2016.06.00203

Download PDF

Abstract

Brain tumors are the second leading cause of cancer-related deaths in children and young adults. About 37% of the brain tumors are malignant, causing brain cancers. Current standard of cure includes surgery, radiotherapy and chemotherapy. Brain surgery is too complex and involves more risks in treatment. In external beam radiation therapy, multiple treatments of standard-dose fractions of radiation are applied to the brain. This process is repeated for a total of 10 to 30 treatments, depending on the type of tumor. Chemotherapy may improve overall survival in patients with the most malignant primary brain tumors. However, it does so in only about 20% of the patients. Chemotherapy is often used in young children instead of radiation, as radiation may have negative effects on the developing brain. Thus, there is a critical need for alternate/new therapies, and Electrochemotheraphy is an attractive alternative. By applying electric field at critical strength across cell membrane, it is shown to make it permeable to anti-cancer drugs. In this study, the electric field distribution of human brain is investigated using commercial, state-of-the-art, COMSOL software that utilizes finite element method, with variable parameters like tissue conductivity and permeability. The electric field distribution on tumor regions at various depths (slices) were examined for two types of electrode geometries, viz., needle and microelectrode configurations. The needle electrodes showed the max intensity of 1600V/cm at the 0 cm slice, while the microelectrodes showed reduced values. These results could be used for clinical applications.

Background

Our body is electric and it is possible to elicit both therapeutic and diagnostic phenomena using electrical pulses of appropriate magnitude and frequency (duration). The use of electrical fields to treat brain cancer is gaining momentum.1,2 Electrical fields were also used in Israel to treat brain tumors successfully, with twice the median time to recur of the cancer and the overall survival period.1 At Yale University, alternating electric fields were used to arrest cell proliferation in human brain tumors.2 This is possible because living cells consist of ions, polar or charged molecules, membranes, and organelles, and they are responsive to electric fields of appropriate magnitudes and frequencies therapeutically.1‒3 The electric activity of cells plays a key role in many essential biological processes, which could be utilized to treat cancer also, using application of electrical pulses. This kind of alternate, novel, and physical therapies are critical to treat brain tumors as 20-40% of all other cancers later develop brain metastasis4 in addition, for children and young adults at ages 20-29, brain cancer is the 2nd leading cause of death.5,6 There are millions of people worldwide suffering from the brain tumor morbidities and mortalities and it necessitates the application of electrical-pulse-mediated therapy known as electrochemotherapy.7‒9 This involves application of high intensity, short duration pulses after injection of a low dose chemo drug to render the plasma membranes permeable to otherwise impermeable drugs. Up to 1000 fold enhancement in drug uptake is observed: hence, electrical-pulse-mediated chemotherapy has a lot of promise as a tool to treat brain tumors effectively and economically. Towards this, we study the electric field distribution in brain tumors under various conditions. It is important to study the electric field distribution, as, if the field intensity is low, the drug could not get into the cells; if the field intensity is high, the cells will be killed, hence for effective drug delivery, the optimum field conditions is necessary and the distribution of this electric field intensity depends on type of electrodes. We studied the electric field distribution of brain tumor using needle and microneedle electrodes, at various depths of the tumor.10

Materials & methods

3D Model

A three-dimentional (3D) model consisting of the human brain, tumor tissue and electrodes was developed to study electric field distribution. Figure 1a shows a brain tumor configuration.11 This is modeled as shown in Figure 1b. Here, the upper electrode (positive) and the lower electrode (ground) are used to apply voltages. The electric field distribution using two needle electrode and microneedle electrode array was studied. The dimensions of model components along with the voltages applied in each case (to obtain a nominal value of 1200V/cm7) are listed in Table 1 and the electrical parameters in Table 2.12

Figure 1 Brain Tumor (a) [11] and the 3D Brain Model (b) created for simulation.

Model Component

Breadth x Height (cm x cm)   

Diameter (cm)

Applied Voltage (V)

Needle electrode
(with 0.25 cm distance between them)    

-

0.1

200

Micro Needle electrode array
(with 1 mm distance between them)

-

0.025

15

Dura matter (Ellipsoid)

-

16 x 20 x 12
(0.5 cm thickness)   

-

Gray matter (Ellipsoid)

-

15 x 19 x 11
(0.5 cm thickness)

-

White matter (Ellipsoid)

-

14 x 18 x 10

-

Tumor (sphere)

-

2

-

Table 1  Dimensions of the geometry used

Material

Relative permittivity, εr

Tumor

260

Dura matter

1.54x103

Gray matter

4.23x104

White matter

2.09x104

Electrodes

1

Table 2 Properties of the materials used inside the Geometry.12

Software

In order to determine the electric field distribution inside the tumor and other layers of brain, a three-dimensional finite element model (Figure 2) was implemented in the COMSOL software package (COMSOL Multiphysics, V.4.3b; Stockholm, Sweden). COMSOL solves a series of equations for solving the given problem. The governing equations that COMSOL uses to compute electric field (E) are given as follows:

D=ρV MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiabgEGirlaads eacqGH9aqpcqaHbpGCcaWGwbaaaa@3CFE@ (1)

D= ε 0 ε r E MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadseacqGH9a qpcqaH1oqzdaWgaaWcbaGaaGimaaqabaGccqaH1oqzdaWgaaWcbaGa amOCaaqabaGccaWGfbaaaa@3F12@ (2)

E=V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=Mj0xXdbba91rFfpec8Eeeu0xXdbba9frFj0=OqFf ea0dXdd9vqaq=JfrVkFHe9pgea0dXdar=Jb9hs0dXdbPYxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaiaadweacqGH9a qpcqGHsislcqGHhis0caWGwbaaaa@3C2C@ (3)

Figure 2 Slices of brain consisting tumor across thickness (Z-Axis).

where, D is electric displacement, r is the electrical charge density, eo is the permittivity of free space, and er is the relative permittivity of the medium. Electrostatics (es) module of AC/DC physics in COMSOL is used to compute electric field while varying geometry, voltage and some other parameters. The electrical boundary conditions along the tissue in contact with the energized electrode were ϕ=V0 (electrode voltage) and ϕ= 0 at the ground electrode. The boundaries where the analyzed domain was not in contact with an electrode were treated as electrically isolative.

Results & discussion

Figure 3 shows a typical distribution of the electric field along the depth using microneedle array electrode. The electric field intensity varies from a lower value at the bottom to a higher value at the top. It is the highest at the tip of the needles. The variation of the electric field distributions for the two electrode configurations, namely, the two needle, and microneedle array electrodes, are shown in Figures 4a-4c, at the 0.25cm, 0cm, and -0.25cm slices respectively (the slices were cut at 0, 0.25, 0.5, 0.75, -0.25, -0.5, and -0.75cm). The electric field intensity is the highest at 121.6V/cm for the two needle electrodes and 7.79V/cm for the microneedle array electrodes, at +0.25 cm slice. At 0cm, the maximum electric field intensities for these electrode configurations are, 1603.1 and 434.49V/cm, respectively. The electric field intensity is highest at the needle configuration due to the sharp edge of the needle. The electric fields at -0.25 cm slice are, 604.86, and 487.64V/cm respectively for the two electrode types. This indicates the variation of the electric field intensity along the tumor thickness, maximum being at the surface where the electrodes are placed (in the center of the tumor, at 0cm) and they vary in magnitudes above and below 0 cm, as shown for 0.25cm and -0.25cm. Figure 5 shows a comparison of average electric field for the two needle electrode at +0.25, 0, and -0.25cm layers along the depth. The variation is highly non-linear from -0.25 to 0 to +0.25 cm depths. It is only 7.6% at +0.25cm, compared to the magnitude at 0 cm, while it is 37.7% at -0.25cm, compared to that at 0 cm. Similar variation is also observed for the microneedle electrode array (Figure 6). In this case it varies to a low 1.8% at the +0.25cm level compared to 0 cm level. However, at -0.2cm level, it is more than at 0 cm level, it is of 112.22% of that of at 0cm level. These variations could be due to the difference in the electrode geometry and their penetration into the tumor thickness at various slice levels. The applied voltage can be varied to obtain desired field intensities.

Figure 3 Electric field distribution along the axes of the brain tumor model.

Figure 4 Electric field distribution at +0.25 (a), 0 (b) and -0.25 cm (c) slices for the two needle and microneedle array electrodes.

Figure 5 Comparison of Maximum Electric Field Magnitudes at various depths inside the tumor (1-at +0.25 cm, 2-at 0 cm, 3-at -0.25 cm levels) for the two needle electrode array.

Figure 6 Comparison of Maximum Electric Field Magnitudes at various depths inside the tumor (1-at +0.25 cm, 2-at 0 cm, 3-at -0.25 cm levels) for microneedle electrode array.

Conclusion

Electric field intensities vary at various depths of the tumors. It also varies based on the electrode geometries, viz. two needle and microneedle array electrode configurations. Maximum field intensity of 1603V/cm is obtained for the two needle electrode at slice at 0 cm. The variation is highly non-uniform, from 0 cm to +0.25 cm and -0.25cm levels. For the given configuration, if higher intensities are needed, two needle electrodes are suitable. However, if lower electric field strengths are needed, the microarray needles are preferred. However, by varying the applied voltage, it is possible to obtain desired electric field intensity at a given depth for a given electrode configuration. This study is useful to select the appropriate intensity at the clinic to do electrochemotherapy for brain tumors.

Acknowledgments

None.

Conflicts of interest

Authors declare there are no conflicts of interest.

References

Creative Commons Attribution License

©2016 Venuturumilli, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.