Simulation of radiofrequency ablation in real human anatomy

Computer simulation-based nanothermal field and tissue damage analysis for cardiac tumor ablation

Radiofrequency ablation is a nominally invasive technique to eradicate cancerous or non-cancerous cells by heating. However, it is still hampered to acquire a successful cell destruction process due to inappropriate RF intensities that will not entirely obliterate tumorous tissues, causing in treatment failure. In this study, we are acquainted with a nanoassisted RF ablation procedure of cardiac tumor to provide better outcomes for long-term survival rate without any recurrences. A three-dimensional thermo-electric energy model is employed to investigate nanothermal field and ablation efficiency into the left atrium tumor. The cell death model is adopted to quantify the degree of tissue injury while injecting the Fe3O4 nanoparticles concentrations up to 20% into the target tissue. The results reveal that when nanothermal field extents as a function of tissue depth (10 mm) from the electrode tip, the increasing thermal rates were approximately 0.54362%, 3.17039%, and 7.27397% for the particle concentration levels of 7%, 10%, and 15% compared with no-particle case. In the 7% Fe3O4 nanoparticles, 100% fractional damage index is achieved after ablation time of 18 s whereas tissue annihilation approach proceeds longer to complete for no-particle case. The outcomes indicate that injecting nanoparticles may lessen ablation time in surgeries and prevent damage to adjacent healthy tissue.

Improving Prediction of the Potential Distribution Induced by Cylindrical Electrodes within a Homogeneous Rectangular Grid during Irreversible Electroporation

Background: Irreversible electroporation (IRE) is an ablation technique based on the application of short, high-voltage pulses between needle electrodes (diameter: ~1.0 × 10−3 m). A Finite Difference-based software simulating IRE treatment generally uses rectangular grids, yielding discretization issues when modeling cylindrical electrodes and potentially affecting the validity of treatment planning simulations. Aim: Develop an Electric-Potential Estimation (EPE) method for accurate prediction of the electric-potential distribution in the vicinity of cylindrical electrodes. Methods: The electric-potential values in the voxels neighboring the cylindrical electrode voxels were corrected based on analytical solutions derived for coaxial/cylindrical electrodes. Simulations at varying grid resolutions were validated using analytical models. Low-resolution heterogeneous simulations at 2.0 × 10−3 m excluding/including EPE were compared with high-resolution results at 0.25 × 10−3 m. Results: EPE significantly reduced maximal errors compared to analytical results for the electric-potential distributions (26.6–71.8%→0.4%) and for the electrical resistance (30%→1–6%) at 3.0 × 10−3 m voxel-size. EPE significantly improved the mean-deviation (43.1–52.8%→13.0–24.3%) and the calculation-time gain (>15,000×) of low-resolution compared to high-resolution heterogeneous simulations. Conclusions: EPE can accurately predict the potential distribution of neighboring cylindrical electrodes, regardless of size, position, and orientation in a rectangular grid. The simulation time of treatment planning can therefore be shortened by using large voxel-sized models without affecting accuracy of the electric-field distribution, enabling real-time clinical IRE treatment planning.

Pre-operative Assessment of Ablation Margins for Variable Blood Perfusion Metrics in a Magnetic Resonance Imaging Based Complex Breast Tumour Anatomy: Simulation Paradigms in Thermal Therapies

BACKGROUND AND OBJECTIVES

Image-guided medical interventions facilitates precise visualization at treatment site. The conformal prediction for sparing healthy tissue fringes precisely in the vicinity of irregular tumour anatomy remains clinically challenging. Pre-clinical image-based computational modelling is imperative as it helps in enhancement of treatment quality, augmenting clinical-decision making, while planning, targeting, controlling, monitoring and assessing treatment response with an effective risk assessment before the onset of treatment in clinical settings. In this study, the influence of heat deposition rate (SAR), exposure duration, and variable blood perfusion metrics for a patient-specific breast tumour is quantified considering the tumour margins thereby suggesting need of geometrically accurate models.

METHODS

A three-dimensional realistic model mimicking dimensions of a female breast, comprising ~1.7 cm irregular tumour, was generated from patient specific two-dimensional DICOM format MRI images through image segmentation tools MIMICS 19.0® and 3-Matic 11.0® which is finally exported to COMSOL Multiphysics 5.2® as a volumetric mesh for finite element analysis. The Pennes bioheat transfer model and Arrhenius thermal damage model of cell-death are integrated to simulate a coupled biophysics problem. A comparative blood perfusion analysis is done to evaluate the response of tumour during heating considering thermal damage extent, including the tumour margins while sparing critical adjoining healthy tissues.

RESULTS

The evaluated thermal damage zones for 1 mm, 2 mm and 3 mm fringe heating region (beyond tumour boundary) reveals 0.09%, 0.21% and 0.34% thermal damage to the healthy tissue (which is <1%) and thus successful necrosis of the tumour. The iterative computational experiments suggests treatment margins < 5 mm are sufficient enough as heating beyond 3 mm fringe layer leads to higher damage surrounding the tumour approximately 1.5 times the tumour volume. Further, the heat-dosage requirements are 22% more for highly perfused tumour as compared to moderately perfused tumour with an approximate double time to ablate the whole tumour volume.

CONCLUSIONS

Depending on the blood perfusion characteristics of a tumour, it is a trade-off between heat-dosage (SAR) and exposure/treatment duration to get desired thermal damage including the irregular tumour boundaries while taking into account, the margin of healthy tissue. The suggested patient-specific integrated multiphysics-model based on MRI-Images may be implemented for pre-treatment planning based on the tumour blood perfusion to evaluate the thermal ablation zone dimensions clinically and thereby avoiding the damage of off-target tissues. Thus, risks involving underestimation or overestimation of thermal coagulation zones may be minimised while preserving the surrounding normal breast parenchyma.

Robot-assisted ultrasound navigation platform for 3D HIFU treatment planning: Initial evaluation for conformal interstitial ablation

Interstitial Ultrasound-guided High Intensity Focused Ultrasound (USgHIFU) therapy has the potential to deliver ablative treatments which conform to the target tumor. In this study, a robot-assisted US-navigation platform has been developed for 3D US guidance and planning of conformal HIFU ablations. The platform was used to evaluate a conformal therapeutic strategy associated with an interstitial dual-mode USgHIFU catheter prototype (64 elements linear-array, measured central frequency f = 6.5 MHz), developed for the treatment of HepatoCellular Carcinoma (HCC). The platform included a 3D navigation environment communicating in real-time with an open research dual-mode US scanner/HIFU generator and a robotic arm, on which the USgHIFU catheter was mounted. 3D US-navigation was evaluated in vitro for guiding and planning conformal HIFU ablations using a tumor-mimic model in porcine liver. Tumor-mimic volumes were then used as targets for evaluating conformal HIFU treatment planning in simulation. Height tumor-mimics (ovoid- or disc-shaped, sizes: 3-29 cm³) were created and visualized in liver using interstitial 2D US imaging. Robot-assisted spatial manipulation of these images and real-time 3D navigation allowed reconstructions of 3D B-mode US images for accurate tumor-mimic volume estimation (relative error: 4 ± 5%). Sectorial and full-revolution HIFU scanning (angular sectors: 88-360°) could both result in conformal ablations of the tumor volumes, as soon as their radii remained ≤ 24 mm. The presented US navigation-guided HIFU procedure demonstrated advantages for developing conformal interstitial therapies in standard operative rooms. Moreover, the modularity of the developed platform makes it potentially useful for developing other HIFU approaches.

Thermal ablation of biological tissues in disease treatment: A review of computational models and future directions

