Seminal Vesicle Treatment for Localized Prostate Cancer Treated with External Beam Radiotherapy.

This study retrospectively reviewed data from men with localized prostate cancer treated with external beam radiotherapy (EBRT). We identified 359 men with localized prostate cancer treated with curative EBRT at the Cross Cancer Institute between 2010-2011. The volume of seminal vesicles (SVs) treated as well as dose values were extracted. These volumes were compared to gold standard contours drawn by a trained expert based on consensus European Society for Radiotherapy and Oncology (ESTRO) contouring guidelines. Patient and tumor characteristics were extracted for these patients. Memorial Sloan Kettering prostate cancer nomogram was used to assign a predicted risk of SV involvement for each patient based on baseline tumor characteristics. In patients with a predicted risk of SV involvement greater than 15% (n = 184), 86.5% (SD = 18.6) of the base of the SVs were treated with EBRT, compared to 66.7% (SD = 32.6) for patients with a predicted risk of SV involvement less than 15% (n = 175, p < 0.0001). Similarly, the mean percentage of proximal and total SV volumes treated with EBRT was 75.6% (SD = 24.4) and 68.7% (SD = 26.0) for patients with a predicted risk of SV involvement of greater than 15%, compared to 50.3% (SD = 31.0, p < 0.0001) and 41.0% (SD = 27.8, p < 0.0001) for patients with a risk of less than 15%. The results indicate that all parts of the SVs are more likely to be contoured in men with >15% risk of SV involvement than those with <15% risk. However, radiation oncologists still contour a high percentage of SVs in men with <15% risk of SV involvement, suggesting that there may be over-treatment of SVs that increases the risk of rectal or bladder toxicity.

Tumor Volume Predicts for Pathologic Complete Response in Rectal Cancer Patients Treated With Neoadjuvant Chemoradiation.

OBJECTIVES: Nonoperative management (NOM) of locally advanced rectal cancer is an emerging approach allowing patients to preserve their anal sphincter. Identifying clinical factors associated with pathologic complete response (pCR) is essential for physicians and patients considering NOM. MATERIALS AND METHODS: In total, 412 locally advanced rectal cancer patients were included in this retrospective analysis. Tumor volumes were derived from pretreatment MRI. Clinical parameters such as tumor volume, stage, and location were analyzed by univariate and multivariate analysis, against pCR. A receiver operator characteristic curve was generated to identify a tumor volume cut-off with the highest clinically relevant Youden index for predicting pCR. RESULTS: Seventy-five of 412 patients (18%) achieved pCR. A tumor volume threshold of 37.3 cm 3 was identified as predictive for pCR. On regression analysis, a tumor volume >37.3 cm 3 was associated with a greater than 78% probability of not achieving pCR. On multivariate analysis, a GTV <37.3 cm 3 [odds ratio (OR)=3.7, P <0.0001] was significantly associated with an increased pCR rate, whereas tumor length > 4.85 cm was associated with pCR on univariate (OR=3.03, P <0.01) but not on multivariate analysis (OR=1.45, P =0.261). Other clinical parameters did not impact pCR rates. CONCLUSIONS: A tumor volume threshold of 37.3 cm 3 was identified as predictive for pCR in locally advanced rectal cancer patients receiving neoadjuvant chemoradiation. Tumors above this volume threshold corresponded to a greater than 78% probability of not achieving pCR. This information will be helpful at diagnosis for clinicians who are considering potential candidates for NOM.

Uncertainties Associated with Clonogenic Assays using a Cs-137 Irradiator and Ir-192 Afterloader: A Comprehensive Compilation for Radiation Researchers.

