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

Dosimetric study of uniform scanning proton therapy planning for prostate cancer patients with a metal hip prosthesis, and comparison with volumetric-modulated arc therapy

The main purposes of this study were to 1) investigate the dosimetric quality of uniform scanning proton therapy planning (USPT) for prostate cancer patients with a metal hip prosthesis, and 2) compare the dosimetric results of USPT with that of volumetric‐modulated arc therapy (VMAT). Proton plans for prostate cancer (four cases) were generated in XiO treatment planning system (TPS). The beam arrangement in each proton plan consisted of three fields (two oblique fields and one lateral or slightly angled field), and the proton beams passing through a metal hip prosthesis was avoided. Dose calculations in proton plans were performed using the pencil beam algorithm. From each proton plan, planning target volume (PTV) coverage value (i.e., relative volume of the PTV receiving the prescription dose of 79.2 CGE) was recorded. The VMAT prostate planning was done using two arcs in the Eclipse TPS utilizing 6 MV X‐rays, and beam entrance through metallic hip prosthesis was avoided. Dose computation in the VMAT plans was done using anisotropic analytical algorithm, and calculated VMAT plans were then normalized such that the PTV coverage in the VMAT plan was the same as in the proton plan of the corresponding case. The dose‐volume histograms of calculated treatment plans were used to evaluate the dosimetric quality of USPT and VMAT. In comparison to the proton plans, on average, the maximum and mean doses to the PTV were higher in the VMAT plans by 1.4% and 0.5%, respectively, whereas the minimum PTV dose was lower in the VMAT plans by 3.4%. The proton plans had lower (or better) average homogeneity index (HI) of 0.03 compared to the one for VMAT (HI = 0.04). The relative rectal volume exposed to radiation was lower in the proton plan, with an average absolute difference ranging from 0.1% to 32.6%. In contrast, using proton planning, the relative bladder volume exposed to radiation was higher at high‐dose region with an average absolute difference ranging from 0.4% to 0.8%, and lower at low‐ and medium‐dose regions with an average absolute difference ranging from 2.7% to 10.1%. The average mean dose to the rectum and bladder was lower in the proton plans by 45.1% and 22.0%, respectively, whereas the mean dose to femoral head was lower in VMAT plans by an average difference of 79.6%. In comparison to the VMAT, the proton planning produced lower equivalent uniform dose (EUD) for the rectum (43.7 CGE vs. 51.4 Gy) and higher EUD for the femoral head (16.7 CGE vs. 9.5 Gy), whereas both the VMAT and proton planning produced comparable EUDs for the prostate tumor (76.2 CGE vs. 76.8 Gy) and bladder (50.3 CGE vs. 51.1 Gy). The results presented in this study show that the combination of lateral and oblique fields in USPT planning could potentially provide dosimetric advantage over the VMAT for prostate cancer involving a metallic hip prosthesis.

PACS number: 87.55.D‐, 87.55.ne, 87.55.dk

Dosimetric effect of tissue heterogeneity for (125)I prostate implants

AIM

To use Monte Carlo (MC) together with voxel phantoms to analyze the tissue heterogeneity effect in the dose distributions and equivalent uniform dose (EUD) for (125)I prostate implants.

BACKGROUND

Dose distribution calculations in low dose-rate brachytherapy are based on the dose deposition around a single source in a water phantom. This formalism does not take into account tissue heterogeneities, interseed attenuation, or finite patient dimensions effects. Tissue composition is especially important due to the photoelectric effect.

MATERIALS AND METHODS

The computed tomographies (CT) of two patients with prostate cancer were used to create voxel phantoms for the MC simulations. An elemental composition and density were assigned to each structure. Densities of the prostate, vesicles, rectum and bladder were determined through the CT electronic densities of 100 patients. The same simulations were performed considering the same phantom as pure water. Results were compared via dose-volume histograms and EUD for the prostate and rectum.

RESULTS

The mean absorbed doses presented deviations of 3.3-4.0% for the prostate and of 2.3-4.9% for the rectum, when comparing calculations in water with calculations in the heterogeneous phantom. In the calculations in water, the prostate D 90 was overestimated by 2.8-3.9% and the rectum D 0.1cc resulted in dose differences of 6-8%. The EUD resulted in an overestimation of 3.5-3.7% for the prostate and of 7.7-8.3% for the rectum.

CONCLUSIONS

The deposited dose was consistently overestimated for the simulation in water. In order to increase the accuracy in the determination of dose distributions, especially around the rectum, the introduction of the model-based algorithms is recommended.

Volumetric modulated arc therapy versus step-and-shoot intensity modulated radiation therapy in the treatment of large nerve perineural spread to the skull base: a comparative dosimetric planning study

Introduction

Cutaneous squamous cell carcinoma with large nerve perineural (LNPN) infiltration of the base of skull is a radiotherapeutic challenge given the complex target volumes to nearby organs at risk (OAR). A comparative planning study was undertaken to evaluate dosimetric differences between volumetric modulated arc therapy (VMAT) versus intensity modulated radiation therapy (IMRT) in the treatment of LNPN.

Methods

Five consecutive patients previously treated with IMRT for LNPN were selected. VMAT plans were generated for each case using the same planning target volumes (PTV), dose prescriptions and OAR constraints as IMRT. Comparative parameters used to assess target volume coverage, conformity and homogeneity included V95 of the PTV (volume encompassed by the 95% isodose), conformity index (CI) and homogeneity index (HI). In addition, OAR maximum point doses, V20, V30, non-target tissue (NTT) point max doses, NTT volume above reference dose, monitor units (MU) were compared.

Results

IMRT and VMAT plans generated were comparable for CI (P = 0.12) and HI (P = 0.89). VMAT plans achieved better V95 (P = < 0.001) and reduced V20 and V30 by 652 cubic centimetres (cc) (28.5%) and 425.7 cc (29.1%), respectively. VMAT increased MU delivered by 18% without a corresponding increase in NTT dose.

Conclusion

Compared with IMRT plans for LNPN, VMAT achieved comparable HI and CI.

Redefining relative biological effectiveness in the context of the EQDX formalism: implications for alpha-particle emitter therapy

Alpha-particle radiopharmaceutical therapy (αRPT) is currently enjoying increasing attention as a viable alternative to chemotherapy for targeting of disseminated micrometastatic disease. In theory, αRPT can be personalized through pre-therapeutic imaging and dosimetry. However, in practice, given the particularities of α-particle emissions, a dosimetric methodology that accurately predicts the thresholds for organ toxicity has not been reported. This is in part due to the fact that the biological effects caused by α-particle radiation differ markedly from the effects caused by traditional external beam (photon or electron) radiation or β-particle emitting radiopharmaceuticals. The concept of relative biological effectiveness (RBE) is used to quantify the ratio of absorbed doses required to achieve a given biological response with alpha particles versus a reference radiation (typically a beta emitter or external beam radiation). However, as conventionally defined, the RBE varies as a function of absorbed dose and therefore a single RBE value is limited in its utility because it cannot be used to predict response over a wide range of absorbed doses. Therefore, efforts are underway to standardize bioeffect modeling for different fractionation schemes and dose rates for both nuclear medicine and external beam radiotherapy. Given the preponderant use of external beams of radiation compared to nuclear medicine in cancer therapy, the more clinically relevant quantity, the 2 Gy equieffective dose, EQD2(α/β), has recently been proposed by the ICRU. In concert with EQD2(α/β), we introduce a new, redefined RBE quantity, named RBE2(α/β), as the ratio of the two linear coefficients that characterize the α particle absorbed dose-response curve and the low-LET megavoltage photon 2 Gy fraction equieffective dose-response curve. The theoretical framework for the proposed new formalism is presented along with its application to experimental data obtained from irradiation of a breast cancer cell line. Radiobiological parameters are obtained using the linear quadratic model to fit cell survival data for MDA-MB-231 human breast cancer cells that were irradiated with either α particles or a single fraction of low-LET (137)Cs γ rays. From these, the linear coefficient for both the biologically effective dose (BED) and the EQD2(α/β) response lines were derived for fractionated irradiation. The standard RBE calculation, using the traditional single fraction reference radiation, gave RBE values that ranged from 2.4 for a surviving fraction of 0.82-6.0 for a surviving fraction of 0.02, while the dose-independent RBE2(4.6) value was 4.5 for all surviving fraction values. Furthermore, bioeffect modeling with RBE2(α/β) and EQD2(α/β) demonstrated the capacity to predict the surviving fraction of cells irradiated with acute and fractionated low-LET radiation, α particles and chronic exponentially decreasing dose rates of low-LET radiation. RBE2(α/β) is independent of absorbed dose for α-particle emitters and it provides a more logical framework for data reporting and conversion to equieffective dose than the conventional dose-dependent definition of RBE. Moreover, it provides a much needed foundation for the ongoing development of an α-particle dosimetry paradigm and will facilitate the use of tolerance dose data available from external beam radiation therapy, thereby helping to develop αRPT as a single modality as well as for combination therapies.

Prostate stereotactic ablative radiation therapy using volumetric modulated arc therapy to dominant intraprostatic lesions

PURPOSE

To investigate boosting dominant intraprostatic lesions (DILs) in the context of stereotactic ablative radiation therapy (SABR) and to examine the impact on tumor control probability (TCP) and normal tissue complication probability (NTCP).

METHODS AND MATERIALS

Ten prostate datasets were selected. DILs were defined using T2-weighted, dynamic contrast-enhanced and diffusion-weighted magnetic resonance imaging. Four plans were produced for each dataset: (1) no boost to DILs; (2) boost to DILs, no seminal vesicles in prescription; (3) boost to DILs, proximal seminal vesicles (proxSV) prescribed intermediate dose; and (4) boost to DILs, proxSV prescribed higher dose. The prostate planning target volume (PTV) prescription was 42.7 Gy in 7 fractions. DILs were initially prescribed 115% of the PTV(Prostate) prescription, and PTV(DIL) prescriptions were increased in 5% increments until organ-at-risk constraints were reached. TCP and NTCP calculations used the LQ-Poisson Marsden, and Lyman-Kutcher-Burman models respectively.

RESULTS

When treating the prostate alone, the median PTV(DIL) prescription was 125% (range: 110%-140%) of the PTV(Prostate) prescription. Median PTV(DIL) D50% was 55.1 Gy (range: 49.6-62.6 Gy). The same PTV(DIL) prescriptions and similar PTV(DIL) median doses were possible when including the proxSV within the prescription. TCP depended on prostate α/β ratio and was highest with an α/β ratio = 1.5 Gy, where the additional TCP benefit of DIL boosting was least. Rectal NTCP increased with DIL boosting and was considered unacceptably high in 5 cases, which, when replanned with an emphasis on reducing maximum dose to 0.5 cm(3) of rectum (Dmax(0.5cc)), as well as meeting existing constraints, resulted in considerable rectal NTCP reductions.

CONCLUSIONS

Boosting DILs in the context of SABR is technically feasible but should be approached with caution. If this therapy is adopted, strict rectal constraints are required including Dmax(0.5cc). If the α/β ratio of prostate cancer is 1.5 Gy or less, then high TCP and low NTCP can be achieved by prescribing SABR to the whole prostate, without the need for DIL boosting.

Nonisocentric treatment strategy for breast radiation therapy: a proof of concept study

PURPOSE

To propose a nonisocentric treatment strategy as a special form of station parameter optimized radiation therapy, to improve sparing of critical structures while preserving target coverage in breast radiation therapy.

METHODS AND MATERIALS

To minimize the volume of exposed lung and heart in breast irradiation, we propose a novel nonisocentric treatment scheme by strategically placing nonconverging beams with multiple isocenters. As its name suggests, the central axes of these beams do not intersect at a single isocenter as in conventional breast treatment planning. Rather, the isocenter locations and beam directions are carefully selected, in that each beam is only responsible for a certain subvolume of the target, so as to minimize the volume of irradiated normal tissue. When put together, the beams will provide an adequate coverage of the target and expose only a minimal amount of normal tissue to radiation. We apply the nonisocentric planning technique to 2 previously treated clinical cases (breast and chest wall).

RESULTS

The proposed nonisocentric technique substantially improved sparing of the ipsilateral lung. Compared with conventional isocentric plans using 2 tangential beams, the mean lung dose was reduced by 38% and 50% using the proposed technique, and the volume of the ipsilateral lung receiving ≥ 20 Gy was reduced by a factor of approximately 2 and 3 for the breast and chest wall cases, respectively. The improvement in lung sparing is even greater compared with volumetric modulated arc therapy.

CONCLUSIONS

A nonisocentric implementation of station parameter optimized radiation therapy has been proposed for breast radiation therapy. The new treatment scheme overcomes the limitations of existing approaches and affords a useful tool for conformal breast radiation therapy, especially in cases with extreme chest wall curvature.

Impact of geometric variations on delivered dose in highly focused single vocal cord IMRT

PURPOSE

To investigate the robustness of single vocal cord intensity modulated radiation therapy (IMRT) treatment plans for set-up errors, respiration, and deformation.

MATERIAL AND METHODS

Four-dimensional computed tomography (4D-CT) scans of 10 early glottic carcinoma patients, previously treated with conventional techniques, were used in this simulation study. For each patient a pre-treatment 4D-CT was used for IMRT planning, generating a reference dose distribution. Prescribed PTV dose was 66 Gy. The impact of systematic set-up errors was simulated by applying shifts of ± 2 mm to the planning CT scans, followed by dose re-calculation with original beam segments, MUs, etc. Effects of respiration and deformation were determined utilizing extreme inhale and exhale CT scans, and repeat scans acquired after 22 Gy, 44 Gy, and 66 Gy, respectively. All doses were calculated using Monte Carlo dose simulations.

RESULTS

Considering all investigated geometrical perturbations, reductions in the clinical target volume (CTV) V95%, D98%, D2%, and generalized equivalent uniform dose (gEUD) were limited to 1.2 ± 2.2%, 2.4 ± 2.9%, 0.2 ± 1.8%, and 0.6 ± 1.1 Gy, respectively. The near minimum dose, D98%, was always higher than 89%, and gEUD always remained higher than 66 Gy. Planned contra-lateral (CL) vocal cord DMean, gEUD, and V40 Gy were 38.2 ± 6.0 Gy, 43.4 ± 5.6 Gy, and 42.7 ± 14.9%. With perturbations these values changed by -0.1 ± 4.3 Gy, 0.1 ± 4.0 Gy, and -1.0 ± 9.6%, respectively.

CONCLUSIONS

On average, CTV dose reductions due to geometrical perturbations were very low, and sparing of the CL vocal cord was maintained. In a few observations (6 of 103 simulated situations), the near-minimum CTV-dose was around 90%, requiring attention in deciding on a future clinical protocol.

Improving IMRT delivery efficiency with reweighted L1-minimization for inverse planning

PURPOSE

This study presents an improved technique to further simplify the fluence-map in intensity modulated radiation therapy (IMRT) inverse planning, thereby reducing plan complexity and improving delivery efficiency, while maintaining the plan quality.

METHODS

First-order total-variation (TV) minimization (min.) based on L1-norm has been proposed to reduce the complexity of fluence-map in IMRT by generating sparse fluence-map variations. However, with stronger dose sparing to the critical structures, the inevitable increase in the fluence-map complexity can lead to inefficient dose delivery. Theoretically, L0-min. is the ideal solution for the sparse signal recovery problem, yet practically intractable due to its nonconvexity of the objective function. As an alternative, the authors use the iteratively reweighted L1-min. technique to incorporate the benefits of the L0-norm into the tractability of L1-min. The weight multiplied to each element is inversely related to the magnitude of the corresponding element, which is iteratively updated by the reweighting process. The proposed penalizing process combined with TV min. further improves sparsity in the fluence-map variations, hence ultimately enhancing the delivery efficiency. To validate the proposed method, this work compares three treatment plans obtained from quadratic min. (generally used in clinic IMRT), conventional TV min., and our proposed reweighted TV min. techniques, implemented by a large-scale L1-solver (template for first-order conic solver), for five patient clinical data. Criteria such as conformation number (CN), modulation index (MI), and estimated treatment time are employed to assess the relationship between the plan quality and delivery efficiency.

RESULTS

The proposed method yields simpler fluence-maps than the quadratic and conventional TV based techniques. To attain a given CN and dose sparing to the critical organs for 5 clinical cases, the proposed method reduces the number of segments by 10-15 and 30-35, relative to TV min. and quadratic min. based plans, while MIs decreases by about 20%-30% and 40%-60% over the plans by two existing techniques, respectively. With such conditions, the total treatment time of the plans obtained from our proposed method can be reduced by 12-30 s and 30-80 s mainly due to greatly shorter multileaf collimator (MLC) traveling time in IMRT step-and-shoot delivery.

CONCLUSIONS

The reweighted L1-minimization technique provides a promising solution to simplify the fluence-map variations in IMRT inverse planning. It improves the delivery efficiency by reducing the entire segments and treatment time, while maintaining the plan quality in terms of target conformity and critical structure sparing.

Liver toxicity prediction with stereotactic body radiation therapy: The impact of accounting for fraction size

PURPOSE

Current liver SBRT protocols rely on the calculation of "effective volume" without accounting for the biologic effect of fraction size to estimate the risk of liver toxicity, which subsequently defines tumor prescription doses. This study compared effective volume and liver toxicity predictions with and without correction for fraction size.

METHODS AND MATERIALS

The effective volume was determined for 18 liver SBRT plans with and without biologic normalization using the linear quadratic formula. Lyman-Kutcher-Burman normal tissue complication probability models estimated the risk of liver toxicity. Effective volumes and corresponding toxicity predictions were compared with and without biologic normalization.

RESULTS

Accounting for the biologic difference of larger fraction size reduced the effective volume in all treatment plans compared with the unadjusted effective volume (median effective volume 0.21 vs 0.32). The lower effective volume with biologic normalization substantially reduced the estimated risk of liver toxicity (average risk of toxicity 32% vs 4.5%).

CONCLUSIONS

This study demonstrates that accounting for the biologic effect of fraction size with effective volume significantly decreases predicted hepatic toxicity, which suggests that the risk of liver toxicity may be overestimated in clinical practice if biologic normalization is omitted. The effective volume toxicity model has proven safe in prospective clinical trials, though safe dose escalation with liver SBRT may be feasible.

Analysis of equivalent uniform dose (EUD) and conventional radiation treatment parameters after primary and re-irradiation of malignant glioma

BACKGROUND

Re-irradiation is a reasonable second treatment option for patients with recurrent malignant glioma (MG) after previous radio(chemo)therapy. However, only limited data is available allowing for a precise selection of patients suitable for re-treatment in regard to safety and efficacy.

METHODS

Using the department database, 58 patients with two courses of percutaneous radiation were identified. Besides classical dose-volume histogram (DVH) parameters equivalent uniform dose (EUD) values were calculated for the tumor and organs at risk (OARs), retrospectively analyzed and correlated to survival outcome parameters. Cumulative EUD values were also calculated in all cases where previous OAR DVHs were available.

RESULTS

Median follow-up was 265 days and no relevant toxicity was observed after re-irradiation in our patient cohort during follow-up. Time interval between first and second irradiation was regularly above 6 months. As a conservative estimation of the cumulative EUD to the OARs, the EUDs of first and second irradiation were added. Median cumulative EUD to the optic chiasm was 48.8 Gy (range, 2.5-76.5 Gy), 57.4 Gy (range, 2.7-75.3 Gy) to the brainstem, 20.9/22.1 Gy (range, 0.0-68.3 Gy) to the right/left optic nerve and 73.8 Gy (range, 64.9-77.3 Gy) to the brain. No correlation between treated volume and survival was seen.

CONCLUSIONS

This study provides retrospective estimates on cumulative doses at the OARs. EUD values are derived and may serve as reference for further studies, including planning studies where specific constraints are needed.

An adaptive planning strategy for station parameter optimized radiation therapy (SPORT): Segmentally boosted VMAT

PURPOSE

Conventional volumetric modulated arc therapy (VMAT) discretizes the angular space into equally spaced control points during planning and then optimizes the apertures and weights of the control points. The aperture at an angle in between two control points is obtained through interpolation. This approach tacitly ignores the differential need for intensity modulation of different angles. As such, multiple arcs are often required, which may oversample some angle(s) and undersample others. The purpose of this work is to develop a segmentally boosted VMAT scheme to eliminate the need for multiple arcs in VMAT treatment with improved dose distribution and∕or delivery efficiency.

METHODS

The essence of the new treatment scheme is how to identify the need of individual angles for intensity modulation and to provide the necessary beam intensity modulation for those beam angles that need it. We introduce a "demand metric" at each control point to decide which station or control points need intensity modulation. To boost the modulation at selected stations, additional segments are added in the vicinity of the selected stations. The added segments are then optimized together with the original set of station or control points as a whole. The authors apply the segmentally boosted planning technique to four previously treated clinical cases: two head and neck (HN) cases, one prostate case, and one liver case. The proposed planning technique is compared with conventional one-arc and two-arc VMAT.

RESULTS

The proposed segmentally boosted VMAT technique achieves better critical structure sparing than one-arc VMAT with similar or better target coverage in all four clinical cases. The segmentally boosted VMAT also outperforms two-arc VMAT for the two complicated HN cases, yet with ∼30% reduction in the machine monitor units (MUs) relative to two-arc VMAT, which leads to less leakage∕scatter dose to the patient and can potentially translate into faster dose delivery. For the less challenging prostate and liver cases, similar critical structure sparing as the two-arc VMAT plans was obtained using the segmentally boosted VMAT. The benefit for the two simpler cases is the reduction of MUs and improvement of treatment delivery efficiency.

CONCLUSIONS

Segmentally boosted VMAT achieves better dose conformality and∕or reduced MUs through effective consideration of the need of individual beam angles for intensity modulation. Elimination of the need for multiple arcs in rotational arc therapy while improving the dose distribution should lead to improved workflow and treatment efficacy, thus may have significant implication to radiation oncology practice.

Maximizing the biological effect of proton dose delivered with scanned beams via inhomogeneous daily dose distributions

PURPOSE

Biological effect of radiation can be enhanced with hypofractionation, localized dose escalation, and, in particle therapy, with optimized distribution of linear energy transfer (LET). The authors describe a method to construct inhomogeneous fractional dose (IFD) distributions, and evaluate the potential gain in the therapeutic effect from their delivery in proton therapy delivered by pencil beam scanning.

METHODS

For 13 cases of prostate cancer, the authors considered hypofractionated courses of 60 Gy delivered in 20 fractions. (All doses denoted in Gy include the proton's mean relative biological effectiveness (RBE) of 1.1.) Two types of plans were optimized using two opposed lateral beams to deliver a uniform dose of 3 Gy per fraction to the target by scanning: (1) in conventional full-target plans (FTP), each beam irradiated the entire gland, (2) in split-target plans (STP), beams irradiated only the respective proximal hemispheres (prostate split sagittally). Inverse planning yielded intensity maps, in which discrete position control points of the scanned beam (spots) were assigned optimized intensity values. FTP plans preferentially required a higher intensity of spots in the distal part of the target, while STP, by design, employed proximal spots. To evaluate the utility of IFD delivery, IFD plans were generated by rearranging the spot intensities from FTP or STP intensity maps, separately as well as combined using a variety of mixing weights. IFD courses were designed so that, in alternating fractions, one of the hemispheres of the prostate would receive a dose boost and the other receive a lower dose, while the total physical dose from the IFD course was roughly uniform across the prostate. IFD plans were normalized so that the equivalent uniform dose (EUD) of rectum and bladder did not increase, compared to the baseline FTP plan, which irradiated the prostate uniformly in every fraction. An EUD-based model was then applied to estimate tumor control probability (TCP) and normal tissue complication probability (NTCP). To assess potential local RBE variations, LET distributions were calculated with Monte Carlo, and compared for different plans. The results were assessed in terms of their sensitivity to uncertainties in model parameters and delivery.

RESULTS

IFD courses included equal number of fractions boosting either hemisphere, thus, the combined physical dose was close to uniform throughout the prostate. However, for the entire course, the prostate EUD in IFD was higher than in conventional FTP by up to 14%, corresponding to the estimated increase in TCP to 96% from 88%. The extent of gain depended on the mixing factor, i.e., relative weights used to combine FTP and STP spot weights. Increased weighting of STP typically yielded a higher target EUD, but also led to increased sensitivity of dose to variations in the proton's range. Rectal and bladder EUD were same or lower (per normalization), and the NTCP for both remained below 1%. The LET distributions in IFD also depended strongly on the mixing weights: plans using higher weight of STP spots yielded higher LET, indicating a potentially higher local RBE.

CONCLUSIONS

In proton therapy delivered by pencil beam scanning, improved therapeutic outcome can potentially be expected with delivery of IFD distributions, while administering the prescribed quasi-uniform dose to the target over the entire course. The biological effectiveness of IFD may be further enhanced by optimizing the LET distributions. IFD distributions are characterized by a dose gradient located in proximity of the prostate's midplane, thus, the fidelity of delivery would depend crucially on the precision with which the proton range could be controlled.

Retrospective evaluation of dosimetric quality for prostate carcinomas treated with 3D conformal, intensity modulated and volumetric modulated arc radiotherapy

INTRODUCTION

This study examines and compares the dosimetric quality of radiotherapy treatment plans for prostate carcinoma across a cohort of 163 patients treated across five centres: 83 treated with three-dimensional conformal radiotherapy (3DCRT), 33 treated with intensity modulated radiotherapy (IMRT) and 47 treated with volumetric modulated arc therapy (VMAT).

METHODS

Treatment plan quality was evaluated in terms of target dose homogeneity and organs at risk (OAR), through the use of a set of dose metrics. These included the mean, maximum and minimum doses; the homogeneity and conformity indices for the target volumes; and a selection of dose coverage values that were relevant to each OAR. Statistical significance was evaluated using two-tailed Welch's T-tests. The Monte Carlo DICOM ToolKit software was adapted to permit the evaluation of dose metrics from DICOM data exported from a commercial radiotherapy treatment planning system.

RESULTS

The 3DCRT treatment plans offered greater planning target volume dose homogeneity than the other two treatment modalities. The IMRT and VMAT plans offered greater dose reduction in the OAR: with increased compliance with recommended OAR dose constraints, compared to conventional 3DCRT treatments. When compared to each other, IMRT and VMAT did not provide significantly different treatment plan quality for like-sized tumour volumes.

CONCLUSIONS

This study indicates that IMRT and VMAT have provided similar dosimetric quality, which is superior to the dosimetric quality achieved with 3DCRT.

Critical dose and toxicity index of organs at risk in radiotherapy: analyzing the calculated effects of modified dose fractionation in non-small cell lung cancer

To increase the efficacy of radiotherapy for non-small cell lung cancer (NSCLC), many schemes of dose fractionation were assessed by a new "toxicity index" (I), which allows one to choose the fractionation schedules that produce less toxic treatments. Thirty-two patients affected by non resectable NSCLC were treated by standard 3-dimensional conformal radiotherapy (3DCRT) with a strategy of limited treated volume. Computed tomography datasets were employed to re plan by simultaneous integrated boost intensity-modulated radiotherapy (IMRT). The dose distributions from plans were used to test various schemes of dose fractionation, in 3DCRT as well as in IMRT, by transforming the dose-volume histogram (DVH) into a biological equivalent DVH (BDVH) and by varying the overall treatment time. The BDVHs were obtained through the toxicity index, which was defined for each of the organs at risk (OAR) by a linear quadratic model keeping an equivalent radiobiological effect on the target volume. The less toxic fractionation consisted in a severe/moderate hyper fractionation for the volume including the primary tumor and lymph nodes, followed by a hypofractionation for the reduced volume of the primary tumor. The 3DCRT and IMRT resulted, respectively, in 4.7% and 4.3% of dose sparing for the spinal cord, without significant changes for the combined-lungs toxicity (p < 0.001). Schedules with reduced overall treatment time (accelerated fractionations) led to a 12.5% dose sparing for the spinal cord (7.5% in IMRT), 8.3% dose sparing for V20 in the combined lungs (5.5% in IMRT), and also significant dose sparing for all the other OARs (p < 0.001). The toxicity index allows to choose fractionation schedules with reduced toxicity for all the OARs and equivalent radiobiological effect for the tumor in 3DCRT, as well as in IMRT, treatments of NSCLC.

Equivalence of Gyn GEC-ESTRO guidelines for image guided cervical brachytherapy with EUD-based dose prescription

Background

To establish a generalized equivalent uniform dose (gEUD) -based prescription method for Image Guided Brachytherapy (IGBT) that reproduces the Gyn GEC-ESTRO WG (GGE) prescription for cervix carcinoma patients on CT images with limited soft tissue resolution.

Methods

The equivalence of two IGBT planning approaches was investigated in 20 patients who received external beam radiotherapy (EBT) and 5 concomitant high dose rate IGBT treatments. The GGE planning strategy based on dose to the most exposed 2 cm³ (D2cc) was used to derive criteria for the gEUD-based planning of the bladder and rectum. The safety of gEUD constraints in terms of GGE criteria was tested by maximizing dose to the gEUD constraints for individual fractions.

Results

The gEUD constraints of 3.55 Gy for the rectum and 5.19 Gy for the bladder were derived. Rectum and bladder gEUD-maximized plans resulted in D2cc averages very similar to the initial GGE criteria. Average D2ccs and EUDs from the full treatment course were comparable for the two techniques within both sets of normal tissue constraints. The same was found for the tumor doses.

Conclusions

The derived gEUD criteria for normal organs result in GGE-equivalent IGBT treatment plans. The gEUD-based planning considers the entire dose distribution of organs in contrast to a single dose-volume-histogram point.

Correlation of dose computed using different algorithms with local control following stereotactic ablative radiotherapy (SABR)-based treatment of non-small-cell lung cancer

PURPOSE

To retrospectively compute dose distributions for lung cancer patients treated with SABR, and to correlate dose distributions with outcome using a tumor control probability (TCP) model.

METHODS

Treatment plans for 133 NSCLC patients treated using 12 Gy/fxn × 4 (BED=106 Gy), and planned using a pencil-beam (1D-equivalent-path-length, EPL-1D) algorithm were retrospectively re-calculated using model-based algorithms (including convolution/superposition, Monte Carlo). 4D imaging was performed to manage motion. TCP was computed using the Marsden model and associations between dose and outcome were inferred.

RESULTS

Mean D95 reductions of 20% (max.=33%) were noted with model-based algorithms (relative to EPL-1D) for the smallest tumors (PTV<20 cm(3)), corresponding to actual delivered D95 BEDs of ≈ 60-85 Gy. For larger tumors (PTV>100 cm(3)), D95 reductions were ≈ 10% (BED>100 Gy). Mean lung doses (MLDs) were 15% lower for model-based algorithms for PTVs<20 cm(3). No correlation between tumor size and 2-year local control rate was observed clinically, consistent with TCP calculations, both of which were ≈ 90% across all PTV bins.

CONCLUSION

Results suggest that similar control rates might be achieved for smaller tumors using lower BEDs relative to larger tumors. However, more studies with larger patient cohorts are necessary to confirm this possible finding.

Dosimetric comparison between IMRT delivery modes: Step-and-shoot, sliding window, and volumetric modulated arc therapy - for whole pelvis radiation therapy of intermediate-to-high risk prostate adenocarcinoma

THIS STUDY WAS PERFORMED TO EVALUATE DOSIMETRIC DIFFERENCES BETWEEN CURRENT INTENSITY MODULATED RADIATION THERAPY (IMRT) DELIVERY MODES: Step-and-shoot (SS), sliding window (SW), and volumetric modulated arc therapy (VMAT). Plans for 15 prostate cancer patients with 10 MV photon beams using each IMRT mode were generated. Patients had three planning target volumes (PTVs) including prostate, prostate plus seminal vesicles, and pelvic lymphatics. Dose volume histograms (DVHs) of PTVs and organs at risk (OARs), tumor control probability (TCP) and normal tissue complication probabilities (NTCPs), conformation number, and monitor units (MUs) used were compared. Statistical analysis was performed using the analysis of variance (ANOVA) technique. The TCPs were < 99% with insignificant differences among modalities (P > 0.99). Doses to all critical structures were higher on average with SW method compared to SS, but insignificant. NTCP values were lowest for VMAT in all structures excepting bladder. Normal tissue volumes receiving doses in the 20-30 Gy range were reduced for VMAT compared to SS. Percentage of MUs required for VMAT to deliver a comparable plan to SS and SW was at least 40% less. In conclusion, similar target coverage and normal tissue doses were found by the three compared modes and the dosimetric differences were small.

Biological effects and equivalent doses in radiotherapy: A software solution

BACKGROUND

The limits of TDF (time, dose, and fractionation) and linear quadratic models have been known for a long time. Medical physicists and physicians are required to provide fast and reliable interpretations regarding delivered doses or any future prescriptions relating to treatment changes.

AIM

We, therefore, propose a calculation interface under the GNU license to be used for equivalent doses, biological doses, and normal tumor complication probability (Lyman model).

MATERIALS AND METHODS

THE METHODOLOGY USED DRAWS FROM SEVERAL SOURCES: the linear-quadratic-linear model of Astrahan, the repopulation effects of Dale, and the prediction of multi-fractionated treatments of Thames.

RESULTS AND CONCLUSIONS

The results are obtained from an algorithm that minimizes an ad-hoc cost function, and then compared to an equivalent dose computed using standard calculators in seven French radiotherapy centers.

Endorectal balloons in the post prostatectomy setting: do gains in stability lead to more predictable dosimetry?

PURPOSE

To perform a comparative study assessing potential benefits of endorectal-balloons (ERB) in post-prostatectomy patients.

METHOD AND MATERIALS

Ten retrospective post-prostatectomy patients treated without ERB and ten prospective patients treated with the ERB in situ were recruited. All patients received IMRT and IGRT using kilovoltage cone-beam computed tomography (kVCBCT). kVCBCT datasets were registered to the planning dataset, recontoured and the original plan recalculated on the kVCBCTs to recreate anatomical conditions during treatment. The imaging, structure and dose data were imported into in-house software for the assessment of geometric variation and cumulative equivalent uniform dose (EUD) in the two groups.

RESULTS

The difference in location (ΔCOV) for the bladder between planning and each CBCT was similar for each group. The range of mean ΔCOV for the rectum was 0.15-0.58 cm and 0.15-0.59 cm for the non-ERB and ERB groups. For superior-CTV and inferior-CTV the difference between planned and delivered D95% (mean ± SD) for the non-ERB group was 2.1 ± 6.0 Gy and -0.04 ± 0.20 Gy. While for the ERB group the difference in D95% was 8.7 ± 12.6 Gy and 0.003 ± 0.104 Gy.

CONCLUSIONS

The use of ERBs in the post-prostatectomy setting did improve geometric reproducibility of the target and surrounding normal tissues, however no improvement in dosimetric stability was observed for the margins employed.

Novel parameter predicting grade 2 rectal bleeding after iodine-125 prostate brachytherapy combined with external beam radiation therapy

PURPOSE

To propose a novel parameter predicting rectal bleeding on the basis of generalized equivalent uniform doses (gEUD) after (125)I prostate brachytherapy combined with external beam radiation therapy and to assess the predictive value of this parameter.

METHODS AND MATERIALS

To account for differences among radiation treatment modalities and fractionation schedules, rectal dose-volume histograms (DVHs) of 369 patients with localized prostate cancer undergoing combined therapy retrieved from corresponding treatment planning systems were converted to equivalent dose-based DVHs. The gEUDs for the rectum were calculated from these converted DVHs. The total gEUD (gEUDsum) was determined by a summation of the brachytherapy and external-beam radiation therapy components.

RESULTS

Thirty-eight patients (10.3%) developed grade 2+ rectal bleeding. The grade 2+ rectal bleeding rate increased as the gEUDsum increased: 2.0% (2 of 102 patients) for <70 Gy, 10.3% (15 of 145 patients) for 70-80 Gy, 15.8% (12 of 76 patients) for 80-90 Gy, and 19.6% (9 of 46 patients) for >90 Gy (P=.002). Multivariate analysis identified age (P=.024) and gEUDsum (P=.000) as risk factors for grade 2+ rectal bleeding.

CONCLUSIONS

Our results demonstrate gEUD to be a potential predictive factor for grade 2+ late rectal bleeding after combined therapy for prostate cancer.

Interplay effects in proton scanning for lung: a 4D Monte Carlo study assessing the impact of tumor and beam delivery parameters

Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4-167.1 cc) and motion amplitude (2.9-30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (τes), time required for magnet settling (τss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior-inferior motion amplitude alone. Larger spot sizes (σ ~ 9-16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (σ ~ 2-4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes respectively. The interplay effect is highly patient specific, depending on the motion amplitude, tumor location and the delivery parameters. Large degradations of the dose distribution in a single fraction were observed, but improved significantly using conventional fractionation.

A quantitative analysis of intensity-modulated radiation therapy plans and comparison of homogeneity indices for the treatment of gynecological cancers

The aim of present study was to evaluate the intensity-modulated radiation therapy (IMRT) plans using different homogeneity and conformity indices in gynecological cancers, as well as to compare and find out the most reliable and accurate measure of the dose homogeneity among the available indices. In this study, a cohort of 12 patients were registered for evaluation, those receiving dynamic IMRT treatment on Clinac-2300C/D linear accelerator with 15-Mega Voltage (MV) photon beam. Dynamic IMRT plans were created on Eclipse treatment planning system with Helios dose volume optimization software. Homogeneity indices (HI) such as H index, modified H index, HI index, modified HI index, and S-index (sigma-index) proposed by M Yoon et al. (2007) were calculated and compared. The values of S-index vary from 1.63 to 2.99. The results indicate that the H and HI indices and their modified versions may not provide the correct dose homogeneity information, but the S-index provides accurate information about the dose homogeneity in the Planning Target Volume (PTV). Each plan was compared with 6-MV photon energy on the basis of S-index and conformity index (CI). Organs at risk (OAR) doses with 6-MV and 15-MV beams were also reported.

Dose mapping sensitivity to deformable registration uncertainties in fractionated radiotherapy - applied to prostate proton treatments

BACKGROUND

Calculation of accumulated dose in fractionated radiotherapy based on spatial mapping of the dose points generally requires deformable image registration (DIR). The accuracy of the accumulated dose thus depends heavily on the DIR quality. This motivates investigations of how the registration uncertainty influences dose planning objectives and treatment outcome predictions.A framework was developed where the dose mapping can be associated with a variable known uncertainty to simulate the DIR uncertainties in a clinical workflow. The framework enabled us to study the dependence of dose planning metrics, and the predicted treatment outcome, on the DIR uncertainty. The additional planning margin needed to compensate for the dose mapping uncertainties can also be determined. We applied the simulation framework to a hypofractionated proton treatment of the prostate using two different scanning beam spot sizes to also study the dose mapping sensitivity to penumbra widths.

RESULTS

The planning parameter most sensitive to the DIR uncertainty was found to be the target D95. We found that the registration mean absolute error needs to be ≤0.20 cm to obtain an uncertainty better than 3% of the calculated D95 for intermediate sized penumbras. Use of larger margins in constructing PTV from CTV relaxed the registration uncertainty requirements to the cost of increased dose burdens to the surrounding organs at risk.

CONCLUSIONS

The DIR uncertainty requirements should be considered in an adaptive radiotherapy workflow since this uncertainty can have significant impact on the accumulated dose. The simulation framework enabled quantification of the accuracy requirement for DIR algorithms to provide satisfactory clinical accuracy in the accumulated dose.

The use and QA of biologically related models for treatment planning: short report of the TG-166 of the therapy physics committee of the AAPM

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

Software for quantitative analysis of radiotherapy: overview, requirement analysis and design solutions

Radiotherapy is a fast-developing discipline which plays a major role in cancer care. Quantitative analysis of radiotherapy data can improve the success of the treatment and support the prediction of outcome. In this paper, we first identify functional, conceptional and general requirements on a software system for quantitative analysis of radiotherapy. Further we present an overview of existing radiotherapy analysis software tools and check them against the stated requirements. As none of them could meet all of the demands presented herein, we analyzed possible conceptional problems and present software design solutions and recommendations to meet the stated requirements (e.g. algorithmic decoupling via dose iterator pattern; analysis database design). As a proof of concept we developed a software library "RTToolbox" following the presented design principles. The RTToolbox is available as open source library and has already been tested in a larger-scale software system for different use cases. These examples demonstrate the benefit of the presented design principles.

A comparison of the dosimetric effects of intrafraction motion on step-and-shoot, compensator, and helical tomotherapy-based IMRT

Intrafraction motion during intensity‐modulated radiation therapy can cause differences between the planned and delivered patient dose. The magnitude of these differences is dependent on a number of variables, including the treatment modality. This study was designed to compare the relative susceptibility of plans generated with three different treatment modalities to intrafraction motion. The dosimetric effects of motion were calculated using computational algorithms for seven lung tumor patients. Three delivery techniques — MLC‐based step‐and‐shoot (SNS), beam attenuating compensators, and helical tomotherapy (HT) — were investigated. In total 840 motion‐encoded dose‐volume histograms (DVHs) were calculated for various combinations of CTV margins and sinusoidal CTV motion including CTV offsets. DVH‐based metrics (e.g., D95% and D05%) were used to score plan degradations. For all three modalities, dosimetric degradations were typically smaller than 3% if the CTV displacement was smaller than the CTV margin. For larger displacements, technique and direction‐specific sensitivities existed. While the HT plans show similar D95% degradations for motion in the SI and AP directions, SNS and compensator plans showed larger D95% degradations for motion in the SI direction than for motion in the AP direction. When averaged over all motion/margin combinations, compensator plans resulted in 0.9% and 0.6% smaller D95% reductions compared to SNS and HT plans, respectively. These differences were statistically significant. No statistically significant differences in D95% degradations were found between SNS and HT for data averaged over all margin and motion track combinations. For CTV motion that is larger than the CTV margin, the dosimetric impact on the CTV varies with treatment technique and the motion direction. For the cases presented here, the effect of motion on CTV dosimetry was statistically smaller for compensator deliveries than SNS and HT, likely due to the absence of the interplay effect which is present for the more dynamic treatment deliveries. The differences between modalities were, however, small and might not be clinically significant. As expected, margins that envelop the CTV motion provide dosimetric protection against motion for all three modalities.

PACS numbers: 87.53.Jw, 87.55.dk, 87.55.de

Dose-response relationships for an atomized symptom of fecal incontinence after gynecological radiotherapy

PURPOSE

The aim of this study was to investigate what bowel organ and delivered dose levels are most relevant for the development of 'emptying of all stools into clothing without forewarning' so that the related dose-responses could be derived as an aid in avoiding this distressing symptom in the future.

MATERIAL AND METHODS

Of the 77 gynecological cancer survivors treated with radiotherapy (RT) for gynecological cancer, 13 developed the symptom. The survivors were treated between 1991 and 2003. The anal-sphincter region, the rectum, the sigmoid and the small intestines were all delineated and the dose-volume histograms were exported for each patient. The dose-volume parameters were estimated fitting the data to the Relative Seriality (RS), the Lyman and the generalized Equivalent Uniform Dose (gEUD) model.

RESULTS

The dose-response parameters for all three models and four organs at risk (OARs) were estimated. The data from the sigmoid fits the studied models best: D50 was 58.8 and 59.5 Gy (RS, Lyman), γ50 was 1.60 and 1.57 (RS, Lyman), s was 0.32, n was 0.13 and a was 7.7 (RS, Lyman, gEUD). The estimated volume parameters indicate that the investigated OARs behave serially for this endpoint. Our results for the three models studied indicate that they have the same predictive power (similar LL values) for the symptom as a function of the dose for all investigated OARs.

CONCLUSIONS

In our study, the anal-sphincter region and sigmoid fit our data best, but all OARs were found to have steep dose-responses for 'emptying of all stools into clothing without forewarning' and thus, the outcome can be predicted with an NTCP model. In addition, the dose to the four studied OARs may be considered when minimizing the risk of the symptom.

Degradation of target coverage due to inter-fraction motion during intensity-modulated proton therapy of prostate and elective targets

Internal target and organ motion during treatment is a challenge in radiotherapy (RT) of the prostate and the involved elective targets, with residual motion being present also following image-guidance strategies. The aim of this study was to investigate organ motion-induced dose degradations for the prostate, seminal vesicle and the pelvic lymph node when treating these targets with proton therapy, using different image-guidance and delivery strategies.

MATERIAL AND METHODS

Four patients were selected from a larger series as they displayed large inter-fractional variation in bladder and rectum volume. Intensity-modulated proton therapy plans were generated using both simultaneous integrated and sequential boost delivery. For each technique, three isotropic margin expansions (in the range of 4-10 mm) were evaluated for the clinical target volume of prostate (CTV-p), seminal vesicles (CTV-sv) and lymph nodes (CTV-ln). Simulation of the dose degradations for all treatment plans were based on dose re-calculations for the 8-9 repeat CTs available for each patient, after applying rigid registrations to reproduce set-up based on either intra-prostatic fiducials or bony anatomy.

RESULTS

The simulated dose received by 99% of the target volume (D(99)) and generalized equivalent dose (gEUD) showed substantial inter-patient variations. For 40% of the investigated scenarios, the patient average simulated D(99) for all targets were within 2 GyE from the planned dose. The largest difference between simulated and planned dose was seen for the CTV-sv when using SIB delivery, with an average relative reduction in D(99) of 13% and 15% for the largest margin expansion, when positioned using fiducials and bony anatomy, respectively.

CONCLUSIONS

The most severe dose degradations were found for CTV-sv, but they were also evident for CTV-ln. The degradations could not be completely resolved, neither by using the largest margin expansion nor with the choice of set-up. With fiducial set-up CTV-p was robust against the inter-fraction changes.

Treatment plan comparison using grading analysis based on clinical judgment

PURPOSE

In this work we explore a method named clinical grading analysis (CGA) which is based on clinical assessments performed by radiation oncologists (ROs). The purpose is to investigate how useful the method is for treatment plan comparisons, and how the CGA results correlate with dosimetric evaluation parameters, traditionally used for treatment plan comparisons.

MATERIAL AND METHODS

Helical tomotherapy (HTT) and seven-beam step-and-shoot intensity modulated radiation therapy (SS-IMRT) plans were compared and assessed by 10 experienced ROs for 23 patient cases. A CGA was performed where the plans were graded based on how the ROs thought they compared to each other. The resulting grades from the CGA were analyzed and compared to dose-volume statistics and equivalent uniform dose (EUD) data.

RESULTS

For eight of the 23 cases the CGA revealed a significant difference between the HTT and the SS-IMRT plans, five cases were in favor of HTT, and three in favor of SS-IMRT. Comparing the dose-volume statistics and EUD-data with the result from the CGA showed that CGA results correlated well with dose-volume statistics for cases regarding difference in target coverage or doses to organs at risk. The CGA results also correlated well with EUD-data for cases with difference in clinical target volume (CTV) coverage but the correlation for cases with difference in planning target volume (PTV) coverage was not as clear.

CONCLUSIONS

This study presents CGA as a useful method of comparing radiotherapy treatment plans. The proposed method offers a formalized way of introducing and evaluating the implementation of new radiotherapy techniques in a clinical setting. The CGA identify patients that have a clinical benefit of one or the other of the advanced treatment techniques available to them, i.e. in this study HTT and SS-IMRT, which facilitates a more optimal use of a clinics' advanced treatment resources.

Maximizing dosimetric benefits of IMRT in the treatment of localized prostate cancer through multicriteria optimization planning

We examine the quality of plans created using multicriteria optimization (MCO) treatment planning in intensity-modulated radiation therapy (IMRT) in treatment of localized prostate cancer. Nine random cases of patients receiving IMRT to the prostate were selected. Each case was associated with a clinically approved plan created using Corvus. The cases were replanned using MCO-based planning in RayStation. Dose-volume histogram data from both planning systems were presented to 2 radiation oncologists in a blinded evaluation, and were compared at a number of dose-volume points. Both physicians rated all 9 MCO plans as superior to the clinically approved plans (p<10(-5)). Target coverage was equivalent (p = 0.81). Maximum doses to the prostate and bladder and the V50 and V70 to the anterior rectum were reduced in all MCO plans (p<0.05). Treatment planning time with MCO took approximately 60 minutes per case. MCO-based planning for prostate IMRT is efficient and produces high-quality plans with good target homogeneity and sparing of the anterior rectum, bladder, and femoral heads, without sacrificing target coverage.

Local tumor control probability to evaluate an applicator-guided volumetric-modulated arc therapy solution as alternative of 3D brachytherapy for the treatment of the vaginal vault in patients affected by gynecological cancer

The purpose of this study was to evaluate the applicator‐guided volumetric‐modulated arc therapy (AGVMAT) solution as an alternative to high‐dose‐rate brachytherapy (HDR‐BRT) treatment of the vaginal vault in patients with gynecological cancer (GC). AGVMAT plans for 51 women were developed. The volumetric scans used for plans were obtained with an implanted CT‐compatible vaginal cylinder which provides spatial registration and immobilization of the gynecologic organs. Dosimetric and radiobiological comparisons for planning target volume (PTV) and organs at risk (OARs) were performed by means of a dose‐volume histogram (DVH), equivalent uniform dose (EUD), and local tumor control probability (LTCP). In addition, the integral dose and the overall delivery time, were evaluated. The HDR‐BRT averages of EUD and minimum LTCP were significantly higher than those of AGVMAT. Doses for the OARs were comparable for the bladder and sigmoid, while, although HDR‐BRT was able to better spare the bowel, AGVMAT provided a significant reduction of d2cc, d1cc, and dmax (p<0.01) for the rectum. AGVMAT integral doses were higher than HDR‐BRT with low values in both cases. Delivery times were about two or three times higher for HDR‐BRT with respect to the single arc technique (AGVMAT1) and dual arc technique (AGVMAT2), respectively. The applicator‐guided volumetric‐modulated arc therapy seems to have the potential of improving rectum avoidance. However, brachytherapy improves performance in terms of PTV coverage, as demonstrated by a greater EUD and better LTCP curves.

PACS number: 87

Changes in radiobiological parameters in 131Cs permanent prostate implants

In prostate permanent implants using 131Cs seeds, the prostatic edema developed during the implantation procedure, increases the separation between the seeds. This leads to a decrease in the prostate coverage and thus causes an edema induced dose reduction, which results in an increase in tumour cell surviving fraction (SF) with a corresponding decrease in tumour control probability (TCP). To investigate the impact of edema on the SF and the TCP, the expression of the SF of the linear quadratic (LQ) model was extended to account for the effects of edema using the exponential nature of edema resolution and the dose delivered to the edematous prostate. The SF and the TCP for edematous prostate implants were calculated for 31 patients who underwent real time 131Cs permanent seed implantation. The dose delivered to the edematous prostate was calculated to compute the SF and the TCP for these patients for edema half lives (EHL) ranging from 4 days to 34 days and for edemas of magnitudes (M0) varying from 5 to 60% of the actual prostate volume.

A reduction in the dose delivered to the edematous prostate was found with the increase of EHL and edema magnitude which results in an increase of the SF, and corresponding decrease in the TCP. The dose reductions in 131Cs implants varied from 1.1% (for EHL = 4 days and M0 = 5%) to 32.3% (for EHL = 34 days and M0 = 60%). These are higher than the dose reduction in 125I implants, which vary from 0.3% (for EHL = 4 days and M0 = 5%) to 17.5% (for EHL = 34 days and M0 = 60%). As edema half life increased from 4 days to 34 days and edema magnitude increased from 5 to 60% the SF increased by 4.57 log, and the TCP decreased by 0.80. Compensation of edema induced increase in the SF and decrease in the TCP in 131Cs seed implants should be carefully done by redefining seed positions with the guidance of post-needle plans. The presented model in this study can be used to estimate the SF or the TCP for pre plan or real time permanent prostate implants using day 0 post-implant CT images.

Radiobiological Evaluation of Breast Cancer Radiotherapy Accounting for the Effects of Patient Positioning and Breathing in Dose Delivery. A Meta Analysis

In breast cancer radiotherapy, significant discrepancies in dose delivery can contribute to underdosage of the tumor or overdosage of normal tissue, which is potentially related to a reduction of local tumor control and an increase of side effects. To study the impact of these factors in breast cancer radiotherapy, a meta analysis of the clinical data reported by Mavroidis et al. (2002) in Acta Oncol (41:471–85), showing the patient setup and breathing uncertainties characterizing three different irradiation techniques, were employed. The uncertainties in dose delivery are simulated based on fifteen breast cancer patients (5 mastectomized, 5 resected with negative node involvement (R-) and 5 resected with positive node involvement (R+)), who were treated by three different irradiation techniques, respectively. The positioning and breathing effects were taken into consideration in the determination of the real dose distributions delivered to the CTV and lung in each patient. The combined frequency distributions of the positioning and breathing distributions were obtained by convolution. For each patient the effectiveness of the dose distribution applied is calculated by the Poisson and relative seriality models and a set of parameters that describe the dose-response relations of the target and lung. The three representative radiation techniques are compared based on radiobiological measures by using the complication-free tumor control probability, P+ and the biologically effective uniform dose, D̿ concepts. For the Mastectomy case, the average P+ values of the planned and delivered dose distributions are 93.8% for a D̿CTV of 51.8 Gy and 85.0% for a D̿CTV of 50.3 Gy, respectively. The respective total control probabilities, PB values are 94.8% and 92.5%, whereas the corresponding total complication probabilities, PI values are 0.9% and 7.4%. For the R- case, the average P+ values are 89.4% for a D̿CTV of 48.9 Gy and 88.6% for a D̿CTV of 49.0 Gy, respectively. The respective PB values are 89.8% and 89.9%, whereas the corresponding PI values are 0.4% and 1.2%. For the R+ case, the average P+ values are 86.1% for a D̿CTV of 49.2 Gy and 85.5% for a D̿CTV of 49.1 Gy, respectively. The respective PB values are 90.2% and 90.1%, whereas the corresponding PI values are 4.1% and 4.6%. The combined effects of positioning uncertainties and breathing can introduce a significant deviation between the planned and delivered dose distributions in lung in breast cancer radiotherapy. The positioning and breathing uncertainties do not affect much the dose distribution to the CTV. The simulated delivered dose distributions show larger lung complication probabilities than the treatment plans. This means that in clinical practice the true expected complications are underestimated. Radiation pneumonitis of Grade 1–2 is more frequent and any radiotherapy optimization should use this as a more clinically relevant endpoint.

A Comparison of Critical Structure Dose and Toxicity Risks in Patients with Low Grade Gliomas Treated with IMRT versus Proton Radiation Therapy

Proton therapy offers dosimetric advantage of decreased dose to non-target tissues. This study explored the potential benefits of proton radiation therapy versus photon based intensity modulated radiation therapy (IMRT) for patients with low grade gliomas (LGG) through dosimetric comparison and biological modeling of potential radiation-induced toxicities. Eleven patients were treated with fractionated proton radiation therapy on a prospective protocol assessing for feasibility and treatment toxicity of proton radiation therapy in patients with LGG. IMRT treatment plans were created for each patient using the same CT planning data set and defined structures. The prescription dose to the clinical target volume (CTV) was 54 Gy(RBE). The toxicity risk of IMRT and protons was estimated based upon equivalent uniform dose (EUD) and normal tissue complication probability (NTCP) modeling. The risk of secondary tumors for each modality was estimated. Proton EUD for most immediate normal tissue structures was between 10–20 Gy lower than the EUD delivered by IMRT. However, the difference in NTCP was negligible for both modalities. The mean excess risk of proton radiation-induced second tumor in the brain per 10,000 cases per year is 47 (range 11–83), while the mean risk for IMRT is 106 (range 70–134). The mean ratio of excess risk IMRT/protons is 2.2 (range 1.6–6.5), demonstrating that the risk of secondary tumors is consistently higher for IMRT. Proton therapy effectively reduces the dose to surrounding normal tissues in LGG patients. IMRT has a twofold higher risk of secondary intracranial tumors as compared to proton therapy. In most cases, NTCP is negligible for both modalities. The benefit of proton therapy over IMRT may be more substantial in patients with tumors in proximity to critical structures.

Dosimetric assessment of prostate cancer patients through principal component analysis (PCA)

The aims of this study were twofold: first, to determine the impact of variance in dose-volume histograms (DVH) on patient-specific toxicity after 2 high-dose fractions in a sample of 22 men with prostate cancer; and second, to compare the effectiveness of traditional DVH analysis and principal component analysis (PCA) in predicting rectum and urethra toxicity. A series of 22 patients diagnosed with prostate adenocarcinoma was treated with 45 Gy external beam and 20 Gy dose rate brachytherapy. Principal component analysis was applied to model the shapes of the rectum and urethra dose-volume histograms. We used logistic regression to measure the correlations between the principal components and the incidence of rectal bleeding and urethra stricture. We also calculated the equivalent uniform dose (EUD) and normal tissue complication probability (NTCP) for the urethra and rectum, and tumor control probability (TCP) for the prostate using BioSuite software. We evaluated their correlations with rectal and urethra toxicity. The rectum DVHs are well described by one principal component (PC1), which accounts for 93.5% of the variance in their shapes. The urethra DVHs are described by two principal components, PC1 and PC2, which account for 94.98% and 3.15% of the variance, respectively. Multivariate exact logistic regression suggests that urethra PC2 is a good predictor of stricture, with Nagelkerke's R2 estimated at 0.798 and a Wald criterion of 5.421 (p &lt; 0.021). The average NTCPs were 0.06% ± 0.04% and 1.25% ± 0.22% for the rectum and urethra, respectively. The average TCP was 85.29% ± 2.28%. This study suggests that principal component analysis can be used to identify the shape variation in dose-volume histograms, and that the principal components can be correlated with the toxicity of a treatment plan based on multivariate analysis. The principal components are also correlated with traditional dosimetric parameters.

Quantification of dose nonuniformities by voxel-based dosimetry in patients receiving 90Y-ibritumomab-tiuxetan

Abstract Objective: To assess the impact of nonuniform dose distribution within lesions and tumor-involved organs of patients receiving Zevalin(®), and to discuss possible implications of equivalent uniform biological effective doses (EU-BED) on treatment efficacy and toxicity. MATLAB™ -based software for voxel-based dosimetry was adopted for this purpose.

METHODS

Eleven lesions from seven patients with either indolent or aggressive non-Hodgkin lymphoma were analyzed, along with four organs with disease. Absorbed doses were estimated by a direct integration of single-voxel kinetic data from serial tomographic images. After proper corrections, differential BED distributions and surviving cell fractions were estimated, allowing for the calculation of EU-BED. To quantify dose uniformity in each target area, a heterogeneity index was defined.

RESULTS

Average doses were below those prescribed by conventional radiotherapy to eradicate lymphoma lesions. Dose heterogeneity and effect on tumor control varied among lesions, with no apparent relation to tumor mass. Although radiation doses to involved organs were safe, unexpected liver toxicity occurred in one patient who presented with a pattern of diffuse infiltration.

CONCLUSION

Voxel-based dosimetry and radiobiologic modeling can be successfully applied to lesions and tumor-involved organs, representing a methodological advance over estimation of mean absorbed doses. However, effects on tumor control and organ toxicity still cannot be easily predicted.

EQUIVALENT UNIFORM DOSE SENSITIVITY TO CHANGES IN ABSORBED DOSE DISTRIBUTION

Treatment planning in external beam radiation therapy (EBRT) utilizes dose volume histograms (DVHs) as optimization and evaluation tools. They present the fraction of planning target volume (PTV) receiving more than a given absorbed dose, against the absorbed dose values, and a number of radiobiological indices can be computed with their help. Equivalent uniform dose (EUD) is the absorbed dose that, uniformly imparted, would yield the same biological effect on a tumor as the dose distribution described by the DVH. Uncertainty and missing information can affect the dose distribution, therefore DVHs can be modeled as samples from a set of possible outcomes. This work studies the sensitivity of the EUD index when a small change in absorbed dose distribution takes place. EUD is treated as a functional on the set of DVHs. Defining a Lévy distance on this set and using a suitable expansion of the functional, a very simple expression for a bound on the variation of EUD when the dose distribution changes is found. This bound is easily interpreted in terms of standard treatment planning practice.

Proton therapy for hepatocellular carcinoma

Proton radiotherapy has seen an increasing role in the treatment of hepatocellular carcinoma (HCC). Historically, external beam radiotherapy has played a very limited role in HCC due to a high incidence of toxicity to surrounding normal structures. The ability to deliver a high dose of radiation to the tumor is a key factor in improving outcomes in HCC. Advances in photon radiotherapy have improved dose conformity and allowed dose escalation to the tumor. However, despite these advances there is still a large volume of normal liver that receives a considerable radiation dose during treatment. Proton beams do not have an exit dose along the beam path once they enter the body. The inherent physical attributes of proton radiotherapy offer a way to maximize tumor control via dose escalation while avoiding excessive radiation to the remaining liver, thus increasing biological effectiveness. In this review we discuss the physical attributes and rationale for proton radiotherapy in HCC. We also review recent literature regarding clinical outcomes of using proton radiotherapy for the treatment of HCC.

Adaptive IMRT using a multiobjective evolutionary algorithm integrated with a diffusion-invasion model of glioblastoma

We demonstrate a patient-specific method of adaptive IMRT treatment for glioblastoma using a multiobjective evolutionary algorithm (MOEA). The MOEA generates spatially optimized dose distributions using an iterative dialogue between the MOEA and a mathematical model of tumor cell proliferation, diffusion and response. Dose distributions optimized on a weekly basis using biological metrics have the potential to substantially improve and individualize treatment outcomes. Optimized dose distributions were generated using three different decision criteria for the tumor and compared with plans utilizing standard dose of 1.8 Gy/fraction to the CTV (T2-visible MRI region plus a 2.5 cm margin). The sets of optimal dose distributions generated using the MOEA approach the Pareto Front (the set of IMRT plans that delineate optimal tradeoffs amongst the clinical goals of tumor control and normal tissue sparing). MOEA optimized doses demonstrated superior performance as judged by three biological metrics according to simulated results. The predicted number of reproductively viable cells 12 weeks after treatment was found to be the best target objective for use in the MOEA.