Simulation techniques in hyperthermia treatment planning

Magnetic Particle Imaging-Guided Thermal Simulations for Magnetic Particle Hyperthermia

Magnetic particle hyperthermia (MPH) enables the direct heating of solid tumors with alternating magnetic fields (AMFs). One challenge with MPH is the unknown particle distribution in tissue after injection. Magnetic particle imaging (MPI) can measure the nanoparticle content and distribution in tissue after delivery. The objective of this study was to develop a clinically translatable protocol that incorporates MPI data into finite element calculations for simulating tissue temperatures during MPH. To verify the protocol, we conducted MPH experiments in tumor-bearing mouse cadavers. Five 8-10-week-old female BALB/c mice bearing subcutaneous 4T1 tumors were anesthetized and received intratumor injections of Synomag®-S90 nanoparticles. Immediately following injection, the mice were euthanized and imaged, and the tumors were heated with an AMF. We used the Mimics Innovation Suite to create a 3D mesh of the tumor from micro-computerized tomography data and spatial index MPI to generate a scaled heating function for the heat transfer calculations. The processed imaging data were incorporated into a finite element solver, COMSOL Multiphysics®. The upper and lower bounds of the simulated tumor temperatures for all five cadavers demonstrated agreement with the experimental temperature measurements, thus verifying the protocol. These results demonstrate the utility of MPI to guide predictive thermal calculations for MPH treatment planning.

Patient-specific temperature distribution prediction in laser interstitial thermal therapy: single-irradiation data-driven method

Laser interstitial thermal therapy (LITT) is popular for treating brain tumours and epilepsy. The strict control of tissue thermal damage extent is crucial for LITT. Temperature prediction is useful for predicting thermal damage extent. Accurately predictingin vivobrain tissue temperature is challenging due to the temperature dependence and the individual variations in tissue properties. Considering these factors is essential for improving the temperature prediction accuracy.Objective. To present a method for predicting patient-specific tissue temperature distribution within a target lesion area in the brain during LITT.Approach. A magnetic resonance temperature imaging (MRTI) data-driven estimation model was constructed and combined with a modified Pennes bioheat transfer equation (PBHE) to predict patient-specific temperature distribution. In the PBHE for temperature prediction, the individual specificity and temperature dependence of thermal tissue properties and blood perfusion, as well as the individual specificity of optical tissue properties were considered. Only MRTI data during one laser irradiation were required in the method. This enables the prediction of patient-specific temperature distribution and the resulting thermal damage region for subsequent ablations.Main results. Patient-specific temperature prediction was evaluated based on clinical data acquired during LITT in the brain, using intraoperative MRTI data as the reference standard. Our method significantly improved the prediction performance of temperature distribution and thermal damage region. The average root mean square error was decreased by 69.54%, the average intraclass correlation coefficient was increased by 37.5%, the average Dice similarity coefficient was increased by 43.14% for thermal damage region prediction.Significance. The proposed method can predict temperature distribution and thermal damage region at an individual patient level during LITT, providing a promising approach to assist in patient-specific treatment planning for LITT in the brain.

A Validated Methodological Approach to Prove the Safety of Clinical Electromagnetic Induction Systems in Magnetic Hyperthermia

The present study focuses on the development of a methodology for evaluating the safety of MNH systems, through the numerical prediction of the induced temperature rise in superficial skin layers due to eddy currents heating under an alternating magnetic field (AMF). The methodology is supported and validated through experimental measurements of the AMF's distribution, as well as temperature data from the torsos of six patients who participated in a clinical trial study. The simulations involved a computational model of the actual coil, a computational model of the cooling system used for the cooling of the patients during treatment, and a detailed human anatomical model from the Virtual Population family. The numerical predictions exhibit strong agreement with the experimental measurements, and the deviations are below the estimated combined uncertainties, confirming the accuracy of computational modeling. This study highlights the crucial role of simulations for translational medicine and paves the way for personalized treatment planning.

Assessment of the thermal tissue models for the head and neck hyperthermia treatment planning

PURPOSE

To compare different thermal tissue models for head and neck hyperthermia treatment planning, and to assess the results using predicted and measured applied power data from clinical treatments.

METHODS

Three commonly used temperature models from literature were analysed: "constant baseline", "constant thermal stress" and "temperature dependent". Power and phase data of 93 treatments of 20 head and neck patients treated with the HYPERcollar3D applicator were used. The impact on predicted median temperature T50 inside the target region was analysed with maximum allowed temperature of 44 °C in healthy tissue. The robustness of predicted T50 for the three models against the influence of blood perfusion, thermal conductivity and the assumed hotspot temperature level was analysed.

RESULTS

We found an average predicted T50 of 41.0 ± 1.3 °C (constant baseline model), 39.9 ± 1.1 °C (constant thermal stress model) and 41.7 ± 1.1 °C (temperature dependent model). The constant thermal stress model resulted in the best agreement between the predicted power (P = 132.7 ± 45.9 W) and the average power measured during the hyperthermia treatments (P = 129.1 ± 83.0 W).

CONCLUSION

The temperature dependent model predicts an unrealistically high T50. The power values for the constant thermal stress model, after scaling simulated maximum temperatures to 44 °C, matched best to the average measured powers. We consider this model to be the most appropriate for temperature predictions using the HYPERcollar3D applicator, however further studies are necessary for developing of robust temperature model for tissues during heat stress.

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.

Three dimensional models of human thermoregulation: A review

Numerous human thermoregulatory models have been developed and widely used in various applications such as aerospace, medicine, public health, and physiology research. This paper is a review of three dimensional (3D) models for human thermoregulation. This review begins with a short introduction of thermoregulatory model development followed by key principles for mathematical description of human thermoregulation systems. Different representations of 3D human bodies are discussed with respect to their detail and prediction capability. The human body was divided into fifteen layered cylinders in early 3D models (cylinder model). Recent 3D models have utilized medical image datasets to develop geometrically correct human models (realistic geometry model). The finite element method is mostly used to solve the governing equations and get numerical solutions. The realistic geometry models provide a high degree of anatomical realism and predict whole-body thermoregulatory responses at high resolution and at organ and tissue levels. Thus, 3D models extend to a wide range of applications where temperature distribution is critical, such as hypothermia/hyperthermia therapy and physiology research. The development of thermoregulatory models will continue with the growth in computational power, advancement in numerical methods and simulation software, advances in modern imaging techniques, and progress in the basic science of thermal physiology.

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.