Clonogenic assays are the gold standard for measuring cell clonogenic survival and enable quantification of a cell line's radiosensitivity through the calculation of the surviving fraction, the ratio of cell clusters (colonies) formed after radiation exposure compared to the number formed without exposure. Such studies regularly utilize Cs-137 irradiators. While uncertainties for specific procedural aspects have been described previously, a comprehensive review has not been completed. We therefore quantified uncertainties associated with clonogenic assays performed using a Cs-137 Shepherd irradiator, and a recently established brachytherapy afterloader in vitro radiation delivery apparatus (BAIRDA), through a series of experiments and a literature review. The clonogenic assay is subject to uncertainties that affect the determination of the surviving fraction (e.g., accuracy of the number of cells seeded, potential effects of hypothermia, and the threshold number of cells for a cluster to be identified as a colony). Furthermore, dose delivery uncertainties related to both the Cs-137 irradiator and BAIRDA were also quantified. The combined standard (k = 1) uncertainty was ± 6.0% in the surviving fraction for the Cs-137 irradiator (±6.3% for BAIRDA), up to ± 2.2% in the dose delivered by the Cs-137 irradiator, and up to ± 4.3% in the dose delivered by BAIRDA. The largest individual uncertainties were associated with the number of cells seeded on a plate (3.4%) and inter-observer variability in counting (4.1%), suggesting that effective reduction of uncertainties in the conduct of the clonogenic assay may provide the greatest relief on the uncertainty budget. Finally, measurable impact on experimental findings was assessed by applying this uncertainty to clonogenic assays of SW756 cells using either a Cs-137 irradiator or BAIRDA, introducing a maximum shift in the reported radiobiological parameters α/ß and T1/2 of 0.3 Gy and 0.4 h, respectively, while the 95% confidence interval increased by 0.5 Gy and decreased by 0.4 h, respectively. Though the overall impact on radiobiological parameter estimation was small, the individual uncertainties could have a significant influence in other applications of in vitro experiments in radiation biology. Hence, better understanding of the uncertainties associated with both clonogenic assays and the radiation source used can improve the accuracy of experimental analysis and reproducibility of the results.

BAIRDA: a novelin vitrosetup to quantify radiobiological parameters for cervical cancer brachytherapy dose estimations.

Objective. Brachytherapy (BT) dose prescriptions for locally advanced cervical cancer are made with account for the radiobiological parameters,α/ßratio and halftime of repair (T1/2). However, a wide range of parameter values has been reported which can challenge commonly held equivalencies between dose prescriptions. This is the first reported study that aims to develop anin vitroexperimental technique using clinical high-dose-rate (HDR) and pulsed-dose-rate (PDR) Ir-192 brachytherapy afterloaders to quantify these parametersin vitroand to contextualize findings within contemporary practice.Approach. To efficiently quantifyα/ßandT1/2,in vitroexperiments more reflective of clinical BT practice than traditional clonogenic survival assays were developed and applied to four squamous cell carcinoma cell lines (CaSki, C-33A, SiHa, and SW756). Radiation was delivered using single acute and fractionated dose treatments with a conventional irradiator and clinical HDR and PDR BT afterloaders. For the latter, a novelbrachytherapyafterloaderin vitroradiationdeliveryapparatus (BAIRDA) was developed.Main Results. Theα/ßandT1/2values determined using BAIRDA and the conventional irradiator showed close agreement, validating the novel apparatus and technique. For CaSki, C-33A, SiHa, and SW756, the BAIRDA-measuredα/ßratios (5.2 [4.6-5.8], 5.6 [4.5-6.6], 6.3 [4.9-7.7], and 5.3 [4.7-6.0] Gy, respectively) were consistently smaller, while theT1/2(3.3 [2.7-3.9], 2.7 [2.0-3.3], 2.8 (2.4-3.1], and 4.8 [4.1-5.4] hours) larger, than the widely accepted values in clinical practice (α/ß= 10 Gy;T1/2 = 1.5 h).Significance.In vitroexperiments using BAIRDA provided evidence for differences between the conventionally selected and experimentally determinedα/ßratio andT1/2. Treatment regimens using HDR-BT and PDR-BT, designed to deliver equivalent radiobiological doses based on conventional values, were shown to differ by up to 27 Gy EQD2 - an effect that could impact treatment outcomes in cervical cancer. Furthermore, with BAIRDA, we have developed a novel method for radiobiological research in BT.

More Is Not Better When It Comes to Treating Rectal Cancer With Multimodal Chemoradiation Beyond the Standard Radiation Dose of 5040 cGy.

BACKGROUND: Radiation dose schedules for neoadjuvant chemoradiation for rectal cancers differ, with the most common dose schedule using 5040 cGy in 28 fractions. OBJECTIVES: The aim of this retrospective study was to assess the benefit of higher radiation doses beyond 5040 cGy in the context of pathological response and follow-up events. SETTING: The database from a provincial tertiary cancer center in Canada was the source of information for this study. PATIENTS: Included in this study were 508 consecutive patients with rectal cancer with locally advanced disease (clinical T3/T4 or N1/N2) who received neoadjuvant chemoradiation followed by surgery. Of the 508 patients, 281 received the standard radiation dose of 4500 to 5040 cGy and 227 received a dose >5040 cGy. MAIN OUTCOME MEASURE: The postsurgical pathology, late toxicities, and follow-up outcomes were analyzed. The outcomes were evaluated in relation to the dose of radiation received. RESULTS: Data regarding the clinical outcomes were comparable between the 4500 to 5040 cGy and >5040 cGy radiation groups with pathological complete response rates of 20.9% and 15.4% (p = 0.104); distant recurrence rates of 17.4% and 19.4% (p = 0.36); local recurrence rates of 3.2% and 3.5% (p = 0.36); and the median overall survival rates of 61 and 60.5 months (p = 0.8). No statistically significant correlation of improvement in outcomes was noted with radiation doses beyond 5040 cGy. LIMITATIONS: This is a retrospective study. CONCLUSION: Our study showed that dose escalation beyond the standard dose of 4500 to 5040cGy failed to achieve meaningful clinical outcomes. See Video Abstract at http://links.lww.com/DCR/B633. MS NO ES MEJOR CUANDO SE TRATA DE TRATAR EL CNCER DE RECTO CON QUIMIORRADIACIN MULTIMODAL MS ALL DE LA DOSIS DE RADIACIN ESTNDAR DE CGY: ANTECEDENTES:En neoadyuvancia de cáncer rectal es posible encontrar muchas variaciones, en radioterapia la dosis más común que usa 5040 cGy en 28 fracciones.OBJETIVOS:El objetivo de este estudio retrospectivo fue evaluar el beneficio de dosis de radiación más altas más allá de 5040cGy en el contexto de la respuesta patológica y en su seguimiento.AJUSTE:Base de datos de un centro de cáncer terciario provincial en Canadá.PACIENTES:Se incluyeron en este estudio quinientos ocho pacientes consecutivos con cáncer de recto y enfermedad localmente avanzada (clínica T3 / T4 o N1 / N2) que recibieron quimiorradiación neoadyuvante seguida de cirugía. De los 508 pacientes, 281 recibieron la dosis de radiación estándar de 4500-5040 cGy y 227 recibieron una dosis > 5040 cGy.PRINCIPAL MEDIDA DE RESULTADO:Se analizo evolucion posquirúrgica, toxicidad tardía y seguimiento. Los resultados se evaluaron en relación con la dosis de radiación recibida.RESULTADOS:Los datos con respecto a los resultados clínicos fueron comparables entre los grupos de radiación de 4500-5040 cGy y> 5040 cGy con tasas de respuesta patológica completa de 20,9% y 15,4% respectivamente (p = 0,104); tasas de recurrencia a distancia de 17,4% y 19,4%, respectivamente (p = 0,36); tasas de recurrencia local de 3,2% y 3,5%, respectivamente (p = 0,36); y la mediana de las tasas de supervivencia global de 61 y 60,5 meses, respectivamente (p = 0,8). No se observó una correlación estadísticamente significativa de mejoría en los resultados con dosis de radiación superiores a 5040 cGy.LIMITACIONES:Este es un estudio retrospectivo.CONCLUSIONES:Nuestro estudio mostró que el aumento de la dosis más allá de la dosis estándar de 4500-5040cGy no logró resultados clínicos significativos. Consulte Video Resumen en http://links.lww.com/DCR/B633. (Traducción-Dr. Gunther Bocic).

