Clinical implementation of hyperthermia treatment planning guided steering: A cross over trial to assess its current contribution to treatment quality

Robust, planning-based targeted locoregional tumour heating in small animals

Objective.To improve hyperthermia in clinical practice, pre-clinical hyperthermia research is essential to investigate hyperthermia effects and assess novel treatment strategies. Translating pre-clinical hyperthermia findings into clinically viable protocols requires laboratory animal treatment techniques similar to clinical hyperthermia techniques. The ALBA micro8 electromagnetic heating system (Med-logix SRL, Rome, Italy) has recently been developed to provide the targeted locoregional tumour heating currently lacking for pre-clinical research. This study evaluates the heat focusing properties of this device and its ability to induce robust locoregional tumour heating under realistic physiological conditions using simulations.Approach.Simulations were performed using the Plan2Heat treatment planning package (Amsterdam UMC, the Netherlands). First, the specific absorption rate (SAR) focus was characterised using a homogeneous phantom. Hereafter, a digital mouse model was used for the characterisation of heating robustness in a mouse. Device settings were optimised for treatment of a pancreas tumour and tested for varying circumstances. The impact of uncertainties in tissue property and perfusion values was evaluated using polynomial chaos expansion. Treatment quality and robustness were evaluated based on SAR and temperature distributions.Main results.The SAR distributions within the phantom are well-focused and can be adjusted to target any specific location. The focus size (full-width half-maximum) is a spheroid with diameters 9 mm (radially) and 20 mm (axially). The mouse model simulations show strong robustness against respiratory motion and intestine and stomach filling (∆T90≤0.14°C).Mouse positioning errors in the cranial-caudal direction lead to∆T90≤0.23°C. Uncertainties in tissue property and perfusion values were found to impact the treatment plan up to 0.56 °C (SD), with a variation onT90of 0.32 °C (1 SD).Significance.Our work shows that the pre-clinical phased-array system can provide adequate and robust locoregional heating of deep-seated target regions in mice. Using our software, robust treatment plans can be generated for pre-clinical hyperthermia research.

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.

Quantification of tissue property and perfusion uncertainties in hyperthermia treatment planning: Multianalysis using polynomial chaos expansion

INTRODUCTION

Hyperthermia treatment planning (HTP) tools can guide treatment delivery, particularly with locoregional radiative phased array systems. Uncertainties in tissue and perfusion property values presently lead to quantitative inaccuracy of HTP, leading to sub-optimal treatment. Assessment of these uncertainties would allow for better judgement of the reliability of treatment plans and improve their value for treatment guidance. However, systematically investigating the impact of all uncertainties on treatment plans is a complex, high-dimensional problem and too computationally expensive for traditional Monte Carlo approaches. This study aims to systematically quantify the treatment-plan impact of tissue property uncertainties by investigating their individual contribution to, and combined impact on predicted temperature distributions.

METHODS

A novel Polynomial Chaos Expansion (PCE)-based HTP uncertainty quantification was developed and applied for locoregional hyperthermia of modelled tumours in the pancreatic head, prostate, rectum, and cervix. Patient models were based on the Duke and Ella digital human models. Using Plan2Heat, treatment plans were created to optimise tumour temperature (represented by T90) for treatment using the Alba4D system. For all 25-34 modelled tissues, the impact of tissue property uncertainties was analysed individually i.e., electrical and thermal conductivity, permittivity, density, specific heat capacity and perfusion. Next, combined analyses were performed on the top 30 uncertainties with the largest impact.

RESULTS

Uncertainties in thermal conductivity and heat capacity were found to have negligible impact on the predicted temperature ( < 1 × 10-10 °C), density and permittivity uncertainties had a small impact (< 0.3 °C). Uncertainties in electrical conductivity and perfusion can lead to large variations in predicted temperature. However, variations in muscle properties result in the largest impact at locations that could limit treatment quality, with a standard deviation up to almost 6 °C (pancreas) and 3.5 °C (prostate) for perfusion and electrical conductivity, respectively. The combined influence of all significant uncertainties leads to large variations with a standard deviation up to 9.0, 3.6, 3.7 and 4.1 °C for the pancreatic, prostate, rectal and cervical cases, respectively.

CONCLUSION

Uncertainties in tissue and perfusion property values can have a large impact on predicted temperatures from hyperthermia treatment planning. PCE-based analysis helps to identify all major uncertainties, their impact and judge the reliability of treatment plans.

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.

POD-Kalman filtering for improving noninvasive 3D temperature monitoring in MR-guided hyperthermia

Background

During resonance frequency (RF) hyperthermia treatment, the temperature of the tumor tissue is elevated to the range of 39–44°C. Accurate temperature monitoring is essential to guide treatments and ensure precise heat delivery and treatment quality. Magnetic resonance (MR) thermometry is currently the only clinical method to measure temperature noninvasively in a volume during treatment. However, several studies have shown that this approach is not always sufficiently accurate for thermal dosimetry in areas with motion, such as the pelvic region. Model‐based temperature estimation is a promising approach to correct and supplement 3D online temperature estimation in regions where MR thermometry is unreliable or cannot be measured. However, complete 3D temperature modeling of the pelvic region is too complex for online usage.

Purpose

This study aimed to evaluate the use of proper orthogonal decomposition (POD) model reduction combined with Kalman filtering to improve temperature estimation using MR thermometry. Furthermore, we assessed the benefit of this method using data from hyperthermia treatment where there were limited and unreliable MR thermometry measurements.

Methods

The performance of POD–Kalman filtering was evaluated in several heating experiments and for data from patients treated for locally advanced cervical cancer. For each method, we evaluated the mean absolute error (MAE) concerning the temperature measurements acquired by the thermal probes, and we assessed the reproducibility and consistency using the standard deviation of error (SDE). Furthermore, three patient groups were defined according to susceptibility artifacts caused by the level of intestinal gas motion to assess if the POD–Kalman filtering could compensate for missing and unreliable MR thermometry measurements.

Results