Percutaneous thermal ablation has proven to be an effective modality for treating both benign and malignant tumours in various tissues. Among these modalities, radiofrequency ablation (RFA) is the most promising and widely adopted approach that has been extensively studied in the past decades. Microwave ablation (MWA) is a newly emerging modality that is gaining rapid momentum due to its capability of inducing rapid heating and attaining larger ablation volumes, and its lesser susceptibility to the heat sink effects as compared to RFA. Although the goal of both these therapies is to attain cell death in the target tissue by virtue of heating above 50°C, their underlying mechanism of action and principles greatly differs. Computational modelling is a powerful tool for studying the effect of electromagnetic interactions within the biological tissues and predicting the treatment outcomes during thermal ablative therapies. Such a priori estimation can assist the clinical practitioners during treatment planning with the goal of attaining successful tumour destruction and preservation of the surrounding healthy tissue and critical structures. This review provides current state-of-the-art developments and associated challenges in the computational modelling of thermal ablative techniques, viz., RFA and MWA, as well as touch upon several promising avenues in the modelling of laser ablation, nanoparticles assisted magnetic hyperthermia and non-invasive RFA. The application of RFA in pain relief has been extensively reviewed from modelling point of view. Additionally, future directions have also been provided to improve these models for their successful translation and integration into the hospital work flow.

Shape-shifting thermal coagulation zone during saline-infused radiofrequency ablation: A computational study on the effects of different infusion location

BACKGROUND AND OBJECTIVE

The majority of the studies on radiofrequency ablation (RFA) have focused on enlarging the size of the coagulation zone. An aspect that is crucial but often overlooked is the shape of the coagulation zone. The shape is crucial because the majority of tumours are irregularly-shaped. In this paper, the ability to manipulate the shape of the coagulation zone following saline-infused RFA by altering the location of saline infusion is explored.

METHODS

A 3D model of the liver tissue was developed. Saline infusion was described using the dual porosity model, while RFA was described using the electrostatic and bioheat transfer equations. Three infusion locations were investigated, namely at the proximal end, the middle and the distal end of the electrode. Investigations were carried out numerically using the finite element method.

RESULTS

Results indicated that greater thermal coagulation was found in the region of tissue occupied by the saline bolus. Infusion at the middle of the electrode led to the largest coagulation volume followed by infusion at the proximal and distal ends. It was also found that the ability to delay roll-off, as commonly associated with saline-infused RFA, was true only for the case when infusion is carried out at the middle. When infused at the proximal and distal ends, the occurrence of roll-off was advanced. This may be due to the rapid and more intense heating experienced by the tissue when infusion is carried out at the electrode ends where Joule heating is dominant.

CONCLUSION

Altering the location of saline infusion can influence the shape of the coagulation zone following saline-infused RFA. The ability to 'shift' the coagulation zone to a desired location opens up great opportunities for the development of more precise saline-infused RFA treatment that targets specific regions within the tissue.

Possible involvement of HSP70 in pancreatic cancer cell proliferation after heat exposure and impact on RFA postoperative patient prognosis

Purpose

As an alleviative treatment measured in patients with unresectable advanced pancreatic cancer, radiofrequency ablation (RFA) needed more clinical data to prove its advantages and to explore limitations in its utilization. This study was determined to observe the efficacy of RFA, and to explore its impact on perioperative periphery carcinoma as well as the normal pancreatic tissues.

Methods

Clinical data of 32 patients with pancreatic cancer accepted RFA surgery were collected. Followed up patients' pain degree and the changes in serum tumor markers CA19-9 and CA 242 before and after surgery. Ex vivo, gave human pancreatic cancer cell line PANC-1 heat treatment to simulate the heat exposure condition periphery carcinoma was experienced during RFA surgery, and to observe the proliferation rate and HSP70 expression change compared with control group.

Results

Of the 32 patients, 1 died of upper gastrointestinal hemorrhage, and 29 survived for more than 5 months, 2 of which for more than 16 months. The average CA19-9 and CA 242 levels of the patients were significantly decreased in 3 months after surgery (t = 9.873, 5.978, P < 0.001). During in vitro experiments, the proliferation rate of PANC-1 cells after heating was significantly increased, accompanied with the increased HSP70 expression. The addition of HSP70 inhibitor can inhibit the rise of proliferation after heat therapy.

Conclusion

Utilizing RFA treat patients with unresectable advanced pancreatic cancer, could effectively relieve the pain, decline jaundice, and deduce tumor marker levels significantly. However, it failed to extend the long-term survival rate of the patients significantly. This study found that a higher proliferative rate accompanied with a higher HSP70 expression level were observed on in vitro cultured pancreatic carcinoma cells after heat treatment, which could be altered by HSP70 inhibitor. And these findings indicated that the heat exposure might impact periphery carcinoma during RFA surgery and HSP70 might play an important role in patients' prognosis.

