An algorithm for maximizing the probability of complication-free tumour control in radiation therapy

NTCP Modeling and Dose-Volume Correlations of Significant Hematocrit Drop 3 Months After Prostate Radiation Therapy

Purpose

Our purpose was to determine and model the dose-response relations of different parts of the pelvis regarding the endpoint of hematocrit level drop after pelvic radiation therapy (RT).

Methods and Materials

Two hundred and twenty-one patients treated with RT for prostate adenocarcinoma between 2014 and 2016 were included. All patients had complete blood counts collected at baseline and 3 months post-RT. The net difference of hematocrit level post-RT versus baseline was calculated, and the level of the 15th percentiles defined the thresholds of response in each case. The doses to 8 different pelvic structures were derived and fitted to the hematocrit levels using the relative seriality normal tissue complication probability model and the biologically equivalent uniform dose (D = ).

Results

Pelvic structures that correlated with significant decreases in hematocrit were the os coxae bilaterally superior to the acetabulum (OCUB), the total os coxae bilaterally, and the bone volume of the whole pelvis. The structure showing the highest correlation was OCUB with a maximum area under the curve (AUC) of 0.74. For V20 Gy < 30% the odds ratio was 9.8 with 95% CI of 2.9 to 32.9. For mean dose (Dmean) to OCUB, an AUC of 0.73 was observed where the dose threshold was 23 Gy and the odds ratio was 2.7 and 95% CI 1.3 to 5.6. The values for the D50, γ, and s parameters of the relative seriality model were 26.9 Gy (25.9-27.9), 1.3 (1.2-2.2), and 0.12 (0.10-0.83), respectively. The AUC of D = was 0.73 and patients with D = to OCUB ≥ 27 Gy had 8.2 times higher rate of significant hematocrit drop versus <27 Gy.

Conclusions

These findings confirm the association of radiation-induced damage to pelvic bone marrow with a drop in hematocrit. A threshold of V20 Gy < 30%, Dmean < 23 Gy, or D = < 27 Gy to OCUB may significantly reduce the risk for this endpoint.

Evaluation of the clinical impact of the differences between planned and delivered dose in prostate cancer radiotherapy based on CT-on-rails IGRT and patient-reported outcome scores

PURPOSE

To estimate the clinical impact of differences between delivered and planned dose using dose metrics and normal tissue complication probability (NTCP) modeling.

METHODS

Forty-six consecutive patients with prostate adenocarcinoma between 2010 and 2015 treated with intensity-modulated radiation therapy (IMRT) and who had undergone computed tomography on rails imaging were included. Delivered doses to bladder and rectum were estimated using a contour-based deformable image registration method. The bladder and rectum NTCP were calculated using dose-response parameters applied to planned and delivered dose distributions. Seven urinary and gastrointestinal symptoms were prospectively collected using the validated prostate cancer symptom indices patient reported outcome (PRO) at pre-treatment, weekly treatment, and post-treatment follow-up visits. Correlations between planned and delivered doses against PRO were evaluated in this study.

RESULTS

Planned mean doses to bladder and rectum were 44.9 ± 13.6 Gy and 42.8 ± 7.3 Gy, while delivered doses were 46.1 ± 13.4 Gy and 41.3 ± 8.7 Gy, respectively. D_10cc for rectum was 64.1 ± 7.6 Gy for planned and 60.1 ± 9.3 Gy for delivered doses. NTCP values of treatment plan were 22.3% ± 8.4% and 12.6% ± 5.9%, while those for delivered doses were 23.2% ± 8.4% and 9.9% ± 8.3% for bladder and rectum, respectively. Seven of 25 patients with follow-up data showed urinary complications (28%) and three had rectal complications (12%). Correlations of NTCP values of planned and delivered doses with PRO follow-up data were random for bladder and moderate for rectum (0.68 and 0.67, respectively).

CONCLUSION

Sensitivity of bladder to clinical variations of dose accumulation indicates that an automated solution based on a DIR that considers inter-fractional organ deformation could recommend intervention. This is intended to achieve additional rectum sparing in cases that indicate higher than expected dose accumulation early during patient treatment in order to prevent acute severity of bowel symptoms.

Spatial descriptions of radiotherapy dose: normal tissue complication models and statistical associations

For decades, dose-volume information for segmented anatomy has provided the essential data for correlating radiotherapy dosimetry with treatment-induced complications. Dose-volume information has formed the basis for modelling those associations via normal tissue complication probability (NTCP) models and for driving treatment planning. Limitations to this approach have been identified. Many studies have emerged demonstrating that the incorporation of information describing the spatial nature of the dose distribution, and potentially its correlation with anatomy, can provide more robust associations with toxicity and seed more general NTCP models. Such approaches are culminating in the application of computationally intensive processes such as machine learning and the application of neural networks. The opportunities these approaches have for individualising treatment, predicting toxicity and expanding the solution space for radiation therapy are substantial and have clearly widespread and disruptive potential. Impediments to reaching that potential include issues associated with data collection, model generalisation and validation. This review examines the role of spatial models of complication and summarises relevant published studies. Sources of data for these studies, appropriate statistical methodology frameworks for processing spatial dose information and extracting relevant features are described. Spatial complication modelling is consolidated as a pathway to guiding future developments towards effective, complication-free radiotherapy treatment.

Radiobiological assessment of nasopharyngeal cancer IMRT using various collimator angles and non-coplanar fields

Aim:

The aim of this study was to evaluate clinical efficacy and radiobiological outcome of intensity-modulated radiation therapy (IMRT) modalities using various collimator angles and non-coplanar fields for nasopharyngeal cancer (NPC).

Materials and methods:

A 70-Gy planning target volume dose was administered for 30 NPC patients referred for IMRT. Standard IMRT plans were constructed based on the target and organs at risk (OARs) volume; and dose constraints recommended by Radiation Therapy Oncology Group (RTOG). Using various collimator angles and non-coplanar fields, 11 different additional IMRT protocols were investigated. Homogeneity indexes (HIs) and conformation numbers (CNs) were calculated. Poisson and relative seriality models were utilised for estimating tumour control probability (TCP) and normal tissue complication probabilities (NTCPs), respectively.

Results:

Various collimator angles and non-coplanar fields had no significant effect on HI, CN and TCP, while significant effects were noted for some OARs, with a maximum mean dose (Dmax). No significant differences were observed among the calculated NTCPs of all the IMRT protocols. However, the protocol with 10° collimator angle (for five fields out of seven) and 8° couch angle had the lowest NTCP. Furthermore, the standard and some of non-coplanar IMRT protocols led to the reduction in OARs Dmax.

Conclusions:

Using appropriate standard/non-coplanar IMRT protocols for NPC treatment could potentially reduce the dose to the OARs and the probability of inducing secondary cancer in patients.

Generalized approach for radiotherapy treatment planning by optimizing projected health outcome: preliminary results for prostate radiotherapy patients

Research in cancer care increasingly focuses on survivorship issues, e.g. managing disease- and treatment-related morbidity and mortality occurring during and after treatment. This necessitates innovative approaches that consider treatment side effects in addition to tumor cure. Current treatment-planning methods rely on constrained iterative optimization of dose distributions as a surrogate for health outcomes. The goal of this study was to develop a generally applicable method to directly optimize projected health outcomes. We developed an outcome-based objective function to guide selection of the number, angle, and relative fluence weight of photon and proton radiotherapy beams in a sample of ten prostate-cancer patients by optimizing the projected health outcome. We tested whether outcome-optimized radiotherapy (OORT) improved the projected longitudinal outcome compared to dose-optimized radiotherapy (DORT) first for a statistically significant majority of patients, then for each individual patient. We assessed whether the results were influenced by the selection of treatment modality, late-risk model, or host factors. The results of this study revealed that OORT was superior to DORT. Namely, OORT maintained or improved the projected health outcome of photon- and proton-therapy treatment plans for all ten patients compared to DORT. Furthermore, the results were qualitatively similar across three treatment modalities, six late-risk models, and 10 patients. The major finding of this work was that it is feasible to directly optimize the longitudinal (i.e. long- and short-term) health outcomes associated with the total (i.e. therapeutic and stray) absorbed dose in all of the tissues (i.e. healthy and diseased) in individual patients. This approach enables consideration of arbitrary treatment factors, host factors, health endpoints, and times of relevance to cancer survivorship. It also provides a simpler, more direct approach to realizing the full beneficial potential of cancer radiotherapy.

A DNA Repair-Based Model of Cell Survival with Important Clinical Consequences

This work provides a description of a new interaction, cross-section-based model for radiation-induced cellular inactivation, sublethal damage, DNA repair and cell survival, with the ability to more accurately elucidate different radiation-response phenomena. The principal goal of this work is to describe the damage-induction cross sections, as well as repair and survival, as Poisson processes with two main types of damage: mild damage that can be rapidly handled by the most basic repair processes; and more complex damage requiring longer repair times and the high-fidelity homologous recombination (HR) repair process to ensure accuracy and safety in the survival. This work is unique in its use of Poisson statistics to quantify the main repairable cell compartments that are exposed to simple and more complex sublethal hits, the cross section of which determines what is homologically and non-homologically repairable. The new method is applied to central radiation damage and survival data, such as in vitro cellular repair and survival with key DNA repair genes knocked out, low-dose hypersensitivity (LDHS), change in survival over the cell cycle, and variation with linear energy transfer (LET) for densely ionizing ions, all results supporting our basic assumptions. Among the results, it was shown that less than 1% of the simple DSBs are lethal at approximately 2 Gy and below for sparsely ionizing radiations, but their δ-electron track ends of between 1.5 and 0.5 keV can deliver 0.5 MGy to a few hundred nm3 volumes, mainly due to multiple scatter detours and multiple secondary electrons. They can cause dual double-strand breaks (DSBs) on the periphery of nucleosomes that are the most common multiply damaged sites, with an average of 1-2 δ-electron track ends per cell nucleus at 2 Gy. LDHS is most likely due to the normal lack of fast, efficient repair of sublethal damage below approximately 0.5 Gy, and requires largely intact key DNA repair genes to achieve significant repair recovery at higher doses. The new repair model describes this phenomenon quite accurately. Cells with key non-homologous end joining (NHEJ) genes knocked-out, lose LDHS but provoke HR repair, and cells with HR genes knocked out may lose some LDHS, but provoke NHEJ repair. The DNA duplication during the S phase results in a direct doubling as well of the total and sublethal hit cross sections. For the lowest LET carbon ions, NHEJ is reduced to where it is almost eliminated at maximum relative biological effectiveness (RBE), while HR is induced more than by X rays, due to complex damage and misrepair of DSBs produced by numerous δ electrons. The use of a lower LET such as electrons or photons during the final week of radiation treatment may potentially maximize complication-free cure. Optimally-designed weekly fractionation schedules are proposed to maximize the DNA repair potential in normal tissues. Additionally, the optimal therapeutic ion species, LET, apoptosis and permanent growth arrest/senescence window is identified with helium, lithium and boron ions and LETs at approximately 15-55 eV/nm, to maximize these quantities in the tumor and minimize them in the normal tissues, resulting in a very high probability of complication-free cure.

Justification and optimization of radiation exposures: a new framework to aggregate arbitrary detriments and benefits

Myriad radiation effects, including benefits and detriments, complicate justifying and optimizing radiation exposures. The purpose of this study was to develop a comprehensive conceptual framework and corresponding quantitative methods to aggregate the detriments and benefits of radiation exposures to individuals, groups, and populations. In this study, concepts from the ICRP for low dose were integrated with clinical techniques focused on high dose to develop a comprehensive figure of merit (FOM) that takes into account arbitrary host- and exposure-related factors, endpoints, and time points. The study built on existing methods with three new capabilities: application to individuals, groups, and populations; extension to arbitrary numbers and types of endpoints; and inclusion of limitation, where relevant. The FOM was applied to three illustrative exposure situations: emergency response, diagnostic imaging, and cancer radiotherapy, to evaluate its utility in diverse settings. The example application to radiation protection revealed the FOM's utility in optimizing the benefits and risks to a population while keeping individual exposures below applicable regulatory limits. Examples in diagnostic imaging and cancer radiotherapy demonstrated the FOM's utility for guiding population- and patient-specific decisions in medical applications. The major finding of this work is that it is possible to quantitatively combine the benefits and detriments of any radiation exposure situation involving an individual or population to perform cost-effectiveness analyses using the ICRP key principles of radiation protection. This FOM fills a chronic gap in the application of radiation-protection theory, i.e., limitations of generalized frameworks to algorithmically justify and optimize radiation exposures. This new framework potentially enhances objective optimization and justification, especially in complex exposure situations.

Proposals of models for new formulations of the current complication-free cure (P+) and uncomplicated tumor control probability (UTCP) concepts, and total normal tissue complication probability of late complications

This study proposes phenomenological models for total normal tissue complication probability (TNTCP) and NTCP0. NTCP0 is a new acronym for reformulating the current complication-free cure (P+) and uncomplicated tumor control probability (UTCP) concepts, and TNTCP will reformulate the current NTCP involving multiple organs at risks. The current probabilistic concepts are incoherently formulated with mathematical operations of tumor control probability (TCP) and normal tissue complication probability (NTCP) that are associated with different stochastic processes and random variables. NTCP0 is equal to NTCP₀ (normal tissue non-complication probability) that is calculated as the ratio of a number of patients of a population without late complications and a total of them. As a cumulative distribution function (CDF) of late complications, TNTCP = sum(NTCP_i), where NTCP_i is the NTCP of the i^th late complication. TNTCP is also a new acronym, and the probabilistic complement of NTCP0, then NTCP0 = 100% - TNTCP. The NTCP0/TNTCP (D(d)) proposing models are based on the relationship between the NTCP0/TNTCP and total dose (D = n×d; where d = dose per fraction, and n = number of fractions). TNTCP(D) model will be correlated with LKB model (the normal CDF) that is an increasing function; and NTCP0(D) model with a decreasing function, which additionally will define clear limits of three possible regions for NTCP0: 0 and 100% deterministic, and a stochastic. These models are function D, which is widely used for characterizing radiation therapies.

Incorporating biological modeling into patient‐specific plan verification

Purpose

Dose–volume histogram (DVH) measurements have been integrated into commercially available quality assurance systems to provide a metric for evaluating accuracy of delivery in addition to gamma analysis. We hypothesize that tumor control probability and normal tissue complication probability calculations can provide additional insight beyond conventional dose delivery verification methods.

Methods

A commercial quality assurance system was used to generate DVHs of treatment plan using the planning CT images and patient‐specific QA measurements on a phantom. Biological modeling was performed on the DVHs produced by both the treatment planning system and the quality assurance system.

Results

The complication‐free tumor control probability, P₊, has been calculated for previously treated intensity modulated radiotherapy (IMRT) patients with diseases in the following sites: brain (−3.9% ± 5.8%), head‐neck (+4.8% ± 8.5%), lung (+7.8% ± 1.3%), pelvis (+7.1% ± 12.1%), and prostate (+0.5% ± 3.6%).

Conclusion

Dose measurements on a phantom can be used for pretreatment estimation of tumor control and normal tissue complication probabilities. Results in this study show how biological modeling can be used to provide additional insight about accuracy of delivery during pretreatment verification.

Dosimetric and Radiobiological Evaluation of Patient Setup Accuracy in Head-and-neck Radiotherapy Using Daily Computed Tomography-on-rails-based Corrections

Introduction:

This study evaluates treatment plans aiming at determining the expected impact of daily patient setup corrections on the delivered dose distribution and plan parameters in head-and-neck radiotherapy.

Materials and Methods:

In this study, 10 head-and-neck cancer patients are evaluated. For the evaluation of daily changes of the patient internal anatomy, image-guided radiation therapy based on computed tomography (CT)-on-rails was used. The daily-acquired CT-on-rails images were deformedly registered to the CT scan that was used during treatment planning. Two approaches were used during data analysis (“cascade” and “one-to-all”). The dosimetric and radiobiological differences of the dose distributions with and without patient setup correction were calculated. The evaluation is performed using dose–volume histograms; the biologically effective uniform dose () and the complication-free tumor control probability (P₊) were also calculated. The dose–response curves of each target and organ at risk (OAR), as well as the corresponding P₊ curves, were calculated.

Results:

The average difference for the “one-to-all” case is 0.6 ± 1.8 Gy and for the “cascade” case is 0.5 ± 1.8 Gy. The value of P₊ was lowest for the cascade case (in 80% of the patients).

Discussion:

Overall, the lowest P_I is observed in the one-to-all cases. Dosimetrically, CT-on-rails data are not worse or better than the planned data.

Conclusions:

The differences between the evaluated “one-to-all” and “cascade” dose distributions were small. Although the differences of those doses against the “planned” dose distributions were small for the majority of the patients, they were large for given patients at risk and OAR.

Three-dimensional dose prediction and validation with the radiobiological gamma index based on a relative seriality model for head-and-neck IMRT

This study proposes a quality assurance (QA) method incorporating radiobiological factors based on the QUANTEC-determined tumor control probability and the normal tissue complication probability (NTCP) of head-and-neck intensity-modulated radiation therapy (HN-IMRT). Per-beam measurements were conducted for 20 cases using a 2D detector array. Three-dimensional predicted dose distributions within targets and organs at risk were reconstructed based on the per-beam QA results derived from differences between planned and measured doses. Under the predicted dose distributions, the differences between the physical and radiobiological gamma indices (PGI and RGI, respectively) based on the relative seriality (RS) model were evaluated. The NTCP values in the RS and Niemierko models were compared. The dose covers 98% (D98%) of the clinical target volume (CTV) decreased by 3.2% (P < 0.001), and the mean dose of the ipsilateral parotid increased by 6.3% (P < 0.001) compared with the original dose. RGI passing rates in the CTV and brain stem were greater than PGI ones by 5.8% (P < 0.001) and 2.0% (P < 0.001), respectively. The RS model's average NTCP values for the ipsilateral and contralateral parotids under the original dose were smaller than those of the Niemierko model by 9.0% (P < 0.001) and 7.0% (P < 0.001), respectively. The 3D predicted dose evaluation with RGI based on the RS model was introduced for QA of HN-IMRT, leading to dose evaluation for each organ with consideration of the radiobiological effect. This method constitutes a rational way to perform QA of HN-IMRT in clinical practice.

Modeling Radiotherapy Induced Normal Tissue Complications: An Overview beyond Phenomenological Models

An overview of radiotherapy (RT) induced normal tissue complication probability (NTCP) models is presented. NTCP models based on empirical and mechanistic approaches that describe a specific radiation induced late effect proposed over time for conventional RT are reviewed with particular emphasis on their basic assumptions and related mathematical translation and their weak and strong points.

Comparison of composite prostate radiotherapy plan doses with dependent and independent boost phases

Prostate cases commonly consist of dual phase planning with a primary plan followed by a boost. Traditionally, the boost phase is planned independently from the primary plan with the risk of generating hot or cold spots in the composite plan. Alternatively, boost phase can be planned taking into account the primary dose. The aim of this study was to compare the composite plans from independently and dependently planned boosts using dosimetric and radiobiological metrics. Ten consecutive prostate patients previously treated at our institution were used to conduct this study on the Raystation™ 4.0 treatment planning system. For each patient, two composite plans were developed: a primary plan with an independently planned boost and a primary plan with a dependently planned boost phase. The primary plan was prescribed to 54 Gy in 30 fractions to the primary planning target volume (PTV1) which includes prostate and seminal vesicles, while the boost phases were prescribed to 24 Gy in 12 fractions to the boost planning target volume (PTV2) that targets only the prostate. PTV coverage, max dose, median dose, target conformity, dose homogeneity, dose to OARs, and probabilities of benefit, injury, and complication-free tumor control (P+) were compared. Statistical significance was tested using either a 2-tailed Student's t-test or Wilcoxon signed-rank test. Dosimetrically, the composite plan with dependent boost phase exhibited smaller hotspots, lower maximum dose to the target without any significant change to normal tissue dose. Radiobiologically, for all but one patient, the percent difference in the P+ values between the two methods was not significant. A large percent difference in P+ value could be attributed to an inferior primary plan. The benefits of considering the dose in primary plan while planning the boost is not significant unless a poor primary plan was achieved.

Disregarding RBE variation in treatment plan comparison may lead to bias in favor of proton plans

PURPOSE

Currently in proton radiation therapy, a constant relative biological effectiveness (RBE) equal to 1.1 is assumed. The purpose of this study is to evaluate the impact of disregarding variations in RBE on the comparison of proton and photon treatment plans.

METHODS

Intensity modulated treatment plans using photons and protons were created for three brain tumor cases with the target situated close to organs at risk. The proton plans were optimized assuming a standard RBE equal to 1.1, and the resulting linear energy transfer (LET) distribution for the plans was calculated. In the plan evaluation, the effect of a variable RBE was studied. The RBE model used considers the RBE variation with dose, LET, and the tissue specific parameter α/β of photons. The plan comparison was based on dose distributions, DVHs and normal tissue complication probabilities (NTCPs).

RESULTS

Under the assumption of RBE=1.1, higher doses to the tumor and lower doses to the normal tissues were obtained for the proton plans compared to the photon plans. In contrast, when accounting for RBE variations, the comparison showed lower doses to the tumor and hot spots in organs at risk in the proton plans. These hot spots resulted in higher estimated NTCPs in the proton plans compared to the photon plans.

CONCLUSIONS

Disregarding RBE variations might lead to suboptimal proton plans giving lower effect in the tumor and higher effect in normal tissues than expected. For cases where the target is situated close to structures sensitive to hot spot doses, this trend may lead to bias in favor of proton plans in treatment plan comparisons.

Comparing conformal, arc radiotherapy and helical tomotherapy in craniospinal irradiation planning

Currently, radiotherapy treatment plan acceptance is based primarily on dosimetric performance measures. However, use of radiobiological analysis to assess benefit in terms of tumor control and harm in terms of injury to normal tissues can be advantageous. For pediatric craniospinal axis irradiation (CSI) patients, in particular, knowing the technique that will optimize the probabilities of benefit versus injury can lead to better long-term outcomes. Twenty-four CSI pediatric patients (median age 10) were retrospectively planned with three techniques: three-dimensional conformal radiation therapy (3D CRT), volumetric-modulated arc therapy (VMAT), and helical tomotherapy (HT). VMAT plans consisted of one superior and one inferior full arc, and tomotherapy plans were created using a 5.02cm field width and helical pitch of 0.287. Each plan was normalized to 95% of target volume (whole brain and spinal cord) receiving prescription dose 23.4Gy in 13 fractions. Using an in-house MATLAB code and DVH data from each plan, the three techniques were evaluated based on biologically effective uniform dose (D=), the complication-free tumor control probability (P+), and the width of the therapeutically beneficial range. Overall, 3D CRT and VMAT plans had similar values of D= (24.1 and 24.2 Gy), while HT had a D= slightly lower (23.6 Gy). The average values of the P+ index were 64.6, 67.4, and 56.6% for 3D CRT, VMAT, and HT plans, respectively, with the VMAT plans having a statistically significant increase in P+. Optimal values of D= were 28.4, 33.0, and 31.9 Gy for 3D CRT, VMAT, and HT plans, respectively. Although P+ values that correspond to the initial dose prescription were lower for HT, after optimizing the D= prescription level, the optimal P+ became 94.1, 99.5, and 99.6% for 3D CRT, VMAT, and HT, respectively, with the VMAT and HT plans having statistically significant increases in P+. If the optimal dose level is prescribed using a radiobiological evaluation method, as opposed to a purely dosimetric one, the two IMRT techniques, VMAT and HT, will yield largest overall benefit to CSI patients by maximizing tumor control and limiting normal tissue injury. Using VMAT or HT may provide these pediatric patients with better long-term outcomes after radiotherapy.

Consequences of anorectal cancer atlas implementation in the cooperative group setting: radiobiologic analysis of a prospective randomized in silico target delineation study

PURPOSE

The aim of this study is to ascertain the subsequent radiobiological impact of using a consensus guideline target volume delineation atlas.

MATERIALS AND METHODS

Using a representative case and target volume delineation instructions derived from a proposed IMRT rectal cancer clinical trial, gross tumor volume (GTV) and clinical/planning target volumes (CTV/PTV) were contoured by 13 physician observers (Phase 1). The observers were then randomly assigned to follow (atlas) or not-follow (control) a consensus guideline/atlas for anorectal cancers, and instructed to re-contour the same case (Phase 2).

RESULTS

The atlas group was found to have increased tumor control probability (TCP) after the atlas intervention for both the CTV (p<0.0001) and PTV1 (p=0.0011) with decreasing normal tissue complication probability (NTCP) for small intestine, while the control group did not. Additionally, the atlas group had reduced variance in TCP for all target volumes and reduced variance in NTCP for the bowel. In Phase 2, the atlas group had increased TCP relative to the control for CTV (p=0.03).

CONCLUSIONS

Visual atlas and consensus treatment guideline usage in the development of rectal cancer IMRT treatment plans reduced the inter-observer radiobiological variation, with clinically relevant TCP alteration for CTV and PTV volumes.

γ+ index: A new evaluation parameter for quantitative quality assurance

PURPOSE

The accuracy of dose delivery and the evaluation of differences between calculated and delivered dose distributions, has been studied by several groups. The aim of this investigation is to extend the gamma index by including radiobiological information and to propose a new index that we will here forth refer to as the gamma plus (γ+). Furthermore, to validate the robustness of this new index in performing a quality control analysis of an IMRT treatment plan using pure radiobiological measures such as the biologically effective uniform dose (D) and complication-free tumor control probability (P+).

MATERIAL AND METHODS

A new quality assurance index, the (γ+), is proposed based on the theoretical concept of gamma index presented by Low et al. (1998). In this study, the dose difference, including the radiobiological dose information (biological effective dose, BED) is used instead of just the physical dose difference when performing the γ+ calculation. An in-house software was developed to compare different dose distributions based on the γ+ concept. A test pattern for a two-dimensional dose comparison was built using the in-house software platform. The γ+ index was tested using planar dose distributions (exported from the treatment planning system) and delivered (film) dose distributions acquired in a solid water phantom using a test pattern and a theoretical clinical case. Furthermore, a lung cancer case for a patient treated with IMRT was also selected for the analysis. The respective planar dose distributions from the treatment plan and the film were compared based on the γ+ index and were evaluated using the radiobiological measures of P+ and D.

RESULTS

The results for the test pattern analysis indicate that the γ+ index distributions differ from those of the gamma index since the former considers radiobiological parameters that may affect treatment outcome. For the theoretical clinical case, it is observed that the γ+ index varies for different treatment parameters (e.g. dose per fraction). The dose area histogram (DAH) from the plan and film dose distributions are associated with P+ values of 50.8% and 49.0%, for a D to the target of 54.0 Gy and 53.3 Gy, respectively.

CONCLUSION

The γ+ index shows advantageous properties in the quantitative evaluation of dose delivery and quality control of IMRT treatments because it includes information about the expected responses and radiobiological doses of the individual tissues.

Analysis of fractionation correction methodologies for multiple phase treatment plans in radiation therapy

PURPOSE

Radiation therapy is often delivered by multiple sequential treatment plans. For an accurate radiobiological evaluation of the overall treatment, fractionation corrections to each dose distribution must be applied before summing the three-dimensional dose matrix of each plan since the simpler approach of performing the fractionation correction to the total dose-volume histograms, obtained by the arithmetical sum of the different plans, becomes inaccurate for more heterogeneous dose patterns. In this study, the differences between these two fractionation correction methods, named here as exact (corrected before) and approximate (after summation), respectively, are assessed for different cancer types.

METHODS

Prostate, breast, and head and neck (HN) tumor patients were selected to quantify the differences between two fractionation correction methods (the exact vs the approximate). For each cancer type, two different treatment plans were developed using uniform (CRT) and intensity modulated beams (IMRT), respectively. The responses of the target and normal tissue were calculated using the Poisson linear-quadratic-time model and the relative seriality model, respectively. All treatments were radiobiologically evaluated and compared using the complication-free tumor control probability (P+), the biologically effective uniform dose (D) together with common dosimetric criteria.

RESULTS

For the prostate cancer patient, an underestimation of around 14%-15% in P+ was obtained when the fractionation correction was applied after summation compared to the exact approach due to significant biological and dosimetric variations obtained between the two fractionation correction methods in the involved lymph nodes. For the breast cancer patient, an underestimation of around 3%-4% in the maximum dose in the heart was obtained. Despite the dosimetric differences in this organ, no significant variations were obtained in treatment outcome. For the HN tumor patient, an underestimation of about 5% in treatment outcome was obtained for the CRT plan as a result of an underestimation of the planning target volume control probability by about 10%. An underestimation of about 6% in the complication probability of the right parotid was also obtained. For all the other organs at risk, dosimetric differences of up to 4% were obtained but with no significant impact in the expected clinical outcome. However, for the IMRT plan, an overestimation in P+ of 4.3% was obtained mainly due to an underestimation of the complication probability of the left and right parotids (2.9% and 5.8%, respectively).

CONCLUSIONS

The use of the exact fractionation correction method, which is applying fractionation correction on the separate dose distributions of a multiple phase treatment before their summation was found to have a significant expected clinical impact. For regions of interest that are irradiated with very heterogeneous dose distributions and significantly different doses per fraction in the different treatment phases, the exact fractionation correction method needs to be applied since a significant underestimation of the true patient outcome can be introduced otherwise.

A gEUD-based inverse planning technique for HDR prostate brachytherapy: feasibility study

PURPOSE

The purpose of this work was to study the feasibility of a new inverse planning technique based on the generalized equivalent uniform dose for image-guided high dose rate (HDR) prostate cancer brachytherapy in comparison to conventional dose-volume based optimization.

METHODS

The quality of 12 clinical HDR brachytherapy implants for prostate utilizing HIPO (Hybrid Inverse Planning Optimization) is compared with alternative plans, which were produced through inverse planning using the generalized equivalent uniform dose (gEUD). All the common dose-volume indices for the prostate and the organs at risk were considered together with radiobiological measures. The clinical effectiveness of the different dose distributions was investigated by comparing dose volume histogram and gEUD evaluators.

RESULTS

Our results demonstrate the feasibility of gEUD-based inverse planning in HDR brachytherapy implants for prostate. A statistically significant decrease in D10 or/and final gEUD values for the organs at risk (urethra, bladder, and rectum) was found while improving dose homogeneity or dose conformity of the target volume.

CONCLUSIONS

Following the promising results of gEUD-based optimization in intensity modulated radiation therapy treatment optimization, as reported in the literature, the implementation of a similar model in HDR brachytherapy treatment plan optimization is suggested by this study. The potential of improved sparing of organs at risk was shown for various gEUD-based optimization parameter protocols, which indicates the ability of this method to adapt to the user's preferences.

Quality assessment for VMAT prostate radiotherapy planning based on data envelopment analysis

The majority of commercial radiotherapy treatment planning systems requires planners to iteratively adjust the plan parameters in order to find a satisfactory plan. This iterative trial-and-error nature of radiotherapy treatment planning results in an inefficient planning process and in order to reduce such inefficiency, plans can be accepted without achieving the best attainable quality. We propose a quality assessment method based on data envelopment analysis (DEA) to address this inefficiency. This method compares a plan of interest to a set of past delivered plans and searches for evidence of potential further improvement. With the assistance of DEA, planners will be able to make informed decisions on whether further planning is required and ensure that a plan is only accepted when the plan quality is close to the best attainable one. We apply the DEA method to 37 prostate plans using two assessment parameters: rectal generalized equivalent uniform dose (gEUD) as the input and D95 (the minimum dose that is received by 95% volume of a structure) of the planning target volume (PTV) as the output. The percentage volume of rectum overlapping PTV is used to account for anatomical variations between patients and is included in the model as a non-discretionary output variable. Five plans that are considered of lesser quality by DEA are re-optimized with the goal to further improve rectal sparing. After re-optimization, all five plans improve in rectal gEUD without clinically considerable deterioration of the PTV D95 value. For the five re-optimized plans, the rectal gEUD is reduced by an average of 1.84 Gray (Gy) with only an average reduction of 0.07 Gy in PTV D95. The results demonstrate that DEA can correctly identify plans with potential improvements in terms of the chosen input and outputs.

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.

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.

The radiobiological P(+) index for pretreatment plan assessment with emphasis on four-dimensional radiotherapy modalities

PURPOSE

Radiation treatment modalities will continue to emerge that promise better clinical outcomes albeit technologically challenging to implement. An important question facing the radiotherapy community then is the need to justify the added technological effort for the clinical return. Mobile tumor radiotherapy is a typical example, where 4D tumor tracking radiotherapy (4DTRT) has been proposed over the simpler conventional modality for better results. The modality choice per patient can depend on a wide variety of factors. In this work, we studied the complication-free tumor control probability (P(+)) index, which combines the physical complexity of the treatment plan with the radiobiological characteristics of the clinical case at hand and therefore found to be useful in evaluating different treatment techniques and estimating the expected clinical effectiveness of different radiation modalities.

METHODS

4DCT volumes of 18 previously treated lung cancer patients with tumor motion and size ranging from 2 mm to 15 mm and from 4 cc to 462 cc, respectively, were used. For each patient, 4D treatment plans were generated to extract the 4D dose distributions, which were subsequently used with clinically derived radiobiological parameters to compute the P(+) index per modality.

RESULTS

The authors observed, on average, a statistically significant increase in P(+) of 3.4% ± 3.8% (p < 0.003) in favor of 4DTRT. There was high variability among the patients with a <0.5% up to 13.4% improvement in P(+).

CONCLUSIONS

The observed variability in the improvement of the clinical effectiveness suggests that the relative benefit of tracking should be evaluated on a per patient basis. Most importantly, this variability could be effectively captured in the computed P(+). The index can thus be useful to discriminate and hence point out the need for a complex modality like 4DTRT over another. Besides tumor mobility, a wide range of other factors, e.g., size, location, fractionation, etc., can affect the relative benefits. Application of the P(+) objective is a simple and effective way to combine these factors in the evaluation of a treatment plan.

Automated volumetric modulated Arc therapy treatment planning for stage III lung cancer: how does it compare with intensity-modulated radio therapy?

PURPOSE

To compare the quality of volumetric modulated arc therapy (VMAT) or intensity-modulated radiation therapy (IMRT) plans generated by an automated inverse planning system with that of dosimetrist-generated IMRT treatment plans for patients with stage III lung cancer.

METHODS AND MATERIALS

Two groups of 8 patients with stage III lung cancer were randomly selected. For group 1, the dosimetrists spent their best effort in designing IMRT plans to compete with the automated inverse planning system (mdaccAutoPlan); for group 2, the dosimetrists were not in competition and spent their regular effort. Five experienced radiation oncologists independently blind-reviewed and ranked the three plans for each patient: a rank of 1 was the best and 3 was the worst. Dosimetric measures were also performed to quantitatively evaluate the three types of plans.

RESULTS

Blind rankings from different oncologists were generally consistent. For group 1, the auto-VMAT, auto-IMRT, and manual IMRT plans received average ranks of 1.6, 2.13, and 2.18, respectively. The auto-VMAT plans in group 1 had 10% higher planning tumor volume (PTV) conformality and 24% lower esophagus V70 (the volume receiving 70 Gy or more) than the manual IMRT plans; they also resulted in more than 20% higher complication-free tumor control probability (P+) than either type of IMRT plans. The auto- and manual IMRT plans in this group yielded generally comparable dosimetric measures. For group 2, the auto-VMAT, auto-IMRT, and manual IMRT plans received average ranks of 1.55, 1.75, and 2.75, respectively. Compared to the manual IMRT plans in this group, the auto-VMAT plans and auto-IMRT plans showed, respectively, 17% and 14% higher PTV dose conformality, 8% and 17% lower mean lung dose, 17% and 26% lower mean heart dose, and 36% and 23% higher P+.

CONCLUSIONS

mdaccAutoPlan is capable of generating high-quality VMAT and IMRT treatment plans for stage III lung cancer. Manual IMRT plans could achieve quality similar to auto-IMRT plans if best effort was spent.

Consideration of the likely benefit from implementation of prostate image-guided radiotherapy using current margin sizes: a radiobiological analysis

OBJECTIVE

To estimate the benefit of introduction of image-guided radiotherapy (IGRT) to prostate radiotherapy practice with current clinical target volume-planning target volume (PTV) margins of 5-10 mm.

METHODS

Systematic error data collected from 50 patients were used together with a random error of σ=3.0 mm to model non-IGRT treatment. IGRT was modelled with residual errors of Σ=σ=1.5 mm. Population tumour control probability (TCP(pop)) was calculated for two three-dimensional conformal radiotherapy techniques: two-phase and concomitant boost. Treatment volumes and dose prescriptions were ostensibly the same. The relative field sizes of the treatment techniques, distribution of systematic errors and correlations between movement axes were examined.

RESULTS

The differences in TCP(pop) between the IGRT and non-IGRT regimes were 0.3% for the two-phase and 1.5% for the concomitant boost techniques. A 2-phase plan, in each phase of which the 95% isodose conformed to its respective PTV, required fields that were 3.5 mm larger than those required for the concomitant boost plan. Despite the larger field sizes, the TCP (without IGRT) in the two-phase plan was only 1.7% higher than the TCP in the concomitant boost plan. The deviation of craniocaudal systematic errors (p=0.02) from a normal distribution, and the correlation of translations in the craniocaudal and anteroposterior directions (p<0.0001) were statistically significant.

CONCLUSIONS

The expected population benefit of IGRT for the modelled situation was too small to be detected by a clinical trial of reasonable size, although there was a significant benefit to individual patients. For IGRT to have an observable population benefit, the trial would need to use smaller margins than those used in this study. Concomitant treatment techniques permit smaller fields and tighter conformality than two phases planned separately.

Response-probability volume histograms and iso-probability of response charts in treatment plan evaluation

PURPOSE

This study aims at demonstrating a new method for treatment plan evaluation and comparison based on the radiobiological response of individual voxels. This is performed by applying them on three different cancer types and treatment plans of different conformalities. Furthermore, their usefulness is examined in conjunction with traditionally applied radiobiological and dosimetric treatment plan evaluation criteria.

METHODS

Three different cancer types (head and neck, breast and prostate) were selected to quantify the benefits of the proposed treatment plan evaluation method. In each case, conventional conformal radiotherapy (CRT) and intensity modulated radiotherapy (IMRT) treatment configurations were planned. Iso-probability of response charts was produced by calculating the response probability in every voxel using the linear-quadratic-Poisson model and the dose-response parameters of the corresponding structure to which this voxel belongs. The overall probabilities of target and normal tissue responses were calculated using the Poisson and the relative seriality models, respectively. The 3D dose distribution converted to a 2 Gy fractionation, D2(GY) and iso-BED distributions are also shown and compared with the proposed methodology. Response-probability volume histograms (RVH) were derived and compared with common dose volume histograms (DVH). The different dose distributions were also compared using the complication-free tumor control probability, P+, the biologically effective uniform dose, D, and common dosimetric criteria.

RESULTS

3D Iso-probability of response distributions is very useful for plan evaluation since their visual information focuses on the doses that are likely to have a larger clinical effect in that particular organ. The graphical display becomes independent of the prescription dose highlighting the local radiation therapy effect in each voxel without the loss of important spatial information. For example, due to the exponential nature of the Poisson distribution, cold spots in the target volumes or hot spots in the normal tissues are much easier to be identified. Response-volume histograms, as DVH, can also be derived and used for plan comparison. RVH are advantageous since by incorporating the radiobiological properties of each voxel they summarize the 3D distribution into 2D without the loss of relevant information. Thus, more clinically relevant radiobiological objectives and constraints could be defined and used in treatment planning optimization. These measures become increasingly important when dose distributions need to be designed according to the microscopic biological properties of tumor and normal tissues.

CONCLUSIONS

The proposed methods do not aim to replace quantifiers like the probabilities of total tissue response, which ultimately are the quantities of interest to evaluate treatment success. However, iso-probability of response charts and response-probability volume histograms illustrates more clearly the difference in effectiveness between different treatment plans than the information provided by alternative dosimetric data. The use of 3D iso-probability of response distributions could serve as a good descriptor of the effectiveness of a dose distribution indicating primarily the regions in a tissue that dominate its response.

A methodology for automatic intensity-modulated radiation treatment planning for lung cancer

In intensity-modulated radiotherapy (IMRT), the quality of the treatment plan, which is highly dependent upon the treatment planner's level of experience, greatly affects the potential benefits of the radiotherapy (RT). Furthermore, the planning process is complicated and requires a great deal of iteration, and is often the most time-consuming aspect of the RT process. In this paper, we describe a methodology to automate the IMRT planning process in lung cancer cases, the goal being to improve the quality and consistency of treatment planning. This methodology (1) automatically sets beam angles based on a beam angle automation algorithm, (2) judiciously designs the planning structures, which were shown to be effective for all the lung cancer cases we studied, and (3) automatically adjusts the objectives of the objective function based on a parameter automation algorithm. We compared treatment plans created in this system (mdaccAutoPlan) based on the overall methodology with plans from a clinical trial of IMRT for lung cancer run at our institution. The 'autoplans' were consistently better, or no worse, than the plans produced by experienced medical dosimetrists in terms of tumor coverage and normal tissue sparing. We conclude that the mdaccAutoPlan system can potentially improve the quality and consistency of treatment planning for lung cancer.

Investigating the Clinical Aspects of Using CT vs. CT-MRI Images during Organ Delineation and Treatment Planning in Prostate Cancer Radiotherapy

In order to apply highly conformal dose distributions, which are characterized by steep dose fall-offs, it is necessary to know the exact target location and extension. This study aims at evaluating the impact of using combined CT-MRI images in organ delineation compared to using CT images alone, on the clinical results. For 10 prostate cancer patients, the respective CT and MRI images at treatment position were acquired. The CTV was delineated using the CT and MRI images, separately, whereas bladder and rectum were delineated using the CT images alone. Based on the CT and MRI images, two CTVs were produced for each patient. The mutual information algorithm was used in the fusion of the two image sets. In this way, the structures drawn on the MRI images were transferred to the CT images in order to produce the treatment plans. For each set of structures of each patient, IMRT and 3D-CRT treatment plans were produced. The individual treatment plans were compared using the biologically effective uniform dose ([Formula: see text]) and the complication-free tumor control probability ( P+) concepts together with the DVHs of the targets and organs at risk and common dosimetric criteria. For the IMRT treatment, at the optimum dose level of the average CT and CT-MRI delineated CTV dose distributions, the P+ values are 74.7% in both cases for a [Formula: see text] of 91.5 Gy and 92.1 Gy, respectively. The respective average total control probabilities, PB are 90.0% and 90.2%, whereas the corresponding average total complication probabilities, PI are 15.3% and 15.4%. Similarly, for the 3D-CRT treatment, the average P+ values are 42.5% and 46.7%, respectively for a [Formula: see text] of 86.4 Gy and 86.7 Gy, respectively. The respective average PB values are 80.0% and 80.6%, whereas the corresponding average PI values are 37.4% and 33.8%, respectively. For both radiation modalities, the improvement mainly stems from the better sparing of rectum. According to these results, the expected clinical effectiveness of IMRT can be increased by a maximum Δ P+ of around 0.9%, whereas of 3D-CRT by about 4.2% when combined CT-MRI delineation is performed instead of using CT images alone. It is apparent that in both IMRT and 3D-CRT radiation modalities, the better knowledge of the CTV extension improved the produced dose distribution. It is shown that the CTV is irradiated more effectively, while the complication probabilities of bladder and rectum, which is the principal organs at risk, are lower in the CT-MRI based treatment plans.

Radiobiological and Dosimetric Analysis of Daily Megavoltage CT Registration on Adaptive Radiotherapy with Helical Tomotherapy

Pre-treatment patient repositioning in highly conformal image-guided radiation therapy modalities is a prerequisite for reducing setup uncertainties. In Helical Tomotherapy (HT) treatment, a megavoltage CT (MVCT) image is usually acquired to evaluate daily changes in the patient's internal anatomy and setup position. This MVCT image is subsequently compared to the kilovoltage CT (kVCT) study that was used for dosimetric planning, by applying a registration process. This study aims at investigating the expected effect of patient setup correction using the Hi-Art tomotherapy system by employing radiobiological measures such as the biologically effective uniform dose ([Formula: see text]) and the complication-free tumor control probability ( P+). A new module of the Tomotherapy software (TomoTherapy, Inc, Madison, WI) called Planned Adaptive is employed in this study. In this process the delivered dose can be calculated by using the sinogram for each delivered fraction and the registered MVCT image set that corresponds to the patient's position and anatomical distribution for that fraction. In this study, patients treated for lung, pancreas and prostate carcinomas are evaluated by this method. For each cancer type, a Helical Tomotherapy plan was developed. In each cancer case, two dose distributions were calculated using the MVCT image sets before and after the patient setup correction. The fractional dose distributions were added and renormalized to the total number of fractions planned. The dosimetric and radiobiological differences of the dose distributions with and without patient setup correction were calculated. By using common statistical measures of the dose distributions and the P+ and [Formula: see text] concepts and plotting the tissue response probabilities vs. [Formula: see text] a more comprehensive comparison was performed based on radiobiological measures. For the lung cancer case, at the clinically prescribed dose levels of the dose distributions, with and without patient setup correction, the complication-free tumor control probabilities, P+ are 48.5% and 48.9% for a [Formula: see text] of 53.3 Gy. The respective total control probabilities, PB are 56.3% and 56.5%, whereas the corresponding total complication probabilities, PI are 7.9% and 7.5%. For the pancreas cancer case, at the prescribed dose levels of the two dose distributions, the P+ values are 53.7% and 45.7% for a [Formula: see text] of 54.7 Gy and 53.8 Gy, respectively. The respective PB values are 53.7% and 45.8%, whereas the corresponding PI values are ~0.0% and 0.1%. For the prostate cancer case, at the prescribed dose levels of the two dose distributions, the P+ values are 10.9% for a [Formula: see text] of 75.2 Gy and 11.9% for a [Formula: see text] of 75.4 Gy, respectively. The respective PB values are 14.5% and 15.3%, whereas the corresponding PI values are 3.6% and 3.4%. Our analysis showed that the very good daily patient setup and dose delivery were very close to the intended ones. With the exception of the pancreas cancer case, the deviations observed between the dose distributions with and without patient setup correction were within ±2% in terms of P+. In the radiobiologically optimized dose distributions, the role of patient setup correction using MVCT images could appear to be more important than in the cases of dosimetrically optimized treatment plans were the individual tissue radiosensitivities are not precisely considered.

Comparison of the helical tomotherapy against the multileaf collimator-based intensity-modulated radiotherapy and 3D conformal radiation modalities in lung cancer radiotherapy

OBJECTIVES

The aim of this study was to compare three-dimensional (3D) conformal radiotherapy and the two different forms of IMRT in lung cancer radiotherapy.

METHODS

Cases of four lung cancer patients were investigated by developing a 3D conformal treatment plan, a linac MLC-based step-and-shoot IMRT plan and an HT plan for each case. With the use of the complication-free tumour control probability (P(+)) index and the uniform dose concept as the common prescription point of the plans, the different treatment plans were compared based on radiobiological measures.

RESULTS

The applied plan evaluation method shows the MLC-based IMRT and the HT treatment plans are almost equivalent over the clinically useful dose prescription range; however, the 3D conformal plan inferior. At the optimal dose levels, the 3D conformal treatment plans give an average P(+) of 48.1% for a effective uniform dose to the internal target volume (ITV) of 62.4 Gy, whereas the corresponding MLC-based IMRT treatment plans are more effective by an average ΔP(+) of 27.0% for a Δ effective uniform dose of 16.3 Gy. Similarly, the HT treatment plans are more effective than the 3D-conformal plans by an average ΔP(+) of 23.8% for a Δ effective uniform dose of 11.6 Gy.

CONCLUSION

A radiobiological treatment plan evaluation can provide a closer association of the delivered treatment with the clinical outcome by taking into account the dose-response relations of the irradiated tumours and normal tissues. The use of P - effective uniform dose diagrams can complement the traditional tools of evaluation to compare and effectively evaluate different treatment plans.

Toolkit for determination of dose-response relations, validation of radiobiological parameters and treatment plan optimization based on radiobiological measures

Accurately determined dose-response relations of the different tumors and normal tissues should be estimated and used in the clinic. The aim of this study is to demonstrate developed tools that are necessary for determining the dose-response parameters of tumors and normal tissues, for clinically verifying already published parameter sets using local patient materials and for making use of all this information in the optimization and comparison of different treatment plans and radiation techniques. One of the software modules (the Parameter Determination Module) is designed to determine the dose-response parameters of tumors and normal tissues. This is accomplished by performing a maximum likelihood fitting to calculate the best estimates and confidence intervals of the parameters used by different radiobiological models. Another module of this software (the Parameter Validation Module) concerns the validation and compatibility of external or reported dose-response parameters describing tumor control and normal tissue complications. This is accomplished by associating the expected response rates, which are calculated using different models and published parameter sets, with the clinical follow-up records of the local patient population. Finally, the last module of the software (the Radiobiological Plan Evaluation Module) is used for estimating and optimizing the effectiveness a treatment plan in terms of complication-free tumor control, P(+). The use of the Parameter Determination Module is demonstrated by deriving the dose-response relation of proximal esophagus from head and neck cancer radiotherapy. The application of the Parameter Validation Module is illustrated by verifying the clinical compatibility of those dose-response parameters with the examined treatment methodologies. The Radiobiological Plan Evaluation Module is demonstrated by evaluating and optimizing the effectiveness of head and neck cancer treatment plans. The results of the radiobiological evaluation are compared against dosimetric criteria. The presented toolkit appears to be very convenient and efficient for clinical implementation of radiobiological modeling. It can also be used for the development of a clinical data and health information database for assisting the performance of epidemiological studies and the collaboration between different institutions within research and clinical frameworks.

Dose-painting IMRT optimization using biological parameters

PURPOSE

Our work on dose-painting based on the possible risk characteristics for local recurrence in tumor subvolumes and the optimization of treatment plans using biological objective functions that are region-specific are reviewed.

MATERIALS AND METHODS

A series of intensity modulated dose-painting techniques are compared to their corresponding intensity modulated plans in which the entire PTV is treated to a single dose level, delivering the same equivalent uniform dose (EUD) to the entire PTV. Iso-TCP and iso-NTCP maps are introduced as a tool to aid the planner in the evaluation of the resulting non-uniform dose distributions. Iso-TCP and iso-NTCP maps are akin to iso-dose maps in 3D conformal radiotherapy. The impact of the currently limited diagnostic accuracy of functional imaging on a series of dose-painting techniques is also discussed.

RESULTS

Utilizing biological parameters (risk-adaptive optimization) in the generation of dose-painting plans results in an increase in the therapeutic ratio as compared to conventional dose-painting plans in which optimization techniques based on physical dose are employed.

CONCLUSION

Dose-painting employing biological parameters appears to be a promising approach for individualized patient- and disease-specific radiotherapy.

Radiobiological evaluation of the influence of dwell time modulation restriction in HIPO optimized HDR prostate brachytherapy implants

Purpose

One of the issues that a planner is often facing in HDR brachytherapy is the selective existence of high dose volumes around some few dominating dwell positions. If there is no information available about its necessity (e.g. location of a GTV), then it is reasonable to investigate whether this can be avoided. This effect can be eliminated by limiting the free modulation of the dwell times. HIPO, an inverse treatment plan optimization algorithm, offers this option. In treatment plan optimization there are various methods that try to regularize the variation of dose non-uniformity using purely dosimetric measures. However, although these methods can help in finding a good dose distribution they do not provide any information regarding the expected treatment outcome as described by radiobiology based indices.

Material and methods

The quality of 12 clinical HDR brachytherapy implants for prostate utilizing HIPO and modulation restriction (MR) has been compared to alternative plans with HIPO and free modulation (without MR). All common dose-volume indices for the prostate and the organs at risk have been considered together with radiobiological measures. The clinical effectiveness of the different dose distributions was investigated by calculating the response probabilities of the tumors and organs-at-risk (OARs) involved in these prostate cancer cases. The radiobiological models used are the Poisson and the relative seriality models. Furthermore, the complication-free tumor control probability, P₊ and the biologically effective uniform dose (D¯¯) were used for treatment plan evaluation and comparison.

Results

Our results demonstrate that HIPO with a modulation restriction value of 0.1-0.2 delivers high quality plans which are practically equivalent to those achieved with free modulation regarding the clinically used dosimetric indices. In the comparison, many of the dosimetric and radiobiological indices showed significantly different results. The modulation restricted clinical plans demonstrated a lower total dwell time by a mean of 1.4% that was proved to be statistically significant (p = 0.002). The HIPO with MR treatment plans produced a higher P₊ by 0.5%, which stemmed from a better sparing of the OARs by 1.0%.

Conclusions

Both the dosimetric and radiobiological comparison shows that the modulation restricted optimization gives on average similar results with the optimization without modulation restriction in the examined clinical cases. Concluding, based on our results, it appears that the applied dwell time regularization technique is expected to introduce a minor improvement in the effectiveness of the optimized HDR dose distributions.

A graphic user interface toolkit for specification, report and comparison of dose-response relations and treatment plans using the biologically effective uniform dose

A toolkit (BEUDcal) has been developed for evaluating the effectiveness and for predicting the outcome of treatment plans by calculating the biologically effective uniform dose (BEUD) and complication-free tumor control probability. The input for the BEUDcal is the differential dose-volume histograms of organs exported from the treatment planning system. A clinical database is built for the dose-response parameters of different tumors and normal tissues. Dose-response probabilities of all the examined organs are illustrated together with the corresponding BEUDs and the P(+) values. Furthermore, BEUDcal is able to generate a report that simultaneously presents the radiobiological evaluation together with the physical dose indices, showing the complementary relation between the physical and radiobiological treatment plan analysis performed by BEUDcal. Comparisons between treatment plans for helical tomotherapy and multileaf collimator-based intensity modulated radiotherapy of a lung patient were demonstrated to show the versatility of BEUDcal in the assessment and report of dose-response relations.

Expected clinical impact of the differences between planned and delivered dose distributions in helical tomotherapy for treating head and neck cancer using helical megavoltage CT images

Helical Tomotherapy (HT) has become increasingly popular over the past few years. However, its clinical efficacy and effectiveness continues to be investigated. Pre-treatment patient repositioning in highly conformal image-guided radiation therapy modalities is a prerequisite for reducing setup uncertainties. A MVCT image set has to be acquired to account for daily changes in the patient's internal anatomy and setup position. Furthermore, a comparison should be performed to the kVCT study used for dosimetric planning, by a registration process which results in repositioning the patient according to specific transitional and rotational shifts. Different image registration techniques may lead to different repositioning of the patient and, as a result, to varying delivered doses. This study aims to investigate the expected effect of patient setup correction using the Hi-Art tomotherapy system by employing radiobiological measures such as the biologically effective uniform dose (BEUD) and the complication-free tumor control probability (P+). In this study, a typical case of lung cancer with metastatic head & neck disease was investigated by developing a Helical Tomotherapy plan. For the Tomotherapy HiArt plan, the dedicated Tomotherapy treatment planning station was used. Three dose distributions (planned and delivered with and without patient setup correction) were compared based on radiobiological measures by using the P+ index and the BEUD concept as the common prescription point of the plans and plotting the tissue response probabilities against the mean target dose for a range of prescription doses. The applied plan evaluation method shows that in this cancer case the planned and delivered dose distributions with and without patient setup correction give a P+ of 81.6%, 80.9% and 72.2%, for a BEUD to the planning target volume (PTV) of 78.0Gy, 77.7Gy and 75.4Gy, respectively. The corresponding tumor control probabilities are 86.3%, 85.1% and 75.1%, whereas the total complication probabilities are 4.64%, 4.20% and 2.89%, respectively. HT can encompass the often large PTV required while minimizing the volume of the organs at risk receiving high dose. However, the effectiveness of a HT treatment plan can be considerably deteriorated if an accurate patient setup system is not available. Taking into account the dose-response relations of the irradiated tumors and normal tissues, a radiobiological treatment plan evaluation can be performed, which may provide a closer association of the delivered treatment with the clinical outcome. In such situations, for effective evaluation and comparison of different treatment plans, traditional dose based evaluation tools can be complemented by the use of P+,BEUD diagrams.

Optimal treatment margins for radiotherapy of prostate cancer based on interfraction imaging

PURPOSE

To present a methodology to estimate optimal treatment margins for radiotherapy of prostate cancer based on interfraction imaging.

MATERIALS AND METHODS

Cone beam CT images of a prostate cancer patient undergoing fractionated radiotherapy were acquired at all treatment sessions. The clinical target volume (CTV) and organs at risk (OARs; bladder and rectum) were delineated in the images. Random sampling from the CTV-OAR library was performed in order to simulate fractionated radiotherapy including intra- and interpatient variability in setup and organ motion/deformation. For each simulated patient, four treatment fields defined by multileaf collimators were automatically generated around the planning CTV. The treatment margin (the distance from the CTV to the field border) was varied between 2.5 and 20 mm. Resulting dose distributions were calculated by a convolution method. Doses to OARs were reconstructed by polynomial warping, while the CTV was assumed to be a rigid body. The equivalent uniform dose (EUD), the tumor control probability (TCP) and the normal tissue complication probability (NTCP) were used to estimate the clinical effect. Patient repositioning strategies at treatment were compared.

RESULTS

The simulations produced population based EUD histograms for the CTV and the OARs. The number of patients receiving an optimal target EUD increased with increasing margins, but at the cost of an increasing number receiving a high EUD to the OARs. Calculations of the probability of complication-free tumor control and subsequent analysis gave an optimal treatment margin of about 10mm for the simulated population, if no correction strategy was undertaken.

CONCLUSIONS

The current work illustrates the principle of optimal treatment margins based on interfraction imaging. Clinically applicable margins may be obtained if a large patient image database is available.

Variation in radiation sensitivity and repair kinetics in different parts of the spinal cord

BACKGROUND

The spinal cord, known for its strongly serial character and high sensitivity to radiation even when a small segment is irradiated, is one of the most critical organs at risk to be spared during radiation therapy. To compare the sensitivity of different parts of the spinal cord, data for radiation myelopathy have been used.

MATERIAL AND METHODS

In the present study, the relative seriality model was fitted to two different datasets of clinical radiation myelitis concerning cervical spinal cord after treating 248 patients for head and neck cancer and thoracic spinal cord after treating 43 patients with lung carcinoma. The maximum likelihood method was applied to fit the clinical data. The model parameters and their 68% confidence intervals were calculated for each dataset. The alpha/beta ratio for the thoracic cord was also was also found to be 0.9 (0-3.0) Gy.

RESULTS

The dose-response curve for the more sensitive cervical myelopathy is well described by the parameters D(50)=55.9 (54.8-57.1) Gy, gamma=6.9 (5.0-9.2), s=0.13 (0.07-0.24), whereas the thoracic myelopathy is described by the parameters D(50)=75.5 (70.5-80.8) Gy, gamma=1.1 (0.6-1.6), s=36 (3.3-infinity).

DISCUSSION AND CONCLUSIONS

Large differences in radiation response between the cervical and thoracic region of spinal cord are thus observed: cervical myelopathy seems to be characterized by medium seriality, while thoracic spinal cord is characterized by a highly serial dose-response. The much steeper dose-response curve for cervical spinal cord myelopathy can be interpreted as a higher number of functional subunits consistent with a higher amount of white matter close to the brain.

Assessing four-dimensional radiotherapy planning and respiratory motion-induced dose difference based on biologically effective uniform dose

Four-dimensional (4D) radiotherapy is considered as a feasible and ideal solution to accommodate intra-fractional respiratory motion during conformal radiation therapy. With explicit inclusion of the temporal changes in anatomy during the imaging, planning, and delivery of radiotherapy, 4D treatment planning in principle provides better dose conformity. However, the clinical benefits of developing 4D treatment plans in terms of tumor control rate and normal tissue complication probability as compared to other treatment plans based on CT images of a fixed respiratory phase remains mostly unproven. The aim of our study is to comprehensively evaluate 4D treatment planning for nine lung tumor cases with both physical and biological measures using biologically effective uniform dose (D =) together with complication-free tumor control probability, P+. Based on the examined lung cancer patients and PTV margin applied, we found similar but not identical curves of DVH, and slightly different mean doses in tumor (up to 1.5%) and normal tissue in all cases when comparing 4D, P0%, and P50% plans. When it comes to biological evaluations, we did not observe definitively PTV size dependence in P+ among these nine lung cancer patients with various sizes of PTV. Moreover, it is not necessary that 4D plans would have better target coverage or higher P+ as compared to a fixed phase IMRT plan. However, on the contrary to significant deviations in P+ (up to 14.7%) observed if delivering the IMRT plan made at end-inhalation incorrectly at end-exhalation phase, we estimated the overall P+, PB, and PI for 4D composite plans that have accounted for intra-fractional respiratory motion.