Skip to main content

Biological modelling of the radiation dose escalation effect of regional hyperthermia in cervical cancer



Locoregional hyperthermia combined with radiotherapy significantly improves locoregional control and overall survival for cervical tumors compared to radiotherapy alone. In this study biological modelling is applied to quantify the effect of radiosensitization for three cervical cancer patients to evaluate the improvement in equivalent dose for the combination treatment with radiotherapy and hyperthermia.


The Linear-Quadratic (LQ) model extended with temperature-dependent LQ-parameters α and β was used to model radiosensitization by hyperthermia and to calculate the conventional radiation dose that is equivalent in biological effect to the combined radiotherapy and hyperthermia treatment. External beam radiotherapy planning was performed based on a prescription dose of 46Gy in 23 fractions of 2Gy. Hyperthermia treatment using the AMC-4 system was simulated based on the actual optimized system settings used during treatment.


The simulated hyperthermia treatments for the 3 patients yielded a T50 of 40.1 °C, 40.5 °C, 41.1 °C and a T90 of 39.2 °C, 39.7 °C, 40.4 °C, respectively. The combined radiotherapy and hyperthermia treatment resulted in a D95 of 52.5Gy, 55.5Gy, 56.9Gy in the GTV, a dose escalation of 7.3–11.9Gy compared to radiotherapy alone (D95 = 45.0–45.5Gy).


This study applied biological modelling to evaluate radiosensitization by hyperthermia as a radiation-dose escalation for cervical cancer patients. This model is very useful to compare the effectiveness of different treatment schedules for combined radiotherapy and hyperthermia treatments and to guide the design of clinical studies on dose escalation using hyperthermia in a multi-modality setting.


Achieving locoregional control is a challenge for locally advanced cervix tumors (FIGO stage IIB-IV). Currently combination of radiotherapy with concurrent cisplatin-based chemotherapy is widely accepted as a standard treatment [15]. An alternative is combining radiotherapy with hyperthermia. Hyperthermia (i.e. increasing tumor temperatures to 40–45 °C) is a very powerful radio and chemosensitizer [69] improving clinical outcome [1012]. Several randomized clinical trials have demonstrated that locoregional hyperthermia combined with radiotherapy (thermoradiotherapy) is at least as effective in improving locoregional control for cervical tumors as chemoradiotherapy [1317]. The Dutch Deep Hyperthermia Trial not only showed an 83 % local control rate for radiotherapy + hyperthermia versus 57 % for radiotherapy alone, but also a 3-year overall survival of 51 % for radiotherapy combined with hyperthermia, compared to 27 % with radiotherapy alone [13]. After 12 years survival is still almost twice as high after treatment with thermoradiotherapy (37 %), versus 20 % with radiotherapy [18]. A major advantage of hyperthermia given sequentially before or after radiotherapy is that the incidence of late toxicity grade 3 or higher is not significantly different from the toxicity reported after treatment with radiotherapy alone. Thermoradiotherapy is therefore often used as an alternative for the large number of stage IIB-IV patients unable to receive the standard therapy of concurrent chemotherapy and radiotherapy, e.g. because of kidney failure. Local tumor control achieved for standard therapy is good, particularly if combined with image guided brachytherapy [1922], the challenge is still to improve control of regional and distant disease and the reduction of side effects. Hyperthermia may be effective in achieving these objectives.

Planning for locoregional hyperthermia has become a standard clinical procedure [2325] but does not take the synergy with radiotherapy into account. The radiosensitizing effect of hyperthermia can be considered as a local increase in tumor dose, which can also be quantified using biological modelling [26, 27]. Hyperthermic radiosensitization can be modelled as a change in the α and ß parameters of the Linear Quadratic model [2830]. Recently a method has been presented to estimate the effect of a hyperthermia treatment in terms of equivalent dose distributions, i.e. the radiation dose that has a biological effect equivalent to that of the combined hyperthermia plus radiotherapy treatment. A conservative estimate of the temperature-dependency of the modelling parameters was applied due to lack of published temperature-dependent radiosensitivity data [31].

In this study biological modelling was applied for three cervical cancer patients treated with radiotherapy (external beam + brachytherapy) and hyperthermia to evaluate the improvement in equivalent dose for this combined treatment. The radiotherapy dose distribution was matched onto the geometry of the hyperthermia data set. Temperature-dependent parameters describing the radiosensitivity were estimated from in vitro experiments with a human cervical cancer cell line [29]. Hyperthermia treatments were simulated to evaluate the difference in equivalent dose. Predicted temperature distributions were compared to the measured tumor temperatures.


