 Research
 Open Access
 Published:
Radiotoxicity in robotic radiosurgery: proposing a new quality index for optimizing the treatment planning of brain metastases
Radiation Oncology volume 12, Article number: 136 (2017)
Abstract
Background
As irradiated brain volume at 12 Gy (V12) is a predictor for radionecrosis, the purpose of the study was to develop a model for Cyberknife (CK) plans that is able to predict the lowest achievable V12 at a given tumor size and prescription dose (PD), and to suggest a new quality index regarding V12 for optimizing the treatment planning of brain metastases.
Method
In our model V12 was approximated as a spherical shell around the tumor volume. The radial distance between tumor surface and the 12 Gy isodose line was calculated using an approximation of the mean dose gradient in that area. Assuming a radially symmetrical irradiation from the upper half space, the dose distribution is given by the superposition of single fields. The dose profiles of a single field were derived by the measured offcenter ratios (OCR) of the CK system. Using the calculated gradients of the sum dose profiles, minimalV12 was estimated for different tumor sizes. The model calculation was tested using a phantom dataset and retrospectively applied on clinical cases.
Results
Our model allows the prediction of a bestcase scenario for V12 at a given tumor size and PD which was confirmed by the results of the isocentric phantom plans. The results of the nonisocentric phantom plans showed that an optimization of coverage caused an increase in V12. This was in accordance with the results of the retrospective analysis. V12 s of the clinical cases were on average twice that of the predicted model calculation. A good agreement was achieved for plans with an optimal conformity index (nCI). Replanning of cases with high V12 showed that lower values could be reached by selecting smaller collimators and by allowing a larger number of total MU and more MU per beam.
Conclusions
V12 is a main parameter for assessing plan quality in terms of radiotoxicity. The index f12 defined as the ratio of V12 from the actual plan with the evaluated V12 from our model describes the conformity of an optimally possible V12 and thus can be used as a new quality index for optimizing treatment plans.
Background
Stereotactic radiosurgery (SRS) is an established and widely used technique for the treatment of single and multiple brain metastases [1,2,3]. For the treatment of brain metastases, different radiosurgical modalities such as gamma knife (GK), LINACSRS or the robotic Cyberknife system (CK) exist, which all show dosimetrically comparable plan quality [4,5,6,7]. CK (Accuray, Sunnyvale, CA, USA) has the advantage of a frameless and imageguided targeting technology, which not only enables wholebody radiosurgery, but also hypofractionated treatment of large tumors (>3 cm) [8, 9]. The concept of hypofractionated treatment and the possibility to treat patients with up to ten brain metastases without wholebrain radiotherapy [10] underlines the attractiveness of SRS as a minimalinvasive, effective and safe therapy option for brain metastases. Nevertheless, side effects may arise whereby the development of brain radionecrosis represents the most common complication after SRS [11,12,13,14,15,16].
The risk of radionecrosis increases with larger doses and sizes of the target volume. Besides risk factors such as tumor biology [17], the volume of healthy brain tissue that is irradiated with a single dose of more than 10 Gy (V10) or 12 Gy (V12) is a predictor of brain radionecrosis [11, 12, 18,19,20,21]. Blonigen et al. and Minniti et al. have shown that the risk of radionecrosis increases gradually with V10 and V12 [20, 21]. Bohoudi et al. have presented a formula for estimating the optimal individual radiosurgery dose for isotoxic treatment planning of brain metastases [22]. In our study we did not aim to optimize the prescription dose in order not to exceed a maximum accepted V12. The objective of our study was to develop a model that allows the prediction of a minimal achievable V12 for the treatment of brain metastases with the CK system. Based on the presented model we suggest a new quality parameter which can be used for assessing V12 in the context of other plan quality criterions.
Methods
Derivation of a model for estimating a minimalV12
Physical and geometrical assumptions have been made to develop a predictive model for estimating a minimalV12 at a given tumor volume and prescription dose (PD). One assumption is that all metastases are spherical and that the prescription isodoseline can be fitted perfectly to the target surface. Thereby V12 is given as spherical shell around the planning tumor volume (PTV). For singular metastases V12 can be evaluated as follows:
where Δr_{12} is the radial distance between the tumor surface and the 12 Gy isodoseline. It depends on the prescribed dose and steepness of the dose gradient \( \overrightarrow{\nabla}D \) and can be derived using an approximation of the mean dose gradient \( \left\langle \left\overrightarrow{\nabla}D\right\right\rangle \) between the tumor surface and the 12 Gy isodoseline:
ΔD describes the dose difference between PD and 12 Gy. For the approximation of the mean dose gradient the measured offcenter ratios (OCR) of the CK at a depth of 10 cm with an 80 cm sourcedetector distance have been used (Fig. 1a). Assuming an isocentric irradiation technique from the upper half space of a sphere that represents the possible beam setup for the CK system, the dose distribution of a spherical metastasis is given by the superposition of single fields, which summed up to a cumulative dose profile [23]:
The radial distance R is given by r/sinθ where r is the radial distance of the central axis. All single fields were equally weighted, the depthdose has been neglected and an optimal coverage of the PTV has been expected. A program for dose calculation was developed and implemented using IDL (Exelis Visual Information Solutions, Boulder, CO, USA). The sum dose profiles were calculated for different collimator sizes using Eq. 3 (Fig. 1b) and the corresponding dose gradients were determined for all dose profiles by numerical differentiation (Fig. 1c). Assuming the availability of a continuously adjustable collimator that would be able to precisely cover any tumor size with optimum conformity, it is possible to find a fit formula which describes the dose falloff of any desired tumor size (Fig. 1d):
Using Eqs. 4 and 2 minimalV12 can be calculated for any desired tumor volume at a given PD and can be compared to that of the CK plan.
For multiple metastases the total minimalV12 can be derived as the sum of the calculated V12 s (Eq. 1) from each individual metastasis:
Another possibility would be only considering the total target volume so that the total minimalV12 is given by:
Model validation and application
A set of virtual CT phantoms, created with MATLAB (MathWorks, Version 9.0) were used for model validation. The phantoms comprise a spherical volume filled with water (HU = 0) and embedded in a bone shell (HU = 1000). The Hounsfield values outside the sphere were set to −1000 (air). Spherical tumors of different size were set in the centre of each phantom. For ten CK collimators (diameter 5 up to 40) an optimal target radius using the radial distance R of the sum dose profiles at 65% (see Fig. 1b) was used for generating the spherical tumors. To simulate treatment cases as close as possible to our model assumptions, isocentric CK plans of the phantom dataset were performed using MultiPlan (Accuray, Version 4.5) (Fig. 2a). A PD of 65% at 20 Gy, commonly used in our department for the treatment of brain metastases, was applied. V12 s of CK plans were compared with the calculated minimalV12 of our model. Additionally we replanned all phantom cases using the sequential planning option of MultiPlan. We chose the same optimization planning strategy as for our clinical cases to analyze the effect of the nonisocentric irradiation technique on V12 compared to our model. To assess the application of our model on clinical cases, we analyzed retrospectively a patient collective comprising the treatment of 176 metastases and 67 patients. The patient collective was chosen randomly (in alphabetical order) from our data pool. The planning strategy of the sequential optimization included coverage (Cov = PTV covered by the prescription isodose volume (PIV)/PTV), three shells for optimizing conformity (CI = Cov * (PTV covered by the PIV)/PIV) [24, 25] and minimal monitor units (MU) as inverse planning objectives. The collimator selection was made manually for each case. Several collimator settings were tested to find the best solution in terms of dosage conformity, minimal dose in target volume, coverage and V10 and V12, respectively. Beam reduction was performed to allow a minimum MU of 10 per beam. MU limits for maximum MU per beam varied between 150 and 200 depending on PTV size and number of metastases. All treatments were performed in one fraction. The deviation of the predicted V12 s from the resulting V12 s of the CK plans was described by index f12, which is given by the following ratio:
The influence of two calculation methods (Eq. 5 versus Eq. 6) on V12 estimation for multiple metastases was compared and the dependence of index f12 on different factors was analyzed in terms of target volume, tumor location, number of metastases, new conformity index (nCI = 1/CI) [26], treatment time and total number of MUs.
Results
Model derivation and validation
The calculated sum dose profiles for the 12 existing collimator sizes showed that the slope of the dose profiles varies depending on the field size (Fig. 1b). Small field sizes showed sharp dose gradients whereas larger field size gradients were flatter. The steepest dose gradient was in the range of 55%  65% for collimator diameters between 5 and 15 mm and drifted to higher dose levels (65%  70%) for collimator diameters between 20 and 60 mm (Fig. 1c). The mean dose gradient decreased exponentially with increasing PTV radius (Fig. 1d). Results of the exponential fits (cf. equation 4) are shown in Table 1.
Resulting V12 s of the isocentric plans in the phantom dataset showed a good agreement with the predicted V12 s of our model (Fig. 2b). Index f12 was independent on tumor size for both irradiation techniques. For isocentric plans, f12 ranged from 1.06 to 1.26 (mean = 1.18 ± 0.06) and showed no dependence on coverage (mean = 95.6% ± 2.0%, range: 91.9%  97.8%) and nCI (mean = 1.13 ± 0.03, range: 1.10–1.20) (Fig. 2c, d). The sequential optimization yielded results with better coverage (mean = 99.4% ± 0.9%, range: 97.3%  100%) and lower nCI (mean = 1.07 ± 0.02, range: 1.04–1.11) compared to isocentric plans, but at the cost of larger f12 values (mean = 1.41 ± 0.17, range: 1.06–1.62). Index f12 and coverage showed an opposite trend while using the sequential optimization technique, whereas in terms of nCI no trend was observable (Fig. 2c, d). A better coverage is per definition associated with a lower nCI; this effect was dominant for isocentric plans, whereas due to the sequential optimization technique both coverage and conformity were on a high level for nonisocentric plans (Fig. 2e).
Model application on patient data
The retrospective analysis of the patient collective included 40 cases of singular and 40 cases of multiple metastases (176 metastases in total). The number of metastases per treatment ranged from 1 to 7 and was on average 2.2. In 72 cases PD was 65% at 20 Gy, in 4 cases 65% at 18 Gy, in 3 cases 65% at 22 Gy and in one case 65% at 16 Gy. The total PTV per treatment ranged from 0.09 ml to 46.07 ml and was on average 4.47 ml. The PTV per metastasis ranged from 0.04 ml to 39.88 ml and was on average 1.78 ml. The mean coverage of the CK plans of the patient data was 99.4% ± 1.2% (range: 92.1% to 100%) and the treatment time per metastasis was less than 60 min.
Our data showed that the resulting V12 s of the CK plans are influenced by the number of metastases. Cases that had a similar total target volume but different number of metastases showed V12 increasing with the number of metastases (Fig. 3a). For V12 determination of multiple metastases, two calculation options (Eq. 5 versus Eq. 6) are possible which causes different results (Fig. 3b). The deviation between the two calculation methods was on average 19.5% ± 13.3% (range: 3.4% to 57.6%). There was a weak tendency towards greater deviations with the increasing number of metastases (Fig. 3b). The deviation between both calculation methods increased when treating small and similarly sized tumor volumes. For example, one case with 3 metastases treated simultaneously (similar tumor volumes: V_{1} = 0.18 ml, V_{2} = 0.14 ml, V_{3} = 1.11 ml) resulted in a deviation of 38%, whereas another case with 3 metastases (one large tumor volume V_{1} = 4.41 ml, two smaller lesions: V_{2} = 0.21 ml, V_{3} = 0.12 ml) resulted in a deviation of only 8%. Index f12 showed no dependence on the number of metastases while using Eq. 5 for the determination of V12 (Fig. 3c) whereas a trend to larger f12 values with the increasing number of metastases was observable when using Eq. 6 (Fig. 3d).
In the following analyses, we used Eq. 5 for V12 estimation of multiple metastases.
For all treatments (40 cases of singular and 40 cases of multiple metastases), f12 ranged from 1.12 (good agreement) up to 4.22 (poor agreement) and was on average 1.91 ± 0.47. Results of index f12 (separated by tumor location and number of metastases) are shown in Table 2.
Tumor location showed little influence on the resulting V12 s of the CK plans (Fig. 4). Peripheral tumors (concerns 11 cases) showed less V12 and lower f12 values compared to tumors located in the middle of the parenchyma (Fig. 4c and d). To avoid a bias when evaluating the dependence of index f12 in terms of target volume, nCI, treatment time and total number of MUs, the 11 cases of peripheral tumors were subsequently excluded.
Analyzing the dependence of index f12 on the factors mentioned above revealed that f12 is independent of PTV size (Fig. 5a), whereas V12 increases with growing tumor volume (Fig. 3a and Fig. 4a). Better plan conformity (nCI close to one) was associated with smaller f12 values (slope of the linear fit = 3, R^{2} = 0.52, Fig. 5b). For the used patient collective, nCI ranged from 1.05 up to 1.78 (mean = 1.21 ± 0.11). Index f12 was independent of treatment time (Fig. 5c) and of the total number of MUs (Fig. 5d).
Discussion
Risk analyses of stereotactic radiosurgery have shown that V10 and V12 act as a predictor for brain radionecrosis and thus representing an important parameter for radiotoxicity [11, 12, 18,19,20,21]. In our study we present a model that allows the prediction of the lowest achievable V12 at a given tumor volume and prescription dose. The model follows optimum conditions as it assumes a radially symmetrical distribution of the beams from the upper half space, optimal target coverage and plan conformity, as well as spherical target volumes. Considering the underlying assumptions, our model provides a bestcase scenario for estimating a minimalV12. The validity of the model was tested by isocentric plans of a phantom dataset that led to a mean f12 value of 1.2. Deviation from unity might be explained by using voxel data for phantom generation and due to deviations from radially symmetrical dose distribution. Nonisocentric irradiation technique includes the solution space of the isocentric beam setup, but also offers the advantages of sequential optimization which can lead to plan improvements concerning coverage and dosage conformity [6].
Plan quality is characterized by different metrics such as the volumetrically defined coverage, the minimal dose in target volume, dosage conformity and the proportion of the healthy brain volume that is irradiated by more than 10 Gy and 12 Gy, respectively. The planning strategy used in our department aimed at achieving a high coverage (Cov > 98%), an optimal dosage conformity (nCI < 1.3) and a moderate treatment time (30 min1 h per target). The plan quality was evaluated individually for each case by a decision board (2 physicians, 1 physicist) having regard to patients history, patients age, tumor localization and the number of metastases. In general a high coverage is required to ensure tumor control and to prevent recurrences [27]. Based on our retrospective analyses there were individual cases (tumors located in eloquent brain regions, for instance brainstem metastases or located close to critical structures) in which the physicians accepted a lower coverage or reduced the PD, since a compromise in coverage or dose reduction results in a lower exposure of the normal brain tissue. A reasonable balance between sufficient PD (coverage) to ensure local control and V12 to reduce the risk of radionecrosis is important. The QUANTEC review of Marks et al. recommends keeping V12 lower than 5–10 ml [27]. Our model allows predicting at what tumor size a threshold value of V12 is reached. For example, the group of Minniti et al. [21] recommends considering lesions with V12 > 8.5 ml for hypofractionated SRS to reduce the risk of radionecrosis. Assuming the bestcase scenario proposed by our model, V12 of 8.5 ml is achieved at a tumor size of 10.9 ml (tumor diameter of 2.75 cm) when using 65% prescription dose at 20 Gy.
Results of the nonisocentric phantom plans suggested that an expected f12 value of 1.4 on average is possible on the basis of the described sequential optimization strategy. In regard to provide a plan with sufficiently high coverage and since the difference between V12 s of the nonisocentric plans and V12 s of the isocentric plans was low (<1 ml) for small tumors (PTV <10 ml, cf. Fig. 2a), this seems to be an acceptable value of f12.
The application of the model to the clinical cases led to an f12 value of on average 1.9. Index f12 improved with better conformity index nCI and showed no dependence on PTV size or on other plan parameters such as treatment time, total number of MUs, and number of metastases. When keeping the coverage constant at a high level, nCI will be improved if the ratio of PIV to the proportion of PTV covered by the prescribed dose is decreased. Since f12 is defined by the ratio of V12 to the optimum V12 calculated by our model, f12 acts like a conformity index for the 12 Gy isodose level and improves with lower nCI values.
Results of the retrospective analysis showed that f12 values of singular (f12_{singular} = 1.80 ± 0.39) and multiple metastases (f12_{multiple} = 2.02 ± 0.52) were comparable within the standard deviation. Depending on the distance between multiple metastases, dose bridges may occur, which can result in higher V12 s compared to singular metastases. Our model takes the number of metastases into account (Eq. 5), but not the effect of dose bridges. Narayanasamy et al. showed that the total treatment volume is a better predictor of whole brain dose from gamma knife based SRS than the number, shape, or location of the lesions [28] while Sahgal et al. found a significant increase in V12 with increasing number of targets for multitarget SRS [29]. In our study, we found that tumor location shows little influence on the resulting V12 s of the CK plans, while V12 is dependent on the number of metastases and a more precise determination of V12 for multiple metastases is achieved when applying the model to each individual metastasis (cf. Fig. 3).
In order to check whether the proposed f12 index can be used as new planning objective for optimizing the treatment planning of brain metastases, we replanned some cases of the retrospective analysis (Additional file 1: Table S1). All plans aimed at achieving a high coverage and dosage conformity. The plan comparison showed that for singular metastases lower f12 (V12) can be reached when using a various number of collimators (combination of small and larger ones) and by allowing a longer treatment time which is connected with a larger number of total MUs. Performing additional time reduction optimization steps resulted in increased f12 (V12) values. For multiple metastases a lower f12 (V12) was reached by changing the preconditions of MUlimits (allowing more MU per beam, no conditions for total MU or max. MU per node). A disadvantage of such a strategy is that it can lead to dose fingers in the lowdose range. The CK offers only two collimation systems, a limited number of circular collimators or a variable circular Iris collimator, both providing the same selection of circular fields. The application of a micromultileafcollimator (MLC) (for instance the new “InCise MLC”, available for CK M6 system), which enabled the production of irregularshaped fields, might be a good alternative [30, 31]. Previous planning studies showed a notable reduction in treatment time and in total MU [32, 33].
For assessing plan quality different decision objectives have to be taken into account.
Moreover, the plan quality depends on planning time (testing different beam configurations and collimator settings) and the experience of the treatment planner, and varies between different sites [34]. The proposed f12 index is an instrument which allows an estimation of how far the actual V12 is from the bestcase scenario at a given tumor size (independently for small or medium/large tumors) and at a given PD. This information may help to save planning time and to reduce the possible risk of radionecrosis.
The model prediction of V12 and the clinical application of the model have some limitations. The model does not consider the location of the metastases. Theoretically an f12 index less than one is possible for peripheral tumors at the margin to bone or for metastases located close to other structures outside the brain contour because this would lead to an overestimation of the V12model. Another issue concerns the shape of the tumors; our model assumes spherical metastases that do not reflect the reality. Any deviation from the radially symmetrical tumor shape could influence the dose gradient and lead to larger V12 s. A possible solution to take this aspect into account could be an estimation of the sphericity of the tumors and to implement it in the model calculation.
Our study is focused on the treatment planning of brain metastases with CK, but it is conceivable to transfer the described methods and model assumptions to other systems as well.
Conclusion
Achieving a minimal V12 for the treatment plans of brain metastases is desirable to lessen the risk of radionecrosis and to spare the healthy surrounding brain. In our study we suggest a new quality parameter (index f12), which is a measure of how far the actual V12 of the plan is from the predicted minimalV12. Thus, f12 describes the conformity of V12 with the minimal possible V12. It is volume independent and can be used for optimizing treatment plans of singular and multiple brain metastases equally.
Abbreviations
 CK:

