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.

Treatment planning tools that use biologically related models for plan optimization and/or evaluation are being introduced for clinical use. A variety of dose-response models and quantities along with a series of organ-specific model parameters are included in these tools. However, due to various limitations, such as the limitations of models and available model parameters, the incomplete understanding of dose responses, and the inadequate clinical data, the use of biologically based treatment planning system (BBTPS) represents a paradigm shift and can be potentially dangerous. There will be a steep learning curve for most planners. The purpose of this task group is to address some of these relevant issues before the use of BBTPS becomes widely spread. In this report, the authors (1) discuss strategies, limitations, conditions, and cautions for using biologically based models and parameters in clinical treatment planning; (2) demonstrate the practical use of the three most commonly used commercially available BBTPS and potential dosimetric differences between biologically model based and dose-volume based treatment plan optimization and evaluation; (3) identify the desirable features and future directions in developing BBTPS; and (4) provide general guidelines and methodology for the acceptance testing, commissioning, and routine quality assurance (QA) of BBTPS.

Improved critical structure sparing with biologically based IMRT optimization.

The impact of using biological models in treatment planning on plan quality is studied by comparing IMRT plans generated using selected commercially available treatment planning systems (TPSs) employing biological models/quantities in IMRT optimization (bIMRT) and the conventional physically (dose-volume) based optimization (pIMRT). A total of 25 IMRT plans, generated for five cases of different anatomic sites (brain, head and neck, lung, pancreas, and prostate) using five TPSs, two bIMRT (CMS Monaco and Phillips Pinnacle3 P3IMRT) and three pIMRT (CMS Xio, Phillips Pinnacle3, and Tomotherapy) systems, were compared. Dose-volume histograms, maximum, minimum, and mean doses, target heterogeneity and conformity indices, equivalent uniform dose (EUD), and an overall plan-ranking index (fEUD) were used in the comparison. It is clear from the comparison that the use of biological models in treatment planning optimization can generate IMRT plans with significantly improved normal tissue sparing with similar or slightly increased dose heterogeneity in the target, as compared to the conventional dose-volume based optimization for the same beam arrangement. For example, the bIMRT plans lead to smaller EUDs in 32 out of 37 normal structures in all five cases combined, as compared to the pIMRT plans. Caution should be exercised in choosing appropriate models and/or model parameters and in evaluating the plan obtained when using the biologically based treatment planning system.

Irradiation of varying volumes of rat lung to same mean lung dose: a little to a lot or a lot to a little?

PURPOSE: To investigate whether irradiating small lung volumes with a large dose or irradiating large lung volumes with a small dose, given the same mean lung dose (MLD), has a different effect on pulmonary function in laboratory animals. METHODS AND MATERIALS: WAG/Rij/MCW male rats were exposed to single fractions of 300 kVp X-rays. Four treatments, in decreasing order of irradiated lung volume, were administered: (1) whole lung irradiation, (2) right lung irradiation, (3) left lung irradiation, and (4) irradiation of a small lung volume with four narrow beams. The irradiation times were chosen to accumulate the same MLD of 10, 12.5, or 15 Gy with each irradiated lung volume. The development of radiation-induced lung injury for < or =20 weeks was evaluated as increased breathing frequency, mortality, and histopathologic changes in the irradiated and control rats. RESULTS: A significant elevation of respiratory rate, which correlated with the lung volume exposed to single small doses (> or =5 Gy), but not with the MLD, was observed. The survival of the rats in the whole-lung-irradiated group was MLD dependent, with all events occurring between 4.5 and 9 weeks after irradiation. No mortality was observed in the partial-volume irradiated rats. CONCLUSIONS: The lung volume irradiated to small doses might be the dominant factor influencing the loss of pulmonary function in the rat model of radiation-induced lung injury. Caution should be used when new radiotherapy techniques that result in irradiation of large volumes of normal tissue are used for the treatment of lung cancer and other tumors in the thorax.

Combined use of Monte Carlo DNA damage simulations and deterministic repair models to examine putative mechanisms of cell killing.

A kinetic repair-misrepair-fixation (RMF) model is developed to better link double-strand break (DSB) induction to reproductive cell death. Formulas linking linear-quadratic (LQ) model radiosensitivity parameters to DSB induction and repair explicitly account for the contribution to cell killing of unrejoinable DSBs, misrepaired and fixed DSBs, and exchanges formed through intra- and intertrack DSB interactions. Information from Monte Carlo simulations is used to determine the initial yields and complexity of DSBs formed by low- and high-LET radiations. Our analysis of published survival data for human kidney cells suggests that intratrack DSB interactions are negligible for low-LET radiations but increase rapidly with increasing LET. The analysis suggests that no class of DSB is intrinsically unrejoinable or that DSB reparability is not strictly determined by the number of lesions forming the DSB. For radiations with LET >110 keV/mum, the model predicts that the relative cell killing efficiency, per unit absorbed dose, should continue to increase, whereas data from published experiments indicate a reduced cell killing efficiency. This observation suggests that the Monte Carlo simulation overestimates the DSB yield beyond 110 keV/microm or that other biological phenomena not included in the model, such as proximity effects, are important. For 200-250 kVp X rays ( approximately 1.9 keV/microm), only about 1% of the one-track killing is attributed to intratrack binary misrepair interactions. The analysis indicates that the remaining 99% of the lethal damage is due to other types of one-track damage, including possible unrepairable, misrepaired and fixed damage. Compared to the analysis of the X-ray results, 48% of the one-track lethal damage caused by 5.1 MeV alpha particles (approximately 88 keV/microm) is due to intratrack DSB interactions while the remainder is due to other forms of one-track damage.

Evaluation of a commercial biologically based IMRT treatment planning system.

A new inverse treatment planning system (TPS) for external beam radiation therapy with high energy photons is commercially available that utilizes both dose-volume-based cost functions and a selection of cost functions which are based on biological models. The purpose of this work is to evaluate quality of intensity-modulated radiation therapy (IMRT) plans resulting from the use of biological cost functions in comparison to plans designed using a traditional TPS employing dose-volume-based optimization. Treatment planning was performed independently at two institutions. For six cancer patients, including head and neck (one case from each institution), prostate, brain, liver, and rectal cases, segmental multileaf collimator IMRT plans were designed using biological cost functions and compared with clinically used dose-based plans for the same patients. Dose-volume histograms and dosimetric indices, such as minimum, maximum, and mean dose, were extracted and compared between the two types of treatment plans. Comparisons of the generalized equivalent uniform dose (EUD), a previously proposed plan quality index (fEUD), target conformity and heterogeneity indices, and the number of segments and monitor units were also performed. The most prominent feature of the biologically based plans was better sparing of organs at risk (OARs). When all plans from both institutions were combined, the biologically based plans resulted in smaller EUD values for 26 out of 33 OARs by an average of 5.6 Gy (range 0.24 to 15 Gy). Owing to more efficient beam segmentation and leaf sequencing tools implemented in the biologically based TPS compared to the dose-based TPS, an estimated treatment delivery time was shorter in most (five out of six) cases with some plans showing up to 50% reduction. The biologically based plans were generally characterized by a smaller conformity index, but greater heterogeneity index compared to the dose-based plans. Overall, compared to plans based on dose-volume optimization, plans with equivalent target coverage obtained using the biologically based TPS demonstrate improved dose distributions for the majority of normal structures.

Effects of oxygen on intrinsic radiation sensitivity: A test of the relationship between aerobic and hypoxic linear-quadratic (LQ) model parameters.

The poor treatment prognosis for tumors with high levels of hypoxia is usually attributed to the decreased sensitivity of hypoxic cells to ionizing radiation. Mechanistic considerations suggest that linear quadratic (LQ) survival model radiosensitivity parameters for hypoxic (H) and aerobic (A) cells are related by alphaH = alphaA/oxygen enhancement ratio (OER) and (alpha/beta)H=OER(alpha/beta)A. The OER parameter may be interpreted as the ratio of the dose to the hypoxic cells to the dose to the aerobic cells required to produce the same number of DSBs per cell. The validity of these expressions is tested against survival data for mammalian cells irradiated in vitro with low- and high-LET radiation. Estimates of hypoxic and aerobic radiosensitivity parameters are derived from independent and simultaneous least-squares fits to the survival data. An external bootstrap procedure is used to test whether independent fits to the survival data give significantly better predictions than simultaneous fits to the aerobic and hypoxic data. For low-LET radiation, estimates of the OER derived from the in vitro data are between 2.3 and 3.3 for extreme levels of hypoxia. The estimated range for the OER is similar to the oxygen enhancement ratios reported in the literature for the initial yield of DSBs. The half-time for sublethal damage repair was found to be independent of oxygen concentration. Analysis of patient survival data for cervix cancer suggests an average OER less than or equal to 1.5, which corresponds to a pO2 of 5 mm Hg (0.66%) in the in vitro experiments. Because the OER derived from the cervix cancer data is averaged over cells at all oxygen levels, cells irradiated in vivo under extreme levels of hypoxia (<0.5 mm Hg) may have an OER substantially higher than 1.5. The reported analyses of in vitro data, as well as mechanistic considerations, provide strong support for the expressions relating hypoxic and aerobic radiosensitivity parameters. The formulas are also useful for the analysis of clinical data because the number of radiosensitivity parameters that need to be determined is reduced from four to three without a substantial decrease in the ability of the LQ to accurately predict the surviving faction. The relationships among radiosensitivity parameters imply that the dose to the hypoxic subvolume of the tumor needs to be escalated by a factor of the OER to achieve the same level of tumor control as in well oxygenated tumor regions.

Validation of normal tissue complication probability predictions in individual patient: late rectal toxicity.

PURPOSE: To perform validation of risk predictions for late rectal toxicity (LRT) in prostate cancer obtained using a new approach to synthesize published normal tissue complication data. METHODS AND MATERIALS: A published study survey was performed to identify the dose-response relationships for LRT derived from nonoverlapping patient populations. To avoid mixing models based on different symptoms, the emphasis was placed on rectal bleeding. The selected models were used to compute the risk estimates of grade 2+ and grade 3+ LRT for an independent validation cohort composed of 269 prostate cancer patients with known toxicity outcomes. Risk estimates from single studies were combined to produce consolidated risk estimates. An agreement between the actuarial toxicity incidence 3 years after radiation therapy completion and single-study or consolidated risk estimates was evaluated using the concordance correlation coefficient. Goodness of fit for the consolidated risk estimates was assessed using the Hosmer-Lemeshow test. RESULTS: A total of 16 studies of grade 2+ and 5 studies of grade 3+ LRT met the inclusion criteria. The consolidated risk estimates of grade 2+ and 3+ LRT were constructed using 3 studies each. For grade 2+ LRT, the concordance correlation coefficient for the consolidated risk estimates was 0.537 compared with 0.431 for the best-fit single study. For grade 3+ LRT, the concordance correlation coefficient for the consolidated risk estimates was 0.477 compared with 0.448 for the best-fit single study. No evidence was found for a lack of fit for the consolidated risk estimates using the Hosmer-Lemeshow test (P=.531 and P=.397 for grade 2+ and 3+ LRT, respectively). CONCLUSIONS: In a large cohort of prostate cancer patients, selected sets of consolidated risk estimates were found to be more accurate predictors of LRT than risk estimates derived from any single study.

Consolidating risk estimates for radiation-induced complications in individual patient: late rectal toxicity.

PURPOSE: To test the feasibility of a new approach to synthesize published normal tissue complication data using late rectal toxicity in prostate cancer as an example. METHODS AND MATERIALS: A data survey was performed to identify the published reports on the dose-response relationships for late rectal toxicity. The risk estimates for Grade 1 or greater, Grade 2 or greater, and Grade 3 or greater toxicity were obtained for a test cohort of patients treated at our institution. The influence of the potential factors that might have affected the reported toxicity levels was investigated. The studies that did not conform to the general data trends were excluded, and single, combined risk estimates were derived for each patient and toxicity level. RESULTS: A total of 21 studies of nonoverlapping patient populations were identified. Three studies provided dose-response models for more than one level of toxicity. Of these 21 studies, 6, 14, and 5 were used to derive the initial risk estimates for Grade 1, 2, and 3 or greater toxicity, respectively. A comparison of risk estimates between the studies reporting rectal bleeding and rectal toxicity (bleeding plus other symptoms) or between studies with follow-up <36 months and ≥36 months did not reveal significant differences (p ≥ .29 for all comparisons). After excluding three reports that did not conform to the general data trends, the combined risk estimates were derived from 5 reports (647 patients), 11 reports (3,369 patients), and 5 reports (1,330 patients) for Grade 1, 2, and 3 or greater toxicity, respectively. CONCLUSIONS: The proposed approach is feasible and allows for more systematic use of published dose-response data to estimate the complication risks for the individual patient.

Elective lymph node irradiation with intensity-modulated radiotherapy: is conventional dose fractionation necessary?

PURPOSE: Intensity-modulated radiation therapy (IMRT) is the standard of care for head-and-neck cancer (HNC). We treated patients with HNC by delivering either a moderate hypofractionation (MHF) schedule (66 Gy at 2.2 Gy per fraction to the gross tumor [primary and nodal]) with standard dose fractionation (54-60 Gy at 1.8-2.0 Gy per fraction) to the elective neck lymphatics or a conventional dose and fractionation (CDF) schedule (70 Gy at 2.0 Gy per fraction) to the gross tumor (primary and nodal) with reduced dose to the elective neck lymphatics. We analyzed these two cohorts for treatment outcomes. METHODS AND MATERIALS: Between November 2001 and February 2009, 89 patients with primary carcinomas of the oral cavity, larynx, oropharynx, hypopharynx, and nasopharynx received definitive IMRT with or without concurrent chemotherapy. Twenty patients were treated using the MHF schedule, while 69 patients were treated with the CDF schedule. Patient characteristics and dosimetry plans were reviewed. Patterns of failure including local recurrence (LR), regional recurrence (RR), distant metastasis (DM), disease-free survival (DFS), overall survival (OS), and toxicities, including rate of feeding tube placement and percentage of weight loss, were reviewed and analyzed. RESULTS: Median follow-up was 31.2 months. Thirty-five percent of patients in the MHF cohort and 77% of patients in the CDF cohort received chemotherapy. No RR was observed in either cohort. OS, DFS, LR, and DM rates for the entire group at 2 years were 89.3%, 81.4%, 7.1%, and 9.4%, respectively. Subgroup analysis showed no significant differences in OS (p = 0.595), DFS (p = 0.863), LR (p = 0.833), or DM (p = 0.917) between these two cohorts. Similarly, no significant differences were observed in rates of feeding tube placement and percentages of weight loss. CONCLUSIONS: Similar treatment outcomes were observed for MHF and CDF cohorts. A dose of 50 Gy at 1.43 Gy per fraction may be sufficient to electively treat low-risk neck lymphatics.

Renin-Angiotensin system suppression mitigates experimental radiation pneumonitis.

PURPOSE: To find the mitigators of pneumonitis induced by moderate doses of thoracic radiation (10-15 Gy). METHODS AND MATERIALS: Unanesthetized WAG/RijCmcr female rats received a single dose of X-irradiation (10, 12, or 15 Gy at 1.615 Gy/min) to the thorax. Captopril (an angiotensin-converting enzyme inhibitor) or losartan (an angiotensin receptor blocker) was administered in the drinking water after irradiation. Pulmonary structure and function were assessed after 8 weeks in randomly selected rats by evaluating the breathing rate, ex vivo vascular reactivity, and histopathologic findings. Survival analysis was undertaken on all animals, except those scheduled for death. RESULTS: Survival after a dose of 10 Gy to the thorax was not different from that of unirradiated rats for