The interest in hybrid nanoparticles for various applications in biomedicine is growing inevitably, stimulating research towards not only more effective, but also more accessible systems. This resulted in the emergence of advanced synthetic protocols with optimized conditions for the production of nanoparticles with high yields and desired morphologies, which ultimately determine their physicochemical and biomedical properties. While these challenges were sufficient for scientists a few decades ago, the sustainability of the synthetic methods is now an important aspect. From this perspective, nanoparticle production methods based on physical principles, such as spark discharge phenomena, could provide an interesting alternative to labor-intensive and environmentally harmful chemical synthesis. The benefits of clean and sustainable physical production routes for various nanomaterials are already recognized in the fields of catalysis and electronics. Biomedicine on the other hand has been reluctant to embrace the new methodologies, as they do not inherently provide nanoparticles dispersed in aqueous media, which is essential for their safe administration and reliable physiological performance. In this work, we investigated the potential of spark discharge as an alternative method to produce hybrid palladium/iron oxide nanoparticles intended for cancer thermo-brachytherapy by leveraging the magnetothermal properties of iron and the favorable radioactive features of the palladium radioisotope. Focusing on the aqueous harvesting of the nanoparticles produced in VSParticle’s spark discharge generator, we determined the optimal settings compatible with the connected bubbling column and identified the pitfalls and possible solutions to the intrinsic challenges, such as low yields and aggregation.

Introduction: Hyperthermia efficacy depends on the temperatures achieved in the target area. Therefore, hyperthermia systems must deliver both controlled and conformal heating. This study presents a comprehensive multi-institutional quality assurance (QA) evaluation of deep hyperthermia devices. Methods: Six European institutions equipped with BSD- Sigma 60 and Sigma Eye deep hyperthermia applicators participated in the study. Up to six measurements per applicator were performed in each institution. The thermal distribution in cylindrical homogeneous phantoms after 10 minutes of heating with a total power delivered of 1000 watts was assessed using the applicator’s integrated mapping thermometry system. Evaluated quality parameters included temperature increase, focus location, and focus symmetry. Results: A total of 54 measurements were conducted, with 43 included in the analysis. All applicators, except one, achieved a temperature increase of 6 °C in 10 minutes. Central heating capabilities were demonstrated, with mean deviations from the intended location of −1.4 ± 1.6 cm for Sigma 60 and 1.5 ± 1.4 cm for Sigma Eye. Symmetry evaluations showed differences in radial temperature profiles of 6.2 ± 4.5 % for the Sigma 60 and 5.9 ± 4.4 % for the Sigma Eye. We propose minimum acceptable values for each quality parameter based on these results. Conclusion: The measurements were reproducible with acceptable values for the various quality parameters. Potential deviations might be attributed to inaccuracies in the mapping thermometry system rather than the heating system. The presented protocol and practical recommendations should be applied for future QA measurements in deep hyperthermia.

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.

Magnetic resonance thermometry (MRT) can measure in-vivo 3D-temperature changes in real-time and noninvasively. However, for the oropharynx region and the entire head and neck, motion potentially introduces large artifacts. Considering long treatment times of 60–90 min, this study aims to evaluate whether MRT around the oropharynx is clinically feasible for hyperthermia treatments and quantify the effects of breathing and swallowing on MRT performance. A 3D-ME-FGRE sequence was used in a phantom cooling down and around the oropharynx of five volunteers over ∼75 min. The imaging protocol consisted of imaging with acceleration (ARC = 2), number of image averages (NEX = 1,2 and 3). For volunteers, the acquisitions included a breath-hold scan and scans with deliberate swallowing. MRT performance was quantified in neck muscle, spinal cord and masseter muscle, using mean average error (MAE), mean error (ME) and spatial standard deviation (SD). In phantom, an increase in NEX leads to a significant decrease in SD, but MAE and ME were unchanged. No significant difference was found in volunteers between the different scans. There was a significant difference between the regions evaluated: neck muscle had the best MAE (=1.96 °C) and SD (=0.82 °C), followed by spinal cord (MAE = 3.17 °C, SD = 0.92 °C) and masseter muscle (MAE = 4.53 °C, SD = 1.16 °C). Concerning the ME, spinal cord did best, then neck muscle and masseter muscle, with values of −0.64 °C, 1.15 °C and −3.05 °C respectively. Breathing, swallowing, and different ways of imaging (acceleration and NEX) do not significantly influence the MRT performance in the oropharynx region. The ROI selected however, leads to significant differences.

While hyperthermia has been shown to induce a variety of cytotoxic and sensitizing effects on cancer tissues, the thermal dose–effect relationship is still not well quantified, and it is still unclear how it can be optimally combined with other treatment modalities. Additionally, it is speculated that different methods of applying hyperthermia, such as water bath heating or electromagnetic energy, may have an effect on the resulting biological mechanisms involved in cell death or in sensitizing tumor cells to other oncological treatments. In order to further quantify and characterize hyperthermia treatments on a cellular level, in vitro experiments shifted towards the use of 3D cell spheroids. These are in fact considered a more representative model of the cell environment when compared to 2D cell cultures. In order to perform radiofrequency (RF)-induced heating in vitro, we have recently developed a dedicated electromagnetic field applicator. In this study, using this applicator, we designed and validated an experimental setup which can heat 3D cell spheroids in a conical polypropylene vial, thus providing a reliable instrument for investigating hyperthermia effects at the cellular scale.

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.

Purpose: This study aimed to assess the quality of the lucite cone applicator (LCA), the standard applicator for superficial hyperthermia at the Erasmus MC Cancer Institute, using the most recent quality assurance guidelines, thus verifying their feasibility. Materials and methods: The assessment was conducted on each of the six LCAs available for clinical treatments. The temperature distribution was evaluated using an infrared camera across different layers of a fat-muscle mimicking phantom. The maximum temperature increase, thermal effective penetration depth (TEPD), and thermal effective field size (TEFS) were used as quality metrics. The experimental results were validated through comparison with simulated results, using a canonical phantom model and a realistic phantom model segmented from CT imaging. Results: A maximum temperature increase above 6 °C at 2 cm depth in the fat-muscle phantom for all the experiments was found. A mean negative difference between simulated and experimental data was of 1.3 °C when using the canonical phantom model. This value decreased to a mean negative difference of 0.4 °C when using the realistic model. Simulated and measured TEPD showed good agreement for both in silico scenarios, while discrepancies were present for TEFS. Conclusions: The LCAs passed all QA guidelines requirements for superficial hyperthermia delivery when used singularly or in an array configuration. A further characterization of parameters such as antenna efficiency and heat transfer coefficients would be beneficial for translating experimental results to simulated values. Implementing the QA guidelines was time-consuming and demanding, requiring careful preparation and correct setup of antenna elements.

Purpose: MR thermometry (MRT) enables noninvasive temperature monitoring during hyperthermia treatments. MRT is already clinically applied for hyperthermia treatments in the abdomen and extremities, and devices for the head are under development. In order to optimally exploit MRT in all anatomical regions, the best sequence setup and post-processing must be selected, and the accuracy needs to be demonstrated. Methods: MRT performance of the traditionally used double-echo gradient-echo sequence (DE-GRE, 2 echoes, 2D) was compared to multi-echo sequences: a 2D fast gradient-echo (ME-FGRE, 11 echoes) and a 3D fast gradient-echo sequence (3D-ME-FGRE, 11 echoes). The different methods were assessed on a 1.5 T MR scanner (GE Healthcare) using a phantom cooling down from 59 °C to 34 °C and unheated brains of 10 volunteers. In-plane motion of volunteers was compensated by rigid body image registration. For the ME sequences, the off-resonance frequency was calculated using a multi-peak fitting tool. To correct for B0 drift, the internal body fat was selected automatically using water/fat density maps. Results: The accuracy of the best performing 3D-ME-FGRE sequence was 0.20 °C in phantom (in the clinical temperature range) and 0.75 °C in volunteers, compared to DE-GRE values of 0.37 °C and 1.96 °C, respectively. Conclusion: For hyperthermia applications, where accuracy is more important than resolution or scan-time, the 3D-ME-FGRE sequence is deemed the most promising candidate. Beyond its convincing MRT performance, the ME nature enables automatic selection of internal body fat for B0 drift correction, an important feature for clinical application.

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.

Biological modeling for anti-cancer treatments using mathematical models can be very supportive in gaining more insight into dynamic processes responsible for cellular response to treatment, and predicting, evaluating and optimizing therapeutic effects of treatment. This review presents an overview of the current status of biological modeling for hyperthermia in combination with radiotherapy (thermoradiotherapy). Various distinct models have been proposed in the literature, with varying complexity; initially aiming to model the effect of hyperthermia alone, and later on to predict the effect of the combined thermoradiotherapy treatment. Most commonly used models are based on an extension of the linear-quadratic (LQ)-model enabling an easy translation to radiotherapy where the LQ model is widely used. Basic predictions of cell survival have further progressed toward 3 D equivalent dose predictions, i.e., the radiation dose that would be needed without hyperthermia to achieve the same biological effect as the combined thermoradiotherapy treatment. This approach, with the use of temperature-dependent model parameters, allows theoretical evaluation of the effectiveness of different treatment strategies in individual patients, as well as in patient cohorts. This review discusses the significant progress that has been made in biological modeling for hyperthermia combined with radiotherapy. In the future, when adequate temperature-dependent LQ-parameters will be available for a large number of tumor sites and normal tissues, biological modeling can be expected to be of great clinical importance to further optimize combined treatments, optimize clinical protocols and guide further clinical studies.

The challenge to explain the diffuse and unconclusive message reported by hyperthermia studies investigating the thermal dose parameter is still to be unravelled. In the present review, we investigated a wide range of technical and clinical parameters characterising hyperthermia treatment to better understand and improve the probability of detecting a thermal dose effect relationship in clinical studies. We performed a systematic literature review to obtain hyperthermia clinical studies investigating the associations of temperature and thermal dose parameters with treatment outcome or acute toxicity. Different hyperthermia characteristics were retrieved, and their influence on temperature and thermal dose parameters was assessed. In the literature, we found forty-eight articles investigating thermal dose effect relationships. These comprised a total of 4107 patients with different tumour pathologies. The association between thermal dose and treatment outcome was the investigated endpoint in 90% of the articles, while the correlation between thermal dose and toxicity was investigated in 50% of the articles. Significant associations between temperature-related parameters and treatment outcome were reported in 63% of the studies, while those between temperature-related parameters and toxicity were reported in 15% of the studies. One clear difficulty for advancement is that studies often omitted fundamental information regarding the clinical treatment, and among the different characteristics investigated, thermometry details were seldom and divergently reported. To overcome this, we propose a clear definition of the terms and characteristics that should be reported in clinical hyperthermia treatments. A consistent report of data will allow their use to further continue the quest for thermal dose effect relationships.

The evaluation of the biological effects of therapeutic hyperthermia in oncology and the precise quantification of thermal dose, when heating is coupled with radiotherapy or chemotherapy, are active fields of research. The reliable measurement of hyperthermia effects on cells and tissues requires a strong control of the delivered power and of the induced temperature rise. To this aim, we have developed a radiofrequency (RF) electromagnetic applicator operating at 434 MHz, specifically engineered for in vitro tests on 3D cell cultures. The applicator has been designed with the aid of an extensive modelling analysis, which combines electromagnetic and thermal simulations. The heating performance of the built prototype has been validated by means of temperature measurements carried out on tissue-mimicking phantoms and aimed at monitoring both spatial and temporal temperature variations. The experimental results demonstrate the capability of the RF applicator to produce a well-focused heating, with the possibility of modulating the duration of the heating transient and controlling the temperature rise in a specific target region, by simply tuning the effectively supplied power.

In high-dose-rate brachytherapy (HDR-BT) for prostate cancer treatment, interstitial hyperthermia (IHT) is applied to sensitize the tumor to the radiation (RT) dose, aiming at a more efficient treatment. Simultaneous application of HDR-BT and IHT is anticipated to provide maximum radiosensitization of the tumor. With this rationale, the ThermoBrachyTherapy applicators have been designed and developed, enabling simultaneous irradiation and heating. In this research, we present a method to optimize the three-dimensional temperature distribution for simultaneous HDR-BT and IHT based on the resulting equivalent physical dose (EQDphys) of the combined treatment. First, the temperature resulting from each electrode is precomputed. Then, for a given set of electrode settings and a precomputed radiation dose, the EQDphys is calculated based on the temperature-dependent linear-quadratic model. Finally, the optimum set of electrode settings is found through an optimization algorithm. The method is applied on implant geometries and anatomical data of 10 previously irradiated patients, using reported thermoradiobiological parameters and physical doses. We found that an equal equivalent dose coverage of the target can be achieved with a physical RT dose reduction of 20% together with a significantly lower EQDphys to the organs at risk (p-value < 0.001), even in the least favorable scenarios. As a result, simultaneous ThermoBrachy-Therapy could lead to a relevant therapeutic benefit for patients with prostate cancer.

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

Objective: The addition of hyperthermia in the treatment of intact breast cancer with the aim to improve local response is currently in a research phase. First, optimal hyperthermia devices need to be developed, for which a diverse, anatomically and pathologically accurate set of patient models is necessary. Methods: To investigate the effects of inter-subject variations on hyperthermia treatment plans, we generated a repository of 22 anatomically and pathologically diverse patient models based on MR images of breast cancer patients. Hyperthermia treatment plans were generated for the 22 models using a generic theoretical phased array hyperthermia applicator. Results: Good temperature coverage was achieved in the vast majority of the models, with median values for T ₁₀ = 43.5°C (41.9–43.8°C), T ₅₀ = 42.5°C (41.3–43.3°C), and T ₉₀ = 41.3°C (39.8–42.6°C) under the condition that the maximum temperature increase in the patient is limited to 44°C. Conclusions: For future development of hyperthermia devices and treatment methods, a repository with a sufficiently large number of representative patient models, such as the one provided in this study, should be used to ensure applicability to a wide variety of patients. This repository is therefore made publicly available.

Temperature monitoring plays a central role in improving clinical effectiveness of adjuvant hyperthermia. The potential of magnetic resonance thermometry for treatment monitoring purposes led to several MR-guided hyperthermia approaches. However, the proposed solutions were sub-optimal due to technological and intrinsic limitations. These hamper achieving target conformal heating possibilities (applicator limitations) and accurate thermometry (inadequate signal-to-noise-ratio (SNR)). In this work, we studied proof of principle of a dual-function hyperthermia approach based on a coil array (64 MHz, 1.5 T) that is integrated in-between a phased array for heating (434 MHz) for maximum signal receive in order to improve thermometry accuracy. Hereto, we designed and fabricated a superficial hyperthermia mimicking planar array setup to study the most challenging interactions of generic phased-array setups in order to validate the integrated approach. Experiments demonstrated that the setup complies with the superficial hyperthermia guidelines for heating and is able to improve SNR at 2-4 cm depth by 17%, as compared to imaging using the body coil. Hence, the results showed the feasibility of our dual-function MR-guided hyperthermia approach as basis for the development of application specific setups.