Cyberknife
 GK:

Gamma knife
 MU:

Monitor unit
 NCI:

new conformity index
 OCR:

Offcenter ratios
 PD:

Prescription dose
 PIV:

Prescription isodose volume
 PTV:

Planning tumor volume
 SRS:

Stereotactic radiosurgery
 V10:

Irradiated brain volume (minus tumor volume) at 10 Gy
 V12:

Irradiated brain volume (minus tumor volume) at 12 Gy
References
 1.
Trifiletti DM, Lee CC, Winardi W, et al. Brainstem metastases treated with stereotactic radiosurgery: safety, efficacy, and dose response. J NeuroOncol. 2015;125:385–92.
 2.
Moraes FY, Taunk NK, Marta GN, Suh JH, Yamada Y. The rationale for targeted therapies and stereotactic Radiosurgery in the treatment of brain metastases. Oncologist. 2016;21:244–51.
 3.
Knoll MA, Oermann EK, Yang AI, et al. Survival of patients with multiple intracranial metastases treated with stereotactic Radiosurgery does the number of tumors matter? Am J Clin Oncol. 2016;00:000.
 4.
Yu C, Jozsef G, Apuzzo ML, Petrovich Z. Dosimetric comparison of CyberKnife with other radiosurgical modalities for an ellipsoidal target. Neurosurgery. 2003;53:1155–62.
 5.
