3D radiobiological evaluation of combined radiotherapy and hyperthermia treatments

A Novel Framework for Thermoradiotherapy Treatment Planning

PURPOSE

Thermoradiotherapy combines radiation therapy with hyperthermia to increase therapeutic effectiveness. Currently, both modalities are optimized separately and in state-of-the-art research the enhanced therapeutic effect is evaluated using equivalent radiation dose in 2-Gy fractions (EQD2). This study proposes a novel thermoradiotherapy treatment planning framework with voxelwise EQD2 radiation therapy optimizing including thermal radiosensitization and direct thermal cytotoxicity.

METHODS AND MATERIALS

To demonstrate proof-of-concept of the planning framework, 3 strategies consisting of 20 radiation therapy fractions were planned for 4 prostate cancer cases with substantially different temperature distributions: (1) Conventional radiation therapy plan of 60 Gy combined with 4 hyperthermia sessions (RT60 + HT), (2) standalone uniform dose escalation to 68 Gy without hyperthermia (RT68), and (3) uniform target EQD2 that maximizes the tumor control probability (TCP) accounting for voxelwise thermal effects of 4 hyperthermia sessions without increasing normal tissue doses (RTHT + HT). Assessment included dose, EQD2, TCP, and rectal normal tissue complication probability (NTCP), alongside robustness analyses for TCP and NTCP against parameter uncertainties.

RESULTS

The estimated TCP of around 76% for RT60 without hyperthermia was increased to an average of 85.9% (range, 81.3%-90.5%) for RT60 + HT, 92.5% (92.4%-92.5%) for RT68, and 94.4% (91.7%-96.6%) for RTHT + HT. The corresponding averaged rectal NTCPs were 8.7% (7.9%-10.0%), 14.9% (13.8%-17.1%), and 8.4% (7.5%-9.7%), respectively. RT68 and RTHT + HT exhibited slightly enhanced TCP robustness against parameter uncertainties compared with RT60 + HT, and RT68 presented higher and less robust rectal NTCP values compared with the other planning strategies.

CONCLUSIONS

This study introduces an innovative thermoradiotherapy planning approach, integrating thermal effects into EQD2-based radiation therapy optimization. Results demonstrate an ability to achieve enhanced and uniform target EQD2 and TCP across various temperature distributions without elevating normal tissue EQD2 or NTCP compared with conventional methods. Although promising for improving clinical outcomes, realizable enhancements depend on accurate tumor- and tissue-specific data and precise quantification of hyperthermic effects, which are seamlessly integrable in the planning framework as they emerge.

Magnetic-guided nanocarriers for ionizing/non-ionizing radiation synergistic treatment against triple-negative breast cancer

BACKGROUND

Triple-negative breast cancer (TNBC) is a subtype of breast cancer with the worst prognosis. Radiotherapy (RT) is one of the core modalities for the disease; however, the ionizing radiation of RT has severe side effects. The consistent development direction of RT is to achieve better therapeutic effect with lower radiation dose. Studies have demonstrated that synergistic effects can be achieved by combining RT with non-ionizing radiation therapies such as light and magnetic therapy, thereby achieving the goal of dose reduction and efficacy enhancement.

METHODS

In this study, we applied FeCo NPs with magneto thermal function and phototherapeutic agent IR-780 to construct an ionizing and non-ionizing radiation synergistic nanoparticle (INS NPs). INS NPs are first subjected to morphology, size, colloidal stability, loading capacity, and photothermal conversion tests. Subsequently, the cell inhibitory and cellular internalization were evaluated using cell lines in vitro. Following comprehensive assessment of the NPs' in vivo biocompatibility, tumor-bearing mouse model was established to evaluate their distribution, targeted delivery, and anti-tumor effects in vivo.

RESULTS

INS NPs have a saturation magnetization exceeding 72 emu/g, a hydrodynamic particle size of approximately 40 nm, a negatively charged surface, and good colloidal stability and encapsulation properties. INS NPs maintain the spectral characteristics of IR-780 at 808 nm. Under laser irradiation, the maximum temperature was 92 °C, INS NPs also achieved the effective heat temperature in vivo. Both in vivo and in vitro tests have proven that INS NPs have good biocompatibility. INS NPs remained effective for more than a week after one injection in vivo, and can also be guided and accumulated in tumors through permanent magnets. Later, the results exhibited that under low-dose RT and laser irradiation, the combined intervention group showed significant synergetic effects, and the ROS production rate was much higher than that of the RT and phototherapy-treated groups. In the mice model, 60% of the tumors were completely eradicated.

CONCLUSIONS

INS NPs effectively overcome many shortcomings of RT for TNBC and provide experimental basis for the development of novel clinical treatment methods for TNBC.

Quantitative analysis of radiosensitizing effect for magnetic hyperthermia-radiation combined therapy on prostate cancer cells

