Impact of time-related factors on biologically accurate radiotherapy treatment planning
Radiation Oncology volume 13, Article number: 30 (2018)
The incomplete repair (IR) model expresses the cell repair effect from radiation-induced damage over time, which is given little consideration in actual treatment planning. By incorporating the IR model into the normal tissue complication probability (NTCP), the accuracy and safety of treatment plan evaluations concerning the effect of repair can be improved. This study aims to evaluate the impact of incorporating the IR model into the NTCP by varying time-related factors such as the repair half-time (T1/2) and the junction-shift sc3hedule in craniospinal irradiation (CSI).
CSI was planned retrospectively, and the NTCP of the spinal cord was calculated with the IR model for values of T1/2 from 1 to 10 h. The NTCP in the case of changing the junction-shift schedule was also examined in the same manner.
The NTCP with the IR model increased with increasing T1/2, which is prominent for the larger T1/2. By changing the junction-shift schedule, the NTCP with the IR model decreased when adjacent fields overlapped.
The IR model is a valuable addition to treatment planning because it enables the NTCP to be evaluated including the effect of repair and differences in scheduling to be reflected in the NTCP. However, these are largely dependent on the value of the T1/2.
Currently, physical dosimetric information such as the dose–volume histogram and dose distribution are used to evaluate treatment planning. In addition, radiobiological treatment planning has great potential to estimate the treatment outcome . Although the normal tissue complication probability (NTCP) is a radiobiological index that can mathematically express the probability of complications, it overlooks several biological effects such as cell repair.
To account for the cell repair effect, Oliver et al.  proposed the time-based incomplete repair (IR) model, which was later modified by Levin-Plotnik et al. , to include successive doses. By incorporating this model into the NTCP formula, the probability of complications, including the cell repair effect, can be estimated. Many authors have studied the cell repair effect, particularly in the spinal cord. Several studies have considered the repair half-time, an important parameter of the IR model, by investigating the kinetics of repair in a rat’s spinal cord [4,5,6]. However, the reported values of repair half-time are wide ranging and contain uncertainties. For example, a repair half-time of 5.0 h with a wide confidence interval (CI) of 0.6–9.3 h has been reported . Moreover, these data were for rats, not humans. In the present study, we focus on the fact that the repair half-time of the spinal cord is comparatively long and contains uncertainties. We also examine the difference in the NTCP of the spinal cord with the IR model for varying values of repair half-time.
A form of radiotherapy that can potentially induce complications with the spinal cord is craniospinal irradiation (CSI), an essential component in the curative treatment of patients with primitive neuroectodermal tumors such as medulloblastoma and other brain tumors that have an increased risk of leptomeningeal spread [7, 8]. Owing to the planning target volume length in the craniocaudal direction, several fields are generally used to cover the whole target volume. The junction-shift technique is used to smooth any over- or underdose arising at a junction between different fields, which results from small errors in patient set up. However, the junction-shift schedule varies among facilities (e.g., change junction every week, two weeks, or one third of treatment.) , and is not covered by any recommendations. Considering that the variations in this schedule have no influence on the dose distribution, the conventional NTCP is calculated using a fixed physical dose distribution. However, we adopt the time-sensitive IR model in the expectation that the NTCP can depend on the junction-shift schedule.
In this study, we performed a planning study to examine the impact on NTCP of the factors related to cell repair, such as the repair half-time and the junction-shift schedule, by applying the IR model to CSI. In addition to a fixed-dose CSI plan, we calculate the NTCP dependence on both repair half-time and dose under the assumption of simple dose distributions.
NTCP with IR model
The cell survival fraction S is given by the linear–quadratic (LQ) model as
where α and β are the LQ model parameters and d k is the dose of fraction k. Therefore, the Poisson LQ NTCP model based on the relative seriality model  and the critical volume model  is expressed as follows:
where M denotes the total number of voxels, d k,i is the dose of fraction k to voxel i, and n is the total number of fractions. N 0 denotes the initial number of cells, v i is the relative volume of voxel i, and s represents the seriality.
To take radiation-induced cell damage into account, we use the IR model , which introduces the concept of “dose equivalent of IR.” The idea is that after a dose of size d, the injury induced by some fraction θ of the dose is still unrepaired by the time the next dose is given. The fraction θ is assumed to decay exponentially with time according to
where Δt is the inter-fractional interval and μ is the repair constant. Using the repair half-time T1/2 as the time needed for half the damaged calls to be repaired, μ is expressed as
By incorporating the cell repair effect on successive doses , the survival fraction derived from the LQ model is given by
where θ q = exp (−μ∆t q ), and Δt q is the time between fractions q and q + 1. Thereafter, the NTCP including cell repair is expressed as
NTCP parameters for the spinal cord were reported as D 50 = 68.6 Gy, γ = 1.90, α/β = 3.00, and a seriality of 4.00 , which were based on radiation tolerance data for clinical myelitis necrosis . However, this NTCP model should have repair term T1/2 potentially because this was derived from clinical fractionated irradiation. On the other study, T1/2 was reported to be 5 h with 95% CI of 0.6 – 9.3 h, which is obtained by irradiation to rat spinal cord with variable intervals . In this study, therefore, we examined the dependence of NTCP on T1/2 varying from 0 h to 10 h to cover 95% CI on the assumption that NTCP reported in Ref.  was given by potentially including T1/2 = 5 h, as following procedure.
Apparent NTCP including repair NTCP repair (T 1/2 ) was calculated by eq. (6) for T1/2 up to 10 h with 1 h intervals.
True NTCP including repair NTCP repair ’(T 1/2 ) was obtained by shifting NTCP repair (T 1/2 ) to NTCP repair (5 h) be equal to original NTCP (NTCP repair (0 h)) by following equation;
This procedure is also described in Fig. 1.
Beam arrangement in CSI
RayStation (ver. 7.2.5, RaySearch Laboratories, Stockholm, Sweden) was used to create retrospective 3D conformal radiation therapy plans for six adult patients receiving CSI from a Siemens Artiste linear accelerator with a 6-MV photon beam, in which the dose distribution is calculated by the collapsed cone convolution superposition method .
For each patient, we prescribed a total dose of 36 Gy (20 fractions of 1.8 Gy each) to the target with the margin of 5 mm to both the spinal cord and the brain, in which the minimum dose was constrained not to be less than 95% of the prescribed dose. The treatment was assumed to be performed in five fractions every week without weekends, as is usual in clinical practice.
We used two opposing lateral fields to treat the entire brain and two abutted spinal fields to treat the upper and lower spinal fields . The brain-field isocenter was located at the volumetric center of the brain, and we used only a longitudinal shift for the spinal-field isocenters. We used the couch angle to match the brain field with the inferior edges of the two lateral beams, and we rotated the collimator to match the upper spinal-field divergence. For the lower spinal field, we rotated the gantry to match the inferior edge divergence of the upper spinal field with a couch rotation of 90°. To facilitate several types of junction shift, we defined three different junction positions in the case of no gap (i.e., a gap size of zero) (Fig. 2). Junction position B was established to match each field exactly. Junction position A was achieved by closing two multi-leaf collimator (MLC) leaves (1 cm) of each anterior and posterior edge of the upper spinal field and by opening two MLC leaves of each posterior edge of the brain field and the anterior edge of the lower spinal field. Junction position C was achieved by moving the MLC leaves in the opposite directions of those for junction position B.
In addition to the case of no gap, gaps of 5 and − 5 mm (i.e., a 5-mm overlap) were created in the same manner. Although these are worst case scenario of systematic errors, it could arise because of setup errors or patient movement. For example, a mistake in counting a patient’s vertebrae could lead to the treatment area being located erroneously, thus resulting in overlapping junction areas .
Schedules of junction shifts
Using the junction positions defined above, some junction-shift schedules were determined for each gap size assuming five fractions per week (Table 1); this process was based on Ref. . Assuming two junction positions (positions A and B), biweekly, weekly, and daily plans refer to changing the junction every two weeks, every week, and every day, respectively. Assuming three junction positions (positions A, B, and C), a 1/3 treatment plan refers to changing the junction every 1/3 of the treatment (i.e., after 7 and 14 fractions). The simplest schedules among each number of junction positions (i.e., the biweekly plan for two junctions and the 1/3 treatment plan for three junctions) were taken as the reference schedules.
Calculation of NTCP in RayStation
First, to estimate the effect of the value of T1/2, the NTCP was calculated with the IR model via the biological evaluation tool RayBiology in RayStation, in which the junction-shift schedule was fixed to the reference schedules for each number of the junction positions.
Second, the NTCP in the case of changing the junction-shift schedule was examined to investigate whether it was affected. For each junction position, the variation in NTCP when changing from the reference schedules to the other schedules was calculated.
Manual calculation of NTCP
To examine the NTCP dependence on both T1/2 and the dose, we calculated the NTCP manually for the spinal cord by assuming many voxels. We began by providing a uniform dose to the entire spinal cord. Thereafter, we gave a uniform dose to 1% of the volume of the spinal cord and no dose to the remaining 99%. To eliminate variations in dose per fraction under a fixed total dose, the total dose was expressed by the equivalent dose in 2 Gy fractions (EQD2) and was varied from 30 to 120 Gy. The fractionation schedule (20 fractions/4 weeks) was the same as those for the CSI planning. Given that we assumed each voxel to be exposed to a fixed dose, no variation in junction-shift schedule was considered in this calculation.
We evaluated the statistical significance of the proposed changes in NTCP. Given the small sample size (n = 6), we used the Shapiro–Wilk test to validate the NTCP normality. Thereafter, we used a paired Student t-test to investigate the statistical significance with a confidence level of 95%.
T1/2 dependence of NTCP considering IR model
Figure 3 shows the NTCP as a function of T1/2 on the reference junction-shift schedules (i.e., the biweekly plan for two junction positions and the 1/3 treatment plan for three junction positions). The NTCP was largest with a gap size of − 5 mm, followed by 0 and 5 mm, although the difference between the latter two was small. With the same gap size, three junction positions resulted in a lower NTCP than two junction positions, particularly for a gap size of − 5 mm.
For all gap sizes, the NTCP with the IR model varied depending on T1/2. In comparison to the NTCP with T1/2 = 5 h, that with the shorter T1/2 was significantly smaller, and that with the longer T1/2 was significantly larger. In addition, the longer T1/2 gave more prominent increase of NTCP, which reached a maximum with T1/2 = 10 h.
Dependence of NTCP on junction-shift schedule with the IR model
Figure 4 shows the NTCP deviation when the junction-shift schedule was changed from the reference schedules to the other schedules for each number of junctions. For two junction positions, no NTCP reduction was observed upon changing to the weekly plan. However, for a gap size of − 5 mm, a statistically significant deviation of NTCP was seen for T1/2 longer than 6 h upon changing to the daily plan; the NTCP deviation increased negatively as T1/2 decreased and reached a maximum of − 0.07 (T1/2 = 10 h). For three junction positions, the NTCP deviation was statistically significant for T1/2 longer than 7 h upon changing to the daily plan for a gap size of − 5 mm; the NTCP deviation increased negatively as T1/2 decreased and reached a maximum of − 0.03 (T1/2 = 10 h).
Manual calculation of NTCP with varying T1/2 and EQD2
We manually calculated the NTCP curves with the IR model by varying EQD2 which is given to 100% or 1% of the volume of spinal cord (Fig. 5). When EQD2 was given to the entire spinal cord (Fig. 5(a)), an EQD2 of 68.6 Gy gave an NTCP of 50% for T1/2 = 0 (i.e., without the IR model). However, the NTCP curve shifted upward upon incorporating the IR model for T1/2 longer than ~ 3 h. Moreover, a larger T1/2 gave a higher NTCP curve. When EQD2 was given to 1% of the volume of the spinal cord (Fig. 5(b)), the NTCP was not saturated with high EQD2 because of partial exposure and continued to increase with increasing EQD2. Although the NTCP was lower than that when the entire spinal cord was irradiated, it increased as T1/2 increased over ~ 3 h, similar to that when the entire spinal cord was irradiated.
Figure 5 also shows the EQD2 that gave an NTCP of 5% for various T1/2, assuming that EQD2 was given to (c) 100% and (d) 1% of the volume of the spinal cord. In Fig. 5(c), an EQD2 of 52 Gy gives an NTCP of 5% for T1/2 = 0, but this value decreases for T1/2 over 3 h and reaches at least 44 Gy (T1/2 = 10 h). In Fig. 5(d), an EQD2 of 58 Gy gives an NTCP of 5% for T1/2 = 0, but this value decreases for T1/2 over 3 h and reaches at least 49 Gy (T1/2 = 10 h).
In this study, the impact of incorporating the IR model on NTCP was evaluated by retrospectively planning CSI for six patients. It should be noted that the study needs more cases to confirm the results because of a limited number of samples. Table 2 summarizes the results in this study.
For CSI planning, NTCP with T1/2 > 6 h was significantly higher than NTCP with T1/2 = 5 h, and it increased as T1/2 increased, particularly with a gap size of − 5 mm. Therefore, the setup of patients for CSI should be conducted carefully, particularly for patients having variations in the number of vertebrae. The increase of T1/2 with the IR model provides a long time for cells to be repaired from radiation damage, thus resulting in high NTCP. Under existing circumstances wherein wide 95% CI of T1/2 are reported, the largest conceivable value of T1/2 should be applied to evaluate the NTCP on the safest side because it largely depends on the value of T1/2. In the future, we need to search for an exact value of T1/2 not for rats but for humans to confirm the NTCP of cell repair. T1/2 for a carbon-ion beam was reported to be shorter than that for a photon beam . Hence, T1/2 should also be investigated for different types of particle radiation beams.
Additionally, we examined the dependence of the NTCP on the junction-shift schedule (Fig. 4). For gap sizes of 5 and 0 mm, the NTCP did not vary with the junction-shift schedule because the maximum dose outside the junction regions was almost unchanged by junction shifts. However, for a gap size of − 5 mm, the NTCP decreased with changing schedules from reference plans to daily plans; we attribute this result to the hot spot in the junction area. For example, with the biweekly plan, an over-dose continued for two weeks and a subsequent mid-dose continued for a further two weeks on a certain point in the junction area, thus restraining cell repair. In contrast, with the daily plan, the over- and mid-doses were alternated every day on a certain point in the junction area, thus fostering cell repair and resulting in a significantly reduced NTCP. In particular, a larger T1/2 resulted in a larger reduction of NTCP because of greater cell repair. A deviation in NTCP associated with different schedules was seen because of tissue with comparatively long T1/2, similar to the case for the spinal cord.
According to the results of manually calculating the NTCP considering the IR model, as shown in Fig. 5(a)-(b), the increase in NTCP due to the increase of T1/2 became prominent for an EQD2 of 60–70 Gy, which corresponds to a dose with a steep normalized gradient of the dose response curve. By contrast, the NTCP hardly increased as T1/2 increased for an EQD2 of 30–40 Gy as irradiated in CSI planning; this result is consistent with the results of CSI planning demonstrated in this study.
The EQD2 with an NTCP of 5% was also calculated with IR model for T1/2 ranging from 1 h to 10 h (Fig. 5(c)-(d)) which can be considered the tolerance dose for normal tissue. Assuming whole-organ irradiation (the entire spinal cord), the tolerance dose of 52 Gy (T1/2 = 0 h) was decreased by 15%, reaching at least 44 Gy (T1/2 = 10 h). Assuming partial organ irradiation (1% volume), the tolerance dose of 58 Gy (T1/2 = 0 h) was decreased by 15%, reaching at least 49 Gy (T1/2 = 10 h). To think on the worst case supposed in this study, the dose should be under 44 Gy for whole-organ irradiation and under 49 Gy for partial (1%) organ irradiation.
Although the NTCP can estimate complications with normal tissue, the parameters of the LQ model and NTCP used in this study contain uncertainties. For example, the α/β value of 0.87 is reported in several papers [18, 19], which is smaller than the value used in this study. In this case, the NTCP could depend on T1/2 more largely, and increase more notably as T1/2 increased, because the value of β is related to the repairable damage as shown in eq. (6). Therefore, the α/β value has an impact on NTCP. Further studies are needed in order to estimate the dependency on the α/β value. Furthermore, other NTCP models have been proposed, such as the Lyman–Kutcher–Burman model [20, 21], and the absolute value of the NTCP depends on those models. However, the relative effects on the NTCP of incorporating the IR model should have the same tendency irrespective of parametric uncertainty or NTCP model.
Although we used a mono-exponential repair model that assumes a single repair constant (i.e., the repair half-time), other studies have reported a bi-exponential repair model that assumes two repair kinetics with long and short repair components [5, 6, 22, 23]. However, Levin-Plotnik et al. detected no difference between the mono- and bi-exponential repair models because the time to complete repair is affected only by the long repair component after an inter-fractional time of 5 h when the short repair component is saturated . Therefore, our results with the mono-exponential model should not be markedly different from those with the bi-exponential repair model for actually used inter-fractional times, such as 24 h. Even if the bi-exponential model were to fit the experimental data, determining the length of the long repair component is more important than that of the short one because the proportion of damage repaired by the long component is larger than that repaired by the short component, as suggested by Ang et al. .
In contrast to physical dosimetric criteria, biological criteria, such as the NTCP, are advantageous to be able to account for not only the cell repair effect but also for several radiobiological phenomena such as clonogenic cell density, radiosensitivity, and hypoxia . By incorporating these biological phenomena into the figures of merit, biologically based treatment planning could increasingly reproduce treatment outcomes toward biologically accurate radiotherapy.
Incorporating an IR model allows the NTCP to be evaluated including the effect of cell repair and to reflect the effect of scheduling. However, the NTCP is significantly affected by the repair half-time T1/2. In the future, using precise values of T1/2 for humans according to the radiation used, biologically based treatment planning using NTCP values that include the of cell repair effect may be clinically performed.
equivalent dose in 2 Gy fractions
normal tissue complication probability
Warkentin B, Stavrev P, Stavreva N, Field C, Fallone BG. A TCP-NTCP estimation module using DVHs and known radiobiological models and parameter sets. J Appl Clin Med Phys. 2004;5:50–63.
Oliver R. A comparison of the effects of acute and protracted gamma-radiation on the growth of seedlings of Visia Faba. II. Theoretical calculations. Int J Radiat Biol Relat Stud Phys Chem Med. 1964;8:475–88.
Levin-Plotnik D, Hamilton RJ, Niemierko A, Akselrod S. A model for optimizing normal tissue complication probability in the spinal cord using a generalized incomplete repair scheme. Radiat Res. 2001;155:593–602.
Kim JJ, Hao Y, Jang D, Wong CS. Lack of influence of sequence of top-up doses on repair kinetics in rat spinal cord. Radiother Oncol. 1997;43:211–7.
Ang KK, Jiang GL, Guttenberger R, Thames HD, Stephens LC, Smith CD, Feng Y. Impact of spinal cord repair kinetics on the practice of altered fractionation schedules. Radiother Oncol. 1992;25:287–94.
Landuyt W, Fowler J, Ruifrok A, Stüben G, van der Kogel A, van der Schueren E. Kinetics of repair in the spinal cord of the rat. Radiother Oncol. 1997;45:55–62.
Jenkin D. The radiation treatment of medulloblastoma. J Neuro-Oncol. 1996;29:45–54.
Paulino AC, Cha DT, Barker JL Jr, Lo S, Manera RB. Patterns of failure in relation to radiotherapy fields in supratentorial primitive neuroectodermal tumor. Int J Radiat Oncol Biol Phys. 2004;58:1171–6.
Urie M, FitzGerald TJ, Followill D, Laurie F, Marcus R, Michalski J. Current calibration, treatment, and treatment planning techniques among institutions participating in the Children's oncology group. Int J Radiat Oncol Biol Phys. 2003;55:245–60.
Källman P, Agren A, Brahme A. Tumour and normal tissue responses to fractionated non-uniform dose delivery. Int J Radiat Biol. 1992;62:249–62.
Stavreva N, Niemierko A, Stavrev P, Goitein M. Modelling the dose-volume response of the spinal cord, based on the idea of damage to contiguous functional subunits. Int J Radiat Biol. 2001;77:695–702.
Ågren Cronqvist A. Quantification of the response of heterogeneous tumors and organized normal tissues to fractionated radiotherapy [Ph.D. thesis]. Stockholm University; 1995.
Emami B, Lyman J, Brown A, Coia L, Goitein M, Munzenrider JE, Shank B, Solin LJ, Wesson M. Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys. 1991;21:109–22.
Mackie TR, Ahnesjö A, Dickof P, Snider A. Development of a convolution/ superposition method for photon beams. Use of Comp in Rad Ther. 1987:107–10.
South M, Chiu JK, Teh BS, Bloch C, Schroeder TM, Paulino AC. Supine craniospinal irradiation using intrafractional junction shifts and field-in-field dose shaping: early experience at Methodist hospital. Int J Radiat Oncol Biol Phys. 2008;71:477–83.
Ibrahim DA, Myung KS, Skaggs DL. Ten percent of patients with adolescent idiopathic scoliosis have variations in the number of thoracic or lumbar vertebrae. J Bone Joint Surg Am. 2013;95:828–33.
Hirayama R, Furusawa Y, Fukawa T, Ando K. Repair kinetics of DNA-DSB induced by X-rays or carbon ions under oxic and hypoxic conditions. J Radiat Res. 2005;46:325–32.
Schultheiss TE. The radiation dose-response of the human spinal cord. Int J Radiat Oncol Biol Phys. 2008;71:1455–9.
Kirkpatrick JP, van der Kogel AJ, Schultheiss TE. Radiation dose-volume effects in the spinal cord. Int J Radiat Oncol Biol Phys. 2010;76:S42–9.
Niemierko A, Goitein M. Calculation of normal tissue complication probability and dose-volume histogram reduction schemes for tissues with a critical element architecture. Radiother Oncol. 1991;20:166–76.
Semenenko VA, Li XA. Lyman-Kutcher-Burman NTCP model parameters for radiation pneumonitis and xerostomia based on combined analysis of published clinical data. Phys Med Biol. 2008;53:737–55.
Millar WT, Canney PA. Derivation and application of equations describing the effects of fractionated protracted irradiation, based on multiple and incomplete repair processes. Part I. Derivation of equations. Int J Radiat Biol. 1993;64:275–91.
Thames HD. Repair kinetics in tissues: alternative models. Radiother Oncol. 1989;14:321–7.
Levin-Plotnik D, Niemierko A, Akselrod S. Effect of incomplete repair on normal tissue complication probability in the spinal cord. Int J Radiat Oncol Biol Phys. 2000;46:631–8.
Allen Li X, Alber M, Deasy JO, Jackson A, Ken Jee KW, Marks LB, Martel MK, Mayo C, Moiseenko V, Nahum AE, Niemierko A, Semenenko VA, Yorke ED. The use and QA of biologically related models for treatment planning: short report of the TG-166 of the therapy physics committee of the AAPM. Med Phys. 2012;39:1386–409.
The authors would like to thank Hiroki Adachi and Jun Ishishita of Hitachi Ltd., Tokyo, Japan, for their help with the RayStation treatment planning system.
This work was supported by the Japan Society for the Promotion of Science KAKENHI Grant.
[Grant-in-Aid for Scientific Research (B) JP17H04119].
Availability of data and materials
Please contact author for data requests.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wakisaka, Y., Yagi, M., Sumida, I. et al. Impact of time-related factors on biologically accurate radiotherapy treatment planning. Radiat Oncol 13, 30 (2018). https://doi.org/10.1186/s13014-018-0973-6