

ORIGINAL ARTICLE 



Year : 2018  Volume
: 43
 Issue : 1  Page : 1622 

Comparison of volumetric modulated arc therapy and intensity modulated radiation therapy for whole brain hippocampal sparing treatment plans based on radiobiological modeling
Ethan Kendall, Ozer Algan, Salahuddin Ahmad
Department of Radiation Oncology, Peggy and Charles Stephenson Cancer Center, University of Oklahoma HSC, Oklahoma City, OK, USA
Date of Submission  20Jul2017 
Date of Decision  29Nov2017 
Date of Acceptance  29Nov2017 
Date of Web Publication  12Mar2018 
Correspondence Address: Dr. Salahuddin Ahmad Department of Radiation Oncology, Peggy and Charles Stephenson Cancer Center, University of Oklahoma HSC, 800 NE. 10^{th} Street, OKCC L100, Oklahoma City, OK 73104 USA
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/jmp.JMP_85_17
Abstract   
Introduction: In this article, we report the results of our investigation on comparison of radiobiological aspects of treatment plans with linear acceleratorbased intensitymodulated radiation therapy and volumetricmodulated arc therapy for patients having hippocampal avoidance wholebrain radiation therapy. Materials and Methods: In this retrospective study using the dosevolume histogram, we calculated and compared biophysical indices of equivalent uniform dose, tumor control probability, and normal tissue complication probability (NTCP) for 15 wholebrain radiotherapy patients. Results and Discussions: Doseresponse models for tumors and critical structures were separated into two groups: mechanistic and empirical. Mechanistic models formulate mathematically with describable relationships while empirical models fit data through empirical observations to appropriately determine parameters giving results agreeable to those given by mechanistic models. Conclusions: Techniques applied in this manuscript could be applied to any other organs or types of cancer to evaluate treatment plans based on radiobiological modeling.
Keywords: Hippocampal sparing treatment, radiobiological modeling, VMAT versus IMRT
How to cite this article: Kendall E, Algan O, Ahmad S. Comparison of volumetric modulated arc therapy and intensity modulated radiation therapy for whole brain hippocampal sparing treatment plans based on radiobiological modeling. J Med Phys 2018;43:1622 
How to cite this URL: Kendall E, Algan O, Ahmad S. Comparison of volumetric modulated arc therapy and intensity modulated radiation therapy for whole brain hippocampal sparing treatment plans based on radiobiological modeling. J Med Phys [serial online] 2018 [cited 2019 May 22];43:1622. Available from: http://www.jmp.org.in/text.asp?2018/43/1/16/227073 
Introduction   
Brain metastases are the most common intracranial tumors in adults and occur in 10%–30% of cancer patients.^{[1]} One of the treatment options for patients with multiple brain metastases is wholebrain radiation therapy (WBRT) which can be combined with radiosensitizers such as motexafin gadolinium, efaproxiral, and temozolomide.^{[1]} Neurological decline may occur in a patient with WBRT treatment causing cerebral edema in the patient within the first few weeks of treatment, symptoms of which include headaches, somnolence, vomiting, nausea, and inability to focus.^{[2]} One to six months after WBRT, subacute encephalopathy may occur, which increases patient headaches, somnolence, and fatigue. Six months after WBRT, patients may experience severe dementia and memory loss. Irradiating hippocampi is thought to be a factor in WBRTassociated cognitive decline because the hippocampi are essential in learning, retrieval of information, and the formation of new memories. One of the ways to reduce cognitive decline is to avoid irradiating the hippocampus when performing WBRT in a process called hippocampalavoidance WBRT (HAWBRT).^{[2]} It is possible to create WBRT treatment plans that spare hippocampi with equivalent dose distributions to conventional WBRT treatment plans. The most common delivery techniques for HAWBRT are intensitymodulated radiation therapy (IMRT) and volumetricmodulated arc therapy (VMAT). VMAT is a subset of IMRT where radiation is delivered in an arc instead of at static angles.
The most common biological metrics used to guide clinical predictions are tumor control probability (TCP) and normal tissue complication probability (NTCP). TCP is the probability of localized tumor control for a given dose distribution, while NTCP is the probability of an undesirable effect from an organ due to radiation. A challenge of radiation therapy is to find a balance of these values to give each patient the optimal treatment by raising TCP as high as possible, while lowering NTCPs as low as possible. Doseresponse models for tumors and critical structures are separated into two groups: mechanistic and empirical. Mechanistic models formulate mathematically with describable relationships while empirical models fit data through empirical observations. We may not be able to accurately formulate some biological processes because of their complexity to make mechanistic models, but this is not a challenge for empirical models which simply need to fit data.^{[3]} Dosevolume histogram (DVH) shows how much dose is being received by a structure. It can be used to view physical parameters such as minimum, mean, and maximum dose to calculate TCP and NTCP. Because of limited information in the setting of HAWBRT, the purpose of our study was to calculate and compare target doses, normal tissue doses, TCP, and NTCP values for treatment plans utilizing IMRT or VMATbased on HAWBRT.
Materials and Methods   
Treatment planning for 15 patients was carried out with Varian Eclipse treatment planning system using 3arc VMAT (Rapid Arc) and 9field stepandshoot (S and S) IMRT, both calculated with Anisotropic Analytical Algorithm. From this point in our manuscript, whenever we refer to IMRT, we are referring to S and S IMRT. The planning target volume (PTV) was defined as the wholebrain excluding a uniform threedimensional 5mm expansion around the hippocampus. Prescribed doses in all plans were 30 Gy delivered over 10 fractions normalized to a minimum of 95% of the target volume receiving 100% of the prescribed dose. Simultaneous integrated boost or stereotactic boost was not utilized for any of the treatment plans. The radiation therapy oncology group (RTOG) trial 0933 study criteria ^{[4]} were followed for contouring and dosevolume constraints. Bilateral hippocampal contours were made on imaging datasets consisting of thinslice spoiled gradient echo (SPGR) magnetic resonance imaging sequences fused to a treatment planning computed tomography (CT) scan. The PTV was defined as the wholebrain parenchyma excluding the bilateral hippocampal contours plus a 5mm margin. The variation acceptable criteria for the protocol were used. In brief, these included PTV D2% ≤40 Gy and D98% ≥25 Gy. The hippocampal constraints included a D100% ≤10 Gy and a maximum hippocampal dose of ≤17 Gy. All treatment plans used identical contour sets and a single radiation oncologist reviewed and evaluated all treatment plans. Calculations of statistical significance were performed using Student's paired ttest to compare VMAT and IMRT plans. Treatment plans were evaluated based on tumor dose fall off, amount of dose received by normal brain tissue, and critical structures as well as the radiobiological comparison which was the purpose of this manuscript. This study was reviewed by our institutional IRB before initiation.
Equivalent uniform dose
While we strive to give tumors a homogenous dose in radiation therapy, in reality, they are inhomogeneous. To be able to better compare dose distributions, the concept of the equivalent uniform dose (EUD) was conceived. The EUD is the dose, which when homogenously given to a tumor gives, on average, the same number of surviving clonogens as a real clinical dose distribution gives. This opens a numerical way to be able to compare treatment plans from the DVH. Niemierko ^{[5]} defined the EUD mechanistically as follows:
where
In this formulation, N_{f} is the number of fractions in the treatment, D_{ref} is the reference dose of 2 Gy, the α/β is the ratio of two radiosensitivity parameters found in the linear quadratic model explaining the linear and quadratic parts of cell killing, SF_{2} is the clonogenic cell survival fraction at a dose of 2 Gy, D_{i} is the dose being received by each partial volume segment, and v_{i} and p_{i} are, respectively, the local volumes and densities of clonogens.
We calculated EUDs for the IMRT and VMAT treatment plans using the formula shown in equation 1. When calculating the EUDs, we varied the values of SF_{2 } as 0.4, 0.5, and 0.6 while varying the values of the α/β ratio to 2, 6, and 10 for each value of SF_{2} to determine appropriate model parameters for the future studies. Our value for N_{f} was 10 and the value for D_{ref} was 2 Gy. The DVH was divided into 10 partial volume segments each covering 10% of the structure volume. D_{i} was taken to be an average of the doses being received by the highest and lowest partial volumes. For example, when calculating the dose for the partial volume segment which spans 80% to 90% of the structure volume, we took D_{i} to be the average of the dose covering 80% of the structure volume and the dose covering 90% of the structure volume.
Niemierko ^{[6],[7]} also came up with the concept of a generalized EUD (gEUD) with an empirical formula:
where v_{i} is the fractional volume irradiated to a dose D_{i,} and “a” is a unitless tissuespecific parameter that describes the dosevolume effect. As the value of “a” becomes a large negative number, the gEUD leans toward the minimum dose received by the tissue. When “a = 1,” the gEUD is equal to the average dose. As the value of “a” becomes a large positive number, the value of the gEUD leans toward the maximum dose. For tissues which experience serial organ complication, the value of “a” should be a large positive number, for tissues that experience parallel organ complication, the value of “a” should be around 1, and for tumors, the value of “a” should be a large negative number. We also calculated gEUDs using formula shown in equation 2. D_{i}'s were reused from our previous EUD calculation. Because the DVH was divided into 10 equal segments, v_{i} was kept constant at 0.10. While calculating gEUDs for PTV, we varied our fitting parameter “a” to integers from −1 to −16. For calculating gEUDs for normal tissues, we used fitting parameter “a” equals +7.
Tumor control probability
The TCP is a quantitative measurement of local tumor control which can be calculated from information obtained from a DVH. While values of TCP calculated from formulas may not be the reflection of reality, the relative values calculated for different DVHs are useful. TCPs are calculated ^{[8]} with the following mechanistic formula:
where N is the number of clonogenic tumor cells (assumed in this study to be 4000, 40000, 400000, and 800000), D is EUD, N_{f}, D_{ref}, SF_{2}, and α/β are terms explained before. Radiation therapy is traditionally given in many fractions. This is because delivering the dose over time allows killing of maximum clonogenic tumor cells while having the least toxic effect on normal tissues. We calculated TCPs for each combination of N, SF_{2}, and α/β ratio using the formula shown in equation 3.
Another common way to calculate TCP is with the following empirical formula ^{[7]} shown in equation 4.
where TCD_{50} is the dose to control the tumor 50% of the time when the tumor is irradiated homogenously, gEUD is the gEUD calculated with formula shown in equation 2, and is a unitless fitting parameter that determines the slope of the doseresponse curve. We also calculated values of TCP using the formula shown in equation 4. We varied the value of parameter “a” in equation 2 from 1 to 16 to get values of gEUD. We used a TCD_{50} value ^{[9]} of 22.17 Gy and varied γ_{50} to determine value of TCP that agreed the best with our mechanistic model prediction.
Normal tissue complication probability
The NTCP is a quantitative measurement of the probability a dose of radiation will have an undesirable effect on an organ. The following mechanistic formula ^{[10]} is used to calculate the NTCP as shown in equation 5.
where
ϑ is the fraction of the organ irradiated, V_{ref} is the reference volume for TD_{50}, TD_{50} is the tolerance dose to the reference volume which would lead to a 50% complication probability, n is a parameter which determines the volume dependence of the complication probability, and m is a parameter which determines the slope of the complication probability versus the dose curve. Doseresponse models are commonly based on linearquadratic model of cell killing which has two adjustable parameters: α (linear component of cell killing) and β (the quadratic component of cell killing). The α/β ratio is the dose at which the linear and quadratic components of cell killing have equal contributions.^{[11]} The isoeffective dose is a method to account for the effect of fractionations in radiation therapy treatment that takes into account the behavior of early and lateresponding tissues. This is calculated with the following formula ^{[12]} shown in equation 6 relative to a standard 2 Gy per fractionation scheme usually used in conventional fractionated radiation therapy.
where D_{f} (=3 Gy) has been the fractionated dose used in the present treatment.
We compared the tissuesparing capabilities of IMRT and VMAT by calculating the NTCPs for various tissues that included among others hippocampus, optical structures, etc., from each of the corresponding treatment plans. We calculated NTCPs using the formula shown in equation 5. Microsoft Excel's normal distribution function was used to calculate NTCPs once tvalues were calculated with the formula shown in equation 5. Conventional fractionated radiation therapy usually uses doses of approximately 2 Gy/fraction and so typical value of the parameters TD_{50} used to calculate tvalues are available in the literature ^{[13]} for 2 Gy/fraction schemes. However, we used 3 Gy/fraction for treatment in the present study. Using the isoeffective dose value of 1.25 as shown in equation 6, TD_{50} values ^{[13]} of brain (since the TD_{50} value of hippocampus is not known) and optical structures of 60 and 65 are modified and values of the parameters TD_{50}, m, and n used to calculate tvalues are given in [Table 1].  Table 1: TD_{50}, m, and n values used to calculate normal tissue complication probability for each critical structure
Click here to view 
Another common method of calculating the NTCP for various tissues is through the following empirical formula ^{[9]} shown in equation 7.
where TD_{50} is the tolerance dose that gives a 50% complication rate, gEUD, and γ_{50} are parameters previously explained. We also calculated NTCPs using the above formula. Modified TD_{50} values shown in [Table 1] were used along with a constant γ_{50} value of 3.^{[14]} Secondary to limitations on the availability of hippocampusspecific parameter estimates, the same calculation end point of necrosis that was used for the brain and optic structures were also used for the hippocampus.
Results   
All VMAT and IMRT treatment plans met the RTOG trial 0933 criteria. The average EUDs calculated using equation 1 with SF_{2} of 0.4, 0.5, and 0.6 and α/β ratios of 2, 6, and 10 for IMRT and VMAT treatment plans ranged from 25.74–29.11 Gy for IMRT to 25.27–28.58 Gy for VMAT treatment plans. In all cases, the average EUD for IMRT was greater than that of VMAT by about 2% (P ≤ 0.02) as shown in [Figure 1]. In general, the average EUD increased with the α/β ratio and the SF_{2} for both IMRT and VMAT. Examples of typical EUD (with SF_{2}= 0.5 and α/β = 2) and gEUD (with a = −11) for the PTV for each patient are shown in [Figure 2]. The average EUD and gEUD for IMRT was again greater by 2% compared to that for VMAT. Furthermore, we found that the calculated gEUDs of the PTV using equation 2 decreased when the fitting parameter, “a,” decreased and the average gEUD was 2% higher in IMRT treatment plans than that of VMAT treatment plans when −16≤ a ≤ −2, and 3% when “a” = −1. In all cases, the difference between average IMRT and VMAT gEUDs were statistically significant (P ≤ 0.02).  Figure 1: Average equivalent uniform dose for various SF_{2}and α/β values for intensitymodulated radiation therapy and volumetricmodulated arc therapy treatment plans. Standard deviation is shown by error bars. All P ≤ 0.02
Click here to view 
 Figure 2: Comparison of intensitymodulated radiation therapy and volumetricmodulated arc therapy treatment plans using formulas for equivalent uniform dose and generalized equivalent uniform dose (generalized equivalent uniform dose). The average equivalent uniform dose for intensitymodulated radiation therapy versus volumetricmodulated arc therapy were 26.61 versus 26.16 and the average generalized equivalent uniform dose was 27.73 versus 27.12, respectively, with P ≤ 0.02. For both the equivalent uniform dose and the generalized equivalent uniform dose, the ratio of intensitymodulated radiation therapy/volumetricmodulated arc therapy was 1.02
Click here to view 
The TCPs for both IMRT and VMAT treatment plans were calculated using equation 3 with various SF_{2}, α/β, and estimated clonogenic cell counts. [Figure 3] shows results of a typical TCP calculation with SF_{2} of 0.5, α/β ratio of 2, and clonogenic cell counts estimated with CCD of 500 using formula 4. The ratio of average IMRT TCPs to average VMAT TCPs were also found to increase with the increase of SF_{2}, α/β ratio, and CCD values. We varied the value of γ_{50} in equation 4 and found that the TCP values were most agreeable with results obtained using the formula in equation 3 when γ_{50}= 3.6. Using equation 4, it was also found that the ratio of average IMRT TCPs over average VMAT TCPs varied from 1.00 to 1.03, as we varied the parameter “a” from −1 to −16. [Figure 3] also shows a typical TCP calculation with “a” equals11, γ_{50}= 3.6 and TCD_{50}= 22.17. Our results show that on the basis of TCP calculations, IMRT performed 2% better than VMAT using equation 3, and 1% better using equation 4; which again indicated better tumor control in favor of IMRT.  Figure 3: Comparison of tumor control probability values of intensitymodulated radiation therapy and volumetricmodulated arc therapy treatment plans. The average mechanistic tumor control probability was 0.92 for intensitymodulated radiation therapy and 0.90 for volumetricmodulated arc therapy, with an average ratio of intensitymodulated radiation therapy/volumetricmodulated arc therapy of 1.02. For the empirical tumor control probability, these values were 0.95 for intensitymodulated radiation therapy and 0.94 for volumetricmodulated arc therapy, respectively, with an average intensitymodulated radiation therapy/volumetricmodulated arc therapy ratio of 1.01
Click here to view 
With the mechanistic formula shown in equation 5, NTCPs were calculated for hippocampus, optic chiasm, right optic nerves, and left optic nerves. Specifically, IMRT lowered the NTCP of optic chiasm by a factor of two, right and left optic nerves by a factor of twentyfive but raised the NTCP of hippocampus by a factor of thirty. [Figure 4] contains the average gEUDs of critical structures for IMRT and VMAT treatment plans calculated with the formula shown in equation 2. The hippocampus in IMRT treatment plans had a higher gEUD than that of VMAT treatment plans (P = 0.19); however, the opposite was true for all the three optical structures (P < 0.01) which supports the previous trend obtained using the mechanistic formula. The average NTCPs calculated from gEUDs using equation 7 were greater for IMRT compared to VMAT for the hippocampus by about 27% (P = 0.13). The opposite was true for the optical structures, where the NTCPs calculated were lower for IMRT compared to VMAT by about 27% for optic chiasm, 62% for right optic nerve and 56% for left optic nerve (P < 0.01). The NTCP values per patient are shown in [Figure 5]. [Figure 6] gives representative examples of dose shown in colorwash superimposed on an axial slice from the treatment planning CT set for both IMRT and VMAT.  Figure 4: Generalized equivalent uniform doses of critical structures using parameter “a” = +7, for intensitymodulated radiation therapy and volumetricmodulated arc therapy treatment plans. The average generalized equivalent uniform dose in Gy for intensitymodulated radiation therapy versus volumetricmodulated arc therapy for the hippocampus, optic chiasm, right optic nerve, and left optic nerve were, respectively, (12.83 vs. 12.51), (27.64 vs. 28.38), (22.89 vs. 24.98), and (22.52 vs. 24.59)
Click here to view 
 Figure 5: Comparison of normal tissue complication probability values for intensitymodulated radiation therapy and volumetricmodulated arc therapy treatment plans for the hippocampus, optic chiasm, right optic nerve, and left optic nerve
Click here to view 
 Figure 6: Comparison of typical dose distribution for a patient planned with hippocampal avoidance wholebrain radiation therapyvolumetricmodulated arc therapy (left) and intensitymodulated radiation therapy (right) superimposed on axial slice of planning computed tomography. The volumetricmodulated arc therapy plan has a lower maximum dose than that of the intensitymodulated radiation therapy treatment plan
Click here to view 
Discussion   
In this manuscript, we estimate TCP, NTCP, and EUD values for IMRT and VMATbased treatment planning for hippocampal sparring WBRT. Our analysis suggests that both of these treatment approaches have their advantages. On the basis of TCP calculations, IMRT performed 2% better than VMAT using equation 3 and 1% better using equation 4. Both of these results indicate improved TCP in favor of IMRTbased treatments. Similarly, IMRTbased treatment plans resulted in lower NTCP values for the optic chiasm and optic nerves while higher NTCP parameters were obtained for the hippocampus. These results are consistent with reports suggesting VMAT may be better at sparing the hippocampus while meeting the treatment goals for the PTV.^{[15],[16]} Even though the relative difference in NTCP was large, the values for the NTCP's were so small that any rational statement based on this comparison is difficult to make.
In this manuscript, the parameters for empirical TCP and NTCP models were derived by fitting our data to mechanistic models. Values for our mechanistic models were taken from values in literature derived from clinical outcomes and tuned to derive optimal results. While these fitting parameters are not the most optimal, we feel they are appropriate for use. TCP calculations have been utilized for the treatment of brain tumors, both in the setting of stereotactic radiation therapy as well as with standard fractionation radiation therapy.^{[17],[18],[19]} The parameters used in our model are consistent with what has been used in these reports. However, there is little information specific to the setting of HAWBRT with IMRT or VMAT, where the overall dose distribution is more heterogeneous, and our results represent best estimates for this setting. In particular, little information specific to the hippocampus is available. For our NTCP calculations, we used similar parameters for the brain and the hippocampus, while the potential toxicity, as well as the dose at which the toxicity occurs, is very different, both clinically and mechanistically. The toxicity of interest for the hippocampus is the development of neurocognitive function decline which is likely to occur at lower doses than what would be expected for brain necrosis and further studies are necessary to establish TD_{5/5} and TD_{50/5} estimates that are specific for the hippocampus.
Conclusions   
HAWBRT is an effective treatment to the brain with its goal being to provide adequate coverage of the brain parenchyma while reducing the radiation doses to the hippocampus. The linacbased IMRT or VMAT are potential techniques to reduce doses to the hippocampus while treating the brain parenchyma. We found that IMRT treatment plans had a higher TCP than VMAT treatment plans and also had a lower NTCP for the optical structures. For the hippocampus, when the end point for the calculation was necrosis, IMRT plans had higher NTCP values when compared to VMAT treatment plans. In this study, we have also found appropriate fitting parameters that can be used in empirical models for EUD, TCP, and NTCP. The concept of isoeffective dose to estimate TD50 for hippocampus and optical structures with 3 Gy per fraction used in this study could be applied to any radiation therapy study for other organs and cancers utilizing different fractionation scheme.
Financial support and sponsorship
Nil.
Conflicts of interest
There are no conflicts of interest.
References   
1.  Khuntia D, Brown P, Li J, Mehta MP. Wholebrain radiotherapy in the management of brain metastasis. J Clin Oncol 2006;24:1295304. 
2.  Marsh JC, Gielda BT, Herskovic AM, Abrams RA. Cognitive sparing during the administration of whole brain radiotherapy and prophylactic cranial irradiation: Current concepts and approaches. J Oncol 2010;2010:198208. 
3.  Allen Li X, Alber M, Deasy JO, Jackson A, Ken Jee KW, Marks LB, et al. The use and QA of biologically related models for treatment planning: Short report of the TG166 of the therapy physics committee of the AAPM. Med Phys 2012;39:1386409. 
4.  Gondi V, Pugh SL, Tome WA, Caine C, Corn B, Kanner A, et al. Preservation of memory with conformal avoidance of the hippocampal neural stemcell compartment during wholebrain radiotherapy for brain metastases (RTOG 0933): A phase II multiinstitutional trial. J Clin Oncol 2014;32:38106. 
5.  Niemierko A. Reporting and analyzing dose distributions: A concept of equivalent uniform dose. Med Phys 1997;24:10310. 
6.  Niemierko A. A generalized concept of equivalent uniform dose (EUD). Med Phys 1999;26:1101. 
7.  Gay HA, Niemierko A. A free program for calculating EUDbased NTCP and TCP in external beam radiotherapy. Phys Med 2007;23:11525. 
8.  Ahmad S, Vogds BJ, McKenna F, Vlachaki MT. Tumor control probability (TCP) in prostate cancer: Role of radiobiological parameters and radiation dose escalation. J Xray Sci Technol 2009;17:34754. 
9.  Okunieff P, Morgan D, Niemierko A, Suit HD. Radiation doseresponse of human tumors. Int J Radiat Oncol Biol Phys 1995;32:122737. 
10.  Lyman JT. Complication probability as assessed from dosevolume histograms. Radiat Res Suppl 1985;8:S139. 
11.  Hall EJ, Giaccia AJ. Radiobiology for the Radiologist. Hagerstone, MD: Lippincott Williams & Wilkins; 2012. 
12.  Metcalfe P, Kron T, Hoban P. The Physics of Radiotherapy Xrays and Electrons. Madison, WI: Medical Physics Publishing Corporation; 2007. p. 450. 
13.  MacDonald SM, Ahmad S, Kachris S, Vogds BJ, DeRouen M, Gittleman AE, et al. Intensity modulated radiation therapy versus threedimensional conformal radiation therapy for the treatment of high grade glioma: A dosimetric comparison. J Appl Clin Med Phys 2007;8:4760. 
14.  ClementeGutiérrez F, PérezVara C, ClavoHerranz MH, LópezCarrizosa C, PérezRegadera J, IbáñezVilloslada C, et al. Assessment of radiobiological metrics applied to patientspecific QA process of VMAT prostate treatments. J Appl Clin Med Phys 2016;17:34167. 
15.  Wang BH, Hua W, Gu X, Wang XL, Li J, Liu LQ, et al. Dosimetric study of different radiotherapy planning approaches for hippocampal avoidance wholebrain radiation therapy (HAWBRT) based on fused CT and MRI imaging. Australas Phys Eng Sci Med 2015;38:76775. 
16.  Rong Y, Evans J, XuWelliver M, Pickett C, Jia G, Chen Q, et al. Dosimetric evaluation of intensitymodulated radiotherapy, volumetric modulated arc therapy, and helical tomotherapy for hippocampalavoidance whole brain radiotherapy. PLoS One 2015;10:e0126222. 
17.  Pedicini P, Fiorentino A, Simeon V, Tini P, Chiumento C, Pirtoli L, et al. Clinical radiobiology of glioblastoma multiforme: Estimation of tumor control probability from various radiotherapy fractionation schemes. Strahlenther Onkol 2014;190:92532. 
18.  TomaDasu I, Sandström H, Barsoum P, Dasu A. To fractionate or not to fractionate? That is the question for the radiosurgery of hypoxic tumors. J Neurosurg 2014;121Suppl:1105. 
19.  Shuryak I, Carlson DJ, Brown JM, Brenner DJ. Highdose and fractionation effects in stereotactic radiation therapy: Analysis of tumor control data from 2965 patients. Radiother Oncol 2015;115:32734. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6]
[Table 1]