The workflow used consists of first computation of the hyperthermia and radiotherapy dose distributions followed by matching the radiotherapy dose distribution onto the geometry of the hyperthermia CT scan and finally computation of the equivalent dose distribution (Fig. 1).

Fig. 1
figure 1

Workflow. The workflow used consists of first computation of the hyperthermia and radiotherapy dose distributions, followed by matching the radiotherapy dose distribution onto the geometry of the hyperthermia CT scan and finally computation of the equivalent dose distribution

Hyperthermia treatment planning

Hyperthermia treatment planning [25] is a standard clinical procedure at the AMC and is based on a CT-scan in hyperthermia treatment position, i.e. in supine position on a water bolus and mattresses. Treatment planning was performed using a non-commercial software package developed at the AMC. Hounsfield Unit based segmentation was applied to segment the CT scan into muscle, fat, bone and air [32], to which literature-based dielectric and thermal tissue properties were assigned (Table 1). The tumor was delineated manually by a radiation oncologist. Hyperthermia treatment was simulated for the AMC-4 locoregional heating device. Electric field distributions in the patient were calculated at a resolution of 2.5x2.5x2.5 mm3 by solving Maxwell’s equations with the Finite Difference Time Domain method [33]. Antenna settings yielding an optimal steady-state tumor temperature distribution were determined using temperature-based optimization [34]. Heat-transfer computations for perfused tissues were based on the Pennes bio heat model [35]. In clinical hyperthermia tumor temperatures are usually reported as T10, T50 and T90, i.e. the temperature at least achieved in 10, 50 and 90 % of the target volume, respectively [36, 37]. A tumor temperature of 43 °C was the objective for optimization, normal tissue temperatures were constrained to 45 °C.

Table 1 Dielectric and thermal tissue properties used in hyperthermia treatment planning

Subsequently five weekly treatment sessions were given during the period external beam irradiation was given, one hour after a radiotherapy session. Hyperthermia was performed with the AMC-4 system consisting of four 70 MHz waveguides around the pelvis of the patient, with water bags cooled to 13 °C between waveguide and skin provide coupling of energy and skin cooling [38]. Thermometry consisted of thermocouple thermometry probes in the bladder, rectum and in a vaginal pelotte for temperature measurements during hyperthermia, yielding measurement points bordering the GTV. Phase settings during treatment are the pre-planned settings, modified after performing phase sweeps at the start of treatment [36], followed by a final phase optimization using three delta T pulses to achieve preferential heating of the cervix compared to bladder and rectum [39]. The simulated tumor temperatures for these phase and amplitude settings were used for biological modelling and were also compared with temperatures measured during the first hyperthermia session. The first session was chosen since the patient anatomy during this session best matches the anatomy as determined during CT and are thus expected to yield more reliable planning results than for later sessions.

Radiotherapy treatment planning

Radiotherapy treatment for locally advanced cervical cancer consists of a combination of external beam and MRI-guided brachytherapy. The prescription dose for external beam irradiation was 46Gy in 23 fractions of 2Gy. External beam radiotherapy planning was performed using Oncentra and based on a CT scan in prone position on a belly board. External beam radiotherapy was followed by a 24Gy brachytherapy boost using PDR brachytherapy. Brachytherapy planning is based on an MRI scan with the applicator in situ. The hyperthermia treatment is unlikely to enhance the brachytherapy dose since hyperthermia is only given during external beam radiotherapy (EBRT) at the AMC and since there is a gap of one week between the EBRT and brachytherapy treatments. Thus we performed our evaluation for the combination of hyperthermia and EBRT without involving the brachytherapy boost.

Image registration/matching

After treatment planning for both modalities, the next step is matching the radiotherapy dose distribution onto the geometry of the hyperthermia data set. To account for organ displacement and deformation due to different patient positioning (prone/supine) in the radiotherapy and hyperthermia CT scans, deformable matching software of Velocity Medical Solutions (Varian Medical systems, Palo Alto) was used (Fig. 2). First, the radiotherapy and hyperthermia CT scans were rigidly matched by visual assessment of the bony anatomy. Next, a deformable registration was made between the EBRT CT and the hyperthermia CT using an intensity-based deformable image registration algorithm. As this study focusses on the effect of hyperthermia on the tumor, accuracy of the deformable match was assessed by warping contours of organs located close to the GTV from the EBRT CT to the hyperthermia CT and verifying that their location was sufficiently accurate. If not, these organs were delineated on the hyperthermia scan as well and a new deformable image registration was performed (starting again from the rigid match), this time using a hybrid algorithm combining both intensity and structure based matching. The resulting deformation vector field was used to warp the radiotherapy dose distribution to the frame of reference of the hyperthermia CT.

