Data on dose–volume effects in the rat spinal cord do not support existing NTCP models

Can dose convolution modelling explain bath and shower effects in rat spinal cord?

Objective. ‘Bath and shower’ effects were first seen in proton irradiations of rat spinal cord, where a low dose ‘bath’ reduced the smaller field ‘shower’ dose needed for limb paralysis giving the appearance of sensitisation of the cord or disproportionate response. This was difficult to reconcile with existing tissue complication models. The purpose of this investigation is to explore a different approach using a dose convolution algorithm to model the 50% isoeffect endpoint. Approach. Bath and shower dose distributions were convolved with Gaussian functions with widths specified by the σ parameter. The hypothesis was that the maximum value from the convolved distributions was constant for isoeffect across the modelled scenarios. A simpler field length dependent relative biological effectiveness (FLRBE) approach was also used for a subset of the data which gave results independent of σ. Main results. The maximum values from the convolved distributions were constant within ±17% across the bath and shower experiments for σ = 3.5 mm, whereas the maximum dose varied by a factor of four. The FLRBE results were also within ±14% confirming the validity of the dose convolution approach. Significance. A simple approach using dose convolution modelling of the 50% isotoxicity gave compelling consistency with the full range of bath and shower results, while the FLRBE approach confirmed the results for the symmetric field data. Convolution modelling and the effect of time interval were consistent with a signalling factor diffusion mechanism such as the ‘bystander effect’. The results suggest biological effectiveness is reduced for very small field sizes, requiring a higher isoeffect dose. By implication, the bath dose does not sensitise the cord to the shower dose; when biological effectiveness is accounted for, a small increase in the bath dose requires a significantly larger reduction in shower dose for isoeffect.

Effects From Nonuniform Dose Distribution in the Spinal Nerves of Pigs: Analysis of Normal Tissue Complication Probability Models

PURPOSE

Our purpose was to evaluate normal tissue complication probability (NTCP) models for their ability to describe the increase in tolerance as the length of irradiated spinal nerve is reduced in a pig.

METHODS AND MATERIALS

Common phenomenological and semimechanistic NTCP models were fit using the maximum likelihood estimate method to dose-response data from spinal nerve irradiation studies in pigs. Statistical analysis was used to compare how well each model fit the data. Model parameters were then applied to a previously published dose distribution used for spinal cord irradiation in rats under the assumption of a similar dose-response.

RESULTS

The Lyman-Kutcher-Burman model, relative seriality, and critical volume model fit the spinal nerve data equally well, but the mean dose logistic and relative seriality models gave the best fit after penalizing for the number of model parameters. The minimum dose logistic regression model was the only model showing a lack of fit. When extrapolated to a 0.5-cm simulated square-wave-like dose distribution, the serial behaving models showed negligible increase in dose-response curve. The Lyman-Kutcher-Burman model and relative seriality models showed significant shifting of NTCP curves due to parallel behaving parameters. The critical volume model gave the closest match to the rat data.

CONCLUSIONS

Several phenomenological and semimechanistic models were observed to adequately describe the increase in the radiation tolerance of the spinal nerves when changing the irradiated length from 1.5 to 0.5 cm. Contrary to common perception, model parameters suggest parallel behaving tissue architecture. Under the assumption that the spinal nerve response to radiation is similar to that of the spinal cord, only the critical volume model was robust when extrapolating to outcome data from a 0.5-cm square-wave-like dose distribution, as was delivered in rodent spinal cord irradiation research.

Regional Responses in Radiation-Induced Normal Tissue Damage

Normal tissue side effects remain a major concern in radiotherapy. The improved precision of radiation dose delivery of recent technological developments in radiotherapy has the potential to reduce the radiation dose to organ regions that contribute the most to the development of side effects. This review discusses the contribution of regional variation in radiation responses in several organs. In the brain, various regions were found to contribute to radiation-induced neurocognitive dysfunction. In the parotid gland, the region containing the major ducts was found to be critical in hyposalivation. The heart and lung were each found to exhibit regional responses while also mutually affecting each other's response to radiation. Sub-structures critical for the development of side effects were identified in the pancreas and bladder. The presence of these regional responses is based on a non-uniform distribution of target cells or sub-structures critical for organ function. These characteristics are common to most organs in the body and we therefore hypothesize that regional responses in radiation-induced normal tissue damage may be a shared occurrence. Further investigations will offer new opportunities to reduce normal tissue side effects of radiotherapy using modern and high-precision technologies.

