Prostate CancerVolume 2011 (2011), Article ID 128360, 9 pagesdoi:10.1155/2011/128360

Research Article

An NTCP Analysis of Urethral Complications from Low Doserate Mono- and Bi-Radionuclide Brachytherapy

V. E. Nuttens,¹ A. E. Nahum,² and S. Lucas¹

¹NAmur Research Institute for LIfe Sciences (NARILIS), Research Center for the Physics of Matter and Radiation (PMR-LARN), University of Namur (FUNDP), Rue de Bruxelles, 61, 5000 Namur, Belgium²Department of Physics, Clatterbridge Center for Oncology, Clatterbridge Road Bebington, Merseyside CH63 4JY, UK

Received 13 January 2011; Accepted 2 May 2011

Academic Editor: M. J. Zelefsky

Copyright © 2011 V. E. Nuttens 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.

Urethral NTCP has been determined for three prostates implanted with seeds based on ¹²⁵I (145 Gy), ¹⁰³Pd (125 Gy), ¹³¹Cs (115 Gy), ¹⁰³Pd-¹²⁵I (145 Gy), or ¹⁰³Pd-¹³¹Cs (115 Gy or 130 Gy). First, , meaning that 20% of the urhral volume receive a dose of at least , is converted into an I-125 LDR equivalent in order to use the urethral NTCP model. Second, the propagation of uncertainties through the steps in the NTCP calculation was assessed in order to identify the parameters responsible for large data uncertainties. Two sets of radiobiological parameters were studied. The NTCP results all fall in the 19%–23% range and are associated with large uncertainties, making the comparison difficult. Depending on the dataset chosen, the ranking of NTCP values among the six seed implants studied changes. Moreover, the large uncertainties on the fitting parameters of the urethral NTCP model result in large uncertainty on the NTCP value. In conclusion, the use of NTCP model for permanent brachytherapy is feasible but it is essential that the uncertainties on the parameters in the model be reduced.

1. Introduction

Radiotherapy treatment planning systems incorporating “biological” models are beginning to make their way into clinical use. The biological models in questions are for Tumor Control Probability (TCP) [3, 4] and Normal Tissue Complication Probability (NTCP) [5, 6]. A common suggestion is that treatment plans be optimized to maximize the TCP while not exceeding a fixed, acceptable NTCP [7, 8]. A considerable amount of work has been done to implement this biological model-based approach in permanent seed prostate brachytherapy [1, 9–12]. Prostate brachytherapy morbidity is generally reported for the urethra [13–17] and rectum [18–21]. Zaider et al. have extended radiotherapy NTCP models for these organs to Low Dose Rate (LDR) permanent brachytherapy [1, 12].

Three radionuclides are generally used in permanent prostate brachytherapy: iodine-125, palladium-103, and cesium-131. Each radionuclide has its advantages but also its drawbacks. ¹²⁵I and ¹³¹Cs have similar emission spectrum whose mean energies are 28.37 and 30.45 keV, respectively. ¹⁰³Pd has a mean emitted energy of 20.74 keV which reduces the dose to the surrounding organs at risk (OAR) such as rectum and urethra but also increases the risk of cold spots (underdosage) in the prostate. By contrast, the ¹²⁵I and ¹³¹Cs dose distributions extend to larger distances, thus reducing the likelihood of cold spots at the cost of delivering more to the OARs for a given prostate dose. Urethral and rectal complications are reported at a similar frequency for ¹²⁵I or ¹³¹Cs [22].

The contrast between the properties of ¹⁰³Pd, ¹²⁵I, and ¹³¹Cs has motivated our research in developing a new kind of seed based on a mixture of two radionuclides, namely, ¹⁰³Pd0.75-¹²⁵I0.25 or ¹⁰³Pd0.25-¹³¹Cs0.75. The subscripts denote the fractions of internal activity of each radionuclide as defined in a previous paper [23]. To avoid a cumbersome style, the ¹⁰³Pd0.75-¹²⁵I0.25 and ¹⁰³Pd0.25-¹³¹Cs0.75 mixture will be referred to in the text as Pd-I and Pd-Cs. The dosimetry characteristics and prescription doses of these sources were derived in previous studies [23, 24].

In our study, we use the Zaider et al. NTCP model for the urethra to compare a monoradionuclide seed implant (¹⁰³Pd, ¹²⁵I or ¹³¹Cs) and bi-radionuclide seed implants (Pd-I or Pd-Cs). A sensitivity analysis on the modeling parameters is also performed.

2. Patients and Methods

2.1. Treatment Plans

A source data file for the Prostate Seed Implant Dosimetry (PSID) Treatment Planning System (TPS) has been generated for InterSource seeds loaded with either ¹⁰³Pd, ¹²⁵I, ¹³¹Cs, Pd-I, or Pd-Cs following the AAPM TG43U1 formalism, which we have adapted for multiple radionuclide brachytherapy sources [23].

Three patients were planned using peripheral seed placement: patient 1 with the smallest prostate, patient 2 with a medium-sized prostate, and patient 3 with the largest prostate. Prostate and urethra volumes of each patient are given in Table 1. The urethral volume is defined as the volume enclosed by the urethra surface. Each patient has been planned with the five types of seeds. The prescription doses used for bi-radionuclide implants have been derived previously [24] and are summarized in Table 2. The method that we used yields a single, fixed value of 142 Gy for the prescription dose of a Pd-I implant whereas it results in two different values for the Pd-Cs mixture: 115 Gy and 128 Gy. However, these values are associated with large uncertainties. They were rounded to 145 Gy, 115 Gy, and 130 Gy, respectively, as safety margins on tumor control. The total number of treatment plans is therefore 18 (6 plans per patient).

Table 1: Prostate and urethra volumes for each planned patient.

Patient
p01
p02
p03

Prostate volume (cc)
24.64
38.9
45.33

Urethra volume (cc)
0.32
0.19
0.28

Table 2: Prescription dose for the different seeds used for planning.

¹⁰³Pd
¹²⁵I
¹³¹Cs
Pd-I
Pd-Cs

Internal activity
100%
100%
100%
75%–25%
25%–75%

Prescription dose (Gy)
125
145
115
145
130 (a)

115 (b)

2.2. Urethra NTCP

To the best of our knowledge, only Zaider et al. [1] have developed a model for the urethra NTCP after I-125 LDR. This model is based on a correlation between the probability of urethral toxicity and the dose received by the “hottest” 20% of the urethral volume after permanent prostate brachytherapy using ¹²⁵I seeds. In order to apply such a model, and because in our study we are dealing with different radionuclides, we have converted the dose distribution for each implant to an equivalent I-125 LDR one. In order to do this, we derive the of an I-125 LDR treatment that would yield the same Biologically Effective Dose (BED) as the implant in question. The BED includes a linear and a quadratic term in dose. In the case of bi-radionuclide implant, its formulation requires the contribution of each radionuclide to the . The models and procedure are described in what follows.

2.2.1. Logistic Model

Based on a logistic regression analysis of patients treated by I-125 LDR, Zaider et al. [1] inferred the probability of unresolved Grade 2 or higher urethral toxicity at 12 months as a function of :

where and are two fitting parameters given with their uncertainties in Table 3. Further in the text, iodine-125 referring to the I-125 LDR equivalent implant will be noted I-125 whereas iodine-125 referring to our implants will be noted ¹²⁵I.

Table 3: Modeling parameters and their absolute uncertainties. is the value of the considered parameter.

Parameter
Units

—
−2.60
−3.10
−2.10

Gy
0.0066
0.0050
0.0082

h⁻¹
1.0
0
1.5

Gy
16.6
5
25

—
1.41
1.00
1.75

—
1.28
1.00
1.45

—
1.28
1.00
1.45

—
Variable

2.2.2. I-125 LDR Equivalent DU₂₀

The parameters of the urethral NTCP model correspond to an I-125 LDR treatment. Hence, the for the different other radionuclides or mixtures has to be converted to an I-125 LDR equivalent . This conversion will be referred to below as “I-125 conversion”. The Biologically Effective Dose (BED) is a powerful tool for this purpose as it allows a comparison of different treatment modalities.

The BED for an LDR implant whose seeds contain one radionuclide is [2, 25] where is the Relative Biological Effectiveness of the radionuclide at very low-dose rate; [2]; is the ratio of radiosensitivity parameters; is the radioactive decay constant of the radionuclide; is the sublethal cell damage repair rate. For the urethra, the commonly used values for and are 3 Gy and 0.5 h⁻¹ respectively. For normal tissues, cell repopulation is not taken into account. The values used in this study are from Wang et al. [26]: 1.41 for ¹⁰³Pd and 1.28 for ¹²⁵I. Due to the lack of experimental data on the of ¹³¹Cs, its value has been set to the same value (i.e., 1.28) as that for ¹²⁵I. This assumption is based on the fact that their emission spectra are similar. Note that (2) differs from the expression in Zaider et al. [1] one. However, as the sublethal cell damage repair rate for the urethra is much larger than the radioactive decay constant, both results fit within 1% in the worst case.

The clinical outcome would be equivalent for two treatments if they yield the same BED: The I-125 LDR equivalent of the resulting from an implant based on radionuclide can be calculated by expressing (4) as a quadratic where

As , and must have positive values, the only valid solution of (5) is

For a bi-radionuclide seed implant, the solution is the same as (9) with same expression for and , but the expression of is where is the one of the bi-radionuclide implant and is the contribution of radionuclide in the seed to the .

2.2.3. Bi-radionuclide Case: Contributions

Equation (10) requires that the contribution of each radionuclide to be known for each point in space that is considered. Treatment plans show that the relative contribution of each radionuclide can be considered as constant throughout the urethral volume. Hence, everywhere in the urethra: Let us include the RBE effect directly in the TPS file that defines the seed data (dose distribution in the adapted AAPM TG43U1 dosimetry formalism [23]). This can be done by modifying the adapted radial dose function and 1D-anisotropy function of the seed (see Appendix). The dose distribution in the organs provided by the TPS would therefore also include the RBE effect:

As the ’s are spatially independent, (13) can be applied directly to the value. Therefore, if the from the treatment plan with and without RBE effect ( and , resp.,) are known, then the contribution of each radionuclide to can be computed from

These values turn out to match within 1% the mean ’s that one would obtain at the middle of the urethra contour on each TRUS image slice. This result is not surprising as the urethra is approximately in the middle of the prostate and therefore in the middle of the seed distribution. All the seeds of the implant contribute with different magnitudes to urethral dose, providing a homogeneous dose distribution throughout the urethra for each radionuclide.

2.2.4. Propagation of Uncertainty

The method described in the above three sections includes many parameters associated with uncertainties. The uncertainties in come from the radiobiological parameters ( and the of each radionuclide. These uncertainties along with those of the NTCP fitting parameters will affect the urethral NTCP. The NTCP fitting parameters’ uncertainties will not affect the relative NTCPs of the different seed implants as the fitting parameters are not radionuclide dependent. By contrast, parameters are radionuclide dependent. Therefore, the repercussion of the parameters’ uncertainties on and the subsequent repercussion on urethral NTCP were studied.

The uncertainty can be expressed as follows:

The uncertainty on RBE values contains three different terms: the uncertainty on I-125 RBE and on the RBE of each radionuclide in the seed. Note that in the case of mono-radionuclide implants, and . If the implant uses ¹²⁵I or ¹³¹Cs, their value has to be the same as the I-125 value for consistency reasons. RBE uncertainties are therefore correlated and the absolute value has to be taken over the whole RBE uncertainties. Moreover, the palladium-103 RBE is correlated to the iodine-125 RBE. Depending on the study [26, 27], if a high value of iodine-125 (1.45 instead of 1.28) is chosen then a high value of palladium-103 is obtained (1.75 instead of 1.41). As a result, the uncertainties on ¹⁰³Pd RBE are correlated to the uncertainties on iodine-125 RBE, which can be rewritten as

The term into bracket has been calculated using the above mentioned value of ¹⁰³Pd and ¹²⁵I values and is equal to 1.65.

The NTCP uncertainty expression includes only three terms: The partial derivatives of and NTCP are given in Appendices B and C, respectively.

3. Results

3.1. Modeling Parameters and Associated Uncertainties

The radiobiological parameters used in this study are limited to the ratio and sublethal cell damage repair rate, . First, we choose the values that Zaider et al. used in their study. These are the one reported in Table 3. This high ratio is typical for tumors (typical range 5–25 Gy) but does not correspond to the commonly used value for late responding normal tissue (typical range 1–5 Gy) [2]. Therefore, we also calculate the and NTCP with  Gy (1–5 Gy) and  h⁻¹ (0–1.5 h⁻¹) reported in the Dale and Jones text book [2]. The values into brackets are the uncertainty range used for the uncertainty calculations.

Relative Biological Effectiveness values are the one published by Wang et al. [26]. Due to the lack of experimental data on the of ¹³¹Cs, its value has been set to the same value as the one of ¹²⁵I. The uncertainty interval is based on RBE values published by other authors. The minimum value is 1 by definition. The maximum values have been set to the ones published by Antipas et al. [27].

Finally, the uncertainty on , the contribution of the first radionuclide to the , is equal to the maximal deviation with respect to the mean value of the obtained in the middle of the urethra on TRUS slices.

3.2. DU₂₀ and I-125 LDR Equivalent DU₂₀, DU20-I

The TPS output for the different patients and radionuclides are presented in Figure 1.

Figure 1: for the different seeds (mono- and bi-radionuclide) and patients. The averaged over the three patients is also plotted.

The mean value is also shown. As discussed above, these values are not comparable as they do not correspond to the same delivery scheme. Nevertheless, it will be interesting to observe how much the is modified for each seed implant after the I-125 conversion. It is emphasized that the does not change significantly from one patient to another one.

The implants can be compared after the I-125 conversion. The of each seed implant has been converted to an I-125 LDR equivalent using (6) to (10) and (14). Figure 2 shows the results and the numerical values are reported in Table 4.

Table 4: Dose received by the hottest 20% of the urethra volume (DU₂₀), 125I LDR equivalent DU₂₀ (DU20-I), and urethra Normal Tissue Complication Probability (NTCP) for the different seed implants. The values are averaged over the three patients. Radiobiological parameters from Zaider et al. are used.

(Gy)
(Gy)
NTCP (%)

¹²⁵I
202.3
202.3
22.0

Pd-I
201.8
207.4
22.6

¹⁰³Pd
176.3
195.8
21.3

Pd-Cs (a)
185.5
194.9
21.2

Pd-Cs (b)
169.0
177.3
19.3

¹³¹Cs
168.4
171.7
18.7

Figure 2: I-125 LDR equivalent results for the different implants (mono- and bi-radionuclide) and patients. The mean over the three patients is also plotted. Zaider et al. radiobiological parameters are used. Error bars were calculated using (15).

The effect of this conversion is to increase the . The smallest increase is observed for the ¹³¹Cs implant with 3.3 Gy and the largest one for the ¹⁰³Pd seed implant with 19.5 Gy. The uncertainty associated with these values will be the subject of a separate point. As expected, ¹²⁵I is not affected by the I-125 conversion. The ¹³¹Cs implant is also not significantly affected by the I-125 conversion which suggests that the quadratic contribution of BED is small.

3.3. Urethral NTCP

Although the figures correspond to the same delivery scheme (I-125 LDR) and they can be used to compare treatments, they do not give information about the probability of urethral complications. These probabilities have been computed using (1) and are plotted in Figure 3. Numerical values are reported in Table 4.

Figure 3: Urethral NTCP results for each patient and each implant (mono- and bi-radionuclide). The mean NTCP over the three patients is also plotted for each implant. Zaider et al. radiobiological parameters are used. Error bars were calculated using (17).

The NTCPs are larger than the one published by Zaider et al. who obtained NTCPs of about 16% for 143 Gy I-125 LDR [1]. However, a more recent study by Zelefsky et al. reports 19% of patient experiencing late Grade 2 urinary symptoms [28].

All the NTCP results fall in the 19 to 23% range. The comparison between these data is therefore very difficult as these results are affected by large uncertainties. Their origin is discussed in the next section.

3.4. Uncertainties

We have calculated the contribution of each parameters uncertainty to the total and NTCP uncertainties. The values are tabulated in Tables 5 and 6, respectively.

Table 5: DU20-I uncertainty produced by each parameter and total DU20-I uncertainty for each implant using Zaider et al. radiobiological parameters.

(Gy)
¹²⁵I
Pd-I
¹⁰³Pd
Pd-Cs (a)
Pd-Cs (b)
¹³¹Cs

0.0
0.1
0.8
1.6
1.3
1.6

0.0
0.3
1.6
3.1
2.6
3.3

0.0
0.1
0.8
1.6
1.3
1.7

0.0
0.2
1.1
2.2
1.8
2.3

0.0
0.4
0.0
0.3
0.3
0.0

0.0
0.4
0.0
0.3
0.3
0.0

(tot)

0.0
5.6
19.3
7.1
6.4
0.4

0.0
7.2
24.8
9.3
8.3
0.7

Total

0.0
6.3
20.9
10.6
9.3
3.7

0.0
8.1
27.5
14.9
13.0
6.3

Table 6: NTCP uncertainty produced by each parameter and total NTCP uncertainty for each implant using Zaider et al. radiobiological parameters.

(%)
¹²⁵I
Pd-I
¹⁰³Pd
Pd-Cs (a)
Pd-Cs (b)
¹³¹Cs

8.6%
8.7%
8.4%
8.4%
7.8%
7.6%

5.6%
5.8%
5.2%
5.2%
4.4%
4.2%

0.0%
0.7%
2.3%
1.2%
1.0%
0.4%

0.0%
0.9%
3.0%
1.6%
1.3%
0.6%

Total

14.1%
15.3%
15.9%
14.7%
13.2%
12.2%

14.1%
15.5%
16.7%
15.2%
13.6%
12.4%

Firstly, the results show that the uncertainty is mainly due to the uncertainty on the value. The uncertainty is the greatest for ¹⁰³Pd and the least for ¹³¹Cs apart from ¹²⁵I which is not affected by the conversion and is not subjected to model-related errors. Sublethal damage repair rate, , as well as ratio have little effect on the value.

Secondly, the urethral NTCP uncertainties are mainly related to the parameter-fitting uncertainties. The error produced by has a secondary importance. The total relative uncertainty is therefore almost constant among the different modalities, ranging from 12 to 17%.

However, these figures also depend on the radiobiological data chosen for the calculations. If the Dale and Jones radiobiological parameters (and ) are used, the uncertainty generated by the sublethal damage repair rate becomes the most important one. The and urethral NTCP obtained with Zaider et al. and Dale and Jones radiobiological parameters are compared in Figures 4 and 5, respectively.

Figure 4: Comparison of for two different values of and  Gy and  h⁻¹ (Zaider et al. [1]) and  Gy and  h⁻¹ (Dale and Jones [2]). Error bars were calculated using (15).

Figure 5: Comparison of urethral NTCP for two different values of and  Gy and  h⁻¹ (Zaider et al. [1]) and  Gy and  h⁻¹ (Dale and Jones [2]). Error bars were calculated using (17).

¹³¹Cs is the most affected radionuclide. RBE uncertainty remains almost unchanged. Uncertainties produced by contributions are still negligible compared to other parameter’s uncertainties. These errors will have repercussions on NTCP uncertainties with the largest effect for ¹³¹Cs. The contribution of the total on the NTCP total uncertainty is now of the same order of magnitude as the one produced by the parameter-fitting error.

4. Discussion and Conclusions

The variation in the due to I-125 conversion is influenced by two major contributions. First, the linear term in the BED expression depends only on the value. Second, the quadratic term is more complex and results from the influence of three separate parameters: and which are tissue dependent and the radionuclide half-life.

Firstly, the ¹²⁵I will not be affected by the I-125 conversion as only one radionuclide is involved. Secondly, as one chooses the same value for ¹³¹Cs and ¹²⁵I, the linear term of the ¹³¹Cs BED is not affected by the ¹²⁵I conversion. Since the value of chosen by Zaider et al. is large, the quadratic contribution to the BED is small and influences the ¹³¹Cs BED by only 2%. If future experiments provide a value for ¹³¹Cs that is different from that of ¹²⁵I, larger variation in I-125 equivalent may be observed. Thirdly, ¹⁰³Pd has a larger value than ¹²⁵I. Its linear contribution to BED will therefore also contribute to a modified value of the . This contribution of the quadratic term accounts for only 1% of the total BED. The difference in value between ¹⁰³Pd and I-125 is therefore the main cause of the large increase in after I-125 conversion. Finally, the changes in the Pd-I and Pd-Cs due to I-125 conversion are related to both the linear and the quadratic contribution of each radionuclide to BED. Pd-Cs is the most affected mixture as both radionuclides differ from the one used for modeling (I-125).

The very similar values obtained for the NTCP and their large uncertainties makes it difficult, if not impossible, to conclude definitively whether, for equal tumour effect, bi-radionuclide brachytherapy would reduce the urethral complication probability relative to mono-radionuclide brachytherapy. Planning more patients would not improve the situation as these uncertainties are mainly due to the fitting parameter uncertainties of the empirical model of Zaider et al.

The large value assumed by Zaider et al. for the ratio reduces the influence of the quadratic term on the total BED considerably. Therefore the sublethal damage repair on the BED will show its effect only for radionuclides with short half-lives (like ¹³¹Cs). This is also proven by the low uncertainty associated with and to the ratio. The uncertainty related to the also indicates that our approximation will not significantly affect the final results.

However, the Zaider et al. radiobiological parameters are not the only ones in common use. The results obtained with the radiobiological parameters from Dale and Jones (Figures 4 and 5) clearly show that such parameters could also have a large impact on the final NTCP value. Moreover, the lower value of ratio increases the quadratic contribution of the BED, leading to larger uncertainties associated with and .

Finally, complication probabilities are certainly dependent on the way the treatment planning is performed. A logistic regression among patients treated by the each institution could also provide different NTCP parameter values. However, this would not affect the results.

We can conclude that the use of the urethra NTCP model for biologically based treatment planning in permanent seed prostate brachytherapy requires either better fitting parameters (with less associated errors) or a different NTCP model (different morbidity indicator than e.g.).

Appendix

A. Inclusion of RBE in TPS Dose Distribution

In order to include the RBE effect in the TPS dose distribution, the radial dose function and the 1D or 2D anisotropy functions have to be modified. This modification is trivial in the case of mono-radionuclide brachytherapy: the radial dose function is multiplied by the radionuclide value and the anisotropy functions are not altered. In the case of bi-radionuclide implants, the value will be different from one radionuclide to another. Therefore,

It can be shown that if (A.1) replaces (11), (14) and (16) in (4) or (5) of Nuttens and Lucas [23], this would give which means that the TPS dose distribution includes well the RBE effect.

B. Partial Derivatives of

is expressed as a function of , , and ((6) to (10)). Its derivative with respect to parameter can also be expressed as a function of , , and and their derivative with respect to parameter noted , and The expressions of , , and are now given for each parameter . , , , , ,

These expressions apply only for the bi-radionuclide implant for which

C. Partial Derivatives of NTCP

The NTCP expression is given by (1) and its partial derivatives with respect to each parameter are

Acknowledgments

V. E. Nuttens was supported by a Grant from the Belgian “Fonds pour la formation à la Recherche dans l’Industrie et dans l’Agriculture” (FRIA). Special thanks are due to H. C. Winter, C. D. Lee (Clatterbridge Center for Oncology) for helpful discussions.

References

M. Zaider, M. J. Zelefsky, G. N. Cohen, C. S. Chui, E. D. Yorke, L. Ben-Porat, and L. Happersett, “Methodology for biologically-based treatment planning for combined low-dose-rate (permanent implant) and high-dose-rate (fractionated) treatment of prostate cancer,” International Journal of Radiation Oncology Biology Physics, vol. 61, no. 3, pp. 702–713, 2005.

R. G. Dale and B. Jones, Radiobiological Modelling in Radiation Oncology, British Institute of Radiology, London, UK, 2007.

A. Niemierko and M. Goitein, “Implementation of a model for estimating tumor control probability for an inhomogeneously irradiated tumor,” Radiotherapy and Oncology, vol. 29, no. 2, pp. 140–147, 1993.

S. Webb and A. E. Nahum, “A model for calculating tumour control probability in radiotherapy including the effects of inhomogeneous distributions of dose and clonogenic cell density,” Physics in Medicine and Biology, vol. 38, no. 6, pp. 653–666, 1993.

B. Sanchez-Nieto and A. E. Nahum, “Bioplan: software for the biological evaluation of radiotherapy treatment plans,” Medical Dosimetry, vol. 25, no. 2, pp. 71–76, 2000.

A. Niemierko, M. Urie, and M. Goitein, “Optimization of 3D radiation therapy with both physical and biological end points and constraints,” International Journal of Radiation Oncology Biology Physics, vol. 23, no. 1, pp. 99–108, 1992.

Y. Kim and W. A. Tomé, “Risk-adaptive optimization: selective boosting of high-risk tumor subvolumes,” International Journal of Radiation Oncology Biology Physics, vol. 66, no. 5, pp. 1528–1542, 2006.

Z. Malik, C. Eswar, J. Dobson, J. Fenwick, and A. Nahum, “Iso-NTCP customization of the prescription dose in lung-tumour radiotherapy,” Radiotherapy and Oncology, vol. 84, pp. S278–S279, 2007.

A. Haworth, M. Ebert, D. Waterhouse, D. Joseph, and G. Duchesne, “Assessment of i-125 prostate implants by tumor bioeffect,” International Journal of Radiation Oncology Biology Physics, vol. 59, no. 5, pp. 1405–1413, 2004.

A. Haworth, M. Ebert, D. Waterhouse, D. Joseph, and G. Duchesne, “Prostate implant evaluation using tumour control probability—the effect of input parameters,” Physics in Medicine and Biology, vol. 49, no. 16, pp. 3649–3664, 2004.

W. D. D'Souza, H. D. Thames, and D. A. Kuban, “Dose-volume conundrum for response of prostate cancer to brachytherapy: summary dosimetric measures and their relationship to tumor control probability,” International Journal of Radiation Oncology Biology Physics, vol. 58, no. 5, pp. 1540–1548, 2004.

M. Zaider, “Permanent-implant brachytherapy in prostate cancer,” in New Technologies in Radiation Oncology, pp. 379–388, Springer, Berlin, Germany, 2006.

S. Mallick, R. Azzouzi, L. Cormier, D. Peiffert, and PH. Mangin, “Urinary morbidity after I brachytherapy of the prostate,” BJU International, vol. 92, no. 6, pp. 555–558, 2003.

Z. A. Allen, G. S. Merrick, W. M. Butler, K. E. Wallner, B. Kurko, R. L. Anderson, B. C. Murray, and R. W. Galbreath, “Detailed urethral dosimetry in the evaluation of prostate brachytherapy-related urinary morbidity,” International Journal of Radiation Oncology Biology Physics, vol. 62, no. 4, pp. 981–987, 2005.

M. Neill, G. Studer, L. Le, M. McLean, I. Yeung, G. Pond, and J. M. Crook, “The nature and extent of urinary morbidity in relation to prostate brachytherapy urethral dosimetry,” Brachytherapy, vol. 6, no. 3, pp. 173–179, 2007.

J. Crook, N. Fleshner, C. Roberts, and G. Pond, “Long-term urinary sequelae following iodine prostate brachytherapy,” Journal of Urology, vol. 179, no. 1, pp. 141–146, 2008.

T. Ohashi, A. Yorozu, K. Toya, S. Saito, and T. Momma, “Acute urinary morbidity following I-125 prostate brachytherapy,” International Journal of Clinical Oncology, vol. 10, no. 4, pp. 262–268, 2005.

S. A. Shah, R. R. Cima, E. Benoit, E. L. Breen, and R. Bleday, “Rectal complications after prostate brachytherapy,” Diseases of the Colon and Rectum, vol. 47, no. 9, pp. 1487–1492, 2004.

J. N. Shah and R. D. Ennis, “Rectal toxicity profile after transperineal interstitial permanent prostate brachytherapy: use of a comprehensive toxicity scoring system and identification of rectal dosimetric toxicity predictors,” International Journal of Radiation Oncology Biology Physics, vol. 64, no. 3, pp. 817–824, 2006.

F. M. Waterman and A. P. Dicker, “Probability of late rectal morbidity in I prostate brachytherapy,” International Journal of Radiation Oncology Biology Physics, vol. 55, no. 2, pp. 342–353, 2003.

T. Ohashi, A. Yorozu, K. Toya, S. Saito, T. Momma, H. Nagata, and M. Kosugi, “Rectal morbidity following I-125 prostate brachytherapy in relation to dosimetry,” Japanese Journal of Clinical Oncology, vol. 37, no. 2, pp. 121–126, 2007.

B. R. Prestidge, W. S. Bice, I. Jurkovic, E. Walker, S. Marianne, and A. Sadeghi, “Cesium-131 permanent prostate brachytherapy: an initial report,” International Journal of Radiation Oncology, Biology, Physics, vol. 63, pp. S336–S337, 2005.

V. E. Nuttens and S. Lucas, “AAPM TG-43U1 formalism adaptation and Monte Carlo dosimetry simulations of multiple-radionuclide brachytherapy sources,” Medical Physics, vol. 33, no. 4, pp. 1101–1107, 2006.

V. E. Nuttens and S. Lucas, “Determination of the prescription dose for bi-radionuclide permanent prostate brachytherapy,” Medical Physics, vol. 35, no. 12, pp. 5451–5462, 2008.

R. G. Dale and B. Jones, “The clinical radiobiology of brachytherapy,” British Journal of Radiology, vol. 71, pp. 465–483, 1998.

J. Z. Wang, N. A. Mayr, S. Nag, J. Montebello, N. Gupta, N. Samsami, and C. Kanellitsas, “Effect of edema, relative biological effectiveness, and dose heterogeneity on prostate brachytherapy,” Medical Physics, vol. 33, no. 4, pp. 1025–1032, 2006.

V. Antipas, R. G. Dale, and I. P. Coles, “A theoretical investigation into the role of tumour radiosensitivity, clonogen repopulation, tumour shrinkage and radionuclide RBE in permanent brachytherapy implants of I and Pd,” Physics in Medicine and Biology, vol. 46, no. 10, pp. 2557–2569, 2001.

M. J. Zelefsky, Y. Yamada, G. N. Cohen, A. Shippy, H. Chan, D. Fridman, and M. Zaider, “Five-year outcome of intraoperative conformal permanent I-125 interstitial implantation for patients with clinically localized prostate cancer,” International Journal of Radiation Oncology Biology Physics, vol. 67, no. 1, pp. 65–70, 2007.