Computational and Mathematical Methods in Medicine
Volume 2016 (2016), Article ID 2796186, 9 pages
http://dx.doi.org/10.1155/2016/2796186

Review Article

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

Marco D’Andrea, Marcello Benassi, and Lidia Strigari

Laboratory of Medical Physics and Expert Systems, Regina Elena National Cancer Institute, Rome, Italy

Received 14 July 2016; Accepted 23 October 2016

Academic Editor: David A. Winkler

Copyright © 2016 Marco D’Andrea et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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.

1. Introduction

Modern radiotherapy techniques allow unprecedented levels of accuracy, precision, and conformity in target localization, patient setup, and dose delivery thanks to the aid of many different imaging modalities.

Contemporary treatment strategies almost always involve delivering higher doses to the targeted tissue with the aim of improving tumor control, but before such approaches can be safely implemented an accurate and reliable knowledge on toxic effects on surrounding tissues has to be secured.

With the aim of normal tissue preservation many models have been proposed to describe radiation induced complications mostly focusing on late complications which, being irreversible, are considered to have the highest impact on the patient quality of life.

In this, as in most overviews [1–5], normal tissue complication probability (NTCP) models have been divided into mechanistic and (semi)empirical, according to the level of detail in tissue structure that is introduced. Ideally models should be able to accommodate the body of knowledge coming from cellular radiobiology and more specifically the linear quadratic (LQ) model of cell kill and its more sophisticated evolutions that incorporate cell proliferation, cycle effects, and repair as well as local environmental effects and vascularization. The mechanistic models which almost invariably rely on the paradigm on viewing the tissue as a cooperative collection of functional subunits that allow preservation of the tissue functionalities are able to integrate the LQ model in a more straightforward way.

This overview focuses on the description of tissue organization without any initial assumption on the subunit response to radiation. Our approach stresses the mathematical translation of the features of the presented models to better expose their versatility and the opportunities for further developments. Indeed, radiobiological modeling needs a quite complex mathematical toolbox, mirroring the complexity of the biological systems. Each organ is not just an agglomerate of cells but it has an underlying architecture/organization that is the very basis of its functional role, enabling many different strategies (renewal/replacement of damaged cells, intricate microscopic repair pathways, etc.) to successfully deal with radiation damage [6].

Many of these radiobiological models have been integrated into Treatment Planning Systems [7–10] sometimes in a simplified form to allow a biological optimization of the dose delivery or a ranking of competing treatment strategies. This can be seen as a mandatory first step towards a patient specific design of a radiotherapy care path.

2. The Lyman-Kutcher-Burman Model

2.1. The Pure Lyman Model

The most widely used NTCP model in clinical radiobiology is the Lyman model [11].

The core of the model is a fit of frequency data collected for chosen clinical toxicity endpoints and for a particular class of dose irradiation patterns, namely, those patterns where a given portion of an organ or tissue volume absorbs a spatially uniform dose and the rest of its volume absorbs (ideally) no dose at all [12]. If we label as the point set within the irradiated organ where the dose takes the constant dose value , the chosen fitting function is assumed to be given by the following expression [13]:

where is the measure of the support set (assuming the set , representing the whole organ, to have measure 1); , , and are fitting parameters, specific for organ and endpoint , and is the error function, defined as

As the only features of the dose map that are included in the model are its constant value and the measure of its support, it is clear that the actual position and shape of the set within the organ are considered irrelevant.

In more technical terms, if we define a random indicator such that flags the event that endpoint is observed in organ when the latter absorbs the dose , we can say that at this stage we assume to know the probability only for dose maps of the form for any , where is a point in the organ under investigation and is the indicator function of the set ; that is,

This knowledge is summarized by a database of triples .

2.2. The Kutcher-Burman Interpolation

Having described the radiobiological consequences of the partial uniform dose distributions, the model needs to be extended to arbitrary distributions; that is, we want to know the probability for arbitrary . We are thus facing an interpolation/extrapolation problem in the domain of dose maps and the missing information must be fed in, through some further assumption. Many solutions have been proposed to this problem [14, 15] the most widely accepted being the so-called “Kutcher-Burman histogram reduction scheme” [16–20]. This interpolation algorithm relies on the concept of Dose Volume Histogram (DVH) [21, 22], which we now briefly introduce. First of all, to be more adherent to the clinical practice and to avoid dimensional issues, we note that a dose map is usually given in the form where is the prescription dose and is a dimensionless function (the “normalized dose distribution”).

The range of can be divided into, say, intervals (bins) of width , each interval defining its inverse image. To a chosen binning we can associate the set of couples where one representative value, (e.g., the lower boundary or the middle point), in each dose interval is matched to the size of the inverse image, , of that interval. This set of couples is called differential (normalized) DVH (dDVH).

The key idea of the Kutcher-Burman interpolation is to consider each couple of the dDVH under investigation as describing an independent partial uniform distribution to which a new, equivalent (here “equivalent” means being mapped to the same real number by ), partial uniform distribution, , is associated. This new equivalent partial uniform distributions is taken to have a support whose size is fixed by the requirement of preserving the complication probability, that is to say, the value of in (1). In formulas we have

from which we derive

The final assumptions are as follows: (a) consider the supports as disjoint (which we can do in accordance with the assumption that the actual position of the dose absorbing regions is irrelevant) so that the partial uniform distribution can be defined and (b) take the initial distribution as being equivalent to the latter partial uniform distribution for which we are allowed to use formula (1) to calculate a value of and a complication probability. Indeed,

We have that

A further trivial manipulation of the last formula yields

And it is clear that the product in the numerator can be interpreted as a uniform dose delivered to the whole organ volume which has the same complication probability of the original distribution; hence it is named Equivalent Uniform Dose (EUD):

It is easily checked that a change in the normalization for does not change the value of , provided that the prescription dose is scaled accordingly (indeed formula (8) is most often written as for arbitrary ; however note that = − = = − .

Note also that in the above reasoning the values are not dependent on the dose map; therefore they are expected to be those listed in the “elementary” database. Continually updated versions of this body of knowledge can be found in the literature [23], most notably the QUANTEC collaboration [24] which is the reference for daily clinical practice in radiotherapy.

2.3. Issues on Fractionation Effects

When the Lyman model was proposed, patient irradiation had a simple ballistician that delivered dose distributions with sharp penumbras. As a consequence the basic assumption that dose delivery was performed on a dichotomic ‎basis, with irradiated tissue absorbing all the dose at the intended fractionation and the rest of the tissues absorbing no dose at all, was fairly accurate. In a few years, however, technological advances allowed irradiations strategies that, in spite of a much better performance in concentrating the high dose region around the tumour tissue, could induce highly nonuniform low dose deposit in surrounding healthy tissues. This motivated, at least partially, the introduction of the Kutcher-Burman reduction scheme, which then needed to accommodate also the fact that different portions of the irradiated volume receive the accumulated dose in different fraction sizes according to their position in the inhomogeneous dose map. The most used approach to solve this issue is a nonlinear rescaling of the dose axis according to the BED formula of the LQ model [25, 26]. This solution is obtained by forcing microscopic radiobiology into a model that has no built-in slot to account for “microscopic resolution.” However, once this step is accepted, additional microscopic details coming from dedicated in vitro experiments can be included such as incomplete repair, track structure effects [27, 28], and high dose linearization of survival (i.e., the LQL model) [29, 30].

3. The Models Based on Functional Subunits

3.1. The Threshold Poisson-Binomial (TPB) Models: Critical Volume Model and Critical Element Model

The Critical Volume Model [31–33] relies on the identification of an “atomic” biological entity, usually denoted as Functional Subunit (FSU), whose response to ionizing radiation can be described as a Bernoulli trial, with an associated dichotomic random variable (random indicator) , taking the value 1 with probability if the FSU has been “inactivated,” that is, it is unable to perform its biological task after the irradiation, and the value 0 with probability if the FSU is still functioning properly.

It is also assumed that a definite collection of of these atomic entities forms a tissue, that is, a macroscopic, composite, biological entity, whose status after irradiation is evaluated from the observation of some chosen biological endpoint.

It is only from this viewpoint that the FSUs are considered as “atomic” in the sense that all the information about their interaction with ionizing radiation is summarized by the inactivation probability . In fact, to arrive at a consistent expression relating absorbed dose to inactivation probability the internal structure of the FSUs must be considered.

If we define as the random indicator associated with the th FSU, then the sum is a random variable which counts the number of FSUs that have been inactivated by radiation.

In the standard formulation of the model, the FSUs are assumed independent, and the behavior of their macroscopic aggregate, that is, the tissue, is described by the variable just introduced, together with a characteristic threshold . These two can be combined in a new, global, random indicator , where is the Heaviside step function. and are related in such a way that if more than FSUs are inactivated by radiation, that is, , then the tissue , as a whole, becomes unable to perform one or more of its functions and its impaired status, , is flagged by the biological endpoint under observation.

From the assumed independence of the FSUs it follows that if is the probability associated with the variable taking the value 1, then the variable has a probability distribution, frequently named “Poisson-binomial,” which reduces to the usual binomial distribution if all the s are identical. In the following, an arbitrary distribution of inactivation probabilities among the FSUs, that is, the -tuple of probabilities , will be denoted as , a uniform probability distribution will be denoted as , where is the indicator -tuple , and a partial uniform probability distribution will be denoted as , where is the indicator -tuple having the first components equal to 1 and the remaining equal to .

Given , the probability that exactly out of FSUs are inactivated, that is, the Poisson-binomial probability mass function, is

where the summation ranges over all the possible subcollections of out of FSUs and, therefore, has terms.

From expression (10) it is clear that a permutation of the indices of the components of the vector leaves the value of unchanged. Since the index locates a particular FSU in space, this reflects the fact that independence of the FSUs wipes out any global spatial feature of the tissue under investigation.

A different expression is often used for (1) by first summing over the FSUs that share the same value of inactivation probability [34, 35]. The components of the vector are therefore clustered in, say, probability bins, with , each bin being referred to FSUs having the common inactivation probability value and, of course, . Expression (10) can thus be rewritten as

that is, as a constrained summation of products of binomial probability mass functions, each describing the inactivation of out of FSUs in the th probability bin. From (11) we see that if each bin refers to exactly one FSU, that is, and , we obtain again expression (10); if, on the other hand, there is just one bin, that is, we have a uniform probability distribution, then we obtain the usual binomial probability mass function:

A closed form expression for (11) can be found, for example, in [34]. Discrete Fourier Transform techniques can be used to obtain a closed form also for the general expression (10) [36].

From (10) or (11) it can be shown that [37]

With the usual notation for the sample mean of an -tuple . Introducing the (biased) sample variance it follows from (11) that

which shows how the variance of increases with increasing uniformity among the inactivation probabilities and attains its maximum value when all the s are identical.

The probability that at most FSUs are inactivated defines the cumulative distribution function:

The complement to 1 of this function gives the probability that at least FSUs are inactivated:

When is equal to the tissue characteristic threshold , we are looking at the event that we identify with the loss of some functionality in the tissue, which shows up in the observation of the related biological endpoint, and the corresponding probability could be called “tissue failure probability” .

The customary interpretation of the described probability distributions and the related nomenclature is such that when we are actually looking at healthy tissue or an organ being damaged by ionizing radiation, the probabilities can span the whole range , and the probability is usually called NTCP. The lower extreme, , describes a tissue where damage to a single FSU is enough to trigger the endpoint. This scenario actually has its own model, namely, the so-called Critical Element model, which, although being a special case of the Critical Volume, was introduced somewhat before the latter [38–41]. At the upper extreme , instead, we are looking at a tumour undergoing therapeutic irradiation, the probabilities are, hopefully, all close to 1, and the probability is called “Tumour Control Probability” .

The number of FSUs is usually very large and approximate expressions for (15) that are easier to handle can be used, the normal and the Poisson one being the most popular. The normal approximation stems from the well-known DeMoivre-Laplace limit theorem generalized to the Poisson-binomial distribution. Indeed, if the vector is such that [34] we have that

where and are given by formulas (13). On the other hand, for the case of therapeutic irradiation of a tumour since , with , the variable which counts the number of FSUs still active has approximately a Poisson distribution whose zero class identifies the event of total FSU inactivation; that is,

where and .

It can be shown [42–44] that in (10) is well approximated also by the following:

which is a Poisson-binomial probability mass function that, from the viewpoint of representation (11), has only bins: , , , and .

and are defined by the system of equations

whose solution, taking into account (13) and the fact that must be an integer, is

with being the floor function which outputs the integer part of a real number. Using the sample mean and the sample variance equations (21) can be rewritten as

This approximation also works for the cumulative distribution (15) for any and therefore, in particular, for . Since is less than 1, a notion of “effective fraction,” , of irradiated FSUs comes naturally into play and its value is fixed by the inactivation probabilities according to formula (22). The “equivalent uniform probability,” , over the effective fraction of irradiated FSUs can be used to define an Equivalent Uniform Dose via the specification of how the inactivation probability is related to dose, provided that appropriate uniformity assumptions are introduced, as will be explained in the following.

If we introduce the random variable describing the relative amount of inactivated FSUs together with the relative threshold the approximation (17) becomes

and . Using (17) we have

This formula can be rewritten using the error function defined as

and the result is

A number of FSUs large enough to enforce the normal approximation (17) turn the Gaussian integrand in (25) into a -like distribution and we see that

therefore tissue failure takes place (with probability one) whenever

that is,

In other words the law of large numbers further erases the details of the distribution of the s and the only remaining feature is the arithmetic mean .

Even if no specification other than monotonic increase is made about the functional dependence of the FSU inactivation probability on radiation dose, it is clear that, upon choice of a given tissue, that is, a value of , formula (28) predicts a sharp transition in as a function of dose which is not observed in clinical data. This could be amended by going back to formula (25)–(27), but the fact that is sufficiently large to strongly concentrate the Gaussian density function around its expectation value seems hardly questionable.

To explain the gentler slope of it is necessary to recognize the sampling character of the clinical toxicity frequency data that we are investigating. Indeed, choosing a tissue does not specify a value of , but only its distribution among the population of patients. Given the symmetrical roles of and in (28), the same idea works also when a distribution of radiobiological parameters is assumed among the patients so that a given radiation dose pattern, , does not yield a definite value of and, therefore, of but, again, a distribution of values. It may be worth pointing out that heterogeneity of radiobiological parameters among the FSUs of a single patient is already taken into account by the assumed heterogeneity of the s and it cannot prevent the sharp transition when the number of FSUs becomes very large.

In formulas, if we rewrite (28) as

where the notation has now been used in replacement of to parallel the treatment of the LKB model, the above reasoning amounts to substitute with a density distribution for the relative threshold and with a (dose pattern dependent) density distribution for . The resulting general expression, therefore, becomes

Going back one step by leaving and focusing only on the distribution the above formula turns into

If we denote by the cumulative distribution of and integrate (31) by parts, we obtain

Therefore, varying the radiation dose pattern actually samples the cumulative distribution of the relative threshold through the mean value of the inactivation probability or, in simpler words, the probability of undergoing a tissue failure, given , is just the probability of having a relative threshold smaller than .

3.2. FSU Inactivation Probability

Having shown the assumptions that regulate the organizational response of tissue in the Critical Volume Model and discussed their consequences we must now specify the dose dependence of the FSU inactivation probabilities .

The simplest approach, which however somehow spoils the mechanistic attitude of the CV model, is to assume a generic sigmoid shaped dose response curve whose parameters may or may not depend on the position of the FSU within the organ. In the former case we have some means to introduce a spatial feature within the organ, so that the position and shape of the irradiated volume become relevant.

A more mechanistically oriented assumption is the replication of the Critical Volume approach at a smaller scale: the FSU is considered as an aggregate of independent stem cells and FSU inactivation occurs when a sufficiently large fraction, , of these stem cells has been killed by radiation [45, 46]. To complete the picture we can invoke the LQ model to describe the cell killing dependence on radiation dose and we may assume that even a single surviving stem cell can allow the FSU to regenerate; that is, that :

where the subscript , which, as already pointed out, locates the FSU in space, has been attached also to the radiobiological parameters to allow for their position dependence within the organ.

3.3. Nesting TPBs: Kallman’s Model

The idea of nesting one Poisson-binomial probability distribution into another allows modeling increasing complexity of tissue organization.

At the level of complexity immediately next to the Critical Volume Model, an organ is viewed as a bundle of “meta-FSUs” that can be damaged according to the CV model with a threshold, and each meta-FSU is itself made of FSUs and again follows the CV model with a threshold.

Single FSU inactivation can be described, for example, by expression (35).

In formulas for the th meta-FSU, borrowing from (10), (15), and (16), this translates into

And for the whole organ under investigation

where is the double-indexed vector of FSU inactivation probabilities and and are the subcollections of FSUs out of and meta-FSUs out of , respectively. This general formula is usually restricted to a pure parallel behavior of the organ with respect to the meta-FSUs, [47, 48], that is, , and a pure serial behavior of each meta-FSU, that is, . Therefore we have

As it can be seen by direct inspection of formulas (36)-(37) or (38) there is no invariance for arbitrary index permutation, which means that this model has some built-in global spatial feature. Indeed even if no spatial dependence is assumed for the parameters that regulate FSU response to radiation dose, different positions and shapes of the irradiated volume with respect to the meta-FSU/FSU arrangement yield different values of NTCP. This means that to compare the model to clinical data one needs to take into account spatial details which are most often not available; therefore a simplified version of the model with forced position invariance is used. Such a model can be inspired by (38) after a few observations. First of all, when the inactivation probabilities are assumed identical among the subunits, formula (38) turns into

Which can be trivially solved for

Now we focus on a subset of meta-FSUs, each one limited to only FSUs and, by analogy with (39), we assume that the failure probability of this subset of FSUs is

Having used (40) in the second step. This “tile” of tissue with relative volume is now by definition our new elementary FSU and we introduce a new parameter termed the “relative seriality” :

which is the ratio of the number of the “old” FSUs contained in a meta-FSU over the total number of “old” FSUs. With these new definitions formula (40) can be rewritten as

and allowing the s to have spatial dependence, for an organ comprising “tiles,” we have

where, as before, labels a position within the organ under investigation. In summary we got rid of the difficulties of formula (38) at the expense of collapsing the tissue complexity into an elementary “tile” which has become the new FSU. The usual binning of the inactivation probabilities yields

where is the relative size of the th dose bin.

4. Conclusions

The main aim of the LKB model is to accommodate the available clinical data into a reasonably manageable function, in order to help the clinician in assessing the odds of safe and successful treatment. As long as the clinical needs are involved, any understanding of the underlying biological phenomena is pursued only as long as it is instrumental to the above purpose. From the viewpoint of detailed radiobiological modeling, however, this model can be taken as a summary of the available clinical knowledge.

Quite general features of a normal tissue complication probability model which are already present in the LKB model are the following.

The observation of endpoint in organ is an event that can be encoded in an indexed family of random indicators , each described by a Bernoulli distribution. The set of dose maps is the index set for this family.

The only parameter of the Bernoulli distribution, that is, the probability of occurrence of the event, is given by a functional defined on . In particular in the LKB model the functional is of the form

with , and . Noteworthy is the composition , that is, the “histogram reduction,” which is responsible for the dramatic reduction of the dimensionality of the problem.

At this stage, then, a rather featureless sample space is enough to describe experimental data. The TPB models add some details to the sample space and related probability structure. The cost of trying to magnify the microscopic biological features is that the experimental noise, that is, the conditions in which data are gathered, has to be taken into account too, and at the same level of detail.

The mechanistic models, here exemplified by the TPB models, have an intrinsic ability to include microscopic biological features of the radiation induced damage; however the parameters needed to describe such details have not been determined for most treatment sites and thus cannot be used in daily clinical practice. Accurate determination of fine structure parameters in a clinical setting is quite a formidable task, even when foreseeable biasing is taken into account in the model framework. In addition, even the QUANTEC parameter database for the widely used LKB model suffers from indeterminacies due to nonuniformity of volume definitions (e.g., hollow organs) and heterogeneity in the quality of data and in the radiobiological assumptions

All these models and the proposed parameters summarizing the features of a specific tissue/organ have originated in the years of two-dimensional and three-dimensional (conformal) radiotherapy. Nevertheless they currently are the only available approaches to predict the expected toxicities deriving from the modern intensity modulated delivery techniques such as IMRT and VMAT.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contributions

Marco D’Andrea, Marcello Benassi, and Lidia Strigari contributed equally to this work.

References

V. Moiseenko, J. Battista, and J. Van Dyk, “Normal tissue complication probabilities: dependence on choice of biological model and dose-volume histogram reduction scheme,” International Journal of Radiation Oncology Biology Physics, vol. 46, no. 4, pp. 983–993, 2000.

E. D. Yorke and D. Ellen, “Modeling the effects of inhomogeneous dose distributions in normal tissues,” Seminars in Radiation Oncology, vol. 11, no. 3, pp. 197–209, 2001.

M. E. P. Philippens, L. A. M. Pop, A. G. Visser, S. A. M. W. Schellekens, and A. J. van der Kogel, “Dose-volume effects in rat thoracolumbar spinal cord: an evaluation of NTCP models,” International Journal of Radiation Oncology Biology Physics, vol. 60, no. 2, pp. 578–590, 2004.

P. Van Luijk, H. P. Bijl, A. W. T. Konings, A. J. Van der Kogel, and J. M. Schippers, “Data on dose-volume effects in the rat spinal cord do not support existing NTCP models,” International Journal of Radiation Oncology Biology Physics, vol. 61, no. 3, pp. 892–900, 2005.

M. Adamus-Górka, P. Mavroidis, B. K. Lind, and A. Brahme, “Comparison of dose response models for predicting normal tissue complications from cancer radiotherapy: application in rat spinal cord,” Cancers, vol. 3, no. 2, pp. 2421–2443, 2011.

E. L. Travis and L. Elizabeth, “Organizational response of normal tissues to irradiation,” Seminars in Radiation Oncology, vol. 11, no. 3, pp. 184–196, 2001.

Q. Wu, D. Djajaputra, Y. Wu, J. Zhou, H. H. Liu, and R. Mohan, “Intensity-modulated radiotherapy optimization with gEUD-guided dose-volume objectives,” Physics in Medicine and Biology, vol. 48, no. 3, pp. 279–291, 2003.

A. Ólafsson, R. Jeraj, and S. J. Wright, “Optimization of intensity-modulated radiation therapy with biological objectives,” Physics in Medicine and Biology, vol. 50, no. 22, pp. 5357–5379, 2005.

S. Das, “A role for biological optimization within the current treatment planning paradigm,” Medical Physics, vol. 36, no. 10, pp. 4672–4682, 2009.

X. Allen Li, M. Alber, J. O. Deasy et al., “The use and QA of biologically related models for treatment planning: short report of the TG-166 of the therapy physics committee of the AAPM,” Medical Physics, vol. 39, no. 3, pp. 1386–1408, 2012.

J. L. Lyman, “Complication probability as assessed from dose-volume histograms,” Radiation Research, vol. 104, pp. S13–S19, 1985.

B. Emami, J. Lyman, A. Brown et al., “Tolerance of normal tissue to therapeutic irradiation,” International Journal of Radiation Oncology, Biology, Physics, vol. 21, no. 1, pp. 109–122, 1991.

C. Burman, G. J. Kutcher, B. Emami, and M. Goitein, “Fitting of normal tissue tolerance data to an analytic function,” International Journal of Radiation Oncology, Biology, Physics, vol. 21, no. 1, pp. 123–135, 1991.

J. T. Lyman and A. B. Wolbarst, “Optimization of radiation therapy, III: a method of assessing complication probabilities from dose-volume histograms,” International Journal of Radiation Oncology, Biology, Physics, vol. 13, no. 1, pp. 103–109, 1987.

J. T. Lyman and A. B. Wolbarst, “Optimization of radiation therapy, IV: a dose-volume histogram reduction algorithm,” International Journal of Radiation Oncology, Biology, Physics, vol. 17, no. 2, pp. 433–436, 1989.

G. Kutcher and C. Burman, “Calculation of complication probability factors for non-uniform normal tissue irradiation: the effective volume method,” International Journal of Radiation Oncology, Biology, Physics, vol. 16, no. 6, pp. 1623–1630, 1989.

G. J. Kutcher, C. Burman, L. Brewster, M. Goitein, and R. Mohan, “Histogram reduction method for calculating complication probabilities for three-dimensional treatment planning evaluations,” International Journal of Radiation Oncology, Biology, Physics, vol. 21, no. 1, pp. 137–146, 1991.

R. Mohan, G. Mageras, L. Brewster, G. Kutcher, S. Leibel, C. Burman et al., “Clinically relevant optimization of 3-D cinformal treatments,” Medical Physics, vol. 19, no. 4, pp. 933–944, 1992.

J. O. Deasy, A. Eisbruch, and R. K. Ten Haken, “Comments on the use of the Lyman-Kutcher-Burman model to describe tissue response to nonuniform irradiation,” International Journal of Radiation Oncology Biology Physics, vol. 47, no. 5, pp. 1458–1460, 2000.

G. Luxton, P. J. Keall, and C. R. King, “A new formula for normal tissue complication probability (NTCP) as a function of equivalent uniform dose (EUD),” Physics in Medicine and Biology, vol. 53, no. 1, pp. 23–36, 2008.

R. E. Drzymala, R. Mohan, L. Brewster et al., “Dose-volume histograms,” International Journal of Radiation Oncology, Biology, Physics, vol. 21, no. 1, pp. 71–78, 1991.

A. Niemierko and M. Goitein, “Dose-volume distributions: a new approach to dose-volume histograms in three-dimensional treatment planning,” Medical Physics, vol. 21, no. 1, pp. 3–11, 1994.

F. Kong, C. Pan, A. Eisbruch, and R. K. T. Haken, “Physical models and simpler dosimetric descriptors of radiation late toxicity,” Seminars in Radiation Oncology, vol. 17, no. 2, pp. 108–120, 2007.

S. M. Bentzen, L. S. Constine, J. O. Deasy et al., “Quantitative Analyses of Normal Tissue Effects in the Clinic (QUANTEC): an Introduction to the Scientific Issues,” International Journal of Radiation Oncology Biology Physics, vol. 76, no. 3, pp. S3–S9, 2010.

T. E. Wheldon, C. Deehan, E. G. Wheldon, and A. Barrett, “The linear-quadratic transformation of dose-volume histograms in fractionated radiotherapy,” Radiotherapy and Oncology, vol. 46, no. 3, pp. 285–295, 1998.

C. S. Park, Y. Kim, K. M. Bucci et al., “Method to account for dose fractionation in analysis of IMRT plans: modified equivalent uniform dose,” International Journal of Radiation Oncology, Biology, Physics, vol. 62, no. 3, pp. 925–932, 2005.

M. Scholz and G. Kraft, “Calculation of heavy ion inactivation probabilities based on track structure, X ray sensitivity and target size,” Radiation Protection Dosimetry, vol. 52, no. 1–4, pp. 29–33, 1994.

R. B. Hawkins, “A statistical theory of cell killing by radiation of varying linear energy transfer,” Radiation Research, vol. 140, no. 3, pp. 366–374, 1994.

M. Guerrero and X. A. Li, “Extending the linear-quadratic model for large fraction doses pertinent to stereotactic radiotherapy,” Physics in Medicine and Biology, vol. 49, no. 20, pp. 4825–4835, 2004.

C. Park, L. Papiez, S. Zhang, M. Story, and R. D. Timmerman, “Universal survival curve and single fraction equivalent dose: useful tools in understanding potency of ablative radiotherapy,” International Journal of Radiation Oncology Biology Physics, vol. 70, no. 3, pp. 847–852, 2008.

A. B. Wolbarst, L. M. Chin, and G. K. Svensson, “Optimization of radiation therapy: integral-response of a model biological system,” International Journal of Radiation Oncology, Biology, Physics, vol. 8, no. 10, pp. 1761–1769, 1982.

A. Niemierko and M. Goitein, “Modeling of normal tissue response to radiation: the critical volume model,” International Journal of Radiation Oncology, Biology, Physics, vol. 25, no. 1, pp. 135–145, 1993.

E. D. Yorke, G. J. Kutcher, A. Jackson, and C. C. Ling, “Probability of radiation-induced complications in normal tissues with parallel architecture under conditions of uniform whole or partial organ irradiation,” Radiotherapy and Oncology, vol. 26, no. 3, pp. 226–237, 1993.

A. Jackson, G. J. Kutcher, and E. D. Yorke, “Probability of radiation-induced complications for normal tissues with parallel architecture subject to non-uniform irradiation,” Medical Physics, vol. 20, no. 3, pp. 613–625, 1993.

P. Stavrev, N. Stavreva, A. Niemierko, and M. Goitein, “Generalization of a model of tissue response to radiation based on the idea of functional subunits and binomial statistics,” Physics in Medicine and Biology, vol. 46, no. 5, pp. 1501–1518, 2001.

M. Fernández and S. Williams, “Closed-form expression for the poisson-binomial probability density function,” IEEE Transactions on Aerospace and Electronic Systems, vol. 46, no. 2, pp. 803–817, 2010.

W. Feller, Introduction to Probability Theory and Its Applications, vol. 1, John Wiley & Sons, New York, NY, USA, 1968.

T. E. Schultheiss, C. G. Orton, and R. A. Peck, “Models in radiotherapy: volume effects,” Medical Physics, vol. 10, no. 4, pp. 410–415, 1983.

A. B. Wolbarst, “Optimization of radiation therapy II: the critical-voxel model,” International Journal of Radiation Oncology, Biology, Physics, vol. 10, no. 5, pp. 741–745, 1984.

H. Rodney Withers, J. M. G. Taylor, and B. Maciejewski, “Treatment volume and tissue tolerance,” International Journal of Radiation Oncology, Biology, Physics, vol. 14, no. 4, pp. 751–759, 1988.

A. Niemierko and M. Goitein, “Calculation of normal tissue complication probability and dose-volume histogram reduction schemes for tissues with a critical element architecture,” Radiotherapy and Oncology, vol. 20, no. 3, pp. 166–176, 1991.

B. Roos, “Binomial approximation to the poisson binomial distribution. The Krawtchouk Expansion,” Theory of Probability and Its Applications, vol. 45, no. 2, pp. 258–272, 2001.

K. P. Choi and A. Xia, “Approximating the number of successes in independent trials: binomial versus poisson,” Annals of Applied Probability, vol. 12, no. 4, pp. 1139–1148, 2002.

S. Y. Soon, “Binomial approximation for dependent indicators,” Statistica Sinica, vol. 6, no. 3, pp. 703–714, 1996.

J. H. Hendry and H. D. Thames, “The tissue-rescuing unit,” The British Journal of Radiology, vol. 59, no. 702, pp. 628–630, 1986.

E. Rutkowska, C. Baker, and A. Nahum, “Mechanistic simulation of normal-tissue damage in radiotherapy—implications for dose-volume analyses,” Physics in Medicine and Biology, vol. 55, no. 8, pp. 2121–2136, 2010.

P. Källman, A. Ågren, and A. Brahme, “Tumour and normal tissue responses to fractionated non-uniform dose delivery,” International Journal of Radiation Biology, vol. 62, no. 2, pp. 249–262, 1992.

P. Kallman, B. K. Lind, and A. Brahme, “An algorithm for maximizing the probability of complication-free tumour control in radiation therapy,” Physics in Medicine and Biology, vol. 37, no. 4, pp. 871–890, 1992.