Dose-Limiting Organs at Risk in Carbon Ion Re-Irradiation of Head and Neck Malignancies: An Individual Risk-Benefit Tradeoff

Background: Carbon ion re-irradiation (CIR) was evaluated to investigate treatment planning and the consequences of individual risk-benefit evaluations concerning dose-limiting organs at risk (OAR). Methods: A total of 115 consecutive patients with recurrent head and neck cancer (HNC) were analyzed after initial radiotherapy and CIR at the same anatomical site. Toxicities were evaluated in line with the Common Terminology Criteria for Adverse Events 4.03. Results: The median maximum cumulative equivalent doses applied in fractions of 2 Gy (EQD2) to the brainstem, optic chiasm, ipsilateral optic nerve, and spinal cord were 56.8 Gy (range 0.94-103.9), 51.4 Gy (range 0-120.3 Gy), 63.6 Gy (range 0-146.1 Gy), and 28.8 Gy (range 0.2-87.7 Gy). The median follow up after CIR was 24.0 months (range 2.5-72.0 months). The cumulative rates of acute and late severe (≥grade III) side effects after CIR were 1.8% and 14.3%. Conclusion: In recurrent HNC, an individual risk-benefit tradeoff is frequently inevitable due to unfavorable location of tumors in close proximity to vital OAR. There are uncertainties about the dose tolerance of OAR after CIR, which warrant increased awareness about the potential treatment toxicity and further studies on heavy ion re-irradiation.

Spatially fractionated proton minibeams

Extraordinary normal tissue response to highly spatially fractionated X-ray beams has been explored for over 25 years. More recently, alternative radiation sources have been developed and utilized with the aim to evoke comparable effects. These include protons, which lend themselves well for this endeavour due to their physical depth dose characteristics as well as corresponding variable biological effectiveness. This paper addresses the motivation for using protons to generate spatially fractionated beams and reviews the technological implementations and experimental results to date. This includes simulation and feasibility studies, collimation and beam characteristics, dosimetry and biological considerations as well as the results of in vivo and in vitro studies. Experimental results are emerging indicating an extraordinary normal tissue sparing effect analogous to what has been observed for synchrotron generated X-ray microbeams. The potential for translational research and feasibility of spatially modulated proton beams in clinical settings is discussed.

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

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

Biological considerations when comparing proton therapy with photon therapy

Owing to the limited availability of data on the outcome of proton therapy, treatments are generally optimized based on broadly available data on photon-based treatments. However, the microscopic pattern of energy deposition of protons differs from that of photons, leading to a different biological effect. Consequently, proton therapy needs a correction factor (relative biological effectiveness) to relate proton doses to photon doses, and currently, a generic value is used. Moreover, the macroscopic distribution of dose in proton therapy differs compared with photon treatments. Although this may offer new opportunities to reduce dose to normal tissues, it raises the question whether data obtained from photon-based treatments offer sufficient information on dose-volume effects to optimally use unique features of protons. In addition, there are potential differences in late effects due to low doses of secondary radiation outside the volume irradiated by the primary beam. This article discusses the controversies associated with these 3 issues when comparing proton and photon therapy.

Novel approaches to improve the therapeutic index of head and neck radiotherapy: an analysis of data from the PARSPORT randomised phase III trial

PURPOSE

Subjective xerostomia is a common side-effect following radiotherapy for the treatment of head-and-neck cancer. Standard mean dose models previously used to model xerostomia only that partially predict the occurrence of xerostomia. Studies in animal models have suggested that there are regional variations in the radiosensitivity of the parotid glands. In this work we tested the hypothesis that this is also true for the human parotid gland.

METHODS

We present novel dose-response models explicitly taking the spatial distribution of the radiation dose into account. We considered dose to the submandibular gland and other clinical factors and used a variable-selection algorithm to select the best dose-response model. This methodology was applied to 63 head and neck cancer patients and validated using two independent patient cohorts of 19 and 29 patients, respectively.

RESULTS

The predictive accuracy of dose-response models improved significantly when including regional variations of radiosensitivity of the parotid glands compared to standard mean-dose models (p = 0.001, t-test). Beneficial dose-pattern analysis demonstrated the importance of minimising dose to the lateral and cranial component of the human parotid gland in order to avoid xerostomia. Furthermore we found an evidence that surgical removal of the sub-mandibular gland significantly increases the risk of radiation-induced xerostomia.

CONCLUSION

Dose-response models which take the shape of the dose-distribution into account predicted xerostomia significantly better than standard mean-dose models. Our novel model could be used to rank potential treatment plans more reliably according to their therapeutic index and may be useful to generate better treatment plans.

External validation of three dimensional conformal radiotherapy based NTCP models for patient-rated xerostomia and sticky saliva among patients treated with intensity modulated radiotherapy

PURPOSE

The purpose of this study was to investigate the ability of predictive models for patient-rated xerostomia (XER(6M)) and sticky saliva (STIC(6M)) at 6 months after completion of primary (chemo)radiation developed in head and neck cancer patients treated with 3D-conformal radiotherapy (3D-CRT) to predict outcome in patients treated with intensity modulated radiotherapy (IMRT).

METHODS AND MATERIALS

Recently, we published the results of a prospective study on predictive models for patient-rated xerostomia and sticky saliva in head and neck cancer patients treated with 3D-CRT (3D-CRT based NTCP models). The 3D-CRT based model for XER(6M) consisted of three factors, including the mean parotid dose, age, and baseline xerostomia (none versus a bit). The 3D-CRT based model for STIC(6M) consisted of the mean submandibular dose, age, the mean sublingual dose, and baseline sticky saliva (none versus a bit). In the current study, a population consisting of 162 patients treated with IMRT was used to test the external validity of these 3D-CRT based models. External validity was described by the explained variation (R(2) Nagelkerke) and the Brier score. The discriminative abilities of the models were calculated using the area under the receiver operating curve (AUC) and calibration (i.e. the agreement between predicted and observed outcome) was assessed with the Hosmer-Lemeshow "goodness-of-fit" test.

RESULTS

Overall model performance of the 3D-CRT based predictive models for XER(6M) and STIC(6M) was significantly worse in terms of the Brier score and R(2) Nagelkerke among patients treated with IMRT. Moreover the AUC for both 3D-CRT based models in the IMRT treated patients were markedly lower. The Hosmer-Lemeshow test showed a significant disagreement for both models between predicted risk and observed outcome.

CONCLUSION

3D-CRT based models for patient-rated xerostomia and sticky saliva among head and neck cancer patients treated with primary radiotherapy or chemoradiation turned out to be less valid for patients treated with IMRT. The main message from these findings is that models developed in a population treated with a specific technique cannot be generalised and extrapolated to a population treated with another technique without external validation.

Comparison of dose response models for predicting normal tissue complications from cancer radiotherapy: application in rat spinal cord

Seven different radiobiological dose-response models have been compared with regard to their ability to describe experimental data. The first four models, namely the critical volume, the relative seriality, the inverse tumor and the critical element models are mainly based on cell survival biology. The other three models: the Lyman (Gaussian distribution), the parallel architecture and the Weibull distribution models are semi-empirical and rather based on statistical distributions. The maximum likelihood estimation was used to fit the models to experimental data and the χ2-distribution, AIC criterion and F-test were applied to compare the goodness-of-fit of the models. The comparison was performed using experimental data for rat spinal cord injury. Both the shape of the dose-response curve and the ability of handling the volume dependence were separately compared for each model. All the models were found to be acceptable in describing the present experimental dataset (p > 0.05). For the white matter necrosis dataset, the Weibull and Lyman models were clearly superior to the other models, whereas for the vascular damage case, the Relative Seriality model seems to have the best performance although the Critical volume, Inverse tumor, Critical element and Parallel architecture models gave similar results. Although the differences between many of the investigated models are rather small, they still may be of importance in indicating the advantages and limitations of each particular model. It appears that most of the models have favorable properties for describing dose-response data, which indicates that they may be suitable to be used in biologically optimized intensity modulated radiation therapy planning, provided a proper estimation of their radiobiological parameters had been performed for every tissue and clinical endpoint.

Spinal cord tolerance in the age of spinal radiosurgery: lessons from preclinical studies

Clinical implementation of spinal radiosurgery has increased rapidly in recent years, but little is known regarding human spinal cord tolerance to single-fraction irradiation. In contrast, preclinical studies in single-fraction spinal cord tolerance have been ongoing since the 1970s. The influences of field length, dose rate, inhomogeneous dose distributions, and reirradiation have all been investigated. This review summarizes literature regarding single-fraction spinal cord tolerance in preclinical models with an emphasis on practical clinical significance. The outcomes of studies that incorporate uniform irradiation are surprisingly consistent among multiple small- and large-animal models. Extensive investigation of inhomogeneous dose distributions in the rat has demonstrated a significant dose-volume effect while preliminary results from one pig study are contradictory. Preclinical spinal cord dose-volume studies indicate that dose distribution is more critical than the volume irradiated suggesting that neither dose-volume histogram analysis nor absolute volume constraints are effective in predicting complications. Reirradiation data are sparse, but results from guinea pig, rat, and pig studies are consistent with the hypothesis that the spinal cord possesses a large capacity for repair. The mechanisms behind the phenomena observed in spinal cord studies are not readily explained and the ability of dose response models to predict outcomes is variable underscoring the need for further investigation. Animal studies provide insight into the phenomena and mechanisms of radiosensitivity but the true significance of animal studies can only be discovered through clinical trials.

The risks to healthy tissues from the use of existing and emerging techniques for radiation therapy

As radical radiotherapy treatments become more effective, more and more cancer patients are becoming cured of their disease and surviving for decades. Damage to exposed healthy tissues that becomes manifest in the medium-to-long-term is becoming a more significant factor in the choice of individual treatment plans and treatment modality. However, currently there are no reliable objective methods for predicting in an individual patient the occurrence of normal tissue complications, or second cancers caused by radiation. This is especially needed as new competing techniques and modalities become available, such as IMRT, protons, carbon ions, etc., all advancing the ability to focus the radiation dose on the target while sparing normal tissue. ALLEGRO is a Euratom-funded project that is currently investigating the current state of knowledge, and attempting to define the priority research areas. Preliminary considerations of the problems to be solved and research priorities are presented.

Quantitative Analyses of Normal Tissue Effects in the Clinic (QUANTEC): an introduction to the scientific issues

Advances in dose-volume/outcome (or normal tissue complication probability, NTCP) modeling since the seminal Emami paper from 1991 are reviewed. There has been some progress with an increasing number of studies on large patient samples with three-dimensional dosimetry. Nevertheless, NTCP models are not ideal. Issues related to the grading of side effects, selection of appropriate statistical methods, testing of internal and external model validity, and quantification of predictive power and statistical uncertainty, all limit the usefulness of much of the published literature. Synthesis (meta-analysis) of data from multiple studies is often impossible because of suboptimal primary analysis, insufficient reporting and variations in the models and predictors analyzed. Clinical limitations to the current knowledge base include the need for more data on the effect of patient-related cofactors, interactions between dose distribution and cytotoxic or molecular targeted agents, and the effect of dose fractions and overall treatment time in relation to nonuniform dose distributions. Research priorities for the next 5-10 years are proposed.

Assessing correlations between the spatial distribution of the dose to the rectal wall and late rectal toxicity after prostate radiotherapy: an analysis of data from the MRC RT01 trial (ISRCTN 47772397)

Many studies have been performed to assess correlations between measures derived from dose-volume histograms and late rectal toxicities for radiotherapy of prostate cancer. The purpose of this study was to quantify correlations between measures describing the shape and location of the dose distribution and different outcomes. The dose to the rectal wall was projected on a two-dimensional map. In order to characterize the dose distribution, its centre of mass, longitudinal and lateral extent, and eccentricity were calculated at different dose levels. Furthermore, the dose-surface histogram (DSH) was determined. Correlations between these measures and seven clinically relevant rectal-toxicity endpoints were quantified by maximally selected standardized Wilcoxon rank statistics. The analysis was performed using data from the RT01 prostate radiotherapy trial. For some endpoints, the shape of the dose distribution is more strongly correlated with the outcome than simple DSHs. Rectal bleeding was most strongly correlated with the lateral extent of the dose distribution. For loose stools, the strongest correlations were found for longitudinal extent; proctitis was most strongly correlated with DSH. For the other endpoints no statistically significant correlations could be found. The strengths of the correlations between the shape of the dose distribution and outcome differed considerably between the different endpoints. Due to these significant correlations, it is desirable to use shape-based tools in order to assess the quality of a dose distribution.

The significance of the choice of radiobiological (NTCP) models in treatment plan objective functions

A Clinician's discrimination between radiation therapy treatment plans is traditionally a subjective process, based on experience and existing protocols. A more objective and quantitative approach to distinguish between treatment plans is to use radiobiological or dosimetric objective functions, based on radiobiological or dosimetric models. The efficacy of models is not well understood, nor is the correlation of the rank of plans resulting from the use of models compared to the traditional subjective approach. One such radiobiological model is the Normal Tissue Complication Probability (NTCP). Dosimetric models or indicators are more accepted in clinical practice. In this study, three radiobiological models, Lyman NTCP, critical volume NTCP and relative seriality NTCP, and three dosimetric models, Mean Lung Dose (MLD) and the Lung volumes irradiated at 10Gy (V10) and 20Gy (V20), were used to rank a series of treatment plans using, harm to normal (Lung) tissue as the objective criterion. None of the models considered in this study showed consistent correlation with the Radiation Oncologists plan ranking. If radiobiological or dosimetric models are to be used in objective functions for lung treatments, based on this study it is recommended that the Lyman NTCP model be used because it will provide most consistency with traditional clinician ranking.

Bath and shower effects in the rat parotid gland explain increased relative risk of parotid gland dysfunction after intensity-modulated radiotherapy

PURPOSE

To assess in a rat model whether adding a subtolerance dose in a region adjacent to a high-dose irradiated subvolume of the parotid gland influences its response (bath-and-shower effect).

METHODS AND MATERIALS

Irradiation of the whole, cranial 50%, and/or the caudal 50% of the parotid glands of Wistar rats was performed using 150-MeV protons. To determine suitable (i.e., subtolerance) dose levels for a bath-dose, both whole parotid glands were irradiated with 5 to 25 Gy. Subsequently groups of Wistar rats received 30 Gy to the caudal 50% (shower) and 0 to 10 Gy to the cranial 50% (bath) of both parotid glands. Stimulated saliva flow rate (function) was measured before and up to 240 days after irradiation.

RESULTS

Irradiation of both glands up to a dose of 10 Gy did not result in late loss of function and is thus regarded subtolerance. Addition of a dose bath of 1 to 10 Gy to a high-dose in the caudal 50% of the glands resulted in enhanced function loss.

CONCLUSION

Similar to the spinal cord, the parotid gland demonstrates a bath and shower effect, which may explain the less-than-expected sparing of function after IMRT.

Dose-volume effects in rat thoracolumbar spinal cord: the effects of nonuniform dose distribution

PURPOSE

To investigate dose-volume effects in rat spinal cord irradiated with nonuniform dose distributions and to assess regional differences in radiosensitivity.

METHODS AND MATERIALS

A total of 106 rats divided into three groups were irradiated with (192)Ir gamma-rays at a high dose rate. The groups were irradiated with one, two, or six catheters distributed around the thoracolumbar spinal cord to create different dose distributions. After irradiation, the animals were tested for motor function for 9 months. The response was defined as motor dysfunction and WM or nerve root necrosis. Dose-response data were analyzed with a probit analysis as function of the dose level at a percentage of the volume (D(%)) and with different normal tissue complication probability models. Additionally, the histologic responses of the individual dose voxels were analyzed after registration with the histologic sections.

RESULTS

The probit analysis at D(24) (24% of the volume) gave the best fit results. In addition, the Lyman Kutcher Burman model and the relative seriality model showed acceptable fits, with volume parameters of 0.17 and 0.53, respectively. The histology-based analysis revealed a lower radiosensitivity for the dorsal (50% isoeffective dose [ED(50)] = 32.3) and lateral WM (ED(50) = 33.7 Gy) compared with the dorsal (ED(50) = 25.9 Gy) and ventral nerve roots (ED(50) = 24.1 Gy).

CONCLUSIONS

For this nonuniform irradiation, the spinal cord did not show typical serial behavior. No migration terms were needed for an acceptable fit of the dose-response curves. A higher radiosensitivity for the lumbar nerve roots than for the thoracic WM was found.

Toward improving the therapeutic ratio in stereotactic radiosurgery: selective modulation of the radiation responses of both normal tissues and tumor

✓A review of the radiobiological factors that influence the response of the brain to radiation is provided in relation to stereotactic radiosurgery (SRS). The prospects for intervention after radiation treatment to selectively modulate the expression of late central nervous system (CNS) injury is considered, as well as an account of recent interest in the use of radiation enhancers to selectively increase the response of tumors to radiation.

Brain necrosis in humans, after conventional irradiation, indicates that the risk of necrosis increases rapidly after an equivalent single dose of 12 or 13 Gy. When single-dose treatments are extended due to 60Co decay or planned extension of treatment times, account should be taken of the effects of the repair of sublethal radiation damage to DNA on the efficacy of treatment. Both repair capacity and repair kinetics will also influence tumor control, but parameters to quantify this effect have not yet been established.

The volume of CNS tissue that has been irradiated affects the tissue response, but this effect is only significant for volumes less than 0.05 cm3. The gain obtained from irradiation of small volumes is reduced, however, when focal irradiation is given within a wider field of irradiation.

Based on a vascular hypothesis explaining the pathogenesis of late CNS damage, approaches designed to selectively modulate the frequency of late CNS damage have been validated. Given the high intrinsic radioresistance of some tumors, as opposed to the presence of hypoxia, an interest has developed in the use of selective radiation enhancers in the treatment of tumors. The compound presently available has proved to be disappointing clinically due to toxicity at effective doses, when repeated administration is required. However, when given at high single doses it is less toxic and may be more effective. Less toxic radiation enhancers need to be developed.

Self-consistent tumor control probability and normal tissue complication probability models based on generalized EUDa)

Traditional methods to compute the tumor control probability (TCP) or normal tissue complication probability (NTCP) typically require a heterogeneous radiation dose distribution to be converted into a simple uniform dose distribution with an equivalent biological effect. Several power-law type dose-volume-histogram reduction schemes, particularly Niemierko's generalized equivalent uniform dose model [Med. Phys. 26, 1000 (1999)], have been proposed to achieve this goal. In this study, we carefully examine the mathematical outcome of these schemes. We demonstrate that (1) for tumors, with each tumor cell independently responding to local radiation dose, a closed-form analytical solution for tumor survival fraction and TCP can be obtained; (2) for serial structured normal tissues, an exponential power-law form relating survival to functional sub-unit (FSU) radiation is required, and a closed-form analytical solution for the related NTCP is provided; (3) in the case of a parallel structured normal tissue, when NTCP is determined solely by the number of the surviving FSUs, a mathematical solution is available only when there is a non-zero threshold dose and/or a finite critical dose defining the radiotherapy response. Some discussion is offered for the partial irradiation effect on normal tissues in this category; (4) for normal tissues with alternative architectures, where the radiation response of FSU is inhomogeneous, there is no exact global mathematical solution for SF or NTCP within the available schemes. Finally, numerical fits of our models to some experimental data are also presented.