First, we showed that this method is beneficial and reproducible in phantom experiments. Second, we demonstrated that the combined method improved the match between temperature prediction and temperature acquired by intraluminal thermometry for patients treated for locally advanced cervical cancer. Considering all patients, the POD–Kalman filter improved MAE by 43% (filtered MR thermometry = 1.29°C, POD–Kalman filtered temperature = 0.74°C). Moreover, the SDE was improved by 47% (filtered MR thermometry = 1.16°C, POD–Kalman filtered temperature = 0.61°C). Specifically, the POD–Kalman filter reduced the MAE by approximately 60% in patients whose MR thermometry was unreliable because of the great amount of susceptibilities caused by the high level of intestinal gas motion.

Conclusions

We showed that the POD–Kalman filter significantly improved the accuracy of temperature monitoring compared to MR thermometry in heating experiments and hyperthermia treatments. The results demonstrated that POD–Kalman filtering can improve thermal dosimetry during RF hyperthermia treatment, especially when MR thermometry is inaccurate.

Adapt2Heat: treatment planning-assisted locoregional hyperthermia by on-line visualization, optimization and re-optimization of SAR and temperature distributions

BACKGROUND

Hyperthermia treatment planning is increasingly used in clinical applications and recommended in quality assurance guidelines. Assistance in phase-amplitude steering during treatment requires dedicated software for on-line visualization of SAR/temperature distributions and fast re-optimization in response to hot spots. As such software tools are not yet commercially available, we developed Adapt2Heat for on-line adaptive hyperthermia treatment planning and illustrate possible application by different relevant real patient examples.

METHODS

Adapt2Heat was developed as a separate module of the treatment planning software Plan2Heat. Adapt2Heat runs on a Linux operating system and was developed in C++, using the open source Qt, Qwt and VTK libraries. A graphical user interface allows interactive and flexible on-line use of hyperthermia treatment planning. Predicted SAR/temperature distributions and statistics for selected phase-amplitude settings can be visualized instantly and settings can be re-optimized manually or automatically in response to hot spots.

RESULTS

Pretreatment planning E-Field, SAR and temperature calculations are performed with Plan2Heat and imported in Adapt2Heat. Examples show that Adapt2Heat can be helpful in assisting with phase-amplitude steering, e.g., by suppressing indicated hot spots. The effects of phase-amplitude adjustments on the tumor and potential hot spot locations are comprehensively visualized, allowing intuitive and flexible assistance by treatment planning during locoregional hyperthermia treatments.

CONCLUSION

Adapt2Heat provides an intuitive and flexible treatment planning tool for on-line treatment planning-assisted hyperthermia. Extensive features for visualization and (re-)optimization during treatment allow practical use in many locoregional hyperthermia applications. This type of tools are indispensable for enhancing the quality of hyperthermia treatment delivery.

Fast Adaptive Temperature-Based Re-Optimization Strategies for On-Line Hot Spot Suppression during Locoregional Hyperthermia

Simple Summary

When treatment limiting hot spots occur during locoregional hyperthermia (i.e., heating tumors to 40–44 °C for ~1 h), system settings are adjusted based on experience. In this study, we developed and evaluated treatment planning with temperature-based re-optimization and compared the predicted effectiveness to clinically applied protocol/experience-based steering. Re-optimization times were typically ~10 s; sufficiently fast for on-line use. Effective hot spot suppression was predicted, while maintaining adequate tumor heating. Inducing new hot spots was avoided. Temperature-based re-optimization to suppress treatment limiting hot spots seemed feasible to match the effectiveness of long-term clinical experience and will be further evaluated in a clinical setting. When numerical algorithms are proven to match long-term experience, the overall treatment quality within hyperthermia centers can significantly improve. Implementing these strategies would then imply that treatments become less dependent on the experience of the center/operator.

Abstract

Background: Experience-based adjustments in phase-amplitude settings are applied to suppress treatment limiting hot spots that occur during locoregional hyperthermia for pelvic tumors. Treatment planning could help to further optimize treatments. The aim of this research was to develop temperature-based re-optimization strategies and compare the predicted effectiveness with clinically applied protocol/experience-based steering. Methods: This study evaluated 22 hot spot suppressions in 16 cervical cancer patients (mean age 67 ± 13 year). As a first step, all potential hot spot locations were represented by a spherical region, with a user-specified diameter. For fast and robust calculations, the hot spot temperature was represented by a user-specified percentage of the voxels with the largest heating potential (HPP). Re-optimization maximized tumor T90, with constraints to suppress the hot spot and avoid any significant increase in other regions. Potential hot spot region diameter and HPP were varied and objective functions with and without penalty terms to prevent and minimize temperature increase at other potential hot spot locations were evaluated. Predicted effectiveness was compared with clinically applied steering results. Results: All strategies showed effective hot spot suppression, without affecting tumor temperatures, similar to clinical steering. To avoid the risk of inducing new hot spots, HPP should not exceed 10%. Adding a penalty term to the objective function to minimize the temperature increase at other potential hot spot locations was most effective. Re-optimization times were typically ~10 s. Conclusion: Fast on-line re-optimization to suppress treatment limiting hot spots seems feasible to match effectiveness of ~30 years clinical experience and will be further evaluated in a clinical setting.

ESHO benchmarks for computational modeling and optimization in hyperthermia therapy

BACKGROUND

The success of cancer hyperthermia (HT) treatments is strongly dependent on the temperatures achieved in the tumor and healthy tissues as it correlates with treatment efficacy and safety, respectively. Hyperthermia treatment planning (HTP) simulations have become pivotal for treatment optimization due to the possibility for pretreatment planning, optimization and decision making, as well as real-time treatment guidance.

MATERIALS AND METHODS

The same computational methods deployed in HTP are also used for in silico studies. These are of great relevance for the development of new HT devices and treatment approaches. To aid this work, 3 D patient models have been recently developed and made available for the HT community. Unfortunately, there is no consensus regarding tissue properties, simulation settings, and benchmark applicators, which significantly influence the clinical relevance of computational outcomes.

RESULTS AND DISCUSSION

