ESHO benchmarks for computational modeling and optimization in hyperthermia therapy

Helmet Radio Frequency Phased Array Applicators Enhance Thermal Magnetic Resonance of Brain Tumors

Thermal Magnetic Resonance (ThermalMR) integrates Magnetic Resonance Imaging (MRI) diagnostics and targeted radio-frequency (RF) heating in a single theranostic device. The requirements for MRI (magnetic field) and targeted RF heating (electric field) govern the design of ThermalMR applicators. We hypothesize that helmet RF applicators (HPA) improve the efficacy of ThermalMR of brain tumors versus an annular phased RF array (APA). An HPA was designed using eight broadband self-grounded bow-tie (SGBT) antennae plus two SGBTs placed on top of the head. An APA of 10 equally spaced SGBTs was used as a reference. Electromagnetic field (EMF) simulations were performed for a test object (phantom) and a human head model. For a clinical scenario, the head model was modified with a tumor volume obtained from a patient with glioblastoma multiforme. To assess performance, we introduced multi-target evaluation (MTE) to ensure whole-brain slice accessibility. We implemented time multiplexed vector field shaping to optimize RF excitation. Our EMF and temperature simulations demonstrate that the HPA improves performance criteria critical to MRI and enhances targeted RF and temperature focusing versus the APA. Our findings are a foundation for the experimental implementation and application of a HPA en route to ThermalMR of brain tumors.

Adapting Temperature Predictions to MR Imaging in Treatment Position to Improve Simulation-Guided Hyperthermia for Cervical Cancer

Hyperthermia treatment consists of elevating the temperature of the tumor to increase the effectiveness of radiotherapy and chemotherapy. Hyperthermia treatment planning (HTP) is an important tool to optimize treatment quality using pre-treatment temperature predictions. The accuracy of these predictions depends on modeling uncertainties such as tissue properties and positioning. In this study, we evaluated if HTP accuracy improves when the patient is imaged inside the applicator at the start of treatment. Because perfusion is a major uncertainty source, the importance of accurate treatment position and anatomy was evaluated using different perfusion values. Volunteers were scanned using MR imaging without ("planning setup") and with the MR-compatible hyperthermia device ("treatment setup"). Temperature-based quality indicators were used to assess the differences between the standard, apparent and the optimized hyperthermia dose. We conclude that pre-treatment imaging can improve HTP predictions accuracy but also, that tissue perfusion modelling is crucial if temperature-based optimization is applied.

Optimization of laser dosimetry based on patient-specific anatomical models for the ablation of pancreatic ductal adenocarcinoma tumor

Laser-induced thermotherapy has shown promising potential for the treatment of unresectable primary pancreatic ductal adenocarcinoma tumors. Nevertheless, heterogeneous tumor environment and complex thermal interaction phenomena that are established under hyperthermic conditions can lead to under/over estimation of laser thermotherapy efficacy. Using numerical modeling, this paper presents an optimized laser setting for Nd:YAG laser delivered by a bare optical fiber (300 µm in diameter) at 1064 nm working in continuous mode within a power range of 2-10 W. For the thermal analysis, patient-specific 3D models were used, consisting of tumors in different portions of the pancreas. The optimized laser power and time for ablating the tumor completely and producing thermal toxic effects on the possible residual tumor cells beyond the tumor margins were found to be 5 W for 550 s, 7 W for 550 s, and 8 W for 550 s for the pancreatic tail, body, and head tumors, respectively. Based on the results, during the laser irradiation at the optimized doses, thermal injury was not evident either in the 15 mm lateral distances from the optical fiber or in the nearby healthy organs. The present computational-based predictions are also in line with the previous ex vivo and in vivo studies, hence, they can assist in the estimation of the therapeutic outcome of laser ablation for pancreatic neoplasms prior to clinical trials.

Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors

Thermal Magnetic Resonance (ThermalMR) is a theranostic concept that combines diagnostic magnetic resonance imaging (MRI) with targeted thermal therapy in the hyperthermia (HT) range using a radiofrequency (RF) applicator in an integrated system. ThermalMR adds a therapeutic dimension to a diagnostic MRI device. Focused, targeted RF heating of deep-seated brain tumors, accurate non-invasive temperature monitoring and high-resolution MRI are specific requirements of ThermalMR that can be addressed with novel concepts in RF applicator design. This work examines hybrid RF applicator arrays combining loop and self-grounded bow-tie (SGBT) dipole antennas for ThermalMR of brain tumors, at magnetic field strengths of 7.0 T, 9.4 T and 10.5 T. These high-density RF arrays improve the feasible transmission channel count, and provide additional degrees of freedom for RF shimming not afforded by using dipole antennas only, for superior thermal therapy and MRI diagnostics. These improvements are especially relevant for ThermalMR theranostics of deep-seated brain tumors because of the small surface area of the head. ThermalMR RF applicators with the hybrid loop+SGBT dipole design outperformed applicators using dipole-only and loop-only designs, with superior MRI performance and targeted RF heating. Array variants with a horse-shoe configuration covering an arc (270°) around the head avoiding the eyes performed better than designs with 360° coverage, with a 1.3 °C higher temperature rise inside the tumor while sparing healthy tissue. Our EMF and temperature simulations performed on a virtual patient with a clinically realistic intracranial tumor provide a technical foundation for implementation of advanced RF applicators tailored for ThermalMR theranostics of brain tumors.

Field and Temperature Shaping for Microwave Hyperthermia: Recent Treatment Planning Tools to Enhance SAR-Based Procedures

The aim of the article is to provide a summary of the work carried out in the framework of a research project funded by the Italian Ministry of Research. The main goal of the activity was to introduce multiple tools for reliable, affordable, and high-performance microwave hyperthermia for cancer therapy. The proposed methodologies and approaches target microwave diagnostics, accurate in vivo electromagnetic parameters estimation, and improvement in treatment planning using a single device. This article provides an overview of the proposed and tested techniques and shows their complementarity and interconnection. To highlight the approach, we also present a novel combination of specific absorption rate optimization via convex programming with a temperature-based refinement method implemented to mitigate the effect of thermal boundary conditions on the final temperature map. To this purpose, numerical tests were carried out for both simple and anatomically detailed 3D scenarios for the head and neck region. These preliminary results show the potential of the combined technique and improvements in the temperature coverage of the tumor target with respect to the case wherein no refinement is adopted.

Antenna Arrangement in UWB Helmet Brain Applicators for Deep Microwave Hyperthermia

Deep microwave hyperthermia applicators are typically designed as narrow-band conformal antenna arrays with equally spaced elements, arranged in one or more rings. This solution, while adequate for most body regions, might be sub-optimal for brain treatments. The introduction of ultra-wide-band semi-spherical applicators, with elements arranged around the head and not necessarily aligned, has the potential to enhance the selective thermal dose delivery in this challenging anatomical region. However, the additional degrees of freedom in this design make the problem non-trivial. We address this by treating the antenna arrangement as a global SAR-based optimization process aiming at maximizing target coverage and hot-spot suppression in a given patient. To enable the quick evaluation of a certain arrangement, we propose a novel E-field interpolation technique which calculates the field generated by an antenna at any location around the scalp from a limited number of initial simulations. We evaluate the approximation error against full array simulations. We demonstrate the design technique in the optimization of a helmet applicator for the treatment of a medulloblastoma in a paediatric patient. The optimized applicator achieves 0.3 °C higher T90 than a conventional ring applicator with the same number of elements.

Hyperhydration of Cancers: A Characteristic Biophysical Trait Strongly Increasing O2, CO2, Glucose and Lactate Diffusivities, and Improving Thermophysical Properties of Solid Malignancies