Application of photogrammetry reconstruction for hyperthermia quality control measurements

PURPOSE

Hyperthermia is a cancer treatment in which the target region is heated to temperatures of 40-44 °C usually applying external electromagnetic field sources. The behavior of the hyperthermia applicators (antennas) in clinical practice should be periodically checked with phantom experiments to verify the applicator's performance over time. The purpose of this study was to investigate the application of photogrammetry reconstructions of 3D applicator position in these quality control procedure measurements.

METHODS

Photogrammetry reconstruction was applied at superficial hyperthermia scenario using the Lucite cone applicator (LCA) and phased-array heating in the head and neck region using the HYPERcollar3D. Wire-frame models of the entire measurement setups were created from multiple-view images and used for recreation of the setup inside 3D electromagnetic field simulation software. We evaluated applicator relation (R_a) between measured and simulated absolute specific absorption rate (SAR) for manually created and photogrammetry reconstructed simulation setups.

RESULTS

We found a displacement of 7.9 mm for the LCA and 8.2 mm for the HYPERcollar3D setups when comparing manually created and photogrammetry reconstructed applicator models placements. R_a improved from 1.24 to 1.18 for the LCA and from 1.17 to 1.07 for the HYPERcollar3D when using photogrammetry reconstructed simulation setups.

CONCLUSION

Photogrammetry reconstruction technique holds promise to improve measurement setup reconstruction and agreement between measured and simulated absolute SAR.

Integrated thermal and magnetic susceptibility modeling for air-motion artifact correction in proton resonance frequency shift thermometry

PURPOSE

Hyperthermia treatments are successful adjuvants to conventional cancer therapies in which the tumor is sensitized by heating. To monitor and guide the hyperthermia treatment, measuring the tumor and healthy tissue temperature is important. The typical clinical practice heavily relies on intraluminal probe measurements that are uncomfortable for the patient and only provide spatially sparse temperature information. A solution may be offered through recent advances in magnetic resonance thermometry, which allows for three-dimensional internal temperature measurements. However, these measurements are not widely used in the pelvic region due to a low signal-to-noise ratio and presence of image artifacts.

METHODS

To advance the clinical integration of magnetic resonance-guided cancer treatments, we consider the problem of removing air-motion-induced image artifacts. Thereto, we propose a new combined thermal and magnetic susceptibility model-based temperature estimation scheme that uses temperature estimates to improve the removal of air-motion-induced image artifacts. The method is experimentally validated using a dedicated phantom that enables the controlled injection of air-motion artifacts and with in vivo thermometry from a clinical hyperthermia treatment.

RESULTS

We showed, using probe measurements in a heated phantom, that our method reduced the mean absolute error (MAE) by 58% compared to the state-of-the-art near a moving air volume. Moreover, with in vivo thermometry our method obtained a MAE reduction between 17% and 95% compared to the state-of-the-art.

CONCLUSION

We expect that the combined thermal and magnetic susceptibility modeling used in model-based temperature estimation can significantly improve the monitoring in hyperthermia treatments and enable feedback strategies to further improve MR-guided hyperthermia cancer treatments.

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.

Spatio-Spectral Ultrasound Characterization of Reflection and Transmission Through Bone With Temperature Dependence

Transcranial focused ultrasound (tFUS) is a promising approach for the treatment of neurological disorders. It has proven useful in several clinical applications, with promising outcomes reported in the recent literature. Furthermore, it is currently being investigated in a range of neuromodulation (NM) and ablative applications, including epilepsy. In this application, tFUS access through the temporal window is the key to optimizing the treatment safety and efficacy. Traditional approaches have utilized transducers with low operating frequencies for tFUS applications. Modern array transducers and driving systems allow for more intelligent use of the temporal window by exploiting the spatio-spectral transmission bandwidth to a specified target or targets within the brain. To demonstrate the feasibility of this approach, we have investigated the ultrasound reflection and transmission characteristics for different access points within the temporal window of human skull samples ex vivo. Different transmit-receive (Rx) configurations are used for characterization of the spatio-spectral variability in reflection and transmission through the temporal window. In this article, we show results from a dual-piston transducer set up in the frequency range of 2-7 MHz. Broadband pulses as well as synthesized orthogonal frequency division multiplexed (OFDM) waveforms were used. The latter was used to improve the magnitude and phase measurements in 100-kHz subbands within the 2-7 MHz spectral window. A temperature-controlled water bath was used to characterize the change in reflection and transmission characteristics with temperature in the 25°C-43°C range. The measured values of the complex reflection and transmission coefficients exhibited significant variations with space, frequency, and temperature. On the other hand, the measured transmission phase varied more with location and frequency, with smaller sensitivity to temperature. A measurement-based hybrid angular spectrum (HAS) simulation through the human temporal bone was used to demonstrate the dependence of focusing gain on the skull profile and spatial distribution of change of speed of sound (SOS) at different skull temperatures.

Influence of Magnetic Scaffold Loading Patterns on their Hyperthermic Potential against Bone Tumors

Magnetic scaffolds have been investigated as promising tools for the interstitial hyperthermia treatment of bone cancers, to control local recurrence by enhancing radio- and chemotherapy effectiveness. The potential of magnetic scaffolds motivates the development of production strategies enabling tunability of the resulting magnetic properties. Within this framework, deposition and drop-casting of magnetic nanoparticles on suitable scaffolds offer advantages such as ease of production and high loading, although these approaches are often associated with a non-uniform final spatial distribution of nanoparticles in the biomaterial. The implications and the influences of nanoparticle distribution on the final therapeutic application have not yet been investigated thoroughly. In this work, poly-caprolactone scaffolds are magnetized by loading them with synthetic magnetic nanoparticles through a drop-casting deposition and tuned to obtain different distributions of magnetic nanoparticles in the biomaterial. The physicochemical properties of the magnetic scaffolds are analyzed. The microstructure and the morphological alterations due to the reworked drop-casting process are evaluated and correlated to static magnetic measurements. THz tomography is used as an innovative investigation technique to derive the spatial distribution of nanoparticles. Finally, multiphysics simulations are used to investigate the influence on the loading patterns on the interstitial bone tumor hyperthermia treatment.

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.

Suitability of eigenvalue beam-forming for discrete multi-frequency hyperthermia treatment planning

PURPOSE