Coupled thermo-electro-mechanical models for thermal ablation of biological tissues and heat relaxation time effects

Thermal ablation is a widely applied electrosurgical process in medical treatment of soft biological tissues. Numerical modeling and simulations play an important role in prediction of temperature distribution and damage volume during the treatment planning stage of associated therapies. In this contribution we report a coupled thermo-electro-mechanical model, accounting for heat relaxation time, for more accurate and precise prediction of the temperature distribution, tissue deformation and damage volume during the thermal ablation of biological tissues. Finite element solutions are obtained for most widely used percutaneous thermal ablative techniques, viz., radiofrequency ablation (RFA) and microwave ablation (MWA). Importantly, both tissue expansion and shrinkage have been considered for modeling the tissue deformation in the coupled model of high temperature thermal ablation. The coupled model takes into account the non-Fourier effects, considering both single-phase-lag (SPL) and dual-phase-lag (DPL) models of bio-heat transfer. The temperature-dependent electrical and thermal parameters, damage-dependent blood perfusion rate and phase change effect accounting for tissue vaporization have been accounted for obtaining more clinically relevant model. The proposed model predictions are found to be in good agreement against the temperature distribution and damage volume reported by previous experimental studies. The numerical simulation results revealed that the non-Fourier effects cause a decrease in the predicted temperature distribution, tissue deformation and damage volume during the high temperature thermal ablative procedures. Furthermore, the effects of different magnitudes of phase lags of the heat flux and temperature gradient on the predicted treatment outcomes of the considered thermal ablative modalities are also quantified and discussed in detail.

Design of a Novel Electrode of Radiofrequency Ablation for Large Tumors: In Vitro Validation and Evaluation

In the present study, a monopolar expandable electrode (MEE) in radiofrequency ablation (RFA) proposed in our previous study was validated and evaluated using the in vitro experiment and computer simulation. Two commercial RF electrodes (conventional electrode, CE and umbrella electrode, UE) was used to compare the ablation results with MEE using the in vitro egg white model (experiment and computer simulation) and in vivo liver tumor model (computer simulation) to verify the efficacy of MEE in the large tumor ablation. The sharp increase in impedance during RFA procedures was taken as the termination of RFA protocols. The volume and sphericity of ablation zone generated by MEE, CE, and UE in the in vitro egg white experiment were 75.3 1.6 cm3, 2.7 0.4 cm3, 12.4 1.8 cm3 (P <0.001), and 88.1 0.9%, 12.9 1.3%, 62.0 3.0% (P <0.001), respectively. Correspondingly, a similar result was obtained in the egg white simulation. In the liver tumor simulation, the volume and sphpericity of ablation zone generated by MEE, CE, and UE were 35.4 cm3 and 86.8%, 3.7 cm3 and 17.7%, and 12.7 cm3 and 59.6%, respectively. In summary, MEE has the potential to achieve complete ablation in the treatment of large tumors (>3 cm in diameter) compared with CE and UE due to the larger electrode-tissue interface and more round shape of hooks.

A Pilot Study Into the Use of FDG-mNP as an Alternative Approach in Neuroblastoma Cell Hyperthermia

Herein, we present a pilot study concerning the use of fluorodeoxy glucose conjugated magnetite nanoparticles (FDG-mNP) as a potential agent in magnetic nanoparticle mediated neuroblastoma cancer cell hyperthermia. This approach makes use of the 'Warburg effect', utilizing the fact that cancer cells have a higher metabolic rate than normal cells. FDG-mNP were synthesized, then applied to the SH-SY5Y neuroblastoma cancer cell line and exposed to an ac magnetic field. 3D Calorimetry was performed on the FDG-mNP compound. Simulations were performed using SEMCAD X software using Thelonious, (an anatomically correct male child model) in order to understand more about the end requirements with respect to cancer cell destruction. We investigated FDG-mNP mediated neuroblastoma cytotoxicity in conjunction with ac magnetic field exposure. Results are presented for 3D FDG-mNP SAR _mnp (10.86 ± 0.99 W/g of particles) using a therapeutic dose of 0.83 mg/ mL. Human model simulations suggest that 43 W/kg SAR _Theo would be required to obtain 42 °C within the centre of a liver tumor (Tumor size, bounding box x = 64, y = 61, z = 65 [mm]), and that the temperature distribution is inhomogeneous within the tumor. Our study suggests that this approach could potentially be used to increase the temperature within cells that would result in cancer cell death due to hyperthermia. Further development of this research will also involve using whole tumors removed from living organisms in conjunction with magnetic resonance imaging and positron emission tomography.

Computational simulation of the predicted dosimetric impact of adjuvant yttrium-90 PET/CT-guided percutaneous ablation following radioembolization

BACKGROUND

90Y PET/CT post-radioembolization imaging has demonstrated that the distribution of 90Y in a tumor can be non-uniform. Using computational modeling, we predicted the dosimetric impact of post-treatment 90Y PET/CT-guided percutaneous ablation of the portions of a tumor receiving the lowest absorbed dose. A cohort of fourteen patients with non-resectable liver cancer previously treated using 90Y radioembolization were included in this retrospective study. Each patient exhibited potentially under-treated areas of tumor following treatment based on quantitative 90Y PET/CT. 90Y PET/CT was used to guide electrode placement for simulated adjuvant radiofrequency ablation in areas of tumor receiving the lowest dose. The finite element method was used to solve Penne's bioheat transport equation, coupled with the Arrhenius thermal cell-death model to determine 3D thermal ablation zones. Tumor and unablated tumor absorbed-dose metrics (average dose, D50, D70, D90, V100) following ablation were compared, where D70 is the minimum dose to 70% of tumor and V100 is the fractional tumor volume receiving more than 100 Gy.

RESULTS

Compared to radioembolization alone, 90Y radioembolization with adjuvant ablation was associated with predicted increases in all tumor dose metrics evaluated. The mean average absorbed dose increased by 11.2 ± 6.9 Gy. Increases in D50, D70, and D90 were 11.0 ± 6.9 Gy, 13.3 ± 10.9 Gy, and 11.8 ± 10.8 Gy, respectively. The mean increase in V100 was 7.2 ± 4.2%. All changes were statistically significant (P < 0.01). A negative correlation between pre-ablation tumor volume and D50, average dose, and V100 was identified (ρ < - 0.5, P < 0.05) suggesting that adjuvant radiofrequency ablation may be less beneficial to patients with large tumor burdens.

CONCLUSIONS

This study has demonstrated that adjuvant 90Y PET/CT-guided radiofrequency ablation may improve tumor absorbed-dose metrics. These data may justify a prospective clinical trial to further evaluate this hybrid approach.

A review of radiofrequency ablation: Large target tissue necrosis and mathematical modelling

Radiofrequency ablation (RFA) is an effective clinical method for tumour ablation with minimum intrusiveness. However, the use of RFA is mostly restricted to small tumours, especially those <3cm in diameter. This paper discusses the state-of-the-art of RFA, drawn from experimental and clinical results, for large tumours (i.e. ⩾3cm in diameter). In particular, the paper analyses clinical results related to target tissue necrosis (TTN) and mathematical modelling of the RFA procedure to understand the mechanism whereby the TTN is limited to under 3cm with RFA. This paper also discusses a strategy of controlling of the temperature of target tissue in the RFA procedure with the state-of-art device, which has the potential to increase the size of TTN. This paper ends with a discussion of some future ideas to solve the so-called 3-cm problem with RFA.

Intrauterine adhesion after transvaginal ultrasound-guided radiofrequency myolysis

Transvaginal ultrasound-guided radiofrequency myolysis offers an alternative to surgery for symptomatic uterine myomas, with encouraging efficacy and safety, suggesting low complication rates. This case study describes the first reported intrauterine adhesion after transvaginal ultrasound-guided radiofrequency myolysis. Women who desire further pregnancy should be warned about the particular risk of intrauterine adhesion of radiofrequency myolysis.