Cancers are complex, heterogeneous, dynamic and aggressive diseases exhibiting a series of characteristic biophysical traits which complement the original biological hallmarks of cancers favouring progressive growth, metastasis, and contributing to immune evasion and treatment resistance. One of the prevalent differences between most solid tumors and their corresponding, healthy tissues is a significantly higher water content (hyperhydration) in cancers. As a consequence, cancers have distinctly higher (Fick's) diffusion coefficients D [cm2 s-1] for the respiratory gases O2 and CO2, the key substrate glucose, and for the oncometabolite lactate. In addition, cancers have (a) clearly increased specific heat capacities cp [J g-1 K-1], thus representing high-capacity-tissues upon therapeutic heating induced by electromagnetic irradiation, and (b) higher thermal conductivities k [W m-1 K-1], i.e., increased abilities to conduct heat. Therefore, in diffusion analyses (e.g., when describing critical O2 and glucose supplies or CO2 removal, and the development of hypoxic subvolumes) and for modeling temperature distributions in hyperthermia treatment planning, these specific cancer-related data must be considered in order to reliably reflect oncologic thermo-radiotherapy settings.

Using MRI to measure position and anatomy changes and assess their impact on the accuracy of hyperthermia treatment planning for cervical cancer

PURPOSE

We studied the differences between planning and treatment position, their impact on the accuracy of hyperthermia treatment planning (HTP) predictions, and the relevance of including true treatment anatomy and position in HTP based on magnetic resonance (MR) images.

MATERIALS AND METHODS

All volunteers were scanned with an MR-compatible hyperthermia device, including a filled waterbolus, to replicate the treatment setup. In the planning setup, the volunteers were scanned without the device to reproduce the imaging in the current HTP. First, we used rigid registration to investigate the patient position displacements between the planning and treatment setup. Second, we performed HTP for the planning anatomy at both positions and the treatment mimicking anatomy to study the effects of positioning and anatomy on the quality of the simulated hyperthermia treatment. Treatment quality was evaluated using SAR-based parameters.

RESULTS

We found an average displacement of 2 cm between planning and treatment positions. These displacements caused average absolute differences of ∼12% for TC25 and 10.4%-15.9% in THQ. Furthermore, we found that including the accurate treatment position and anatomy in treatment planning led to an improvement of 2% in TC25 and 4.6%-10.6% in THQ.

CONCLUSIONS

This study showed that precise patient position and anatomy are relevant since these affect the accuracy of HTP predictions. The major part of improved accuracy is related to implementing the correct position of the patient in the applicator. Hence, our study shows a clear incentive to accurately match the patient position in HTP with the actual treatment.

Rotationally Adjustable Hyperthermia Applicators: A Computational Comparative Study of Circular and Linear Array Applicators

Microwave breast hyperthermia (MH) aims to increase the temperature at the tumor location with minimal change in the healthy tissue. To this end, the specific absorption rate (SAR) inside the breast is optimized. The choice of the MH applicator design is important for a superior energy focus on the target. Although hyperthermia treatment planning (HTP) changes for every patient, the MH applicator is required to be effective for different breast models and tumor types. The linear applicator (LA) is one of the previously proposed applicator designs with linearly arranged antennas; however, it suffers from low focusing ability in certain breast regions due to its unsymmetrical geometrical features. In this paper, we propose to radially adjust the LA to obtain alternative excitation schemes without actually changing the applicator. Antipodal Vivaldi antennas were utilized, and the antenna excitations were optimized with particle swarm optimization (PSO). The comparison of the rotated and the fixed linear applicator, between 12-antenna circular and linear applicators, and finally, between a 24-antenna circular applicator are provided. Within the 12 rotation angles and two target locations that were analyzed, the 135° axially rotated linear applicator gave a 35% to 84% higher target-to-breast SAR ratio (TBRS) and a 21% to 28% higher target-to-breast temperature ratio (TBRT) than the fixed linear applicator. For the deep-seated target, the 135° rotated linear applicator had an 80% higher TBRS and a 59% higher TBRT than the 12-antenna circular applicator, while the results were comparable to the 24-antenna circular applicator.

Review of Thermal and Physiological Properties of Human Breast Tissue

Electromagnetic thermal therapies for cancer treatment, such as microwave hyperthermia, aim to heat up a targeted tumour site to temperatures within 40 and 44 °C. Computational simulations used to investigate such heating systems employ the Pennes' bioheat equation to model the heat exchange within the tissue, which accounts for several tissue properties: density, specific heat capacity, thermal conductivity, metabolic heat generation rate, and blood perfusion rate. We present a review of these thermal and physiological properties relevant for hyperthermia treatments of breast including fibroglandular breast, fatty breast, and breast tumours. The data included in this review were obtained from both experimental measurement studies and estimated properties of human breast tissues. The latter were used in computational studies of breast thermal treatments. The measurement methods, where available, are discussed together with the estimations and approximations considered for values where measurements were unavailable. The review concludes that measurement data for the thermal and physiological properties of breast and tumour tissue are limited. Fibroglandular and fatty breast tissue properties are often approximated from those of generic muscle or fat tissue. Tumour tissue properties are mostly obtained from approximating equations or assumed to be the same as those of glandular tissue. We also present a set of reliable data, which can be used for more accurate modelling and simulation studies to better treat breast cancer using thermal therapies.

Accurate Three-Dimensional Thermal Dosimetry and Assessment of Physiologic Response Are Essential for Optimizing Thermoradiotherapy

Numerous randomized trials have revealed that hyperthermia (HT) + radiotherapy or chemotherapy improves local tumor control, progression free and overall survival vs. radiotherapy or chemotherapy alone. Despite these successes, however, some individuals fail combination therapy; not every patient will obtain maximal benefit from HT. There are many potential reasons for failure. In this paper, we focus on how HT influences tumor hypoxia, since hypoxia negatively influences radiotherapy and chemotherapy response as well as immune surveillance. Pre-clinically, it is well established that reoxygenation of tumors in response to HT is related to the time and temperature of exposure. In most pre-clinical studies, reoxygenation occurs only during or shortly after a HT treatment. If this were the case clinically, then it would be challenging to take advantage of HT induced reoxygenation. An important question, therefore, is whether HT induced reoxygenation occurs in the clinic that is of radiobiological significance. In this review, we will discuss the influence of thermal history on reoxygenation in both human and canine cancers treated with thermoradiotherapy. Results of several clinical series show that reoxygenation is observed and persists for 24-48 h after HT. Further, reoxygenation is associated with treatment outcome in thermoradiotherapy trials as assessed by: (1) a doubling of pathologic complete response (pCR) in human soft tissue sarcomas, (2) a 14 mmHg increase in pO2 of locally advanced breast cancers achieving a clinical response vs. a 9 mmHg decrease in pO2 of locally advanced breast cancers that did not respond and (3) a significant correlation between extent of reoxygenation (as assessed by pO2 probes and hypoxia marker drug immunohistochemistry) and duration of local tumor control in canine soft tissue sarcomas. The persistence of reoxygenation out to 24-48 h post HT is distinctly different from most reported rodent studies. In these clinical series, comparison of thermal data with physiologic response shows that within the same tumor, temperatures at the higher end of the temperature distribution likely kill cells, resulting in reduced oxygen consumption rate, while lower temperatures in the same tumor improve perfusion. However, reoxygenation does not occur in all subjects, leading to significant uncertainty about the thermal-physiologic relationship. This uncertainty stems from limited knowledge about the spatiotemporal characteristics of temperature and physiologic response. We conclude with recommendations for future research with emphasis on retrieving co-registered thermal and physiologic data before and after HT in order to begin to unravel complex thermophysiologic interactions that appear to occur with thermoradiotherapy.

Simultaneous ThermoBrachytherapy: Electromagnetic Simulation Methods for Fast and Accurate Adaptive Treatment Planning

The combination of interstitial hyperthermia treatment (IHT) with high dose rate brachytherapy (HDR-BT) can improve clinical outcomes since it highly enhances the efficiency of cell kill, especially when applied simultaneously. Therefore, we have developed the ThermoBrachy applicators. To effectively apply optimal targeted IHT, treatment planning is considered essential. However, treatment planning in IHT is rarely applied as it is regarded as difficult to accurately calculate the deposited energy in the tissue in a short enough time for clinical practice. In this study, we investigated various time-efficient methods for fast computation of the electromagnetic (EM) energy deposition resulting from the ThermoBrachy applicators. Initially, we investigated the use of an electro-quasistatic solver. Next, we extended our investigation to the application of geometric simplifications. Furthermore, we investigated the validity of the superpositioning principle, which can enable adaptive treatment plan optimization without the need for continuous recomputation of the EM field. Finally, we evaluated the accuracy of the methods by comparing them to the golden standard Finite-Difference Time-Domain calculation method using gamma-index analysis. The simplifications considerably reduced the computation time needed, improving from >12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a >99% passing rate with 1%/0.5 mm Dose Difference and Distance-to-Agreement criteria. These results allow the proposed electromagnetic simulation method to be used for fast and accurate adaptive treatment planning.