Thermal dose delivery in microwave hyperthermia for cancer treatment is expected to benefit from the introduction of ultra-wideband (UWB)-phased array applicators. A full exploitation of the combination of different frequencies to improve the deposition pattern is, however, a nontrivial problem. It is unclear whether the cost functions used for hyperthermia treatment planning (HTP) optimization in the single-frequency setting can be meaningfully extended to the UWB case.

METHOD

We discuss the ability of the eigenvalue (EV) and a novel implementation of iterative-EV (i-EV) beam-forming methods to fully exploit the available frequency spectrum when a discrete set of simultaneous operating frequencies is available for treatment. We show that the quadratic power deposition ratio solved by the methods can be maximized by only one frequency in the set, therefore rendering EV inadequate for UWB treatment planning. We further investigate whether this represents a limitation in two realistic test cases, comparing the thermal distributions resulting from EV and i-EV to those obtained by optimizing for other nonlinear cost functions that allow for multi-frequency.

RESULTS

The classical EV-based single-frequency HTP yields systematically lower target SAR deposition and temperature values than nonlinear HTP. In a larynx target, the proposed single-frequency i-EV scheme is able to compensate for this and reach temperatures comparable to those given by global nonlinear optimization. In a meninges target, the multi-frequency setting outperforms the single-frequency one, achieving better target coverage and 0 . 5 ∘ C higher T 90 in the tumor than single-frequency-based HTP.

CONCLUSIONS

Classical EV performs poorly in terms of resulting target temperatures. The proposed single-frequency i-EV scheme can be a viable option depending on the patient and tumor to be treated, as long as the proper operating frequency can be selected across a UWB range. Multi-frequency HTP can bring a considerable benefit in regions typically difficult to treat such as the brain.

Improved patient-specific hyperthermia planning based on parametrized electromagnetic and thermal models for the SIGMA-30 applicator

OBJECTIVE

To create an improved planning method for pediatric regional hyperthermia (RHT) using the SIGMA-30 applicator (SIGMA-30).

MATERIALS AND METHODS

An electromagnetic model of SIGMA-30 was generated for use with the finite-difference time-domain (FDTD) method. Applying special MATLAB-based algorithms, voxel models of a pediatric patient with pelvic rhabdomyosarcoma were created from Computed-Tomography (CT) contours for use with the FDTD method and the finite-difference (FD) method capable of using either temperature-independent or temperature-dependent perfusion models for solving the Bioheat Transfer Equation (BHTE). Patient models were parametrized regarding, first, the positioning in the applicator, second, the absorbed power range and, third, different perfusion models, resulting in the so-called Parametrized Treatment Models (PTMs). A novel dedicated optimization procedure was developed based on quantitative comparison of numerical calculations against temperature and power measurements from two RHT therapies.

RESULTS

Using measured data, a realistic absorbed power range in the patient model was estimated. Within this range, several FDTD and BHTE runs were performed and, applying the aforementioned optimization scheme, the best PTMs and perfusion models were identified for each therapy via a retrospective comparison with measurements in 14 temperature sensor positions: 5 in the tumor, 8 in rectum and one in bladder.

CONCLUSION

A novel dedicated optimization procedure for identification of suitable patient-specific electromagnetic and thermal models, which can be used for improved patient planning, was developed and evaluated by comparison with treatment-derived measurements using SIGMA-30. The optimization procedure can be extended to other hyperthermia applicators and to other patient types, including adults.

In silico assessment of collateral eddy current heating in biocompatible implants subjected to magnetic hyperthermia treatments

Purpose: Bearing partially or fully metallic passive implants represents an exclusion criterion for patients undergoing a magnetic hyperthermia procedure, but there are no specific studies backing this restrictive decision. This work assesses how the secondary magnetic field generated at the surface of two common types of prostheses affects the safety and efficiency of magnetic hyperthermia treatments of localized tumors. The paper also proposes the combination of a multi-criteria decision analysis and a graphical representation of calculated data as an initial screening during the preclinical risk assessment for each patient.Materials and methods: Heating of a hip joint and a dental implant during the treatment of prostate, colorectal and head and neck tumors have been assessed considering different external field conditions and exposure times. The Maxwell equations including the secondary field produced by metallic prostheses have been solved numerically in a discretized computable human model. The heat exchange problem has been solved through a modified version of the Pennes' bioheat equation assuming a temperature dependency of blood perfusion and metabolic heat, i.e. thermorregulation. The degree of risk has been assessed using a risk index with parameters coming from custom graphs plotting the specific absorption rate (SAR) vs temperature increase, and coefficients derived from a multi-criteria decision analysis performed following the MACBETH approach.Results: The comparison of two common biomaterials for passive implants - Ti6Al4V and CoCrMo - shows that both specific absorption rate (SAR) and local temperature increase are found to be higher for the hip prosthesis made by Ti6Al4V despite its lower electrical and thermal conductivity. By tracking the time evolution of temperature upon field application, it has been established that there is a 30 s delay between the time point for which the thermal equilibrium is reached at prostheses and tissues. Likewise, damage may appear in those tissues adjacent to the prostheses at initial stages of treatment, since recommended thermal thresholds are soon surpassed for higher field intensities. However, it has also been found that under some operational conditions the typical safety rule of staying below or attain a maximum temperature increase or SAR value is met.Conclusion: The current exclusion criterion for implant-bearing patients in magnetic hyperthermia should be revised, since it may be too restrictive for a range of the typical field conditions used. Systematic in silico treatment planning using the proposed methodology after a well-focused diagnostic procedure can aid the clinical staff to find the appropriate limits for a safe treatment window.

Towards personalized and versatile monitoring of temperature fields within heterogeneous tissues during laser therapies

Advancements in medical laser technology have paved the way for its widespread acceptance in a variety of treatments and procedures. Selectively targeting particular tissue structures with minimally invasive procedures limits the damage to surrounding tissue and allows for reduced post-procedural downtime. In many treatments that are hyperthermia-based, the efficiency depends on the achieved temperature within the targeted tissues. Current approaches for monitoring subdermal temperature distributions are either invasive, complex, or offer inadequate spatial resolution. Numerical studies are often therapy-tailored and source tissue parameters from the literature, lacking versatility and a tissue-specific approach. Here, we show a protocol that estimates the temperature distribution within the tissue based on a thermographic recording of its surface temperature evolution. It couples a time-dependent matching algorithm and thermal-diffusion-based model, while recognizing tissue-specific characteristics yielded by a fast calibration process. The protocol was employed during hyperthermic laser treatment performed ex-vivo on a heterogeneous porcine tissue, and in-vivo on a human subject. In both cases the calibrated thermal parameters correlate with the range of values reported by other studies. The matching algorithm sufficiently reproduced the temperature dynamics of heterogeneous tissue. The estimated temperature distributions within ex-vivo tissue were validated by simultaneous reference measurements, and the ones estimated in-vivo reveal a distribution trend that correlates well with similar studies. The presented method is versatile, supported by the protocol for tissue-specific tailoring, and can readily be implemented for temperature monitoring of various hyperthermia-based procedures by means of recording the surface temperature evolution with a miniature thermal camera implemented within a handheld laser scanner or similar.