BACKGROUND

Magnetic hyperthermia (MHT) has emerged as a promising therapeutic approach in the field of radiation oncology due to its superior precision in controlling temperature and managing the heating area compared to conventional hyperthermia. Recent studies have proposed solutions to address clinical safety concerns associated with MHT, which arise from the use of highly concentrated magnetic nanoparticles and the strong magnetic field needed to induce hyperthermic effects. Despite these efforts, challenges remain in quantifying therapeutic outcomes and developing treatment plan systems for combining MHT with radiation therapy (RT).

PURPOSE

This study aims to quantitatively measure the therapeutic effect, including radiation dose enhancement (RDE) in the magnetic hyperthermia-radiation combined therapy (MHRT), using the equivalent radiation dose (EQD) estimation method.

METHODS

To conduct EQD estimation for MHRT, we compared the therapeutic effects between the conventional hyperthermia-radiation combined therapy (HTRT) and MHRT in human prostate cancer cell lines, PC3 and LNCaP. We adopted a clonogenic assay to validate RDE and the radiosensitizing effect induced by MHT. The data on survival fractions were analyzed using both the linear-quadradic model and Arrhenius model to estimate the biological parameters describing RDE and radiosensitizing effect of MHRT for both cell lines through maximum likelihood estimation. Based on these parameters, a new survival fraction model was suggested for EQD estimation of MHRT.

RESULTS

The newly designed model describing the MHRT effect, effectively captures the variations in thermal and radiation dose for both cell lines (R2 > 0.95), and its suitability was confirmed through the normality test of residuals. This model appropriately describes the survival fractions up to 10 Gy for PC3 cells and 8 Gy for LNCaP cells under RT-only conditions. Furthermore, using the newly defined parameter r, the RDE effect was calculated as 29% in PC3 cells and 23% in LNCaP cells. EQDMHRT calculated through this model was 9.47 Gy for PC3 and 4.71 Gy for LNCaP when given 2 Gy and MHT for 30 min. Compared to EQDHTRT, EQDMHRT showed a 26% increase for PC3 and a 20% increase for LNCaP.

CONCLUSIONS

The proposed model effectively describes the changes of the survival fraction induced by MHRT in both cell lines and adequately represents actual data values through residual analysis. Newly suggested parameter r for RDE effect shows potential for quantitative comparisons between HTRT and MHRT, and optimizing therapeutic outcomes in MHRT for prostate cancer.

Biological treatment evaluation in thermoradiotherapy: application in cervical cancer patients

BACKGROUND

Hyperthermia treatment quality is usually evaluated by thermal (dose) parameters, though hyperthermic radiosensitization effects are also influenced by the time interval between the two modalities. This work applies biological modelling for clinical treatment evaluation of cervical cancer patients treated with radiotherapy plus hyperthermia by calculating the equivalent radiation dose (EQDRT, i.e., the dose needed for the same effect with radiation alone). Subsequent analyses evaluate the impact of logistics.

METHODS

Biological treatment evaluation was performed for 58 patients treated with 23-28 fractions of 1.8-2 Gy plus 4-5 weekly hyperthermia sessions. Measured temperatures (T50) and recorded time intervals between the radiotherapy and hyperthermia sessions were used to calculate the EQDRT using an extended linear quadratic (LQ) model with hyperthermic LQ parameters based on extensive experimental data. Next, the impact of a 30-min time interval (optimized logistics) as well as a 4‑h time interval (suboptimal logistics) was evaluated.

RESULTS

Median average measured T50 and recorded time intervals were 41.2 °C (range 39.7-42.5 °C) and 79 min (range 34-125 min), respectively, resulting in a median total dose enhancement (D50) of 5.5 Gy (interquartile range [IQR] 4.0-6.6 Gy). For 30-min time intervals, the enhancement would increase by ~30% to 7.1 Gy (IQR 5.5-8.1 Gy; p < 0.001). In case of 4‑h time intervals, an ~ 40% decrease in dose enhancement could be expected: 3.2 Gy (IQR 2.3-3.8 Gy; p < 0.001). Normal tissue enhancement was negligible (< 0.3 Gy), even for short time intervals.

CONCLUSION

Biological treatment evaluation is a useful addition to standard thermal (dose) evaluation of hyperthermia treatments. Optimizing logistics to shorten time intervals seems worthwhile to improve treatment efficacy.

Theoretical evaluation of the impact of diverse treatment conditions by calculation of the tumor control probability (TCP) of simulated cervical cancer Hyperthermia-Radiotherapy (HT-RT) treatments in-silico

INTRODUCTION

Hyperthermia (HT) induces various cellular biological processes, such as repair impairment and direct HT cell killing. In this context, in-silico biophysical models that translate deviations in the treatment conditions into clinical outcome variations may be used to study the extent of such processes and their influence on combined hyperthermia plus radiotherapy (HT + RT) treatments under varying conditions.

