Search

For many years the linear-quadratic (LQ) model has been widely used to describe the effects of total dose and dose per fraction at low-to-intermediate doses in conventional fractionated radiotherapy. Recent advances in stereotactic radiosurgery (SRS) and stereotactic radiotherapy (SRT) have increased the interest in finding a reliable cell survival model, which will be accurate at high doses, as well. Different models have been proposed for improving descriptions of high dose survival responses, such as the Universal Survival Curve (USC), the Kavanagh-Newman (KN) and several generalizations of the LQ model, e.g. the Linear-Quadratic-Linear (LQL) model and the Padé Linear Quadratic (PLQ) model. The purpose of the present study is to compare a number of models in order to find the best option(s) which could successfully be used as a fractionation correction method in SRT. In this work, six independent experimental data sets were used: CHOAA8 (Chinese hamster fibroblast), H460 (non-small cell lung cancer, NSLC), NCI-H841 (small cell lung cancer, SCLC), CP3 and DU145 (human prostate carcinoma cell lines) and U1690 (SCLC). By detailed comparisons with these measurements, the performance of nine different radiobiological models was examined for the entire dose range, including high doses beyond the shoulder of the survival curves. Using the computed and measured cell surviving fractions, comparison of the goodness-of-fit for all the models was performed by means of the reduced χ2-test with a 95% confidence interval. The obtained results indicate that models with dose-independent final slopes and extrapolation numbers generally represent better choices for SRT. This is especially important at high doses where the final slope and extrapolation numbers are presently found to play a major role. The PLQ, USC and LQL models have the least number of shortcomings at all doses. The extrapolation numbers and final slopes of these models do not depend on dose. Their asymptotes for the cell surviving fractions are exponentials at low as well as high doses, and this is in agreement with the behaviour of the corresponding experimental data. This is an important improvement over the LQ model which predicts a Gaussian at high doses. Overall and for the highlighted reasons, it was concluded that the PLQ, USC and LQL models are theoretically well-founded. They could prove useful compared to the other proposed radiobiological models in clinical applications for obtaining uniformly accurate cell surviving fractions encountered in stereotactic high-dose radiotherapy as well as at medium and low doses.

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 aim of this study was to investigate what bowel organ and delivered dose levels are most relevant for the development of 'emptying of all stools into clothing without forewarning' so that the related dose-responses could be derived as an aid in avoiding this distressing symptom in the future.Of the 77 gynecological cancer survivors treated with radiotherapy (RT) for gynecological cancer, 13 developed the symptom. The survivors were treated between 1991 and 2003. The anal-sphincter region, the rectum, the sigmoid and the small intestines were all delineated and the dose-volume histograms were exported for each patient. The dose-volume parameters were estimated fitting the data to the Relative Seriality (RS), the Lyman and the generalized Equivalent Uniform Dose (gEUD) model.The dose-response parameters for all three models and four organs at risk (OARs) were estimated. The data from the sigmoid fits the studied models best: D50 was 58.8 and 59.5 Gy (RS, Lyman), γ50 was 1.60 and 1.57 (RS, Lyman), s was 0.32, n was 0.13 and a was 7.7 (RS, Lyman, gEUD). The estimated volume parameters indicate that the investigated OARs behave serially for this endpoint. Our results for the three models studied indicate that they have the same predictive power (similar LL values) for the symptom as a function of the dose for all investigated OARs.In our study, the anal-sphincter region and sigmoid fit our data best, but all OARs were found to have steep dose-responses for 'emptying of all stools into clothing without forewarning' and thus, the outcome can be predicted with an NTCP model. In addition, the dose to the four studied OARs may be considered when minimizing the risk of the symptom.

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 DP+ of 27.0% for a D effective uniform dose of 16.3 Gy. Similarly, the HT treatment plans are more effective than the 3D-conformal plans by an average DP+ of 23.8% for a D 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.

Determination of the dose-response relations for oesophageal stricture after radiotherapy of the head and neck.In this study 33 patients who developed oesophageal stricture and 39 patients as controls are included. The patients received radiation therapy for head and neck cancer at Karolinska University Hospital, Stockholm, Sweden. For each patient the 3D dose distribution delivered to the upper 5 cm of the oesophagus was analysed. The analysis was conducted for two periods, 1992-2000 and 2001-2005, due to the different irradiation techniques used. The fitting has been done using the relative seriality model.For the treatment period 1992-2005, the mean doses were 49.8 and 33.4 Gy, respectively, for the cases and the controls. For the period 1992-2000, the mean doses for the cases and the controls were 49.9 and 45.9 Gy and for the period 2001-2005 were 49.8 and 21.4 Gy. For the period 2001-2005 the best estimates of the dose-response parameters are D(50)=61.5 Gy (52.9-84.9 Gy), γ=1.4 (0.8-2.6) and s=0.1 (0.01-0.3).Radiation-induced strictures were found to have a dose response relation and volume dependence (low relative seriality) for the treatment period 2001-2005. However, no dose response relation was found for the complete material.

Radiation treatment of arteriovenous malformations (AVMs) has a slow and progressive vaso-occlusive effect. Some studies suggested the possible role of vascular structure in this process. A detailed biomathematical model has been used, where the morphological, biophysical and hemodynamic characteristics of intracranial AVM vessels are faithfully reproduced. The effect of radiation on plexiform and fistulous AVM nidus vessels was simulated using this theoretical model. The similarities between vascular and electrical networks were used to construct this biomathematical AVM model and provide an accurate rendering of transnidal and intranidal hemodynamics. The response of different vessels to radiation and their obliteration probability as a function of different angiostructures were simulated and total obliteration was defined as the probability of obliteration of all possible vascular pathways. The dose response of the whole AVM is observed to depend on the vascular structure of the intra-nidus AVM. Furthermore, a plexiform AVM appears to be more prone to obliteration compared with an AVM of the same size but having more arteriovenous fistulas. Finally, a binomial model was introduced, which considers the number of crucial vessels and is able to predict the dose response behavior of AVMs with a complex vascular structure.

During the last 20 years, the field of cellular and not least molecular radiation biology has been developed substantially and can today describe the response of heterogeneous tumors and organized normal tissues to radiation therapy quite well. An increased understanding of the sub-cellular and molecular response is leading to a more general systems biological approach to radiation therapy and treatment optimization. It is interesting that most of the characteristics of the tissue infrastructure, such as the vascular system and the degree of hypoxia, have to be considered to get an accurate description of tumor and normal tissue responses to ionizing radiation. In the limited space available, only a brief description of some of the most important concepts and processes is possible, starting from the key functional genomics pathways of the cell that are not only responsible for tumor development but also responsible for the response of the cells to radiation therapy. The key mechanisms for cellular damage and damage repair are described. It is further more discussed how these processes can be brought to inactivate the tumor without severely damaging surrounding normal tissues using suitable radiation modalities like intensity-modulated radiation therapy (IMRT) or light ions. The use of such methods may lead to a truly scientific approach to radiation therapy optimization, particularly when invivo predictive assays of radiation responsiveness becomes clinically available at a larger scale. Brief examples of the efficiency of IMRT are also given showing how sensitive normal tissues can be spared at the same time as highly curative doses are delivered to a tumor that is often radiation resistant and located near organs at risk. This new approach maximizes the probability to eradicate the tumor, while at the same time, adverse reactions in sensitive normal tissues are as far as possible minimized using IMRT with photons and light ions.

The investigation of the clinical efficacy and effectiveness of Intensity Modulated Radiotherapy (IMRT) using Multileaf Collimators (MLC) and Helical Tomotherapy (HT) has been an issue of increasing interest over the past few years. In order to assess the suitability of a treatment plan, dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean, and standard deviation of the dose distribution are typically used. Nevertheless, the radiobiological parameters of the different tumors and normal tissues are often not taken into account. The use of the biologically effective uniform dose (D̿) together with the complication-free tumor control probability ( P+) were applied to evaluate the two radiation modalities. Two different clinical cases of brain and cranio-spinal axis cancers have been investigated by developing a linac MLC-based step-and-shoot IMRT plan and a Helical Tomotherapy plan. The treatment plans of the MLC-based IMRT were developed on the Philips treatment planning station using the Pinnacle 7.6 software release while the dedicated Tomotherapy treatment planning station was used for the HT plan. With the use of the P+ index and the D̿ concept as the common prescription point, the different treatment plans were compared based on radiobiological measures. The tissue response probabilities were plotted against D̿ for a range of prescription doses. The applied plan evaluation method shows that in the brain cancer, the HT treatment gives slightly better results than the MLC-based IMRT in terms of optimum expected clinical outcome ( P+ of 66.1% and 63.5% for a D̿ to the PTV of 63.0 Gy and 62.0 Gy, respectively). In the cranio-spinal axis cancer, the HT plan is significantly better compared to the MLC-based IMRT plan over the clinically useful dose prescription range ( P+ of 84.1% and 28.3% for a D̿ to the PTV of 50.6 Gy and 44.0 Gy, respectively). If a higher than 5% risk for complications could be allowed, the complication-free tumor control could be increased by almost 30% compared to the initial dose prescription. In comparison to MLC based-IMRT, HT can better encompass the often large PTV while minimizing the volume of the OARs receiving high dose. 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 tumors and normal tissues. The use of P — D̿ diagrams can complement the traditional tools of evaluation such as DVHs, in order to compare and effectively evaluate different treatment plans.

The full potential of biologically optimized radiation therapy can only be maximized with the prediction of individual patient radiosensitivity prior to treatment. Unfortunately, the available biological parameters, derived from clinical trials, reflect an average radiosensitivity of the examined populations. In the present study, a breast cancer patient of stage I-II with positive lymph nodes was chosen in order to analyse the effect of the variation of individual radiosensitivity on the optimal dose distribution. Thus, deviations from the average biological parameters, describing tumour, heart and lung response, were introduced covering the range of patient radiosensitivity reported in the literature. Two treatment configurations of three and seven biologically optimized intensity-modulated beams were employed. The different dose distributions were analysed using biological and physical parameters such as the complication-free tumour control probability (P(+)), the biologically effective uniform dose (D), dose volume histograms, mean doses, standard deviations, maximum and minimum doses. In the three-beam plan, the difference in P(+) between the optimal dose distribution (when the individual patient radiosensitivity is known) and the reference dose distribution, which is optimal for the average patient biology, ranges up to 13.9% when varying the radiosensitivity of the target volume, up to 0.9% when varying the radiosensitivity of the heart and up to 1.3% when varying the radiosensitivity of the lung. Similarly, in the seven-beam plan, the differences in P(+) are up to 13.1% for the target, up to 1.6% for the heart and up to 0.9% for the left lung. When the radiosensitivity of the most important tissues in breast cancer radiation therapy was simultaneously changed, the maximum gain in outcome was as high as 7.7%. The impact of the dose-response uncertainties on the treatment outcome was clinically insignificant for the majority of the simulated patients. However, the jump from generalized to individualized radiation therapy may significantly increase the therapeutic window for patients with extreme radio sensitivity or radioresistance, provided that these are identified. Even for radiosensitive patients a simple treatment technique is sufficient to maximize the outcome, since no significant benefits were obtained with a more complex technique using seven intensity-modulated beams portals.

Haptic technology enables us to incorporate the sense of touch into computer applications, providing an additional input/ output channel. The purpose of this study was to examine if haptic feedback can help physicians and other practitioners to interact with medical imaging and treatment planning systems. A haptic application for outlining target areas (a key task in radiation therapy treatment planning) was implemented and then evaluated via a controlled experiment with ten subjects. Even though the sample size was small, and the application only a prototype, results showed that haptic feedback can significantly increase (p0.05) the speed of outlining target volumes and organs at risk. No significant differences were found regarding precision or perceived usability. This promising result warrants further development of a full haptic application for this task. Improvements to the usability of the application as well as to the forces generated have been implemented and an experiment with more subjects is planned.

This work aims to evaluate the predictive strength of the relative seriality, parallel and Lyman-Kutcher-Burman (LKB) normal tissue complication probability (NTCP) models regarding the incidence of radiation pneumonitis (RP), in a group of patients following lung cancer radiotherapy and also to examine their correlation with pulmonary function tests (PFTs). The study was based on 47 patients who received radiation therapy for stage III non-small-cell lung cancer. For each patient, lung dose volume histograms (DVHs) and the clinical treatment outcome were available. Clinical symptoms, radiological findings and pulmonary function tests incorporated in a post-treatment follow-up period of 18 months were used to assess the manifestation of radiation induced complications. Thirteen of the 47 patients were scored as having radiation induced pneumonitis, with RTOG criteria grade 3 and 28 of the 47 with RTOG criteria grade 2. Using this material, different methods of estimating the likelihood of radiation effects were evaluated, by analysing patient data based on their full dose distributions and associating the calculated complication rates with the clinical follow-up records. Lungs were evaluated as a paired organ as well as individual lungs. Of the NTCP models examined in the overall group considering the dose distribution in the ipsilateral lung, all models were able to predict radiation induced pneumonitis only in the case of grade 2 radiation pneumonitis score, with the LKB model giving the best results (chi2-test: probability of agreement between the observed and predicted results Pchi(chi2)=0.524 using the 0.05 significance level). The NTCP modelling considering lungs as a paired organ did not give statistically acceptable results. In the case of lung cancer radiotherapy, the application of different published radiobiological parameters alters the NTCP results, but not excessively as in the case of breast cancer radiotherapy. In this relatively small group of lung cancer patients, no positive statistical correlation could be established between the incidence of radiation pneumonitis as estimated by NTCP models and the pulmonary function test evaluation. However, the use of PFTs as markers or predictors for the incidence or severity of radiation induced pneumonitis must be investigated further.

Background and purpose To investigate the dose-response relation between the dose to the vagina and the patient-reported symptom ‘absence of vaginal elasticity’ and how time to follow-up influences this relation. Materials and methods The study included 78 long-term gynecological cancer survivors treated between 1991 and 2003 with external beam radiation therapy. Of those, 24 experienced absence of vaginal elasticity. A normal tissue complication model is introduced that takes into account the influence of time to follow-up on the dose-response relation and the patient’s age. The best estimates of the dose-response parameters were calculated using Probit, Probit-Relative Seriality (RS) and Probit-time models. Log likelihood (LL) values and the Akaike Information Criterion (AIC) were used to evaluate the model fit. Results The dose-response parameters for ‘absence of vaginal elasticity’ according to the Probit and Probit-time models with the 68% Confidence Intervals (CI) were: LL=−39.8, D50 = 49.7 (47.2–52.4) Gy, γ 50 = 1.40 (1.12–1.70) and LL = −37.4, D50 = 46.9 (43.5–50.9) Gy, γ 50 = 1.81(1.17–2.51) respectively. Conclusions The proposed model, which describes the influence of time to follow-up on the dose-response relation, fits our data best. Our data indicate that the steepness of the dose-response curve of the dose to the vagina and the symptom ‘absence of vaginal elasticity’ increases with time to follow-up, while D50decreases. • Alevronta E, Avall-Lundqvist E, Al-Abany M, et al. Time-dependent dose-response relation for absence of vaginal elasticity after gynecological radiation therapy. Radiother Oncol 2016. http://dx.doi.org/10.1016/j.radonc.2016.02.013 .

The choice of the appropriate model and parameter set in determining the relation between the incidence of radiation pneumonitis and dose distribution in the lung is of great importance, especially in the case of breast radiotherapy where the observed incidence is fairly low. From our previous study based on 150 breast cancer patients, where the fits of dose-volume models to clinical data were estimated (Tsougos et al 2005 Evaluation of dose-response models and parameters predicting radiation induced pneumonitis using clinical data from breast cancer radiotherapy Phys. Med. Biol. 50 3535-54), one could get the impression that the relative seriality is significantly better than the LKB NTCP model. However, the estimation of the different NTCP models was based on their goodness-of-fit on clinical data, using various sets of published parameters from other groups, and this fact may provisionally justify the results. Hence, we sought to investigate further the LKB model, by applying different published parameter sets for the very same group of patients, in order to be able to compare the results. It was shown that, depending on the parameter set applied, the LKB model is able to predict the incidence of radiation pneumonitis with acceptable accuracy, especially when implemented on a sub-group of patients (120) receiving [see text]|EUD higher than 8 Gy. In conclusion, the goodness-of-fit of a certain radiobiological model on a given clinical case is closely related to the selection of the proper scoring criteria and parameter set as well as to the compatibility of the clinical case from which the data were derived.

Comprehensive atlases of lymph node topography are necessary tools to provide a detailed description of the lymphatic distribution in relation to other organs and structures. Despite the recent developments of atlases and guidelines focusing on definitions of lymphatic regions, a comprehensive and detailed description of the three-dimensional (3D) nodal distribution is lacking. This article describes a new 3D atlas of lymph node topography based on the digital images of the Visible Human Male Anatomical (VHMA) data set. About 1,200 lymph nodes were localized in the data set and their distribution was compared with data from current cross-sectional lymphatic atlases. The identified nodes were delineated and then labeled with different colors that corresponded to their anatomical locations. A series of 2D illustrations, showing discrete locations, description, and distribution of major lymph nodes, was compiled to form a cross-sectional atlas. The resultant contours of all localized nodes in the VHMA data set were superimposed to develop a volumetric model. A 3D reconstruction was generated for the lymph nodes and surrounding structures. The volumetric lymph node topography was also integrated into the existing VOXEL-MAN digital atlas to obtain an interactive and photo-realistic visualization of the lymph nodes showing their proximity to blood vessels and surrounding organs. The lymph node topography forms part of our whole body atlas database, which includes organs, definitions, and parameters that are related to radiation therapy. The lymph node topography atlas could be utilized for visualization and exploration of the 3D lymphatic distribution to assist in defining the target volume for treatment based on the lymphatic spread surrounding the primary tumor.

The aim was to investigate and compare the influence of linear energy transfer (LET), dose and time on the induction of apoptosis in a human melanoma cell line exposed to accelerated light boron ((10)B) ions and photons. Cells were exposed in vitro to doses up to 6 Gy accelerated boron ions (40, 80, 125 and 160 eV nm(-1)) and up to 12 Gy photons (0.2 eV nm(-1)). The induction of apoptosis was measured up to 9 days after irradiation using morphological characterization of apoptotic cells and bodies. In parallel, measurements of cell-cycle distribution, monitored by DNA flow cytometry, and cell survival based on the clonogenic cell survival assay, were performed. In addition, the induction and repair of DNA double-strand breaks (DSB), using pulsed-field gel electrophoresis (PFGE) were studied. Accelerated boron ions induced a significant increase in apoptosis as compared with photons at all time points studied. At 1-5 h the percentage of radiation-induced apoptotic cells increased with both dose and LET. At the later time points (24-216 h) the apoptotic response was more complex and did not increase in a strictly LET-dependent manner. The early premitotic apoptotic cells disappeared at 24 h following exposure to the highest LET (160 eV nm(-1)). A postmitotic apoptotic response was seen after release of the dose-, time- and LET-dependent G2/M accumulations. The loss of clonogenic ability was dose- and LET-dependent and the fraction of un-rejoined DSB increased with increasing LET. Despite the LET-dependent clonogenic cell killing, it was not possible to measure quantitatively a LET-dependent apoptotic response. This was due to the different time course of appearance and disappearance of apoptotic cells.

The purpose of this work is to provide some statistical methods for evaluating the predictive strength of radiobiological models and the validity of dose-response parameters for tumour control and normal tissue complications. This is accomplished by associating the expected complication rates, which are calculated using different models, with the clinical follow-up records. These methods are applied to 77 patients who received radiation treatment for head and neck cancer and 85 patients who were treated for arteriovenous malformation (AVM). The three-dimensional dose distribution delivered to esophagus and AVM nidus and the clinical follow-up results were available for each patient. Dose-response parameters derived by a maximum likelihood fitting were used as a reference to evaluate their compatibility with the examined treatment methodologies. The impact of the parameter uncertainties on the dose-response curves is demonstrated. The clinical utilization of the radiobiological parameters is illustrated. The radiobiological models (relative seriality and linear Poisson) and the reference parameters are validated to prove their suitability in reproducing the treatment outcome pattern of the patient material studied (through the probability of finding a worse fit, area under the ROC curve and chi2 test). The analysis was carried out for the upper 5 cm of the esophagus (proximal esophagus) where all the strictures are formed, and the total volume of AVM. The estimated confidence intervals of the dose-response curves appear to have a significant supporting role on their clinical implementation and use.

The advent of intensity-modulated radiation therapy makes it increasingly important to model the response accurately when large volumes of normal tissues are irradiated by controlled graded dose distributions aimed at maximizing tumor cure and minimizing normal tissue toxicity. The cell survival model proposed here is very useful and flexible for accurate description of the response of healthy tissues as well as tumors in classical and truly radiobiologically optimized radiation therapy. The repairable-conditionally repairable (RCR) model distinguishes between two different types of damage, namely the potentially repairable, which may also be lethal, i.e. if unrepaired or misrepaired, and the conditionally repairable, which may be repaired or may lead to apoptosis if it has not been repaired correctly. When potentially repairable damage is being repaired, for example by nonhomologous end joining, conditionally repairable damage may require in addition a high-fidelity correction by homologous repair. The induction of both types of damage is assumed to be described by Poisson statistics. The resultant cell survival expression has the unique ability to fit most experimental data well at low doses (the initial hypersensitive range), intermediate doses (on the shoulder of the survival curve), and high doses (on the quasi-exponential region of the survival curve). The complete Poisson expression can be approximated well by a simple bi-exponential cell survival expression, S(D) = e(-aD) + bDe(-cD), where the first term describes the survival of undamaged cells and the last term represents survival after complete repair of sublethal damage. The bi-exponential expression makes it easy to derive D(0), D(q), n and alpha, beta values to facilitate comparison with classical cell survival models.

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

A generalization of the standard dose-response gradient to arbitrarily heterogeneous dose distributions has been developed. The generalized dose-response gradient is the scalar product of the vector representing the dose distribution and the gradient of the dose-response relation with respect to that dose vector. It is shown that, for a tumor, the individual &#110 -values for each portion of the tumor divided by the corresponding local tumor control probability should be added to get the total value for the heterogeneously irradiated tumor. This corresponds to summing up the contributions of all tumor volumes so that the total value of the gradient is related to the logarithm of the total tumor clonogen number. General expressions are also derived for the change in the dose-response relation as a function of a change in the delivered dose distribution.

The objective of the dynamic radiotherapy project 'Dynarad' within the European Community has been to compare and grade treatment techniques that are currently applied or being developed at the participating institutions. Cervical cancer was selected as the tumour site on the grounds that the involved organs at risk, mainly the rectum and the bladder, are very close to the tumour and partly located inside the internal target volume. In this work, a solid phantom simulating the pelvic anatomy was used by institutions in Belgium, France, Greece, Holland, Italy, Sweden and the United Kingdom. The results were evaluated using both biological and physical criteria. The main purpose of this parallel evaluation is to test the value of biological and physical evaluations in comparing treatment techniques. It is demonstrated that the biological objective functions allow a much higher conformality and a more clinically relevant scoring of the outcome. Often external beam treatment techniques have to be combined with intracavitary therapy to give clinically acceptable results. However, recent developments can reduce or even eliminate this need by delivering more conformal dose distributions using intensity modulated external dose delivery. In these cases the reliability of the patient set-up procedure becomes critical for the effectiveness of the treatment.

Purpose : (a) To assess the increased risk of death due to ischemic heart disease (IHD) in a group of patients treated for Hodgkin's disease (HD) with radiation therapy (RT) as the primary treatment. (b) To quantify the dose response of IHD using a biophysical model. Materials and methods : Patient material consisted of 157 patients diagnosed for HD between 1972 and 1985 who received RT as the primary treatment at Radiumhemmet, Karolinska Hospital. The general population formed the control group. The RT treatments were reconstructed based on the individual treatment data and simulator films. Individual clinical and dosimetrical data were analyzed with the relative seriality model. The material was also analyzed grouping the material according to dose-volume constraints. Results : Of the 157 patients, 13 (8.3%) died due to IHD. The standardized mortality ratio (SMR) was 5.0 (95% CI, 2.7–8.6). Analysis of dose-volume histograms (DVH) showed an increasing risk with increasing dose to a larger volume fraction. The observed individual clinical complication data could not be modeled unambiguously. The group analysis resulted in the dose-response parameters: D 50 =71 Gy, γ =0.96 and s =1.0. Conclusions : A significantly increased risk of death due to IHD following RT for HD was found. The risk was found to increase with higher dose and larger volume fraction irradiated.

The aim of the present work is to better account for the known uncertainties in radiobiological response parameters when optimizing radiation therapy. The radiation sensitivity of a specific patient is usually unknown beyond the expectation value and possibly the standard deviation that may be derived from studies on groups of patients. Instead of trying to find the treatment with the highest possible probability of a desirable outcome for a patient of average sensitivity, it is more desirable to maximize the expectation value of the probability for the desirable outcome over the possible range of variation of the radiation sensitivity of the patient. Such a stochastic optimization will also have to consider the distribution function of the radiation sensitivity and the larger steepness of the response for the individual patient. The results of stochastic optimization are also compared with simpler methods such as using biological response `margins' to account for the range of sensitivity variation. By using stochastic optimization, the absolute gain will typically be of the order of a few per cent and the relative improvement compared with non-stochastic optimization is generally less than about 10 per cent. The extent of this gain varies with the level of interpatient variability as well as with the difficulty and complexity of the case studied. Although the dose changes are rather small (

A new general beam optimization algorithm for inverse treatment planning is presented. It utilizes a new formulation of the probability to achieve complication-free tumour control. The new formulation explicitly describes the dependence of the treatment outcome on the incident fluence distribution, the patient geometry, the radiobiological properties of the patient and the fractionation schedule. In order to account for both measured and non-measured positioning uncertainties, the algorithm is based on a combination of dynamic and stochastic optimization techniques. Because of the difficulty in measuring all aspects of the intra- and interfractional variations in the patient geometry, such as internal organ displacements and deformations, these uncertainties are primarily accounted for in the treatment planning process by intensity modulation using stochastic optimization. The information about the deviations from the nominal fluence profiles and the nominal position of the patient relative to the beam that is obtained by portal imaging during treatment delivery, is used in a feedback loop to automatically adjust the profiles and the location of the patient for all subsequent treatments. Based on the treatment delivered in previous fractions, the algorithm furnishes optimal corrections for the remaining dose delivery both with regard to the fluence profile and its position relative to the patient. By dynamically refining the beam configuration from fraction to fraction, the algorithm generates an optimal sequence of treatments that very effectively reduces the influence of systematic and random set-up uncertainties to minimize and almost eliminate their overall effect on the treatment. Computer simulations have shown that the present algorithm leads to a significant increase in the probability of uncomplicated tumour control compared with the simple classical approach of adding fixed set-up margins to the internal target volume.

Sublethal damage after radiation exposure may become lethal or be repaired according to repair kinetics. This is a well-established concept in conventional radiotherapy. It also plays an important role in single-dose stereotactic radiotherapy treatments, often called stereotactic radiosurgery, when duration of treatment is extended due to source decay or treatment planning protocol. The purpose of this study is to look into the radiobiological characteristics of normal brain tissue and treatment protocols and find a way to optimize the time course of these protocols. The general problem is nonlinear and can be solved numerically. For numerical optimization of the time course of radiation protocol, a biexponential repair model with slow and fast components was considered. With the clinically imposed constraints of a fixed total dose and total treatment time, three parameters for each fraction (dose-rate, fraction duration, time of each fraction) were simultaneously optimized. A biological optimization can be performed by maximizing the therapeutic difference between tumor control probability and normal tissue complication probability. Specifically, for gamma knife radiosurgery, this approach can be implemented for normal brain tissue or tumor voxels separately in a treatment plan. Differences in repair kinetics of normal tissue and tumors can be used to find clinically optimized protocols. Thus, in addition to considering the physical dose in tumor and normal tissue, we also account for repair of sublethal damage in both these tissues.

An efficient algorithm has been developed to estimate the relative biological effectiveness (RBE) from experimental RBE data and measured single-event energy deposition spectra. As benchmarks for the calculations well established RBE data for crypt cell survival and associated microdosimetric distributions from a number of therapeutic high- and low-LET beams was used. As a by-product of the calculations the RBE for this system was determined for the lineal energy range 0.1-5000 keV microns-1 in very good agreement with experimental data. These data have been applied to estimate the biological effectiveness of therapeutic high-energy bremsstrahlung beams characterized by microdosimetric measurements with a wall-less proportional counter. The absorbed dose component due to photoneutrons and charged particles from photonuclear reactions in scanned 50 MV bremsstrahlung beams was measured to be about 2% of the total absorbed dose. The RBE of 50 MV scanned bremsstrahlung beams was estimated to be 1.09 for jejunum crypt cell survival as the biological endpoint at absorbed doses of the order of 10 Gy, in fair agreement with reported results based on radiobiological experiments. The RBE is fairly independent of the bremsstrahlung target used. Using the estimated increase in the RBE of neutrons compared to photons for absorbed dose levels applied in clinical fractionation schedules increases the RBE further to 1.13 at absorbed doses per fraction of the order of 2 Gy.

The purpose of this study was to quantify and model various types of cell death for a small-cell lung cancer (SCLC) cell line (U1690) after exposure to a 137Cs source and as well as to compare the linear-quadratic (LQ) and repairable-conditionally repairable model (RCR). This study is based on four different experiments that were taken place at Cancer Centrum Karolinska (CCK). A human small-cell lung cancer (SCLC) cell line after the exposure to a 137Cs source was used for the extraction of the clonogenic cell survival curve. Additionally, for the determination and quantification of various modes of cell death the method of fluorescence staining was implemented, where the cell deaths were categorized based on morphological characteristics. The percentage of cells in each phase of the cell cycle was investigated with flow cytometry analysis. The quantification of senescent cells was performed by staining the samples with senescence-associated β-galactosidase (SA-β-Gal) solution and then scoring as senescent cells those that had incorporated the substance. These data were introduced into a maximum likelihood fitting to calculate the best estimates of the parameters used by the examined model. In this model, the modes of cell death are divided into three categories: apoptotic, senescent and other types of cell death (necrotic/apoptotic, necrotic, micronuclei and giant). In the clonogenic cell survival assay, the fitting of the RCR model gives a χ(2)-value of 6.10 whereas for the LQ model became 9.61. In the fluorescence microscopy and senescence assay, the probability of the three different modes of cell death on day 2 seems to increases with a dose up to about 10 Gy where there is saturation. On day 7 a significant induction of apoptosis in a dose- and time-dependent manner was evident, whereas senescence was slightly increased in response to dose but not to time. As for the 'other types of cell death' mode on day 7 showed a higher probability than the one on day 2 and as well as a prominent dose-dependence. The RCR model fits better to the experimental data than the LQ model. On day 2 there is a slight increase of the apoptotic and senescent probability with dose. On the other hand, on day 7 the shape of the curve of apoptosis differs and a sigmoidal increase with dose is observed. At both time-points, the present model fits the data reasonably well. Due to the fact that the clonogenic survival does not coincide with the one extracted from the fluorescence microscopy, a more accurate way to quantify cell death needs to be used, e.g. computerized video time-lapse (CVTL).

The biological effects of particles are often expressed in relation to that of photons through the concept of relative biological effectiveness, RBE. In proton radiotherapy, a constant RBE of 1.1 is usually assumed. However, there is experimental evidence that RBE depends on various factors. The aim of this study is to develop a model to predict the RBE based on linear energy transfer (LET), dose, and the tissue specific parameter α/β of the linear-quadratic model for the reference radiation. Moreover, the model should capture the basic features of the RBE using a minimum of assumptions, each supported by experimental data.The α and β parameters for protons were studied with respect to their dependence on LET. An RBE model was proposed where the dependence of LET is affected by the (α/β)phot ratio of photons. Published cell survival data with a range of well-defined LETs and cell types were selected for model evaluation rendering a total of 10 cell lines and 24 RBE values.A statistically significant relation was found between α for protons and LET. Moreover, the strength of that relation varied significantly with (α/β)phot. In contrast, no significant relation between β and LET was found. On the whole, the resulting RBE model provided a significantly improved fit (p-value0.01) to the experimental data compared to the standard constant RBE. By accounting for the α/β ratio of photons, clearer trends between RBE and LET of protons were found, and our results suggest that late responding tissues are more sensitive to LET changes than early responding tissues and most tumors. An advantage with the proposed RBE model in optimization and evaluation of treatment plans is that it only requires dose, LET, and (α/β)phot as input parameters. Hence, no proton specific biological parameters are needed.

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.

The often large uncertainties that exist in beam patient alignment during radiation therapy may require modification of the incident beams to ensure an optimal delivered dose distribution to the target volume. This problem becomes increasingly severe when the required dose distribution of the incident beams becomes more heterogeneous. A simple analytical formula is derived for the case when the fraction number is high, and the desired relative dose variations are small. This formula adjusts the fluence distribution of the incident beam so that the resultant dose distribution will be as close as possible to the desired one considering the uncertainties in beam patient alignment. When sharp dose gradients are important, for instance at the border of the target volume, the problem is much more difficult. It is shown here that, if the tumor is surrounded by organs at risk, it is generally best to open up the field by about one standard deviation of the positional uncertainty--that is sigma/2 on each side of the target volume. In principle it is simultaneously desirable to increase the prescribed dose by a few per cent compared to the case where the positional uncertainty is negligible, in order to compensate for the rounded shoulders of the delivered dose distribution. When the tissues surrounding the tumor no longer are dose limiting even larger increases in field size may be advantageous. For more critical clinical situations the positional uncertainty may even limit the success of radiotherapy. In such cases one generally wants to create a steeper dose distribution than the underlying random Gaussian displacement process allows. The problem is then best handled by quantifying the treatment outcome under the influence of the stochastic process of patient misalignment. Either the coincidence with the desired dose distribution, or the expectation value of the probability of achieving complication-free tumor control is maximized under the influence of this stochastic process. It is shown that the most advantageous treatment is to apply beams that are either considerably widened or slightly widened and over flattened near the field edges for small and large fraction numbers respectively.

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 & 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 & 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.

In light-ion radiation therapy, both the dose and the local energy spectrum, which is often characterized with the linear energy transfer (LET), must be considered. In treatment optimization, it is advantageous to use a radiobiological model that analytically accounts for both dose and LET for the ion type of interest. With such a model the biological effect can also be estimated for dose and LET combinations for which there are no observations in the underlying experimental data. In this study, the repairable-conditionally repairable (RCR) damage model was extended by expressing its parameters as functions of LET to provide a radiobiological model that accounts for both the dose and the LET for a given ion type and cell line. This LET-parameterized RCR model was fitted to published cell survival data for HSG and V79 cells irradiated with carbon ions and for T1 cells irradiated with helium ions. To test the robustness of the model, fittings to only a subset of the data were performed. Good agreement with the cell survival data was obtained, including survival data for LET values not used for model fitting, opening up the possibility of using the model in treatment planning for light ions.

New radiobiological models are used to describe tumour and normal tissue reactions and to account for their dependence on the irradiated volume and inhomogeneities of the delivered dose distribution and cell sensitivity. The probability of accomplishing complication-free tumour control is maximized by an iterative algorithm. The algorithm is demonstrated by applying it to a one-dimensional (1D) tumour model but also to a more clinically relevant 2D case. The new algorithm is n-dimensional so it could simultaneously optimize the dose delivery in a 3D volume and in principle also select the ideal beam orientations, beam modalities (photons, electrons, neutrons, etc) and optimal spectral distributions of the corresponding modalities. To make calculation time reasonable, 2D-3D problems are most practical, and suitable beam orientations are preselected by the choice of irradiation kernel. The energy deposition kernel should therefore be selected in order to avoid irradiation through organs at risk. Clinically established dose response parameters for the tissues of interest are used to make the optimization as relevant as possible to the clinical problems at hand. The algorithm can be used even with a poorly selected kernel because it will always, as far as possible, avoid irradiating organs at risk. The generated dose distribution will be optimal with respect to the spatial distribution and assumed radiobiological properties of the tumour and normal tissues at risk for the kernel chosen. More specifically the probability of achieving tumour control without fatal complications in normal tissues is maximized. In the clinical examples a reduced tumour dose is seen at the border to sensitive organs at risk, but instead an increased dose just inside the tumour border is generated. The increased tumour dose has the effect that the dose fall-off is as steep as possible at the border to organs at risk.

The problem of choosing radiation quantities and units for energy deposition kernels and their associated kernel densities is treated with the aim of making them consistent with related classical radiation quantities and units such as restricted mass stopping powers and mass attenuation coefficients. It is shown that it is very useful to define the kernels h(r), in terms of the quotient of the mean specific energy imparted to the medium by the radiant energy incident on a volume element centred at the origin of the kernel. The basic building block used to generate these kernels is the point energy deposition kernel, h(p), describing the spatial distribution of the energy imparted by a photon interacting at a point in a medium. This will allow the kernels to be regarded as generalizations of the traditional mass stopping and attenuation coefficients, which in detail describe the spatial distribution of the mean energy deposition around an interaction site. As a consequence, the irradiation or kernel density, f(r) should be expressed in terms of the radiant energy incident per unit volume of the medium. It is shown that the kernel density is equal to minus the divergence of the incident unattenuated vectorial energy fluence, and it therefore acts as an irradiation density for the incident vectorial energy fluence. The microscopic kernels or the irradiation density may thus be viewed as a perfect 'sink' distribution to the required incident photon energy fluence which is totally absorbed at f(r), and instead replaced by the kernels which describe the detailed energy deposition in the medium in coordinates centred at the sinks. From these definitions the required incident energy fluence from an external radiation source used for treatment realization can be determined directly by projecting the irradiation density on the relevant positions of the radiation source. This procedure has the valuable property that maximal calculational accuracy is achieved in the tumour because the irradiation density has non-zero values only in the tumour, and the accuracy of the kernel is highest at its origin.

Object Radiation treatment of large arteriovenous malformations (AVMs) remains difficult and not very effective, even though seemingly promising methods such as staged volume treatments have been proposed by some radiation treatment centers. In symptomatic patients harboring large intracranial AVMs not amenable to embolization or resection, single-session high-dose stereotactic radiation therapy is a viable option, and the special characteristics of high-ionization-density light-ion beams offer several treatment advantages over photon and proton beams. These advantages include a more favorable depth-dose distribution in tissue, an almost negligible lateral scatter of the beam, a sharper penumbra, a steep dose falloff beyond the Bragg peak, and a higher probability of vascular response due to high ionization density and associated induction of endothelial cell proliferation and/or apoptosis. Carbon ions were recently shown to be an effective treatment for skull-base tumors. Bearing that in mind, the authors postulate that the unique physical and biological characteristics of light-ion beams should convey considerable clinical advantages in the treatment of large AVMs. In the present meta-analysis the authors present a comparison between light-ion beam therapy and more conventional modalities of radiation treatment with respect to these lesions. Methods Dose-volume histograms and data on peripheral radiation doses for treatment of large AVMs were collected from various radiation treatment centers. Dose-response parameters were then derived by applying a maximum likelihood fitting of a binomial model to these data. The present binomial model was needed because the effective number of crucial blood vessels in AVMs (the number of vessels that must be obliterated to effect a cure, such as large fistulous nidus vessels) is low, making the Poisson model less suitable. In this study the authors also focused on radiobiological differences between various radiation treatments. Results Light-ion Bragg-peak dose delivery has the precision required for treating very large AVMs as well as for delivering extremely sharp, focused beams to irregular lesions. Stereotactic light-ion radiosurgery resulted in better angiographically defined obliteration rates, less white-matter necrosis, lower complication rates, and more favorable clinical outcomes. In addition, in patients treated by He ion beams, a sharper dose-response gradient was observed, probably due to a more homogeneous radiosensitivity of the AVM nidus to light-ion beam radiation than that seen when low-ionization-density radiation modalities, such as photons and protons, are used. Conclusions Bragg-peak radiosurgery can be recommended for most large and irregular AVMs and for the treatment of lesions located in front of or adjacent to sensitive and functionally important brain structures. The unique physical and biological characteristics of light-ion beams are of considerable advantage for the treatment of AVMs: the densely ionizing beams of light ions create a better dose and biological effect distribution than conventional radiation modalities such as photons and protons. Using light ions, greater flexibility can be achieved while avoiding healthy critical structures such as diencephalic and brainstem nuclei and tracts. Treatment with the light ion He or Li is more suitable for AVMs ≤ 10 cm3, whereas treatment with the light ion Li, Be, or C may be more appropriate for larger AVMs. A binomial model based on the effective number of crucial vessels in the AVM may be used quite well to predict AVM obliteration probabilities for both small and large AVMs when therapies involving either photons or light ions are used.

Risk factors for development of a stricture of the upper esophagus after radiotherapy for head and neck cancer are poorly defined.This was a retrospective case-control study of patients diagnosed and treated for esophageal stricture after radiotherapy for head and neck cancer.The incidence of esophageal stricture after external beam radiation therapy (EBRT) was 3.3%. Seventy patients with stricture and 66 patients without stricture were identified. A multivariate analysis showed that there was increased risk of stricture in receiving enteral feeding during EBRT or in receiving a mean dose of45 Gy to the upper esophagus.Enteral feeding during EBRT is strongly associated with the development of stricture of the esophagus, as is a mean dose of45 Gy to the upper esophagus. Treatment of the stricture with Savary-Gilliard bougienage or through scope balloon dilatation is safe and successful but often has to be repeated.

In IMRT treatment plan optimization there are various methods that try to regularize the variation of dose nonuniformity 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. When a treatment plan optimization is performed using biological measures, the final goal should be some indication about the expected tumor control or normal tissue complications, which is the primary goal of treatment planning (the association of treatment configurations and dose prescription with the treatment outcome). In this study, this issue is analyzed distinguishing the dose-oriented treatment plan optimization from the response-oriented optimization. Three different dose distributions were obtained by using a dose-based optimization technique, an EUD-based optimization without applying any technique for regularizing the nonuniformity of the dose distribution, and an EUD-based optimization using a variational regularization technique, which controls dose nonuniformity. The clinical effectiveness of the three dose distributions was investigated by calculating the response probabilities of the tumors and organs-at-risk (OARs) involved in two head and neck and prostate cancer cases. The radiobiological models used are the linear-quadratic-Poisson and the Relative Seriality models. Furthermore, the complication-free tumor control probability and the biologically effective uniform dose (D) were used for treatment plan evaluation and comparison. The radiobiological comparison shows that the EUD-based optimization using L-curve regularization gives better results than the EUD-based optimization without regularization and dose-based optimization in both clinical cases. Concluding, it appears that the applied dose nonuniformity regularization technique is expected to improve the effectiveness of the optimized IMRT dose distributions. However, more patient cases are needed to validate the statistical significance of the results and conclusions presented in this paper.

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.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.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).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.

Introduction. In SRT/SRS, dedicated treatment planning systems are used for the calculation of the dose distribution. The majority of these systems utilize the standard TMR/OAR formalism for dose calculation as well as they usually neglect any perturbation due to head heterogeneities. The aim of this study is to examine the errors due to head heterogeneities for both absolute and relative dose distributions in stereotactic radiotherapy. Materials and methods. Dosimetric measurements in phantoms have been made for linac stereotactic irradiation. CT-based phantoms have been used for Monte Carlo simulations for both linac-based stereotactic system and Gamma Knife unit. Absolute and relative dose distributions have been compared between homogeneous and heterogeneous media. DVH and TCP results are presented for all cases. Results. The maximum absolute dose difference at the isocenter was 2.2% and 6.9% for the linac and Gamma Knife respectively. The impact of heterogeneity in the target DVH was minor for the linac technique whereas considerable difference was observed for the Gamma Knife treatment. This was reflected also to the radiobiological evaluation, where the maximum TCP difference for the linac system was 2.7% and for the Gamma Knife was 4%. Discussion and conclusions. The errors rising from the existence of head heterogeneities are not negligible especially for the Gamma Knife which uses lower energy beams. The errors of the absolute dose calculation could be easily eliminated by implementing a simple heterogeneity correction algorithm at the TPS. Nevertheless, the errors for not taking into account the lateral electron transport would require a more sophisticated approach and even direct Monte Carlo calculation.

The rapid implementation of advanced treatment planning and delivery technologies for radiation therapy has brought new challenges in evaluating the most effective treatment modality. Intensity-modulated radiotherapy (IMRT) using multi-leaf collimators (MLC) and helical tomotherapy (HT) are becoming popular modes of treatment delivery and their application and effectiveness continues to be investigated. Presently, there are several treatment planning systems (TPS) that can generate and optimize IMRT plans based on user-defined objective functions for the internal target volume (ITV) and organs at risk (OAR). However, the radiobiological parameters of the different tumours and normal tissues are typically not taken into account during dose prescription and optimization of a treatment plan or during plan evaluation. The suitability of a treatment plan is typically decided based on dosimetric criteria such as dose-volume histograms (DVH), maximum, minimum, mean and standard deviation of the dose distribution. For a more comprehensive treatment plan evaluation, the biologically effective uniform dose (D) is applied together with the complication-free tumour control probability (P(+)). Its utilization is demonstrated using three clinical cases that were planned with two different forms of IMRT. In this study, three different cancer types at different anatomical sites were investigated: head and neck, lung and prostate cancers. For each cancer type, a linac MLC-based step-and-shoot IMRT plan and a HT plan were developed. The MLC-based IMRT treatment plans were developed on the Philips treatment-planning platform, using the Pinnacle 7.6 software release. For the tomotherapy HiArt plans, the dedicated tomotherapy treatment planning station was used, running version 2.1.2. By using D as the common prescription point of the treatment plans and plotting the tissue response probabilities versus D for a range of prescription doses, a number of plan trials can be compared based on radiobiological measures. The applied plan evaluation method shows that in the head and neck cancer case the HT treatment gives better results than MLC-based IMRT in terms of expected clinical outcome P(+) of 62.2% and 46.0%, D to the ITV of 72.3 Gy and 70.7 Gy, respectively). In the lung cancer and prostate cancer cases, the MLC-based IMRT plans are better over the clinically useful dose prescription range. For the lung cancer case, the HT and MLC-based IMRT plans give a P(+) of 66.9% and 72.9%, D to the ITV of 64.0 Gy and 66.9 Gy, respectively. Similarly, for the prostate cancer case, the two radiation modalities give a P(+) of 68.7% and 72.2%, D to the ITV of 86.0 Gy and 85.9 Gy, respectively. If a higher risk of complications (higher than 5%) could be allowed, the complication-free tumour control could increase by over 40%, 2% and 30% compared to the initial dose prescription for the three cancer cases, respectively. Both MLC-based IMRT and HT can encompass the often-large ITV required while they minimize the volume of the organs at risk receiving high doses. Radiobiological evaluation of treatment plans may provide an improved correlation of the delivered treatment with the clinical outcome by taking into account the dose-response characteristics of the irradiated targets and normal tissues. There may exist clinical cases, which may look dosimetrically similar but in radiobiological terms may be quite different. In such situations, traditional dose-based evaluation tools can be complemented by the use of P(+)--D diagrams to effectively evaluate and compare treatment plans.

In a recently published paper (Nioutsikou et al 2005 Phys. Med. Biol. 50 L17) the authors showed that the use of the dose-mass histogram (DMH) concept is a more accurate descriptor of the dose delivered to lung than the traditionally used dose-volume histogram (DVH) concept. Furthermore, they state that if a functional imaging modality could also be registered to the anatomical imaging modality providing a functional weighting across the organ (functional mass) then the more general and realistic concept of the dose-functioning mass histogram (D[F]MH) could be an even more appropriate descriptor. The comments of the present letter to the editor are in line with the basic arguments of that work since their general conclusions appear to be supported by the comparison of the DMH and DVH concepts using radiobiological measures. In this study, it is examined whether the dose-mass histogram (DMH) concept deviated significantly from the widely used dose-volume histogram (DVH) concept regarding the expected lung complications and if there are clinical indications supporting these results. The problem was investigated theoretically by applying two hypothetical dose distributions (Gaussian and semi-Gaussian shaped) on two lungs of uniform and varying densities. The influence of the deviation between DVHs and DMHs on the treatment outcome was estimated by using the relative seriality and LKB models using the Gagliardi et al (2000 Int. J. Radiat. Oncol. Biol. Phys. 46 373) and Seppenwoolde et al (2003 Int. J. Radiat. Oncol. Biol. Phys. 55 724) parameter sets for radiation pneumonitis, respectively. Furthermore, the biological equivalent of their difference was estimated by the biologically effective uniform dose (D) and equivalent uniform dose (EUD) concepts, respectively. It is shown that the relation between the DVHs and DMHs varies depending on the underlying cell density distribution and the applied dose distribution. However, the range of their deviation in terms of the expected clinical outcome was proven to be very large. Concluding, the effectiveness of the dose distribution delivered to the patients seems to be more closely related to the radiation effects when using the DMH concept.

Aims Because of the highly conformal distributions that can be obtained with intensity-modulated radiotherapy (IMRT), any discrepancy between the intended and delivered distributions would probably affect the clinical outcome. Consequently, there is a need for a measure that would quantify those differences in terms of a change in the expected clinical outcome. Materials and methods To evaluate such a measure, cancer of the cervix was used, where the bladder and rectum are proximal and partially overlapping with the internal target volume. A solid phantom simulating the pelvic anatomy was fabricated and a treatment plan was developed to deliver the prescribed dose to the phantom. The phantom was then irradiated with films positioned in several transverse planes. The racetrack microtron at 50MV was used in the treatment planning and delivery processes. The dose distribution delivered was analysed based on the film measurements and compared against the treatment plan. The differences in the measurements were evaluated using both physical and biological criteria. Whereas the physical comparison of dose distributions can assess the geometric accuracy of delivery, it does not reflect the clinical effect of any measured dose discrepancies. Results It is shown how small inaccuracies in delivered dose can affect the treatment outcome in terms of complication-free tumour cure. Conclusions With highly conformal IMRT, the accuracy of the patient set-up and treatment delivery are critical for the success of the treatment. A method is proposed to evaluate the precision of the delivered plan based on changes in complication and control rates as they relate to uncertainties in dose delivery.

The purpose of this work is to evaluate the predictive strength of the relative seriality, parallel and LKB normal tissue complication probability (NTCP) models regarding the incidence of radiation pneumonitis, in a large group of patients following breast cancer radiotherapy, and furthermore, to illustrate statistical methods for examining whether certain published radiobiological parameters are compatible with a clinical treatment methodology and patient group characteristics. The study is based on 150 consecutive patients who received radiation therapy for breast cancer. For each patient, the 3D dose distribution delivered to lung and the clinical treatment outcome were available. Clinical symptoms and radiological findings, along with a patient questionnaire, were used to assess the manifestation of radiation-induced complications. Using this material, different methods of estimating the likelihood of radiation effects were evaluated. This was attempted by analysing patient data based on their full dose distributions and associating the calculated complication rates with the clinical follow-up records. Additionally, the need for an update of the criteria that are being used in the current clinical practice was also examined. The patient material was selected without any conscious bias regarding the radiotherapy treatment technique used. The treatment data of each patient were applied to the relative seriality, LKB and parallel NTCP models, using published parameter sets. Of the 150 patients, 15 experienced radiation-induced pneumonitis (grade 2) according to the radiation pneumonitis scoring criteria used. Of the NTCP models examined, the relative seriality model was able to predict the incidence of radiation pneumonitis with acceptable accuracy, although radiation pneumonitis was developed by only a few patients. In the case of modern breast radiotherapy, radiobiological modelling appears to be very sensitive to model and parameter selection giving clinically acceptable results in certain cases selectively (relative seriality model with Seppenwoolde et al and Gagliardi et al parameter sets). The use of published parameters should be considered as safe only after their examination using local clinical data. The variation of inter-patient radiosensitivity seems to play a significant role in the prediction of such low incidence rate complications. Scoring grades were combined to give stronger evidence of radiation pneumonitis since their differences could not be strictly associated with dose. This obviously reveals a weakness of the scoring related to this endpoint, and implies that the probability of radiation pneumonitis induction may be too low to be statistically analysed with high accuracy, at least with the latest advances of dose delivery in breast radiotherapy.

The aim of the study was to examine the influence of overall treatment time (OTT) on the value of calculated biological effective doses (BEDs) for different biological variables. These variables were: tumour proliferation rate, different cell radiosensitivity (alpha=0.2, 0.3, and 0.4 /Gy), and different start time for repopulation (Tk=21, 28, and 35 days). Also the influence of age (/= 50 years), Hb level (/= 116 g/l), tumor proliferation rate (bromodeoxyuridine labelling index; BrdUrdLI), and DNA ploidy on survival after shorter (/= 60 days) or longer (60 days) OTT was investigated. The study included 229 patients with cervix carcinoma treated entirely by standard radiotherapy (RT) (external beam RT plus low-medium dose-rate (LDR/MDR) brachytherapy (BT) at the Center of Oncology in Krakow. The linear quadratic equation was used to calculate BED, which is proportional to log cell kill. BEDs 10 (for tumours) were calculated with consideration of OTT for each patient and tumour proliferation rate (standardized potential doubling time; standardized Tpot) based on BrdUrdLI assessed on biopsy material before RT. Median OTT was 90 days (range 30-210). The mean calculated total BED for point A for tumour and 'early reactions' was equal to 103.0 Gy10. The longest median survival time--52 months--was seen for patients treated with OTT/= 60 days. If OTT exceeded 90 days to more than 120 days, loss in BED10 for relatively radiosensitive tumours (alpha=0.3-0.4/Gy and Tk=28 days) was equal to 0.37-0.26 Gy/day. However, for radioresistant tumours (alpha=0.2/Gy) it was 0.6 Gy/day. For fast proliferating tumours (BrdUrdLI8.8%) BED loss was 1.4 Gy/day and for slowly proliferating tumours (BrdUrdLI/= 8.8%) it was 0.2 Gy/day. Assuming shorter (21 days) or longer (35 days) periods for Tk and relatively radiosensitive tumours similar BED loss of 0.38 Gy/day was observed. Kaplan-Meier analysis revealed that OTT/= 60 days was a significant prognostic factor for overall survival (OS) (p=0.019), disease-free survival (DFS) (p=0.0173), and local control (LC) (p=0.011). BED10 had significant influence on survival (p=0.047). Cox multivariate analysis revealed that for OTT shorter than 60 days the only favourable significant parameters were: age50 years (p=0.003) and high Hb level (116 g/l) (p=0.041). For longer treatments (OTT60 days) the unfavourable parameters were: age/= 50 years (p=0.037), BrdUrdLI/= 8.8% (p=0.003), tumour aneuploidy (p=0.043), and BED10/= 103 Gy (p=0.017). The examined tumour biological parameters should be taken into account for RT and provide a basis for adjuvant RT.

Background: The estimation of the parameters that describe the dose-response relations of anal sphincter regarding the clinical endpoints of fecal leakage and blood or phlegm in stools is important in the optimization of prostate cancer radiotherapy. Also, the validity of the relative seriality model for this clinical case needs to be examined by associating the clinical follow-up results with the predicted complication rates. Patients and Methods: In this study, 65 patients who received radiation therapy for clinically localized prostate adenocarcinoma are analyzed. The clinical treatment outcome and the three-dimensional dose distribution delivered to anal sphincter were available for each patient. A questionnaire was used for assessing the clinical bowel and urinary symptoms. A maximum likelihood fitting was performed to calculate the best estimates of the parameters used by the relative seriality model. The clinical utilization of the calculated parameters in predicting anal sphincter complication probabilities was illustrated by applying the best estimate of the parameters to a subset of the patient population. Results: The estimated values of the parameters for the two clinical endpoints are D 50 = 70.2 Gy, γ = 1.22, s = 0.35 for fecal leakage and D 50 = 74.0 Gy, γ = 0.75, s ≈ 0 for blood or phlegm in stools. The standard deviations of the parameters were also calculated together with the confidence intervals of the dose-response curves. The analysis proved that the treatment outcome pattern of the patient material can suitably be reproduced by the relative seriality model (probability of finding a worse fit = 60.2%, the area under the receiver operating characteristic curve = 0.72 and 0.69 and χ2-test = 0.97 and 0.86, respectively). Conclusion: Fecal leakage is characterized by a medium relative seriality whereas blood or phlegm in stools was found to have strong volume dependence (low relative seriality). Diminishing the biologically effective uniform dose to anal sphincter < 40–45 Gy may significantly reduce the risk of fecal leakage or blood or phlegm in stools for patients irradiated for prostate cancer.

The purpose of this work is to determine the parameters and evaluate the predictive strength of the relative seriality model. This is accomplished by associating the calculated complication rates with the clinical follow-up records. The study is based on 82 patients who received radiation treatment for head and neck cancer. For each patient the 3D dose distribution delivered to the esophagus and the clinical treatment outcome were available. Clinical symptoms and radiological findings were used to assess the manifestation of radiation-induced esophageal strictures. These data were introduced into a maximum likelihood fitting to calculate the best estimates of the parameters used by the relative seriality model (D50 = 68.4 Gy, gamma = 6.55, s = 0.22). The uncertainties of these parameters were also calculated and their individual influence on the dose-response curve was demonstrated. The best estimate of the parameters was applied to 58 patients of the study material and their esophageal stricture induction probabilities were calculated to illustrate the clinical utilization of the calculated parameters. The calculation of the biological effective dose (BED) appeared to be significantly sensitive to the applied fractionation correction for complex treatment plans. The relative seriality model was proved suitable in reproducing the treatment outcome pattern of the patient material studied (probability of finding a worse fit = 61.0%, the area under the ROC curve = 0.84 and chi2 test = 0.95). The analysis was carried out for the upper 5 cm of the esophagus (proximal esophagus) where all the strictures are formed. Radiation-induced strictures were found to have a strong volume dependence (low relative seriality). The uncertainties of the parameters appear to have a significant supporting role on the estimated dose-response curve.

To investigate dysfunction caused by unwanted radiation to the anal-sphincter region and the rectum.A questionnaire assessing bowel symptoms, sexual function, and urinary symptoms was sent to 72 patients with clinically localized prostatic adenocarcinoma treated by external beam radiation therapy at the Radiumhemmet, Karolinska Hospital, in Stockholm, Sweden, 2-4 years after treatment. The mean percentage dose-volume histograms for patients with and without the specific symptom were calculated.Of the 65 patients providing information, 9 reported fecal leakage, 10 blood and mucus in stools, 10 defecation urgency, and 7 diarrhea or loose stools. None of the 19 and 13 patients who received, respectively, a dose ofor =35 Gy toor =60% oror =40 Gy toor =40% of the anal-sphincter region volume reported fecal leakage (p0.05). In dose-volume histograms, a statistically significant correlation was found between radiation to the anal-sphincter region and the risk of fecal leakage in the interval 45-55 Gy. There was also a statistically significant correlation between radiation to the rectum and the risk of defecation urgency and diarrhea or loose stools in the interval 25-42 Gy. No relationship was found between anatomic rectal wall volume and the investigated late effects.Although the limited data in this study prevent the definition of a conclusive threshold regarding volume and dose to the anal-sphincter region and untoward morbidity, it seems that careful monitoring of unnecessary irradiation to this area should be done because it can potentially help reduce the risk of adverse effects, such as fecal leakage. Future studies should pay more attention to the anal-sphincter region and help to more rigorously define its radiotherapeutic tolerance.

BACKGROUND It is well recognized that many patients with head and neck carcinoma have problems with food intake and malnutrition. The objective of the current study was to determine the clinical pattern of patients with nonneoplastic stricture of the upper esophagus after radiotherapy for head and neck carcinoma. METHODS A retrospective chart study of 22 patients with stricture of the proximal esophagus diagnosed between 1993 and 1999 at Karolinska Hospital was performed. The dose volume histograms of the first 2 cm and 5 cm, respectively, of the proximal esophagus were calculated. RESULTS Five of the patients (23%) had total obliteration. The first 2 cm of the esophagus received at least 60 grays (Gy) in > 80% of the volume. Radiation injury was not reported to occur at doses < 60 Gy. There was a correlation found between dysphagia during radiotherapy and the development of proximal esophageal stricture. Stricture was diagnosed 1–60 months (median, 6 months) after radiotherapy. In 18 patients, the stricture was treated with single or repeated endoscopic dilation. These treatments allowed a nearly normal diet in 78% of the patients. CONCLUSIONS Stricture of the upper esophagus is one deglutition disorder that is reported to occur after radiotherapy for head and neck carcinoma. In the current study, the authors emphasize the importance of knowing the tolerance of the normal esophagus to irradiation as well as early diagnosis of stricture of the proximal esophagus because this condition may lead to physical and emotional distress. Cancer 2003;97:1693–700. © 2003 American Cancer Society. DOI 10.1002/cncr.11236

The quality of the radiation therapy delivered in the treatment of breast cancer is susceptible to setup errors and organ motion uncertainties. For 60 breast cancer patients (24 resected with negative node involvement, 13 resected with positive node involvement and 23 ablated) who were treated with three different irradiation techniques. these uncertainties are simulated. The delivered dose distributions in the lung were recalculated taking positioning uncertainty and breathing effects into account. In this way the real dose distributions delivered to the patients are more closely determined. The positioning uncertainties in the anteroposterior (AP) and the craniocaudal (CC) directions are approximated by Gaussian distributions based on the fact that setup errors are random. Breathing is assumed to have a linear behavior because of the chest wall movement during expiration and inspiration. The combined frequency distribution of the positioning and breathing distributions is obtained by convolution. By integrating the convolved distribution over a number of intervals, the positions and the weights of the fields that simulate the original 'effective fields' are calculated. Opposed tangential fields are simulated by a set of 5 pairs of fields in the AP direction and 3 such sets in the CC direction. Opposed AP + PA fields are simulated by a set of 3 pairs of fields in the AP direction and 3 such sets in the CC direction. Single frontal fields are simulated by a set of 5 fields. In radiotherapy for breast cancer, the lung is often partly within the irradiated volume even though it is a sensitive organ at risk. The influence of the deviation in the dose delivered by the original and the adjusted treatment plans on the clinical outcome is estimated by using the relative seriality model and the biologically effective uniform dose concept. Radiation pneumonitis is used as the clinical endpoint for lung complications. The adjusted treatment plans show larger lung complication probabilities than the original plans. This means that the true expected complications are often underestimated in clinical practice. The lung density variation during breathing is calculated from the maximal change in average density during tidal breathing. The change in density in the lung due to breathing is shown to have almost no influence on the dose distribution in the lung. The proposed treatment-plan adjustments taking positioning uncertainty and breathing effects into account indicate significant deviations in the dose delivery and the predicted lung complications.

The prognostic significance of intratumour microvessel density (IMD) and haemoglobin (Hb) level was studied in 152 (FIGO stage IB-IIIB) cervical cancer patients before radiotherapy. Patients' age and tumour stage, grade and degree of keratinization were also studied. IMD measurement expressed as the mean vessel count per 1 mm2 was performed on formalin-fixed, paraffin-embedded tumour biopsies stained with anti-factor VIII antibody (DAKO Ltd.) using immunohistochemistry and the vascular hot-spot technique. The median age of patients was 55 years (29-80). Median values for IMD and Hb level were: 142.5 (range 56.3-476.6) vessels/mm2 and 129 (range 81-160) gil, respectively. The median time of follow-up was 26 months, with a range of 2-145 months. Tumour stage (p = 0.7009), grade (p = 0.6660) and degree of keratinization (0.2669) were not significant in the Kaplan-Meier univariate analysis. However, patients' age50 years (p = 0.0079). high vascularity (IMD190.0 vessels/mm2, minimal cut-off point, p = 0.0503) and Hb concentration116 g/l (p = 0.0213) were favourable prognostic factors for cervical cancer treated with radiotherapy. In the Cox multivariate analysis only higher vascularity (IMD190/mm2 and Hb concentration116 g/l were favourable prognostic factors in terms of patients' survival. However, when a Cox analysis was done separately for keratinizing and non-keratinizing tumours, it was found that higher vascularity was significant only for keratinizing, and higher Hb level only for non-keratinizing cancers.