Feasibility of Temperature Control by Electrical Impedance Tomography in Hyperthermia

We present a simulation study investigating the feasibility of electrical impedance tomography (EIT) as a low cost, noninvasive technique for hyperthermia (HT) treatment monitoring and adaptation. Temperature rise in tissues leads to perfusion and tissue conductivity changes that can be reconstructed in 3D by EIT to noninvasively map temperature and perfusion. In this study, we developed reconstruction methods and investigated the achievable accuracy of EIT by simulating HT treatmentlike scenarios, using detailed anatomical models with heterogeneous conductivity distributions. The impact of the size and location of the heated region, the voltage measurement signal-to-noise ratio, and the reference model personalization and accuracy were studied. Results showed that by introducing an iterative reconstruction approach, combined with adaptive prior regions and tissue-dependent penalties, planning-based reference models, measurement-based reweighting, and physics-based constraints, it is possible to map conductivity-changes throughout the heated domain, with an accuracy of around 5% and cm-scale spatial resolution. An initial exploration of the use of multifrequency EIT to separate temperature and perfusion effects yielded promising results, indicating that temperature reconstruction accuracy can be in the order of 1 ∘C. Our results suggest that EIT can provide valuable real-time HT monitoring capabilities. Experimental confirmation in real-world conditions is the next step.

Standardization of patient modeling in hyperthermia simulation studies: introducing the Erasmus Virtual Patient Repository

Purpose: Thermal dose-effect relations have demonstrated that clinical effectiveness of hyperthermia would benefit from more controlled heating of the tumor. Hyperthermia treatment planning (HTP) is a potent tool to study strategies enabling target conformal heating, but its accuracy is affected by patient modeling approximations. Homogeneous phantoms models are being used that do not match the body shape of patients in treatment position and often have unrealistic target volumes. As a consequence, simulation accuracy is affected, and performance comparisons are difficult. The aim of this study is to provide the first step toward standardization of HTP simulation studies in terms of patient modeling by introducing the Erasmus Virtual Patient Repository (EVPR): a virtual patient model database.Methods: Four patients with a tumor in the head and neck or the pelvis region were selected, and corresponding models were created using a clinical segmentation procedure. Using the Erasmus University Medical Center standard procedure, HTP was applied to these models and compared to HTP for commonly used surrogate models.Results: Although this study was aimed at presenting the EVPR database, our study illustrates that there is a non-negligible difference in the predicted SAR patterns between patient models and homogeneous phantom-based surrogate models. We further demonstrate the difference between actual and simplified target volumes being used today.Conclusion: Our study describes the EVPR for the research community as a first step toward standardization of hyperthermia simulation studies.

Dielectric Properties of Ovine Heart at Microwave Frequencies

Accurate knowledge of the dielectric properties of biological tissues is important in dosimetry studies and for medical diagnostic, monitoring and therapeutic technologies. In particular, the dielectric properties of the heart are used in numerical simulations of radiofrequency and microwave heart ablation. In one recent study, it was demonstrated that the dielectric properties of different components of the heart can vary considerably, contrary to previous literature that treated the heart as a homogeneous organ with measurements that ignored the anatomical location. Therefore, in this study, we record and report the dielectric properties of the heart as a heterogeneous organ. We measured the dielectric properties at different locations inside and outside of the heart over the 500 MHz to 20 GHz frequency range. Different parts of the heart were identified based on the anatomy of the heart and their function; they include the epicardium, endocardium, myocardium, exterior and interior surfaces of atrial appendage, and the luminal surface of the great vessels. The measured dielectric properties for each part of the heart are reported at both a single frequency (2.4 GHz), which is of interest in microwave medical applications, and as parameters of a broadband Debye model. The results show that in terms of dielectric properties, different parts of the heart should not be considered the same, with more than 25% difference in dielectric properties between some parts. The specific Debye models and single frequency dielectric properties from this study can be used to develop more detailed models of the heart to be used in electromagnetic modeling.

Simulation guided design of the MRcollar: a MR compatible applicator for deep heating in the head and neck region

PURPOSE

To develop a head and neck hyperthermia phased array system compatible with a 1.5 T magnetic resonance (MR) scanner for noninvasive thermometry.

METHODS

We designed a dielectric-parabolic-reflector antenna (DiPRA) based on a printed reflector backed dipole antenna and studied its predicted and measured performance in a flat configuration (30 mm thick water bolus and muscle equivalent layer). Thereafter, we designed a phased array applicator model ('MRcollar') consisting of 12 DiPRA modules placed on a radius of 180 mm. Theoretical heating performance of the MRcollar model was benchmarked against the current clinical applicator (HYPERcollar3D) using specific (3D) head and neck models of 28 treated patients. Lastly, we assessed the influence of the DiPRA modules on MR scanning quality.

RESULTS

The predicted and measured reflection coefficients (S₁₁) of the DiPRA module are below -20 dB. The maximum specific absorption rate (SAR) in the area under the antenna was 47% higher than for the antenna without encasing. Compared to the HYPERcollar3D, the MRcollar design incorporates 31% less demineralized water (-2.5 L), improves the predicted TC25 (target volume enclosed by 25% iso-SAR contour) by 4.1% and TC50 by 8.5%, while the target-to-hotspot quotient (THQ) is minimally affected (-1.6%). MR experiments showed that the DiPRA modules do not affect MR transmit/receive performance.

CONCLUSION

Our results suggest that head and neck hyperthermia delivery quality with the MRcollar can be maintained, while facilitating simultaneous noninvasive MR thermometry for treatment monitoring and control.

Optimization of the Clinical Setting Using Numerical Simulations of the Electromagnetic Field in an Obese Patient Model for Deep Regional Hyperthermia of an 8 MHz Radiofrequency Capacitively Coupled Device in the Pelvis

Simple Summary

