Causes and consequences of inhomogeneous dose distributions in radiation therapy

Rethinking the potential role of dose painting in personalized ultra-fractionated stereotactic adaptive radiotherapy

Fractionated radiotherapy was established in the 1920s based upon two principles: (1) delivering daily treatments of equal quantity, unless the clinical situation requires adjustment, and (2) defining a specific treatment period to deliver a total dosage. Modern fractionated radiotherapy continues to adhere to these century-old principles, despite significant advancements in our understanding of radiobiology. At UT Southwestern, we are exploring a novel treatment approach called PULSAR (Personalized Ultra-Fractionated Stereotactic Adaptive Radiotherapy). This method involves administering tumoricidal doses in a pulse mode with extended intervals, typically spanning weeks or even a month. Extended intervals permit substantial recovery of normal tissues and afford the tumor and tumor microenvironment ample time to undergo significant changes, enabling more meaningful adaptation in response to the evolving characteristics of the tumor. The notion of dose painting in the realm of radiation therapy has long been a subject of contention. The debate primarily revolves around its clinical effectiveness and optimal methods of implementation. In this perspective, we discuss two facets concerning the potential integration of dose painting with PULSAR, along with several practical considerations. If successful, the combination of the two may not only provide another level of personal adaptation ("adaptive dose painting"), but also contribute to the establishment of a timely feedback loop throughout the treatment process. To substantiate our perspective, we conducted a fundamental modeling study focusing on PET-guided dose painting, incorporating tumor heterogeneity and tumor control probability (TCP).

Improving organ at risk sparing in oropharyngeal treatment planning by increasing target dose heterogeneity: A feasibility study

Target dose homogeneity has historically been a priority in radiotherapy treatment planning. However, in an era of more advanced modulated techniques, there is now greater flexibility in shaping dose distributions suggesting that allowing controlled target dose heterogeneity may consequently improve organ at risk (OAR) sparing. This study sought to determine the feasibility of allowing an increase in target dose heterogeneity in oropharyngeal VMAT plans, and to examine the dosimetric impact this has on target coverage and OARs such as the parotid glands, spinal cord, brainstem and mandible. Nineteen oropharyngeal patients' plans were created with homogeneous dose distributions specified in the London Cancer Head and Neck Radiotherapy Protocol. The upper dose constraint (UDC) objective of the primary planning target volumes (PTV) for each plan were increased in increments of 10% until a maximum of 150% of the prescribed dose was reached. These plans were dosimetrically compared to plans with a uniform dose distribution in terms of OAR sparing and target coverage. Minimal coverage was not compromised, with the largest median changes being a 0.81% decrease [98.6 to 97.8%] to the PTV_70Gy D98% and a 2.86% decrease [99.81 to 96.96%] to the PTV_54Gy D98% at a UDC of 150% of the prescription dose. An OAR sparing effect was observed for the parotid glands, spinal cord and oral cavity sub PTV. Mandible and brainstem Dmax values increased as the PTV UDC increased. Changes in brainstem dose were not statistically significant. All other differences were statistically significant for UDC's above 130%. Target coverage was not compromised as a result of increased target dose heterogeneity. The OAR sparing effect was promising for most organs, however further research with a larger dataset is necessary surrounding the effect on organs that overlap with the PTV.

Impact of target dose inhomogeneity on BED and EUD in lung SBRT

PURPOSE

To evaluate the effect of dose heterogeneity in the treatment target on biologically effective dose (BED) for frequently used hypofractionation regimens in stereotactic body radiation therapy (SBRT).

METHODS

In the case of non-uniform target dose, BED in the planning target volume (PTV) is determined by using the linear-quadratic model. An expression for BED is obtained for an arbitrary dose distribution in the PTV in the case of small variance of the target dose. Another analytical expression for BED is obtained by assuming a Gaussian dose distribution in the target.

RESULTS

Analytical expressions for BED as a function of the variance of the target dose have been derived. It is shown that a relatively small dose inhomogeneity (<5%-6%) can cause a significant reduction (i.e. >10%) in the corresponding BED and equivalent uniform dose (EUD) compared to the case of uniform target dose.

CONCLUSIONS

Small variations in the absorbed dose can significantly reduce BED and EUD in the PTV. The effect of dose non-uniformity on BED increases with increasing dose per fraction. The observed reduction in BED compared to that for uniform target dose can be several times greater for SBRT than for standard fractionation with dose per fraction varying between 1.8 and 2 Gy.

Simultaneous Integrated Boost Intensity-Modulated Radiotherapy for Locally Advanced Drug-Resistant Gastrointestinal Stromal Tumors: A Feasibility Study

Background

As an emerging clinical problem, locally advanced drug-resistant gastrointestinal stromal tumors (LADRGISTs) has relatively few therapeutic schemes. Although radiotherapy is not often considered for GISTs, it could be a valuable contributing modality. The aim of our study is to explore a safe and effective radiation regimen for LADR-GISTs.

Methods

Three patients with LADR-GISTs were treated with simultaneous integrated boost intensity-modulated radiation therapy (SIB-IMRT) plans. In the SIB-IMRT plans, gross target volume (GTV) was divided into GTV-outer, GTV-mid, and GTV-center. And the prescribed dose of planning gross target volume (PGTV) and GTV-outer were both set to 50.4 Gy in 28 fractions. GTV-mid and GTV-center were simultaneously boosted to 60–62 Gy and 62–64 Gy respectively. For comparison purposes, conventional IMRT (Con-IMRT) plans with uniform dose distribution were generated for same optimization objectives without a dose boost to GTV-mid and GTV-center. All plans were optimized to make sure that deliver at least 95% of the prescription dose was delivered to PGTV. Isodose distribution, dose profiles, conformity indexes (CIs), monitor units (MUs), and dose volume histogram (DVH) was evaluated for each individual patient. After the three patients were treated with SIB-IMRT plans, the relative changes in the tumor size and CT values by CT scanning were also tracked.

Results

Compared with Con-IMRT plans, SIB-IMRT plans saw a significant increase from D₉₅ to D₂ of the GTV. With steeper dose gradients in the dose profiles, SIB-IMRT plans had GTV-mid and GTV-center accumulated with higher dose mainly by delivering extra 93 MUs in average. However, there was no significant difference in CIs and organs at risks (OARs) DVH. The relative changes in tumor size and CT values of the three patients in follow up were up to the Choi criteria and the three patients were all assessed as partial response.

Conclusions

The proposed SIB-IMRT may be a potential technique for achieving objective response and prolonging survival of selected GISTs patients.

A planning study of focal dose escalations to multiparametric MRI-defined dominant intraprostatic lesions in prostate proton radiation therapy

OBJECTIVES

The purpose of this study is to investigate the dosimetric effect and clinical impact of delivering a focal radiotherapy boost dose to multiparametric MRI (mp-MRI)-defined dominant intraprostatic lesions (DILs) in prostate cancer using proton therapy.

METHODS

We retrospectively investigated 36 patients with pre-treatment mp-MRI and CT images who were treated using pencil beam scanning (PBS) proton radiation therapy to the whole prostate. DILs were contoured on co-registered mp-MRIs. Simultaneous integrated boost (SIB) plans using intensity-modulated proton therapy (IMPT) were created based on conventional whole-prostate-irradiation for each patient and optimized with additional DIL coverage goals and urethral constraints. DIL dose coverage and organ-at-risk (OAR) sparing were compared between conventional and SIB plans. Tumor control probability (TCP) and normal tissue complication probability (NTCP) were estimated to evaluate the clinical impact of the SIB plans.

RESULTS

Optimized SIB plans significantly escalated the dose to DILs while meeting OAR constraints. SIB plans were able to achieve 125, 150 and 175% of prescription dose coverage in 74, 54 and 17% of 36 patients, respectively. This was modeled to result in an increase in DIL TCP by 7.3-13.3% depending on α / β and DIL risk level.

CONCLUSION

The proposed mp-MRI-guided DIL boost using proton radiation therapy is feasible without violating OAR constraints and demonstrates a potential clinical benefit by improving DIL TCP. This retrospective study suggested the use of IMPT-based DIL SIB may represent a strategy to improve tumor control.

ADVANCES IN KNOWLEDGE

This study investigated the planning of mp-MRI-guided DIL boost in prostate proton radiation therapy and estimated its clinical impact with respect to TCP and NTCP.

A new homogeneity index definition for evaluation of radiotherapy plans

Purpose

The goal of this study was to define a new homogeneity index (HI) to evaluate dose homogeneity within a target volume.

Materials and Methods

The new HI is based on the area under an ideal dose‐volume histogram curve (IA), the area under the achieved dose‐volume histogram curve (AA), and the overlapping area between the IA and AA (OA). It is defined as the ratio of the square of OA to the product of the IA and AA. To evaluate the performance of the new HI, 88 cases were selected and two plans were designed for each case. The homogeneity of the two plans was first evaluated by three physicists, with their judgments forming the evaluation standard and then evaluated by the new HI and other HIs of D_max/D_p, D₅/D₉₅, (D₂ − D₉₈)/D_p, (D₂ − D₉₈)/D₅₀ and S‐index. An evaluation was determined to be accurate if its result was agreed upon by physicists. The percentage accuracy of evaluation was calculated as the ratio of the number of accurate evaluations to the total number of evaluations. Pearson's chi‐square test was performed for statistical analysis.

Results

The percentage accuracies of the new HI, D_max/D_p, D₅/D₉₅, (D₂ − D₉₈)/D_p, (D₂ − D₉₈)/D₅₀, and S‐index were 98.51%, 88.80%, 94.78%, 94.78%, 96.27%, and 97.01%, respectively. The newly defined HI had the highest accuracy of all the HIs, with the difference being statistically significant (P < 0.05).

Conclusions

The newly defined HI was shown to be effective in the evaluation of dose homogeneity, and we recommended it for evaluating the homogeneity of radiotherapy plans.

Trade-off between the conflicting planning goals in correlation with patient’s anatomical parameters for intensity-modulated radiotherapy of prostate cancer patients

Aim

To quantify the relationship between the planning target volume (PTV) dose homogeneity and organs at risk (OARs) sparing in correlation with anatomical parameters in prostate intensity-modulated radiotherapy (IMRT).

Materials and methods

Nine IMRT plans with various target dose constraints’ priorities were created for 15 prostate cancer patients. Selected PTV and OARs parameters were calculated for the patients. A trade-off was assessed between homogeneity index (HI) and OAR sparing. Several anatomical parameters were evaluated to investigate their effects on the OAR sparing and HI.

Results

Inverse exponential relationships were found between the OAR sparing and HI (average R 2 of 0·983 and 0·994 for bladder and rectum, respectively). Decreasing the priority led to more OARs sparing (normal tissue complication probability reduction: 97·6 and 74·5%; mean dose reduction: 16·3 and 11·3% for bladder and rectum, respectively) and worsening of the HI (0·095–0·322) but with no significant effect on tumour control probability. Furthermore, OARs volumes, distances between OARs and PTV and their joint volumes had stronger correlations with OARs’ mean doses.

Conclusion

Enforcement of target dose constraints was more effective on the improvement of HIs for the patients with initial high HI values at low dose constraints’ priorities. Reducing the priority had more effects on the OARs sparing compared to HI, especially for the patients with high OAR doses in high priority plans. This can be attributed to smaller distances or greater joint volumes between the OARs and PTV.

A quantitative assessment of the consequences of allowing dose heterogeneity in prostate radiation therapy planning

Target dose uniformity has been historically an aim of volumetric modulated arc therapy (VMAT) planning. However, for some sites, this may not be strictly necessary and removing this constraint could theoretically improve organ‐at‐risk (OAR) sparing and tumor control probability (TCP). This study systematically investigates the consequences of PTV dose uniformity that results from the application or removal of an upper dose constraint (UDC) in the inverse planning process for prostate VMAT treatments. OAR sparing, target coverage, hotspots, and plan complexity were compared between prostate VMAT plans with and without the PTV UDC optimized using the progressive resolution optimizer (PRO, Varian Medical Systems, Palo Alto, CA). Removing the PTV UDC, the median D1cc reached 144.6% for the CTV and the PTV, and an average increase of 3.2% TCP was demonstrated, while CTV and PTV coverage evaluated by D99% was decreased by less than 0.6% with statistical significance. Moreover, systematic improvement in the rectum dose volume histograms was shown (a 5–10% decrease in the volume receiving 50% to 75% prescribed dose), resulting in an average decrease of 1.3% (P < 0.01) in the rectum normal tissue complication probability. Additional consequences included potentially increased dose to the urethra as evaluated by PTV D0.035cc (median: 153.4%), delivering 283 extra monitor units (MUs), and slightly higher degrees of modulation. In general, the results were consistent when a different optimizer (Photon Optimizer, Varian Medical Systems) was used. In conclusion, removing the PTV UDC is acceptable for localized prostate cases given the systematic improvement of rectal dose and TCP. It can be particularly useful for cases that do not meet the rectum dose constraints with the PTV UDC on. This comes with the foreseeable consequences of increased dose heterogeneity in the PTV and an increase in MUs and plan complexity. It also has a higher requirement for reproducing the position and size of the target and OARs during treatment. Finally, with the PTV UDC completely removed, in some cases the maximum doses within the PTV did approach levels that may be of concern for urethral toxicity and therefore in clinical implementation it may still be necessary to include a PTV UDC, but one based on limiting toxicity rather than enforcing dose homogeneity.

A novel voxel based homogeneity index: Rationale and clinical implications for whole-brain radiation therapy

PURPOSE OR OBJECTIVE

A homogeneity index (HI) measures the uniformity of a dose distribution within a given target volume. Traditional HIs only use a limited number of dose-volume histogram data-points for calculation. A voxel-based homogeneity index (VHI) is proposed which utilizes the entire information of the three-dimensional dose distribution. We compared the VHI with existing HIs and analyzed if VHI results were associated with treatment outcomes in patients who underwent therapeutic WBRT.

MATERIAL AND METHODS

The VHI analyzes deviations from the prescribed dose in each voxel of the target volume. We retrospectively analyzed WBRT treatment plans. Overall survival (OS), CNS progression-free-survival (CNS PFS) and hazard rates were compared for tertile-split levels of the VHI using the Kaplan-Meier methods and multivariable Cox-regression analysis.

RESULTS

WBRT treatment plans (n = 770) were used for HIs comparison. OS and CNS PFS were assessed for 430 patients. The VHI showed a higher sensitivity for dose inhomogeneities. Lower OS and CNS PFS were observed for higher levels of VHI_Underdosage, particularly in patients with good performance status (KPS >70%) (OS: Log-rank P = .007, HR = 1.37 95%CI [1.09, 1.72]).

CONCLUSION

Higher sensitivity and feasibility to assess treatment plan quality using the VHI were demonstrated. First clinical implications were found in terms of compromised OS/CNS PFS for WBRT with radiation underdosage.

Potential implications of the bystander effect on TCP and EUD when considering target volume dose heterogeneity

PURPOSE

In light of in vitro evidence suggesting that radiation-induced bystander effects may enhance non-local cell killing, there is potential for impact on radiotherapy treatment planning paradigms such as the goal of delivering a uniform dose throughout the clinical target volume (CTV). This work applies a bystander effect model to calculate equivalent uniform dose (EUD) and tumor control probability (TCP) for external beam prostate treatment and compares the results with a more common model where local response is dictated exclusively by local absorbed dose. The broad assumptions applied in the bystander effect model are intended to place an upper limit on the extent of the results in a clinical context.

MATERIALS AND METHODS

EUD and TCP of a prostate cancer target volume under conditions of increasing dose heterogeneity were calculated using two models: One incorporating bystander effects derived from previously published in vitro bystander data ( McMahon et al. 2012 , 2013a); and one using a common linear-quadratic (LQ) response that relies exclusively on local absorbed dose. Dose through the CTV was modelled as a normal distribution, where the degree of heterogeneity was then dictated by changing the standard deviation (SD). Also, a representative clinical dose distribution was examined as cold (low dose) sub-volumes were systematically introduced.

RESULTS

The bystander model suggests a moderate degree of dose heterogeneity throughout a target volume will yield as good or better outcome compared to a uniform dose in terms of EUD and TCP. For a typical intermediate risk prostate prescription of 78 Gy over 39 fractions maxima in EUD and TCP as a function of increasing SD occurred at SD ∼ 5 Gy. The plots only dropped below the uniform dose values for SD ∼ 10 Gy, almost 13% of the prescribed dose. Small, but potentially significant differences in the outcome metrics between the models were identified in the clinically-derived dose distribution as cold sub-volumes were introduced.

CONCLUSIONS

In terms of EUD and TCP, the bystander model demonstrates the potential to deviate from the common local LQ model predictions as dose heterogeneity through a prostate CTV varies. The results suggest, at least in a limiting sense, the potential for allowing some degree of dose heterogeneity within a CTV, although further investigation of the assumptions of the bystander model are warranted.

Impact of dose and sensitivity heterogeneity on TCP

This present paper presents an analytical description and numerical simulations of the influence of macroscopic intercell dose variations and intercell sensitivity variations on the probability of controlling the tumour. Computer simulations of tumour control probability accounting for heterogeneity in dose and radiation sensitivity were performed. An analytical expression for tumor control probability accounting for heterogeneity in sensitivity was also proposed and validated against simulations. The results show good agreement between numerical simulations and the calculated TCP using the proposed analytical expression for the case of a heterogeneous dose and sensitivity distributions. When the intratumour variations of dose and sensitivity are taken into account, the total dose required for achieving the same level of control as for the case of homogeneous distribution is only slightly higher, the influence of the variations in the two factors taken into account being additive. The results of this study show that the interplay between cell or tumour variation in the sensitivity to radiation and the inherent heterogeneity in dose distribution is highly complex and therefore should be taken into account when predicting the outcome of a given treatment in terms of tumor control probability.

Dosimetric comparison of intensity-modulated radiosurgery and helical tomotherapy for the treatment of multiple intracranial metastases

The purpose of this study was to evaluate the dosimetry of single fraction, single-isocenter intensity-modulated radiosurgery (IMRS) plans for multiple intracranial metastases and to compare Helical Tomotherapy (HT). Ten treatment plans with 3-6 brain metastases treated with IMRS were re-planned with HT. The mean number of lesions was 5 and mean PTV 22 cm(3). The prescribed dose was 16-20 Gy. The mean V100% was similar for IMRS and HT, and the mean conformity index was 1.4, mean Paddick confirmity index was 0.7, and mean MDPD was 1.1 for both. The mean gradient index was similar for both. The mean 50% _isodose volume was 179.2 cm(3) for IMRS and 277.0 cm(3) for HT (p=0.01). The mean maximum doses to organs at risk were lower for IMRS except brainstem and right optic nerve. For brain, the integral dose was 5.1 and 6.8 Gy-kg (p<0.001) and mean dose 4.0 and 5.4 Gy (p<0.001) for IMRS and HT, respectively. The mean treatment times were 23 (IMRS) and 41 (HT) minutes. Conformity and homogeneity indices were equivalent and sparing of the organs at risk was clinically acceptable for both IMRS and HT. Though the gradient index was similar for IMRS and HT, the mean 50% isodose volume and integral dose to normal brain were lower for IMRS as was treatment time.

A heterogeneous dose distribution in simultaneous integrated boost: the role of the clonogenic cell density on the tumor control probability

IMRT with inverse planning allows simultaneous integrated boost strategies that exploit the heterogeneous dose distribution within the planning target volumes (PTVs). In this scenario, the location of cold spots within the target becomes a crucial issue and has to be related to the distribution of the clonogenic cell density (CCD). The main aim of this work is to provide the means to calculate the optimal prescription dose in a relative inhomogeneous dose distribution. To achieve this, the prescription dose has to be assigned to obtain the same tumor control probability (TCP) as the ideal homogeneous distribution, taking into account different CCDs in different PTVs (i.e. visible and subclinical regions). An adapted formulation of the linear-quadratic model, within the F-factor formalism, has been derived to preserve a chosen TCP value for the whole target volume. The F-factor has been investigated to show its potential applications in clinical practice.

Conformity of LINAC-based stereotactic radiosurgery using dynamic conformal arcs and micro-multileaf collimator

PURPOSE

To assess the conformity of dynamic conformal arc linear accelerator-based stereotactic radiosurgery and to describe a standardized method of isodose surface (IDS) selection.

METHODS AND MATERIALS

In 174 targets, the conformity index (CI) at the prescription IDS used for treatment was calculated as CI = (PIV/PVTV)/(PVTV/TV), where TV is the target volume, PIV (prescription isodose volume) is the total volume encompassed by the prescription IDS, and PVTV is the TV encompassed by the IDS. In addition, a "standardized" prescription IDS (sIDS) was chosen according to the following criteria: 95% of the TV was encompassed by the PIV and 99% of TV was covered by 95% of the prescription dose. The CIs at the sIDS were also calculated.

RESULTS

The median CI at the prescription IDS and sIDS was 1.63 and 1.47, respectively (p < 0.001). In 132 of 174 cases, the volume of normal tissue in the PIV was reduced by the prescription to the sIDS compared with the prescription IDS, in 20 cases it remained unchanged, and in 22 cases it was increased.

CONCLUSION

The CIs obtained with linear accelerator-based stereotactic radiosurgery are comparable to those previously reported for gamma knife stereotactic radiosurgery. Using a uniform method to select the sIDS, adequate target coverage was usually achievable with prescription to an IDS greater than that chosen by the treating physician (prescription IDS), providing sparing of normal tissue. Thus, the sIDS might aid physicians in identifying a prescription IDS that balances coverage and conformity.

Optimization of tumour control probability in hypoxic tumours by radiation dose redistribution: a modelling study

Tumour hypoxia is a known cause of clinical resistance to radiation therapy. The purpose of this work was to model the effects on tumour control probability (TCP) of selectively boosting the dose to hypoxic regions in a tumour, while keeping the mean tumour dose constant. A tumour model with a continuous oxygen distribution, incorporating pO(2) histograms published for head and neck patients, was developed. Temporal and spatial variations in the oxygen distribution, non-uniform cell density and cell proliferation during treatment were included in the tumour modelling. Non-uniform dose prescriptions were made based on a segmentation of the tumours into four compartments. The main findings were: (1) Dose redistribution considerably improved TCP for all tumours. (2) The effect on TCP depended on the degree of reoxygenation during treatment, with a maximum relative increase in TCP for tumours with poor or no reoxygenation. (3) Acute hypoxia reduced TCP moderately, while underdosing chronic hypoxic cells gave large reductions in TCP. (4) Restricted dose redistribution still gave a substantial increase in TCP as compared to uniform dose boosts. In conclusion, redistributing dose according to tumour oxygenation status might increase TCP when the tumour response to radiotherapy is limited by chronic hypoxia. This could potentially improve treatment outcome in a subpopulation of patients who respond poorly to conventional radiotherapy.

Biological optimization of heterogeneous dose distributions in systemic radiotherapy

The standard computational method developed for internal radiation dosimetry is the MIRD (medical internal radiation dose) formalism, based on the assumption that tumor control is given by uniform dose and activity distributions. In modern systemic radiotherapy, however, the need for full 3D dose calculations that take into account the heterogeneous distribution of activity in the patient is now understood. When information on nonuniform distribution of activity becomes available from functional imaging, a more patient specific 3D dosimetry can be performed. Application of radiobiological models can be useful to correlate the calculated heterogeneous dose distributions to the current knowledge on tumor control probability of a homogeneous dose distribution. Our contribution to this field is the introduction of a parameter, the F factor, already used by our group in studying external beam radiotherapy treatments. This parameter allows one to write a simplified expression for tumor control probability (TCP) based on the standard linear quadratic (LQ) model and Poisson statistics. The LQ model was extended to include different treatment regimes involving source decay, incorporating the repair "micro" of sublethal radiation damage, the relative biological effectiveness and the effective "waste" of dose delivered when repopulation occurs. The sensitivity of the F factor against radiobiological parameters (alpha, beta, micro) and the influence of the dose volume distribution was evaluated. Some test examples for 131I and 90Y labeled pharmaceuticals are described to further explain the properties of the F factor and its potential applications. To demonstrate dosimetric feasibility and advantages of the proposed F factor formalism in systemic radiotherapy, we have performed a retrospective planning study on selected patient case. F factor formalism helps to assess the total activity to be administered to the patient taking into account the heterogeneity in activity uptake and dose distribution, giving the same TCP of a homogeneous prescribed dose distribution. Animal studies and collection of standardized clinical data are needed to ascertain the effects of nonuniform dose distributions and to better assess the radiobiological input parameters of the model based on LQ model.

Large deformation three-dimensional image registration in image-guided radiation therapy

In this paper, we present and validate a framework, based on deformable image registration, for automatic processing of serial three-dimensional CT images used in image-guided radiation therapy. A major assumption in deformable image registration has been that, if two images are being registered, every point of one image corresponds appropriately to some point in the other. For intra-treatment images of the prostate, however, this assumption is violated by the variable presence of bowel gas. The framework presented here explicitly extends previous deformable image registration algorithms to accommodate such regions in the image for which no correspondence exists. We show how to use our registration technique as a tool for organ segmentation, and present a statistical analysis of this segmentation method, validating it by comparison with multiple human raters. We also show how the deformable registration technique can be used to determine the dosimetric effect of a given plan in the presence of non-rigid tissue motion. In addition to dose accumulation, we describe a method for estimating the biological effects of tissue motion using a linear-quadratic model. This work is described in the context of a prostate treatment protocol, but it is of general applicability.

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.

A study of the effects of internal organ motion on dose escalation in conformal prostate treatments

BACKGROUND AND PURPOSE

To assess the effect of internal organ motion on the dose distributions and biological indices for the target and non-target organs for three different conformal prostate treatment techniques.

MATERIALS AND METHODS

We examined three types of treatment plans in 20 patients: (1) a six field plan, with a prescribed dose of 75.6 Gy; (2) the same six field plan to 72 Gy followed by a boost to 81 Gy; and (3) a five field plan with intensity modulated beams delivering 81 Gy. Treatment plans were designed using an initial CT data set (planning) and applied to three subsequent CT scans (treatment). The treatment CT contours were used to represent patient specific organ displacement; in addition, the dose distribution was convolved with a Gaussian distribution to model random setup error. Dose-volume histograms were calculated using an organ deformation model in which the movement between scans of individual points interior to the organs was tracked and the dose accumulated. The tumor control probability (TCP) for the prostate and proximal half of seminal vesicles (clinical target volume, CTV), normal tissue complication probability (NTCP) for the rectum and the percent volume of bladder wall receiving at least 75 Gy were calculated.

RESULTS

The patient averaged increase in the planned TCP between plan types 2 and 1 and types 3 and 1 was 9.8% (range 4.9-12.5%) for both, whereas the corresponding increases in treatment TCP were 9.0% (1.3-16%) and 8.1% (-1.3-13.8%). In all patients, plans 2 and 3 (81 Gy) exhibited equal or higher treatment TCP than plan 1 (75.6 Gy). The maximum treatment NTCP for rectum never exceeded the planning constraint and percent volume of bladder wall receiving at least 75 Gy was similar in the planning and treatment scans for all three plans.

CONCLUSION

For plans that deliver a uniform prescribed dose to the planning target volume (PTV) (plan 1), current margins are adequate. In plans that further escalate the dose to part of the PTV (plans 2 and 3), in a fraction of the cases the CTV dose increase is less than planned, yet in all cases the TCP values are higher relative to the uniform dose PTV (plan 1). Doses to critical organs remain within the planning criteria.

Should single or distributed parameters be used to explain the steepness of tumour control probability curves?

Linear quadratic (LQ) modelling allows easy comparison of different fractionation schedules in radiotherapy. However, estimating the radiation effect of a single fractionated treatment introduces many questions with respect to the parameters to be used in the modelling process. Several studies have used tumour control probability (TCP) curves in order to derive the values for the LQ parameters that may be used further for the analysis and ranking of treatment plans. Unfortunately, little attention has been paid to the biological relevance of these derived parameters, either for the initial number of cells or their intrinsic radiosensitivity, or both. This paper investigates the relationship between single values for the TCP parameters and the resulting dose-response curve. The results of this modelling study show how clinical observations for the position and steepness of the TCP curve can be explained only by the choice of extreme values for the parameters, if they are single values. These extreme values are in contradiction with experimental observations. This contradiction suggests that single values for the parameters are not likely to explain reasonably the clinical observations and that some distributions of input parameters should be taken into consideration.

Tumor control probability for selective boosting of hypoxic subvolumes, including the effect of reoxygenation

PURPOSE

To study the effect on tumor control probability of selectively boosting the dose to hypoxic subvolumes.

METHODS AND MATERIALS

A Monte Carlo model was developed that separates the tumor into two compartments, one of which receives a primary dose, and one of which receives a higher boost dose. During radiation delivery, each compartment consists of three clonogen subpopulations: those that are well oxygenated, those that are temporarily hypoxic (geometrically transient hypoxia), and those that are permanently hypoxic (geometrically stable hypoxia). The spatial location of temporary hypoxia within the tumor volume varies over time, whereas, the spatial location of permanent hypoxia does not. The effect of reoxygenation was included. Clonogen proliferation was not included in the model.

RESULTS

A modest boost dose (120%-150% of the primary dose) increases tumor control probability to that found in the absence of permanent hypoxia. The entire hypoxic subvolume need not be included to obtain a significant benefit. However, only tumors with a geometrically stable hypoxic volume will have an improved control rate.

CONCLUSIONS

Tumors with an identifiable geometrically stable hypoxic volume will have an improved control rate if the dose to the hypoxic volume is escalated. Further work is required to determine the spatiotemporal evolution of the hypoxic volumes before and during the course of radiotherapy.

Is uniform target dose possible in IMRT plans in the head and neck?

PURPOSE

Various published reports involving intensity-modulated radiotherapy (IMRT) plans developed using automated optimization (inverse planning) have demonstrated highly conformal plans. These reported conformal IMRT plans involve significant target dose inhomogeneity, including both overdosage and underdosage within the target volume. In this study, we demonstrate the development of optimized beamlet IMRT plans that satisfy rigorous dose homogeneity requirements for all target volumes (e.g., +/-5%), while also sparing the parotids and other normal structures.

METHODS AND MATERIALS

The treatment plans of 15 patients with oropharyngeal cancer who were previously treated with forward-planned multisegmental IMRT were planned again using an automated optimization system developed in-house. The optimization system allows for variable sized beamlets computed using a three-dimensional convolution/superposition dose calculation and flexible cost functions derived from combinations of clinically relevant factors (costlets) that can include dose, dose-volume, and biologic model-based costlets. The current study compared optimized IMRT plans designed to treat the various planning target volumes to doses of 66, 60, and 54 Gy with varying target dose homogeneity while using a flexible optimization cost function to minimize the dose to the parotids, spinal cord, oral cavity, brainstem, submandibular nodes, and other structures.

RESULTS

In all cases, target dose uniformity was achieved through steeply varying dose-based costs. Differences in clinical plan evaluation metrics were evaluated for individual cases (eight different target homogeneity costlets), and for the entire cohort of plans. Highly conformal plans were achieved, with significant sparing of both the contralateral and ipsilateral parotid glands. As the homogeneity of the target dose distributions was allowed to decrease, increased sparing of the parotids (and other normal tissues) may be achieved. However, it was shown that relatively few patients would benefit from the use of increased target inhomogeneity, because the range of improvement in the parotid dose is relatively limited. Hot spots in the target volumes are shown to be unnecessary and do not assist in normal tissue sparing.

CONCLUSION

Sparing of both parotids in patients receiving bilateral neck radiation can be achieved without compromising strict target dose homogeneity criteria. The geometry of the normal tissue and target anatomy are shown to be the major factor necessary to predict the parotid sparing that will be possible for any particular case.

Generalization of a model of tissue response to radiation based on the idea of functional subunits and binomial statistics

This work investigates the existing biological models describing the response of tumours and normal tissues to radiation, with the purpose of developing a general biological model of the response of tissue to radiation. Two different types of normal tissue behaviour have been postulated with respect to its response to radiation, namely critical element and critical volume behaviour. Based on the idea that an organ is composed of functional subunits, models have been developed describing these behaviours. However, these models describe the response of an individual, a particular patient or experimental animal, while the clinically or experimentally observed quantity is the population response. There is a need to extend the models to address the population response, based on the ideas we have about the individual response. We have attempted here to summarize and unify the existing individual models. Finally, the population models are investigated by fitting to pseudoexperimental sets of data and comparing them with each other in terms of goodness-of-fit and in terms of their power to recover the values of the population parameters.

Partially wedged beams improve radiotherapy treatment of urinary bladder cancer

BACKGROUND AND PURPOSE

Partially wedged beams (PWBs) having wedge in one part of the field only, can be shaped using dynamic jaw intensity modulation. The possible clinical benefit of PWBs was tested in treatment plans for muscle-infiltrating bladder cancer.

MATERIAL AND METHODS

Three-dimensional treatment plans for 25 bladder cancer patients were analyzed. The originally prescribed standard conformal four-field box technique, which includes the use of lateral ordinary wedge beams, was compared to a modified conformal treatment using customized lateral PWBs. In these modified treatment plans, only the anterior parts of the two lateral beams had a wedge. To analyze the potential clinical benefit of treatment with PWBs, treatment plans were scored and compared using both physical parameters and biological dose response models. One tumour control probability model and two normal tissue complication probability (NTCP) models were applied. Different parameters for normal tissue radiation tolerance presented in the literature were used.

RESULTS

By PWBs the dose homogeneity throughout the target volume was improved for all patients, reducing the average relative standard deviation of the target dose distribution from 2.3 to 1.8%. A consistent reduction in the maximum doses to surrounding normal tissue volumes was also found. The most notable improvement was demonstrated in the rectum where the volume receiving more than the prescribed tumour dose was halved. Treatment with PWBs would permit a target dose escalation of 2-6 Gy in several of the patients analyzed, without increasing the overall risk for complications. The number of patients suitable for dose escalation ranged from 3 to 15, depending on whether support from all or only one of the five applied NTCP model/parameter combinations were required in each case to recommend dose escalation.

CONCLUSION

PWBs represent a simple dose conformation tool that may allow radiation dose escalation in the treatment of muscle-infiltrating urinary bladder tumours.

The delta-TCP concept: a clinically useful measure of tumor control probability

PURPOSE

The aim of this article is to provide a quantitative tool to evaluate the influence of the different dose regions in a non-uniformly irradiated tumour upon the probability of controlling that tumor.

METHODS AND MATERIALS

First, a method to generate a distribution of the probability of controlling the cells in a voxel (VCP) is explored and found not to be useful. Second, we introduce the concept of delta-TCP, which represents the gain or loss in the overall TCP as a result of each particular bin in a DVH not receiving the prescribed dose (the same concept is applicable to dose cubes or to a fraction of the bin). The delta-TCP method presented here is based on the Poisson TCP model, but any other model could also be used. Third, using this tool, with parameters appropriate to Stage C prostate tumors, the consequences of "cold" and "hot" dose regions have been explored.

RESULTS

We show that TCP is affected by the minimum dose, even if it is delivered to a very small volume (20% dose deficit to 5% of the volume makes the TCP decrease by 18%), and that a hot region may be "wasted" unless the boost is to the bulk of the volume. An example of the application of the delta-TCP concept to a prostate radiotherapy plan is also given.

CONCLUSION

The delta-TCP distribution adds more objective information to the original DVH by enabling the clinician or planner to directly evaluate the effects of a non-uniform dose distribution on local control.

Variations of tumor control and rectum complication probabilities due to random set-up errors during conformal radiation therapy of prostate cancer

BACKGROUND AND PURPOSE

The effect of random set-up errors on tumor control probability (TCP) and rectum complication probability (NTCP) on 3D conformal treatment planning of prostate cancer has been investigated by applying the convolution method originally proposed by Leong (Leong, J. Implementation of random positioning error in computerized radiation treatment planning systems as a result of fractionation. Phys. Med. Biol. 32: 327-334, 1987).

MATERIALS AND METHODS

The combined influence of the standard deviation of the random shifts probability distribution (sigma) of the dose and of the Beam's-eye-view margin (M) between the clinical target volume (CTV) and the edge of the blocks have been investigated in two patients.

RESULTS AND CONCLUSIONS

Random set-up error has been found to decrease TCP (for a typical 70 Gy CTV mean dose) by up to 6% for a 1 cm margin (sigma = 7 mm). When M is equal to or larger than 1.5 cm, no relevant effects on TCP are obtained. Maximum acceptable TCP values (corresponding to a rectum NTCP equal to 5%) have been derived and the dependence on sigma and M has been investigated.

Late GI and GU complications in the treatment of prostate cancer

PURPOSE

To assess the factors that predict late GI and GU morbidity in radiation treatment of the prostate.

METHODS AND MATERIALS

Seven hundred twelve consecutive prostate cancer patients treated at this institution between 1986 and 1994 (inclusive) with conformal or conventional techniques were included in the analysis. Patients had at least 3 months follow-up and received at least 65 Gy. Late GI Grade 3 morbidity was rectal bleeding (requiring three or more procedures) or proctitis. Late Grade 3 GU morbidity was cystitis or stricture. Multivariate analysis (MVA) was used to assess factors related to the complication-free survival. The factors assessed were age, occurrence of side effects > or = Grade 2 during treatment, irradiated volume parameters (use of pelvic fields, treatment of seminal vesicles to full dose or 57 Gy, and use of additional rectal shielding), dose, comorbidities, and other treatments (hormonal manipulation, TURP).

RESULTS

Acute GI and GU side effects (Grade 2 or higher) were noted in 246 and 201 patients, respectively; 67 of these patients exhibited both. GI side effects were not correlated with GU side effects acutely. Late and acute morbidities were correlated (both GI and GU). Fifteen of the 712 patients expressed Grade 3 or 4 GI injuries 3 to 32 months after the end of treatment, with a mean of 14.3 months. One hundred fifteen patients expressed Grade 2 or higher GI morbidity (mean: 13.7 months). The 43 Grade 2 or higher GU morbidities occurred significantly later (mean: 22.7 months). Central axis dose was the only independent variable significantly related to the incidence of late GI morbidity on MVA. No treatment volume parameters were significant for Grade 3. The following parameters were significantly related (by MVA) to Grade 2 GI morbidity: central axis dose, use of the increased rectal shielding, androgen deprivation therapy starting before RT. Acute and late GI morbidities were highly correlated. History of diabetes, treatment of pelvic nodes, and age less than 60 years were significantly related to acute GI side effects. The parameters significantly related to late Grade 2 or higher GU morbidity were central axis dose, androgen deprivation therapy (Zoladex or Lupron) prior to radiation therapy (RT), history of obstructive symptoms, and acute GU side effects. There were too few late Grade 3 GU morbidities to perform multivariate analysis. Acute GU side effects were highly correlated with late GU injury. The following were correlated with acute GU side effects: history of diabetes (+), treatment with conformal fields (-), TURP before RT (-), presentation with urinary obstructive symptoms.

CONCLUSION

Both late GI and GU morbidity demonstrate a dose dependence, but only the volume dependence observed is a reduction in late Grade 2-4 GI morbidity by increasing the rectal shielding in the lateral fields for the final 10 Gy. Moreover, both late GI and GU morbidity was increased in patients treated with hormone manipulation prior to RT. GI and GU injuries were correlated with their corresponding acute side effects. GI and GU complications must not be combined for analysis to determine the factors related to their occurrence.

Adaptive radiation therapy

Adaptive radiation therapy is a closed-loop radiation treatment process where the treatment plan can be modified using a systematic feedback of measurements. Adaptive radiation therapy intends to improve radiation treatment by systematically monitoring treatment variations and incorporating them to re-optimize the treatment plan early on during the course of treatment. In this process, field margin and treatment dose can be routinely customized to each individual patient to achieve a safe dose escalation.

Characterization of dose distribution in radiation therapy plans

As a method of considering only significant radiation doses to different tissues, the ICRU Report 50 recommends taking the dose given to a significant tissue volume (minimum diameter greater then 15 mm) instead of choosing a single, potentially insignificant, voxel value. In order to find this significant volume, we have adapted an emission imaging analysis method to radiation therapy planning. The resulting method finds and characterizes the dose distribution in the volumes of interest in a way that includes spatial arrangement. The data can be used to signal significant hot or cold volumes in the dose plan and to score the plans based on significant dose to the tissues.

Study on the reference dose level in radiotherapy treatment planning

PURPOSE

The reference dose level of the dose distribution in the tumor volume is studied.

METHODS AND MATERIALS

The study is performed using a formula based on the Linear Quadratic (LQ) model. The calculated reference dose level to which the prescribed dose must be referred, for the eradication of a homogeneous tumor, is investigated by varying the dose distribution, that is, the dose volume histogram shape, its range, the prescribed total dose, the fraction size and the linear quadratic model parameters, alpha and beta.

RESULTS

For all the simulated dose volume histograms the calculated reference dose level is lower than the mean dose level, depending on the range of dose variation and the considered tumor sensitivity. When the dose nonuniformity is not too great the reference dose level is very near to the mean dose level; when the inhomogeneity of dose distribution is high the reference level is clearly lower than the mean level but not necessarily equal to the minimum level in the tumor. For the dose volume histograms derived from the actual dose distributions obtained from a two tangential beams technique, a four beams technique and a moving beam technique, the reference levels are calculated and compared with the ICRU 29 reference point dose level. In two cases the reference levels are lower than the level at the ICRU 29 reference point. In the case of the four beams technique, the two levels are equal.

CONCLUSION

These theoretical results show the possibility of administering the prescribed dose to a dose level higher than the minimum in the tumor, with the same value of Tumor Control Probability (TCP) as the one corresponding to a uniform tumor irradiation. The application of the proposed study can offer a general support to the choice of the reference dose level, based on the actual dose distribution in the tumor volume.

Rate of relapse following treatment for localized prostate cancer: a critical analysis of retrospective reports

PURPOSE

Controversy exists over the optimal treatment for patients with clinically localized prostate cancer. Almost all of the treatment results are from non-randomized trials and interseries comparison is difficult since the apparent success of a treatment, as judged by the actuarial freedom from relapse and survival data, depends on patient selection criteria and post-treatment evaluation, in addition to the efficacy of the therapeutic intervention. In this report the calculation of a hazard function is used to estimate and compare the rate of relapse for the different treatments.

METHODS AND MATERIALS

Clinical reports from major surgery and radiation oncology treatment institutions were analyzed. The actuarial recurrence data were used to calculate the annual rate of recurrence within each series.

RESULTS

For all but the lowest volume tumors, patients continue to be at risk of relapse for as long as these series have been followed. Despite the heterogeneity of patient populations, the recurrence rates by stage are similar for patients treated with surgery or irradiation. This result is consistent with pathologic data from prostatectomy specimens which indicate that for lesions > 12 cm3 (approx. 3 cm in diameter) there is high likelihood of extraprostatic disease.

CONCLUSION

Treatment outcome for patients with localized prostate cancer may be more dependent on the inherent tumor biology than the particular type of treatment. Accordingly, the expectation and recommendation of a treatment must take into consideration the continued risk of relapse with either radiation therapy or surgery. There are, as yet, insufficient data regarding the impact of screening and earlier diagnosis on the curability of patients with localized prostate cancer.

Optimization of 3D radiation therapy with both physical and biological end points and constraints

A new optimization model is described and its clinical usefulness is demonstrated. The optimization technique was developed to allow computer optimization of 3-dimensional radiation therapy plans with biological models of tumor and normal tissue response to radiation as well as with scores based on physical dose. The emphasis was placed on the optimization model, which should describe, as closely as possible, the goal of the radiation treatment, which is eradication of the tumor while sparing normal tissues. Since the statement of the goals may vary from case to case, a technique that allows a variety of objective functions and types of constraints was developed. The optimization algorithm is capable of handling nonlinear and even discrete score (objective) functions and constraints and effectively explores the vast space of feasible solutions in a relatively short time (minutes of MicroVax 3200 CPU time). An example of computer optimization of radiation therapy of a chordoma of the sphenoid bone using x-ray and proton beams is shown and compared with the best plans achieved by an experienced planner. Directions for future development of the algorithm, allowing optimization of beam orientation, are presented.

Quantification of doses to mediastinal lymph nodes in Hodgkin's disease

Hodgkin's disease is highly curable today. Radiotherapy (RT) is the treatment of choice in the early stages. A mantle field is often used in the RT of Hodgkin's disease, and the technique and dosimetry are quite complex. We used computerized tomography (CT)-based dosimetry to determine doses delivered to different mediastinal nodes with the commonly used technique in Hodgkin's disease that was originally described by Kaplan. We used dose-volume histograms to determine doses to various groups of nodes in nine patients. Significant inhomogeneity (30%, 30%, 35%, 35%, 30%, 40%, 35%, 35%, and 30% in the nine patients) in dose distribution was found within the mediastinum. With the advent of 3-dimensional CT-based treatment planning, we are able to quantify such inhomogeneities. The question arises whether a homogeneous, lesser dose can achieve equal results. Average doses and "effective doses" were also calculated. The "effective doses" in eight patients (for a prescribed dose of 44 Gy) with a midline posterior spinal cord block added at 20 Gy were 37.3 Gy, 34.3 Gy, 36.0 Gy, 38.4 Gy, 35.8 Gy, 38.1 Gy, 36.7 Gy, and 36.7 Gy, respectively. A homogeneous dose equivalent to effective dose may achieve the same control as an inhomogeneous dose delivery. Prospective 3-D dosimetric studies are required to confirm this concept.

A dose response analysis of injury to cranial nerves and/or nuclei following proton beam radiation therapy

The low tolerance of the central nervous system (CNS) limits the radiation dose which can be delivered in the treatment of many patients with brain and head and neck tumors. Although there are many reports concerning the tolerance of the CNS, few have examined individual substructures of the brain and fewer still have had detailed dose information. This study has both. A three dimensional planning system was used to develop the combined proton beam/photon beam treatments for 27 patients with skull-base tumors. The cranial nerves and their related nuclei were delineated on the planning CT scans and the radiation dose to each was determined from three dimensional dose distributions. In the 594 CNS structures (22 structures/patient in 27 patients), there have been 17 structures (in 5 patients) with clinically manifest radiation injury, after a mean follow-up time of 74 months (range 40-110 months). From statistical analyses, dose is found to be a significant predictor of injury. Using logistic regression analysis, we find that, for each cranial nerve, at 60 Cobalt Gray Equivalent (CGE) the complication rate is 1% (0.5-3% with 95% confidence) and that the 5% complication rate occurs at 70 CGE (64-81 CGE with 95% confidence). The slope of the dose response curve (at 50%) is 3.2 (2.2-5.4 with 95% confidence). No significant relationship between dose and latency period for nerve injury was found.

Beams eye view-based photon radiotherapy I

Geographic miss, dosimetric miss (underdosing), and proximity of the tumor to sensitive normal tissues are some of the causes of inadequate radiation dose delivery; this is one of many causes of failure after radiotherapy. In the past decade, computerized tomography (CT)-based treatment planning has helped to overcome some of these problems. Beam's eye view (BEV)-based radiotherapy planning is an improvement over CT-based treatment planning that may further increase the therapeutic ratio. Since January 1988, we have treated 198 patients with BEV-based photon radiotherapy. About 40% of our patients treated with radical radiotherapy undergo BEV-based treatment, and about 70% of patients who undergo planning CT in the treatment position receive BEV-based radiotherapy. Our findings are as follows: (a) routine use of BEV-based RT (BEVRT) is possible in a busy radiation oncology department; (b) BEVRT improves geometric coverage of tumors; (c) BEVRT is extremely useful in the design of oblique portals; (d) time commitments for various members of the RT treatment-planning team are reasonable; (e) BEVRT helps individualize RT technique; (f) preliminary data suggest decreased acute toxicity with the use of BEVRT for prostate cancer patients. Whether these advantages will help to improve the outcome (i.e., improve local control and survival) and/or decrease the long-term toxicity is not yet known.

The effect on minimum tumor dose of restricting target-dose inhomogeneity in optimized three-dimensional treatment of lung cancer

An examination was made of the effect upon the minimum tumor dose of a limit placed on the variation of dose across target. If the required level of target dose uniformity is slightly relaxed, a substantial improvement in the minimum tumor dose might appear. It was conjectured that this effect could be seen with treatments optimally planned and evaluated in three-dimensions. A model of advanced carcinoma of the lung treated with a computer controlled accelerator was used to test this hypothesis. A mathematical program for optimizing beam weights was used to determine the largest minimum tumor dose possible. In the six cases tested, a minimum tumor dose of greater than 80 Gy could be delivered if a 20% inhomogeneity limit was accepted. The minimum tumor dose fell to the range 44-64 Gy when the inhomogeneity limit was tightened to 13-17%. The results imply a need to examine the choice of a required level of dose uniformity from the range of values suggested in the 2-dimensional planning literature. If a strict bound-on-dose uniformity is preserved, mechanisms--such as formal optimization--which can reduce target dose inhomogeneity will be valuable.

Improved dose distributions for 3D conformal boost treatments in carcinoma of the nasopharynx

This study was designed to demonstrate the feasibility of 3-dimensional (3D) treatment planning in patients with carcinoma of the nasopharynx, and to explore its potential therapeutic advantage over the traditional 2-dimensional (2D) approach in this disease. Qualitative and quantitative comparisons between the two techniques were made for the boost portion of the treatment (19.8 Gy of a total 70.2 Gy treatment schedule) in 10 previously untreated patients and for the entire treatment in 5 patients with locally recurrent disease. The 2D and 3D plans were compared in each patient using dose-volume histograms (DVH's), tumor control probabilities (TCP's), normal tissue complication probabilities (NTCP's), and a new biologic figure of merit that describes the probability of uncomplicated control. Although there was no attempt to optimize the 3D treatment approach by using this method throughout the total treatment course (rather than for the boost only), it was still found that for each of the endpoints examined the 3D approach resulted in improved plans. An average of 22% of the target volume was underdosed at the 95% isodose level with the 2D plans compared to 7% with the 3D plans. The improved treatment planning by 3D increased the mean dose to the tumor volume by an average of 13% over 2D planning. The dose to normal structures such as the mandible and parotid glands was reduced with the 3D plans while the brain stem and spinal cord remained within tolerance limits. The probability of uncomplicated tumor control was increased by an average of 15% with 3D treatment planning compared to the 2D approach. Our findings demonstrate the potential of 3D planning for improving the treatment of carcinoma of the nasopharynx, but prospective studies are required to define the true clinical advantages of this methodology.

The biological effect of inhomogeneous dose distributions in fractionated radiotherapy

The linear quadratic (LQ) model is applied to an organ receiving a fractionated course of radiotherapy with an inhomogeneous dose distribution. It is shown that the gradient in the extrapolated response dose (ERD) will be steeper than the gradient in the physical dose. This effect will be greatest for an organ with a small alpha/beta ratio treated with large dose fractions. Clinical implications are discussed with an emphasis on radiation myelitis.

Optimisation of conformal radiotherapy dose distributions by simulated annealing

A method of computing the beam profiles for multi-element multiple-beam radiotherapy treatment is presented. This is applicable to the treatment planning of conformal radiotherapy using an isocentric rotation technique. The method begins with the treatment dose prescription and calculates the beam profiles rather than the reverse 'conventional' planning technique. The method is the iterative optimisation method of simulated annealing. It is shown using published 'difficult' clinical treatment planning problems that the high-dose region can be very precisely tailored to the tumour volume even when this has a re-entrant (concave) periphery. Simultaneously the dose can be constrained in other sensitive regions. The mathematics of the method is explained together with this implementation. The analogy of this optimisation with certain reconstruction problems in medical imaging is drawn and by way of experiment dose distributions are presented in image form and beam profiles as sinograms.

The promise of a new technology: knowledge-based systems in radiation oncology and diagnostic radiology

The revolutionary changes in computer capabilities in the last decade, both in software and hardware, have opened new doorways for the uses of computers in radiation oncology and diagnostic radiology. Knowledge-based systems offer the potential to function as aids, consultants and advisors in the differential diagnosis of disease, staging, selection of therapy and treatment management and delivery for cancer patients. These computer-based systems can also provide for the training and teaching of radiotherapy and diagnostic radiology residents, and act as advisors and teachers to the medical physicists, dosimetrists and technicians. Following a brief history of the development of knowledge-based systems, the general capabilities of computer-based physician workstations in a department of radiation oncology are described.

Accuracy in radiation field alignment in head and neck cancer: a prospective study

A prospective study has been performed to determine the accuracy of radiation field alignment for a group of 22 patients with tumors in the head and neck. The accuracy was assessed by an analysis of 138 megavolt portal films in comparison to 55 simulation films. The distance (at the patient midplane) between corresponding points at the field edges on verification film and simulation film appeared to be 5 mm on the average and the standard deviation 5 mm. The analysis was extended by translational and rotational matching of the fields in order to separate each error in a translation error of the field with respect to the patient and an error in field size or shape. Translation errors appear to be somewhat larger than field size or shape errors. From an analysis of a series of megavolt films taken every third radiotherapy session, it was concluded that treatment-to-treatment variations are as large as the errors due to the transition from simulation to treatment situation. Further analysis showed that variation of the patient's position within the cast is clearly one of the error sources.

Some implications of the Linear Quadratic model for tumor control probability

To define an optimal radiation therapy strategy, the dependence of the probability of cure and of significant complications, on the parameters controlled by the radiotherapist must be determined. The recent success of the Linear Quadratic (LQ) model in constructing isoeffect relations for normal tissue damage and in describing in vitro cell survival curves, indicate that this model could be used to determine this dependence. The problem of tumor control is addressed. Using LQ model parameters obtained from human tumor cell lines, the sigmoid dose-response curves for controlling tumors of fixed size but of several histologies, with a fractionated course of radiotherapy is obtained. Except for squamous cell carcinoma, the calculated average tumor control doses (TCD37 or TCD50) are unrealistically low, but the model can be made more realistic by including inhomogeneities in the spatial dose distribution and heterogeneous tumor cell populations. The slope of the dose response curves are determined and the significance of the "relative slope" parameter rho as a measure of the number of cells in a tumor's most radioresistant clone is noted. The relation of the model's predictions to qualitative features of the experimental animal data for both tumor control and for normal tissue damage is discussed. Experiments to test the validity of this type of model are suggested.