Disassembly of microtubules by intense terahertz pulses.

The biological effects of terahertz (THz) radiation have been observed across multiple levels of biological organization, however the sub-cellular mechanisms underlying the phenotypic changes remain to be elucidated. Filamentous protein complexes such as microtubules are essential cytoskeletal structures that regulate diverse biological functions, and these may be an important target for THz interactions underlying THz-induced effects observed at the cellular or tissue level. Here, we show disassembly of microtubules within minutes of exposure to extended trains of intense, picosecond-duration THz pulses. Further, the rate of disassembly depends on THz intensity and spectral content. As inhibition of microtubule dynamics is a mechanism of clinically-utilized anti-cancer agents, disruption of microtubule networks may indicate a potential therapeutic mechanism of intense THz pulses.

Dosimetric Parameters Predicting Late Small Bowel Toxicity in Patients With Rectal Cancer Receiving Neoadjuvant Chemoradiation.

PURPOSE: The aim of this study was to identify dosimetric parameters that predict late small bowel (SB) toxicity after neoadjuvant long course chemoradiation (CRT) for rectal cancer. METHODS AND MATERIALS: Four hundred eighty-six consecutive patients with locally advanced rectal cancers (clinical T3/T4 or N1/N2) who received CRT followed by surgery and had dosimetric data available for analysis were included in this study. The dose-volume relationship between small bowel irradiation and late small bowel toxicity was evaluated and a mathematical model to predict for late SB toxicity was derived. RESULTS: Among the 486 patients with a median follow-up of 60 months from completion of radiation, 36 (7.4%) patients experienced ≥ grade 2 and 21 (4.3%) developed ≥ grade 3 late SB toxicity. A statistically significant association between the development of grade ≥3 late small bowel toxicity and the volume of small bowel irradiated was found at each dose level from 5 to 40 Gy (P < .001 for all dose volumes) in 5 Gy intervals. The average SB volume for patients who experienced grade ≥2 SB toxicity was 2149.9 cm3 and the average SB volume for patients who experienced grade ≥3 SB toxicity was 2179.9 cm3. The predicted V30 for a 5% risk for grade ≥2 SB toxicity was 101.5 cm3 and for grade ≥3 SB toxicity was 201.5 cm3. The volume of small bowel receiving at least 30 Gy (V30) was most strongly associated with grade ≥3 SB toxicity. CONCLUSIONS: This study demonstrates the significant dose-volume relationship between volume of small bowel receiving 30 Gy (V30 Gy) and late grade ≥3 SB toxicity. When planning CRT for patients with rectal cancer, restricting V30 to <200 cm3 will be a useful guideline to minimize the 5 year grade ≥3 late SB toxicity to <5%.

Radiobiological dose calculation parameters for cervix cancer brachytherapy: A systematic review.

The GEC-ESTRO recommendation in cervical cancer treatment planning, including external beam radiotherapy and brachytherapy boosts, is to use radiobiological dose calculations. Such calculations utilize the linear-quadratic model to estimate the effect of multiple cellular response factors and dose delivery parameters. The radiobiological parameters utilized in these calculations are literature values estimated based on clinical and experimental results. However, the impact of the uncertainties associated with these parameters is often not fully appreciated. This review includes a summary of the radiobiological dose calculation (for both high-dose-rate and pulsed-dose-rate brachytherapy boost treatments) for cervical cancer and a compilation of the reported values of the associated parameters. As discrepancies exist between conventionally recommended and published values, equivalencies between current brachytherapy boosts may be imprecise and could create underappreciated uncertainties in the radiobiological dose calculations. This review highlights these uncertainties by calculating the radiobiological dose delivered by the brachytherapy boost when assuming different radiobiological parameter values (within the range reported by previous research). Furthermore, conventional treatment planning does not consider the effects of proliferation of the tumor over the treatment time, which can significantly decrease its radiobiological dose and can introduce an additional variance of over 7 Gy10. Further investigation of uncertainties in parameter values and modifications of current dose models could improve the accuracy of radiobiological dose calculation.

Cardiac-sparing radiation therapy using positioning breast shell for patients with left-sided breast cancer who are ineligible for breath-hold techniques.

PURPOSE: Patients with left-sided breast cancer (LSBC) are at increased risk of cardiac morbidity from adjuvant breast radiation therapy (ABRT). Breath-hold (BH) techniques substantially reduce the radiation received by heart during radiation therapy for LSBC. However, a subset of patients with LSBC are ineligible for BH techniques due to an inability to breath-hold or because of other comorbidities. To reduce radiation to the heart, we routinely use a custom-made breast shell for the treatment of patients with LSBC who are ineligible for BH techniques. This study evaluates the dosimetric impact of using a breast shell for patients with LSBC undergoing ABRT. METHODS AND MATERIALS: Sixteen consecutive patients with LSBC who failed BH and underwent ABRT using a breast shell during the period of 2014 to 2016 were identified. Treatment was planned using field-in-field tangents with a prescribed dose of 42.5 Gy in 16 fractions. Comparisons between plans with and without a shell were made for each patient using a paired t test to quantify the sparing of organs at risk (OARs) and target coverage. RESULTS: There was no statistically significant difference in the planning target volume of breast coverage. A statistically significant improvement was observed in sparing the heart, left ventricle (LV), and ipsilateral lung (P-value < .001). Plans with the shell spared OARs better than the no-shell plans with a mean dose of 2.15 Gy versus 5.15 Gy (58.2% reduction) to the heart, 3.27 Gy versus 9.00 Gy (63.7% reduction) to the LV, and 5.16 Gy versus 7.95 Gy (35% reduction) to the ipsilateral lung. The irradiated volumes of OARs for plans with and without shell are 13.3 cc versus 59.5 cc (77.6% reduction) for the heart, 6.2 cc versus 33.2 cc (81.2% reduction) for the LV, and 92.8 cc versus 162.5 cc (42.9% reduction) for the ipsilateral lung. CONCLUSIONS: A positioning breast shell offers significant benefit in terms of sparing the heart for patients with LSBC who are ineligible for BH techniques. It also can be used as a simple cardiac-sparing alternative in centers without BH capability.

Intensity-Modulated Radiotherapy (IMRT) vs Helical Tomotherapy (HT) in Concurrent Chemoradiotherapy (CRT) for Patients with Anal Canal Carcinoma (ACC): an analysis of dose distribution and toxicities.

PURPOSE: Intensity-modulated radiotherapy (IMRT) and helical tomotherapy (HT) have been adopted for radiotherapy treatment of anal canal carcinoma (ACC) due to better conformality, dose homogeneity and normal-tissue sparing compared to 3D-CRT. To date, only one published study compares dosimetric parameters of IMRT vs HT in ACC, but there are no published data comparing toxicities. Our objectives were to compare dosimetry and toxicities between these modalities. METHODS AND MATERIALS: This is a retrospective study of 35 ACC patients treated with radical chemoradiotherapy at two tertiary cancer institutions from 2008-2010. The use of IMRT vs HT was primarily based on center availability. The majority of patients received fluorouracil (5-FU) and 1-2 cycles of mitomycin C (MMC); 2 received 5-FU and cisplatin. Primary tumor and elective nodes were prescribed to ≥54Gy and ≥45Gy, respectively. Patients were grouped into two cohorts: IMRT vs HT. The primary endpoint was a dosimetric comparison between the cohorts; the secondary endpoint was comparison of toxicities. RESULTS: 18 patients were treated with IMRT and 17 with HT. Most IMRT patients received 5-FU and 1 MMC cycle, while most HT patients received 2 MMC cycles (p<0.01), based on center policy. HT achieved more homogenous coverage of the primary tumor (HT homogeneity and uniformity index 0.14 and 1.02 vs 0.29 and 1.06 for IMRT, p=0.01 and p<0.01). Elective nodal coverage did not differ. IMRT achieved better bladder, femoral head and peritoneal space sparing (V30 and V40, p ≤ 0.01), and lower mean skin dose (p<0.01). HT delivered lower bone marrow (V10, p<0.01) and external genitalia dose (V20 and V30, p<0.01). Grade 2+ hematological and non-hematological toxicities were similar. Febrile neutropenia and unscheduled treatment breaks did not differ (both p=0.13), nor did 3-year overall and disease-free survival (p=0.13, p=0.68). CONCLUSIONS: Chemoradiotherapy treatment of ACC using IMRT vs HT results in differences in dose homogenity and normal-tissue sparing, but no significant differences in toxicities.

Spatial and temporal distribution of γH2AX fluorescence in human cell cultures following synchrotron-generated X-ray microbeams: lack of correlation between persistent γH2AX foci and apoptosis.

Formation of γH2AX foci (a marker of DNA double-strand breaks), rates of foci clearance and apoptosis were investigated in cultured normal human fibroblasts and p53 wild-type malignant glioma cells after exposure to high-dose synchrotron-generated microbeams. Doses up to 283 Gy were delivered using beam geometries that included a microbeam array (50 µm wide, 400 µm spacing), single microbeams (60-570 µm wide) and a broad beam (32 mm wide). The two cell types exhibited similar trends with respect to the initial formation and time-dependent clearance of γH2AX foci after irradiation. High levels of γH2AX foci persisted as late as 72 h post-irradiation in the majority of cells within cultures of both cell types. Levels of persistent foci after irradiation via the 570 µm microbeam or broad beam were higher when compared with those observed after exposure to the 60 µm microbeam or microbeam array. Despite persistence of γH2AX foci, these irradiation conditions triggered apoptosis in only a small proportion (<5%) of cells within cultures of both cell types. These results contribute to the understanding of the fundamental biological consequences of high-dose microbeam irradiations, and implicate the importance of non-apoptotic responses such as p53-mediated growth arrest (premature senescence).

Modulation factors calculated with an EPID-derived MLC fluence model to streamline IMRT/VMAT second checks.

This work outlines the development of a robust method of calculating modulation factors used for the independent verification of MUs for IMRT and VMAT treatments, to replace onerous ion chamber measurements. Two-dimensional fluence maps were calculated for dynamic MLC fields that include MLC interleaf leakage, transmission, and tongue-and-groove effects, as characterized from EPID-acquired images. Monte Carlo-generated dose kernels were then used to calculate doses for a modulated field and that field with the modulation removed at a depth specific to the calculation point in the patient using in-house written software, Mod_Calc. The ratio of these two doses was taken to calculate modulation factors. Comparison between Mod_Calc calculation and ion chamber measurement of modulation factors for 121 IMRT fields yielded excellent agreement, where the mean difference between the two was -0.3% ± 1.2%. This validated use of Mod_Calc clinically. Analysis of 5,271 dynamic fields from clinical use of Mod_Calc gave a mean difference of 0.3% ± 1.0% between Mod_Calc and Eclipse-generated factors. In addition, 99.3% and 96.5% fields pass 5% and 2% criteria, respectively, for agreement between these two predictions. The development and use of Mod_Calc at our clinic has considerably streamlined our QA process for IMRT and RapidArc fields, compared to our previous method based on ion chamber measurements. As a result, it has made it feasible to maintain our established and trusted current in-house method of MU verification, without resorting to commercial software alternatives.

Evaluation of dose-volume metrics for microbeam radiation therapy dose distributions in head phantoms of various sizes using Monte Carlo simulations.

This work evaluates four dose-volume metrics applied to microbeam radiation therapy (MRT) using simulated dosimetric data as input. We seek to improve upon the most frequently used MRT metric, the peak-to-valley dose ratio (PVDR), by analyzing MRT dose distributions from a more volumetric perspective. Monte Carlo simulations were used to calculate dose distributions in three cubic head phantoms: a 2 cm mouse head, an 8 cm cat head and a 16 cm dog head. The dose distribution was calculated for a 4 × 4 mm² microbeam array in each phantom, as well as a 16 × 16 mm² array in the 8 cm cat head, and a 32 × 32 mm² array in the 16 cm dog head. Microbeam widths of 25, 50 and 75 µm and center-to-center spacings of 100, 200 and 400 µm were considered. The metrics calculated for each simulation were the conventional PVDR, the peak-to-mean valley dose ratio (PMVDR), the mean dose and the percentage volume below a threshold dose. The PVDR ranged between 3 and 230 for the 2 cm mouse phantom, and between 2 and 186 for the 16 cm dog phantom depending on geometry. The corresponding ranges for the PMVDR were much smaller, being 2-49 (mouse) and 2-46 (dog), and showed a slightly weaker dependence on phantom size and array size. The ratio of the PMVDR to the PVDR varied from 0.21 to 0.79 for the different collimation configurations, indicating a difference between the geometric dependence on outcome that would be predicted by these two metrics. For unidirectional irradiation, the mean lesion dose was 102%, 79% and 42% of the mean skin dose for the 2 cm mouse, 8 cm cat and 16 cm dog head phantoms, respectively. However, the mean lesion dose recovered to 83% of the mean skin dose in the 16 cm dog phantom in intersecting cross-firing regions. The percentage volume below a 10% dose threshold was highly dependent on geometry, with ranges for the different collimation configurations of 2-87% and 33-96% for the 2 cm mouse and 16 cm dog heads, respectively. The results of this study illustrate that different dose-volume metrics exhibit different functional dependences on MRT geometry parameters, and suggest that reliance on the PVDR as a predictor of therapeutic outcome may be insufficient.

Population TCP estimators in case of heterogeneous irradiation: a new discussion of an old problem.

PURPOSE: To investigate the capacity of two phenomenological expressions to describe the population tumor response in case of a heterogeneous irradiation of the tumor. The generalization of the individual tumor control probability (TCP) models to include the case of a heterogeneous irradiation is a trivial problem. However, an analytical solution that results in a closed form population TCP formula for the heterogeneous case is, unfortunately, a very complex mathematical problem. Therefore we applied a numerical approach to the problem. METHOD: Pseudo-experimental data sets are constructed through the generation of dose distributions and population TCP data obtained by a numerical solution of a multi-dimensional integral over an individual TCP model. The capacity of the following two phenomenological - Poisson and equivalent uniform dose (EUD) based - TCP expressions: [Figure: see text] to describe the population tumor response in case of heterogeneous irradiation is investigated through their fitting to the psuedo-experimental data sets. RESULTS AND CONCLUSIONS. While both expressions produce statistically acceptable fits to the pseudo-experimental data within 2% TCP error band, the use of the second expression is preferable since it produces considerably better fits to the data sets.

Analytic investigation into effect of population heterogeneity on parameter ratio estimates.

PURPOSE: A homogeneous tumor control probability (TCP) model has previously been used to estimate the alpha/beta ratio for prostate cancer from clinical dose-response data. For the ratio to be meaningful, it must be assumed that parameter ratios are not sensitive to the type of tumor control model used. We investigated the validity of this assumption by deriving analytic relationships between the alpha/beta estimates from a homogeneous TCP model, ignoring interpatient heterogeneity, and those of the corresponding heterogeneous (population-averaged) model that incorporated heterogeneity. METHODS AND MATERIALS: The homogeneous and heterogeneous TCP models can both be written in terms of the geometric parameters D(50) and gamma(50). We show that the functional forms of these models are similar. This similarity was used to develop an expression relating the homogeneous and heterogeneous estimates for the alpha/beta ratio. The expression was verified numerically by generating pseudo-data from a TCP curve with known parameters and then using the homogeneous and heterogeneous TCP models to estimate the alpha/beta ratio for the pseudo-data. RESULTS: When the dominant form of interpatient heterogeneity is that of radiosensitivity, the homogeneous and heterogeneous alpha/beta estimates differ. This indicates that the presence of this heterogeneity affects the value of the alpha/beta ratio derived from analysis of TCP curves. CONCLUSIONS: The alpha/beta ratio estimated from clinical dose-response data is model dependent--a heterogeneous TCP model that accounts for heterogeneity in radiosensitivity will produce a greater alpha/beta estimate than that resulting from a homogeneous TCP model.

Fundamental form of a population TCP model in the limit of large heterogeneity.

A population tumor control probability (TCP) model for fractionated external beam radiotherapy, based on Poisson statistics and in the limit of large parameter heterogeneity, is studied. A reduction of a general eight-parameter TCP equation, which incorporates heterogeneity in parameters characterizing linear-quadratic radiosensitivity, repopulation, and clonogen number, to an equation with four parameters is obtained. The four parameters represent the mean and standard deviation for both clonogen number and a generalized radiosensitivity that includes linear-quadratic and repopulation descriptors. Further, owing to parameter inter-relationship, it is possible to express these four parameters as three ratios of parameters in the large heterogeneity limit. These ratios can be directly linked to two defining features of the TCP dose response: D50 and gamma50. In the general case, the TCP model can be written in terms of D50, gamma50 and a third parameter indicating the ratio of the levels of heterogeneity in clonogen number and generalized radiosensitivity; however, the third parameter is unnecessary when either of these two sources of heterogeneity is dominant. It is shown that heterogeneity in clonogen number will have little impact on the TCP formula for clinical scenarios, and thus it will generally be the case that the fundamental form of the Poisson-based population TCP model can be specified completely in terms of D50 and gamma50: TCP= 1/2 erfc[square root of pi(gamma50)(D50/D-1)]. This implies that limited radiobiological information can be determined by the analysis of dose response data: information about parameter ratios can be ascertained, but knowledge of absolute values for the fundamental radiobiological parameters will require independent auxiliary measurements.

A TCP-NTCP estimation module using DVHs and known radiobiological models and parameter sets.

Radiotherapy treatment plan evaluation relies on an implicit estimation of the tumor control probability (TCP) and normal tissue complication probability (NTCP) arising from a given dose distribution. A potential application of radiobiological modeling to radiotherapy is the ranking of treatment plans via a more explicit determination of TCP and NTCP values. Although the limited predictive capabilities of current radiobiological models prevent their use as a primary evaluative tool, radiobiological modeling predictions may still be a valuable complement to clinical experience. A convenient computational module has been developed for estimating the TCP and the NTCP arising from a dose distribution calculated by a treatment planning system, and characterized by differential (frequency) dose-volume histograms (DDVHs). The radiobiological models included in the module are sigmoidal dose response and Critical Volume NTCP models, a Poisson TCP model, and a TCP model incorporating radiobiological parameters describing linear-quadratic cell kill and repopulation. A number of sets of parameter values for the different models have been gathered in databases. The estimated parameters characterize the radiation response of several different normal tissues and tumor types. The system also allows input and storage of parameters by the user, which is particularly useful because of the rapidly increasing number of parameter estimates available in the literature. Potential applications of the system include the following: comparing radiobiological predictions of outcome for different treatment plans or types of treatment; comparing the number of observed outcomes for a cohort of patient DVHs to the predicted number of outcomes based on different models/parameter sets; and testing of the sensitivity of model predictions to uncertainties in the parameter values. The module thus helps to amalgamate and make more accessible current radiobiological modeling knowledge, and may serve as a useful aid in the prospective and retrospective analysis of radiotherapy treatment plans.