Several randomized controlled trials have shown that the concurrent use of hyperthermia with radiotherapy results in a significant increase in local control rates. However, in studies that analyzed thermal parameters, the radiosensitizing effect required a favorable increase in tumor temperature. A good temperature increase is more difficult to achieve in deep-seated tumors than in superficial tumors. In this study, the reduction of subcutaneous fat overheating, which is a weak point of the deep heating of a capacitively coupled heating system, and the optimization of temperature distribution in the deep regional hyperthermia in the pelvis were investigated using electromagnetic field numerical simulations. In conclusion, the use of overlay boluses, the salt solution concentration in the overlay boluses, and the intergluteal cleft bolus insertion were found to be important for optimizing the temperature distribution. Further studies with numerical simulations based on the patient’s body shape and tumor location are expected.

Abstract

Background: The purpose of this study was to evaluate the effectiveness of the clinical setting for deep regional hyperthermia of an 8 MHz radiofrequency (RF) capacitively coupled device in the pelvis by using numerical simulations of the electromagnetic field. Methods: A three-dimensional patient model of cervical cancer of the uterus in an obese patient was reconstructed with computed tomography data. The specific absorption rate (SAR) and temperature distributions among the various heating settings were evaluated using numerical simulations. Results: The averaged SAR value of the deep target tumor was similar between with or without overlay boluses (OBs), and that of the subcutaneous fat (SF) at the edges of cooling boluses with OBs was lower than that of the SF without OBs. The use of OBs reduced the overheating of the SF. The 0.5% salt solution in the OB produced the least overheated areas outside the deep target tumor compared with the other concentrations. The insertion of the intergluteal cleft (IGC) bolus could improve the temperature concentration of the deep target tumor. Conclusions: The use of OBs and the salt solution concentration in the OB were important to optimize the temperature distribution. IGC bolus might contribute to temperature optimization. Further studies with individualized numerical simulations in each patient are expected.

Iterative time-reversal for multi-frequency hyperthermia

Time-reversal (TR) is a known wideband array beam-forming technique that has been suggested as a treatment planning alternative in deep microwave hyperthermia for cancer treatment. While the aim in classic TR is to focus the energy at a specific point within the target, no assumptions are made on secondary lobes that might arise in the healthy tissues. These secondary lobes, together with tissue heterogeneity, may result in hot-spots (HSs), which are known to limit the efficiency of the thermal dose delivery to the tumor. This paper proposes a novel wideband TR focusing method that iteratively shifts the focus away from HSs and towards cold-spots from an initial TR solution, a procedure that improves tumor coverage and reduces HSs. We verify this method on two different applicator topologies and several target volume configurations. The algorithm is deterministic and runs within seconds, enabling its use for real-time applications. At the same time, it yields results comparable to those obtained with global stochastic optimizers such as Particle Swarm.

Whither Magnetic Hyperthermia? A Tentative Roadmap

The scientific community has made great efforts in advancing magnetic hyperthermia for the last two decades after going through a sizeable research lapse from its establishment. All the progress made in various topics ranging from nanoparticle synthesis to biocompatibilization and in vivo testing have been seeking to push the forefront towards some new clinical trials. As many, they did not go at the expected pace. Today, fruitful international cooperation and the wisdom gain after a careful analysis of the lessons learned from seminal clinical trials allow us to have a future with better guarantees for a more definitive takeoff of this genuine nanotherapy against cancer. Deliberately giving prominence to a number of critical aspects, this opinion review offers a blend of state-of-the-art hints and glimpses into the future of the therapy, considering the expected evolution of science and technology behind magnetic hyperthermia.

Laser-induced optothermal response of gold nanoparticles: From a physical viewpoint to cancer treatment application

Gold nanoparticles (GNPs)-based photothermal therapy (PTT) is a promising minimally invasive thermal therapy for the treatment of focal malignancies. Although GNPs-based PTT has been known for over two decades and GNPs possess unique properties as therapeutic agents, the delivery of a safe and effective therapy is still an open question. This review aims at providing relevant and recent information on the usage of GNPs in combination with the laser to treat cancers, pointing out the practical aspects that bear on the therapy outcome. Emphasis is given to the assessment of the GNPs' properties and the physical mechanisms underlying the laser-induced heat generation in GNPs-loaded tissues. The main techniques available for temperature measurement and the current theoretical simulation approaches predicting the therapeutic outcome are reviewed. Topical challenges in delivering safe thermal dosage are also presented with the aim to discuss the state-of-the-art and the future perspective in the field of GNPs-mediated PTT.

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.

Numerical simulation of the effect of necrosis area in systemic delivery of magnetic nanoparticles in hyperthermia cancer treatment

In a magnetic hyperthermia treatment, malignant cancerous cells are ablated by the heat production of magnetic nanoparticles (MNP) under an external magnetic field. This novel approach is a promising tool to eliminate the tumor cells by a higher temperature inside the tumor microenvironment. MNPs are needed inside the tumor microenvironment to increase the heat, and this could be possible with intravenous drug injection. However, tumors with necrosis regions are more resistant to drug penetration, and this can cause inadequate and non-homogeneous temperature distribution in the tumor. Hence, in this study, we used numerical methods to investigate the Spatio-temporal temperature field distribution in the necrotic tumor and its surrounding tissue. To this end, an intravenous bolus injection is used to simulate the effect of systemic drug delivery in tumors with necrosis region. Results show that the temperature field with the necrosis region with 10% of the tumor radius is more prone to higher temperature values. The hypoxia region is affected by the high temperature despite the necrosis region in the tumor. However, a broader necrosis region impedes drug penetration inside the inner layers of tumors, which leads to a lower heat generation by the MNPs. Results also demonstrate that only 15.5% of MNP concentration distributed to the necrosis with 50% of tumor radius, leading a temperature of 42^∘C in the necrosis region, which is not sufficient for the tumor ablation. Therefore, the temperature distribution is dependant on the sizes of necrosis regions in tumors, and tumors with a larger necrotic region (over 20% of tumor radius) are challenging to treat with hyperthermia treatment. This study could help the future in vitro and in vivo studies of hyperthermia treatment in necrotic tumors.

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.

Design, Implementation, Evaluation and Application of a 32-Channel Radio Frequency Signal Generator for Thermal Magnetic Resonance Based Anti-Cancer Treatment

Thermal Magnetic Resonance (ThermalMR) leverages radio frequency (RF)-induced heating to examine the role of temperature in biological systems and disease. To advance RF heating with multi-channel RF antenna arrays and overcome the shortcomings of current RF signal sources, this work reports on a 32-channel modular signal generator (SG_PLL). The SG_PLL was designed around phase-locked loop (PLL) chips and a field-programmable gate array chip. To examine the system properties, switching/settling times, accuracy of RF power level and phase shifting were characterized. Electric field manipulation was successfully demonstrated in deionized water. RF heating was conducted in a phantom setup using self-grounded bow-tie RF antennae driven by the SG_PLL. Commercial signal generators limited to a lower number of RF channels were used for comparison. RF heating was evaluated with numerical temperature simulations and experimentally validated with MR thermometry. Numerical temperature simulations and heating experiments controlled by the SG_PLL revealed the same RF interference patterns. Upon RF heating similar temperature changes across the phantom were observed for the SG_PLL and for the commercial devices. To conclude, this work presents the first 32-channel modular signal source for RF heating. The large number of coherent RF channels, wide frequency range and accurate phase shift provided by the SG_PLL form a technological basis for ThermalMR controlled hyperthermia anti-cancer treatment.

Multipoint Temperature Monitoring of Microwave Thermal Ablation in Bones through Fiber Bragg Grating Sensor Arrays

Bones are a frequent site of metastases that cause intolerable cancer-related pain in 90% of patients, making their quality of life poor. In this scenario, being able to treat bone oncology patients by means of minimally invasive techniques can be crucial to avoid surgery-related risks and decrease hospitalization times. The use of microwave ablation (MWA) is gaining broad clinical acceptance to treat bone tumors. It is worth investigating temperature variations in bone tissue undergoing MWA because the clinical outcomes can be inferred from this parameter. Several feasibility studies have been performed, but an experimental analysis of the temperature trends reached into the bone during the MWA has not yet been assessed. In this work, a multi-point temperature study along the bone structure during such treatment is presented. The study has been carried out on ex vivo bovine femur and tibia, subjected to MWA. An overall of 40 measurement points covering a large sensing area was obtained for each configuration. Temperature monitoring was performed by using 40 fiber Bragg grating (FBGs) sensors (four arrays each housing 10 FBGs), inserted into the bones at specific distances to the microwave antenna. As result, the ability of this experimental multi-point monitoring approach in tracking temperature variations within bone tissue during MWA treatments was shown. This study lays the foundations for the design of a novel approach to study the effects of MWA on bone tumors. As consequence, the MWA treatment settings could be optimized in order to maximize the treatment effects of such a promising clinical application, but also customized for the specific tumor and patient.

A moderate thermal dose is sufficient for effective free and TSL based thermochemotherapy

Hyperthermia, i.e. heating the tumor to a temperature of 40-43 °C is considered by many a valuable treatment to sensitize tumor cells to radiotherapy and chemotherapy. In recent randomized trials the great potential of adding hyperthermia to chemotherapy was demonstrated for treatment of high risk soft tissue sarcoma: +11.4% 5 yrs. overall survival (OS) and for ovarian cancer with peritoneal involvement nearly +12 months OS gain. As a result interest in combining chemotherapy with hyperthermia, i.e. thermochemotherapy, is growing. Extensive biological research has revealed that hyperthermia causes multiple effects, from direct cell kill to improved oxygenation, whereby each effect has a specific temperature range. Thermal sensitization of the tumor cell for chemotherapy occurs for many drugs at temperatures ranging from 40 to 42 °C with little additional increase of sensitization at higher temperatures. Increasing perfusion/oxygenation and increased extravasation are two other important hyperthermia induced mechanisms. The combination of free drug and hyperthermia has not been found to increase tumor drug concentration. Hence, enhanced effectiveness of free drug will depend on the thermal sensitization of the tumor cells for the applied drug. In contrast to free drugs, experimental animal studies combining hyperthermia and thermo-sensitive liposomal (TSL) drugs delivery have demonstrated to result in a substantial increase of the drug concentration in the tumor. For TSL based chemotherapy, hyperthermia is critical to both increase perfusion and extravasation as well as to trigger TSL drug release, whereby the temperature controlled induction of a local high drug concentration in a highly permeable vessel is driving the enhanced drug uptake in the tumor. Increased drug concentrations up to 26 times have been reported in rodents. Good control of the tissue temperature is required to keep temperatures below 43 °C to prevent vascular stasis. Further, careful timing of the drug application relative to the start of heating is required to benefit optimal from the combined treatment. From the available experimental data it follows that irrespective whether chemotherapy is applied as free drug or using a thermal sensitive liposomal carrier, the optimal thermal dose for thermochemotherapy should be 40-42 °C for 30-60 min, i.e. equivalent to a CEM43 of 1-15 min. Timing is critical: most free drug should be applied simultaneous with heating, whereas TSL drugs should be applied 20-30 min after the start of hyperthermia.

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.

Recent technological advancements in radiofrequency- andmicrowave-mediated hyperthermia for enhancing drug delivery

Hyperthermia therapy is a potent enhancer of chemotherapy and radiotherapy. In particular, microwave (MW) and radiofrequency (RF) hyperthermia devices provide a variety of heating approaches that can treat most cancers regardless the size. This review introduces the physics of MW/RF hyperthermia, the current state-of-the-art systems for both localized and regional heating, and recent advancements in hyperthermia treatment guidance using real-time computational simulations and magnetic resonance thermometry. Clinical trials involving RF/MW hyperthermia as adjuvant for chemotherapy are also presented per anatomical site. These studies favor the use of adjuvant hyperthermia since it significantly improves curative and palliative clinical outcomes. The main challenge of hyperthermia is the distribution of state-of-the-art heating systems. Nevertheless, we anticipate that recent technology advances will expand the use of hyperthermia to chemotherapy centers for enhanced drug delivery. These new technologies hold great promise not only for (image-guided) perfusion modulation and sensitization for cytotoxic drugs, but also for local delivery of various compounds using thermosensitive liposomes.

Field drift correction of proton resonance frequency shift temperature mapping with multichannel fast alternating nonselective free induction decay readouts

PURPOSE

To demonstrate that proton resonance frequency shift MR thermometry (PRFS-MRT) acquisition with nonselective free induction decay (FID), combined with coil sensitivity profiles, allows spatially resolved B₀ drift-corrected thermometry.

METHODS

Phantom experiments were performed at 1.5T and 3T. Acquisition of PRFS-MRT and FID were performed during MR-guided high-intensity focused ultrasound heating. The phase of the FIDs was used to estimate the change in angular frequency δω_drift per coil element. Two correction methods were investigated: (1) using the average δω_drift over all coil elements (0th-order) and (2) using coil sensitivity profiles for spatially resolved correction. Optical probes were used for independent temperature verification. In-vivo feasibility of the methods was evaluated in the leg of 1 healthy volunteer at 1.5T.

RESULTS

In 30 minutes, B₀ drift led to an apparent temperature change of up to -18°C and -98°C at 1.5T and 3T, respectively. In the sonicated area, both corrections had a median error of 0.19°C at 1.5T and -0.54°C at 3T. At 1.5T, the measured median error with respect to the optical probe was -1.28°C with the 0th-order correction and improved to 0.43°C with the spatially resolved correction. In vivo, without correction the spatiotemporal median of the apparent temperature was at -4.3°C and interquartile range (IQR) of 9.31°C. The 0th-order correction had a median of 0.75°C and IQR of 0.96°C. The spatially resolved method had the lowest median at 0.33°C and IQR of 0.80°C.

CONCLUSION

FID phase information from individual receive coil elements allows spatially resolved B₀ drift correction in PRFS-based MRT.

Preliminary analysis of ultrasound elastography imaging-based thermometry on non-perfused ex vivo swine liver

AIMS

Real-time monitoring of tissue temperature during percutaneous tumor ablation improves treatment efficacy, leading clinicians in adjustment of treatment settings. This study aims at assessing feasibility of ultrasound thermometry during laser ablation of biological tissue using a specific ultrasound imaging techniques based on elastography acoustic radiation force impulse (ARFI).

METHODS

ARFI uses high-intensity focused ultrasound pulses to generate 'radiation force' in tissue; this provokes tissue displacements trackable using correlation-based ultrasound methods: the sensitivity of shear waves velocity is able to detect temperature changes. Experiments were carried out using a Nd:YAG laser (power: 5 W) in three non-perfused ex vivo pig livers. In each organ, a thermocouple was placed close to the applicator tip (distance range 1.5-2.5 cm) used to record a reference temperature. Positioning of laser applicator and thermocouple was eco-guided. The organ was scanned by an echography system equipped with ARFI; propagation velocity was measured in a region of interest of 1 × 0.5 cm located close to thermocouple, to investigate influence of tissue temperature on shear waves velocity.

RESULTS

Shear wave velocity has a very low sensitivity to temperature up to 55-60 °C, and in all cases, velocity is < 5 m s^-1; for temperature > 55-60 °C, velocity shows a steep increment. The system measures a value "over limit", meaning a velocity > 5 m s^-1.

CONCLUSIONS

Ultrasound thermometry during laser ablation of biological tissue based on elastography shows an abrupt output change at temperatures > 55-60 °C. This issue can have a relevant clinical impact, considering tumor necrosis when temperature crosses 55 °C to define the boundary of damaged volume.

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.

Predictive value of SAR based quality indicators for head and neck hyperthermia treatment quality

PURPOSE

Hyperthermia treatment quality determines treatment effectiveness as shown by the clinically derived thermal-dose effect relations. SAR based optimization factors are used as possible surrogate for temperature, since they are not affected by thermal tissue properties uncertainty and variations. Previously, target coverage (TC) at the 25% and 50% iso-SAR level was shown predictive for treatment outcome in superficial hyperthermia and the target-to-hot-spot-quotient (THQ) was shown to highly correlate with predictive temperature in deep pelvic hyperthermia. Here, we investigate the correlation with temperature for THQ and TC using an 'intermediate' scenario: semi-deep hyperthermia in the head & neck region using the HYPERcollar3D.

METHODS

Fifteen patient-specific models and two different planning approaches were used, including random perturbations to circumvent optimization bias. The predicted SAR indicators were compared to predicted target temperature distribution indicators T50 and T90, i.e., the median and 90th percentile temperature respectively.

RESULTS

The intra-patient analysis identified THQ, TC25 and TC50 as good temperature surrogates: with a mean correlation coefficient R²_T50 = 0.72 and R²_T90=0.66. The inter-patient analysis identified the highest correlation with TC25 (R²_T50 = 0.76, R²_T90=0.54) and TC50 (R²_T50 = 0.74, R²_T90 = 0.56).

CONCLUSION

Our investigation confirmed the validity of our current strategy for deep hyperthermia in the head & neck based on a combination of THQ and TC25. TC50 was identified as the best surrogate since it enables optimization and patient inclusion decision making using one single parameter.

The potential of constrained SAR focusing for hyperthermia treatment planning: analysis for the head & neck region

Clinical trials have shown that hyperthermia is a potent adjuvant to conventional cancer treatments, but the temperatures currently achieved in the clinic are still suboptimal. Hyperthermia treatment planning simulations have potential to improve the heating profile of phased-array applicators. An important open challenge is the development of an effective optimization procedure that enables uniform heating of the target region while keeping temperature below a threshold in healthy tissues. In this work, we analyzed the effectiveness and efficiency of a recently proposed optimization approach, i.e. focusing via constrained power optimization (FOCO), using 3D simulations of twelve clinical patient specific models. FOCO performance was compared against a clinically used particle swarm based optimization approach. Evaluation metrics were target coverage at the 25% iso-SAR level, target hotspot quotient, median target temperature (T50) and computational requirements. Our results show that, on average, constrained power focusing performs slightly better than the clinical benchmark ([Formula: see text]T50 [Formula: see text] °C), but outperforms this clinical benchmark for large target volumes ([Formula: see text]40 cm[Formula: see text], [Formula: see text]T50 [Formula: see text] °C). In addition, the results are achieved in a shorter time ([Formula: see text]%) and are repeatable because the approach is formulated as a convex optimization problem.

Functionalized Anatomical Models for Computational Life Sciences

The advent of detailed computational anatomical models has opened new avenues for computational life sciences (CLS). To date, static models representing the anatomical environment have been used in many applications but are insufficient when the dynamics of the body prevents separation of anatomical geometrical variability from physics and physiology. Obvious examples include the assessment of thermal risks in magnetic resonance imaging and planning for radiofrequency and acoustic cancer treatment, where posture and physiology-related changes in shape (e.g., breathing) or tissue behavior (e.g., thermoregulation) affect the impact. Advanced functionalized anatomical models can overcome these limitations and dramatically broaden the applicability of CLS in basic research, the development of novel devices/therapies, and the assessment of their safety and efficacy. Various forms of functionalization are discussed in this paper: (i) shape parametrization (e.g., heartbeat, population variability), (ii) physical property distributions (e.g., image-based inhomogeneity), (iii) physiological dynamics (e.g., tissue and organ behavior), and (iv) integration of simulation/measurement data (e.g., exposure conditions, “validation evidence” supporting model tuning and validation). Although current model functionalization may only represent a small part of the physiology, it already facilitates the next level of realism by (i) driving consistency among anatomy and different functionalization layers and highlighting dependencies, (ii) enabling third-party use of validated functionalization layers as established simulation tools, and (iii) therefore facilitating their application as building blocks in network or multi-scale computational models. Integration in functionalized anatomical models thus leverages and potentiates the value of sub-models and simulation/measurement data toward ever-increasing simulation realism. In our o²S²PARC platform, we propose to expand the concept of functionalized anatomical models to establish an integration and sharing service for heterogeneous computational models, ranging from the molecular to the organ level. The objective of o²S²PARC is to integrate all models developed within the National Institutes of Health SPARC initiative in a unified anatomical and computational environment, to study the role of the peripheral nervous system in controlling organ physiology. The functionalization concept, as outlined for the o²S²PARC platform, could form the basis for many other application areas of CLS. The relationship to other ongoing initiatives, such as the Physiome Project, is also presented.

Clinical validation of a novel thermophysical bladder model designed to improve the accuracy of hyperthermia treatment planning in the pelvic region

PURPOSE

Hyperthermia treatment planning for deep locoregional hyperthermia treatment may assist in phase and amplitude steering to optimize the temperature distribution. This study aims to incorporate a physically correct description of bladder properties in treatment planning, notably the presence of convection and absence of perfusion within the bladder lumen, and to assess accuracy and clinical implications for non muscle invasive bladder cancer patients treated with locoregional hyperthermia.

METHODS

We implemented a convective thermophysical fluid model based on the Boussinesq approximation to the Navier-Stokes equations using the (finite element) OpenFOAM toolkit. A clinician delineated the bladder on CT scans obtained from 14 bladder cancer patients. We performed (1) conventional treatment planning with a perfused muscle-like solid bladder, (2) with bladder content properties without and (3) with flow dynamics. Finally, we compared temperature distributions predicted by the three models with temperature measurements obtained during treatment.

RESULTS

Much higher and more uniform bladder temperatures are predicted with physically accurate fluid modeling compared to previously employed muscle-like models. The differences reflect the homogenizing effect of convection, and the absence of perfusion. Median steady state temperatures simulated with the novel convective model (3) deviated on average -0.6 °C (-12%) from values measured during treatment, compared to -3.7 °C (-71%) and +1.5 °C (+29%) deviation for the muscle-like (1) and static (2) models, respectively. The Grashof number was 3.2 ± 1.5 × 10⁵ (mean ± SD).

CONCLUSIONS

Incorporating fluid modeling in hyperthermia treatment planning yields significantly improved predictions of the temperature distribution in the bladder lumen during hyperthermia treatment.

Estimation of the temperature field in laser-induced hyperthermia experiments with a phantom

BACKGROUND

One of the challenges faced during the hyperthermia treatment of cancer is to monitor the temperature distribution in the region of interest. The main objective of this work was to accurately estimate the transient temperature distribution in the heated region, by using a stochastic heat transfer model and temperature measurements.

METHODS

Experiments involved the laser heating of a cylindrical phantom, partially loaded with iron oxide nanoparticles. The nanoparticles were manufactured and characterized in this work. The solution of the state estimation problem was obtained with an algorithm of the Particle Filter method, which allowed for simultaneous estimation of state variables and model parameters. Measurements of one single sensor were used for the estimation procedure, which is highly desirable for practical applications in order to avoid patient discomfort.

RESULTS

Despite the large uncertainties assumed for the model parameters and for the coupled radiation-conduction model, discrepancies between estimated temperatures and internal measurements were smaller than 0.7 °C. In addition, the estimated fluence rate distribution was physically meaningful. Maximum discrepancies between the prior means and the estimated means were of 2% for thermal conductivity and heat transfer coefficient, 4% for the volumetric heat capacity and 3% for the irradiance.

CONCLUSIONS

This article demonstrated that the Particle Filter method can be used to accurately predict the temperatures in regions where measurements are not available. The present technique has potential applications in hyperthermia treatments as an observer for active control strategies, as well as to plan personalized heating protocols.

Groupwise Multichannel Image Registration

Multichannel image registration is an important challenge in medical image analysis. Multichannel images result from modalities such as dual-energy CT or multispectral microscopy. Besides, multichannel feature images can be derived from acquired images, for instance, by applying multiscale feature banks to the original images to register. Multichannel registration techniques have been proposed, but most of them are applicable to only two multichannel images at a time. In the present study, we propose to formulate multichannel registration as a groupwise image registration problem. In this way, we derive a method that allows the registration of two or more multichannel images in a fully symmetric manner (i.e., all images play the same role in the registration procedure), and therefore, has transitive consistency by definition. The method that we introduce is applicable to any number of multichannel images, any number of channels per image, and it allows to take into account correlation between any pair of images and not just corresponding channels. In addition, it is fully modular in terms of dissimilarity measure, transformation model, regularisation method, and optimisation strategy. For two multimodal datasets, we computed feature images from the initially acquired images, and applied the proposed registration technique to the newly created sets of multichannel images. MIND descriptors were used as feature images, and we chose total correlation as groupwise dissimilarity measure. Results show that groupwise multichannel image registration is a competitive alternative to the pairwise multichannel scheme, in terms of registration accuracy and insensitivity towards registration reference spaces.

Magnetic Hyperthermia and Radiation Therapy: Radiobiological Principles and Current Practice †

Hyperthermia, though by itself generally non-curative for cancer, can significantly increase the efficacy of radiation therapy, as demonstrated by in vitro, in vivo, and clinical results. Its limited use in the clinic is mainly due to various practical implementation difficulties, the most important being how to adequately heat the tumor, especially deep-seated ones. In this work, we first review the effects of hyperthermia on tissue, the limitations of radiation therapy and the radiobiological rationale for combining the two treatment modalities. Subsequently, we review the theory and evidence for magnetic hyperthermia that is based on magnetic nanoparticles, its advantages compared with other methods of hyperthermia, and how it can be used to overcome the problems associated with traditional techniques of hyperthermia.