Herein, we propose a comprehensive set of applicator benchmarks, efficacy and safety optimization algorithms, simulation settings and clinical parameters, to establish benchmarks for method comparison and code verification, to provide guidance, and in view of the 2021 ESHO Grand Challenge (Details on the ESHO grand challenge on HTP will be provided at https://www.esho.info/).

CONCLUSION

We aim to establish guidelines to promote standardization within the hyperthermia community such that novel approaches can quickly prove their benefit as quickly as possible in clinically relevant simulation scenarios. This paper is primarily focused on radiofrequency and microwave hyperthermia but, since 3 D simulation studies on heating with ultrasound are now a reality, guidance as well as a benchmark for ultrasound-based hyperthermia are also included.

Treatment planning facilitates clinical decision making for hyperthermia treatments

Background: Treatment quality is important in clinical hyperthermia. Guideline-based treatment protocols are used to determine system settings and treatment strategies to ensure effective tumor heating and prevent unwanted treatment-limiting normal tissue hot spots. Realizing both these goals can prove challenging using generic guideline-based and operator-dependent treatment strategies. Hyperthermia treatment planning (HTP) can be very useful to support treatment strategies. Although HTP is increasingly integrated into the standard clinical workflow, active clinical application is still limited to a small number of hyperthermia centers and should be further stimulated.Purpose: This paper aims to serve as a practical guide, demonstrating how HTP can be applied in clinical decision making for both superficial and locoregional hyperthermia treatments.HTP in clinical decision making: Seven problems that occur in daily clinical practice are described and we show how HTP can enhance insight to formulate an adequate treatment strategy. Examples use representative commercially available hyperthermia devices and cover all stages during the clinical workflow. Problems include selecting adequate phase settings, heating ability analysis, hot spot suppression, applicator selection, evaluation of target coverage and heating depth, and predicting possible thermal toxicity in case of an implant. Since we aim to promote a general use of HTP in daily practice, basic simulation strategies are used in these problems, avoiding a need for the application of dedicated advanced optimization routines that are not generally available.Conclusion: Even fairly basic HTP can facilitate clinical decision making, providing a meaningful and clinically relevant contribution to maintaining and improving treatment quality.

Advanced patient-specific hyperthermia treatment planning

Hyperthermia treatment planning (HTP) is valuable to optimize tumor heating during thermal therapy delivery. Yet, clinical hyperthermia treatment plans lack quantitative accuracy due to uncertainties in tissue properties and modeling, and report tumor absorbed power and temperature distributions which cannot be linked directly to treatment outcome. Over the last decade, considerable progress has been made to address these inaccuracies and therefore improve the reliability of hyperthermia treatment planning. Patient-specific electrical tissue conductivity derived from MR measurements has been introduced to accurately model the power deposition in the patient. Thermodynamic fluid modeling has been developed to account for the convective heat transport in fluids such as urine in the bladder. Moreover, discrete vasculature trees have been included in thermal models to account for the impact of thermally significant large blood vessels. Computationally efficient optimization strategies based on SAR and temperature distributions have been established to calculate the phase-amplitude settings that provide the best tumor thermal dose while avoiding hot spots in normal tissue. Finally, biological modeling has been developed to quantify the hyperthermic radiosensitization effect in terms of equivalent radiation dose of the combined radiotherapy and hyperthermia treatment. In this paper, we review the present status of these developments and illustrate the most relevant advanced elements within a single treatment planning example of a cervical cancer patient. The resulting advanced HTP workflow paves the way for a clinically feasible and more reliable patient-specific hyperthermia treatment planning.

Heating technology for malignant tumors: a review

The therapeutic application of heat is very effective in cancer treatment. Both hyperthermia, i.e., heating to 39-45 °C to induce sensitization to radiotherapy and chemotherapy, and thermal ablation, where temperatures beyond 50 °C destroy tumor cells directly are frequently applied in the clinic. Achievement of an effective treatment requires high quality heating equipment, precise thermal dosimetry, and adequate quality assurance. Several types of devices, antennas and heating or power delivery systems have been proposed and developed in recent decades. These vary considerably in technique, heating depth, ability to focus, and in the size of the heating focus. Clinically used heating techniques involve electromagnetic and ultrasonic heating, hyperthermic perfusion and conductive heating. Depending on clinical objectives and available technology, thermal therapies can be subdivided into three broad categories: local, locoregional, or whole body heating. Clinically used local heating techniques include interstitial hyperthermia and ablation, high intensity focused ultrasound (HIFU), scanned focused ultrasound (SFUS), electroporation, nanoparticle heating, intraluminal heating and superficial heating. Locoregional heating techniques include phased array systems, capacitive systems and isolated perfusion. Whole body techniques focus on prevention of heat loss supplemented with energy deposition in the body, e.g., by infrared radiation. This review presents an overview of clinical hyperthermia and ablation devices used for local, locoregional, and whole body therapy. Proven and experimental clinical applications of thermal ablation and hyperthermia are listed. Methods for temperature measurement and the role of treatment planning to control treatments are discussed briefly, as well as future perspectives for heating technology for the treatment of tumors.

Effect of gastrointestinal gas on the temperature distribution in pancreatic cancer hyperthermia treatment planning

PURPOSE

In pancreatic cancer treatment, hyperthermia can be added to increase efficacy of chemo- and/or radiotherapy. Gas in stomach, intestines and colon is often in close proximity to the target volume. We investigated the impact of variations in gastrointestinal gas (GG) on temperature distributions during simulated hyperthermia treatment (HT).

METHODS

We used sets of one CT and eight cone-beam CT (CBCT) scans obtained prior to/during fractionated image-guided radiotherapy in four pancreatic cancer patients. In Plan2Heat, we simulated locoregional heating by an ALBA-4D phased array radiofrequency system and calculated temperature distributions for (i) the segmented CT (sCT), (ii) sCT with GG replaced by muscle (sCT₀), (iii) sCT₀ with eight different GG distributions as visible on CBCT inserted (sCT_CBCT). We calculated cumulative temperature-volume histograms for the clinical target volume (CTV) for all ten temperature distributions for each patient and investigated the relationship between GG volume and change in ΔT₅₀ (temperature increase at 50% of CTV volume). We determined location and volume of normal tissue receiving a high thermal dose.

RESULTS

GG volume on CBCT varied greatly (9-991 cm³). ΔT₅₀ increased for increasing GG volume; maximum ΔT₅₀ difference per patient was 0.4-0.6 °C. The risk for GG-associated treatment-limiting hot spots appeared low. Normal tissue high-temperature regions mostly occurred anteriorly; their volume and maximum temperature showed moderate positive correlations with GG volume, while fat-muscle interfaces were associated with higher risks for hot spots.

CONCLUSIONS

Considerable changes in volume and position of gastrointestinal gas can occur and are associated with clinically relevant tumor temperature differences.

Impact of Number of Segmented Tissues on SAR Prediction Accuracy in Deep Pelvic Hyperthermia Treatment Planning

Simple Summary

Hyperthermia treatment planning is the process of optimizing treatment quality using pre-treatment simulations. Although it has become a powerful tool, prediction accuracy is strongly dependent on the patient model. For deep hyperthermia in the pelvis, it is common that only four tissue categories are discriminated (bone, fat, muscle-like, and tumor). For the head and neck region, more tissues have been shown to be required for good prediction accuracy. Delineating is a labor-intensive and difficult process. Hence, it is important to find the optimum between accuracy and labor, but for deep pelvic hyperthermia, there are no published studies showing the impact of the number of tissues. We studied the trade-off between the segmentation detail needed and segmentation feasibility. Our findings indicate that including high water content tissues can impact simulation accuracy. Although our results, in general, underline the suitability of our current clinical protocol, they help to prioritize improvements for specific cases.

Abstract

In hyperthermia, the general opinion is that pre-treatment optimization of treatment settings requires a patient-specific model. For deep pelvic hyperthermia treatment planning (HTP), tissue models comprising four tissue categories are currently discriminated. For head and neck HTP, we found that more tissues are required for increasing accuracy. In this work, we evaluated the impact of the number of segmented tissues on the predicted specific absorption rate (SAR) for the pelvic region. Highly detailed anatomical models of five healthy volunteers were selected from a virtual database. For each model, seven lists with varying levels of segmentation detail were defined and used as an input for a modeling study. SAR changes were quantified using the change in target-to-hotspot-quotient and maximum SAR relative differences, with respect to the most detailed patient model. The main finding of this study was that the inclusion of high water content tissues in the segmentation may result in a clinically relevant impact on the SAR distribution and on the predicted hyperthermia treatment quality when considering our pre-established thresholds. In general, our results underline the current clinical segmentation protocol and help to prioritize any improvements.

Integrating Loco-Regional Hyperthermia Into the Current Oncology Practice: SWOT and TOWS Analyses

Moderate hyperthermia at temperatures between 40 and 44°C is a multifaceted therapeutic modality. It is a potent radiosensitizer, interacts favorably with a host of chemotherapeutic agents, and, in combination with radiotherapy, enforces immunomodulation akin to “in situ tumor vaccination.” By sensitizing hypoxic tumor cells and inhibiting repair of radiotherapy-induced DNA damage, the properties of hyperthermia delivered together with photons might provide a tumor-selective therapeutic advantage analogous to high linear energy transfer (LET) neutrons, but with less normal tissue toxicity. Furthermore, the high LET attributes of hyperthermia thermoradiobiologically are likely to enhance low LET protons; thus, proton thermoradiotherapy would mimic ¹²C ion therapy. Hyperthermia with radiotherapy and/or chemotherapy substantially improves therapeutic outcomes without enhancing normal tissue morbidities, yielding level I evidence reported in several randomized clinical trials, systematic reviews, and meta-analyses for various tumor sites. Technological advancements in hyperthermia delivery, advancements in hyperthermia treatment planning, online invasive and non-invasive MR-guided thermometry, and adherence to quality assurance guidelines have ensured safe and effective delivery of hyperthermia to the target region. Novel biological modeling permits integration of hyperthermia and radiotherapy treatment plans. Further, hyperthermia along with immune checkpoint inhibitors and DNA damage repair inhibitors could further augment the therapeutic efficacy resulting in synthetic lethality. Additionally, hyperthermia induced by magnetic nanoparticles coupled to selective payloads, namely, tumor-specific radiotheranostics (for both tumor imaging and radionuclide therapy), chemotherapeutic drugs, immunotherapeutic agents, and gene silencing, could provide a comprehensive tumor-specific theranostic modality akin to “magic (nano)bullets.” To get a realistic overview of the strength (S), weakness (W), opportunities (O), and threats (T) of hyperthermia, a SWOT analysis has been undertaken. Additionally, a TOWS analysis categorizes future strategies to facilitate further integration of hyperthermia with the current treatment modalities. These could gainfully accomplish a safe, versatile, and cost-effective enhancement of the existing therapeutic armamentarium to improve outcomes in clinical oncology.

Role of Simulations in the Treatment Planning of Radiofrequency Hyperthermia Therapy in Clinics

Hyperthermia therapy is a treatment modality in which tumor temperatures are elevated to higher temperatures to cause damage to cancerous tissues. Numerical simulations are integral in the development of hyperthermia treatment systems and in clinical treatment planning. In this study, simulations in radiofrequency hyperthermia therapy are reviewed in terms of their technical development and clinical aspects for effective clinical use. This review offers an overview of mathematical models and the importance of tissue properties; locoregional mild hyperthermia therapy, including phantom and realistic human anatomy models; phase array systems; tissue damage; thermal dose analysis; and thermoradiotherapy planning. This review details the improvements in numerical approaches in treatment planning and their application for effective clinical use. Furthermore, the modeling of thermoradiotherapy planning, which can be integrated with radiotherapy to provide combined hyperthermia and radiotherapy treatment planning strategies, are also discussed. This review may contribute to the effective development of thermoradiotherapy planning in clinics.

Confirmation of thermal dose as a predictor of local control in cervical carcinoma patients treated with state-of-the-art radiation therapy and hyperthermia

BACKGROUND

Addition of deep hyperthermia results in improved local control (LC) and overall survival (OS) compared to radiotherapy alone in patients with cervical carcinoma. Previously, we showed that the thermal dose of hyperthermia significantly correlates with LC and disease specific survival (DSS). Over the last decade, new radiation techniques were introduced resulting in improved LC.

AIM

To validate the effect of thermal dose in a more recent cohort of patients treated with modern radiotherapy techniques, including image guided brachytherapy (IGBT).

METHODS

We analyzed primary cervical carcinoma patients treated with a combination of radiotherapy and deep hyperthermia between 2005 and 2016 at our institute. Data on patient, tumor and treatment were collected including the thermal dose parameters TRISE and CEM43T90. Follow-up data on LC, disease free survival, DSS, OS as well as late toxicity data were collected. Data were analyzed using the Cox proportional hazard and Kaplan-Meier analyses.

RESULTS

227 patients were included. In multivariate analysis, histology, FIGO stage, lymphadenopathy, TRISE, CEM43T90 and IGBT had a significant effect on LC. In the patients treated with IGBT, the thermal dose parameter TRISE remained to have a significant effect on LC in univariate analysis.

CONCLUSIONS

The positive association between thermal dose and clinical outcome is replicated in an independent, recent cohort of cervical carcinoma patients. Importantly, in patients receiving IGBT, the effect of thermal dose on clinical outcome is still observed.

Systematic quality assurance of the BSD2000-3D MR-compatible hyperthermia applicator performance using MR temperature imaging

BACKGROUND

Radiofrequency (RF) mild hyperthermia (40 °C-44 °C for 60 minutes) is an effective adjuvant treatment for several types of cancer. To ensure treatment efficacy, quality assurance (QA) is necessary. This study presents the first systematic 3D characterisation of the heating performance of the commonly used Pyrexar BSD2000-3D MR-compatible hyperthermia applicator using magnetic resonance temperature imaging (MRTI).

METHODS

A reproducibly positioned phantom was heated with a power of 1000 watts during the 12.4 min needed to measure eight temperature distributions using MRTI. The target heating location was systematically varied between experiments. We analysed focus shape characteristics, steering accuracy, focus deformation due to steering, presence of off-target heating and reproducibility.

RESULTS

The mean maximum temperature increase was 5.9 ± 0.4 °C. The mean full width half maximum (FWHM) was 14.4 ± 0.5 cm in the XY plane and 24.5 ± 0.8 cm in Z-direction. The mean steering error was 0.4 ± 0.2 cm. The focus shape slightly varied between experiments, depending on steering distance in Y-direction. Off-target heating was not detected. Reproducibility of the focus amplitude and shape was determined by comparing the mean deviation from the mean temperature in the central slice was 0.3 ± 0.2 °C.

CONCLUSION

The Pyrexar BSD2000-3D MR-compatible applicator provides robust and reproducible heating. The upper boundary of the 95% confidence interval of the spatial steering accuracy is 0.9 cm, i.e. sufficient to fulfil the criterion of ≤0.2 °C temperature variation due to positioning errors as defined by Canters et al.

Online Adaptive Hyperthermia Treatment Planning During Locoregional Heating to Suppress Treatment-Limiting Hot Spots

BACKGROUND

Adequate tumor temperatures during hyperthermia are essential for good clinical response, but excessive heating of normal tissue should be avoided. This makes locoregional heating using phased array systems technically challenging. Online application of hyperthermia treatment planning could help to improve the heating quality. The aim of this study was to evaluate the clinical benefit of online treatment planning during treatment of pelvic tumors heated with the AMC-8 locoregional hyperthermia system.

METHODS

For online adaptive hyperthermia treatment planning, a graphical user interface was developed. Electric fields were calculated in a preprocessing step using our in-house-developed finite-difference-based treatment planning system. This allows instant calculation of the temperature distribution for user-selected phase-amplitude settings during treatment and projection onto the patient's computed tomographic scan for online visualization. Online treatment planning was used for 14 treatment sessions in 8 patients to reduce the patients' reports of hot spots while maintaining the same level of tumor heating. The predicted decrease in hot spot temperature should be at least 0.5°C, and the tumor temperature should decrease less than 0.2°C. These predictions were compared with clinical data: patient feedback about the hot spot and temperature measurements in the tumor region.

RESULTS

In total, 17 hot spot reports occurred during the 14 sessions, and the alternative settings predicted the hot spot temperature to decrease by at least 0.5°C, which was confirmed by the disappearance of all 17 hot spot reports. At the same time, the average tumor temperature was predicted to change on average -0.01°C (range, -0.19°C to 0.34°C). The measured tumor temperature change was on average only -0.02°C (range, -0.26°C to 0.31°C). In only 2 cases the temperature decrease was slightly larger than 0.2°C, but at most it was 0.26°C.

CONCLUSIONS

Online application of hyperthermia treatment planning is reliable and very useful to reduce hot spots without affecting tumor temperatures.

A three-dimensional visualization preoperative treatment planning system for microwave ablation in liver cancer: a simulated experimental study

PURPOSE

To evaluate the application value of three-dimensional (3D) visualization preoperative treatment planning system (VPTPS) for microwave ablation (MWA) in liver cancer.

MATERIALS AND METHODS

The study was a simulated experimental study using the CT imaging data of patients in DICOM format in a model. Three students (who learn to interventional ultrasound for less than 1 year) and three experts (who have more than 5 years of experience in ablation techniques) in MWA performed the preoperative planning for 39 lesions (mean diameter 3.75 ± 1.73 cm) of 32 patients using two-dimensional (2D) image planning method and 3D VPTPS, respectively. The number of planning insertions, planning ablation rate, and damage rate to surrounding structures were compared between2D image planning group and 3D VPTPS group.

RESULTS

There were fewer planning insertions, lower ablation rate and higher damage rate to surrounding structures in 2D image planning group than 3D VPTPS group for both students and experts. When using the 2D ultrasound planning method, students could carry out fewer planning insertions and had a lower ablation rate than the experts (p < 0.001). However, there was no significant difference in planning insertions, the ablation rate, and the incidence of damage to the surrounding structures between students and experts using 3D VPTPS.

CONCLUSIONS

3DVPTPS enables inexperienced physicians to have similar preoperative planning results to experts, and enhances students' preoperative planning capacity, which may improve the therapeutic efficacy and reduce the complication of MWA.

A printed Yagi-Uda antenna for application in magnetic resonance thermometry guided microwave hyperthermia applicators

Biological studies and clinical trials show that addition of hyperthermia stimulates conventional cancer treatment modalities and significantly improves treatment outcome. This supra-additive stimulation can be optimized by adaptive hyperthermia to counteract strong and dynamic thermoregulation. The only clinically proven method for the 3D non-invasive temperature monitoring required is by magnetic resonance (MR) temperature imaging, but the currently available set of MR compatible hyperthermia applicators lack the degree of heat control required. In this work, we present the design and validation of a high-frequency (433 MHz ISM band) printed circuit board antenna with a very low MR-footprint. This design is ideally suited for use in a range of hyperthermia applicator configurations. Experiments emulating the clinical situation show excellent matching properties of the antenna over a 7.2% bandwidth (S ₁₁  <  -15 dB). Its strongly directional radiation properties minimize inter-element coupling for typical array configurations (S ₂₁  <  -23 dB). MR imaging distortion by the antenna was found negligible and MR temperature imaging in a homogeneous muscle phantom was highly correlated with gold-standard probe measurements (root mean square error: RMSE  =  0.51 °C and R ²  =  0.99). This work paves the way for tailored MR imaging guided hyperthermia devices ranging from single antenna or incoherent antenna-arrays, to real-time adaptive hyperthermia with phased-arrays.

Absence of acute ocular damage in humans after prolonged exposure to intense RF EMF

The eye is considered to be a critical organ when determining safety standards for radio frequency (RF) radiation. Experimental data obtained using animals showed that RF heating of the eye, particularly over a specific threshold, can induce cataracts. During the treatment of cancer in the head and neck by hyperthermia, the eyes receive a considerable dose of RF radiation due to stray radiation from the prolonged (60 min) and intense exposure at 434 MHz of this region. In the current study, we verified the exposure guidelines for humans by determining the association between the electromagnetic and thermal dose in the eyes with the reported ocular effects. We performed a simulation study to retrospectively assess the specific absorption rate (SAR) and temperature increase in the eyes of 16 selected patients (encompassing a total of 74 treatment sessions) whose treatment involved high power delivery as well as a minimal distance between the tumor site and the eye. Our results show that the basic restrictions on the peak 10 g spatial-averaged SAR (10 W kg(-1)) and peak tissue temperature increase (1 °C) are exceeded by up to 10.4 and 4.6 times, on average, and by at least 6.2 and 1.8 times when considering the lower limit of the 95% confidence interval. Evaluation of the acute effects according to patients' feedback (all patients), the common toxicity criteria scores (all patients) and an ophthalmology investigation (one patient with the highest exposure) revealed no indication of any serious acute ocular effect, even though the eyes were exposed to high electromagnetic fields, leading to a high thermal dose. We also found that, although there is a strong correlation (R (2) =  0.88) between the predicted induced SAR and temperature in the eye, there are large uncertainties regarding the temperature-SAR relationship. Given this large uncertainty (129%) compared with the uncertainty of 3D temperature simulations (61%), we recommend using temperature simulations as a dosimetric measure in electromagnetic exposure risk assessments.

Molecular Imaging-Guided Interventional Hyperthermia in Treatment of Breast Cancer

Breast cancer is the most frequent malignancy in women worldwide. Although it is commonly treated via chemotherapy, responses vary among its subtypes, some of which are relatively insensitive to chemotherapeutic drugs. Recent studies have shown that hyperthermia can enhance the effects of chemotherapy in patients with refractory breast cancer or without surgical indications. Recent advances in molecular imaging may not only improve early diagnosis but may also facilitate the development and response assessment of targeted therapies. Combining advanced techniques such as molecular imaging and hyperthermia-integrated chemotherapy should open new avenues for effective management of breast cancer.

Current state of the art of regional hyperthermia treatment planning: a review

Locoregional hyperthermia, i.e. increasing the tumor temperature to 40–45 °C using an external heating device, is a very effective radio and chemosensitizer, which significantly improves clinical outcome. There is a clear thermal dose-effect relation, but the pursued optimal thermal dose of 43 °C for 1 h can often not be realized due to treatment limiting hot spots in normal tissue. Modern heating devices have a large number of independent antennas, which provides flexible power steering to optimize tumor heating and minimize hot spots, but manual selection of optimal settings is difficult. Treatment planning is a very valuable tool to improve locoregional heating. This paper reviews the developments in treatment planning software for tissue segmentation, electromagnetic field calculations, thermal modeling and optimization techniques. Over the last decade, simulation tools have become more advanced. On-line use has become possible by implementing algorithms on the graphical processing unit, which allows real-time computations. The number of applications using treatment planning is increasing rapidly and moving on from retrospective analyses towards assisting prospective clinical treatment strategies. Some clinically relevant applications will be discussed.

Local hyperthermia combined with radiotherapy and-/or chemotherapy: recent advances and promises for the future

Hyperthermia, one of the oldest forms of cancer treatment involves selective heating of tumor tissues to temperatures ranging between 39 and 45°C. Recent developments based on the thermoradiobiological rationale of hyperthermia indicate it to be a potent radio- and chemosensitizer. This has been further corroborated through positive clinical outcomes in various tumor sites using thermoradiotherapy or thermoradiochemotherapy approaches. Moreover, being devoid of any additional significant toxicity, hyperthermia has been safely used with low or moderate doses of reirradiation for retreatment of previously treated and recurrent tumors, resulting in significant tumor regression. Recent in vitro and in vivo studies also indicate a unique immunomodulating prospect of hyperthermia, especially when combined with radiotherapy. In addition, the technological advances over the last decade both in hardware and software have led to potent and even safer loco-regional hyperthermia treatment delivery, thermal treatment planning, thermal dose monitoring through noninvasive thermometry and online adaptive temperature modulation. The review summarizes the outcomes from various clinical studies (both randomized and nonrandomized) where hyperthermia is used as a thermal sensitizer of radiotherapy and-/or chemotherapy in various solid tumors and presents an overview of the progresses in loco-regional hyperthermia. These recent developments, supported by positive clinical outcomes should merit hyperthermia to be incorporated in the therapeutic armamentarium as a safe and an effective addendum to the existing oncological treatment modalities.

Tissue segmentation of head and neck CT images for treatment planning: a multiatlas approach combined with intensity modeling

PURPOSE

Hyperthermia treatment of head and neck tumors requires accurate treatment planning, based on 3D patient models that are derived from segmented 3D images. These segmentations are currently obtained by manual outlining of the relevant tissue regions, which is a tedious and time-consuming procedure (≈ 8 h) limiting the clinical applicability of hyperthermia treatment. In this context, the authors present and evaluate an automatic segmentation algorithm for CT images of the head and neck.

METHODS

The proposed method combines anatomical information, based on atlas registration, with local intensity information in a graph cut framework. The method is evaluated with respect to ground truth manual delineation and compared with multiatlas-based segmentation on a dataset of 18 labeled CT images using the Dice similarity coefficient (DSC), the mean surface distance (MSD), and the Hausdorff surface distance (HSD) as evaluation measures. On a subset of 13 labeled images, the influence of different labelers on the method's accuracy is quantified and compared with the interobserver variability.

RESULTS

For the DSC, the proposed method performs significantly better for the segmentation of all the tissues, except brain stem and spinal cord. The MSD shows a significant improvement for optical nerve, eye vitreous humor, lens, and thyroid. For the HSD, the proposed method performs significantly better for eye vitreous humor and brainstem. The proposed method has a significantly better score for DSC, MSD, and HSD than the multiatlas-based method for the eye vitreous humor. For the majority of the tissues (8/11) the segmentation accuracy of the proposed method is approaching the interobserver agreement. The authors' method showed better robustness to variations in atlas labeling compared with multiatlas segmentation. Moreover, the method improved the segmentation reproducibility compared with human observer's segmentations.

CONCLUSIONS

In conclusion, the proposed framework provides in an accurate automatic segmentation of head and neck tissues in CT images for the generation of 3D patient models, which improves reproducibility, and substantially reduces labor involved in therapy planning.

Thermal modelling using discrete vasculature for thermal therapy: A review

Reliable temperature information during clinical hyperthermia and thermal ablation is essential for adequate treatment control, but conventional temperature measurements do not provide 3D temperature information. Treatment planning is a very useful tool to improve treatment quality, and substantial progress has been made over the last decade. Thermal modelling is a very important and challenging aspect of hyperthermia treatment planning. Various thermal models have been developed for this purpose, with varying complexity. Since blood perfusion is such an important factor in thermal redistribution of energy in in vivo tissue, thermal simulations are most accurately performed by modelling discrete vasculature. This review describes the progress in thermal modelling with discrete vasculature for the purpose of hyperthermia treatment planning and thermal ablation. There has been significant progress in thermal modelling with discrete vasculature. Recent developments have made real-time simulations possible, which can provide feedback during treatment for improved therapy. Future clinical application of thermal modelling with discrete vasculature in hyperthermia treatment planning is expected to further improve treatment quality.

Targeted thermosensitive liposomes: an attractive novel approach for increased drug delivery to solid tumors

INTRODUCTION

Currently available chemotherapy is hampered by a lack in tumor specificity and resulting toxicity. Small and long-circulating liposomes can preferentially deliver chemotherapeutic drugs to tumors upon extravasation from tumor vasculature. Although clinically used liposomal formulations demonstrated significant reduction in toxicity, enhancement of therapeutic activity has not fully met expectations.

AREAS COVERED

Low drug bioavailability from liposomal formulations and limited tumor accumulation remain major challenges to further improve therapeutic activity of liposomal chemotherapy. The aim of this review is to highlight strategies addressing these challenges. A first strategy uses hyperthermia and thermosensitive liposomes to improve tumor accumulation and trigger liposomal drug bioavailability. Image-guidance can aid online monitoring of heat and drug delivery and further personalize the treatment. A second strategy involves tumor-specific targeting to enhance drug delivery specificity and drug internalization. In addition, we review the potential of combinations of the two in one targeted thermosensitive-triggered drug delivery system.

EXPERT OPINION

Heat-triggered drug delivery using thermosensitive liposomes as well as the use of tumor vasculature or tumor cell-targeted liposomes are both promising strategies to improve liposomal chemotherapy. Preclinical evidence has been encouraging and both strategies are currently undergoing clinical evaluation. A combination of both strategies rendering targeted thermosensitive liposomes (TTSL) may appear as a new and attractive approach promoting tumor drug delivery.

Simulation techniques in hyperthermia treatment planning

Abstract Clinical trials have shown that hyperthermia (HT), i.e. an increase of tissue temperature to 39-44 °C, significantly enhance radiotherapy and chemotherapy effectiveness [1]. Driven by the developments in computational techniques and computing power, personalised hyperthermia treatment planning (HTP) has matured and has become a powerful tool for optimising treatment quality. Electromagnetic, ultrasound, and thermal simulations using realistic clinical set-ups are now being performed to achieve patient-specific treatment optimisation. In addition, extensive studies aimed to properly implement novel HT tools and techniques, and to assess the quality of HT, are becoming more common. In this paper, we review the simulation tools and techniques developed for clinical hyperthermia, and evaluate their current status on the path from 'model' to 'clinic'. In addition, we illustrate the major techniques employed for validation and optimisation. HTP has become an essential tool for improvement, control, and assessment of HT treatment quality. As such, it plays a pivotal role in the quest to establish HT as an efficacious addition to multi-modality treatment of cancer.

Figures of merit and their bounds in radiofrequency heating by phased arrays

PURPOSE

The problem of effective power delivery to a semi-deep target by a phased array has been addressed for application to hyperthermia treatment of some tumours in the thorax.

METHODS

Three efficiencies have been introduced, which estimate system ability in power transfer from generators to body, from body to tumour, and from generators to tumour. They are formulated in terms of a dissipation matrix and an interference matrix. Bounds to achievable efficiencies are obtained. Further figures of merit have also been introduced. The necessary mathematics has been developed.

RESULTS

A numerical analysis has been carried out for a partially interdigitated planar array of resonant dipoles. Results show how the new parameters can be exploited for optimal selection of the array's degrees of freedom.

CONCLUSION

The figures of merit and their bounds allow comparisons between RF heating devices and provide guidelines to phased array design.

Computational modeling of temperature elevation and thermoregulatory response in the brains of anesthetized rats locally exposed at 1.5 GHz

The dominant effect of human exposures to microwaves is caused by temperature elevation ('thermal effect'). In the safety guidelines/standards, the specific absorption rate averaged over a specific volume is used as a metric for human protection from localized exposure. Further investigation on the use of this metric is required, especially in terms of thermophysiology. The World Health Organization (2006 RF research agenda) has given high priority to research into the extent and consequences of microwave-induced temperature elevation in children. In this study, an electromagnetic-thermal computational code was developed to model electromagnetic power absorption and resulting temperature elevation leading to changes in active blood flow in response to localized 1.457 GHz exposure in rat heads. Both juvenile (4 week old) and young adult (8 week old) rats were considered. The computational code was validated against measurements for 4 and 8 week old rats. Our computational results suggest that the blood flow rate depends on both brain and core temperature elevations. No significant difference was observed between thermophysiological responses in 4 and 8 week old rats under these exposure conditions. The computational model developed herein is thus applicable to set exposure conditions for rats in laboratory investigations, as well as in planning treatment protocols in the thermal therapy.

Procedure for creating a three-dimensional (3D) model for superficial hyperthermia treatment planning

PURPOSE

To make a patient- and treatment-specific computed tomography (CT) scan and to create a three-dimensional (3D) patient model for superficial hyperthermia treatment planning (SHTP).

PATIENTS, MATERIALS, AND METHODS

Patients with recurrent breast adenocarcinoma in previously irradiated areas referred for radiotherapy (RT) and hyperthermia (HT) treatment and giving informed consent were included. After insertion of the thermometry catheters in the treatment area, a CT scan in the treatment position was made.

RESULTS

A total of 26 patients have been, thus far, included in the study. During the study period, five types of adjustments were made to the procedure: (1) marking the RT field with radioopaque markers, (2) making the CT scan after the first HT treatment instead of before, (3) using an air- and foam-filled (dummy) water bolus, (4) a change to radiolucent catheters for which radioopaque markers were needed, and (5) marking the visible/palpable extent of the tumor with radioopaque markers, if necessary. With these adjustments, all necessary information is visible on the CT scan. Each CT slice was automatically segmented into muscle, fat, bone, and air. RT field, catheters, applicators, and tumor lesions, if indicated, were outlined manually using the segmentation program iSeg. Next the model was imported into SEMCAD X, a 3D electromagnetic field simulator.

CONCLUSION

Using the final procedure to obtain a patient- and treatment-specific CT scan, it is possible to create a 3D model for SHTP.

Improved power steering with double and triple ring waveguide systems: the impact of the operating frequency

INTRODUCTION

Regional hyperthermia systems with 3D power steering have been introduced to improve tumour temperatures. The 3D 70-MHz AMC-8 system has two rings of four waveguides. The aim of this study is to evaluate whether T(90) will improve by using a higher operating frequency and whether further improvement is possible by adding a third ring.

METHODS

Optimised specific absorption rate (SAR) distributions were evaluated for a centrally located target in tissue-equivalent phantoms, and temperature optimisation was performed for five cervical carcinoma patients with constraints to normal tissue temperatures. The resulting T(90) and the thermal iso-effect dose (i.e. the number of equivalent min at 43°C) were evaluated and compared to the 2D 70-MHz AMC-4 system with a single ring of four waveguides. FDTD simulations were performed at 2.5 × 2.5 × 5 mm(3) resolution. The applied frequencies were 70, 100, 120, 130, 140 and 150 MHz.

RESULTS

Optimised SAR distributions in phantoms showed an optimal SAR distribution at 140 MHz. For the patient simulations, an optimal increase in T(90) was observed at 130 MHz. For a two-ring system at 70 MHz the gain in T(90) was about 0.5°C compared to the AMC-4 system, averaged over the five patients. At 130 MHz the average gain in T(90) was ~1.5°C and ~2°C for a two and three-ring system, respectively. This implies an improvement of the thermal iso-effect dose with a factor ~12 and ~30, respectively.

CONCLUSION

Simulations showed that a 130-MHz two-ring waveguide system yields significantly higher tumour temperatures compared to 70-MHz single-ring and double-ring waveguide systems. Temperatures were further improved with a 130-MHz triple-ring system.

Children and adults exposed to electromagnetic fields at the ICNIRP reference levels: theoretical assessment of the induced peak temperature increase

To avoid potentially adverse health effects of electromagnetic fields (EMF), the International Commission on Non-Ionizing Radiation Protection (ICNIRP) has defined EMF reference levels. Restrictions on induced whole-body-averaged specific absorption rate (SAR(wb)) are provided to keep the whole-body temperature increase (T(body, incr)) under 1 °C during 30 min. Additional restrictions on the peak 10 g spatial-averaged SAR (SAR(10g)) are provided to prevent excessive localized tissue heating. The objective of this study is to assess the localized peak temperature increase (T(incr, max)) in children upon exposure at the reference levels. Finite-difference time-domain modeling was used to calculate T(incr, max) in six children and two adults exposed to orthogonal plane-wave configurations. We performed a sensitivity study and Monte Carlo analysis to assess the uncertainty of the results. Considering the uncertainties in the model parameters, we found that a peak temperature increase as high as 1 °C can occur for worst-case scenarios at the ICNIRP reference levels. Since the guidelines are deduced from temperature increase, we used T(incr, max) as being a better metric to prevent excessive localized tissue heating instead of localized peak SAR. However, we note that the exposure time should also be considered in future guidelines. Hence, we advise defining limits on T(incr, max) for specified durations of exposure.