Fig. 2
figure 2

Image registration/matching. Radiotherapy and hyperthermia CT scans are rigidly matched by visual assessment of the bony anatomy, followed by deformable registration using intensity-based deformable image registration software of Velocity Medical Solutions (Varian Medical systems, Palo Alto). Left: overlay showing the excellent match of the RT-CT onto the HT-CT in the Region of Interest (ROI) indicated with the dotted rectangle. Right: Outlines of rectum, cervix, bladder and part of the bony anatomy are shown in the RT-CT and the HT-CT

Biological modelling

To quantify the therapeutic effect of radiosensitization by hyperthermia, the radiotherapy dose distributions combined with hyperthermia were converted to equivalent radiation dose distributions (without hyperthermia in 2Gy/fraction (EQD2)). Calculation of equivalent radiation dose was based on the LQ-model [27], extended with temperature-dependent LQ-parameters α(T) and β(T) derived for cervical cancer cells in the hyperthermic temperature range as is described in the next paragraph.


Franken et al [29] reported α and β values for a human cervical cancer cell line (SiHa) at 37, 41 and 43 °C listed in Table 2. These values were used to define the temperature-dependent LQ-parameters α(T) and β(T) for the entire hyperthermic range needed for the LQ-model. The most significant change observed between 37 and 41 °C is an increase in β, whereas α increases strongly between 41 and 43 °C. We therefore applied a piecewise-linear interpolation for α(T) and β(T), consisting of one linear segment valid between 37 and 41 °C and another linear segment valid between 41 and 43 °C. To evaluate the impact of using realistic temperature-dependent LQ-parameters the outcome predicted by the radiobiological model using this piecewise-linear interpolation was then compared to the outcome using a linear interpolation using only the data for α and β at 37 and 43 °C.

Table 2 LQ-parameters α and β used in the equivalent dose calculations

Equivalent dose calculations

Equivalent dose distributions for combined radiotherapy-hyperthermia treatments were calculated using the LQ-model [27]. The expression used for EQD2 calculation is:

$$ EQ{D}_2=\frac{\alpha (T)\cdot D+\beta (T)\cdot d\cdot D}{\alpha (37)+\beta (37)\cdot 2} $$

with fraction dose d and total dose D and using the temperature-dependent LQ-parameters α and β as described in the previous section on LQ-parameters. Most of the procedure is similar to the procedure used in the simulation study of Kok et al. for prostate cancer [31]. Hyperthermic radiosensitization is assumed to be tumor selective since radiotherapy and hyperthermia are given sequentially with a time interval of at least 1 h [4042]. The temperature-dependent LQ-parameters α and β were thus only applied within the Gross Tumor Volume (GTV), elsewhere α and β were assumed to remain unchanged.


Temperature and Dose Volume Histogram

The temperature volume histograms indicating the tumor temperature distribution are shown in Fig. 3. The corresponding T10, T50 and T90 values are listed in Table 3 as well as the temperature measured during the first treatment session. The effect of adding hyperthermia is shown in the cumulative Dose Volume Histograms (DVH) in Fig. 4. These are computed assuming either a linear or a piecewise-linear temperature dependency for α(T) and β(T), yielding a remarkable difference in predicted dose escalation with a much more pronounced effect for the linear temperature dependency: D95 ranges between 68.4Gy and 82.9Gy for the linear interpolation compared to 52.5–56.9Gy for the piecewise linear interpolation and 45.0-45.5Gy for radiotherapy alone (Table 3). The piecewise linear interpolation results in a bending point in the DVH, separating the sections of the curve representing temperatures lower or higher than 41 °C. Temperatures higher than 41 °C result in far more significant dose escalation than temperatures below 41 °C.

Fig. 3
figure 3

Temperature Volume Histogram. Temperature Volume Histogram (TVH) representing the simulated temperature distribution within the GTV for patient 1, 2 and 3

Table 3 Radiotherapy and hyperthermia dose distribution
Fig. 4
figure 4

Dose Volume Histogram. Dose Volume Histogram reflecting the radiotherapy dose distribution within the GTV for patient 1, 2 and 3 comparing 3 different cases: Radiotherapy alone (RT only), radiotherapy and hyperthermia using a linear interpolation (RT + HT lin) or a piecewise linear interpolation (RT + HT pw lin) to the α and β values from Franken et al. [29]. EBRT only, brachytherapy boost not taken into account

Equivalent Dose distributions

Radiotherapy dose distributions for patient 1 with and without hyperthermia are shown in Fig. 5. Temperature sensitive LQ-parameters α(T) and β(T) were applied to the GTV using the piecewise-linear interpolation to the LQ-data listed in Table 2. The radiotherapy isodose curves for the treatment of radiotherapy alone are therefore nearly identical to the equivalent radiotherapy isodose curves accounting for the combination of radiotherapy and hyperthermia, with exception for the region inside the outlined GTV where a clear dose escalation is evident. The image with the hyperthermia temperature distribution overlaid clearly shows that the gradual increase in temperature from the dorsal to the ventral side of the GTV is associated with a matching gradual increase in effective radiation dose.

Fig. 5
figure 5

Radiotherapy and hyperthermia dose distributions. Radiotherapy isodose curves (top left), equivalent radiotherapy isodose curves for radiotherapy + hyperthermia (top right), idem with the hyperthermia temperature distribution overlaid as a color wash (bottom right) for patient #1 with cervical cancer. Radiotherapy isodose curves overlaid on the CT scan of the patient made in hyperthermia position lying on a hyperthermia water bolus. Gross tumour volume (GTV) is indicated by the bold red contour. The contribution of hyperthermia to equivalent radiotherapy isodose is visible within the GTV and increases with increasing temperature in ventral direction. Vaginal pelotte for temperature measurements during hyperthermia is visible adjacent to GTV


Effect of Hyperthermia on Radiotherapy dose distribution

This study shows a substantial dose escalation of 7-11Gy when radiotherapy is combined with hyperthermia (Table 3 and Fig. 4). Simulated temperatures match fairly well with the measured temperatures. Radiosensitization is not uniformly distributed, but follows the temperature distribution of the hyperthermia treatment (Fig. 5). As a result, the effective radiotherapy dose is no longer homogeneous within the GTV. The absence of any effect on normal tissue outside the GTV is a result of our assumption that hyperthermic radiosensitization is tumor selective when radiotherapy and hyperthermia are given sequentially with a time interval of 1 h or longer. This assumption is based on preclinical and clinical data. Overgaard determined the effect of sequence and time interval between radiotherapy and hyperthermia on the thermal enhancement ratio (TER) in tumor and normal skin in an in vivo tumor model. He found that the TER in skin equaled that of tumor when radiotherapy and hyperthermia are given simultaneously, but TER in skin quickly decreases to 1 when the time interval is 1 h or longer as currently applied in the clinic [41]. Clinical data for cervical cancer confirm that hyperthermia does not enhance radiation associated toxicity when the time interval between hyperthermia and radiotherapy exceeds 1 h [13].

Temperature dependence of radiosensitization

The temperature dependence of the LQ-parameters accounted for the temperature dependence of radiosensitization. The large difference between the DVHs based on the linear versus the piecewise-linear interpolation for the temperature dependent α(T) and β(T) values can be explained by overestimation of the effect of moderate temperature elevations below 41 °C on the LQ-parameters when using a linear interpolation. This overestimation of the effectiveness of hyperthermia results in a very high computed effective dose escalation that is not consistent with the known clinical results for hyperthermia from randomized trials. The dose escalation computed for the piecewise-linear interpolation is more realistic and closer to clinical reality. The piecewise linear interpolation results in a dose escalation that increments much stronger for temperatures above 41 °C, which would imply that temperatures exceeding 41 °C should be pursued in clinical hyperthermia. This demonstrates the need to establish more exact data on the temperature dependent LQ-parameters for the entire clinically relevant hyperthermic temperature range (39–43 °C).

The thermal dose effect relationship used in our simulations is strong and was based on cell survival curves of cervical cancer cell lines [29]. Clinical results for cervical cancer confirm there is a correlation between thermal dose and clinical results. Dinges et al. reported a correlation between local tumor control and CEM43T90 (cumulative equivalent minutes of T90 above 43 °C) for the 4 hyperthermia treatment sessions in a group of 18 patients with advanced carcinomas of the uterine cervix treated with the combination of radiotherapy and hyperthermia (p = 0.019) [43]. Franckena et al. analysed treatment outcome for 420 patients with locally advanced cervical cancer treated with radiotherapy and hyperthermia and found a significant correlation between tumor control and survival and CEM43T90 [44]. Franckena et al. found an even more impressive correlation between tumor control and survival and TRISE, a thermal dose parameter based on the product of the median temperature rise and the duration of heating [44].

Future perspectives

The temperature dependent LQ parameters used are based on in vitro data and represent mainly the direct radiosensitizing effect of hyperthermia on DNA repair pathways [8]. A next step is inclusion of other relevant factors and hyperthermia mechanisms like tissue reoxygenation, direct tumor cell kill, tumor hypoxia and repopulation. Most of these factors can be modelled as changes in the LQ parameters. A recent review discusses how these factors and mechanisms can be incorporated [45].

We did not combine brachytherapy with hyperthermia as combination with the individual 24 PDR pulses given in a period of 24 h is technically impossible. The synergy of hyperthermia given before or after the boost would be limited to a small fraction of the boost. It was therefore sufficient to determine the dose escalation of the EBRT only to estimate the dose escalation of the full radiation treatment due to hyperthermia. Inclusion of the brachytherapy boost in the simulation is logical when an HDR brachytherapy boost is given and combined with hyperthermia, although this is not straightforward. Matching the tissue deformations of the brachytherapy planning MRI images with the external beam radiotherapy CT planning scans is a very challenging task and will require a precise form of deformable image matching. Because of the steep gradients in the brachy dose distribution small errors in geometry due to the matching procedure would lead to relatively large deviations in the brachy dose distribution reconstructed in the EBRT scan. The most robust solution minimizing the risk of introducing errors might be to match the EBRT and hyperthermia datasets onto the brachytherapy image set as the dose gradients in the latter surpass dose gradients in the other sets.

An important step is to move from evaluation to optimization of the combined treatment. The present model evaluated the combined effect of hyperthermia and radiotherapy, where both modalities were individually optimized without taking advantage of the synergistic and complementary effects. In the second step the combined treatment could be optimized. This may result in higher hyperthermia doses being used close to organs at risk with regard to radiotherapy.

The change in Tumor Control Probability (TCP) as a result of adding hyperthermia would be an interesting parameter to analyse. Computation of TCP was left out of the present analysis as the brachytherapy given was not included in the current workflow. For tumors treated by EBRT only, conventional TCP models for radiotherapy [4648] can be used to calculate TCP based on the equivalent radiation dose. Okunieff et al [47] estimate that TCP increases with 1–3 %/Gy near TCD50 based on cervical cancer treatment data of Perez et al [49, 50]. This corresponds to a 10–30 % increase in TCP due to hyperthermia based on the median dose escalation of 10Gy predicted for our 3 patients (Table 3). This number is in line with the clinical data of cervical cancer trials, e.g. the 26 % increase in tumor control by adding hyperthermia to radiotherapy reported by van der Zee et al. [13].

A very useful application of this model is to compare the effectiveness of different treatment schedules and to guide the design of clinical studies on dose escalation using hyperthermia in a multi-modality setting.


This study demonstrated that radiosensitization by hyperthermia can be converted to radiation-dose escalation for cervical cancer patients using biological modelling.


  1. Keys HM, Bundy BN, Stehman FB, Muderspach LI, Chafe WE, Suggs 3rd CL, et al. Cisplatin, radiation, and adjuvant hysterectomy compared with radiation and adjuvant hysterectomy for bulky stage IB cervical carcinoma. N Engl J Med. 1999;340:1154–61.

    Article  CAS  PubMed  Google Scholar 

  2. Morris M, Eifel PJ, Lu J, Grigsby PW, Levenback C, Stevens RE, et al. Pelvic radiation with concurrent chemotherapy compared with pelvic and para-aortic radiation for high-risk cervical cancer. N Engl J Med. 1999;340:1137–43.

    Article  CAS  PubMed  Google Scholar 

  3. Peters III WA, Liu PY, Barrett II RJ, Stock RJ, Monk BJ, Berek JS, et al. Concurrent chemotherapy and pelvic radiation therapy compared with pelvic radiation therapy alone as adjuvant therapy after radical surgery in high-risk early-stage cancer of the cervix. J Clin Oncol. 2000;18:1606–13.

    CAS  PubMed  Google Scholar 

  4. Whitney CW, Sause W, Bundy BN, Malfetano JH, Hannigan EV, Fowler Jr WC, et al. Randomized comparison of fluorouracil plus cisplatin versus hydroxyurea as an adjunct to radiation therapy in stage IIB-IVA carcinoma of the cervix with negative para-aortic lymph nodes: a Gynecologic Oncology Group and Southwest Oncology Group study. J Clin Oncol. 1999;17:1339–48.

    CAS  PubMed  Google Scholar 

  5. Li Z, Yang S, Liu L, Han S. A comparison of concurrent chemoradiotherapy and radiotherapy in Chinese patients with locally advanced cervical carcinoma: a multi-center study. Radiat Oncol. 2014;9:212.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Horsman MR, Overgaard J. Hyperthermia: a potent enhancer of radiotherapy. Clin Oncol (R Coll Radiol). 2007;19:418–26.

    Article  CAS  Google Scholar 

  7. Kampinga HH. Cell biological effects of hyperthermia alone or combined with radiation or drugs: a short introduction to newcomers in the field. Int J Hyperthermia. 2006;22:191–6.

    Article  CAS  PubMed  Google Scholar 

  8. Oei AL, Vriend LE, Crezee J, Franken NA, Krawczyk PM. Effects of hyperthermia on DNA repair pathways: one treatment to inhibit them all. Radiat Oncol. 2015;10:165.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Oei AL, van Leeuwen CM, ten Cate R, Rodermond H, Buist MR, Stalpers LJ, et al. Hyperthermia Selectively Targets Human Papillomavirus in Cervical Tumors via p53-Dependent Apoptosis. Cancer Res. 2015;75:5120–9.

    Article  CAS  PubMed  Google Scholar 

  10. Wust P, Hildebrandt B, Sreenivasa G, Rau B, Gellermann J, Riess H, et al. Hyperthermia in combined treatment of cancer. Lancet Oncol. 2002;3:487–97.

    Article  CAS  PubMed  Google Scholar 

  11. Datta NR, Gomez Ordonez S, Gaipl US, Paulides MM, Crezee H, Gellermann J, et al. Local hyperthermia combined with radiotherapy and-/ or chemotherapy: Recent advances and promises for the future. Cancer Treat Rev. 2015;41:742–53.

    Article  CAS  PubMed  Google Scholar 

  12. Cihoric N, Tsikkinis A, van Rhoon G, Crezee H, Aebersold DM, Bodis S, et al. Hyperthermia-related clinical trials on cancer treatment within the registry. Int J Hyperthermia. 2015;31:609–14.

    Article  PubMed  Google Scholar 

  13. van der Zee J, Gonzalez Gonzalez D, van Rhoon GC, van Dijk JD, van Putten WL, Hart AA. Comparison of radiotherapy alone with radiotherapy plus hyperthermia in locally advanced pelvic tumours: a prospective, randomised, multicentre trial. Dutch Deep Hyperthermia Group. Lancet. 2000;355:1119–25.

    Article  PubMed  Google Scholar 

  14. Datta NR, Bose AK, Kapoor HK. Thermoradiotherapy in the management of carcinoma cervix (stage IIIB): A controlled clinical study. Indian Med Gaz. 1987;121:68–71.

    Google Scholar 

  15. Sharma S, Patel FD, Sandhu APS, Gupta BD, Yadav NS. A prospective randomised study of local hyperthermia as a supplement and radiosensitizer in the treatment of carcinoma of the cervix with radiotherapy. Endocuriether/Hyperthermia Oncol. 1989;5:151–9.

    Google Scholar 

  16. Chen HW, Jun-Jie F, Wei L. A randomized trial of hyper-thermo-radiochemotherapy for uterine cervix cancer. Chin J Clin Oncol. 1997;24:249–51.

    Google Scholar 

  17. Harima Y, Nagata K, Harima K, Ostapenko VV, Tanaka Y, Sawada S. A randomized clinical trial of radiation therapy versus thermoradiotherapy in stage IIIB cervical carcinoma. Int J Hyperthermia. 2001;17:97–105.

    Article  CAS  PubMed  Google Scholar 

  18. Franckena M, Stalpers LJ, Koper PC, Wiggenraad RG, Hoogenraad WJ, van Dijk JD, et al. Long-term improvement in treatment outcome after radiotherapy and hyperthermia in locoregionally advanced cervix cancer: an update of the Dutch Deep Hyperthermia Trial. Int J Radiat Oncol Biol Phys. 2008;70:1176–82.

    Article  PubMed  Google Scholar 

  19. Schmid MP, Kirisits C, Nesvacil N, Dimopoulos JC, Berger D, Pötter R. Local recurrences in cervical cancer patients in the setting of image-guided brachytherapy: a comparison of spatial dose distribution within a matched-pair analysis. Radiother Oncol. 2011;100:468–72.

    Article  PubMed Central  PubMed  Google Scholar 

  20. Nomden CN, de Leeuw AA, Roesink JM, Tersteeg RJ, Moerland MA, Witteveen PO, et al. Clinical outcome and dosimetric parameters of chemo-radiation including MRI guided adaptive brachytherapy with tandem-ovoid applicators for cervical cancer patients: a single institution experience. Radiother Oncol. 2013;107:69–74.

    Article  PubMed  Google Scholar 

  21. Rijkmans EC, Nout RA, Rutten IH, Ketelaars M, Neelis KJ, Laman MS, et al. Improved survival of patients with cervical cancer treated with image-guided brachytherapy compared with conventional brachytherapy. Gynecol Oncol. 2014;135:231–8.

    Article  CAS  PubMed  Google Scholar 

  22. Mazeron R, Maroun P, Castelnau-Marchand P, Dumas I, Del Campo ER, Cao K, et al. Pulsed-dose rate image-guided adaptive brachytherapy in cervical cancer: Dose-volume effect relationships for the rectum and bladder. Radiother Oncol 2015. [at press]

  23. Kok HP, Gellermann J, Van den Berg CA, Stauffer PR, Hand JW, Crezee J. Thermal modelling using discrete vasculature for thermal therapy: a review. Int J Hyperthermia. 2013;29:336–45.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Paulides MM, Stauffer PR, Neufeld E, Maccarini PF, Kyriakou A, Canters RA, et al. Simulation techniques in hyperthermia treatment planning. Int J Hyperthermia. 2013;29:346–57.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Kok HP, Wust P, Stauffer PR, Bardati F, van Rhoon GC, Crezee J. Current state of the art of regional hyperthermia treatment planning: a review. Radiat Oncol. 2015;10:196.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Barendsen GW. Dose fractionation, dose rate and iso-effect relationships for normal tissue responses. Int J Radiat Oncol Biol Phys. 1982;8:1981–97.

    Article  CAS  PubMed  Google Scholar 

  27. Fowler JF. The linear-quadratic formula and progress in fractionated radiotherapy. Br J Radiol. 1989;62:679–94.

    Article  CAS  PubMed  Google Scholar 

  28. Xu M, Myerson RJ, Straube WL, Moros EG, Lagroye I, Wang LL. Radiosensitization of heat resistant human tumour cells by 1 h at 41.1° C and its effect on DNA repair. Int J Hyperthermia. 2002;18:385–403.

    Article  CAS  PubMed  Google Scholar 

  29. Franken NA, Oei AL, Kok HP, Rodermond HM, Sminia P, Crezee J, et al. Cell survival and radiosensitisation: modulation of the linear and quadratic parameters of the LQ model (Review). Int J Oncol. 2013;42:1501–15.

    CAS  PubMed  Google Scholar 

  30. Myerson RJ, Roti Roti JL, Moros EG, Straube WL, Xu M. Modelling heat-induced radiosensitization: clinical implications. Int J Hyperthermia. 2004;20:201–12.

    Article  CAS  PubMed  Google Scholar 

  31. Kok HP, Crezee J, Franken NA, Stalpers LJ, Barendsen GW, Bel A. Quantifying the Combined Effect of Radiation Therapy and Hyperthermia in Terms of Equivalent Dose Distributions. Int J Radiat Oncol Biol Phys. 2014;88:739–45.

    Article  PubMed  Google Scholar 

  32. Hornsleth SN, Mella O, Dahl O. A new segmentation algorithm for finite difference based treatment planning systems. In: Franconi C, Arcangeli G, Cavaliere R, editors. Hyperthermic Oncology 1996 vol. 2. Rome, Italy Tor Vergata: 1996. p. p. 521-3.

  33. Taflove A. The Finite-Difference Time-Domain Method. (Boston, USA: Artech House): 1995.

  34. Das SK, Clegg ST, Samulski TV. Computational techniques for fast hyperthermia temperature optimization. Med Phys. 1999;26:319–28.

    Article  CAS  PubMed  Google Scholar 

  35. Pennes HH. Analysis of tissue and arterial blood temperatures in the resting human forearm. J Appl Physiol. 1948;1:93–122.

    CAS  PubMed  Google Scholar 

  36. Bruggmoser G, Bauchowitz S, Canters R, Crezee H, Ehmann M, Gellermann J, et al. Quality Assurance for Clinical Studies in Regional Deep Hyperthermia. Strahlenther Onkol. 2011;187:605–10.

    Article  PubMed  Google Scholar 

  37. van Rhoon GC. Why high quality hyperthermia is important, lessons to be learned. Radiat Oncol 2015. [at press]

  38. Crezee J, van Haaren PM, Westendorp H, De Greef M, Kok HP, Wiersma J, et al. Improving locoregional hyperthermia delivery using the 3-D controlled AMC-8 phased array hyperthermia system: a preclinical study. Int J Hyperthermia. 2009;25:581–92.

    Article  CAS  PubMed  Google Scholar 

  39. Kok HP, Ciampa S, de Kroon-Oldenhof R, Steggerda-Carvalho EJ3, van Stam G3, Zum Vörde Sive Vörding PJ, et al. Towards on-line adaptive hyperthermia treatment planning: correlation between measured and simulated SAR changes caused by phase steering in patients. Int J Radiat Oncol Biol Phys. 2014;90:438–45.

    Article  PubMed  Google Scholar 

  40. Field SB, Bleehen NM. Hyperthermia in the treatment of cancer. Cancer Treat Rev. 1979;6:63–94.

    Article  CAS  PubMed  Google Scholar 

  41. Overgaard J. Simultaneous and sequential hyperthermia and radiation treatment of an experimental tumor and its surrounding normal tissue in vivo. Int J Radiat Oncol Biol Phys. 1980;6:1507–17.

    Article  CAS  PubMed  Google Scholar 

  42. Myers R, Field SB. The response of the rat tail to combined heat and x rays. Br J Radiol. 1977;50:581–6.

    Article  CAS  PubMed  Google Scholar 

  43. Dinges S, Harder C, Wurm R, Buchali A, Blohmer J, Gellermann J, et al. Combined treatment of inoperable carcinomas of the uterine cervix with radiotherapy and regional hyperthermia. Results of a phase II trial. Strahlenther Onkol. 1998;174:517–21.

    Article  CAS  PubMed  Google Scholar 

  44. Franckena M, Fatehi D, de Bruijne M, Canters RA, van Norden Y, Mens JW, et al. Hyperthermia dose-effect relationship in 420 patients with cervical cancer treated with combined radiotherapy and hyperthermia. Eur J Cancer. 2009;45:1969–78.

    Article  PubMed  Google Scholar 

  45. Crezee H, van Leeuwen CM, Oei AL, Stalpers LJ, Bel A, Franken NA, Kok HP. Thermoradiotherapy planning: integration in routine clinical practice. Int J Hyperthermia 2016 (in press)

  46. Nahum AE, Tait DM. Maximizing local control by customized dose prescription for pelvic tumors: tumor response modelling and treatment planning. In: Advanced Radiation Therapy: Tumor Response Monitoring and Treatment Planning, ed A Breit (Berlin: Springer) 1992, p 425-31

  47. Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation dose-response of human tumors. Int J Radiat Oncol Biol Phys. 1995;32:1227–37.

    Article  CAS  PubMed  Google Scholar 

  48. Gay HA, Niemierko A. A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys Med. 2007;23:115–25.

    Article  PubMed  Google Scholar 

  49. Perez CA, Kao MS. Radiation therapy alone or combined with surgery in the treatment of barrel-shaped carcinoma of the uterine cervix (stages IB, IIA, IIB). Int J Radiat Oncol Biol Phys. 1985;11:1903–9.

    Article  CAS  PubMed  Google Scholar 

  50. Perez CA, Breaux S, Madoc-Jones H, Bedwinek JM, Camel HM, Purdy JA, et al. Radiation therapy alone in the treatment of carcinoma of uterine cervix. I. Analysis of tumor recurrence. Cancer. 1983;51:1393–402.

    Article  CAS  PubMed  Google Scholar 

Download references


This work was financially supported by the Dutch Cancer Society KWF, grant UVA 2012-5393 and UVA 2012-5540.

Author information

Authors and Affiliations


Corresponding author

Correspondence to J. Crezee.

Additional information

Competing interests

The authors have no competing interests to declare.

Authors’ contributions

JC analyzed the literature and carried out the writing of the manuscript. ALO and NAPF provided biological parameters. AB, HPK and CML wrote software. CML performed the planning and prepared the figures. PG assisted in the design of the manuscript. HPK provided general supervision and assisted in the design of the manuscript. All authors participated in the drafting and revising of the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Crezee, J., van Leeuwen, C.M., Oei, A.L. et al. Biological modelling of the radiation dose escalation effect of regional hyperthermia in cervical cancer. Radiat Oncol 11, 14 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: