Reporting and analyzing dose distributions: a concept of equivalent uniform dose

The impact of inter-fraction dose variations on biological equivalent dose (BED): the concept of equivalent constant dose

Inter-fraction dose fluctuations, which appear as a result of setup errors, organ motion and treatment machine output variations, may influence the radiobiological effect of the treatment even when the total delivered physical dose remains constant. The effect of these inter-fraction dose fluctuations on the biological effective dose (BED) has been investigated. Analytical expressions for the BED accounting for the dose fluctuations have been derived. The concept of biological effective constant dose (BECD) has been introduced. The equivalent constant dose (ECD), representing the constant physical dose that provides the same cell survival fraction as the fluctuating dose, has also been introduced. The dose fluctuations with Gaussian as well as exponential probability density functions were investigated. The values of BECD and ECD calculated analytically were compared with those derived from Monte Carlo modelling. The agreement between Monte Carlo modelled and analytical values was excellent (within 1%) for a range of dose standard deviations (0-100% of the dose) and the number of fractions (2 to 37) used in the comparison. The ECDs have also been calculated for conventional radiotherapy fields. The analytical expression for the BECD shows that BECD increases linearly with the variance of the dose. The effect is relatively small, and in the flat regions of the field it results in less than 1% increase of ECD. In the penumbra region of the 6 MV single radiotherapy beam the ECD exceeded the physical dose by up to 35%, when the standard deviation of combined patient setup/organ motion uncertainty was 5 mm. Equivalently, the ECD field was approximately 2 mm wider than the physical dose field. The difference between ECD and the physical dose is greater for normal tissues than for tumours.

The simultaneous integrated boost with proton beams in head and neck patients

The potential benefit of the simultaneous integrated boost (SIB) concept in proton therapy was investigated in a planning study. The proposed fractionation strategy consisted of a SIB treatment following a first phase of conventional fractionation to the elective volume (named SEQ/SIB). The novel method was compared to a conventional sequential fractionation and to a full SIB approach. Treatment plans were designed for five patients on the proton planning system developed for spot scanning at the Paul Scherrer Institute of Villigen (CH). Three to five beams were applied for all plans and fractionations. All effective dose distributions were corrected using biological models to take into account repopulation and time at repopulation onset. Corrected and uncorrected plans were compared on the basis of dosimetry criteria and dose-volume histograms. The results showed a dosimetric advantage for the SEQ/SIB approach in terms of target coverage, without significant disadvantages for risk structures and healthy tissue. Considering the high logistic impact and the limited availability of radiation facilities, the clinical exploitation of accelerated SIB treatment with protons appears promising.

A dynamic supraclavicular field-matching technique for head-and-neck cancer patients treated with IMRT

PURPOSE

The conventional single-isocenter and half-beam (SIHB) technique for matching supraclavicular fields with head-and-neck (HN) intensity-modulated radiotherapy (IMRT) fields is subject to substantial dose inhomogeneities from imperfect accelerator jaw/MLC calibration. It also limits the isocenter location and restricts the useful field size for IMRT. We propose a dynamic field-matching technique to overcome these limitations.

METHODS AND MATERIALS

The proposed dynamic field-matching technique makes use of wedge junctions for the abutment of supraclavicular and HN IMRT fields. The supraclavicular field was shaped with a multileaf collimator (MLC), which was orientated such that the leaves traveled along the superoinferior direction. The leaves that defined the superior field border moved continuously during treatment from 1.5 cm below to 1.5 cm above the conventional match line to generate a 3-cm-wide wedge-shaped junction. The HN IMRT fields were optimized by taking into account the dose contribution from the supraclavicular field to the junction area, which generates a complementary wedge to produce a smooth junction in the abutment region. This technique was evaluated on a polystyrene phantom and 10 HN cancer patients. Treatment plans were generated for the phantom and the 10 patients. Dose profiles across the abutment region were measured in the phantom on films. For patient plans, dose profiles that passed through the center of the neck lymph nodes were calculated using the proposed technique and the SIHB technique, and dose uniformity in the abutment region was compared. Field mismatches of +/- 1 mm and +/- 2 mm because of imperfect jaw/MLC calibration were simulated, and the resulting dose inhomogeneities were studied for the two techniques with film measurements and patient plans. Three-dimensional volumetric doses were analyzed, and equivalent uniform doses (EUD) were computed. The effect of field mismatches on EUD was compared for the two match techniques.

RESULTS

For a perfect jaw/MLC calibration, dose profiles for the 10 patients in the 3-cm match zone had an average inhomogeneity range of -1.6% to +1.6% using the dynamic-matching technique and -3.7% to +3.8% according to the SIHB technique. Measurements showed that dose inhomogeneities that resulted from 1-mm and 2-mm jaw/MLC calibration errors were reduced from as large as 27% and 45% with the SIHB technique to less than 2% and 5.7% with the dynamic technique, respectively. For -1-mm, -2-mm, +1-mm, and +2-mm jaw/MLC calibration errors, respectively, treatment plans for the 10 patients yielded average dose inhomogeneities of -5.9%, -3.0%, +2.7%, and +5.8% with the dynamic technique as compared to -22.8%, -11.1%, +9.8%, and +22.1% with the SIHB technique. Calculation based on a dose-volume histogram (DVH) showed that the SIHB technique resulted in larger changes in EUD of the PTV in the junction area than did the dynamic technique.

CONCLUSION

Compared with the conventional SIHB technique, the dynamic field-matching technique provides superior dose homogeneity in the abutment region between the supraclavicular and HN IMRT fields. The dynamic feathering mechanism substantially reduces dose inhomogeneities that result from imperfect jaw/MLC calibration. In addition, isocenter location in the dynamic field-matching technique can be chosen for reproducible patient setup and for adequate IMRT field size rather than being dictated by the match position. It also allows angling of the supraclavicular field to reduce the volume of healthy lung irradiated, which is impractical with the SIHB technique. In principle, this technique should be applicable to any treatment site that requires the abutment of static and intensity-modulated fields.

Sensitivity of treatment plan optimisation for prostate cancer using the equivalent uniform dose (EUD) with respect to the rectal wall volume parameter

BACKGROUND AND PURPOSE

To analyse the sensitivity of plan optimisation of prostate cancer treatments with respect to changes in the volume parameter (n), when the EUD is used to control the dose in the rectal wall.

PATIENTS AND METHODS

A series of plans was defined, by varying n over a range between 0.08 and 1, and testing different cost functions and beam arrangements. In all cases, the aim was to minimise the EUD in the rectal wall, while ensuring specific dose coverage of the PTV, and limiting the dose in the other OARs. The results were evaluated in terms of 3-D dose distribution and with respect to the current clinical knowledge about late rectal toxicity after irradiation.

RESULTS

Different values of n lead to very similar dose distributions over the PTV (differences in mean dose < 1 Gy, differences in dose given to 99% of the volume < 1%). For the rectal wall, the following observations were made: (a) all cumulative DVH curves crossed each other around 60 Gy; (b) the rectal wall volume receiving doses between 30 and 45 Gy could change by 45 and 30%, respectively, depending on the value of n; (c) for doses higher than 70Gy the differences were typically within 5%. Different values of n also affected the position of isodose surfaces. The distance between the 70 and the 30 Gy isodose curves changed in the AP direction by a factor of 3 when n decreased from 1 to 0.08. High values of n were associated with less dose conformity and a larger volume (at least 20%) of normal tissues receiving 50 Gy or more. All DVHs for the rectal wall were below published dose toxicity thresholds except when the prescribed dose was escalated up to 86 Gy.

CONCLUSIONS

In most cases, the solutions associated with n values up to 0.25 produced similar dose distribution in the rectal wall for doses above 45 Gy, complying with the dose-toxicity thresholds we analysed. The choice of a specific value of n in the optimisation requires an analysis of its effects on the dose distribution for the rectal wall, but also on other aspects, such as the value of the dose to the non-involved normal tissues.

Limitations of a convolution method for modeling geometric uncertainties in radiation therapy: the radiobiological dose-per-fraction effect

The convolution method can be used to model the effect of random geometric uncertainties into planned dose distributions used in radiation treatment planning. This is effectively done by linearly adding infinitesimally small doses, each with a particular geometric offset, over an assumed infinite number of fractions. However, this process inherently ignores the radiobiological dose-per-fraction effect since only the summed physical dose distribution is generated. The resultant potential error on predicted radiobiological outcome [quantified in this work with tumor control probability (TCP), equivalent uniform dose (EUD), normal tissue complication probability (NTCP), and generalized equivalent uniform dose (gEUD)] has yet to be thoroughly quantified. In this work, the results of a Monte Carlo simulation of geometric displacements are compared to those of the convolution method for random geometric uncertainties of 0, 1, 2, 3, 4, and 5 mm (standard deviation). The alpha/betaCTV ratios of 0.8, 1.5, 3, 5, and 10 Gy are used to represent the range of radiation responses for different tumors, whereas a single alpha/betaOAR ratio of 3 Gy is used to represent all the organs at risk (OAR). The analysis is performed on a four-field prostate treatment plan of 18 MV x rays. The fraction numbers are varied from 1-50, with isoeffective adjustments of the corresponding dose-per-fractions to maintain a constant tumor control, using the linear-quadratic cell survival model. The average differences in TCP and EUD of the target, and in NTCP and gEUD of the OAR calculated from the convolution and Monte Carlo methods reduced asymptotically as the total fraction number increased, with the differences reaching negligible levels beyond the treatment fraction number of > or =20. The convolution method generally overestimates the radiobiological indices, as compared to the Monte Carlo method, for the target volume, and underestimates those for the OAR. These effects are interconnected and attributed to assuming an infinite number of fractions inherent in the implementation of the convolution technique, irrespective of the uniqueness of each treatment schedule. Based on the fraction numbers analyzed (1-50), and the range of fraction numbers normally used clinically (> or =20), the convolution method can be used safely to estimate the effects of random geometric uncertainties on prostate treatment radiobiological outcomes, for both the target and the OAR. Although the results of this study is likely to apply to other clinical sites and treatment techniques other than the four-field, further validation similar to those done in this study may be necessary prior to clinical implementation.

A dose-volume-based tool for evaluating and ranking IMRT treatment plans

External beam radiotherapy is commonly used for patients with cancer. While tumor shrinkage and palliation are frequently achieved, local control and cure remain elusive for many cancers. With regard to local control, the fundamental problem is that radiotherapy‐induced normal tissue injury limits the dose that can be delivered to the tumor. While intensity‐modulated radiation therapy (IMRT) allows for the delivery of higher tumor doses and the sparing of proximal critical structures, multiple competing plans can be generated based on dosimetric and/or biological constraints that need to be considered/compared. In this work, an IMRT treatment plan evaluation and ranking tool, based on dosimetric criteria, is presented. The treatment plan with the highest uncomplicated target conformity index (TCI+) is ranked at the top. The TCI+ is a dose‐volume‐based index that considers both a target conformity index (TCI) and a normal tissue‐sparing index (NTSI). TCI+ is designed to assist in the process of judging the merit of a clinical treatment plan. To demonstrate the utility of this tool, several competing lung and prostate IMRT treatment plans are compared. Results show that the plan with the highest TCI+ values accomplished the competing goals of tumor coverage and critical structures sparing best, among rival treatment plans for both treatment sites. The study demonstrates, first, that dose‐volume‐based indices, which summarize complex dose distributions through a single index, can be used to automatically select the optimal plan among competing plans, and second, that this dose‐volume‐based index may be appropriate for ranking IMRT dose distributions.

PACS numbers: 87.53.‐j, 87.53.Tf

Inverse treatment planning with adaptively evolving voxel-dependent penalty scheme

In current inverse planning algorithms it is common to treat all voxels within a target or sensitive structure equally and use structure specific prescriptions and weighting factors as system parameters. In reality, the voxels within a structure are not identical in complying with their dosimetric goals and there exists strong intrastructural competition. Inverse planning objective function should not only balance the competing objectives of different structures but also that of the individual voxels in various structures. In this work we propose to model the intrastructural tradeoff through the modulation of voxel-dependent importance factors and deal with the challenging problem of how to obtain a sensible set of importance factors with a manageable amount of computing. Instead of letting the values of voxel-dependent importance to vary freely during the search process, an adaptive algorithm, in which the importance factors were tied to the local radiation doses through a heuristically constructed relation, was developed. It is shown that the approach is quite general and the EUD-based optimization is a special case of the proposed framework. The new planning tool was applied to study a hypothetical phantom case and a prostate case. Comparison of the results with that obtained using conventional inverse planning technique with structure specific importance factors indicated that the dose distributions from the conventional inverse planning are at best suboptimal and can be significantly improved with the help of the proposed nonuniform penalty scheme.

The Leksell gamma knife Model U versus Model C: a quantitative comparison of radiosurgical treatment parameters

OBJECTIVE

We present a quantitative comparison of radiosurgery treatments for cavernous sinus tumors using the Leksell gamma knife Model U versus the Model C with automatic positioning system (APS) (Elekta Instruments, Norcross, GA).

METHODS

At our medical center from August 1994 through May 2000, the Model U was used to treat 96 patients (37 men [39%] and 59 women [61%]; median age, 54.5 yr) with benign cavernous sinus tumors: 43 meningiomas (45%), 48 pituitary tumors (50%), and 5 others (5%). From June 2000 through April 2002, the Model C with APS treated 45 patients (20 men [44%] and 25 women [56%]; median age, 51.4 yr) with 15 meningiomas (33%), 29 pituitary tumors (65%), and 1 schwannoma (2%). The two groups had similar treated tumor volumes (Model U mean, 4.3 cm(3); Model C mean, 4.2 cm(3)), equivalent tumor distances from critical structures (optic nerve, chiasm, and pons), comparable distributions in Sekhar tumor grades, and the same median prescribed dose of 15 Gy to the 50% isodose line at the tumor periphery. All planning and treatments were performed by the same radiosurgery team to minimize dosage to adjacent critical tissues and to optimize conformity index.

RESULTS

Analysis of multiple treatment parameters showed that the Model C plans were superior. Model C treatments had an improved conformity index (Model U mean, 1.7; Model C mean, 1.6; P < 0.02) and a lower underdosed tumor volume (Model U mean, 0.4 cm(3); Model C mean, 0.1 cm(3); P < 0.004). The total treated volume and the excess treated volume were similar. The Model C group had a reduction in optic chiasm dose (Model C mean dose, 3.8 Gy; Model U mean dose, 5.3 Gy; P < 0.0001). The average number of isocenters was slightly higher for the Model C group (6.7 versus 6 for the Model U), but with a lower mean number of collimator sizes (1 versus 2 for the Model U). Model C plans required a mean of 93 fewer plugs per treatment, thus contributing to an estimated 67.6 minutes saved per treatment session.

CONCLUSION

Comparison of radiosurgery treatments using the Leksell gamma knife Model U versus the Model C with APS was performed by quantitative analysis of treatment parameters on a cohort of benign cavernous sinus tumors. Treatment plans using the Model C resulted in better tumor coverage (improved conformity, less underdosed tumor volume) and decreased optic chiasm dose. An estimated average of 1 hour was saved per treatment when using the Model C with APS.

Incorporating model parameter uncertainty into inverse treatment planning

Radiobiological treatment planning depends not only on the accuracy of the models describing the dose-response relation of different tumors and normal tissues but also on the accuracy of tissue specific radiobiological parameters in these models. Whereas the general formalism remains the same, different sets of model parameters lead to different solutions and thus critically determine the final plan. Here we describe an inverse planning formalism with inclusion of model parameter uncertainties. This is made possible by using a statistical analysis-based frameset developed by our group. In this formalism, the uncertainties of model parameters, such as the parameter a that describes tissue-specific effect in the equivalent uniform dose (EUD) model, are expressed by probability density function and are included in the dose optimization process. We found that the final solution strongly depends on distribution functions of the model parameters. Considering that currently available models for computing biological effects of radiation are simplistic, and the clinical data used to derive the models are sparse and of questionable quality, the proposed technique provides us with an effective tool to minimize the effect caused by the uncertainties in a statistical sense. With the incorporation of the uncertainties, the technique has potential for us to maximally utilize the available radiobiology knowledge for better IMRT treatment.

Comparison between manual and automatic segment generation in step-and-shoot IMRT of prostate cancer

PURPOSE

To compare two methods to generate treatment plans for intensity-modulated radiotherapy (IMRT) of prostate cancer, delivered in a step-and-shoot mode. The first method uses fluence optimization (inverse planning) followed by conversion of the fluence weight map into a limited number of segments. In the second method, segments are manually assigned using a class solution (forward planning), followed by computer optimization of the segment weights.

METHODS

Treatment plans for IMRT, utilizing a simultaneous integrated boost, were created. Plans comprise a five-field technique to deliver 78 Gy to the prostate plus seminal vesicles. Five patients were evaluated. Optimization objectives of both planning approaches concerned dose coverage of the target volumes and the dose distribution in the rectal wall. The two methods were evaluated by comparing dose distributions, the complexity of the resulting plan and the time expenditure to generate and to deliver the plan.

RESULTS

For both planning approaches 99% of the target volumes received 95% of the prescribed dose, which complies with our planning objectives. Inverse planning resulted in more conformal dose distributions than forward planning (conformity index: 1.37 versus 1.51). Inverse planning reduced the dose to the rectal wall compared to a manually designed plan, albeit to a small extent. The theoretical probability of severe rectal proctitis and/or stenosis was reduced on average by 1.9% with inverse planning. Maximal sparing of the rectal wall was achieved with inverse planning for a patient whose target volume was partly wrapped around the rectum. The number of segments generated with inverse planning ranged between 33 and 52, and between 9 and 13 segments for manually created segments.

CONCLUSION

Dose coverage of the planning target volumes is adequate for both approaches of planning. Inverse planning results in slightly better dose distributions with respect to the rectal wall compared to manual planning, at the cost of an increase of the number of segments by a factor of 3.

Simultaneous integrated boost for breast cancer using imrt: a radiobiological and treatment planning study

PURPOSE

The purpose of this work is to explore the possibility of using intensity-modulated radiation therapy (IMRT) to deliver the boost dose to the tumor bed simultaneously with the whole-breast IMRT to reduce the radiation treatment time by 1-2 weeks.

METHODS AND MATERIALS

The biologically effective dose (BED) for different treatments was calculated using the linear-quadratic (LQ) model with parameters previously derived for breast cancer from clinical data (alpha/beta = 10Gy, alpha = 0.3Gy(-1)). A potential doubling time of 15 days (from in vitro measurements) for breast cancer and a generic alpha/beta ratio of 3 Gy for normal tissues were used. A series of regimens that use IMRT as initial treatment and an IMRT simultaneous integrated boost (SIB) were derived using biologic equivalence to conventional schedules. Possible treatment plans with IMRT SIB to the tumor bed were generated for 2 selected breast patients, 1 with a shallow tumor and 1 with a deep-seated tumor. Plans with a simultaneous integrated electron boost were also generated for comparison. Dosimetric merits of these plans were evaluated based on dose volume histograms.

RESULTS

A commonly used conventional treatment of 45 Gy (1.8 Gy x 25) to the whole breast and then a boost of 20 Gy (2 Gy x 10) is biologically equivalent to an alternative plan of 1.8 Gy x 25 to the whole breast with a 2.4 Gy x 25 SIB to the tumor bed. The new regime reduces treatment time from 7 to 5 weeks. For the patient with a deep-seated tumor, the IMRT plans reduce the volume of the breast that receives high doses (compared with the conventional photon boost plan) and provides good coverage of the target volumes.

CONCLUSION

It is biologically and dosimetrically feasible to reduce the overall treatment time for breast radiotherapy by using an IMRT simultaneous integrated boost. For selected patient groups, IMRT plans with a new regimen can be equal to or better than conventional plans.

Les contraintes aux organes à risque en radiothérapie par modulation d'intensité des cancers ORL

Constraint definitions in intensity modulated radiation therapy is a key point factor during the treatment planning process. In literature some data are available about dose constraints and volumes according to the tissue architectures. Following ICRU recommendations, organs at risk organized in a parallel structure could receive an acceptably small proportion of high dose component. Mean dose and dose volume histogram is a most convenient tool for incorporating such constraints. Organs described as a serial structure are supposed to receive less than the given maximum dose, directly linked to the occurrence of complications. Dmax is the best way to describe such events. These constraints are new tools in radiation therapy, available for optimizing the dose distribution in target volume, sparing the organs at risk to protect the organ function or at least decreasing the late functional damages like xerostomia. It is necessary to define with accuracy gross target volumes and clinical target volume with available radio-anatomical guidelines before introducing current constraints on each volume in the inverse dosimetry. The management of these constraints remains under the responsibility of the clinicians. A permanent compromise has to be chosen between homogeneity of the dose distribution in the target volume and the probability of preserving functions of organs at risk.

Prostate implant evaluation using tumour control probability—the effect of input parameters

In this paper, we examine the effect of treatment parameters in a model used to evaluate permanent prostate implants. The model considers the prostate to be composed of 12 sub-sections, each sub-section is assigned a cell density based on the probability of finding cancer foci in that sub-section. Wasted dose as a result of the dose rate from the implant falling below a level adequate to counteract repopulation was found to vary by 2-16% over the range of radiosensitivity and repopulation rates considered. Within the model, applied to five dose distributions, the uncertainty in the tumour control probability (TCP) values calculated for each sub-section as a result of differences in the model parameters, was found to be less than 12% in most cases for the good quality implants. The difference in TCP values was much larger for the poor quality implant. Substituting a heterogeneous distribution of alpha for a single mean value resulted in generally lower TCP values though introducing a cutoff value with a Gaussian distribution had a profound effect on the calculated values. Despite uncertainties in the parameters, the model was able to identify sub-sections at risk of local recurrence but as a result of these uncertainties, the TCP values can only be considered in the relative rather than absolute sense.

Biological effective dose for comparison and combination of external beam and low-dose rate interstitial brachytherapy prostate cancer treatment plans

We report a methodology for comparing and combining dose information from external beam radiotherapy (EBRT) and interstitial brachytherapy (IB) components of prostate cancer treatment using the biological effective dose (BED). On a prototype early-stage prostate cancer patient treated with EBRT and low-dose rate I-125 brachytherapy, a 3-dimensional dose distribution was calculated for each of the EBRT and IB portions of treatment. For each component of treatment, the BED was calculated on a point-by-point basis to produce a BED distribution. These individual BED distributions could then be summed for combined therapies. BED dose-volume histograms (DVHs) of the prostate, urethra, rectum, and bladder were produced and compared for various combinations of EBRT and IB. Transformation to BED enabled computation of the relative contribution of each modality to the prostate dose, as the relative weighting of EBRT and IB was varied. The BED-DVHs of the prostate and urethra demonstrated dramatically increased inhomogeneity with the introduction of even a small component of IB. However, increasing the IB portion relative to the EBRT component resulted in lower dose to the surrounding normal structures, as evidenced by the BED-DVHs of the bladder and rectum. Conformal EBRT and low-dose rate IB conventional dose distributions were successfully transformed to the common "language" of BED distributions for comparison and for merging prostate cancer radiation treatment plans. The results of this analysis can assist physicians in quantitatively determining the best combination and weighting of radiation treatment modalities for individual patients.

Penalized likelihood fluence optimization with evolutionary components for intensity modulated radiation therapy treatment planning

A novel iterative penalized likelihood algorithm with evolutionary components for the optimization of beamlet fluences for intensity modulated radiation therapy (IMRT) is presented. This algorithm is designed to be flexible in terms of the objective function and automatically escalates dose, as long as the objective function increases and all constraints are met. For this study, the objective function employed was the product of target equivalent uniform dose (EUD) and fraction of target tissue within set homogeneity constraints. The likelihood component of the algorithm iteratively attempts to minimize the mean squared error between a homogeneous dose prescription and the actual target dose distribution. The updated beamlet fluences are then adjusted via a quadratic penalty function that is based on the dose-volume histogram (DVH) constraints of the organs at risk. The evolutionary components were included to prevent the algorithm from converging to a local maximum. The algorithm was applied to a prostate cancer dataset, with especially difficult DVH constraints on bladder, rectum, and femoral heads. Dose distributions were generated for manually selected sets of three-, four-, five-, and seven-field treatment plans. Additionally, a global search was performed to find the optimal orientations for an axial three-beam plan. The results from this optimal orientation set were compared to results for manually selected orientation (gantry angle) sets of 3- (0 degrees, 90 degrees, 270 degrees), 4- (0 degrees, 90 degrees, 180 degrees, 270 degrees), 5- (0 degrees, 50 degrees, 130 degrees, 230 degrees, 310 degrees), and 7- (0 degrees, 40 degrees, 90 degrees, 140 degrees, 230 degrees, 270 degrees, 320 degrees) field axial treatment plans. For all the plans generated, all DVH constraints were met and average optimization computation time was approximately 30 seconds. For the manually selected orientations, the algorithm was successful in providing a relatively homogeneous target dose distribution, while simultaneously satisfying dose-volume limits by diverting dose away from proximal critical structures. The global search for an optimal three-beam orientation set yielded gantry angles of 70 degrees, 170 degrees, and 320 degrees. The EUD for this orientation set was 58 Gy, with 96% of the target within the set upper and lower limits. In comparison, optimized EUDs for the manually selected orientation sets of three, four, five and seven beams were 52.3, 52.6, 56.9, and 61.3 Gy, respectively. The orientation optimized three-beam plan yielded higher EUDs than the manually selected three-, four-, and five-beam plans, but lower EUDs than the seven-beam plan. In conclusion, a novel penalized likelihood algorithm with evolutionary components has successfully been implemented to optimize beamlet fluences for IMRT. Initial results are promising for dose conformity and uniformity of dose to target. When combined with optimal beam orientation selection for prostate cancer treatment planning, the results indicate that plans with a small number of optimized beam orientations achieve results comparable to those with a larger number of conventionally oriented beams.

Spinal cord tolerance to high-dose fractionated 3D conformal proton-photon irradiation as evaluated by equivalent uniform dose and dose volume histogram analysis

PURPOSE

To evaluate cervical spinal cord tolerance using equivalent uniform dose (EUD) and dose volume histogram (DVH) analysis after proton-photon radiotherapy.

METHODS AND MATERIAL

The 3D dose distributions were analyzed in 85 patients with cervical vertebral tumors. Mean follow-up was 41.3 months. The mean prescribed dose was 76.3 Cobalt Gray Equivalent (CGE = proton dose x RBE 1.1). Dose constraints to the center and the surface of the cervical cord were 55-58 CGE and 67-70 CGE, respectively. Dose parameters, DVH and EUD, were calculated for each patient. The spinal cord toxicity was graded using the European Organization for Research and Treatment of Cancer (EORTC) and Radiation Therapy Oncology Group (RTOG) late effects scoring system.

RESULTS

Thirteen patients experienced Grade 1-2 toxicity. Four patients had Grade 3 toxicity. For the dose range used in this study, none of the dosimetric parameters was found to be associated with the observed distribution of cord toxicities. The only factor significantly associated with cord toxicity was the number of surgeries before irradiation.

CONCLUSION

The data and our analysis suggest that the integrity of the normal musculoskeletal supportive tissues and vascular supply may be important confounding factors of toxicity at these dose levels. The results also indicate that the cervical spinal cord dose constraints used in treating these patients are appropriate for conformal proton-photon radiotherapy.

A preliminary study of the role of modulated electron beams in intensity modulated radiotherapy, using automated beam orientation and modality selection

PURPOSE

To develop an algorithm for optimal beam arrangement selection in intensity-modulated radiotherapy (IMRT) of mixed photon and electron beams. To apply this algorithm to study the utility of modulated electron beams in the context of IMRT planning.

METHODS AND MATERIALS

The optimization algorithm selects, for a user-specified number of beams, the optimal IMRT arrangement (beam orientations, and photon/electron modality for each orientation) using a novel fast heuristic intensity modulation procedure. The algorithm was employed to select optimal beam arrangements for breast (two, four, and six axial beams) and head-and-neck (three, four, five, and seven nonaxial beams) cases.

RESULTS

For the two cases, increasing the number of selected beams: (1) increased the number of electron beams for the breast case, but not more than one electron beam was selected for the head-and-neck case; (2) decreased critical structure doses for both cases; and (3) decreased target homogeneity for the breast case, but improved it for the head-and-neck case.

CONCLUSIONS

In the two cases analyzed using the selection algorithm, the primary role of modulated electrons differs based on treatment site-normal tissue dose reduction in breast and target homogeneity improvement in head and neck. Although this preliminary study with two cases appears to suggest that the role of intensity-modulated electrons differs based on treatment site, further investigation of large numbers of cases and varied treatment sites are required to establish a definitive conclusion.

Dose–volume response analyses of late rectal bleeding after radiotherapy for prostate cancer

PURPOSE

To compare the fits of various normal tissue complication probability (NTCP) models to a common set of late rectal toxicity data, with the aim of identifying the best model for predicting late rectal injury after irradiation.

METHODS AND MATERIALS

Late toxicity data from 128 prostate cancer patients treated on protocol with three-dimensional conformal radiotherapy at The University of Texas M.D. Anderson Cancer Center (UTMDACC) were analyzed. The dose-volume histogram for total rectal volume, including contents, was obtained for each patient, and the presence or absence of Grade 2 or worse rectal bleeding within 2 years of treatment was scored. Five different NTCP models were fitted to the data using maximum likelihood analysis: the Lyman model, the mean dose model, a parallel architecture model, and models based on either a cutoff dose or a cutoff volume.

RESULTS

All five of the NTCP models considered provided very similar fits to the UTMDACC rectal bleeding data. In particular, none of the more highly parameterized models (the four-parameter parallel model, three-parameter Lyman model, or three-parameter cutoff dose and volume models) provided a better fit than the simplest of the models, the two-parameter NTCP model describing rectal bleeding as a probit function of mean dose to rectum.

CONCLUSION

No dose-volume response model has yet been identified that provides a better description of the UTMDACC rectal toxicity data than the mean dose model. Because this model has relatively low predictive accuracy, the need to identify a better model remains.

Feasibility of optimizing the dose distribution in lung tumors using fluorine-18-fluorodeoxyglucose positron emission tomography and single photon emission computed tomography guided dose prescriptions

The information provided by functional images may be used to guide radiotherapy planning by identifying regions that require higher radiation dose. In this work we investigate the dosimetric feasibility of delivering dose to lung tumors in proportion to the fluorine-18-fluorodeoxyglucose activity distribution from positron emission tomography (FDG-PET). The rationale for delivering dose in proportion to the tumor FDG-PET activity distribution is based on studies showing that FDG uptake is correlated to tumor cell proliferation rate, which is shown to imply that this dose delivery strategy is theoretically capable of providing the same duration of local control at all voxels in tumor. Target dose delivery was constrained by single photon emission computed tomography (SPECT) maps of normal lung perfusion, which restricted irradiation of highly perfused lung and imposed dose-function constraints. Dose-volume constraints were imposed on all other critical structures. All dose-volume/function constraints were considered to be soft, i.e., critical structure doses corresponding to volume/function constraint levels were minimized while satisfying the target prescription, thus permitting critical structure doses to minimally exceed dose constraint levels. An intensity modulation optimization methodology was developed to deliver this radiation, and applied to two lung cancer patients. Dosimetric feasibility was assessed by comparing spatially normalized dose-volume histograms from the nonuniform dose prescription (FDG-PET proportional) to those from a uniform dose prescription with equivalent tumor integral dose. In both patients, the optimization was capable of delivering the nonuniform target prescription with the same ease as the uniform target prescription, despite SPECT restrictions that effectively diverted dose from high to low perfused normal lung. In one patient, both prescriptions incurred similar critical structure dosages, below dose-volume/function limits. However, in the other patient, critical structure dosage from the nonuniform dose prescription exceeded dose-volume/function limits, and greatly exceeded that from the uniform dose prescription. Strict compliance to dose-volume/ function limits would entail reducing dose proportionality to the FDG-PET activity distribution, thereby theoretically reducing the duration of local control. Thus, even though it appears feasible to tailor lung tumor dose to the FDG-PET activity distribution, despite SPECT restrictions, strict adherence to dose-volume/function limits could compromise the effectiveness of functional image guided radiotherapy.

Dose escalation in permanent brachytherapy for prostate cancer: dosimetric and biological considerations

No prospective dose escalation study for prostate brachytherapy (PB) with permanent implants has been reported. In this work, we have performed a dosimetric and biological analysis to explore the implications of dose escalation in PB using 125I and 103Pd implants. The concept of equivalent uniform dose (EUD), proposed originally for external-beam radiotherapy (EBRT), is applied to low dose rate brachytherapy. For a given 125I or 103Pd PB, the EUD for tumour that corresponds to a dose distribution delivered by EBRT is calculated based on the linear quadratic model. The EUD calculation is based on the dose volume histogram (DVH) obtained retrospectively from representative actual patient data. Tumour control probabilities (TCPs) are also determined in order to compare the relative effectiveness of different dose levels. The EUD for normal tissue is computed using the Lyman model. A commercial inverse treatment planning algorithm is used to investigate the feasibility of escalating the dose to prostate with acceptable dose increases in the rectum and urethra. The dosimetric calculation is performed for five representative patients with different prostate sizes. A series of PB dose levels are considered for each patient using 125I and 103Pd seeds. It is found that the PB prescribed doses (minimum peripheral dose) that give an equivalent EBRT dose of 64.8, 70.2, 75.6 and 81 Gy with a fraction size of 1.8 Gy are 129, 139, 150 and 161 Gy for 125I and 103, 112, 122 and 132 Gy for 103Pd implants, respectively. Estimates of the EUD and TCP for a series of possible prescribed dose levels (e.g., 145, 160, 170 and 180 Gy for 125I and 125, 135, 145 and 155 for 103Pd implants) are tabulated. The EUD calculation was found to depend strongly on DVHs and radiobiological parameters. The dosimetric calculations suggest that the dose to prostate can be escalated without a substantial increase in both rectal and urethral dose. For example, increasing the PB prescribed dose from 145 to 180 Gy increases EUD for the rectum by only 3%. Our studies indicate that the dose to urethra can be kept within 100-120% of the prescription dose for all the dose levels studied. In conclusion, dose escalation in permanent implant for localized prostate cancer may be advantageous. It is dosimetrically possible to increase dose to prostate without a substantial increase in the dose to the rectum and urethra. Based on the results of our studies, a prospective dose escalation trial for prostate permanent implants has been initiated at our institution.

A unifying framework for multi-criteria fluence map optimization models

Models for finding treatment plans for intensity modulated radiation therapy are usually based on a number of structure-based treatment plan evaluation criteria, which are often conflicting. Rather than formulating a model that a priori quantifies the trade-offs between these criteria, we consider a multi-criteria optimization approach that aims at finding the so-called undominated treatment plans. We present a unifying framework for studying multi-criteria optimization problems for treatment planning that establishes conditions under which treatment plan evaluation criteria can be transformed into convex criteria while preserving the set of undominated treatment plans. Such transformations are identified for many of the criteria that have been proposed to date, establishing equivalences between these criteria. In addition, it is shown that the use of a nonconvex criterion can often be avoided by transformation to an equivalent convex criterion. In particular, we show that models employing criteria such as tumour control probability, normal tissue complication probability, probability of uncomplicated tumour control, as well as sigmoidal transformations of (generalized) equivalent uniform dose are equivalent to models formulated in terms of separable voxel-based criteria that penalize dose in individual voxels.

Comparative dosimetric evaluation of the simultaneous integrated boost with photon intensity modulation in head and neck cancer patients

BACKGROUND AND PURPOSE

The objective of this study is to evaluate, at planning and dosimetric level, the potential benefits of the simultaneous integrated boost (SIB) concept with intensity-modulated radiation therapy (IMRT), using a comparative analysis on physical dose distributions corrected for radiobiological models. The concept of SIB at the end of the treatment has been analysed as an alternative acceleration scheme.

PATIENTS AND METHODS

Physical dose distributions were computed on a commercial planning system (Varian Cadplan-Helios) for five patients presenting with advanced head and neck carcinomas. Treatment plans were designed using five IMRT beams. Three fractionation strategies were compared in the study: the standard sequential irradiation SEQ of elective and boost volumes, the pure SIB, and a modified SIB (SEQ/SIB), where the actual SIB follows a first phase of conventional fractionation to the elective volume. All physical dose distributions were corrected using a linear quadratic biological model, taking into account also repopulation and time at repopulation onset. Objective quantities, derived from biological dose volume histograms, were used for the analysis.

RESULTS

Physical doses equivalent to 50 and 80 Gy (in fractions of 2 Gy) to elective volume and boost were calculated for the SIB and SEQ/SIB regimes. With SIB 54 and 72 Gy dose levels have to be delivered in 30 fractions, while in the SEQ/SIB scheme 36 Gy are delivered in 20 sessions to the elective volume, and further 18 and 35.5 Gy during the last 10 fractions are delivered to elective volume and boost, respectively (for a total physical dose of 71.5 Gy). The comparison showed: (1) the boost target homogeneity resulted in generally acceptable and comparable among sequential and modified SIB schemes, while it was statistically worse for the pure SIB approach; (2) the fraction of elective target volume not included in the boost volume was characterised by a higher level of dose heterogeneity; (3) the spinal cord never reached tolerance levels and maximum point dose was on average below 38 Gy (biologically corrected to 2 Gy/fraction); and (4) sparing of parotid glands strongly depends on their eventual inclusion in the target volumes: for glands not included or only partially included, it was possible on average to keep the dose to 2/3 of the volume below 29 Gy for all regimes (32 Gy as physical dose).

CONCLUSIONS

Feasibility of SIB techniques and in particular of the modified SIB appears to be dosimetrically proven and the results reported here justify the activation of a phase I protocol to verify clinically their impact using IMRT photon-based techniques.

Incorporation of functional imaging data in the evaluation of dose distributions using the generalized concept of equivalent uniform dose

Advances in the fields of IMRT and functional imaging have greatly increased the prospect of escalating the dose to highly active or hypoxic tumour sub-volumes and steering the dose away from highly functional critical structure regions. However, current clinical treatment planning and evaluation tools assume homogeneous activity/function status in the tumour/critical structures. A method was developed to incorporate tumour/critical structure heterogeneous functionality in the generalized concept of equivalent uniform dose (EUD). The tumour and critical structures functional EUD (FEUD) values were calculated from the dose-function histogram (DFH), which relates dose to the fraction of total function value at that dose. The DFH incorporates flouro-deoxyglucose positron emission tomography (FDG-PET) functional data for tumour, which describes the distribution of metabolically active tumour clonogens, and single photon emission computed tomography (SPECT) perfusion data for critical structures. To demonstrate the utility of the method, the lung dose distributions of two non-small cell lung cancer patients, who received 3D conformal external beam radiotherapy treatment with curative intent, were evaluated. Differences between the calculated lungs EUD and FEUD values of up to 50% were observed in the 3D conformal plans. In addition, a non-small cell lung cancer patient was inversely planned with a target dose prescription of 76 Gy. Two IMRT plans (plan-A and plan-B) were generated for the patient based on the CT, FDG-PET and SPECT treatment planning images using dose-volume objective functions. The IMRT plans were generated with the goal of achieving more critical structures sparing in plan-B than plan-A. Results show the target volume EUD in plan-B is lower than plan-A by 5% with a value of 73.31 Gy, and the FEUD in plan-B is lower than plan-A by 2.6% with a value of 75.77 Gy. The FEUD plan-B values for heart and lungs were lower than plan-A by 22% and 18%, respectively. While EUD values show plan-A is marginally better than plan-B in terms of target volumetric coverage, the FEUD plan-B values show adequate target function coverage with significant critical structure function sparing. In conclusion, incorporating functional data in the calculation of EUD is important in evaluating the biological merit of treatment plans.

Application of influence diagrams to prostate intensity-modulated radiation therapy plan selection

The purpose is to incorporate clinically relevant factors such as patient-specific and dosimetric information as well as data from clinical trials in the decision-making process for the selection of prostate intensity-modulated radiation therapy (IMRT) plans. The approach is to incorporate the decision theoretic concept of an influence diagram into the solution of the multiobjective optimization inverse planning problem. A set of candidate IMRT plans was obtained by varying the importance factors for the planning target volume (PTV) and the organ-at-risk (OAR) in combination with simulated annealing to explore a large part of the solution space. The Pareto set for the PTV and OAR was analysed to demonstrate how the selection of the weighting factors influenced which part of the solution space was explored. An influence diagram based on a Bayesian network with 18 nodes was designed to model the decision process for plan selection. The model possessed nodes for clinical laboratory results, tumour grading, staging information, patient-specific information, dosimetric information, complications and survival statistics from clinical studies. A utility node was utilized for the decision-making process. The influence diagram successfully ranked the plans based on the available information. Sensitivity analyses were used to judge the reasonableness of the diagram and the results. In conclusion, influence diagrams lend themselves well to modelling the decision processes for IMRT plan selection. They provide an excellent means to incorporate the probabilistic nature of data and beliefs into one model. They also provide a means for introducing evidence-based medicine, in the form of results of clinical trials, into the decision-making process.

Therapeutic advantage of grid irradiation for large single fractions

PURPOSE

In the present work, we used model calculations of cell survival to compare the effects of single fraction high-dose grid therapy with those of uniform dose delivery on tumor and normal tissues.

METHODS AND MATERIALS

The grid consisted of a hexagonal pattern of divergent holes in a Cerrobend block. A linear-quadratic model was used to find the surviving fraction of tumor and normal tissue cells after high-dose irradiation. Equivalent uniform doses were determined according to the tumor cell kill. The ratio of the normal tissue surviving fraction under grid irradiation to that obtained under equivalent uniform dose irradiation was taken as a measure of therapeutic gain.

RESULTS

The therapeutic ratio varied from 0.80 to 13.22 for the range of cell sensitivities investigated, with single fraction doses of 10.0-20.0 Gy. Optimization studies showed no significant dependence of therapeutic gain on hole spacing.

CONCLUSION

With high, single-fraction doses, grid irradiation revealed a therapeutic advantage over uniform dose irradiation whenever the tumor and surrounding normal tissues cells were equally radiosensitive, or, particularly, if the tumor cells were more radioresistant than the normal cells. The therapeutic gain did not appear to be a strong function of grid design.

Three-dimensional conformal vs. intensity-modulated radiotherapy in head-and-neck cancer patients: comparative analysis of dosimetric and technical parameters

BACKGROUND AND PURPOSE

The use of intensity-modulated radiotherapy (IMRT) is now widely advocated for the treatment of head-and-neck cancers, to increase the therapeutic ratio of radiotherapy used as sole modality of treatment or in combination with chemotherapy. This report aims to summarize the technical and dosimetric factors to be taken into consideration to assess the respective advantages of the various high conformality treatments in radiotherapy, especially in the framework of quality assurance procedures.

MATERIALS AND METHODS

Twenty-six head-and-neck cancer patients were irradiated following a feasibility internal protocol with IMRT. Treatments were performed with either the static step-and-shoot (20) or the dynamic sliding window (6) techniques on a 6 MV Varian Clinac equipped with a multileaf collimator with 80 leaves. Dose plans were computed using commercial treatment planning systems: MDS-Nordion Helax-TMS for static cases and Varian Eclipse for dynamic cases. Dose plans were evaluated in terms of physical quantities based on dose-volume histograms and isodose distributions. Each IMRT plan was also compared to a reference 3D conformal therapy plan (3DCRT).

RESULTS

Elective target volumes ranged from 530 to 1151 cm(3) with a mean of 780 +/- 141 cm(3). Boost volumes ranged from 248 to 832 cm(3) with a mean of 537 +/- 165 cm(3). Thirty-two dose plans were generated with static technique and 10 with dynamic. In the static mode, 6.8 +/- 3.4 fields were applied on average with 12.5 +/- 1.3 segments per field. In the static mode, 264 +/- 56 MU per Gy were erogated, whereas in the dynamic mode, 387 +/- 126 MU per Gy were erogated, to be compared to 147 +/- 20 computed for reference 3DCRT plans. For all target volumes in general, conformity was improved compared to 3DCRT (e.g. V(95) increased from 85% to 93% with p < 0.001, or equivalent uniform dose normalized to prescribed dose increased from 0.86 to 0.96 with p = 0.002). Irradiation of parotid glands or spinal cord improved, as well: For parotids, D(2/3V) reduced from 59 Gy to 41 Gy (p < 0.001). For spinal cord, D(max) reduced from about 40 Gy to about 30 Gy (p < 0.001).

A novel linear programming approach to fluence map optimization for intensity modulated radiation therapy treatment planning

We present a novel linear programming (LP) based approach for efficiently solving the intensity modulated radiation therapy (IMRT) fluence-map optimization (FMO) problem to global optimality. Our model overcomes the apparent limitations of a linear-programming approach by approximating any convex objective function by a piecewise linear convex function. This approach allows us to retain the flexibility offered by general convex objective functions, while allowing us to formulate the FMO problem as a LP problem. In addition, a novel type of partial-volume constraint that bounds the tail averages of the differential dose-volume histograms of structures is imposed while retaining linearity as an alternative approach to improve dose homogeneity in the target volumes, and to attempt to spare as many critical structures as possible. The goal of this work is to develop a very rapid global optimization approach that finds high quality dose distributions. Implementation of this model has demonstrated excellent results. We found globally optimal solutions for eight 7-beam head-and-neck cases in less than 3 min of computational time on a single processor personal computer without the use of partial-volume constraints. Adding such constraints increased the running times by a factor of 2-3, but improved the sparing of critical structures. All cases demonstrated excellent target coverage (> 95%), target homogeneity (< 10% overdosing and < 7% underdosing) and organ sparing using at least one of the two models.

Analysis of a large number of clinical studies for breast cancer radiotherapy: estimation of radiobiological parameters for treatment planning

Numerous studies of early-stage breast cancer treated with breast conserving surgery (BCS) and radiotherapy (RT) have been published in recent years. Both external beam radiotherapy (EBRT) and/or brachytherapy (BT) with different fractionation schemes are currently used. The present RT practice is largely based on empirical experience and it lacks a reliable modelling tool to compare different RT modalities or to design new treatment strategies. The purpose of this work is to derive a plausible set of radiobiological parameters that can be used for RT treatment planning. The derivation is based on existing clinical data and is consistent with the analysis of a large number of published clinical studies on early-stage breast cancer. A large number of published clinical studies on the treatment of early breast cancer with BCS plus RT (including whole breast EBRT with or without a boost to the tumour bed, whole breast EBRT alone, brachytherapy alone) and RT alone are compiled and analysed. The linear quadratic (LQ) model is used in the analysis. Three of these clinical studies are selected to derive a plausible set of LQ parameters. The potential doubling time is set a priori in the derivation according to in vitro measurements from the literature. The impact of considering lower or higher T(pot) is investigated. The effects of inhomogeneous dose distributions are considered using clinically representative dose volume histograms. The derived LQ parameters are used to compare a large number of clinical studies using different regimes (e.g., RT modality and/or different fractionation schemes with different prescribed dose) in order to validate their applicability. The values of the equivalent uniform dose (EUD) and biologically effective dose (BED) are used as a common metric to compare the biological effectiveness of each treatment regime. We have obtained a plausible set of radiobiological parameters for breast cancer: alpha = 0.3 Gy(-1), alpha/beta = 10 Gy and sub-lethal damage repair time T(rep) = 1 h (mono-exponential behaviour is assumed). This set of parameters is consistent with in vitro experiments and with previously reported analyses. Using this set of parameters, we have found that most of the studies, using BCS plus whole breast RT and a boost to the tumour bed, have EUDs ranging from 60-70 Gy. No correlation is found between BED and the local recurrence rate. The treatments of BCS plus brachytherapy alone have a wide range of EUD (30-50 Gy), which is significantly lower than the treatments with whole breast EBRT plus a boost of the tumour bed. The studies with different fractionation schemes for whole breast EBRT also show a significant variation of EUD. Carefully designed clinical studies with large numbers of patients are required to determine clinically the relative effectiveness of these treatment variations. The derived LQ parameter set based on clinical data is consistent with in vitro experiments and previous studies. As demonstrated in the present work, these radiobiological parameters can be potentially useful in radiotherapy treatment planning for early breast cancer, e.g., in comparing biological effectiveness of different radiotherapy modalities, different fractionation schemes and in designing new treatment strategies.

Impact of prolonged fraction delivery times on tumor control: a note of caution for intensity-modulated radiation therapy (IMRT)

PURPOSE

Intensity-modulated radiation therapy (IMRT) allows greater dose conformity to the tumor target. However, IMRT, especially static delivery, usually requires more time to deliver a dose fraction than conventional external beam radiotherapy (EBRT). The purpose of this work is to explore the potential impact of such prolonged fraction delivery times on treatment outcome.

METHODS AND MATERIALS

The generalized linear-quadratic (LQ) model, which accounts for sublethal damage repair and clonogen proliferation, was used to calculate the cell-killing efficiency of various simulated and clinical IMRT plans. LQ parameters derived from compiled clinical data for prostate cancer (alpha = 0.15 Gy(-1), alpha/beta = 3.1 Gy, and a 16-min repair half-time) were used to compute changes in the equivalent uniform dose (EUD) and tumor control probability (TCP) due to prolonged delivery time of IMRT as compared with conventional EBRT. EUD and TCP calculations were also evaluated for a wide range of radiosensitivity parameters. The effects of fraction delivery times ranging from 0 to 45 min on cell killing were studied.

RESULTS

Our calculations indicate that fraction delivery times in the range of 15-45 min may significantly decrease cell killing. For a prescription dose of 81 Gy in 1.8 Gy fractions, the EUD for prostate cancer decreases from 78 Gy for a conventional EBRT to 69 Gy for an IMRT with a fraction delivery time of 30 min. The values of EUD are sensitive to the alpha/beta ratio, the repair half-time, and the fraction delivery time. The instantaneous dose-rate, beam-on time, number of leaf shapes (segments), and leaf-sequencing patterns given the same overall fraction delivery time were found to have negligible effect on cell killing.

CONCLUSIONS

The total time to deliver a single fraction may have a significant impact on IMRT treatment outcome for tumors with a low alpha/beta ratio and a short repair half-time, such as prostate cancer. These effects, if confirmed by clinical studies, should be considered in designing IMRT treatments.

Adapting inverse planning to patient and organ geometrical variation: algorithm and implementation

Image guided radiotherapy has the potential to improve both tumour control and normal tissue sparing by including temporal patient specific geometry information into the adaptive planning process. In this study we present a practical method of image guided adaptive inverse planning based on computed tomography (CT) and portal image feedback during the treatment course. The method is based on a general description of the radiotherapy optimization problem subject to dynamic geometrical variations of the patient/organs. We will demonstrate the feasibility of off-line image feedback into the inverse planning process with the example of three prostate cancer patients. CT and portal images acquired during the early course of the treatment are used to predict the geometrical variation distribution of a patient and to re-optimize the treatment plan accordingly. We will study the convergence of the optimization problem with respect to the number of image measurements and adaptive feedback loops.

From physical dose constraints to equivalent uniform dose constraints in inverse radiotherapy planning

Optimization algorithms in inverse radiotherapy planning need information about the desired dose distribution. Usually the planner defines physical dose constraints for each structure of the treatment plan, either in form of minimum and maximum doses or as dose-volume constraints. The concept of equivalent uniform dose (EUD) was designed to describe dose distributions with a higher clinical relevance. In this paper, we present a method to consider the EUD as an optimization constraint by using the method of projections onto convex sets (POCS). In each iteration of the optimization loop, for the actual dose distribution of an organ that violates an EUD constraint a new dose distribution is calculated that satisfies the EUD constraint, leading to voxel-based physical dose constraints. The new dose distribution is found by projecting the current one onto the convex set of all dose distributions fulfilling the EUD constraint. The algorithm is easy to integrate into existing inverse planning systems, and it allows the planner to choose between physical and EUD constraints separately for each structure. A clinical case of a head and neck tumor is optimized using three different sets of constraints: physical constraints for all structures, physical constraints for the target and EUD constraints for the organs at risk, and EUD constraints for all structures. The results show that the POCS method converges stable and given EUD constraints are reached closely.

[Optimization criteria in intensity-modulated radiotherapy]

The present paper provides an overview on the inverse treatment planning for the assessment of intensity-modulated fields. The problem is to find the optimal dose distribution for given attributes of the irradiated tissue. The attributes of the optimal dose distribution are delineated by an objective function. In practice, models are used that evaluate the physical dose distribution, either directly or through their radiobiological effects. In the simplest case, the squared deviation of the achieved dose distribution is minimized to the prescribed dose distribution. For organs structured in parallel, it is common to introduce dose-volume constraints. Another approach is to optimize a value for the probability of complication-free tumor control. The complication probability for normal tissue, in turn, is a rather complex function. However, using the relative seriality, a simple model can be devised with a certain approximation. Other models of "effective dose" are also presented.

Optimizer convergence and local minima errors and their clinical importance

Two of the errors common in the inverse treatment planning optimization have been investigated. The first error is the optimizer convergence error, which appears because of non-perfect convergence to the global or local solution, usually caused by a non-zero stopping criterion. The second error is the local minima error, which occurs when the objective function is not convex and/or the feasible solution space is not convex. The magnitude of the errors, their relative importance in comparison to other errors as well as their clinical significance in terms of tumour control probability (TCP) and normal tissue complication probability (NTCP) were investigated. Two inherently different optimizers, a stochastic simulated annealing and deterministic gradient method were compared on a clinical example. It was found that for typical optimization the optimizer convergence errors are rather small, especially compared to other convergence errors, e.g., convergence errors due to inaccuracy of the current dose calculation algorithms. This indicates that stopping criteria could often be relaxed leading into optimization speed-ups. The local minima errors were also found to be relatively small and typically in the range of the dose calculation convergence errors. Even for the cases where significantly higher objective function scores were obtained the local minima errors were not significantly higher. Clinical evaluation of the optimizer convergence error showed good correlation between the convergence of the clinical TCP or NTCP measures and convergence of the physical dose distribution. On the other hand, the local minima errors resulted in significantly different TCP or NTCP values (up to a factor of 2) indicating clinical importance of the local minima produced by physical optimization.

The design and testing of novel clinical parameters for dose comparison

PURPOSE

New multidimensional dose comparison parameters, normalized agreement test (NAT) values and the NAT index, are introduced and compared with an ideal dose comparison parameter. In this article, we analyze a clinically based two-dimensional (2D) quantitative dose comparison case using a wide range of new and old comparison tools. In doing so, we address the benefits and limitations of many common dose comparison tools.

METHODS AND MATERIALS

An in-house software program was developed using the MATLAB 6.5 programming language. Using this software, several 2D quantitative dose comparison parameters were calculated for the computed and measured dose distributions in an intensity-modulated radiotherapy (IMRT) prostate cancer treatment. The experiences gained in the design and testing of this software program form the basis of the dose comparison tool analysis.

RESULTS

Each dose comparison tool has unique strengths and weaknesses. The underlying assumptions of the NAT values and NAT index lead to acceptable generalized behavior, but are not always valid.

CONCLUSION

A thorough 2D quantitative dose comparison analysis can only be accomplished through the use of many dose comparison tools. The introduction of the NAT index allows a 2D dose comparison to be reduced to a single value, and is thus ideal for setting clinical acceptance criteria for IMRT verifications.

[Prostate brachytherapy: current states and future prospects]

The paper presents the characteristics, the place and the limits of brachytherapy in prostate radiotherapy. While sparing the rectal wall, erectile function as well as urinary continence, I(125) and Pd(103) permanent implants represent interesting approaches for good prognosis tumours in comparison to surgery or conformal external beam radiotherapy with similar cure rates. Overcoming easily the problems of organ motion and patient positioning while allowing doses per fraction as high as 10 Gy, brachytherapy is an excellent boosting method in the treatment of intermediate or unfavourable prognosis tumours of which alpha/beta is 1,5 Gy. Encouraging biological control rates of 80-90% have been published in phase II trials. Compared to external beam radiotherapy, the heterogeneity of irradiation inside the clinical target volume should increase the probability of cure as for a specific dose, a significant part will be overdosed. So far, 120-130% of the prescribed doses are delivered to the peripheral zone at the origin of 70% of tumours. On the opposite, this heterogeneity is inducing an overdosage of the urethral bed at the price of higher toxicity levels in situations of previous obstructive syndrome and urethral stenosis. A better integration of the therapeutic modalities available, brachytherapy included, should increase our curative possibilities in the radiation treatment of prostatic cancer.

Radiotherapy of small intracranial tumours with different advanced techniques using photon and proton beams: a treatment planning study

BACKGROUND AND PURPOSE

The potential benefits and limitations of five different radiation techniques, 3D conformal radiotherapy (3DCRT), stereotactic arc therapy (SRS/T), intensity modulated radiotherapy with photons (IMRT), and radiotherapy with protons (spot scanning (SSp) or passive scattering (PSp)), have been assessed using comparative treatment planning methods in a cohort of patients presenting with 'benign' brain tumours.

MATERIAL AND METHODS

Plans for five acoustic neurinomas, five meningiomas, and two pituitary adenomas were computed for all modalities using computed tomography (CT) scans to delineate planning target volume and organs at risk (OARs) and to predict dose distributions. Dose-volume histograms were used for physical and simple biological evaluation.

RESULTS

Proton techniques were shown to be superior to all photon approaches for the irradiation of small brain lesions in terms of target dose uniformity and conformity and in terms of sparing OARs. No major differences were observed between the results of the photon techniques, which were generally good for target coverage. Minimum target doses ranged from 81% with SRS/T to 93% with IMRT. The volume receiving more than 95% of the dose ranged from 95% (SRS/T) to 99% (PSp). No clear patterns of coverage dependence upon target shape were observed. Maximum brain stem irradiation ranged from 60% with IMRT to 26% with protons and the conformity index from 4.4 with IMRT to 2.5 with protons. Considering the rather long life expectancy of the patients suffering from meningiomas, neurinomas, and pituitary adenomas, the most important aspect to be considered, other than target coverage, is toxicity and in the long term, the possibility of secondary tumour induction. Considering these aspects, proton irradiation should be the irradiation technique of choice, when available. If not, IMRT, or even 3DCRT, techniques can provide an acceptable compromise, even without recurring to unconventional treatments like SRS/T, which require complex installations and high machine occupancy.

Biological response to radiation therapy

In an investigation by the Swedish Cancer Society, the present status, critical issues and future aspects and potentials were described by an expert group for each of nine major areas of radiation therapy research. This article deals with biological response to radiation. Separate sections deal with molecular responses to radiation, the stem cell and clonogenic cell concepts and the importance of cell proliferation, cell and tissue responses to doses above and below 1 Gy, respectively, the potential role of intercellular signalling pathways, the so-called bystander effect and radiation biology-based therapy planning and treatment optimization.

Radiation therapy dose delivery

In an investigation by the Swedish Cancer Society, the present status critical issues and future aspects and potentials in each of nine major areas of radiation therapy research were described by an expert group. The report presented here deals with radiation therapy dose delivery. dose distributions, beam shaping and intensity modulation.

[Fluence-modulated radiotherapy with an optimization-integrated sequencer]

On the basis of two clinical cases, we present fluence-modulated radiotherapy with a sequencer integrated into the optimization of our treatment-planning software HYPERION. In each case, we achieved simple relations for the dependence of the total number of segments on the complexity of the sequencing, as well as for the dependence of the dose-distribution quality on the number of segments. For both clinical cases, it was possible to obtain treatment plans that complied with the clinical demands on dose distribution and number of segments. Also, compared to the widespread concept of equidistant steps, our method of sequencing with fluence steps of variable size led to a significant reduction of the number of segments, while maintaining the quality of the dose distribution. Our findings substantiate the value of the integration of the sequencer into the optimization for the clinical efficiency of IMRT.

The method of intercepts in parameter space for the analysis of local minima caused by dose-volume constraints

The local minima problem in radiotherapy optimization has been a concern for both researchers and physicians. In this work, local minima induced by dose-volume histogram (DVH) constraints are discussed. The non-convex property of the feasible set formed by DVH constraints is discussed in beam weight space. An intuitive explanation of the origin of this type of local minima is given by a two-beam model setup. Some interesting properties and insights about the DVH-induced local minima are found. Based on this, a heuristic non-random initial guess sampling method is proposed and applied to a clinical nasopharyngeal case, where some significantly different local minima are located.

The theoretical benefit of beam fringe compensation and field size reduction for iso-normal tissue complication probability dose escalation in radiotherapy of lung cancer

To assess the benefit of beam fringe (50%-90% dose level) sharpening for lung tumors, we performed a numerical simulation in which all geometrical errors (breathing motion, random and systematic errors) are included. A 50 mm diameter lung tumor, located centrally in a lung-equivalent phantom was modeled. Treatment plans were designed with varying number and direction of beams, both with and without the use of intensity modulation to sharpen the beam fringe. Field size and prescribed dose were varied under the constraint of a constant mean lung dose of 20 Gy, which yields a predicted complication probability of about 10%. After numerical simulation of the effect of setup errors and breathing, the resulting dose distribution was evaluated using the minimum dose and the equivalent uniform dose (EUD) in the moving clinical target volume (CTV). When the dose in the CTV was constrained between 95% and 107% of the prescribed dose, the maximum attainable EUD was 71 Gy for a four-field noncoplanar technique with simple conformal beams. When penumbra sharpening was applied using a single beam segment at the edge of the open field, this EUD could be raised to 87 Gy. For a hypothetical infinitely steep penumbra, further escalation to an EUD of 104 Gy was possible. When the dose in the CTV was not constrained, a large escalation of the EUD was possible compared to the constrained case. In this case, the maximum attainable EUD for open fields was 115 Gy, using the four-field noncoplanar technique. The benefit of penumbra sharpening was only modest, with no increase of the EUD for the single-segment technique and a small increase to 125 Gy for the infinitely steep penumbra. From these results we conclude that beam fringe sharpening in combination with field-size reduction leads to a large increase in EUD when a homogeneous target dose is pursued. Further escalation of the EUD is possible when the homogeneity constrained is relaxed, but the relative benefit of beam-fringe sharpening then decreases.

Simultaneous integrated boost intensity-modulated radiotherapy for locally advanced head-and-neck squamous cell carcinomas. I: dosimetric results

PURPOSE

This report describes the dosimetric analyses of a Phase I/II protocol, designed to examine the capabilities of an institutionally developed intensity-modulated radiotherapy (IMRT) system with respect to dose escalation. The protocol employed stringent dosimetric guidelines in the treatment of locally advanced head-and-neck squamous cell carcinomas (HNSCC) with radiotherapy alone using IMRT and the simultaneous integrated boost (SIB) technique.

METHODS AND MATERIALS

The first 14 patients enrolled on the protocol were included in this analysis. Escalating doses of 68.1 Gy (6 patients), 70.8 Gy (6 patients), and 73.8 Gy (2 patients) were delivered to the gross tumor volume (GTV) in 30 fractions. Simultaneously, constant dose coverage was given to the subclinical disease and the electively treated nodal regions, which received 60 Gy and 54 Gy, respectively, in all three cohorts. Parotid glands were spared to the degree possible without compromising target coverage. The following indices are reported for the GTV: (1) dose to specified percent volumes (e.g., D(98) and D(2)); (2) homogeneity index defined as the ratio (D(2) - D(98))/D(prescription); (3) biologically equivalent uniform dose (EUD); and (4) an index of conformality, PITV, defined as the ratio of volume enclosed within the prescribed isodose surface to the target volume. Treatments were planned and delivered with nine 6-MV photon beams using the multileaf collimator (MLC) "sliding window" technique.

RESULTS

Mean doses to 98% of GTV were 68.4 Gy, 70.5 Gy, and 70.8 Gy, and average GTV dose homogeneity was 6.7%, 7.6%, and 8.8% for the three cohorts. The average doses to the parotid gland proximal to and distant from GTV were 41.3 Gy and 25.7 Gy, respectively. Dose distributions measured in phantom showed good agreement with calculations.

CONCLUSIONS

Treatment of locally advanced HNSCC using SIB-IMRT as described is feasible. Treatment planning and delivery are safer and more efficient than with conventional three-dimensional processes. Predicted dose distributions can be accurately delivered with excellent conformality using dynamic MLC. At least one of the parotid glands can be adequately spared. Patient follow-up continues and will allow eventual quantitative correlation of delivered dose distributions with clinical outcomes.

A software tool for specifying voxel models for dosimetry estimation

Many dose estimation problems can be conveniently formulated in terms of finding the energy emitted and absorbed by a set of homogeneous volume elements (voxels) arranged in a rectilinear grid. The solution of these problems requires an accurate model of the source and target geometry to be established, whereupon conventional Monte Carlo simulation of radiation transport can be employed to determine energy deposition. A software application ("MrVoxel") has been developed to assist in the specification of the source and target models. This application includes tools for image segmentation and image registration (2D and 3D, intra- and inter-modality, interactive, and automatic). It employs a plug-in architecture to facilitate customization and future expansion: plug-ins can be written to perform image import and export as well as to implement specialized image processing routines. Using plug-ins, the package can, for example, import DICOM 3.10 files and export input files for a voxel-based Monte Carlo package. Standard dosimetric tools such as the geometric mean method, transmission based attenuation correction, and MIRD-style voxel dose kernel convolution are also implemented as plug-ins. MrVoxel was implemented on a Macintosh computer using a commercial software framework to produce a conventional document-centric application. Hence it includes useful features such as the ability to undo an operation or to save a processed image set at any point. This latter feature enables the production of a processing trail, to allow post-hoc auditing of the analysis process. This paper describes the MrVoxel application and its role in the analysis of a particular dosimetry problem.

Biological-effective versus conventional dose volume histograms correlated with late genitourinary and gastrointestinal toxicity after external beam radiotherapy for prostate cancer: a matched pair analysis

BACKGROUND

To determine whether the dose-volume histograms (DVH's) for the rectum and bladder constructed using biological-effective dose (BED-DVH's) better correlate with late gastrointestinal (GI) and genitourinary (GU) toxicity after treatment with external beam radiotherapy for prostate cancer than conventional DVH's (C-DVH's).

METHODS

The charts of 190 patients treated with external beam radiotherapy with a minimum follow-up of 2 years were reviewed. Six patients (3.2%) were found to have RTOG grade 3 GI toxicity, and similarly 6 patients (3.2%) were found to have RTOG grade 3 GU toxicity. Average late C-DVH's and BED-DVH's of the bladder and rectum were computed for these patients as well as for matched-pair control patients. For each matched pair the following measures of normalized difference in the DVH's were computed: (a) deltaAUC = (Area Under Curve [AUC] in grade 3 patient--AUC in grade 0 patient)/(AUC in grade 0 patient) and (b) deltaV60 = (Percent volume receiving = 60 Gy [V60] in grade 3 patient--V60 in grade 0 patient)/(V60 in grade 0 patient).

RESULTS

As expected, the grade 3 curve is to the right of and above the grade 0 curve for all four sets of average DVH's--suggesting that both the C-DVH and the BED-DVH can be used for predicting late toxicity. deltaAUC was higher for the BED-DVH's than for the C-DVH's--0.27 vs 0.23 (p = 0.036) for the rectum and 0.24 vs 0.20 (p = 0.065) for the bladder. deltaV60 was also higher for the BED-DVH's than for the C-DVH's--2.73 vs 1.49 for the rectum (p = 0.021) and 1.64 vs 0.71 (p = 0.021) for the bladder.

CONCLUSIONS

When considering well-established dosimetric endpoints used in evaluating treatment plans, BED-DVH's for the rectum and bladder correlate better with late toxicity than C-DVH's and should be considered when attempting to minimize late GI and GU toxicity after external beam radiotherapy for prostate cancer.

How should we describe the radioblologic effect of extracranial stereotactic radlosurgery: equivalent uniform dose or tumor control probability?

Extracranial stereotactic radiosurgery (ESR) is now undergoing clinical investigation at numerous institutions as a treatment for solitary malignant lesions. Because there is no standard ESR technique, the same minimum dose might be applied through widely variable target dose-volume histograms. For multicenter trials of ESR or interinstitutional comparisons, a reliable index of radiobiological dose equivalency might facilitate the evaluation of dose-response relationships. Equivalent uniform dose (EUD) and tumor control probability (TCP) were considered for this application. While EUD appears more robust for the prospective description of ESR, TCP is expected to remain more valuable for a post hoc estimation of radiosensitivity parameters.

Intensity-modulated radiation therapy (IMRT) for locally advanced paranasal sinus tumors: incorporating clinical decisions in the optimization process

PURPOSE

Intensity-modulated radiotherapy (IMRT) plans require decisions about priorities and tradeoffs among competing goals. This study evaluates the incorporation of various clinical decisions into the optimization system, using locally advanced paranasal sinus tumors as a model.

METHODS AND MATERIALS

Thirteen patients with locally advanced paranasal sinus tumors were retrospectively replanned using inverse planning. Two clinical decisions were assumed: (1) Spare both optic pathways (OP), or (2) Spare only the contralateral OP. In each case, adequate tumor coverage (treated to 70 Gy in 35 fractions) was required. Two beamlet IMRT plans were thus developed for each patient using a class solution cost function. By altering one key variable at a time, different levels of risk of OP toxicity and planning target volume (PTV) compromise were compared in a systematic manner. The resulting clinical tradeoffs were analyzed using dosimetric criteria, tumor control probability (TCP), equivalent uniform dose (EUD), and normal tissue complication probability.

RESULTS

Plan comparisons representing the two clinical decisions (sparing both OP and sparing only the contralateral OP), with respect to minimum dose, TCP, V(95), and EUD, demonstrated small, yet statistically significant, differences. However, when individual cases were analyzed further, significant PTV underdosage (>5%) was present in most cases for plans sparing both OP. In 6/13 cases (46%), PTV underdosage was between 5% and 15%, and in 3 cases (23%) was greater than 15%. By comparison, adequate PTV coverage was present in 8/13 cases (62%) for plans sparing only the contralateral OP. Mean target EUD comparisons between the two plans (including 9 cases where a clinical tradeoff between PTV coverage and OP sparing was required) were similar: 68.6 Gy and 69.1 Gy, respectively (p = 0.02). Mean TCP values for those 9 cases were 56.5 vs. 61.7, respectively (p = 0.006).

CONCLUSIONS

In IMRT plans for paranasal sinus tumors, tradeoff values between OP toxicity and PTV coverage can be compared for different clinical decisions. The information derived can then be used to individualize the parameters within the optimization system. This process of determining clinical tradeoffs associated with different clinical decisions may be a useful tool in other sites.

The influence of brachytherapy dose heterogeneity on estimates of alpha/beta for prostate cancer

The sensitivity of estimates of alpha/beta for prostate tumours to dose heterogeneity in 125I brachytherapy implants, as well as to variation in selected radiobiological parameters, is analysed. The tumour control probabilities of brachytherapy and external beam radiotherapy are equated for ranges of alpha, Tpot, RBE and external beam dose. For each combination of parameters, the equality is used to derive the value of alpha/beta. Different clinical (non-uniform) brachytherapy dose distributions, and three uniform brachytherapy dose distributions (120, 144 and 160 Gy) are used. For 'nominal' input parameter values of Tpot = 45 days, alpha = 0.2 Gy(-1), RBE = 1.4, and an external beam dose of 70 Gy, the values obtained for alpha/beta ranged between 2.1 and 12.3 Gy for all of the clinical DVHs, between 2.1 and 3.8 Gy for the better quality clinical implants and between 1.0 and 1.8 Gy for the uniform brachytherapy doses. When only 2% of the volume receiving the lowest dose is omitted from the clinical DVHs, the estimated alpha/beta values ranged between 1.4 and 2.1 Gy. When ranges of input parameters were also considered, the overall range of alpha/beta values for the clinical brachytherapy dose distributions lay between 1.1 and 12.3 Gy for the three best clinical implants, and between 0.7 and 6.3 Gy for uniform doses. We conclude that estimation of alpha/beta without taking into account dose heterogeneity and inter-patient variation may underestimate the actual value alpha/beta.

Intensity-modulated radiotherapy optimization with gEUD-guided dose-volume objectives

Currently, most intensity-modulated radiation therapy systems use dose-volume (DV)-based objectives. Although acceptable plans can be generated using these objectives, much trial and error is necessary to plan complex cases with many structures because numerous parameters need to be adjusted. An objective function that makes use of a generalized equivalent uniform dose (gEUD) was developed recently that has the advantage of involving simple formulae and fewer parameters. In addition, not only does the gEUD-based optimization provide the same coverage of the target, it provides significantly better protection of critical structures. However, gEUD-based optimization may not be superior once dose distributions and dose-volume histograms (DVHs) are used to evaluate the plan. Moreover, it is difficult to fine-tune the DVH with gEUD-based optimization. In this paper, we propose a method for combining the gEUD-based and DV-based optimization approaches to overcome these limitations. In this method, the gEUD optimization is performed initially to search for a solution that meets or exceeds most of the treatment objectives. Depending on the requirements, DV-based optimization with a gradient technique is then used to fine-tune the DVHs. The DV constraints are specified according to the gEUD plan, and the initial intensities are obtained from the gEUD plan as well. We demonstrated this technique in two clinical cases: aprostate cancer and ahead and neck cancer case. Compared with the DV-optimized plan, the gEUD plan provided better protection of critical structures and the target coverage was similar. However, homogeneities were slightly poorer. The gEUD plan was then fine-tuned with DV constraints, and the resulting plan was superior to the other plans in terms of the dose distributions. The planning time was significantly reduced as well. This technique is an effective means of optimizing individualized treatment plans.

Reduction in radiation dose to lung and other normal tissues using helical tomotherapy to treat lung cancer, in comparison to conventional field arrangements

The purpose of this study was to determine whether the use of tomotherapy in the treatment of non-small-cell lung cancer (NSCLC) has the potential to reduce radiation dose to normal tissues, in particular, the lungs, esophagus, and spinal cord, as compared with standard radiotherapy. Five patients with anatomically or physiologically inoperable stage III NSCLC were studied, representing a variety of tumor sizes and locations. For each patient, two treatment plans were generated. One was developed using conventional field arrangements (CFA), and the other for tomotherapy. Using dose-volume histogram reduction techniques, including mean normalized dose (NTDmean), V20, and effective uniform dose (EUD), the normal tissue doses for CFA and tomotherapy plans for a given fixed tumor dose were compared. In addition, the maximum tumor doses possible for a given level of mean normalized lung dose were computed and compared for the CFA and tomotherapy plans. The gross tumor volumes in the five patients studied ranged from 13.5 to 87.1 cm. The tumor dose distributions, determined by EUD and minimum dose, were similar for both CFA and tomotherapy plans, as intended. In all cases, the NTDmean of both lungs was significantly reduced using tomotherapy planning (range: 10-53% reduction, mean: 31%). The volume of lung receiving more than 20 Gy was also reduced in all cases using tomotherapy (range: 17-37% reduction, mean: 22%). For a constant lung NTDmean, it is shown that it should be possible to increase tumor dose to up to 160 Gy in certain patients with tomotherapy. The dose to the spinal cord and esophagus was also reduced in all cases with tomotherapy planning, compared with plans generated using conventional field arrangements. Both tomotherapy, and to a lesser extent conventional three-dimensional conformal radiotherapy, have the potential to significantly decrease radiation dose to lung and other normal structures in the treatment of NSCLC. This has important implications for dose escalation strategies in the future.