IsoBED: a tool for automatic calculation of biologically equivalent fractionation schedules in radiotherapy using IMRT with a simultaneous integrated boost (SIB) technique
© Bruzzaniti et al; licensee BioMed Central Ltd. 2011
Received: 24 January 2011
Accepted: 9 May 2011
Published: 9 May 2011
An advantage of the Intensity Modulated Radiotherapy (IMRT) technique is the feasibility to deliver different therapeutic dose levels to PTVs in a single treatment session using the Simultaneous Integrated Boost (SIB) technique. The paper aims to describe an automated tool to calculate the dose to be delivered with the SIB-IMRT technique in different anatomical regions that have the same Biological Equivalent Dose (BED), i.e. IsoBED, compared to the standard fractionation.
Based on the Linear Quadratic Model (LQM), we developed software that allows treatment schedules, biologically equivalent to standard fractionations, to be calculated. The main radiobiological parameters from literature are included in a database inside the software, which can be updated according to the clinical experience of each Institute. In particular, the BED to each target volume will be computed based on the alpha/beta ratio, total dose and the dose per fraction (generally 2 Gy for a standard fractionation). Then, after selecting the reference target, i.e. the PTV that controls the fractionation, a new total dose and dose per fraction providing the same isoBED will be calculated for each target volume.
The IsoBED Software developed allows: 1) the calculation of new IsoBED treatment schedules derived from standard prescriptions and based on LQM, 2) the conversion of the dose-volume histograms (DVHs) for each Target and OAR to a nominal standard dose at 2Gy per fraction in order to be shown together with the DV-constraints from literature, based on the LQM and radiobiological parameters, and 3) the calculation of Tumor Control Probability (TCP) and Normal Tissue Complication Probability (NTCP) curve versus the prescribed dose to the reference target.
Irradiation techniques with Intensity Modulated Radiotherapy (IMRT) allow doses to be delivered to the target with a high conformation of prescribed isodose, sparing Organs at Risk (OARs), compared to conventional 3D-CRT techniques. Another advantage of the IMRT technique is the possibility to achieve the so-called Simultaneous Integrated Boost (SIB), which provides different levels of therapeutic doses to different target volumes during the same treatment session, once the fraction number has been set [1–5].
Historically, to obtain the desired tumor control, the doses were determined using a conventional fractionation that ranged between 50 to 70 Gy at 2 Gy per fraction.
Whereas, in order to obtain Tumor Control Probability (TCP), equivalent to that of a conventional fractionation, the total dose simultaneously delivered to the targets have to be determined according to the Linear Quadratic Model (LQM) to be used with the SIB technique . Thus, the dose per fraction to PTVs and/or boost may differ by 2 Gy per fraction.
The paper aims to: 1) describe home-made software, based on the IsoBED formula, able to calculate the total dose and the dose per fraction with the same TCP as the conventional fractionation, that will be used with the SIB technique, 2) import the DVHs from different TPSs or different plans, convert them into a normalized 2 Gy-fraction-Volume Histogram (NTD2-VH) and compare these amongst themselves and with the Dose-Volume constraints (DV- constraints), 3) calculate and compare the TCPs and the Normal Tissue Complication Probabilities (NTCPs) obtained from different DVHs.
and the (α/β) ratio is a parameter which takes into account the radiobiological effect of fractionation in tumor or OARs.
Equation (2) is the basis on which a comparison of different treatment strategies is performed.
Where d2 is the new dose per fraction delivered in n fractions, resulting in a new total dose D2 = d2n,
Equation (7) is valid for both PTVs and OARs (following the LQM).
The IsoBED software
The software has been developed using the Microsoft Visual Basic 6.0. The main form - the IsoBED Calculator- gives a choice between IsoBED calculation and DVHs analysis modules.
The software allows the anatomical district to be selected. The user has to introduce the total dose, dose per fraction (generally 2 Gy per fraction) for each target (up to 3) and, the (α/β) ratio of investigated tumor must be inserted to calculate the corresponding BED.
Then the software requires the selection of the reference target (which determines the fractions number in the SIB treatment), in order to calculate the new fractionation for the remaining targets, based on equation (7). Furthermore, the software permits a comparison of the biologically equivalent schedules using hyper/hypo-fractionated as well as conventional regimes. It also includes a database with the main DV- constraints at 2 Gy per fraction for different OARs derived from literature and clinical experience in the radiotherapy department of our Institute [9–20] which may be upgraded by the user.
The DV-constraints are converted to those of the new schedule (i.e. hypo or hyper-fractionated) calculated by IsoBED.
Then the converted constraints for OARs can be printed and used as constraints for IMRT optimization.
DVH import and radiobiological analysis
After the IMRT optimization using commercial TPSs (such as: BrainScan, Eclipse, Pinnacle), the obtained DVHs can be imported to our software and can be used to compare techniques and/or dose distributions from the same or different TPSs.
The software automatically recognizes the DVH file format exported from each TPS source and imports it into the patient directory without any changes. In particular, import procedures consist of copying DVH files into a subfolder with the patient's name, contained in a directory where the IsoBED.exe file is held.
Then, a specific window permits the analysis of DVHs to be carried-out. Cumulative or differential DVHs can be visualized after setting dose per fraction and fraction number. In this window up to five plans imported from BrainScan, Eclipse and Pinnacle can be compared. The volumes and the minimum, mean, median, modal and maximum doses can be visualized for OARs and PTVs.
For each volume the software calculates NTD2VH (Appendix 1 equation 1.6) by using the appropriate (α/β)ratio, which may be changed by the user.
Finally, the TCP, NTCP and Therapeutic Gain (P+) curves can be calculated from the DVHs based on radiobiological parameter sets, derived from literature but upgraded by the user, according to the formulas reported in Appendix 1 [21–27].
To illustrate this user friendly IsoBED software some case examples are shown.
The following test cases were considered in order to illustrate the usefulness of the home made software for comparing sequential versus SIB plans for three clinical treatments in this paper.
The first case regards irradiation using IMRT of prostate and pelvic lymph nodes.
The comparison was made between the sum of 2 sequential IMRT plans (50 Gy to the lymph nodes and prostate at 2 Gy per fraction followed by another 30 Gy at 2 Gy per fraction only on the prostate for a total of 40 fractions) and an SIB IMRT plan .
Assuming the same fractionation for prostate, the total dose and dose per fraction of pelvic lymph nodes were calculated with the IsoBED software, using an (α/β)ratio = 1.5 Gy for both targets [28, 29].
The treatment plans were developed using Helios module of Eclipse TPS (Varian Medical System). All 3 treatment plans were performed with the same geometry using 5 coplanar fields (angles: 0, 75, 135, 225 and 285 degrees) with the patient in prone position.
Head and Neck Case
The second case regards the treatment of a rinopharynx cancer patient.
The prescribed dose was 53 Gy at 2.12 Gy per fraction to the Planning Elective Tumor Volume (PETV, i.e. PTV54), 59.36 Gy at 2.12 Gy per fraction to the Planning Clinical Target Volume (PCTV, i.e. PTV60) and 69.96 Gy at 2.12 Gy per fraction to the Planning Gross Target Volume (PGTV, i.e. PTV70).
The first plan, the sequential treatment, was calculated to deliver 53 Gy in 25 fractions to PETV followed by 6.36 Gy in 3 fractions to the PCTV and another 10.6 Gy in 5 fractions to the PGTV, for a total of 33 fractions.
For the SIB plan, the IsoBED doses derived from prescription and the calculated doses from our software were considered in order to deliver 69.96 Gy in 33 fractions to the PGTV.
In a lung cancer patient two volumes had to be irradiated in a hypofractionaction regime . The prescription of the sequential technique was: PTV to receive 40 Gy at 10 Gy per fraction and for the boost an additional fraction of 10 Gy. The SIB technique consisted of an IMRT plan, for which the dose were calculated by IsoBED software, so that the boost received 50 Gy in 5 fractions.
The plan sum was created from the sequential IMRT plans which had to be compared with the IMRT SIB plan. All plans were exported from TPSs and imported into the IsoBED software to calculate and compare NTD2VH, TCP, NTCP and P+.
Then the SIB schedule is calculated by selecting the control box "IsoBED Calculation". The results of such evaluations are visualized in the "IsoBED DOSES" area. The dose limits are visualized in the "OAR CONSTRAINTS" area.
Import procedures consist of copying DVH files, exported from TPS, in a folder with the patient's name contained in a directory where an IsoBED.exe file is installed. DVH files are different depending on the TPS source. IsoBED can import DHV data files from Eclipse, Pinnacle and Brainscan.
Dose distribution and radiobiological analysis
Figures 5a, 5b and 5c show the DVHs imported from TPSs calculated with different modalities (SIB and sequential). The user can choose which volume of interest to view by selecting them from a list visualized at the lower-left corner of the screen. Furthermore, in the same area, the total volume or one between, the minimum, maximum, average, median and modal dose percentage for each plan and each structure shown in the histogram is displayed.
In order to perform radiobiological calculations the (α/β)values can be set for each structure by choosing a dropdown menu in which the list of parameters incorporated in a dedicated database appears. These values are derived from literature data and from experience at our Institute [9–20]. The "NTD2" button transforms every DVH into the NTD2VH (Figures 6a, 6b and 6c).
Finally, the TCP, NTCP and P+ curves against the dose prescribed to the reference target can be calculated with the "TCP-NTCP" button and their values are shown in the lower area of the screen (Figures 7a, 7b and 7c).
All the outcomes from IsoBED software were compared with an automatic excel spreadsheet specially designed for this purpose. In particular, the outcomes from IsoBED calculation and from DVH import and radiobiological analysis modules were tested. The results obtained from the comparison made it possible to validate the software.
The introduction of the IMRT technique in clinical practice, including the SIB approach, requires new treatment schedules able to guarantee the same BED of conventional fractionations to be drawn up. Automatic software that does this is a useful tool when making these estimates, particularly with regard to evaluations and for comparing different forms of DVHs and radiobiological parameters [30–35].
The software, described in this paper, is based on the BED calculation and on LQM. Unlike other software, it allows fractionation schedules to be calculated in SIB-IMRT treatment techniques with both conventional and hypo-fractionation regimes, after setting the desired dose per fraction.
Similar to Bioplan , the IsoBED software is an analysis tool used to compare DVHs with different TPSs or different irradiation techniques.
In addition, this software allows a comparison between plans using NTD2VH. This is a very interesting and useful aspect as it is possible to take into consideration simultaneously the end-points of different OARs.
Moreover, the import of DVHs enables dosimetric and radiobiological comparisons between different TPSs, which is an important issue because this may be used as quality control for treatment planning systems when simple geometry of phantoms are assumed [36, 37].
In addition, the TCP and NTCP curves can be calculated to select the best treatment plans to be discussed with physicians. In fact, the P+ curve can be used to confirm the dose prescription to reference target. In particular, the maximum peak of the P+ curve indicates the dose per fraction to reference target giving the maximum TCP value with the lowest combination of NTCPs.
Furthermore, the possibility of changing the (α/β)value while designing the fractionation scheme might aid the prediction of different effects (such as acute and late effect) related to clinical trials.
Finally, the possibility of updating the radiobiological parameters for OARs stored in the internal database permits us to take into consideration the proven clinical experience of users. The software calculates the radiobiological DV-constrains for different fractionations as shown in the case examples (Figure 1, 2 and 3).
An issue to be considered regards the use of the LQM adopted by IsoBED. In fact, this model is strictly applicable with intermediate doses while its applicability with doses higher than 18-20 Gy per fraction is under debate [38, 39]. Nevertheless, the use of simple analytic models may provide useful suggestions in clinical radiotherapy.
IsoBED software based on LQM allows one to design treatment schedules by using the SIB approach, importing DVHs from different TPSs for dosimetric and radiobiological comparison. It also allows to select and evaluate the best approach able to guarantee maximum TCP and at the same time the minimum NTCP to the organs at risk.
where N* is the total initial number of tumor clonogenic cells and sf is the surviving fraction.
where m and TD50 (v eff ) are the slope of the NTCP curve versus the dose and the tolerance dose at 2 Gy per fraction to a fraction v eff of the organ, respectively.
In order to generalize the LBK method each DVH has been converted into a single value using a DVH reduction method.
where Di is the dose delivered to the volume fraction v i , K is the number of points of the differential DVH, D max is the maximum dose and n is a parameter related to organ response to radiation (n = 0,1 for serial and parallel organs, respectively). By Eq. (1.4), an inhomogeneous dose distribution is converted into an equivalent uniform irradiation of a fraction v eff of the organ treated at the maximum dose (D max ).
where TD50(1) is the tolerance dose to the whole organ, leading to a 50% complication probability.
By using this formula, each dose step in the DVHs was corrected separately. This formalism presumes complete cellular repair between treatment fractions and neglects the role of cellular re-population. The latter assumption is valid for late-responding normal tissues but is inaccurate for acute-responding tissues and tumors. This limitation may be important when using the LQM to compare treatment schedules differing in overall treatment times in terms of their acute effects (for which time-dependent repopulation may be important). For late effects, time factors are generally thought to be of minor importance.
The Authors wish to thank Mrs. Paula Franke for the English revision of the manuscript.
- Ang KK, Peters LJ: Concomitant boost radiotherapy in the treatment of head and neck cancer. Semin Radiat Oncol. 1992, 2: 31-33. 10.1016/S1053-4296(05)80049-9.View ArticleGoogle Scholar
- Ang KK, Peters LJ, Weber RS: Concomitant boost radiotherapy schedules in the treatment of carcinoma of the oropharynx and nasopharynx. Int J Radiat Oncol Biol Phys. 1990, 19: 1339-1345.View ArticlePubMedGoogle Scholar
- Mohan R, Wu Q, Manning M, Schmidt-Ullrich R: Radiobiological considerations in the design of fractionation strategies for intensity-modulated radiation therapy of head and neck cancers. Int J Radiat Oncol Biol Phys. 2000, 46 (3): 619-630.View ArticlePubMedGoogle Scholar
- Dogan N, King S, Emami B, Mohideen N, Mirkovic N, Leybovich LB, Sethi A: Assessment of different IMRT boost delivery methods on target coverage and normal-tissue sparing. Int J Radiat Oncol Biol Phys. 2003, 57: 1480-1491.View ArticlePubMedGoogle Scholar
- Fogliata A, Bolsi A, Cozzi L, Bernier J: Comparative dosimetric evaluation of the simultaneous integrated boost with photon intensity modulation in head and neck cancer patients. Radiother Oncol. 2003, 69: 267-275.View ArticlePubMedGoogle Scholar
- Strigari L, D'Andrea M, Abate A, Benassi M: A heterogeneous dose distribution in simultaneous integrated boost: the role of the clonogenic cell density on the tumor control probability. Phys Med Biol. 2008, 53: 5257-5273.View ArticlePubMedGoogle Scholar
- Stavrev P, Hristov D: Prostate IMRT fractionation strategies: two-phase treatment versus simultaneous integrated boost. Radiol Oncol. 2003, 37: 115-126.Google Scholar
- Mohan R, Wu Q, Manning M, Schmidt-Ullrich R: Radiobiological considerations in the design of fractionation strategies for intensity-modulated radiation therapy of head and neck cancers. Int J Radiat Oncol Biol Phys. 2000, 46: 619-630.View ArticlePubMedGoogle Scholar
- 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-122.View ArticlePubMedGoogle Scholar
- Strigari L, Arcangeli G, Arcangeli S, Benassi M: Mathematical model for evaluating incidence of acute rectal toxicity during conventional or hypofractionated radiotherapy courses for prostate cancer. Int J Radiat Oncol Biol Phys. 2009, 73: 1454-1460.View ArticlePubMedGoogle Scholar
- Marzi S, Arcangeli G, Saracino B, Petrongari MG, Bruzzaniti V, Iaccarino G, Landoni V, Soriani A, Benassi M: Relationships between rectal wall dose-volume constraints and radiobiologic indices of toxicity for patients with prostate cancer. Int J Radiat Oncol Biol Phys. 2007, 68: 41-49.View ArticlePubMedGoogle Scholar
- Rancati T, Fiorino C, Gagliardi G, Cattaneo GM, Sanguineti G, Borca VC, Cozzarini C, Fellin G, Foppiano F, Girelli G, Menegotti L, Piazzolla A, Vavassori V, Valdagni R: Fitting late rectal bleeding data using different NTCP models: results from an Italian multi-centric study (AIROPROS0101). Radiother Oncol. 2004, 73: 21-32.View ArticlePubMedGoogle Scholar
- Abate A, Pressello MC, Benassi M, Strigari L: Comparison of IMRT planning with two-step and one-step optimization: a strategy for improving therapeutic gain and reducing the integral dose. Phys Med Biol. 2009, 54 (23): 7183-98.View ArticlePubMedGoogle Scholar
- Strigari L, Benassi M, Arcangeli G, Bruzzaniti V, Giovinazzo G, Marucci L: A novel dose constraint to reduce xerostomia in head-and-neck cancer patients treated with intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys. 2010, 77: 269-276.View ArticlePubMedGoogle Scholar
- Marzi S, Iaccarino G, Pasciuti K, Soriani A, Benassi M, Arcangeli G, Giovinazzo G, Benassi M, Marucci L: Analysis of salivary flow and dose-volume modeling of complication incidence in patients with head-and-neck cancer receiving intensity-modulated radiotherapy. Int J Radiat Oncol Biol Phys. 2009, 73: 1252-1259.View ArticlePubMedGoogle Scholar
- Eisbruch A, Ten Haken RK, Kim HM, Marsh LH, Ship JA: Dose, volume, and function relationships in parotid salivary glands following conformal and intensity-modulated irradiation of head and neck cancer. Int J Radiat Oncol Biol Phys. 1999, 45: 577-587.View ArticlePubMedGoogle Scholar
- Chao KS, Deasy JO, Markman J, Haynie J, Perez CA, Purdy JA, Low DA: A prospective study of salivary function sparing in patients with head-and-neck cancers receiving intensity-modulated or three-dimensional radiation therapy: initial results. Int J Radiat Oncol Biol Phys. 2001, 49: 907-916.View ArticlePubMedGoogle Scholar
- Mirri MA, Arcangeli G, Benassi M, d'Angelo A, Pinzi V, Caterino M, Rinaldi M, Ceribelli A, Strigari L: Hypofractionated Conformal Radiotherapy (HCRT) for Primary and Metastatic Lung Cancers with Small Dimension. Strahlenther Onkol. 2009, 185: 27-33.View ArticlePubMedGoogle Scholar
- Theuws JC, Kwa SL, Wagenaar AC, Seppenwoolde Y, Boersma LJ, Damen EM, Muller SH, Baas P, Lebesque JV: Prediction of overall pulmonary function loss in relation to the 3-D dose distribution for patients with breast cancer and malignant lymphoma. Radiother Oncol. 1998, 49: 233-243.View ArticlePubMedGoogle Scholar
- Kwa SL, Lebesque JV, Theuws JC, Marks LB, Munley MT, Bentel G, Oetzel D, Spahn U, Graham MV, Drzymala RE, Purdy JA, Lichter AS, Martel MK, Ten Haken RK: Radiation pneumonitis as a function of mean lung dose: an analysis of pooled data of 540 patients. Int J Radiat Oncol Biol Phys. 1998, 42: 1-9.View ArticlePubMedGoogle Scholar
- Marks Lawrence, Yorke Ellen, Jackson Andrew, Ten Haken Randall, Constine Louis, Eisbruch Avraham, Bentzen Søren, Nam Jiho, Deasy Joseph: Use of Normal Tissue Complication Probability Models in the Clinic. Int J Radiat Oncol Biol Phys. 2010, 76 (3): Supplement 1: S10-S19.Google Scholar
- Deasy J: Poisson formulas for tumor control probability with clonogenic proliferation. Radiat Res. 1996, 145: 382-384.View ArticlePubMedGoogle Scholar
- Lyman JT: Complication probability as assessed from dose-volume histograms. Radiat Res Suppl. 1985, 8: S13-19.View ArticlePubMedGoogle Scholar
- Kutcher GJ, Burman C: Calculation of complication probability factors for non-uniform normal tissue irradiation: the effective volume method. Int J Radiat Oncol Biol Phys. 1989, 16: 1623-1630.View ArticlePubMedGoogle Scholar
- Burman C, Kutcher GJ, Emami B, Goitein M: Fitting of normal tissue tolerance data to an analytic function. Int J Radiat Oncol Biol Phys. 1991, 21: 123-135.View ArticlePubMedGoogle Scholar
- Ågren A, Brahme A, Turesson I: Optimization of uncomplicated control for head and neck tumors. Int J Radiat Oncol Biol Phys. 1990, 19: 1077-1085.View ArticlePubMedGoogle Scholar
- Kallman P, Agren A, Brahme A: Tumour and normal tissue responses to fractionated non-uniform dose delivery. Int J Radiat Biol. 1992, 62: 249-262.View ArticlePubMedGoogle Scholar
- Fowler J: The radiobiology of prostate cancer including new aspects of fractionated radiotherapy. Acta Oncol. 2005, 44: 265-276.View ArticlePubMedGoogle Scholar
- Fowler JF, Chappell RJ, Ritter MA: Is α/β for prostate tumors really low?. Int J Radiat Oncol Biol Phys. 2001, 50: 1021-1031.View ArticlePubMedGoogle Scholar
- Sanchez-Nieto B, Nahum AE: BIOPLAN: software for the biological evaluation of radiotherapy treatment plans. Med Dosim. 2000, 25: 71-76.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- El Naqa I, Suneja G, Lindsay PE, Hope AJ, Alaly JR, Vicic M, Bradley JD, Apte A, Deasy JO: Dose response explorer: an integrated open-source tool for exploring and modelling radiotherapy dose-volume outcome relationships. Phys Med Biol. 2006, 51: 5719-5735.View ArticlePubMedGoogle Scholar
- Deasy JO, Blanco AI, Clark VH: CERR: a computational environment for radiotherapy research. Med Phys. 2003, 30: 979-985.View ArticlePubMedGoogle Scholar
- Gay HA, Niemierko A: A free program for calculating EUD-based NTCP and TCP in external beam radiotherapy. Phys Med. 2007, 23: 115-125.View ArticlePubMedGoogle Scholar
- Pyakuryal A, Myint WK, Gopalakrishnan M, Jang S, Logemann JA, Mittal BB: A computational tool for the efficient analysis of dose-volume histograms for radiation therapy treatment plans. J Appl Clin Med Phys. 2010, 11: 137-157.Google Scholar
- Ezzell GA, Galvin JM, Low D, Palta JR, Rosen I, Sharpe MB, Xia P, Xiao Y, Xing L, Yu CX: Guidance document on delivery, treatment planning, and clinical implementation of IMRT: Report of the IMRT subcommittee of the AAPM radiation therapy committee. Med Phys. 2003, 30: 2089-2115.View ArticlePubMedGoogle Scholar
- Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, Van Dyke J: American Association of Physicists in Medicine Radiation Therapy Committee Task Group 53: Quality assurance for clinical radiotherapy treatment planning. Med Phys. 1998, 25: 1773-1829.View ArticlePubMedGoogle Scholar
- Park C, Papiez L, Zhang S, Story M, Timmerman RD: Universal survival curve and single fraction equivalent dose: useful tools in understanding potency of ablative radiotherapy. Int J Radiat Oncol Biol Phys. 2008, 70: 847-52.View ArticlePubMedGoogle Scholar
- Fowler JF: Linear quadratics is alive and well: in regard to Park et al. (Int J Radiat Oncol Biol Phys 2008;70:847-852. Int J Radiat Oncol Biol PhysPhys. 2008, 72: 957-10.1016/j.ijrobp.2008.06.1929.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.