Sio TT, Jang S, Lee SW, Curran B, Pyakuryal AP, Sternick ES. Comparing gamma knife and cyberknife in patients with brain metastases. J Appl Clin Med Phys. 2014;15:4095.
 6.
Treuer H, Hoevels M, Luyken K, et al. Intracranial stereotactic radiosurgery with an adapted linear accelerator vs. robotic radiosurgery: comparison of dosimetric treatment plan quality. Strahlenther Onkol. 2015;191:470–6.
 7.
Cho YH, Lee JM, Lee D, et al. Experiences on two different stereotactic radiosurgery modalities of gamma knife and Cyberknife in treating brain metastases. Acta Neurochir. 2015;157:2003–9.
 8.
Adler JR Jr, Chang SD, Murphy MJ, Doty J, Geis P, Hancock SL. The CyberKnife: a frameless robotic system for radiosurgery. Stereotact Funct Neurosurg. 1997;69:124–8.
 9.
Wowra B, Muacevic A, Tonn JC. CyberKnife Radiosurgery for brain metastases. Prog Neurol Surg Basel Karger. 2012;25:201–9.
 10.
Yamamoto M, Serizawa T, Shuto T, et al. Stereotactic radiosurgery for patients with multiple brain metastases (JLGK0901): a multiinstitutional prospective observational study. Lancet Oncol. 2014;15:387–95.
 11.
Voges J, Treuer H, Sturm V, et al. Risk analysis of linear accelerator radiosurgery. Int J Radiat Oncol Biol Phys. 1996;36:1055–63.
 12.
Flickinger JC, Kondziolka D, Pollock BE, Maitz AH, Lunsford LD. Complications from arteriovenous malformation radiosurgery: multivariate analysis and risk modeling. Int J Radiat Oncol Biol Phys. 1997;38:485–90.
 13.
Chin LS, Ma L, DiBiase S. Radiation necrosis following gamma knife surgery: a casecontrolled comparison of treatment parameters and longterm clinical follow up. J Neurosurg. 2001;94:899–904.
 14.
Williams BJ, Suki D, Fox BD, et al. Stereotactic radiosurgery for metastatic brain tumors: a comprehensive review of complications. J Neurosurg. 2009;111:439–48.
 15.
Kohutek ZA, Yamada Y, Chan TA, et al. Longterm risk of radionecrosis and imaging changes after stereotactic radiosurgery for brain metastases. J NeuroOncol. 2015;125:149–56.
 16.
Le Rhun E, Dhermain F, Vogin G, Reyns N, Metellus P. Radionecrosis after stereotactic radiotherapy for brain metastases. Expert Rev Neurother. 2016;16:903–14.
 17.
Miller JA, Bennett EE, Xiao R, et al. Association between radiation necrosis and tumor biology following stereotactic radiosurgery for brain metastasis. Neurosurgery. 2016;63(Suppl 1):196–7.
 18.
Korytko T, Radivoyevitch T, Colussi V, et al. 12 Gy gamma knife radiosurgical volume is a predictor for radiation necrosis in nonAVM intracranial tumors. Int J Radiat Oncol Biol Phys. 2006;64:419–24.
 19.
Lawrence YR, Li XA, el Naqa I, et al. Radiation dosevolume effects in the brain. Int J Radiat Oncol Biol Phys. 2010;76:20–7.
 20.
Blonigen BJ, Steinmetz RD, Levin L, Lamba MA, Warnick RE, Breneman JC. Irradiated volume as a predictor of brain radionecrosis after linear accelerator stereotactic radiosurgery. Int J Radiat Oncol Biol Phys. 2010;77:996–1001.
 21.
Minniti G, Clarke E, Lanzetta G, et al. Stereotactic radiosurgery for brain metastases: analysis of outcome and risk of brain radionecrosis. Radiat Oncol. 2011;6:48.
 22.
Bohoudi O, Bruynzeel AM, Lagerwaard FJ, Cuijpers JP, Slotman PJ, Palacios MA. Isotoxic radiosurgery planning for brain metastases. Radiother Oncol. 2016;120(2):253–7.
 23.
Treuer U, Treuer H, Hoevels M, Müller RP, Sturm V. Computerized optimization of multiple isocentres in stereotactic convergent beam irradiation. Phys Med Biol. 1998;43:49–64.
 24.
Paddick I. A simple scoring ratio to index the conformity of radiosurgical treatment plans. J Neurosurg. 2000;93(Suppl 3):219–22.
 25.
van’t Riet A, Mak AC, Moerland MA, Elders LH, van der Zee W. A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. Int J Radiat Oncol Biol Phys. 1997;37:731–6.
 26.
Nakamura JL, Verhey LJ, Smith V, et al. Doseconformity of gamma knife radiosurgery and risk factors for complications. Int J Radiat Oncol Biol Phys. 2001;51:1313–9.
 27.
Marks LB, Yorke ED, Jackson A, et al. The use of normal tissue complication probability (NTCP) models in the clinic. Int J Radiat Oncol Biol Phys. 2010;76:10–9.
 28.
Narayanasamy G, Smith A, Van Meter E, Mc Garry R, Molloy JA. Total target volume is a better predictor of whole brain dose from gamma stereotactic radiosurgery than the number, shape, or location of the lesions. Med Phys. 2013;40(9):091714.
 29.
Sahgal A, Barani IJ, Novotny J, et al. Prescription dose guideline based on physical criterion for multiple metastatic brain tumors treated with stereotactic Radiosurgery. Int J Radiat Oncol Biol Phys. 2010;78:605–8.
 30.
Asmerom G, Bourne D, Chappelow J, et al. The design and physical characterization of a multileaf collimator for robotic radiosurgery. Biomed Phys Eng Express. 2016;2:017003.
 31.
Fürweger C, Prins P, Coskan H, Heijmen BJ. Characteristics and performance of the first commercial multileaf collimator for a robotic radiosurgery system. Med Phys. 2016;43(5):2063–71.
 32.
McGuinness CM, Gottschalk AR, Lessard E, et al. Investigating the clinical advantages of a robotic linac equipped with a multileaf collimator in the treatment of brain and prostate cancer patients. J Appl Clin Med Phys. 2015;16(5):284–95.
 33.
Goggin LM, Descovich M, McGuinness C, Shiao S, Pouliot J, Park C. Doisimetric comparison between 3dimensional conformal and robotic SBRT treatment plans for accelerated partial breast radiotherapy. Technol Cancer Res Treat. 2016;15(3):437–45.
 34.
Blanck O, Wang L, Baus W, et al. Inverse treatment planning for spinal robotic radiosurgery: an international multiinstitutional benchmark trial. J Appl Clin Med Phys. 2016;17(3):6151.
Acknowledgements
Not applicable.
Funding
No outside funding was obtained for this study.
Availability of data and materials
The phantom dataset can be requested from the corresponding author. The clinical raw data is provided in the Departement of Stereotaxy and Functional Neurosurgery, Centre of Neurosurgery, University Hospital of Cologne.
Author information
Affiliations
Contributions
HT and AH conceived the study and developed the model. AH created the phantom dataset, carried out the analysis and interpretation of data and wrote the manuscript. KL, MH and AG were involved in data acquisition and analysis. DR and MIR participated in the clinical, radiation and image data collection. WWB, MK, and HT helped in interpreting the data and critically reviewed the manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Alexandra Hellerbach.
Ethics declarations
Ethics approval and consent to participate
In this retrospective study, ethical approval was waived by our institutional review board. Due to the professional code of conduct of the Medical Association’s of Nordrhein from 15th of November 2015 in there is no need for ethics approval in retrospective researches (§15, article 1).
Consent for publication
Not applicable.
Competing interests
MIR reports financial project related support from Accuray, outside the submitted work. All other authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1:
Plan comparison and indicators for low or high f12 indices. (DOCX 29 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Hellerbach, A., Luyken, K., Hoevels, M. et al. Radiotoxicity in robotic radiosurgery: proposing a new quality index for optimizing the treatment planning of brain metastases. Radiat Oncol 12, 136 (2017) doi:10.1186/s130140170867z
Received
Accepted
Published
DOI
Keywords
 Stereotactic radiosurgery
 Robotic radiosurgery
 Brain metastases
 Radiotoxicity
 Dosimetry
 Dose falloff prediction
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.