METHODS

An extended linear-quadratic model calibrated for SiHa and HeLa cell lines (cervical cancer) was used to theoretically study the impact of varying HT treatment conditions on radiosensitization and direct HT cell killing effect. Simulated patients were generated to compute the Tumor Control Probability (TCP) under different HT conditions (number of HT sessions, temperature and time interval), which were randomly selected within margins based on reported patient data.

RESULTS

Under the studied conditions, model-based simulations suggested a treatment improvement with a total CEM43 thermal dose of approximately 10 min. Additionally, for a given thermal dose, TCP increased with the number of HT sessions. Furthermore, in the simulations, we showed that the TCP dependence on the temperature/time interval is more correlated with the mean value than with the minimum/maximum value and that comparing the treatment outcome with the mean temperature can be an excellent strategy for studying the time interval effect.

CONCLUSION

The use of thermoradiobiological models allows us to theoretically study the impact of varying thermal conditions on HT + RT treatment outcomes. This approach can be used to optimize HT treatments, design clinical trials, and interpret patient data.

Radiosensitization by Hyperthermia Critically Depends on the Time Interval

PURPOSE

Hyperthermia is a potent sensitizer of radiation therapy that improves both tumor control and survival in women with locally advanced cervical cancer (LACC). The optimal sequence and interval between hyperthermia and radiation therapy are still under debate.

METHODS AND MATERIALS

We investigated the interval and sequence in vitro in cervical cancer cell lines, patient-derived organoids, and SiHa cervical cancer hind leg xenografts in athymic nude mice and compared the results with retrospective results from 58 women with LACC treated with thermoradiotherapy.

RESULTS

All 3 approaches confirmed that shortening the interval between hyperthermia and radiation therapy enhanced hyperthermic radiosensitization by 2 to 8 times more DNA double-strand breaks and apoptosis and 10 to 100 times lower cell survival, delayed tumor growth in mice, and increased the 5-year survival rate of women with LACC from 22% (interval ≥80 minutes) to 54% (interval <80 minutes). In vitro and in vivo results showed that the sequence of hyperthermia and radiation therapy did not affect the outcome.

CONCLUSIONS

Shortening the interval between hyperthermia and radiation therapy significantly improves treatment outcomes. The sequence of hyperthermia and radiation therapy (before or after) does not seem to matter.

The Relevance of High Temperatures and Short Time Intervals Between Radiation Therapy and Hyperthermia: Insights in Terms of Predicted Equivalent Enhanced Radiation Dose

PURPOSE

The radiosensitization effect of hyperthermia can be considered and quantified as an enhanced equivalent radiation dose (EQD_RT), that is, the dose needed to achieve the same effect without hyperthermia. EQD_RT can be predicted using an extended linear quadratic model, with temperature-dependent parameters. Clinical data show that both the achieved temperature and time interval between radiation therapy and hyperthermia correlate with clinical outcome, but their effect on expected EQD_RT is unknown and was therefore evaluated in this study.

METHODS AND MATERIALS

Biological modeling was performed using our in-house developed software (X-Term), considering a 23- × 2-Gy external beam radiation scheme, as applied for patients with locally advanced cervical cancer. First, the EQD_RT was calculated for homogeneous temperature levels, evaluating time intervals between 0 and 4 hours. Next, realistic heterogeneous hyperthermia treatment plans were combined with radiation therapy plans and the EQD_RT was calculated for 10 patients. Furthermore, the effect of achieving 0.5°C to 1°C lower or higher temperatures was evaluated.

RESULTS

EQD_RT increases substantially with both increasing temperature and decreasing time interval. The effect of the time interval is most pronounced at higher temperatures (>41°C). At a typical hyperthermic temperature level of 41.5°C, an enhancement of ∼10 Gy can be realized with a 0-hour time interval, which is decreased to only ∼4 Gy enhancement with a 4-hour time interval. Most enhancement is already lost after 1 hour. Evaluation in patients predicted an average additional EQD_RT (D95%) of 2.2 and 6.3 Gy for 4- and 0-hour time intervals, respectively. The effect of 0.5°C to 1°C lower or higher temperatures is most pronounced at high temperature levels and short time intervals. The additional EQD_RT (D95%) ranged between 1.5 and 3.3 Gy and between 4.5 and 8.5 Gy for 4- and 0-hour time intervals, respectively.

CONCLUSIONS

Biological modeling provides relevant insight into the relationship between treatment parameters and expected EQD_RT. Both high temperatures and short time intervals are essential to maximize EQD_RT.

Biological modeling in thermoradiotherapy: present status and ongoing developments toward routine clinical use

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

Design and Characterization of an RF Applicator for In Vitro Tests of Electromagnetic Hyperthermia

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

A Novel Framework for the Optimization of Simultaneous ThermoBrachyTherapy

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

Radiobiological Evaluation of Combined Gamma Knife Radiosurgery and Hyperthermia for Pediatric Neuro-Oncology

Simple Summary

This study proposes a novel strategy in brain cancer management. Stereotactic radiosurgery delivered by the Gamma Knife was combined with hyperthermia. For the radiobiological modelling of this synergistic treatment modality, we used the linear-quadratic model with temperature-dependent parameters to assess the potential enhancement of the therapeutic outcome. The results indicate that focused intracranial heating can be used to boost the dose to the target. Alternatively, one can conclude that for the same therapeutic effect, hyperthermia can help to minimize the dose undesirably delivered to healthy tissues. This study is also the first to advocate a combination of stereotactic radiosurgery with focused heating and motivates the future development of hyperthermia systems for brain cancer treatment.

Abstract

Combining radiotherapy (RT) with hyperthermia (HT) has been proven effective in the treatment of a wide range of tumours, but the combination of externally delivered, focused heat and stereotactic radiosurgery has never been investigated. We explore the potential of such treatment enhancement via radiobiological modelling, specifically via the linear-quadratic (LQ) model adapted to thermoradiotherapy through modulating the radiosensitivity of temperature-dependent parameters. We extend this well-established model by incorporating oxygenation effects. To illustrate the methodology, we present a clinically relevant application in pediatric oncology, which is novel in two ways. First, it deals with medulloblastoma, the most common malignant brain tumour in children, a type of brain tumour not previously reported in the literature of thermoradiotherapy studies. Second, it makes use of the Gamma Knife for the radiotherapy part, thereby being the first of its kind in this context. Quantitative metrics like the biologically effective dose (BED) and the tumour control probability (TCP) are used to assess the efficacy of the combined plan.

Advanced patient-specific hyperthermia treatment planning

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

Mathematical model for the thermal enhancement of radiation response: thermodynamic approach

Radiotherapy can effectively kill malignant cells, but the doses required to cure cancer patients may inflict severe collateral damage to adjacent healthy tissues. Recent technological advances in the clinical application has revitalized hyperthermia treatment (HT) as an option to improve radiotherapy (RT) outcomes. Understanding the synergistic effect of simultaneous thermoradiotherapy via mathematical modelling is essential for treatment planning. We here propose a theoretical model in which the thermal enhancement ratio (TER) relates to the cell fraction being radiosensitised by the infliction of sublethal damage through HT. Further damage finally kills the cell or abrogates its proliferative capacity in a non-reversible process. We suggest the TER to be proportional to the energy invested in the sensitisation, which is modelled as a simple rate process. Assuming protein denaturation as the main driver of HT-induced sublethal damage and considering the temperature dependence of the heat capacity of cellular proteins, the sensitisation rates were found to depend exponentially on temperature; in agreement with previous empirical observations. Our findings point towards an improved definition of thermal dose in concordance with the thermodynamics of protein denaturation. Our predictions well reproduce experimental in vitro and in vivo data, explaining the thermal modulation of cellular radioresponse for simultaneous thermoradiotherapy.

Evaluating the Propagation of Uncertainties in Biologically Based Treatment Planning Parameters

Biologically based treatment planning is a broad term used to cover any instance in radiotherapy treatment planning where some form of biological input has been used. This is wide ranging, and the simpler forms (e.g., fractionation modification/optimization) have been in use for many years. However, there is a reluctance to use more sophisticated methods that incorporate biological models either for plan evaluation purposes or for driving plan optimizations. This is due to limited data available regarding the uncertainties in these model parameters and what impact these have clinically. This work aims to address some of these issues and to explore the role that uncertainties in individual model parameters have on the overall tumor control probability (TCP)/normal tissue complication probability (NTCP) calculated, those parameters that have the largest influence and situations where extra care must be taken. In order to achieve this, a software tool was developed, which can import individual clinical DVH's for analysis using a range of different TCP/NTCP models. On inputting individual model parameters, an uncertainty can be applied. Using a normally distributed random number generator, distributions of parameters can be generated, from which TCP/NTCP values can be calculated for each parameter set for the DVH in question. These represent the spread in TCP/NTCP parameters that would be observed for a simulated population of patients all being treated with that particular dose distribution. A selection of clinical DVHs was assessed using published parameters and their associated uncertainties. A range of studies was carried out to determine the impact of individual parameter uncertainties including reduction of uncertainties and assessment of what impact fractionation and dose have on these probabilities.

Holistic View on Cell Survival and DNA Damage: How Model-Based Data Analysis Supports Exploration of Dynamics in Biological Systems

In this work, a method is established to calibrate a model that describes the basic dynamics of DNA damage and repair. The model can be used to extend planning for radiotherapy and hyperthermia in order to include the biological effects. In contrast to “syntactic” models (e.g., describing molecular kinetics), the model used here describes radiobiological semantics, resulting in a more powerful model but also in a far more challenging calibration. Model calibration is attempted from clonogenic assay data (doses of 0–6 Gy) and from time-resolved comet assay data obtained within 6 h after irradiation with 6 Gy. It is demonstrated that either of those two sources of information alone is insufficient for successful model calibration, and that both sources of information combined in a holistic approach are necessary to find viable model parameters. Approximate Bayesian computation (ABC) with simulated annealing is used for parameter search, revealing two aspects that are beneficial to resolving the calibration problem: (1) assessing posterior parameter distributions instead of point-estimates and (2) combining calibration runs from different assays by joining posterior distributions instead of running a single calibration run with a combined, computationally very expensive objective function.

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

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

3D tumour spheroids for the prediction of the effects of radiation and hyperthermia treatments

For multimodality therapies such as the combination of hyperthermia and radiation, quantification of biological effects is key for dose prescription and response prediction. Tumour spheroids have a microenvironment that more closely resembles that of tumours in vivo and may thus be a superior in vitro cancer model than monolayer cultures. Here, the response of tumour spheroids formed from two established human cancer cell lines (HCT116 and CAL27) to single and combination treatments of radiation (0-20 Gy), and hyperthermia at 47 °C (0-780 CEM43) has been evaluated. Response was analysed in terms of spheroid growth, cell viability and the distribution of live/dead cells. Time-lapse imaging was used to evaluate mechanisms of cell death and cell detachment. It was found that sensitivity to heat in spheroids was significantly less than that seen in monolayer cultures. Spheroids showed different patterns of shrinkage and regrowth when exposed to heat or radiation: heated spheroids shed dead cells within four days of heating and displayed faster growth post-exposure than samples that received radiation or no treatment. Irradiated spheroids maintained a dense structure and exhibited a longer growth delay than spheroids receiving hyperthermia or combination treatment at (thermal) doses that yielded equivalent levels of clonogenic cell survival. We suggest that, unlike radiation, which kills dividing cells, hyperthermia-induced cell death affects cells independent of their proliferation status. This induces microenvironmental changes that promote spheroid growth. In conclusion, 3D tumour spheroid growth studies reveal differences in response to heat and/or radiation that were not apparent in 2D clonogenic assays but that may significantly influence treatment efficacy.

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

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

Hyperthermia Treatment Planning Including Convective Flow in Cerebrospinal Fluid for Brain Tumour Hyperthermia Treatment Using a Novel Dedicated Paediatric Brain Applicator

Hyperthermia therapy (40-44 °C) is a promising option to increase efficacy of radiotherapy/chemotherapy for brain tumours, in particular paediatric brain tumours. The Chalmers Hyperthermia Helmet is developed for this purpose. Hyperthermia treatment planning is required for treatment optimisation, but current planning systems do not involve a physically correct model of cerebrospinal fluid (CSF). This study investigates the necessity of fluid modelling for treatment planning. We made treatments plans using the Helmet for both pre-operative and post-operative cases, comparing temperature distributions predicted with three CSF models: a convective "fluid" model, a non-convective "solid" CSF model, and CSF models with increased effective thermal conductivity ("high-k"). Treatment plans were evaluated by T90, T50 and T10 target temperatures and treatment-limiting hot spots. Adequate heating is possible with the helmet. In the pre-operative case, treatment plan quality was comparable for all three models. In the post-operative case, the high-k models were more accurate than the solid model. Predictions to within ±1 °C were obtained by a 10-20-fold increased effective thermal conductivity. Accurate modelling of the temperature in CSF requires fluid dynamics, but modelling CSF as a solid with enhanced effective thermal conductivity might be a practical alternative for a convective fluid model for many applications.

Quantitative Estimation of the Equivalent Radiation Dose Escalation using Radiofrequency Hyperthermia in Mouse Xenograft Models of Human Lung Cancer

Hyperthermia is a potent radiosensitizer, and its effect varies according to the different types of cancer cells. In the present study, the radiosensitizing effect of hyperthermia on lung cancer cell lines A549 and NCI-H1299 was determined based on the equivalent radiation dose escalation. In vitro cell experiments were conducted using lung cancer cell lines A549 and NCI-H1299 to determine thermal radiosensitivity. In vivo experiments were conducted using mouse heterotopic xenograft models to determine the treatment response and increase in the temperature of tumors using a 13.56 MHz radiofrequency (RF) hyperthermia device. Using the α and β values of the linear–quadratic equations of cell survival curves, numerical simulations were performed to calculate the equivalent radiation dose escalations. The dielectric properties of tumors were measured, and their effect on the calculated equivalent radiation dose was analyzed. Hyperthermia increased the equivalent radiation dose of lung cancer xenografts and a higher escalation was found in NCI-H1299 cells compared with that observed in A549 cells. An underestimation of the calculated equivalent radiation dose was observed when the dielectric property of the tumor was varied. This study may contribute to the effective planning of thermoradiotherapy in clinics.

Re‑irradiation plus hyperthermia for recurrent pediatric sarcoma; a simulation study to investigate feasibility

Recurrent pediatric tumors pose a challenge since treatment options may be limited, particularly after previous irradiation. Positive results have been reported for chemotherapy and hyperthermia, but the combination of re‑irradiation and hyperthermia has not been investigated thus far, although it is a proven treatment strategy in adults. The theoretical feasibility of re‑irradiation plus hyperthermia was investigated for infield recurrent pediatric sarcoma in the pelvic region and the extremities. A total of 46 recurrent pediatric sarcoma cases diagnosed at the Academic Medical Center (Amsterdam, The Netherlands) between 2002 and 2017 were evaluated. Patients not previously irradiated, outfield recurrences and locations other than the pelvis and extremities were excluded, ultimately yielding four eligible patients: Two with sarcomas in the pelvis and two in an extremity. Re‑irradiation and hyperthermia treatment plans were simulated for 23x2 Gy treatment schedules and weekly hyperthermia. The radiosensitizing effect of hyperthermia was quantified using biological modelling with a temperature‑dependent change in the parameters of the linear‑quadratic model. The possible effectiveness of re‑irradiation plus hyperthermia was estimated by calculating the equivalent radiotherapy dose distribution. Treatment planning revealed that tumors located in the pelvis and the extremities can be effectively heated in children. Equivalent dose distributions indicated that hyperthermic radiosensitization can be quantified as a target‑selective additional D95% of typically 10 Gy, thereby delivering a possibly curative dose of 54 Gy, without substantially increasing the equivalent dose to the organs at risk. Therefore, re‑irradiation plus hyperthermia is a theoretically feasible and possibly effective treatment option for recurrent pediatric sarcoma in the pelvic region and the extremities, and its clinical feasibility is worthy of evaluation.

A comprehensive model for heat-induced radio-sensitisation

Combined radiotherapy (RT) and hyperthermia (HT) treatments may improve treatment outcome by heat induced radio-sensitisation. We propose an empirical cell survival model (AlphaR model) to describe this multimodality therapy. The model is motivated by the observation that heat induced radio-sensitisation may be explained by a reduction in the DNA damage repair capacity of heated cells. We assume that this repair is only possible up to a threshold level above which survival will decrease exponentially with dose. Experimental cell survival data from two cell lines (HCT116, Cal27) were considered along with that taken from the literature (baby hamster kidney [BHK] and Chinese hamster ovary cells [CHO]) for HT and combined RT-HT. The AlphaR model was used to study the dependence of clonogenic survival on treatment temperature, and thermal dose R2 ≥ 0.95 for all fits). For HT survival curves (0-80 CEM43 at 43.5-57 °C), the number of free fit AlphaR model parameters could be reduced to two. Both parameters increased exponentially with temperature. We derived the relative biological effectiveness (RBE) or HT treatments at different temperatures, to provide an alternative description of thermal dose, based on our AlphaR model. For combined RT-HT, our analysis is restricted to the linear quadratic arm of the model. We show that, for the range used (20-80 CEM43, 0-12 Gy), thermal dose is a valid indicator of heat induced radio-sensitisation, and that the model parameters can be described as a function thereof. Overall, the proposed model provides a flexible framework for describing cell survival curves, and may contribute to better quantification of heat induced radio-sensitisation, and thermal dose in general.

The alfa and beta of tumours: a review of parameters of the linear-quadratic model, derived from clinical radiotherapy studies

BACKGROUND

Prediction of radiobiological response is a major challenge in radiotherapy. Of several radiobiological models, the linear-quadratic (LQ) model has been best validated by experimental and clinical data. Clinically, the LQ model is mainly used to estimate equivalent radiotherapy schedules (e.g. calculate the equivalent dose in 2 Gy fractions, EQD₂), but increasingly also to predict tumour control probability (TCP) and normal tissue complication probability (NTCP) using logistic models. The selection of accurate LQ parameters α, β and α/β is pivotal for a reliable estimate of radiation response. The aim of this review is to provide an overview of published values for the LQ parameters of human tumours as a guideline for radiation oncologists and radiation researchers to select appropriate radiobiological parameter values for LQ modelling in clinical radiotherapy.

METHODS AND MATERIALS

We performed a systematic literature search and found sixty-four clinical studies reporting α, β and α/β for tumours. Tumour site, histology, stage, number of patients, type of LQ model, radiation type, TCP model, clinical endpoint and radiobiological parameter estimates were extracted. Next, we stratified by tumour site and by tumour histology. Study heterogeneity was expressed by the I² statistic, i.e. the percentage of variance in reported values not explained by chance.

RESULTS

A large heterogeneity in LQ parameters was found within and between studies (I² > 75%). For the same tumour site, differences in histology partially explain differences in the LQ parameters: epithelial tumours have higher α/β values than adenocarcinomas. For tumour sites with different histologies, such as in oesophageal cancer, the α/β estimates correlate well with histology. However, many other factors contribute to the study heterogeneity of LQ parameters, e.g. tumour stage, type of LQ model, TCP model and clinical endpoint (i.e. survival, tumour control and biochemical control).

CONCLUSIONS

The value of LQ parameters for tumours as published in clinical radiotherapy studies depends on many clinical and methodological factors. Therefore, for clinical use of the LQ model, LQ parameters for tumour should be selected carefully, based on tumour site, histology and the applied LQ model. To account for uncertainties in LQ parameter estimates, exploring a range of values is recommended.

The alfa and beta of tumours: a review of parameters of the linear-quadratic model, derived from clinical radiotherapy studies

Background

Prediction of radiobiological response is a major challenge in radiotherapy. Of several radiobiological models, the linear-quadratic (LQ) model has been best validated by experimental and clinical data. Clinically, the LQ model is mainly used to estimate equivalent radiotherapy schedules (e.g. calculate the equivalent dose in 2 Gy fractions, EQD₂), but increasingly also to predict tumour control probability (TCP) and normal tissue complication probability (NTCP) using logistic models. The selection of accurate LQ parameters α, β and α/β is pivotal for a reliable estimate of radiation response. The aim of this review is to provide an overview of published values for the LQ parameters of human tumours as a guideline for radiation oncologists and radiation researchers to select appropriate radiobiological parameter values for LQ modelling in clinical radiotherapy.

Methods and materials

We performed a systematic literature search and found sixty-four clinical studies reporting α, β and α/β for tumours. Tumour site, histology, stage, number of patients, type of LQ model, radiation type, TCP model, clinical endpoint and radiobiological parameter estimates were extracted. Next, we stratified by tumour site and by tumour histology. Study heterogeneity was expressed by the I² statistic, i.e. the percentage of variance in reported values not explained by chance.

Results

A large heterogeneity in LQ parameters was found within and between studies (I² > 75%). For the same tumour site, differences in histology partially explain differences in the LQ parameters: epithelial tumours have higher α/β values than adenocarcinomas. For tumour sites with different histologies, such as in oesophageal cancer, the α/β estimates correlate well with histology. However, many other factors contribute to the study heterogeneity of LQ parameters, e.g. tumour stage, type of LQ model, TCP model and clinical endpoint (i.e. survival, tumour control and biochemical control).

Conclusions

The value of LQ parameters for tumours as published in clinical radiotherapy studies depends on many clinical and methodological factors. Therefore, for clinical use of the LQ model, LQ parameters for tumour should be selected carefully, based on tumour site, histology and the applied LQ model. To account for uncertainties in LQ parameter estimates, exploring a range of values is recommended.

Electronic supplementary material

The online version of this article (10.1186/s13014-018-1040-z) contains supplementary material, which is available to authorized users.

The effect of time interval between radiotherapy and hyperthermia on planned equivalent radiation dose

PURPOSE

Thermoradiotherapy is an effective treatment for locally advanced cervical cancer. However, the optimal time interval between radiotherapy and hyperthermia, resulting in the highest therapeutic gain, remains unclear. This study aims to evaluate the effect of time interval on the therapeutic gain using biological treatment planning.

METHODS

Radiotherapy and hyperthermia treatment plans were created for 15 cervical cancer patients. Biological modeling was used to calculate the equivalent radiation dose, that is, the radiation dose that results in the same biological effect as the thermoradiotherapy treatment, for different time intervals ranging from 0-4 h. Subsequently, the thermal enhancement ratio (TER, i.e. the ratio of the dose for the thermoradiotherapy and the radiotherapy-only plan) was calculated for the gross tumor volume (GTV) and the organs at risk (OARs: bladder, rectum, bowel), for each time interval. Finally, the therapeutic gain factor (TGF, i.e. TER_GTV/TER_OAR) was calculated for each OAR.

RESULTS

The median TER_GTV ranged from 1.05 to 1.16 for 4 h and 0 h time interval, respectively. Similarly, for bladder, rectum and bowel, TER_OARs ranged from 1-1.03, 1-1.04 and 1-1.03, respectively. Radiosensitization in the OARs was much less than in the GTV, because temperatures were lower, fractionation sensitivity was higher (lower α/β) and direct cytotoxicity was assumed negligible in normal tissue. TGFs for the three OARs were similar, and were highest (around 1.12) at 0 h time interval.

CONCLUSION

This planning study indicates that the largest therapeutic gain for thermoradiotherapy in cervical cancer patients can be obtained when hyperthermia is delivered immediately before or after radiotherapy.

Combining radiation with hyperthermia: a multiscale model informed by in vitro experiments

Combined radiotherapy and hyperthermia offer great potential for the successful treatment of radio-resistant tumours through thermo-radiosensitization. Tumour response heterogeneity, due to intrinsic, or micro-environmentally induced factors, may greatly influence treatment outcome, but is difficult to account for using traditional treatment planning approaches. Systems oncology simulation, using mathematical models designed to predict tumour growth and treatment response, provides a powerful tool for analysis and optimization of combined treatments. We present a framework that simulates such combination treatments on a cellular level. This multiscale hybrid cellular automaton simulates large cell populations (up to 107 cells) in vitro, while allowing individual cell-cycle progression, and treatment response by modelling radiation-induced mitotic cell death, and immediate cell kill in response to heating. Based on a calibration using a number of experimental growth, cell cycle and survival datasets for HCT116 cells, model predictions agreed well (R2 > 0.95) with experimental data within the range of (thermal and radiation) doses tested (0-40 CEM43, 0-5 Gy). The proposed framework offers flexibility for modelling multimodality treatment combinations in different scenarios. It may therefore provide an important step towards the modelling of personalized therapies using a virtual patient tumour.

Design and characterisation of a phased antenna array for intact breast hyperthermia

PURPOSE

Currently available hyperthermia technology is not well suited to treating cancer malignancies in the intact breast. This study investigates a microwave applicator incorporating multiple patch antennas, with the goal of facilitating controllable power deposition profiles for treating lesions at diverse locations within the intact breast.

MATERIALS AND METHODS

A 3D-computational model was implemented to assess power deposition profiles with 915 MHz applicators incorporating a hemispheric groundplane and configurations of 2, 4, 8, 12, 16 and 20 antennas. Hemispheric breast models of 90 mm and 150 mm diameter were considered, where cuboid target volumes of 10 mm edge length (1 cm³) and 30 mm edge length (27 cm³) were positioned at the centre of the breast, and also located 15 mm from the chest wall. The average power absorption (αPA) ratio expressed as the ratio of the PA in the target volume and in the full breast was evaluated. A 4-antenna proof-of-concept array was fabricated and experimentally evaluated.

RESULTS

Computational models identified an optimal inter-antenna spacing of 22.5° along the applicator circumference. Applicators with 8 and 12 antennas excited with constant phase presented the highest αPA at centrally located and deep-seated targets, respectively. Experimental measurements with a 4-antenna proof-of-concept array illustrated the potential for electrically steering power deposition profiles by adjusting the relative phase of the signal at antenna inputs.

CONCLUSIONS

Computational models and experimental results suggest that the proposed applicator may have potential for delivering conformal thermal therapy in the intact breast.

Measurement and analysis of the impact of time-interval, temperature and radiation dose on tumour cell survival and its application in thermoradiotherapy plan evaluation

PURPOSE

Biological modelling of thermoradiotherapy may further improve patient selection and treatment plan optimisation, but requires a model that describes the biological effect as a function of variables that affect treatment outcome (e.g. temperature, radiation dose). This study aimed to establish such a model and its parameters. Additionally, a clinical example was presented to illustrate the application.

METHODS

Cell survival assays were performed at various combinations of radiation dose (0-8 Gy), temperature (37-42 °C), time interval (0-4 h) and treatment sequence (radiotherapy before/after hyperthermia) for two cervical cancer cell lines (SiHa and HeLa). An extended linear-quadratic model was fitted to the data using maximum likelihood estimation. As an example application, a thermoradiotherapy plan (23 × 2 Gy + weekly hyperthermia) was compared with a radiotherapy-only plan (23 × 2 Gy) for a cervical cancer patient. The equivalent uniform radiation dose (EUD) in the tumour, including confidence intervals, was estimated using the SiHa parameters. Additionally, the difference in tumour control probability (TCP) was estimated.

RESULTS

Our model described the dependency of cell survival on dose, temperature and time interval well for both SiHa and HeLa data (R²=0.90 and R²=0.91, respectively), making it suitable for biological modelling. In the patient example, the thermoradiotherapy plan showed an increase in EUD of 9.8 Gy that was robust (95% CI: 7.7-14.3 Gy) against propagation of the uncertainty in radiobiological parameters. This corresponded to a 20% (95% CI: 15-29%) increase in TCP.

CONCLUSIONS

This study presents a model that describes the cell survival as a function of radiation dose, temperature and time interval, which is essential for biological modelling of thermoradiotherapy treatments.