Next Article in Journal
Erector Spinae Plane Block versus Paravertebral Block after Thoracic Surgery for Lung Cancer: A Propensity Score Study
Next Article in Special Issue
Holographic Focused Ultrasound Hyperthermia System for Uniform Simultaneous Thermal Exposure of Multiple Tumor Spheroids
Previous Article in Journal
Evaluation of Response to Atezolizumab Plus Bevacizumab in Patients with Advanced Hepatocellular Carcinoma Using the Combination of Response Evaluation Criteria in Solid Tumors and Alpha-Fetoprotein
Previous Article in Special Issue
Multimodal Prehabilitation in Patients Undergoing Complex Colorectal Surgery, Liver Resection, and Hyperthermic Intraperitoneal Chemotherapy (HIPEC): A Pilot Study on Feasibility and Potential Efficacy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Max-Delbrück-Center for Molecular Medicine in the Helmholtz Association (MDC), Berlin Ultrahigh Field Facility (B.U.F.F.), 13125 Berlin, Germany
2
Charité—Universitätsmedizin Berlin, Experimental and Clinical Research Center (ECRC), A Joint Cooperation between the Charité Medical Faculty and the Max-Delbrück Center for Molecular Medicine in the Helmholtz Association, 13125 Berlin, Germany
3
MRI.TOOLS GmbH, 13125 Berlin, Germany
4
Brightmind.AI GmbH, 1010 Vienna, Austria
5
Hasso Plattner Institute for Digital Engineering, University of Potsdam, 14482 Potsdam, Germany
*
Authors to whom correspondence should be addressed.
Cancers 2023, 15(8), 2303; https://doi.org/10.3390/cancers15082303
Submission received: 6 March 2023 / Revised: 9 April 2023 / Accepted: 11 April 2023 / Published: 14 April 2023
(This article belongs to the Collection Hyperthermia in Cancer Therapy)

Abstract

:

Simple Summary

Localized thermal therapy has been reported to have clinical benefits as a potent sensitizer of chemo and radiotherapy for various cancers, and for facilitating targeted drug delivery. Thermal magnetic resonance (ThermalMR) integrates targeted radiofrequency (RF) induced heating in the hyperthermia (HT) range, together with diagnostic MRI and in vivo non-invasive temperature mapping within a single RF applicator for thermal theranostics. The potential of dipole antenna arrays for ThermalMR and hyperthermia RF applicators is well recognized, but the additional value of loop elements remains to be investigated. Therefore, we designed circular and elliptical ThermalMR RF applicator arrays with circular 360° and horse-shoe shaped (arc = 270°) coverage of the human head, combining loop antennas and self-grounded bow-tie (SGBT) dipole antennas in a hybrid design. We investigated the performance of these designs in terms of the MRI transmission B1+ field, RF power deposition to optimize targeted RF heating, using electromagnetic field (EMF) and temperature simulations performed on a virtual patient with a clinically realistic intracranial brain tumor. ThermalMR RF applicators with the hybrid loop+SGBT dipole design showed superior MRI performance and targeted RF heating inside the tumor, while preserving healthy tissue.

Abstract

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.

Graphical Abstract

1. Introduction

Temperature is a physical parameter with a multitude of biological implications and a subject of great clinical interest. Fighting fire with fire, therapeutic hyperthermia (HT) (T = 40–43 °C) is an adjunct treatment of cancer therapy that can inhibit tumor growth and enhance the efficacy of other anti-cancer treatments [1,2,3,4,5,6,7,8]. Localized thermal therapy as a potent sensitizer of chemo-, radio- or immunotherapy of cancers has been reported to show clinical benefits [1,2,3,4,5,6,7,8,9,10]. Technical developments and research directions of thermal therapy have followed several trajectories, including magnetic nanoparticle hyperthermia (MNH), focused ultrasound-based hyperthermia, microwave-induced hyperthermia and radio frequency (RF)-driven hyperthermia [10,11,12,13,14,15,16,17,18,19,20,21]. With an ever-increasing clinical interest in thermal therapy, non-invasive in vivo approaches that facilitate diagnostic imaging-guided temperature therapy, along with the characterization of its subsequent effects, are imperative. An integrated platform that combines diagnosis and thermal therapy can help better define the role of temperature in biological systems and disease, and use this insight for improved thermal theranostics. There are several cutting-edge applications of thermal theranostics, including boosting the efficacy of immune-cell therapy and increasing permeability of the blood-brain barrier to allow crucial drugs to better enter the central nervous system and monitor the distribution of molecules in the interstitium [9,22,23,24,25,26]. Thermal therapy can support individualized medicine with temperature-triggered targeted drug delivery, using thermo-responsive nano-carriers for ‘smart’ theranostics [25,27,28,29]. Molecular therapy with heat-activated gene expression is a further potential application [30]. Using thermal modulation for selectively controlling neural activity patterns or genetically sensitized neurons (thermogenetics) is another compelling possibility [31].
Glioblastoma multiforme (GBM) represents 77–81% of most aggressive primary malignant tumors of the central nervous system (CNS) classified by the World Health Organization (WHO), with a median survival rate of approximately 14.6 months and no long-term cure [26,32,33,34,35,36]. Clinical studies suggest that patient survival improves when HT therapy is used as an adjunct to improve the efficacy of radio- and chemotherapy [7,8]. Adjuvant interstitial brain HT applied before and after brachytherapy was reported to significantly improve GBM patient survival (2-year survival 31% vs. 15%) with acceptable toxicity [7]. Notwithstanding the body of evidence on the effectiveness of radiotherapy + HT, thermal therapy is still not clinically widespread [1,37]. Achieving well-controlled tumor temperatures with optimal timing, sequence and dosage management constitutes major challenges of current HT applications, which limit realization of the full clinical potential of HT, especially for deep-seated tumors in the brain [1,17,38]. Current standalone RF HT approaches are constrained by limits in depth penetration, which preclude targeting deep tissue, and may be restricted to regional rather than focal heating [1,17,38]. Moreover, standalone HT approaches lack inherent in vivo non-invasive temperature mapping, real-time dose management and therapy detection capabilities [1,17,38].
Any in vivo HT modality strongly benefits from guidance by diagnostic imaging. Magnetic resonance (MR) is a mainstay of diagnostic imaging and an essential tool in fundamental biomedical research and clinical science. MR imaging (MRI) provides exquisite anatomical reference, facilitates functional contrast and supports non-invasive temperature mapping and therapy detection [39,40,41]. Multi-modal hybrid setups that combine MRI instruments and separate HT devices for targeted RF-induced heating using two independent RF chains have been established [42,43]. This approach requires extra hardware, is not cost effective and comes at the cost of mutual interferences between the two separate RF chains. This approach also suffers from compatibility issues and practical obstacles because each sub-device is provided by a different vendor. Thermal magnetic resonance (ThermalMR) integrates RF-induced heating in the HT range, in vivo non-invasive temperature mapping using MR thermometry (MRTh), and anatomical and functional MRI (and even physiometabolic MRI) in a single, multi-purpose RF applicator permitting supervised thermal theranostics [37,41,44,45,46,47,48,49,50,51,52,53,54]. ThermalMR adds a therapeutic dimension to a diagnostic imaging device and has great promise to address the urgent need for novel personalized patient care.
For ThermalMR, RF antenna arrays are used to non-invasively and selectively increase the temperature of a target region [37,44,46]. ThermalMR RF applicator development and optimization is governed by the operating frequency, type and size of RF antenna, and number of antenna elements. It is also guided by the spatial arrangement of the antennas in a phased array, with the objective to ensure uniform magnetic transmission fields (B1+) for diagnostic MRI and MR thermometry, and facilitate targeted control of electric fields (E) for thermal therapy [45,46]. MRI at ultrahigh magnetic field strengths (UHF-MRI) with B₀ ≥ 7.0 T uses higher RF frequencies than conventional MRI. This permits shorter wavelengths in tissue and offers ample potential for ThermalMR-based theranostics [45]. The potential of dipole antenna arrays for ThermalMR RF applicators, as well as microwave (MW) HT applicators, is well recognized [37,44,45,46,55,56,57,58,59,60]. Trefná et al. recently documented their work on SGBT dipole antenna arrangement over the head for MW-HT [61]. Improved MRI performance of loop antennas combined with dipole antennas was shown for UHF-MRI [62,63,64,65,66,67,68,69,70]. Until now, the potential of the additional value of integrating loop elements with dipole antennas in a ThermalMR RF applicator remained to be investigated. This approach is conceptually appealing because it supports the development of high-density arrays. This could benefit ThermalMR theranostics, because an increase in the number of RF elements per unit area increases the degrees of freedom for improving transmission field (B1+) uniformity for diagnostic MRI, and enhancing targeted RF shimming and RF peak power deposition for thermal therapy. This is especially relevant for ThermalMR of GBM or other deep-seated brain tumors, since the head presents a small surface area that limits the number of RF transmission elements that can be arranged in an array to cover the head.
Recognizing the opportunities of ThermalMR theranostics of brain tumors, we hypothesize that advanced RF applicators comprising combined loop+dipole elements will outperform the dipole-only counterparts. To test this hypothesis, we examine the suitability of hybrid loop+dipole RF applicator concepts for ThermalMR at 300 MHz (7.0 T), 400 MHz (9.4 T) and 450 MHz (10.5 T). For this purpose, we use a compact self-grounded bow-tie (SGBT) dipole antenna integrated with a loop element to form a loop+dipole building block. Eight hybrid loop+SGBT dipole building blocks are arranged in arrays to facilitate MRI and constructive E-field focusing in the target region of the brain. ThermalMR RF applicators with different spatial arrangements of the loop+SGBT building blocks are examined including circular and elliptical array configurations using 360° or horse-shoe-like 270° coverage of the brain. We hypothesize that the horse-shoe configurations are superior to the 360°-coverage configurations, and would improve patient safety and comfort by eliminating RF building blocks close to the nose, chin and eyes. To assess the performance of the ThermalMR RF applicators in a realistic setup, we performed electromagnetic field (EMF) and temperature simulations in a human voxel model, modified by including a small brain tumor replicating a clinical scenario.

2. Materials and Methods

Numerical EMF simulations using optimization algorithms in virtual patient models provide a springboard and technical foundation for the development of advanced ThermalMR RF applicator concepts customized for simultaneous MRI and RF hyperthermia for ThermalMR theranostics [71]. Realistic virtual patient models incorporated into the EMF simulations provide essential information about the performance of ThermalMR RF applicators with respect to the MRI characteristics and RF power deposition, i.e., specific absorption rate (SAR) and temperature distributions in the target regions of the brain. The workflow of ThermalMR theranostics including diagnostic imaging, treatment planning (TMTP), thermal therapy and therapy detection is presented in Figure 1.

2.1. Electromagnetic Field Simulations and Patient Model

The human voxel model ‘Duke’ from the virtual family (IT’IS Foundation Zürich, Switzerland) was used for EMF simulations and TMTP [72]. TMTP is based on EMF results simulated with CST Microwave Studio Suite 2020 (Dassault Systèmes Darmstadt, Germany) using the time-domain solver based on the finite integration technique (FIT). To reduce the computational effort of the EMF simulations, the human voxel model (maximum resolution of 2.0 × 2.0 × 2.0 mm3) was truncated at the level of the neck (Figure 2A).
To mimic a realistic clinical scenario, the Duke model was modified with an intracranial sphere (radius = 2 cm) that represents a small tumor in the right parietal region of the brain, with a target volume (TV) of 33.5 mL [37]. The dielectric and thermal properties of the tumor model were assigned based on the literature reports on human glioma where tumor perfusion is decreased by a factor of 0.7 [37,73,74]. The dielectric and thermal tissue properties of the remaining healthy tissues were assigned according to the database provided by the IT’IS Foundation [75]. To accelerate the simulation speed of broadband simulation by not modeling materials as dispersive, tissue dielectric properties were taken from the IT’IS database at 350 MHz as tissue dielectric properties do not vary strongly in the frequency range of interest. The dielectric and thermal properties of the tumor and healthy tissues used for the human voxel model Duke are summarized in Supplementary Table S1.

2.2. Integration of Loop and Dipole Antenna in a Hybrid RF Building Block

Loop antennas are single-frequency resonant, show asymmetric B1+ field distribution and provide good surface E-fields [62,63,64,65,66]. Dipole antenna arrays tailored for MRI offer uniform B1+ fields, can operate over a broadband frequency range and provide a central E-field distribution compared to loop arrays [62,63,64,65,66,67]. The current distribution of a dipole antenna is symmetric along its long axis, whereas a loop antenna demonstrates an antisymmetric current distribution [64,65,66]. Given the distinct current distribution patterns of loops and dipoles, mutual element decoupling can be achieved with loop+dipole combinations [62,63,64,65,66]. This provides an opportunity to develop advanced ThermalMR RF applicators.
Recognizing the opportunity, we implemented combined loop+dipole RF building blocks in EMF simulations. For this purpose, we used an ultra-wideband (range 250–650 MHz) self-grounded bow-tie (SGBT) dipole antenna (size: 42.3 × 46.3 × 2.5 mm3) adapted for ThermalMR (Figure 2B,C) [76]. The SGBT antennas were fed at the center. Each SGBT antenna uses a water bolus placed between the radiating element and the surface of the object under investigation to enhance the efficiency and directivity for targeted RF heating [77,78,79]. The water bolus was modeled as the shape of the outer surface of the SGBT antenna resonator and positioned between the SGBT and Duke’s head. To form a combined loop+SGBT dipole building block, the SGBT antenna was integrated with a rectangular loop antenna, illustrated in Figure 2B,C. Rectangular loop elements were modelled using 5 mm wide copper strips (size: 75 × 125 × 1 mm3). The loops were fed with a capacitive matching (Cm) network and decoupled from each other with a transformer decoupling Ld (Figure 2B,C) [80]. Eight fixed capacitors were distributed over the length of the loop to serve as tuning capacitors (Ct). The loops and SGBT dipoles were individually tuned and matched to support RF-induced target heating at multiple discrete frequencies: 250, 300, 350, 400, 450 MHz. The loop+SGBT building block will be further referred to as LD.
Figure 2. (A) Human voxel model ‘Duke’ from the virtual family IT’IS modified with a small brain tumor (red) in the right parietal region with coronal, sagittal, axial views [72]. (B) Compact SGBT dipole and loop building block. The wideband self-grounded bow-tie (SGBT) dipole antenna has a range of 250–600 MHz operating frequencies. (C) Schematic of loop+SGBT building blocks configured in an array. Loops were tuned to be multi-resonant to work with wideband SGBT dipoles at frequencies of 250, 300, 350, 400, 450 MHz. The loops are decoupled from each other with a transformer decoupling Ld [80]. Cm and Ct are matching and tuning capacitors of loops.
Figure 2. (A) Human voxel model ‘Duke’ from the virtual family IT’IS modified with a small brain tumor (red) in the right parietal region with coronal, sagittal, axial views [72]. (B) Compact SGBT dipole and loop building block. The wideband self-grounded bow-tie (SGBT) dipole antenna has a range of 250–600 MHz operating frequencies. (C) Schematic of loop+SGBT building blocks configured in an array. Loops were tuned to be multi-resonant to work with wideband SGBT dipoles at frequencies of 250, 300, 350, 400, 450 MHz. The loops are decoupled from each other with a transformer decoupling Ld [80]. Cm and Ct are matching and tuning capacitors of loops.
Cancers 15 02303 g002

2.3. ThermalMR RF Applicators

ThermalMR RF applicators were designed as an array of eight LD building blocks (hereafter referred to as ‘16LD’: 8 loops and 8 dipoles) positioned around the head of the human voxel model Duke (Figure 3A–D). The 16LD array is fed by 16 separate RF transmission channels. Two 16LD arrays with 360° coverage around the head were implemented: circular (cir, radius = 130 mm) and elliptical (ellip, radiusmajor axis = 130 mm, radiusminor axis = 100 mm) configurations (Figure 3A–C). Alternatively, two horse-shoe shaped circular (cir_HS) and elliptical (ellip_HS) arrays were designed, in an arc of 270° to ensure ample brain coverage while sparing the high conductivity regions of the eyes from RF exposure during targeted heating (Figure 3B–D).
All ThermalMR RF applicators support MRI at 300 MHz (B₀ = 7.0 T), 400 MHz (B₀ = 9.4 T) and 450 MHz (B₀ = 10.5 T) and RF-induced targeted heating using multiple discrete frequencies of 250, 300, 350, 400, 450 MHz. For best MRI performance, the magnetic transmission field B1+ should be perpendicular to the static magnetic field (B₀). Therefore, the LD building blocks were arranged with their long axis parallel to the head-feet (cranial-caudal) direction. Each LD was positioned around the head of the human voxel model so that the brain is centered in the x–y-plane (left–right, anterior–posterior) of the array and the TV is centered in the z-direction (superior–inferior) of the array. For comparison, ThermalMR RF applicator arrays with identical geometry and configuration to that of the 16LD, but comprising only eight SGBT dipole antennas (referred to as 8D) or only eight loop elements (referred to as 8L) were also implemented and evaluated in the EMF simulations (Figure 3A–D).

2.4. Transmission Field and RF Power Deposition for MRI

To assess the MRI capability of the designed ThermalMR RF applicators, we examined the B1+ distribution and RF power deposition (SAR10g). Post-processing of B1+ and RF power deposition was performed in MATLAB 2020 (The Mathworks, Natick, CA, USA). The B1+ distribution was calculated for the circularly polarized (CP) mode at magnetic field strengths of B₀ = 7.0 T, 9.4 T and 10.5 T, for 1 W forward power [81]. In the CP mode, each channel of the RF applicator is driven by a phase corresponding to its angular position relative to the center of the head in the transverse plane [82]. The transmission field (B1+) efficiency (given in µT/√kW) and SAR10g were compared for all ThermalMR RF applicators within a spherical region of interest (ROI, radius = 2 cm) centered in the human voxel model Duke’s head, for all three magnetic field strengths. RF power deposition was measured based on the specific absorption rate normalized to 1 W input power and averaged over 10 g of tissue (SAR10g). Safety guidelines governed by the IEC standard 60601-2-33 restrict the maximum SAR10g in 1st level operating mode to 3.2 W/kg (volume coil) and 20 W/kg (local transmit coil) for the whole head [83].

2.5. Targeted RF Heating for Thermal Therapy

Given the great complexity of thermoregulation in vivo, assessment of SAR10g as a common metric of local RF power deposition is essential before proceeding with any thermal intervention [71,73,84]. We used a time and frequency multiplexed vector field shaping (MVFS) algorithm to optimize local RF power deposition for targeted RF heating (Figure 4) [84]. The optimization algorithm provided globally optimal excitation vectors, by defining the phase and amplitude settings for each RF channel of the ThermalMR RF applicator. The MVFS algorithm automatically selects the appropriate intervention frequencies and time-interleaved excitations to best deliver RF power at the desired target shape and location. The resulting SAR distribution of the incident electric fields interference is tailored to focus heating of the TV while minimizing local peak RF exposure to healthy and remote tissue below a defined threshold (constraint SAR). Excitation vectors for each frequency are termed as the excitation mode (M). If M individual modes belong to the same excitation frequency, then time multiplexing is done with each solution vector scaled by √M, and the excitations were played out sequentially. Excitations at different frequencies can be played out concurrently as their electromagnetic fields do not interact coherently. The final resultant target shape is created by superimposing the individual patterns for each frequency.
Multiple discrete frequencies of 250, 300, 350, 400, 450 MHz were used to optimize RF-targeted heating using time and frequency multiplexing. The desired target SAR was set to 100 W/kg to maximize RF power deposition inside the tumor TV. SAR was constrained for regions outside the TV, i.e., the safe limit for healthy tissue: SAR10g,max(healthy tissue) = SARconstraint = 40 W/kg.
C ρ T t = · k T C b W b T T a + P L D
The SARconstraint limit was decided based on simplifying the temperature rise models using the Pennes’ bioheat transfer equation [85]. In Equation (1), Cp is the specific heat capacity (J/kg/°C); ρ is the tissue density (kg/m3); k is the thermal conductivity (W/m/°C); Cb is the specific heat capacity of blood (J/kg/°C); Wb is the volumetric perfusion rate (kg/m3/s); Ta is the local arterial or body core temperature; PLD is power loss density (W/m3) deposited by the heating system.
In steady-state formulation of the Pennes’ bioheat Equation (1), the time derivative on the left-hand side can be set to zero. The first transient term on the right-hand side can be neglected, as the characteristic time of temperature adaptation is short compared to the duration of the thermal therapy. Absent any time or temperature-dependent properties, the remaining terms (T − Ta) govern the temperature (T) increase. With a basal blood perfusion (Wb) of 50 mL/100 g/min, an average SAR of 50 W/kg results in a temperature change (T − Ta) of ~∆T ± 1.5 °C [37] without considering the response of a thermoregulatory perfusion (Cb) increase. Therefore, we set the SARconstraint at 40 W/kg to be in accordance with previous studies that used a SAR of approx. 20–40 W/kg, and reported a temperature increase to ~42 °C in the target region [86,87].

2.6. Temperature Simulations

Temperature distribution for 16LD ThermalMR RF applicators were calculated using the thermal transient solver of CST Microwave Studio Suite 2020, (Dassault Systèmes, Darmstadt, Germany). At first, the steady state temperature without any RF heating, i.e., basal temperature level or initial body temperature, was calculated using the thermal steady state solver at an ambient temperature T = 25 °C [87,88,89,90]. The initial temperature of the water boluses was fixed at T = 25 °C. Subsequently, the transient heating phase was evaluated for an intervention period of 30 min using the thermal transient solver considering isothermal boundaries with the goal to evaluate the relation between total achieved SAR10g and overall temperature rise in the TV from the basal body temperature level [87,88,89,90]. To simulate the effect of RF heating, the temperature was calculated from the power loss density W/m3 (local SAR10g) based on the phase and amplitude setting provided by the MVFS algorithm of the targeted RF heating optimization. The total power applied was the sum of the individual excitation mode power obtained from the targeted RF heating optimization for all 16LD ThermalMR RF applicators. The total power used in the temperature simulation is summarized in Table 1.

2.7. Thermal Therapy Quality Assessment

We evaluated the performance of the ThermalMR RF applicators using clinically relevant quality indicators provided by the European Society for Hyperthermic Oncology (ESHO) according to the criteria of characterization of the SAR10g, index temperature distribution inside TV and hotspots in healthy tissue [73,91,92,93,94,95,96]. These quality indicators allow the investigation of localized and averaged RF power deposition in the target volume, as well as in remote healthy tissue.
SARTV: SARTV (mean SAR10g within the TV) has been shown to be predictive of therapeutic outcomes of thermal treatment [92,93].
SARmax(TV): SARmax(TV) reflects the maximum absolute RF power deposition, and is directly related to the quality of the heating intervention [92,93].
Tumor coverage (TCx): Tumor coverage (TCx), or iso-SAR target coverage, is defined by the following Equation (2):
T C x = ( V T V ) x % S A R m a x V T V
representing the fraction (percentage) of the target volume (VTV) enclosed by x% peak SARmax isolines, i.e., TC25, TC50, TC80 and TC100 indicate the fraction of the tumor volume enclosed within the 25, 50, 80 and 100% isolines of peak SARmax, respectively [73,91,93,94]. Tumor coverage reflects the homogeneity of RF power deposition inside a tumor, and TC25 > 75% is typically considered for clinical treatment [73,91,93,94].
Hotspots-to-target quotient (HTQ): The hotspots-to-target quotient (HTQ) is the ratio of average SAR10g in the first percentile (P1) of healthy voxels exposed to the highest SAR10g versus the average SAR10g in the TV:
H T Q = P 1 , m e a n S A R 10 g H e a l t h y m e a n S A R 10 g T V
The HTQ reflects hotspot sizes, focusing on local maxima creation relative to the TV [37,73,91,93,94]. Values of HTQ ≤ 1 are typically considered acceptable for clinical treatment [73].
SAR amplification factor (SAF): The SAR amplification factor is the ratio of average local SAR10g in the TV to the average local SAR10g of healthy tissue [37,92,93,94]:
S A F = m e a n S A R 10 g T V m e a n S A R 10 g H e a l t h y
A higher SAF indicates less unwanted exposure of healthy tissue to RF power deposition, but does not provide information about specific local maxima or absolute RF power deposition levels [37]. In clinical routines, a SAF > 1 is considered desirable [92,93].
Tmean, T40°C, T41°C, T42°C: Tmean is the mean temperature in the tumor TV, and the cumulative minutes (cumin) for Tmean > TBasal (the basal body temperature) is an indicator of the heating efficiency [87,93,94]. Index temperature coverage (TX°C) is evaluated as the percent of the tumor volume covered with 40 °C (T40°C), 41 °C (T41°C), 42 °C (T42°C) [87,91,93].

2.8. Statistical Analysis

Statistical analysis was done for the MRI transmission field B1+ and optimized targeted RF heating results inside the tumor TV for all ThermalMR RF applicators. Data were evaluated for Gaussian distribution using the Shapiro-Wilk test. Differences in the metrics among the ThermalMR RF applicators were analyzed using the non-parametric Kruskal-Wallis ANOVA, followed by the Dunn’s post-hoc test with correction for multiple comparisons. Data were analyzed using the statistical environment R v.3.6.3; p-values < 0.05 were considered significant.

3. Results

3.1. Transmission Field and RF Power Deposition for MRI

The B1+ transmission field distribution is shown for a central-axial slice through the head of the human voxel model Duke for all ThermalMR RF applicators at B₀ = 7.0 T (Figure 5A), 9.4 T (Figure 5B) and 10.5 T (Figure 5C). For all arrays, the 16LD ThermalMR RF applicators had higher B1+ efficiency than the 8D and 8L arrays, indicating superior MRI performance.
At 7.0 T the 16LDcir, 16LDellip, 16LDcir_HS and 16LDellip_HS variants showed an improvement in mean B1+ efficiency of 51, 28, 53 and 30% over the 8D variants, and 58, 58, 31 and 54% over the 8L variants. The 16LDcir and 16LDellip variants yielded improved B1+ homogeneity in terms of maximum and mean B1+ versus the 16LDcir_HS and 16LDellip_HS counterparts. The 16LDellip variant provided the best overall maximum B1+ efficiency. At 9.4 T, the 16LDcir, 16LDellip, 16LDcir_HS and 16LDellip_HS variants showed an improvement in mean B1+ efficiency of 29, 20, 25 and 20% over the 8D variants, and 46, 71, 61, and 72% increase over 8L variants. At 10.5 T the 16LDcir, 16LDellip, 16LDcir_HS and 16LDellip_HS variants yielded an improvement in mean B1+ efficiency of 22%, 19%, 19% and 19% over the 8D variants, and 35%, 74%, 73% and 74% over the 8L variants. Mean, maximum and minimum values for each ThermalMR RF applicator are listed in Table 2.
Analysis of the B1+ efficiency achieved inside the ROI shows that the 16LD had statistically significant higher B1+ efficiency compared to the 8D and 8L variants, for all configurations at B₀ = 7.0 T (Figure 5D), 9.4 T (Figure 5E) and 10.5 T (Figure 5F) (p-values are listed in Table 2).
Aside from the intentional power deposition desired for targeted RF heating, the suitability of the ThermalMR applicators for MRI requires control of undesired power deposition. The SAR10g deposition maps obtained for a central-axial slice of Duke’s head for all ThermalMR RF applicators at B₀ = 7.0 T (Figure 6A,D), 9.4 T (Figure 6B,E) and 10.5 T (Figure 6C,F) show that in every case the SAR10g is well within the safe limits defined by the IEC standard. The 8Lcir_HS variant showed the highest maximum local SAR10g = 1.8 W/kg.
To summarize, the EMF simulations show that all 16LD variants provide MRI transmission fields and a safe limit of RF power deposition (SAR10g), which are very well suited for MRI at 7.0 T 9.4 T and 10.5 T.

3.2. Targeted RF Heating for Thermal Therapy

Results of targeted RF heating optimization with MVFS algorithm of the spherical tumor of Duke’s head with the individual excitation modes (small maps) and their combined results (large maps) are shown in Figure 7. For 16LDcir and 16LDellip, the MVFS algorithm determined that four excitation modes (M) at three distinct frequencies (300, 350, 450 MHz) was the optimum solution to create the resultant SAR pattern (Figure 7A,C). For 16LDcir_HS and 16LDellip_HS , three modes at two distinct frequencies (350, 450 MHz) were found as the optimum solution (Figure 7B,D). In cases where two modes were found at the same frequency (i.e., 450 MHz), these were added to the resultant total SAR10g by frequency multiplexing. For the 8D RF applicators, the algorithm found two excitation modes at 400 and 450 MHz (8Dcir, 8Dellip and 8Dellip_HS, Figure 7A,B,D) and 350 and 450 MHz (8Dcir_HS, Figure 7C). For the 8L variants 8Lcir, 8Lcir_HS and 8Lellip, three modes at three distinct frequencies were found (Figure 7A–C), and for 8Lellip_HS , two modes at 400 MHz, 450 MHz were found (Figure 7D).
The targeted RF heating optimization local SAR10g maps show that the 16LD and 8D variants are clearly able to deposit RF power inside the tumor TV (Figure 7A–D). Nevertheless, the 8D variants show lower RF power deposition inside the TV and less tumor coverage than the 16LD variants. The 8L variants show even worse performance, with substantially less power deposition, which is mainly concentrated at the periphery of the tumor and head regions due to the limited penetration depth of loop antennas.

3.3. Temperature Simulations

The four 16LD ThermalMR RF applicators are shown in Figure 8A. A detailed view of the SAR10g distribution maps obtained for the 16LD ThermalMR RF applicators is shown for sagittal, coronal and axial orientations (Figure 8B). The steady-state temperature distribution maps derived from the temperature simulations for the four 16LD ThermalMR RF applicators are shown in Figure 8C. The temperature maps show that both horse-shoe variants (16LDcir_HS and 16LDellip_HS) achieved the most uniform focal tumor TV heating, with maximum temperatures (Tmax) of 42.3 °C and 42.2 °C, respectively. The 16LDcir and 16LDellip variants achieved Tmax = 41.5 °C and Tmax = 41.7 °C, respectively, in the tumor TV. However, the 16LDcir and 16LDellip RF applicators also produced undesired temperature hotspots (Tmax = 42.4 °C and Tmax = 43 °C) in the remote facial muscles located close to the right eye. The best heating performance was provided by the 16LDellip_HS RF applicator, which showed the least facial muscle heating (Tmax < 40 °C) and the best tumor TV heating (Tmax = 42.2 °C and Tmean = 41.3 °C).
The transient heating paradigm, which illustrates the kinetics of the maximum tumor TV temperature as a function of the RF power exposure time is illustrated in Figure 9. For all 16LD RF applicators, Tmax of the tumor TV reached an equilibrium at t = 20.5 min. The 16LDcir_HS and 16LDellip_HS arrays showed a temperature rise in the tumor TV from 37.4 °C basal body temperature to a maximum ~42.2 °C with SAR10g of ~40 W/kg. For the same heating paradigm, a temperature increase from basal body temperature (37.4 °C) to Tmax of ~41.5 °C in the tumor TV was observed for the 16LDcir and 16LDellip arrays. This approach suits thermal dose management using the cumulative equivalent minutes (CEM) 43 °C method [17,65,89].

3.4. Thermal Therapy Quality Assessment

SARTV evaluation showed that the 16LD ThermalMR RF applicators had significantly greater SARTV compared with the 8D and 8L variants (Figure 10A). The SARTV obtained for the 16LDcir, 16LDcir_HS, 16LDellip, 16LDellip_HS ThermalMR RF applicators was ~34 W/kg, representing an improvement of 24, 12, 14 and 34% versus the 8D counterparts. The 8Lcir, 8Lcir_HS, 8Lellip, 8Lellip_HS RF applicators achieved substantially less SARTV compared to the 16LD and 8D variants. The values of the mean, maximum, minimum SARTV, and p-values for the statistical comparison are listed in Table 3.
SARmax(TV) analysis demonstrated that 16LDellip and 16LDellip_HS ThermalMR RF applicators afforded a maximum SARmax(TV) = 41 W/kg. This was slightly superior to that of the 16LDcir and 16LDcir_HS variants, which showed maximum SARmax(TV) = 40 W/kg. Notwithstanding this minor difference, all 16LD ThermalMR RF applicators outperformed the 8D and the 8L variants except ellip_HS, which achieved SARmax(TV) = 41 W/kg for the 16LD and 8D horse-shoe configurations. Values are listed in Table 3.
Results of the SAR amplification factor (SAF) are depicted in a color-coded heatmap for all the ThermalMR RF applicators (Figure 10B, left), and spider plots comparing the three designs (16LD, 8D, 8L) across all four configurations (Figure 10B, right). The 16LDellip_HS variant showed the highest SAF (4.78), followed by 16LDcir_HS (SAF = 4.69); the 16LDcir and 16LDellip variants had SAF = 3.95 and SAF = 4.22 (Figure 10B, left). The 8D RF applicators also showed superior SAF values for the horse-shoe configurations, compared to the 360° variants. This observation does not come as a surprise; unlike the 360° variants most of the healthy tissue, e.g., eyes, nose, chin, is not exposed to maximum SAR for the horse-shoe arrays. The overall SAF gain of the 8D variants (8Dcir = 5.39, 8Dcir_HS = 6.25, 8Dellip = 5.36, 8Dellip_HS = 6.02) is plausible because of the overall reduced mean, max and min achieved SAR10g compared with the 16LD variants (Figure 10B, left). All 8L RF applicators showed SAF~2, indicating that there was no greater power deposition in the TV compared with the healthy tissue. By contrast, the 16LD and 8D arrays achieved higher SAF, indicating better SAR deposition in the TV with preservation of remote healthy tissue from SAR exposure.
Quantification of the hotspots-to-target quotient (HTQ) revealed that all the 16LD and 8D RF applicators showed an HTQ of ~0.88 (Figure 10C, left). This indicates that the maximum SAR10g in the hotspot of the healthy tissue does not exceed the mean SAR10g of the TV. The 16LDcir and 16LDcir_hs RF applicators showed HTQ = 0.83, 0.84 and HTQ = 0.8, 0.79 for the 16LDellip and 16LDellip_hs variants. The 8D counterparts showed HTQ = 0.74, 0.72 for the horse-shoe variants 8Dcir_HS and 8Dellip_HS and HTQ = 0.84, 0.88 for 8Dcir, 8Dellip. All 8L variants showed an HTQ > 1: 8Lcir = 2.11, 8Lcir_HS = 2.42, 8Lellip = 2.15, 8Lellip_HS = 2.09. The ellip_HS configuration showed the best results for the 16LD, 8D and 8L variants (Figure 10C, left).
A magnified view of the tumor TV is shown in Figure 11A–D, to better highlight differences in tumor coverage (TCx) of RF power deposition among all ThermalMR RF applicators, along with spider plots showing the TCx quantification in the last column. All 16LD and 8D variants showed T25 = 100%, i.e., 100% of the tumor TV received 25% of peak SAR10g. The 8L variants showed less TC25: 8Lcir = 77%, 8Lcir_HS = 76%, 8Lellip = 78%, 8Lellip_HS = 77%. The 16LD RF applicators outperformed the 8D variants in TC50, TC80 and TC100. All 16LD RF applicators achieved T50 = 100%; the 8Dcir, 8Dcir_HS, 8Dellip and 8Dellip_HS achieved T50 = 87, 97, 83, 95%, respectively; the 8Lcir, 8Lcir_HS, 8Lellip and 8Lellip_HS achieved T50 = 38, 45, 48, 38%, respectively. 16LDcir, 16LDcir_HS, 16LDellip and 16LDellip_HS achieved TC80 = 64, 70, 59, 63%, respectively, which is an 93, 63, 94 and 62% TC80 improvement over the 8D counterparts of TC80 = 33, 43, 31, 39% for 8Dcir, 8Dcir_HS, 8Dellip and 8Dellip_HS , respectively. The 8Lcir, 8Lcir_HS, 8Lellip and 8Lellip_HS achieved much less TC80 = 7, 8, 9, 7%, respectively. 16LDcir, 16LDcir_HS, 16LDellip and 16LDellip_HS showed an 93%, 63%, 94% and 62% TC80 improvement over the 8D counterparts. No 8D or 8L RF applicators achieved any TC100, and 8L showed overall much lower values. But 16LDcir, 16LDcir_HS, 16LDellip and 16LDellip_HS yielded TC100 of 36%, 38%, 30% and 35%, respectively. Among the 16LD RF applicators, the horse-shoe configurations achieved better T80 and TC100 than the 360° variants.
This is also illustrated in the iso-contour temperature maps comparing the four 16LD ThermalMR RF applicators (Figure 12A–D, top), with a magnified view of the temperature achieved within the tumor TV (Figure 12A–D, bottom). The horse-shoe variants 16LDcir_HS and 16LDellip_HS achieved Tmean = 41.5 °C, Tmean = 41.3 °C, respectively, while the 16LDcir and 16LDellip variants showed Tmean = 40.2 °C and Tmean = 40.4 °C inside the tumor. Detailed temperature coverage results show the T40°C, T41°C, T42°C iso-lines within the tumor region. The horse-shoe variants showed superior performance, with T42°C = ~20%, T41°C = ~70%, T40°C = ~10% in the TV (Figure 12B,D, bottom). The 16LDcir and 16LDellip showed T40°C = ~75% and T41°C = ~25% but no T42°C (Figure 12A,C, bottom).
An overview of all the performance metrics, comparing the 16LD ThermalMR applicators among the four configurations is shown in the spider plots in Figure 13. Both horse-shoe-shaped 16LD variants demonstrated superior tumor coverage versus the 360° counterparts. The greater TC50, TC80 and TC100 performance of the 16LD horse-shoe variants facilitated a higher temperature increase from basal body temperature T = 37.4 °C, 16LDcir_HS: ∆T = 4.9 °C, 16LDellip_HS: ∆T = 4.8 °C versus the circular and elliptical 16LD variants.

4. Discussion

Our EMF and temperature simulations in realistic human voxel models of brain tumors add to the literature and demonstrate the efficacy of advanced RF applicators that integrate both loop elements and compact SGBT dipole antennas into hybrid building blocks for ThermalMR theranostics. Our results confirm the hypothesis that ThermalMR RF applicators comprising hybrid loop+dipole building blocks outperform dipole-only and loop-only counterparts in both MRI performance and targeted thermal intervention. Unlike our advanced RF applicators equipped with hybrid loop+SGBT building blocks, loop-only antenna arrays cannot deposit adequate RF power to TV located deep inside the head due to their limited penetration depth in lossy tissue. However, loop antennas support deposition of RF power in peripheral regions of the brain or head, benefiting thermal therapy of TVs located close to the surface. These simulation results provide a crucial technical foundation for the implementation and application of 16 channel LD RF applicators and a springboard for ThermalMR-based theranostics of brain tumors.
Of the 16LD RF applicators, the horse-shoe configurations showed superior tumor coverage and temperature rise inside the tumor compared with the circular 360°-coverage configurations inside the tumor. The horse-shoe configurations have the additional advantages of improved patient comfort by eliminating loop+SGBT building blocks placed close to the eyes, nose and chin, and reduced off-target RF power deposition in the eyes and facial orbit, which benefits RF safety. The antenna arrangement of the horse-shoe variants resulted in a higher SAF and lower HTQ compared with the 360° configurations. Our temperature simulations revealed hotspots in less well-perfused head regions such as the right-side facial muscle tissue, with ~42 °C (16LDcir) and ~43 °C (16LDellip). Nevertheless, these hotspots showed temperatures below the standard upper limit for muscle, fat, and bone tissues (Tlim < 44 °C) [92], which saves healthy tissue from high-temperature exposure. These findings confirm our second hypothesis, that the horse-shoe shaped ThermalMR RF applicators provide better support for MRI and targeted RF heating performance than the 360°-coverage configurations. These advantages suit the clinical needs of ThermalMR theranostics of brain tumors.
We applied a time and frequency multiplexed vector field shaping (MVFS) algorithm to optimize the constructive interference of the electric fields used for targeted RF heating [84]. This algorithm provides multiple excitation modes for appropriate frequencies that do not need to be a priori defined or selected, but rather can adapt to the specific tumor TV for each patient. This is useful for RF applicators, like the self-grounded bow-tie and multi-resonant loop used in our study that are capable of delivering RF power over multiple or broadband frequencies. Although the theoretical maximum number of excitation modes is equal to the total number of antenna elements in the RF applicator, the actual number of excitation modes used in real-world applications of ThermalMR is lower than the maximum. Depending on the ThermalMR RF applicator variants, the algorithm identified specific time and frequency-multiplexed individual modes at different distinct frequencies (250–450 MHz) as the optimum solution to create the overall resultant SAR pattern in the target region. The strength of the MVFS optimization approach is in the utilization of excitations over different frequencies to deposit the maximum power uniformly in the tumor target volume at the desired location. Furthermore, the excitations can be adapted to the size and geometry of tumor TVs to suit the clinical needs of individual patients, and support a ‘personalized medicine’ approach without the need for patient-specific RF applicator hardware.
In this study, we used a single-ring array of LD building blocks, which might constrain MRI to limited anatomical coverage of the head. This potential constraint could be addressed by adding a second ring of LD building blocks along the cranial-caudal axis, which would improve the head coverage of the ThermalMR RF applicator. However, this is not necessarily required since the single-ring RF applicator can be moved along the cranial-caudal axis so that the tumor position is located in the center of the ring, providing a viable and cost-effective approach adapted to the needs of the patient. Incorporating novel metamaterial-surfaces or high dielectric pads offers another alternative approach and research direction to improve MRI performance and heating efficacy [97,98,99,100,101]. Metamaterial planar surfaces (metasurface) are conceptually appealing for pursuing the development of RF applicators for ThermalMR-based theranostics of brain tumors due to the extra degrees of freedom for shaping electromagnetic fields they can provide [97,98,99,100,101]. This benefit could be exploited for further SAR10g reduction and B1+ uniformity improvements for MRI and targeted RF power deposition in deep-seated brain tumors [101,102].
Our thermal intervention results reveal the relationship between temperature change, SAR and the effective perfusion. If the tumor TV includes large blood vessels or anatomical regions with high perfusion, more RF power deposition would be required to achieve the desired temperature increase. To increase the temperature in the brain beyond the (Tmax = ~42.3 °C) achieved in the current study would require SAR > 40 W/kg. Z. Rijnen et al. reported the use of RF power deposition up to 96–178 W/kg inside the tumor region in head and neck cancer treatment [3,22]. However, additional SAR deposition might create more unwanted hotspots. Implementation of a water bolus with cooling circulation placed between the head and the RF array to cool down the less well-perfused facial muscle regions could help limit hotspot generation. Several studies have suggested the use of a water bolus to improve the quality of thermal therapy [71,77,78,92,103,104]. Another way to reduce SAR in unwanted hotspots is to incorporate novel metamaterial absorbers between the RF applicator and head to absorb SAR from hotspots in healthy tissue [102]. Our advanced hybrid loop-dipole ThermalMR RF applicators are also compatible with these approaches.
The physical processes of E-field interference and heat distribution inside the body and brain are complex and heterogeneous. This limits the accuracy with which temperature distributions in heterogeneous perfused tissue can be predicted in thermal models. Nevertheless, numerous applications can benefit from thermal modeling [29,103]. The limitations of existing RF hyperthermia systems can be overcome by developing new simulation tools for optimizing clinical treatments, considering temperature-dependent effects on blood perfusion, by using noninvasive temperature measurement with MR-thermometry. This measurement is also facilitated by an appropriately designed ThermalMR RF applicator. Developing a realistic tumor model remains a challenge for numerical simulations, due to the complex pathophysiology of tumors. The structure of tumor vessels exhibits irregular branching and twisting where blood flow does not follow a constant, unidirectional path and not all vessels are perfused continuously [73]. These complexities limit the accuracy of blood perfusion coefficients used in the temperature simulations of the tumor model. To advance thermal simulations of tumor models, patient-specific vasculature networks can be derived from MR angiography and incorporated into the model [3,85,103]. Also, patient-specific perfusion maps can be derived from perfusion-weighted MRI. Tumors can also exhibit disrupted cellular and extracellular composition with altered dielectric properties that can be heterogeneous even within the tumor itself [74]. MR-based tomography of electric properties can be employed to better model the temperature-dependent conductivity and permittivity of tumors in the simulations [105]. Notwithstanding these potential future advances, numerical EMF and temperature simulations have become essential for assessing the performance of ThermalMR RF applicators and for temperature dose management. The strength of our numerical thermal modeling is that prior to a treatment session, a selection of patient-specific ThermalMR RF applicator settings can be determined, together with specific RF power optimization planning and thermal dose management. This can be used to maximize temperature in the target volume using SAR as a primary input for the planning, evaluation, and optimization of ThermalMR-based theranostics of brain tumors [71,85,94].

5. Conclusions

This work demonstrates the feasibility and applicability of 16-channel hybrid ThermalMR RF applicators combining compact SGBT dipole and loop antennas for ThermalMR theranostics of deep-seated brain tumors at ultrahigh magnetic field strengths of 7.0 T, 9.4 T and 10.5 T. Targeted RF heating of brain tumors can be adapted in a patient-specific manner without changing the RF applicator hardware, by utilizing the multi-resonant and broadband characteristics of the loop+SGBT dipole building blocks in conjunction with the time and frequency multiplexed vector field shaping algorithm. The enhanced RF power deposition and temperature rise up to ~42.3 °C inside the tumor target volume meet the temperature criteria for adjunct hyperthermia therapy of glioblastoma multiforme. Our EMF and temperature simulations establish a solid technical foundation for the development and construction of hybrid loop+SGBT dipole ThermalMR RF applicators. Such numerical simulations are a mandatory precursor for future in vivo studies and provide a rigorous framework en route to clinical ThermalMR-based theranostics of brain tumors. This also has implications for investigating fundamental questions in biology, molecular medicine and molecular imaging. ThermalMR can potentially be used to investigate (patho)physiological processes, opening an entirely new research field of thermal phenotyping of cancer. “Can stages of tumor progression and therapy response be characterized by thermal profiles?”; “What are the links between MRI biomarkers, thermal profiles and the molecular signatures of brain tumors?” ThermalMR has the potential to offer an answer to such questions.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15082303/s1, Table S1: The dielectric and thermal properties of the tumor and healthy tissues used for the human voxel model Duke.

Author Contributions

Conceptualization: N.S., A.K. and T.N.; methodology: N.S., A.K. and J.M.M.; investigation: N.S.; visualization: N.S., T.N. and S.W.; supervision: T.N. and A.K.; writing original draft: N.S.; writing review and editing: N.S., T.N., J.M.M., S.W., L.S. and T.W.E. All authors have read and agreed to the published version of the manuscript.

Funding

This project is funded by the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Program, grant agreement No. 743077 (ThermalMR) and from Hyperboost, an Innovative Training Network (ITN) H2020-MSCA-ITN-2020-955625 of the Marie Skłodowska-Curie Actions of the European Union.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data can be found in the text.

Conflicts of Interest

Thoralf Niendorf is founder and CEO of MRI.TOOLS GmbH, Berlin, Germany.

Abbreviations

MRmagnetic resonance
MRImagnetic resonance imaging
ThermalMRthermal magnetic resonance
MRThMR thermometry
RFradiofrequency
UHFultrahigh field
HThyperthermia
RTradiotherapy
WHOWorld Health Organization
GBMglioblastoma multiforme
CNScentral nervous system
TMTPthermalMR treatment planning
E-fieldelectric Field
B1+RF transmission field used for MRI
B₀main static magnetic field of MRI scanner
SARspecific absorption rate
SAR10g10 g averaged SAR distribution
EMFelectromagnetic field
SGBTself-grounded bow-tie
LDloop+SGBT dipole building block
Lloop building block
DSGBT dipole building block
16LDhybrid 16-channel loop+SGBT dipole ThermalMR RF applicator
8D8-channel SGBT dipole ThermalMR RF applicator
8L8-channel loop ThermalMR RF applicator
circircular array
cir_HScircular horse-shoe array
ellipelliptical array
ellip_HSelliptical horse-shoe array
16LDcir circular array variant of 16LD
16LDcir_HS circular horse-shoe array variant of 16LD
16LDellip elliptical array variant of 16LD
16LDellip_HS elliptical horse-shoe array variant of 16LD
8Dcir circular array variant of 8D
8Dcir_HS circular horse-shoe array variant of 8D
8Dellip elliptical array variant of 8D
8Dellip_HS elliptical horse-shoe array variant of 8D
8Lcir circular array variant of 8L
8Lcir_HS circular horse-shoe array variant of 8L
8Lellip elliptical array variant of 8L
8Lellip_HS elliptical horse-shoe array variant of 8L
MVFSmultiplexed vector field shaping
Mexcitation mode provided by MVFS algorithm
PLDpower loss density
ROIregion of interest
TVtarget volume
TCtumor coverage
TC25fraction of the tumor volume enclosed by 25% peak SARmax isolines
TC50fraction of the tumor volume enclosed by 50% peak SARmax isolines
TC80fraction of the tumor volume enclosed by 80% peak SARmax isolines
TC100fraction of the tumor volume enclosed by 100% peak SARmax isolines
HTQHotspot-to-target quotient
SAFSAR amplification factor
ESHOEuropean Society for Hyperthermic Oncology
Ttemperature
Tmean mean temperature inside TV
Tmax maximum temperature inside TV
Tx°Cindex temperature coverage inside TV
T40°Cfraction of the tumor volume enclosed by 40 °C
T41°Cfraction of the tumor volume enclosed by 41 °C
T42°Cfraction of the tumor volume enclosed by 42 °C
Maxmaximum
Minminimum
CEMcumulative equivalent minutes
FoVfield-of-view

References

  1. Lee Titsworth, W.; Murad, G.J.; Hoh, B.L.; Rahman, M. Fighting fire with fire: The revival of thermotherapy for gliomas. Anticancer Res. 2014, 34, 565–574. [Google Scholar] [PubMed]
  2. Wust, P.; Hildebrandt, B.; Sreenivasa, G.; Rau, B.; Gellermann, J.; Riess, H.; Felix, R.; Schlag, P.M. Hyperthermia in combined treatment of cancer. Lancet Oncol. 2002, 3, 487–497. [Google Scholar] [CrossRef] [PubMed]
  3. Kok, H.P.; Crezee, J. Hyperthermia treatment planning: Clinical application and ongoing research. In Proceedings of the 2020 14th European Conference on Antennas and Propagation (EuCAP), Copenhagen, Denmark, 15–20 March 2020; pp. 1–5. [Google Scholar]
  4. Horsman, M.R.; Overgaard, J. Hyperthermia: A potent enhancer of radiotherapy. Clin. Oncol. 2007, 19, 418–426. [Google Scholar] [CrossRef] [PubMed]
  5. Issels, R.D.; Lindner, L.H.; Verweij, J.; Wust, P.; Reichardt, P.; Schem, B.C.; Abdel-Rahman, S.; Daugaard, S.; Salat, C.; Wendtner, C.M.; et al. Neo-adjuvant chemotherapy alone or with regional hyperthermia for localised high-risk soft-tissue sarcoma: A randomised phase 3 multicentre study. Lancet Oncol. 2010, 11, 561–570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sminia, P.; van der Zee, J.; Wondergem, J.; Haveman, J. Effect of hyperthermia on the central nervous system: A review. Int. J. Hyperth. 1994, 10, 1–30. [Google Scholar] [CrossRef]
  7. Sneed, P.K.; Stauffer, P.R.; McDermott, M.W.; Diederich, C.J.; Lamborn, K.R.; Prados, M.D.; Chang, S.; Weaver, K.A.; Spry, L.; Malec, M.K.; et al. Survival benefit of hyperthermia in a prospective randomized trial of brachytherapy boost +/- hyperthermia for glioblastoma multiforme. Int. J. Radiat. Oncol. Biol. Phys. 1998, 40, 287–295. [Google Scholar] [CrossRef]
  8. Seegenschmiedt, M.H.; Karlsson, U.L.; Black, P.; Brady, L.W. Thermoradiotherapy for brain tumors. Three cases of recurrent malignant astrocytoma and review of clinical experience. Am. J. Clin. Oncol. 1995, 18, 510–518. [Google Scholar] [CrossRef]
  9. Yang, Z.; Gao, D.; Zhao, J.; Yang, G.; Guo, M.; Wang, Y.; Ren, X.; Kim, J.S.; Jin, L.; Tian, Z.; et al. Thermal immuno-nanomedicine in cancer. Nat. Rev. Clin. Oncol. 2023, 20, 116–134. [Google Scholar] [CrossRef]
  10. Jordan, A.; Scholz, R.; Maier-Hauff, K.; Johannsen, M.; Wust, P.; Nadobny, J.; Schirra, H.; Schmidt, H.; Deger, S.; Loening, S.; et al. Presentation of a new magnetic field therapy system for the treatment of human solid tumors with magnetic fluid hyperthermia. J. Magn. Magn. Mater. 2001, 225, 118–126. [Google Scholar] [CrossRef] [Green Version]
  11. Kok, H.P.; Crezee, J. Adapt2Heat: Treatment planning-assisted locoregional hyperthermia by on-line visualization, optimization and re-optimization of SAR and temperature distributions. Int. J. Hyperth. 2022, 39, 265–277. [Google Scholar] [CrossRef]
  12. Hersh, A.M.; Bhimreddy, M.; Weber-Levine, C.; Jiang, K.; Alomari, S.; Theodore, N.; Manbachi, A.; Tyler, B.M. Applications of Focused Ultrasound for the Treatment of Glioblastoma: A New Frontier. Cancers 2022, 14, 4920. [Google Scholar] [CrossRef] [PubMed]
  13. Huang, P.C.; Chaney, E.J.; Aksamitiene, E.; Barkalifa, R.; Spillman, D.R., Jr.; Bogan, B.J.; Boppart, S.A. Biomechanical sensing of in vivo magnetic nanoparticle hyperthermia-treated melanoma using magnetomotive optical coherence elastography. Theranostics 2021, 11, 5620–5633. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, X.; Zhang, Y.; Wang, Y.; Zhu, W.; Li, G.; Ma, X.; Zhang, Y.; Chen, S.; Tiwari, S.; Shi, K.; et al. Comprehensive understanding of magnetic hyperthermia for improving antitumor therapeutic efficacy. Theranostics 2020, 10, 3793–3815. [Google Scholar] [CrossRef]
  15. Kok, H.P.; Korshuize-van Straten, L.; Bakker, A.; de Kroon-Oldenhof, R.; Westerveld, G.H.; Versteijne, E.; Stalpers, L.J.A.; Crezee, J. Feasibility of on-line temperature-based hyperthermia treatment planning to improve tumour temperatures during locoregional hyperthermia. Int. J. Hyperth. 2018, 34, 1082–1091. [Google Scholar] [CrossRef] [Green Version]
  16. Sheybani, N.D.; Batts, A.J.; Mathew, A.S.; Thim, E.A.; Price, R.J. Focused Ultrasound Hyperthermia Augments Release of Glioma-derived Extracellular Vesicles with Differential Immunomodulatory Capacity. Theranostics 2020, 10, 7436–7447. [Google Scholar] [CrossRef]
  17. Kok, H.P.; Cressman, E.N.K.; Ceelen, W.; Brace, C.L.; Ivkov, R.; Grull, H.; Ter Haar, G.; Wust, P.; Crezee, J. Heating technology for malignant tumors: A review. Int. J. Hyperth. 2020, 37, 711–741. [Google Scholar] [CrossRef] [PubMed]
  18. Canters, R.A.; Paulides, M.M.; Franckena, M.F.; van der Zee, J.; van Rhoon, G.C. Implementation of treatment planning in the routine clinical procedure of regional hyperthermia treatment of cervical cancer: An overview and the Rotterdam experience. Int. J. Hyperth. 2012, 28, 570–581. [Google Scholar] [CrossRef]
  19. Fang, Y.; He, Y.; Wu, C.; Zhang, M.; Gu, Z.; Zhang, J.; Liu, E.; Xu, Q.; Asrorov, A.M.; Huang, Y. Magnetism-mediated targeting hyperthermia-immunotherapy in "cold" tumor with CSF1R inhibitor. Theranostics 2021, 11, 6860–6872. [Google Scholar] [CrossRef]
  20. Kato, H.; Ishida, T. Present and future status of noninvasive selective deep heating using RF in hyperthermia. Med. Biol. Eng. Comput. 1993, 31, S2–S11. [Google Scholar] [CrossRef]
  21. Sullivan, D.M.; Ben-Yosef, R.; Kapp, D.S. Stanford 3D hyperthermia treatment planning system. Technical review and clinical summary. Int. J. Hyperth. 1993, 9, 627–643. [Google Scholar] [CrossRef]
  22. Rijnen, Z.; Bakker, J.F.; Canters, R.A.; Togni, P.; Verduijn, G.M.; Levendag, P.C.; Van Rhoon, G.C.; Paulides, M.M. Clinical integration of software tool VEDO for adaptive and quantitative application of phased array hyperthermia in the head and neck. Int. J. Hyperth. 2013, 29, 181–193. [Google Scholar] [CrossRef] [PubMed]
  23. Kiyatkin, E.A.; Sharma, H.S. Permeability of the blood-brain barrier depends on brain temperature. Neuroscience 2009, 161, 926–939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Tabatabaei, S.N.; Girouard, H.; Carret, A.S.; Martel, S. Remote control of the permeability of the blood-brain barrier by magnetic heating of nanoparticles: A proof of concept for brain drug delivery. J. Control. Release 2015, 206, 49–57. [Google Scholar] [CrossRef] [PubMed]
  25. Brito, B.; Price, T.W.; Gallo, J.; Banobre-Lopez, M.; Stasiuk, G.J. Smart magnetic resonance imaging-based theranostics for cancer. Theranostics 2021, 11, 8706–8737. [Google Scholar] [CrossRef]
  26. Rabha, B.; Bharadwaj, K.K.; Pati, S.; Choudhury, B.K.; Sarkar, T.; Kari, Z.A.; Edinur, H.A.; Baishya, D.; Atanase, L.I. Development of Polymer-Based Nanoformulations for Glioblastoma Brain Cancer Therapy and Diagnosis: An Update. Polymers 2021, 13, 4114. [Google Scholar] [CrossRef]
  27. Ji, Y.; Winter, L.; Navarro, L.; Ku, M.C.; Periquito, J.S.; Pham, M.; Hoffmann, W.; Theune, L.E.; Calderon, M.; Niendorf, T. Controlled Release of Therapeutics from Thermoresponsive Nanogels: A Thermal Magnetic Resonance Feasibility Study. Cancers 2020, 12, 1380. [Google Scholar] [CrossRef]
  28. Kneidl, B.; Peller, M.; Winter, G.; Lindner, L.H.; Hossann, M. Thermosensitive liposomal drug delivery systems: State of the art review. Int. J. Nanomed. 2014, 9, 4387–4398. [Google Scholar] [CrossRef] [Green Version]
  29. Paulides, M.M.; Dobsicek Trefna, H.; Curto, S.; Rodrigues, D.B. Recent technological advancements in radiofrequency- andmicrowave-mediated hyperthermia for enhancing drug delivery. Adv. Drug Deliv. Rev. 2020, 163–164, 3–18. [Google Scholar] [CrossRef]
  30. Huang, Q.; Hu, J.K.; Lohr, F.; Zhang, L.; Braun, R.; Lanzen, J.; Little, J.B.; Dewhirst, M.W.; Li, C.Y. Heat-induced gene expression as a novel targeted cancer gene therapy strategy. Cancer Res. 2000, 60, 3435–3439. [Google Scholar]
  31. Bernstein, J.G.; Garrity, P.A.; Boyden, E.S. Optogenetics and thermogenetics: Technologies for controlling the activity of targeted cells within intact neural circuits. Curr. Opin. Neurobiol. 2012, 22, 61–71. [Google Scholar] [CrossRef] [Green Version]
  32. Alexandrescu, S.; Korshunov, A.; Lai, S.H.; Dabiri, S.; Patil, S.; Li, R.; Shih, C.S.; Bonnin, J.M.; Baker, J.A.; Du, E.; et al. Epithelioid Glioblastomas and Anaplastic Epithelioid Pleomorphic Xanthoastrocytomas--Same Entity or First Cousins? Brain Pathol. 2016, 26, 215–223. [Google Scholar] [CrossRef] [PubMed]
  33. Louis, D.N.; Perry, A.; Reifenberger, G.; von Deimling, A.; Figarella-Branger, D.; Cavenee, W.K.; Ohgaki, H.; Wiestler, O.D.; Kleihues, P.; Ellison, D.W. The 2016 World Health Organization Classification of Tumors of the Central Nervous System: A summary. Acta Neuropathol. 2016, 131, 803–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Roussakow, S.V. Clinical and economic evaluation of modulated electrohyperthermia concurrent to dose-dense temozolomide 21/28 days regimen in the treatment of recurrent glioblastoma: A retrospective analysis of a two-centre German cohort trial with systematic comparison and effect-to-treatment analysis. BMJ Open 2017, 7, e017387. [Google Scholar] [CrossRef] [Green Version]
  35. Grech, N.; Dalli, T.; Mizzi, S.; Meilak, L.; Calleja, N.; Zrinzo, A. Rising Incidence of Glioblastoma Multiforme in a Well-Defined Population. Cureus 2020, 12, e8195. [Google Scholar] [CrossRef]
  36. Tan, A.C.; Ashley, D.M.; Lopez, G.Y.; Malinzak, M.; Friedman, H.S.; Khasraw, M. Management of glioblastoma: State of the art and future directions. CA Cancer J. Clin. 2020, 70, 299–312. [Google Scholar] [CrossRef]
  37. Oberacker, E.; Kuehne, A.; Oezerdem, C.; Nadobny, J.; Weihrauch, M.; Beck, M.; Zschaeck, S.; Diesch, C.; Eigentler, T.W.; Waiczies, H.; et al. Radiofrequency applicator concepts for thermal magnetic resonance of brain tumors at 297 MHz (7.0 Tesla). Int. J. Hyperth. 2020, 37, 549–563. [Google Scholar] [CrossRef]
  38. Seebass, M.; Beck, R.; Gellermann, J.; Nadobny, J.; Wust, P. Electromagnetic phased arrays for regional hyperthermia: Optimal frequency and antenna arrangement. Int. J. Hyperth. 2001, 17, 321–336. [Google Scholar] [CrossRef]
  39. Paulides, M.M.; Curto, S.; Wu, M.; Winter, L.; Rhoon, G.C.v.; Yeo, D.T.B. Advances in magnetic resonance guided radiofrequency hyperthermia. In Proceedings of the 2017 11th European Conference on Antennas and Propagation (EUCAP), Paris, France, 19–24 March 2017; pp. 3692–3696. [Google Scholar]
  40. Prinz, C.; Starke, L.; Millward, J.M.; Fillmer, A.; Delgado, P.R.; Waiczies, H.; Pohlmann, A.; Rothe, M.; Nazaré, M.; Paul, F.; et al. In vivo detection of teriflunomide-derived fluorine signal during neuroinflammation using fluorine MR spectroscopy. Theranostics 2021, 11, 2490–2504. [Google Scholar] [CrossRef] [PubMed]
  41. Starke, L.; Millward, J.M.; Prinz, C.; Sherazi, F.; Waiczies, H.; Lippert, C.; Nazaré, M.; Paul, F.; Niendorf, T.; Waiczies, S. First in vivo fluorine-19 magnetic resonance imaging of the multiple sclerosis drug siponimod. Theranostics 2023, 13, 1217–1234. [Google Scholar] [CrossRef]
  42. Paulides, M.M.; Bakker, J.F.; Neufeld, E.; van der Zee, J.; Jansen, P.P.; Levendag, P.C.; van Rhoon, G.C. Winner of the “New Investigator Award” at the European Society of Hyperthermia Oncology Meeting 2007. The HYPERcollar: A novel applicator for hyperthermia in the head and neck. Int. J. Hyperth. 2007, 23, 567–576. [Google Scholar] [CrossRef]
  43. Togni, P.; Rijnen, Z.; Numan, W.C.; Verhaart, R.F.; Bakker, J.F.; van Rhoon, G.C.; Paulides, M.M. Electromagnetic redesign of the HYPERcollar applicator: Toward improved deep local head-and-neck hyperthermia. Phys. Med. Biol. 2013, 58, 5997–6009. [Google Scholar] [CrossRef] [PubMed]
  44. Oberacker, E.; Diesch, C.; Nadobny, J.; Kuehne, A.; Wust, P.; Ghadjar, P.; Niendorf, T. Patient-Specific Planning for Thermal Magnetic Resonance of Glioblastoma Multiforme. Cancers 2021, 13, 1867. [Google Scholar] [CrossRef] [PubMed]
  45. Winter, L.; Niendorf, T. Electrodynamics and radiofrequency antenna concepts for human magnetic resonance at 23.5 T (1 GHz) and beyond. MAGMA 2016, 29, 641–656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Winter, L.; Ozerdem, C.; Hoffmann, W.; Santoro, D.; Muller, A.; Waiczies, H.; Seemann, R.; Graessl, A.; Wust, P.; Niendorf, T. Design and evaluation of a hybrid radiofrequency applicator for magnetic resonance imaging and RF induced hyperthermia: Electromagnetic field simulations up to 14.0 Tesla and proof-of-concept at 7.0 Tesla. PLoS ONE 2013, 8, e61661. [Google Scholar] [CrossRef]
  47. Niendorf, T.; Oezerdem, C.; Ji, Y.; Oberacker, E.; Kuehne, A.; Waiczies, H.; Winter, L. Radiative RF antenna arrays for cardiac, brain and thermal magnetic resonance at ultrahigh and extreme magnetic field strengths: Concepts, electromagnetic field simulations and applications. In Proceedings of the 2017 International Conference on Electromagnetics in Advanced Applications (ICEAA), Verona, Italy, 11–15 September 2017; pp. 1567–1570. [Google Scholar]
  48. Han, H.; Eigentler, T.W.; Wang, S.; Kretov, E.; Winter, L.; Hoffmann, W.; Grass, E.; Niendorf, T. Design, Implementation, Evaluation and Application of a 32-Channel Radio Frequency Signal Generator for Thermal Magnetic Resonance Based Anti-Cancer Treatment. Cancers 2020, 12, 1720. [Google Scholar] [CrossRef]
  49. Han, H.; Oberacker, E.; Kuehne, A.; Wang, S.; Eigentler, T.W.; Grass, E.; Niendorf, T. Multi-Channel RF Supervision Module for Thermal Magnetic Resonance Based Cancer Therapy. Cancers 2021, 13, 1001. [Google Scholar] [CrossRef]
  50. Winter, L.; Oberacker, E.; Paul, K.; Ji, Y.; Oezerdem, C.; Ghadjar, P.; Thieme, A.; Budach, V.; Wust, P.; Niendorf, T. Magnetic resonance thermometry: Methodology, pitfalls and practical solutions. Int. J. Hyperth. 2016, 32, 63–75. [Google Scholar] [CrossRef] [Green Version]
  51. Davis, R.M.; Viglianti, B.L.; Yarmolenko, P.; Park, J.Y.; Stauffer, P.; Needham, D.; Dewhirst, M.W. A method to convert MRI images of temperature change into images of absolute temperature in solid tumours. Int. J. Hyperth. 2013, 29, 569–581. [Google Scholar] [CrossRef] [Green Version]
  52. Rieke, V.; Butts Pauly, K. MR thermometry. J. Magn. Reson. Imaging 2008, 27, 376–390. [Google Scholar] [CrossRef]
  53. Odeen, H.; Parker, D.L. Magnetic resonance thermometry and its biological applications-Physical principles and practical considerations. Prog. Nucl. Magn. Reson. Spectrosc. 2019, 110, 34–61. [Google Scholar] [CrossRef]
  54. Denis de Senneville, B.; Quesson, B.; Moonen, C.T. Magnetic resonance temperature imaging. Int. J. Hyperth. 2005, 21, 515–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Pegah Takook, H.D.T.; Zeng, X.; Fhager, A.; Persson, M. A Computational Study Using Time Reversal Focusing for Hyperthermia Treatment Planning. Prog. Electromagn. Res. B 2017, 73, 117–130. [Google Scholar] [CrossRef] [Green Version]
  56. Takook, P.; Persson, M.; Gellermann, J.; Trefna, H.D. Compact self-grounded Bow-Tie antenna design for an UWB phased-array hyperthermia applicator. Int. J. Hyperth. 2017, 33, 387–400. [Google Scholar] [CrossRef]
  57. Takook, P.; Persson, M.; Trefná, H.D. Performance Evaluation of Hyperthermia Applicators to Heat Deep-Seated Brain Tumors. IEEE J. Electromagn. RF Microw. Med. Biol. 2018, 2, 18–24. [Google Scholar] [CrossRef]
  58. Takook, P.; Shafiemehr, M.; Persson, M.; Trefná, H.D. Experimental evaluation of UWB applicator prototype for head and neck hyperthermia. In Proceedings of the 2017 11th European Conference on Antennas and Propagation (EUCAP), Paris, France, 19–24 March 2017; pp. 3619–3620. [Google Scholar]
  59. Takook, P.; Trefná, H.D.; Persson, M. Performance evaluation of 2 hyperthermia applicators for deep-seated brain tumors. In Proceedings of the 2017 First IEEE MTT-S International Microwave Bio Conference (IMBIOC), Gothenburg, Sweden, 15–17 May 2017; pp. 1–3. [Google Scholar]
  60. Trefná, H.D.; Imtiaz, A.; Lui, H.S.; Rubæk, T.; Persson, M. Evolution of an UWB antenna for hyperthermia array applicator. In Proceedings of the 2012 6th European Conference on Antennas and Propagation (EUCAP), Prague, Czech Republic, 26–30 March 2012; pp. 1046–1048. [Google Scholar]
  61. Zanoli, M.; Ek, E.; Dobšíček Trefná, H. Antenna Arrangement in UWB Helmet Brain Applicators for Deep Microwave Hyperthermia. Cancers 2023, 15, 1447. [Google Scholar] [CrossRef] [PubMed]
  62. Erturk, M.A.; Raaijmakers, A.J.; Adriany, G.; Ugurbil, K.; Metzger, G.J. A 16-channel combined loop-dipole transceiver array for 7 Tesla body MRI. Magn. Reson. Med. 2017, 77, 884–894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Woo, M.K.; Delabarre, L.; Lee, B.Y.; Waks, M.; Lagore, R.L.; Radder, J.W.; Eryaman, Y.; Ugurbil, K.; Adriany, G. Evaluation of a 16-Channel Transceiver Loop + Dipole Antenna Array for Human Head Imaging at 10.5 Tesla. IEEE Access 2020, 8, 203555–203563. [Google Scholar] [CrossRef]
  64. Raaijmakers, A.J.; Luijten, P.R.; van den Berg, C.A. Dipole antennas for ultrahigh-field body imaging: A comparison with loop coils. NMR Biomed. 2016, 29, 1122–1130. [Google Scholar] [CrossRef]
  65. Lattanzi, R.; Wiggins, G.C.; Zhang, B.; Duan, Q.; Brown, R.; Sodickson, D.K. Approaching ultimate intrinsic signal-to-noise ratio with loop and dipole antennas. Magn. Reson. Med. 2018, 79, 1789–1803. [Google Scholar] [CrossRef]
  66. Pfrommer, A.; Henning, A. The ultimate intrinsic signal-to-noise ratio of loop- and dipole-like current patterns in a realistic human head model. Magn. Reson. Med. 2018, 80, 2122–2138. [Google Scholar] [CrossRef]
  67. Hurshkainen, A.; Steensma, B.R.; Glybovski, S.B.; Voogt, I.J.; Melchakova, I.V.; Belov, P.A.; Berg, C.A.T.v.d.; Raaijmakers, A.J.E. On Optimization of Radiative Dipole Body Array Coils for 7 Tesla MRI. arXiv 2017, arXiv:1710.02399. [Google Scholar]
  68. Erturk, M.A.; Wu, X.; Eryaman, Y.; Van de Moortele, P.F.; Auerbach, E.J.; Lagore, R.L.; DelaBarre, L.; Vaughan, J.T.; Ugurbil, K.; Adriany, G.; et al. Toward imaging the body at 10.5 tesla. Magn. Reson. Med. 2017, 77, 434–443. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Eryaman, Y.; Guerin, B.; Keil, B.; Mareyam, A.; Herraiz, J.L.; Kosior, R.K.; Martin, A.; Torrado-Carvajal, A.; Malpica, N.; Hernandez-Tamames, J.A.; et al. SAR reduction in 7T C-spine imaging using a ”dark modes“ transmit array strategy. Magn. Reson. Med. 2015, 73, 1533–1539. [Google Scholar] [CrossRef] [Green Version]
  70. Lagore, R.L.; Moeller, S.; Zimmermann, J.; DelaBarre, L.; Radder, J.; Grant, A.; Ugurbil, K.; Yacoub, E.; Harel, N.; Adriany, G. An 8-dipole transceive and 24-loop receive array for non-human primate head imaging at 10.5 T. NMR Biomed. 2021, 34, e4472. [Google Scholar] [CrossRef]
  71. Prasad, B.; Kim, J.K.; Kim, S. Role of Simulations in the Treatment Planning of Radiofrequency Hyperthermia Therapy in Clinics. J. Oncol. 2019, 2019, 9685476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Foundation, I.I. Virtual Population (ViP) Models:DUKE. Available online: https://itis.swiss/virtual-population/virtual-population/overview/ (accessed on 1 January 2022).
  73. Zanoli, M.; Trefna, H.D. Iterative time-reversal for multi-frequency hyperthermia. Phys. Med. Biol. 2021, 66, 045027. [Google Scholar] [CrossRef]
  74. Lu, Y.; Li, B.; Xu, J.; Yu, J. Dielectric properties of human glioma and surrounding tissue. Int. J. Hyperth. 1992, 8, 755–760. [Google Scholar] [CrossRef]
  75. Hasgall, P.A.; Di Gennaro, F.; Baumgartner, C.; Neufeld, E.; Lloyd, B.; Gosselin, M.C.; Payne, D.; Klingenböck, A.; Kuster, N. IT’IS Database for Thermal and Electromagnetic Parameters of Biological Tissues; IT’IS Foundation: Zurich, Switzerland, 2022. [Google Scholar] [CrossRef]
  76. Eigentler, T.W.; Winter, L.; Han, H.; Oberacker, E.; Kuehne, A.; Waiczies, H.; Schmitter, S.; Boehmert, L.; Prinz, C.; Trefna, H.D.; et al. Wideband Self-Grounded Bow-Tie Antenna for Thermal MR. NMR Biomed. 2020, 33, e4274. [Google Scholar] [CrossRef] [Green Version]
  77. Rijnen, Z.; Togni, P.; Roskam, R.; van de Geer, S.G.; Goossens, R.H.; Paulides, M.M. Quality and comfort in head and neck hyperthermia: A redesign according to clinical experience and simulation studies. Int. J. Hyperth. 2015, 31, 823–830. [Google Scholar] [CrossRef] [Green Version]
  78. Van der Gaag, M.L.; De Bruijne, M.; Samaras, T.; Van der Zee, J.; Van Rhoon, G.C. Development of a guideline for the water bolus temperature in superficial hyperthermia. Int. J. Hyperth. 2006, 22, 637–656. [Google Scholar] [CrossRef]
  79. Jin-Yuan, C.; Gandhi, O.P. Numerical simulation of annular-phased arrays of dipoles for hyperthermia of deep-seated tumors. IEEE Trans. Biomed. Eng. 1992, 39, 209–216. [Google Scholar] [CrossRef] [PubMed]
  80. Roemer, P.B.; Edelstein, W.A.; Hayes, C.E.; Souza, S.P.; Mueller, O.M. The NMR phased array. Magn. Reson. Med. 1990, 16, 192–225. [Google Scholar] [CrossRef]
  81. Kuehne, A.; Goluch, S.; Waxmann, P.; Seifert, F.; Ittermann, B.; Moser, E.; Laistler, E. Power balance and loss mechanism analysis in RF transmit coil arrays. Magn. Reson. Med. 2015, 74, 1165–1176. [Google Scholar] [CrossRef]
  82. Graessl, A.; Renz, W.; Hezel, F.; Dieringer, M.A.; Winter, L.; Oezerdem, C.; Rieger, J.; Kellman, P.; Santoro, D.; Lindel, T.D.; et al. Modular 32-channel transceiver coil array for cardiac MRI at 7.0T. Magn. Reson. Med. 2014, 72, 276–290. [Google Scholar] [CrossRef] [PubMed]
  83. 60601-2-33:2010; Medical Electrical Equipment-Part 2-33: Particular Requirements for the Basic Safety and Essential Performance of Magnetic Resonance Equipment for Medical Diagnosis. IEC: Genève, Switzerland, 2015. Available online: https://mriquestions.com/uploads/3/4/5/7/34572113/safety_iec_60601-2-33previews_1897819_pre.pdf (accessed on 30 November 2022).
  84. Kuehne, A.; Oberacker, E.; Waiczies, H.; Niendorf, T. Solving the Time- and Frequency-Multiplexed Problem of Constrained Radiofrequency Induced Hyperthermia. Cancers 2020, 12, 1072. [Google Scholar] [CrossRef] [PubMed]
  85. Kok, H.P.; Kotte, A.; Crezee, J. Planning, optimisation and evaluation of hyperthermia treatments. Int. J. Hyperth. 2017, 33, 593–607. [Google Scholar] [CrossRef] [Green Version]
  86. Jha, S.; Sharma, P.K.; Malviya, R. Hyperthermia: Role and Risk Factor for Cancer Treatment. Achiev. Life Sci. 2016, 10, 161–167. [Google Scholar] [CrossRef] [Green Version]
  87. Tilly, W.; Wust, P.; Rau, B.; Harder, C.; Gellermann, J.; Schlag, P.; Budach, V.; Felix, R. Temperature data and specific absorption rates in pelvic tumours: Predictive factors and correlations. Int. J. Hyperth. 2001, 17, 172–188. [Google Scholar] [CrossRef]
  88. Oh, S.; Ryu, Y.C.; Carluccio, G.; Sica, C.T.; Collins, C.M. Measurement of SAR-induced temperature increase in a phantom and in vivo with comparison to numerical simulation. Magn. Reson. Med. 2014, 71, 1923–1931. [Google Scholar] [CrossRef] [Green Version]
  89. Li, Z.; Vogel, M.; Maccarini, P.F.; Stakhursky, V.; Soher, B.J.; Craciunescu, O.I.; Das, S.; Arabe, O.A.; Joines, W.T.; Stauffer, P.R. Improved hyperthermia treatment control using SAR/temperature simulation and PRFS magnetic resonance thermal imaging. Int. J. Hyperth. 2011, 27, 86–99. [Google Scholar] [CrossRef] [Green Version]
  90. Kok, H.P.; Schooneveldt, G.; Bakker, A.; de Kroon-Oldenhof, R.; Korshuize-van Straten, L.; de Jong, C.E.; Steggerda-Carvalho, E.; Geijsen, E.D.; Stalpers, L.J.A.; Crezee, J. Predictive value of simulated SAR and temperature for changes in measured temperature after phase-amplitude steering during locoregional hyperthermia treatments. Int. J. Hyperth. 2018, 35, 330–339. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Paulides, M.M.; Rodrigues, D.B.; Bellizzi, G.G.; Sumser, K.; Curto, S.; Neufeld, E.; Montanaro, H.; Kok, H.P.; Dobsicek Trefna, H. ESHO benchmarks for computational modeling and optimization in hyperthermia therapy. Int. J. Hyperth. 2021, 38, 1425–1442. [Google Scholar] [CrossRef] [PubMed]
  92. Bellizzi, G.G.; Drizdal, T.; van Rhoon, G.C.; Crocco, L.; Isernia, T.; Paulides, M.M. Predictive value of SAR based quality indicators for head and neck hyperthermia treatment quality. Int. J. Hyperth. 2019, 36, 456–465. [Google Scholar] [CrossRef] [Green Version]
  93. Canters, R.A.; Wust, P.; Bakker, J.F.; Van Rhoon, G.C. A literature survey on indicators for characterisation and optimisation of SAR distributions in deep hyperthermia, a plea for standardisation. Int. J. Hyperth. 2009, 25, 593–608. [Google Scholar] [CrossRef] [PubMed]
  94. Kok, H.P.; van der Zee, J.; Guirado, F.N.; Bakker, A.; Datta, N.R.; Abdel-Rahman, S.; Schmidt, M.; Wust, P.; Crezee, J. Treatment planning facilitates clinical decision making for hyperthermia treatments. Int. J. Hyperth. 2021, 38, 532–551. [Google Scholar] [CrossRef]
  95. Zanoli, M.; Trefna, H.D. Suitability of eigenvalue beam-forming for discrete multi-frequency hyperthermia treatment planning. Med. Phys. 2021, 48, 7410–7426. [Google Scholar] [CrossRef] [PubMed]
  96. Bruggmoser, G. Some aspects of quality management in deep regional hyperthermia. Int. J. Hyperth. 2012, 28, 562–569. [Google Scholar] [CrossRef]
  97. Jaffar, N.A.; Lias, K.B.; Madzhi, N.K.; Buniyamin, N. An overview of metamaterials used in applicators in hyperthermia cancer treatment procedure. In Proceedings of the 2017 International Conference on Electrical, Electronics and System Engineering (ICEESE), Kanazawa, Japan, 9–10 November 2017; pp. 32–36. [Google Scholar]
  98. Wang, G.; Gong, Y. Metamaterial lens applicator for microwave hyperthermia of breast cancer. Int. J. Hyperth. 2009, 25, 434–445. [Google Scholar] [CrossRef]
  99. Velazquez-Ahumada, M.C.; Freire, M.J.; Marques, R. Metamaterial applicator for microwave hyperthermia. In Proceedings of the 2011 XXXth URSI General Assembly and Scientific Symposium, Istanbul, Turkey, 13–20 August 2011; pp. 1–4. [Google Scholar]
  100. Schmidt, R.; Webb, A. Metamaterial Combining Electric- and Magnetic-Dipole-Based Configurations for Unique Dual-Band Signal Enhancement in Ultrahigh-Field Magnetic Resonance Imaging. ACS Appl. Mater. Interfaces 2017, 9, 34618–34624. [Google Scholar] [CrossRef]
  101. Webb, A.; Shchelokova, A.; Slobozhanyuk, A.; Zivkovic, I.; Schmidt, R. Novel materials in magnetic resonance imaging: High permittivity ceramics, metamaterials, metasurfaces and artificial dielectrics. MAGMA 2022, 35, 875–894. [Google Scholar] [CrossRef]
  102. Chen, H.; Guo, L.; Li, M.; Destruel, A.; Liu, C.; Weber, E.; Liu, F.; Crozier, S. Metamaterial-Inspired Radiofrequency (RF) Shield With Reduced Specific Absorption Rate (SAR) and Improved Transmit Efficiency for UHF MRI. IEEE Trans. Biomed. Eng. 2021, 68, 1178–1189. [Google Scholar] [CrossRef] [PubMed]
  103. Paulides, M.M.; Stauffer, P.R.; Neufeld, E.; Maccarini, P.F.; Kyriakou, A.; Canters, R.A.; Diederich, C.J.; Bakker, J.F.; Van Rhoon, G.C. Simulation techniques in hyperthermia treatment planning. Int. J. Hyperth. 2013, 29, 346–357. [Google Scholar] [CrossRef] [PubMed]
  104. Drizdal, T.; van Rhoon, G.C.; Verhaart, R.F.; Fiser, O.; Paulides, M.M. A Guide for Water Bolus Temperature Selection for Semi-Deep Head and Neck Hyperthermia Treatments Using the HYPERcollar3D Applicator. Cancers 2021, 13, 6126. [Google Scholar] [CrossRef] [PubMed]
  105. Katscher, U.; van den Berg, C.A.T. Electric properties tomography: Biochemical, physical and technical background, evaluation and clinical applications. NMR Biomed. 2017, 30, e3729. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Scheme of thermal MR treatment planning (TMTP). B1+ = transmission field used for MRI; TC = tumor coverage; T = temperature; SAR = specific absorption rate; SAF = SAR amplification factor; HTQ = hotspots-to-target quotient.
Figure 1. Scheme of thermal MR treatment planning (TMTP). B1+ = transmission field used for MRI; TC = tumor coverage; T = temperature; SAR = specific absorption rate; SAF = SAR amplification factor; HTQ = hotspots-to-target quotient.
Cancers 15 02303 g001
Figure 3. Overview of ThermalMR RF applicator configurations with 16LD, 8D and 8L variants around the head of Duke. Circular and elliptical array configurations were designed using 360° coverage, and the horse-shoe shape using arc = 270° coverage of the head. Arrangements of the loop+SGBT dipole building blocks are illustrated on the left column: (A) circular array (cir), (B) circular horse-shoe array (cir_HS), (C) elliptical array (ellip), (D) elliptical horse-shoe array (ellip_HS).
Figure 3. Overview of ThermalMR RF applicator configurations with 16LD, 8D and 8L variants around the head of Duke. Circular and elliptical array configurations were designed using 360° coverage, and the horse-shoe shape using arc = 270° coverage of the head. Arrangements of the loop+SGBT dipole building blocks are illustrated on the left column: (A) circular array (cir), (B) circular horse-shoe array (cir_HS), (C) elliptical array (ellip), (D) elliptical horse-shoe array (ellip_HS).
Cancers 15 02303 g003
Figure 4. Multiplexed vector field shaping (MVFS) optimization [84]. The optimization process is a convex formulation of the time-frequency-multiplexed constrained RF heating problem, and provides an iterative algorithm to efficiently find its globally optimum solution. The MVFS algorithm provides globally optimal excitation vectors by defining the appropriate phase and amplitude settings for each RF channel of the ThermalMR RF applicator, and automatically selects the appropriate intervention frequencies and time-interleaved excitations to match the constructive interference pattern to best deliver RF power at the desired target location and shape. The resulting SAR distribution of the incident E-field interference is tailored to focus heating of the target volume while minimizing local peak RF exposure to healthy tissue below a safe threshold. If M individual excitation modes belong to the same excitation frequency, then time multiplexing is done with each solution vector scaled by √M, and the excitations were played out sequentially [84]. Excitations at different frequencies can be played out concurrently as their electromagnetic fields do not interact coherently [84]. The final resultant target shape is created by superimposing these individual patterns for each frequency.
Figure 4. Multiplexed vector field shaping (MVFS) optimization [84]. The optimization process is a convex formulation of the time-frequency-multiplexed constrained RF heating problem, and provides an iterative algorithm to efficiently find its globally optimum solution. The MVFS algorithm provides globally optimal excitation vectors by defining the appropriate phase and amplitude settings for each RF channel of the ThermalMR RF applicator, and automatically selects the appropriate intervention frequencies and time-interleaved excitations to match the constructive interference pattern to best deliver RF power at the desired target location and shape. The resulting SAR distribution of the incident E-field interference is tailored to focus heating of the target volume while minimizing local peak RF exposure to healthy tissue below a safe threshold. If M individual excitation modes belong to the same excitation frequency, then time multiplexing is done with each solution vector scaled by √M, and the excitations were played out sequentially [84]. Excitations at different frequencies can be played out concurrently as their electromagnetic fields do not interact coherently [84]. The final resultant target shape is created by superimposing these individual patterns for each frequency.
Cancers 15 02303 g004
Figure 5. B1+ distribution maps (central axial slice of the head of the human voxel model Duke) of 16LD, 8D and 8L ThermalMR RF applicators cir, cir_HS, ellip, ellip_HS at magnetic field strength of (A) 7.0 T, (B) 9.4 T, (C) 10.5 T in circular polarization mode for a spherical region of interest (ROI) (radius = 2 cm, red circle) centered in the head of the human voxel model Duke. All four 16LD ThermalMR RF applicators show enhanced B1+ distribution inside the ROI compared with the corresponding 8D and 8L ThermalMR RF applicators. Comparison of the B1+ efficiency in the ROI for the 16LD, 8D and 8L ThermalMR RF applicators variants cir, cir_HS, ellip, ellip_HS at field strengths (D) 7.0 T, (E) 9.4 T, (F) 10.5 T shows increasing mean and maximum B1+ with increasing magnetic field strength, with reduced B1+ homogeneity. The 16LD variants had significantly higher B1+ compared to the 8D and 8L variants (see Table 2).
Figure 5. B1+ distribution maps (central axial slice of the head of the human voxel model Duke) of 16LD, 8D and 8L ThermalMR RF applicators cir, cir_HS, ellip, ellip_HS at magnetic field strength of (A) 7.0 T, (B) 9.4 T, (C) 10.5 T in circular polarization mode for a spherical region of interest (ROI) (radius = 2 cm, red circle) centered in the head of the human voxel model Duke. All four 16LD ThermalMR RF applicators show enhanced B1+ distribution inside the ROI compared with the corresponding 8D and 8L ThermalMR RF applicators. Comparison of the B1+ efficiency in the ROI for the 16LD, 8D and 8L ThermalMR RF applicators variants cir, cir_HS, ellip, ellip_HS at field strengths (D) 7.0 T, (E) 9.4 T, (F) 10.5 T shows increasing mean and maximum B1+ with increasing magnetic field strength, with reduced B1+ homogeneity. The 16LD variants had significantly higher B1+ compared to the 8D and 8L variants (see Table 2).
Cancers 15 02303 g005
Figure 6. Maximum local SAR10g does not exceed the IEC limits for all ThermalMR RF applicators for MR imaging. SAR10g maps of the head of the human voxel model Duke (central axial slice) of ThermalMR RF applicators 16LD, 8D and 8L with cir, cir_HS, ellip, ellip_HS variants at magnetic field strength (A) 7.0 T, (B) 9.4 T, (C) 10.5 T. Quantification of maximum local SAR10g at (D) 7.0 T, (E) 9.4 T, (F) 10.5 T.
Figure 6. Maximum local SAR10g does not exceed the IEC limits for all ThermalMR RF applicators for MR imaging. SAR10g maps of the head of the human voxel model Duke (central axial slice) of ThermalMR RF applicators 16LD, 8D and 8L with cir, cir_HS, ellip, ellip_HS variants at magnetic field strength (A) 7.0 T, (B) 9.4 T, (C) 10.5 T. Quantification of maximum local SAR10g at (D) 7.0 T, (E) 9.4 T, (F) 10.5 T.
Cancers 15 02303 g006
Figure 7. Optimization results obtained for targeted RF heating using the MVFS algorithm. Results are shown for frequencies of 250, 300, 350, 400, 450 MHz on the head of the human voxel model Duke (axial slice) for each ThermalMR RF applicator (16LD, 8D, 8L) with the four variants: (A) cir, (B) cir_HS, (C) ellip, (D) ellip_HS. Each panel shows a central axial slice through the tumor center; the targeted tumor region (TV) is depicted with a white line. The larger plots on the left of each panel shows the total achieved SAR10g in the tumor TV (the mean ± SD and range (min–max value) is indicated above each plot. The resultant target pattern was created from the individual contributing time and frequency-multiplexed modes, shown in the smaller plots on the right sides of each panel (scaled to their individual maxima), with the respective peak contribution inside the tumor TV indicated above.
Figure 7. Optimization results obtained for targeted RF heating using the MVFS algorithm. Results are shown for frequencies of 250, 300, 350, 400, 450 MHz on the head of the human voxel model Duke (axial slice) for each ThermalMR RF applicator (16LD, 8D, 8L) with the four variants: (A) cir, (B) cir_HS, (C) ellip, (D) ellip_HS. Each panel shows a central axial slice through the tumor center; the targeted tumor region (TV) is depicted with a white line. The larger plots on the left of each panel shows the total achieved SAR10g in the tumor TV (the mean ± SD and range (min–max value) is indicated above each plot. The resultant target pattern was created from the individual contributing time and frequency-multiplexed modes, shown in the smaller plots on the right sides of each panel (scaled to their individual maxima), with the respective peak contribution inside the tumor TV indicated above.
Cancers 15 02303 g007
Figure 8. SAR10g of targeted RF heating optimization and temperature distribution maps obtained for the 16LD ThermalMR RF applicators. (A) Diagram showing the placement of the 16LD arrays around the head. (B) Maps of total local SAR10g achieved within the tumor target volume (white line) from targeted RF heating optimization; sagittal, coronal, axial views. The maximum (max) and mean SAR10g in the target volume is indicated above the maps. (C) Temperature distribution maps based on total achieved SAR10g results obtained from targeted RF heating optimization. T = 37.4 °C was used as a baseline body temperature. Maximum (max) and mean temperature (°C) is indicated above the maps.
Figure 8. SAR10g of targeted RF heating optimization and temperature distribution maps obtained for the 16LD ThermalMR RF applicators. (A) Diagram showing the placement of the 16LD arrays around the head. (B) Maps of total local SAR10g achieved within the tumor target volume (white line) from targeted RF heating optimization; sagittal, coronal, axial views. The maximum (max) and mean SAR10g in the target volume is indicated above the maps. (C) Temperature distribution maps based on total achieved SAR10g results obtained from targeted RF heating optimization. T = 37.4 °C was used as a baseline body temperature. Maximum (max) and mean temperature (°C) is indicated above the maps.
Cancers 15 02303 g008
Figure 9. Maximum temperature (Tmax) obtained in tumor TV with respect to RF exposure time with deposited RF power (SAR10g) calculated from targeted RF heating optimization. For all 16LD ThermalMR RF applicators, Tmax of the tumor TV reached a steady state after ~20.5 min. The 16LDcir and 16LDellip yielded a steady state Tmax~41.7 °C and the horse-shoe 16LDcir_HS and 16LDellip_HS achieved a steady state Tmax~42.3 °C.
Figure 9. Maximum temperature (Tmax) obtained in tumor TV with respect to RF exposure time with deposited RF power (SAR10g) calculated from targeted RF heating optimization. For all 16LD ThermalMR RF applicators, Tmax of the tumor TV reached a steady state after ~20.5 min. The 16LDcir and 16LDellip yielded a steady state Tmax~41.7 °C and the horse-shoe 16LDcir_HS and 16LDellip_HS achieved a steady state Tmax~42.3 °C.
Cancers 15 02303 g009
Figure 10. Quantification of SAR10g of targeted RF heating optimization, SAR amplification factor (SAF), and hotspot-to-target quotient (HTQ). (A) The 16LD ThermalMR RF applicators showed significantly greater total SAR10g inside the tumor TV compared with the 8D and 8L RF applicators. The 16LDcir, 16LDcir_HS, 16LDellip, 16LDellip_HS ThermalMR RF applicators show a 24%, 12%, 14% and 34% increase of mean SAR10g over 8D variants. Statistical analysis shows for the cir_HS ThermalMR RF applicator variant that SAR10g was 34.7 ± 6.57 (median ± IQR) compared with 30.9 ± 8.47 for 8D, and 5.57 ± 5.89 for 8L. This result was consistent for the cir, ellip and ellip_HS (Table 3). (B) SAR amplification factor (SAF) presents the ratio of mean SAR10g deposited in the target volume (TV) versus the surrounding healthy tissue. A heatmap shows the SAF across all ThermalMR RF applicators (left). Color ranges from blue (cooler) to yellow (warmer); blue indicates better preservation of remote healthy tissue from SAR exposure. The 16LD and 8D ThermalMR RF applicators achieved higher SAF than the 8L ThermalMR RF applicators, indicating better RF power deposition (SAR10g) inside the TV and better sparing of remote healthy tissue (right). All 8L ThermalMR RF applicators showed SAF~2. Comparing among variants of the 16LD, 8D, 8L ThermalMR RF applicators show that both horse-shoe ThermalMR RF applicators achieved higher SAF than the 360° RF applicators (right). (C) Heatmap of hotspot-to-target quotient (HTQ) of all ThermalMR RF applicators showing that the ellip_HS variant had better HTQ (left). Color ranges from blue (cooler) to yellow (warmer); blue indicates better performance with regards to hotspots. Comparison among all ThermalMR RF applicators shows that both the 16LD and 8D ThermalMR RF applicator have lower values (~0.88) compared to the 8L (left).
Figure 10. Quantification of SAR10g of targeted RF heating optimization, SAR amplification factor (SAF), and hotspot-to-target quotient (HTQ). (A) The 16LD ThermalMR RF applicators showed significantly greater total SAR10g inside the tumor TV compared with the 8D and 8L RF applicators. The 16LDcir, 16LDcir_HS, 16LDellip, 16LDellip_HS ThermalMR RF applicators show a 24%, 12%, 14% and 34% increase of mean SAR10g over 8D variants. Statistical analysis shows for the cir_HS ThermalMR RF applicator variant that SAR10g was 34.7 ± 6.57 (median ± IQR) compared with 30.9 ± 8.47 for 8D, and 5.57 ± 5.89 for 8L. This result was consistent for the cir, ellip and ellip_HS (Table 3). (B) SAR amplification factor (SAF) presents the ratio of mean SAR10g deposited in the target volume (TV) versus the surrounding healthy tissue. A heatmap shows the SAF across all ThermalMR RF applicators (left). Color ranges from blue (cooler) to yellow (warmer); blue indicates better preservation of remote healthy tissue from SAR exposure. The 16LD and 8D ThermalMR RF applicators achieved higher SAF than the 8L ThermalMR RF applicators, indicating better RF power deposition (SAR10g) inside the TV and better sparing of remote healthy tissue (right). All 8L ThermalMR RF applicators showed SAF~2. Comparing among variants of the 16LD, 8D, 8L ThermalMR RF applicators show that both horse-shoe ThermalMR RF applicators achieved higher SAF than the 360° RF applicators (right). (C) Heatmap of hotspot-to-target quotient (HTQ) of all ThermalMR RF applicators showing that the ellip_HS variant had better HTQ (left). Color ranges from blue (cooler) to yellow (warmer); blue indicates better performance with regards to hotspots. Comparison among all ThermalMR RF applicators shows that both the 16LD and 8D ThermalMR RF applicator have lower values (~0.88) compared to the 8L (left).
Cancers 15 02303 g010
Figure 11. Magnified field of view (FoV) of RF power deposition (SAR10g) obtained from targeted RF heating optimization in the tumor target volume. The peak SAR10g within the tumor TV (white line) is shown for 16LD, 8D, 8L ThermalMR RF applicators for the (A) cir, (B) cir_HS, (C) ellip, (D) ellip_HS variants. Quantification of tumor coverage (TC) is shown in the spider plots on the right, which show the fraction of the tumor enclosed within the 25, 50, 80 and 100% isolines of peak SAR10g (TC25, TC50, TC80, TC100, respectively) for each ThermalMR RF applicator. The 16LD ThermalMR RF applicator achieved much greater TC50, TC80 and TC100 than the 8D. The 8L ThermalMR RF applicators failed to achieve TC25 = 100%, and in general had much lower TC values compared to 16LD and 8D ThermalMR RF applicators.
Figure 11. Magnified field of view (FoV) of RF power deposition (SAR10g) obtained from targeted RF heating optimization in the tumor target volume. The peak SAR10g within the tumor TV (white line) is shown for 16LD, 8D, 8L ThermalMR RF applicators for the (A) cir, (B) cir_HS, (C) ellip, (D) ellip_HS variants. Quantification of tumor coverage (TC) is shown in the spider plots on the right, which show the fraction of the tumor enclosed within the 25, 50, 80 and 100% isolines of peak SAR10g (TC25, TC50, TC80, TC100, respectively) for each ThermalMR RF applicator. The 16LD ThermalMR RF applicator achieved much greater TC50, TC80 and TC100 than the 8D. The 8L ThermalMR RF applicators failed to achieve TC25 = 100%, and in general had much lower TC values compared to 16LD and 8D ThermalMR RF applicators.
Cancers 15 02303 g011
Figure 12. Iso-contour temperature distributions maps of Duke’s head (central axial slice) based on SAR10g deposition from targeted RF heating. Top row of iso-contour temperature maps of 16LD ThermalMR RF applicators (A) 16LDcir, (B) 16LDcir_HS, (C) 16LDellip, (D) 16LDellip_HS, showing the temperature increase from 37.4 °C basal body temperature. The 16LDcir (A) and 16LDellip (C) ThermalMR RF applicators show unwanted hotspots of ~42 °C and ~43 °C in the right facial muscles. Correspondingly, magnified field of view (FoV) (bottom row) of highlighted tumor region showing the detail temperature coverage results of T40°C, T41°C, T42°C. Magnified FoV of the tumor region shows that the horse-shoe ThermalMR RF applicators yielded a higher temperature increase (B) 16LDcir_HS: ΔT = 4.9 °C, (D) 16LDellip_HS: ΔT = 4.8 °C from 37.4 °C basal body temperature than the (A) 16LDcir and (C) 16LDellip ThermalMR RF applicators, both of which failed to achieve any target volume with T42°C in contrast to the horse-shoe ThermalMR RF applicators (T42°C~20%).
Figure 12. Iso-contour temperature distributions maps of Duke’s head (central axial slice) based on SAR10g deposition from targeted RF heating. Top row of iso-contour temperature maps of 16LD ThermalMR RF applicators (A) 16LDcir, (B) 16LDcir_HS, (C) 16LDellip, (D) 16LDellip_HS, showing the temperature increase from 37.4 °C basal body temperature. The 16LDcir (A) and 16LDellip (C) ThermalMR RF applicators show unwanted hotspots of ~42 °C and ~43 °C in the right facial muscles. Correspondingly, magnified field of view (FoV) (bottom row) of highlighted tumor region showing the detail temperature coverage results of T40°C, T41°C, T42°C. Magnified FoV of the tumor region shows that the horse-shoe ThermalMR RF applicators yielded a higher temperature increase (B) 16LDcir_HS: ΔT = 4.9 °C, (D) 16LDellip_HS: ΔT = 4.8 °C from 37.4 °C basal body temperature than the (A) 16LDcir and (C) 16LDellip ThermalMR RF applicators, both of which failed to achieve any target volume with T42°C in contrast to the horse-shoe ThermalMR RF applicators (T42°C~20%).
Cancers 15 02303 g012
Figure 13. Overall performance evaluation of ThermalMR RF applicators 16LDcir, 16LDellip, 16LDcir_HS and 16LDellip_HS using the metrics maximum and mean SAR (SARmax(TV), SARmean(TV)) and temperature (Tmax, Tmean) inside tumor, tumor coverage (TC80, TC100), SAR amplification factor (SAF) and hotspot-to-target quotient (HTQ). Both horse-shoe ThermalMR RF applicators 16LDcir_HS and 16LDellip_HS showed superior TC80 and TC100, higher temperature rise, enhanced SAF and a lower HTQ, compared with the 360° ThermalMR RF applicators 16LDcir, 16LDellip.
Figure 13. Overall performance evaluation of ThermalMR RF applicators 16LDcir, 16LDellip, 16LDcir_HS and 16LDellip_HS using the metrics maximum and mean SAR (SARmax(TV), SARmean(TV)) and temperature (Tmax, Tmean) inside tumor, tumor coverage (TC80, TC100), SAR amplification factor (SAF) and hotspot-to-target quotient (HTQ). Both horse-shoe ThermalMR RF applicators 16LDcir_HS and 16LDellip_HS showed superior TC80 and TC100, higher temperature rise, enhanced SAF and a lower HTQ, compared with the 360° ThermalMR RF applicators 16LDcir, 16LDellip.
Cancers 15 02303 g013
Table 1. Input power for temperature simulation from individual excitation mode (M) power obtained from targeted RF heating optimization for 16LD thermalMR RF applicators.
Table 1. Input power for temperature simulation from individual excitation mode (M) power obtained from targeted RF heating optimization for 16LD thermalMR RF applicators.
ThermalMR
RF Applicator
Individual Excitation Mode (M) Power
16LDcir300 MHz350 MHz450 MHz450 MHzTotal
0.65 W31.48 W24.01 W15.12 W71.26 W
16LDcir_HS 350 MHz450 MHz450 MHzTotal
21.33 W34.21 W15.23 W70.77 W
16LDellip300 MHz350 MHz450 MHz450 MHzTotal
14.55 W15.84 W22.72 W18.21 W71.32 W
16LDellip_HS 350 MHz450 MHz450 MHzTotal
18.89 W32.66 W18.75 W70.3 W
Table 2. Summary of statistical analysis of B1+ transmission field for MRI performance of the ThermalMR RF applicators.
Table 2. Summary of statistical analysis of B1+ transmission field for MRI performance of the ThermalMR RF applicators.
ThermalMR
RF Applicator
Mean
7 T
Max
7 T
Min
7 T
* p-Value (vs. 16LD)
7 T
Mean
9.4 T
Max
9.4 T
Min
9.4 T
* p-Value (vs. 16LD) 9.4 TMean
10.5 T
Max
10.5 T
Min
10.5 T
* p-Value (vs. 16LD) 10.5 T
cir16LD19.4221.39 16.14 20.8724.74 15.26 21.2526.1813.90
8D12.8914.63 10.271.05 × 10−8316.1819.0710.838.42 × 10−7017.2621.3010.23 1.56 × 10−43
8L14.6316.26 12.281.09 × 10−8313.416.109.723.33 × 10−8312.6915.868.79 5.70 × 10−82
cir_HS16LD18.33 20.22 14.69 20.0124.0614.02 20.3324.7012.94
8D14.3616.03 10.80 3.23 × 10−8016.8220.1411.083.61 × 10−4616.8420.9610.16 9.48 × 10−34
8L11.5712.69 9.90 1.03 × 10−8311.8614.068.631.07 × 10−8311.5114.117.78 3.80 × 10−83
ellip16LD19.7721.95 16.66 20.0323.8214.80 20.4025.2013.53
8D12.9514.54 10.71 1.05 × 10−8315.9918.5211.30 9.25 × 10−6517.04 20.64 10.66 1.71 × 10−36
8L15.0716.78 12.171.31 × 10−8312.4115.228.342.04 × 10−8311.7015.03 7.16 2.34 × 10−82
ellip_HS16LD18.37 20.47 14.48 19.74 23.4413.85 20.6625.54 13.47
8D14.11 15.52 10.93 8.13 × 10−8216.4319.1810.99 1.31 × 10−4917.2921.20 10.61 3.08 × 10−34
8L11.89 13.83 8.72 1.05 × 10−8311.5113.997.441.11 × 10−8311.7714.92 7.21 4.94 × 10−82
* p-values are from comparison to 16LD.
Table 3. Summary of statistical analysis of the targeted RF heating performance of the ThermalMR RF applicators using the metric total SAR10g inside the tumor volume.
Table 3. Summary of statistical analysis of the targeted RF heating performance of the ThermalMR RF applicators using the metric total SAR10g inside the tumor volume.
ThermalMR
RF Applicator
MeanMaxMin* p-Value (vs. 16LD)
cir16LD34.140.422.9
8D27.138.812.45.32 × 10−5
8L10.824.21.333.96 × 10−8
cir_HS16LD34.740.420.3
8D30.940.317.36.64 × 10−5
8L5.5712.90.7121.18 × 10−7
ellip16LD33.640.620.5
8D25.035.410.55.85 × 10−5
8L9.1419.31.473.52 × 10−7
ellip_HS16LD34.440.920.3
8D30.140.617.21.9 × 10−3
8L7.6918.11.151.56 × 10−7
* p-values are from comparison to 16LD.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saha, N.; Kuehne, A.; Millward, J.M.; Eigentler, T.W.; Starke, L.; Waiczies, S.; Niendorf, T. Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors. Cancers 2023, 15, 2303. https://doi.org/10.3390/cancers15082303

AMA Style

Saha N, Kuehne A, Millward JM, Eigentler TW, Starke L, Waiczies S, Niendorf T. Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors. Cancers. 2023; 15(8):2303. https://doi.org/10.3390/cancers15082303

Chicago/Turabian Style

Saha, Nandita, Andre Kuehne, Jason M. Millward, Thomas Wilhelm Eigentler, Ludger Starke, Sonia Waiczies, and Thoralf Niendorf. 2023. "Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors" Cancers 15, no. 8: 2303. https://doi.org/10.3390/cancers15082303

APA Style

Saha, N., Kuehne, A., Millward, J. M., Eigentler, T. W., Starke, L., Waiczies, S., & Niendorf, T. (2023). Advanced Radio Frequency Applicators for Thermal Magnetic Resonance Theranostics of Brain Tumors. Cancers, 15(8), 2303. https://doi.org/10.3390/cancers15082303

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop