Radiobiologically based treatment plan evaluation for prostate seed implants.

PURPOSE: Accurate prostate low dose-rate brachytherapy treatment plan evaluation is important for future care decisions. Presently, an evaluation is based on dosimetric quantifiers for the tumor and organs at risk. However, these do not account for effects of varying dose-rate, tumor repopulation and other biological effects. In this work, incorporation of the biological response is used to obtain more clinically relevant treatment plan evaluation.
MATERIAL AND METHODS: Eleven patients were evaluated. Each patient received a 145 Gy implant. Iodine-125 seeds were used and the treatment plans were created on the Prowess system. Based on CT images the post-implant plan was created. In the post-plan, the tumor, urethra, bladder and rectum were contoured. The biologically effective dose was used to determine the tumor control probability and the normal tissue complication probabilities for the urethra, bladder, rectum and surrounding tissue.
RESULTS: The average tumor control probability and complication probabilities for the urethra, bladder, rectum and surrounding tissue were 99%, 29%, 0%, 12% and 6%, respectively. These measures provide a simpler means for evaluation and since they include radiobiological factors, they provide more reliable estimation of the treatment outcome.
CONCLUSIONS: The goal of this work was to create more clinically relevant prostate seed-implant evaluation by incorporating radiobiological measures. This resulted in a simpler descriptor of treatment plan quality and was consistent with patient outcomes.

Purpose

Prostate cancer is one of the most common cancers in men, and is responsible for 32 000 deaths in the United States in 2010 [1]. Low dose-rate (LDR) brachytherapy is a common radiation therapy treatment used in the management of prostate cancer. The two most common radioactive seeds used in prostate LDR brachytherapy are 125I and 103Pd. 125I has a half-life of 59.4 days and emits a 28 keV gamma ray. Due to its relatively long half-life and penetrating gamma ray, 125I is often preferred for slow growing tumors [2]. The recommended prescription dose for LDR prostate monotherapy is 145 Gy using 125I seeds [3]. Post-implant treatment evaluation is a necessary step in LDR prostate brachytherapy. Due to the effects of seed migration and tissue edema, the final position of the seeds may be significantly different from their planned position. Since the seeds emits low energy short range radiation any disturbance in the seed distribution may result in unfavorable hot or cold spots. Thus, accurate post-implant evaluation is important to determine if the implant is satisfactory or if further treatment is needed to achieve the desired outcome.

Currently, the treatment plan evaluation is based on isodose distributions, dose-volume histograms (DVH) as well as mean, maximum and minimum doses in the target and organs at risk (OAR). Most of the evaluation tools currently used in clinical practice are dose-based. These do not account for the radiobiological characteristics of the tumor or OAR that are important in the biological outcome of the therapy. Under some circumstances, plans with similar dose distributions have been shown to have different estimated radiobiological outcomes [4]. Inclusion of radiobiological factors such as tissue radiosensitivity, tissue seriality, treatment time and non-linear response of tissues to radiation may improve the estimation of tumor control and normal tissue complications rates.

In this work, a new 3D treatment plan evaluation tool, based on the radiobiological response of each voxel is used. The goal of this new evaluation method is to simplify post implant evaluation while providing more clinically relevant information. Here, the probability of injury to the tumor and OAR is given as a simple probability. Using spatial information is preferable to the more common DVH-based radiobiological analysis. Using spatial information, areas of inadequate treatment can be easily visualized and used in creating additional treatment plans that target these areas.

Material and methods

This study is a retrospective analysis of eleven patients receiving prostate seed implant. Pre-implant prostate volumes were determined using a B-K Medical Leopard 2001 (Mileparken, Denmark) transrectal ultrasound (TRUS) machine with the patients in the lithotomy position. TRUS images were acquired during a continuous sagittal sweep. From the volumetric data acquisition, axial images were reconstructed with 2.5 mm spacing.

On the axial TRUS images, the prostate and urethra were contoured by a physician. The planning target volume (PTV) was defined to be the prostate volume plus 0.5 cm margin. The treatment plan was produced using the Prowess Panther 3D Brachy Pro (Prowess, Concord, California, USA). The treatment planning system uses the TRUS images to map out the placement of each seed with respect to a template. The template is simply a metal plate with holes drilled in a regular grid pattern. The template is then used in the operating room to guide each needle to the planned location. The seeds were implanted with precision under TRUS guidance to ensure their placement in the proper locations. All the patients were treated with 125I seeds (BARD, BrachySource model, Covington, GA, USA).

One month after the implant, the patient received a CT scan. Patients were scanned in the supine position and transverse slices were acquired with 2.5 mm slice thickness. On these CT images the prostate, urethra, bladder and rectum were contoured by a physician. Additionally, the location of each implanted seed was registered using these CT mages. The contoured CT images and seed locations were then used by the in-house software to perform the radiobiological evaluation.

The recommendation given in TG-64 is to perform a post-operative scan 28 days after the implant [5]. The American Brachytherapy Society considers the ideal scan time to be within 4 weeks after the implant [6]. The half-life of prostate edema has been reported to range from 4-25 days, with a mean of 9.3 days [7, 8]. Scanning several weeks after the procedure allows adequate time for most prostate edema to resolve [6]. This is important for plan evaluation. Due to the use of 125I seeds, most of the dose will be delivered after the prostate edema has resolved.

Using the seed positions from the treatment planning system the physical dose distribution was calculated using the TG-43 protocol [9]. Based on the physical dose and type of tissue present in each voxel, the biologically effective dose (BED) and the TCP were determined [4]. The value of the TCP is in part based on calculation of the biologically effective dose (BED), which for the tumor and normal tissues is calculated using equation (1a) and (1b) respectively [5, 10, 11].

where

In equation (1), R0 is the initial dose rate and λ is the decay constant (for 125I, λ = 0.01166 day–1). The sublethal damage repair constant (µ) was calculated using the following lowing expression:

This factor accounts for the decrease in cell kill as the cell repairs damage. Here, a general repair half-life of 15 minutes was assumed for both tumor and normal tissues, making µ = 2.8 hour-1 [12, 13].

The tumor repopulation factor (K) accounts for the growth of new tumor cells during treatment and is calculated as follows:

A potential doubling time (T) of 42 days was used in this analysis resulting in a repopulation factor K of 0.11 Gy × day-1 [12, 14].

Equation (4) was used for the calculation of the effective dose (Deff). The effective treatment time (T) was determined from equation (5). The endpoint for brachytherapy has been defined as the point where the rate of cell kill is equal to the tumor repopulation factor [15]. For normal tissues it is assumed that T = ∞, hence the effective dose is taken to be equal to the total physical dose accumulated over the lifetime of the seeds.

The relative biological effectiveness (RBE) for 125I has been reported to be between 1.4-1.5 [16, 17]. For this study, a middle value of RBE = 1.45 was used. The specific radiobiological parameters α/β, D50 and γ used for each tissue are shown in Table 1. D50 is the dose which gives a 50% response and γ is the maximum normalized dose-response gradient. The D50 and γ parameters are derived from clinical materials and describe the shape of the dose-response curve [18]. These parameters and the radiobiological model used are similar to those recommended in TG-137 [19]. Biological parameters are the subject of some uncertainty due to intra-patient radiosensitivity variations [20-23]. The response probability (P) for each voxel was then determined using equation (6) [24-26].

Table 1

Dose-response parameters for each tissue evaluated

Tissue Type	γ	α/β	D50	s	End point
PTV	4	2	50	1.0	Control
Bladder	3	3	80	0.3	Symptomatic Contracture
Urethra	3	3	230	0.01	RTOG Grade 2
Urethra	3	3	275	0.01	RTOG Grade 3
Rectum	2.2	3	80	0.7	Proctitis, Fistula, Stenosis
Other	4	3	80	0.5	Necrosis

The overall response probability for the tumor and normal tissues is calculated using equations (7a) and (7b), respectively.

Where N is the total number of voxels in the organ, s is the tissue-specific relative seriality parameter and ∆υi is the fractional subvolume of the organ irradiated. The overall probability of injury to the involved normal tissues (PI) and the complication-free tumor control probability (P+) for the treatment were calculated using equations (8) and (9), respectively. For the calculation of P +, both RTOG grade 2 and grade 3 complications to the urethra were considered. The biologically effective uniform dose () which is the uniform dose that causes the same tumor control as actual dose distribution for a given treatment was calculated from equation (10).

The biologically effective uniform dose () calculates the uniform dose that would provide the same clinical outcome as the inhomogeneous dose distribution. It is a function of physical dose and tissue specific radiobiological parameters. The general expression of () is derived numerically from the first part of the following equation, where for a tissue of uniform radiosensitivity, () is given from the analytical formula of the second part of equation (11).

where denotes the 3-dimensional dose distribution delivered to the tissue and P() is the response probability of the tissue. The second part of the equation has been derived using the Poisson model [4].

Results

For all of the patients evaluated in this study, the calculated tumor control probability was very high, 99.6% on average. This result agrees with the dosimetric evaluation. Using the criteria prostate D 90 > 145 Gy for a successful implant, all the implants were considered to deliver an adequate dose to the tumor. The results of the radiobiological evaluation are presented in Table 2. The mean dose to each organ for each patient is given in Table 3. A dosimetric evaluation for each patient is shown in Table 4.

Table 2

Results of radiobiological evaluation. The P+ G2 and P+ G3 columns give the P+ calculation for urethra RTOG grade 2 and 3 injury, respectively

Patient	TCP (%)	NTCP Bladder (%)	NTCP Rectum (%)	NTCP Urethra G2 (%)	NTCP Other (%)	P+ G2 (%)	NTCP Urethra G3 (%)	P+ G3 (%)
1	99.8	0.2	2.2	0.1	7.3	90.3	0.0	90.3
2	100.0	0.0	6.9	20.1	11.6	65.8	0.4	82.0
3	99.5	0.0	8.5	2.3	8.1	81.7	0.0	83.6
4	99.1	0.0	21.4	2.2	6.4	71.1	0.0	72.8
5	99.9	0.0	7.7	30.3	7.8	59.3	1.4	83.9
6	99.8	0.0	16.3	27.7	3.9	58.0	1.0	79.5
7	99.9	0.0	6.0	37.4	1.5	57.9	3.1	89.7
8	99.8	0.0	6.5	66.6	4.3	29.7	17.4	73.7
9	100.0	0.0	13.7	19.9	3.4	66.9	0.0	83.4
10	100.0	0.0	13.3	28.3	8.5	56.9	1.0	78.5
11	99.9	0.0	25.5	85.0	6.3	10.4	48.2	36.1
Avg ± SD	99.8 ± 0.3	0.0 ± 0.0	11.6 ± 7.1	29.1 ± 26.6	6.3 ± 2.8	58.9 ± 22.3	5.0 ± 10.0	77.6 ± 14.9

Table 3

Mean dose to each organ for each patient

Patient	Prostate (Gy)	Urethra (Gy)	Bladder (Gy)	Rectum (Gy)
1	309	182	25.4	29.0
2	343	241	25.2	42.1
3	342	213	39.4	34.0
4	319	220	14.9	51.6
5	335	259	27.4	36.7
6	300	247	18.6	49.6
7	297	261	15.4	33.3
8	348	287	15.5	34.0
9	321	242	19.6	29.9
10	353	225	17.2	32.7
11	357	320	18.0	50.9
Avg ± SD	329 ± 21	245 ± 37	21.5 ± 7.4	38.5 ± 8.5

Table 4

Dosimetric evaluation based on the D90 and V150 from the prostate and maximum rectal dose. The bracketed organ indicated the target for each dose quantifier

Patient	PTV D90 (Gy)	PTV V150 (%)	Rectum Dmax (Gy)	D90 > 145 Gy [tumor]	D90 < 180 Gy [urethra]	V150 < 60% [urethra]	RDmax < 178 Gy [rectum]
1	164	63.9	135.5	yes	yes	no	yes
2	184	78.4	151	yes	no	no	yes
3	161	71.8	183	yes	yes	no	no
4	160	71.4	231	yes	yes	no	no
5	184	78.4	165	yes	no	no	yes
6	169	72.2	174	yes	yes	no	yes
7	154	61.9	155	yes	yes	no	yes
8	171	77.6	161	yes	yes	no	yes
9	180	76	304	yes	no	no	no
10	188	78.4	213	yes	no	no	no
11	190	82.3	210	yes	no	no	no

Considering the rectum, Merrick et al. recommend that the maximum dose to the anterior mucosal wall should be less than 120% of the reference dose [27]. For an implant reference dose of 145 Gy, this 120% dose is 178 Gy. Comparing the radiobiological results with the dosimetric parameters, patients that had higher rectum complication probabilities tended to have larger maximum rectal doses. However, the patient with the highest maximum rectal dose did not have the highest rectum complication probability. This indicates the weakness of simple dosimetric quantifiers and the advantage of radiobiological evaluation which considers biological factors in relation with the dose distribution. Making an allowance for the bladder, the average complication probability was 0.03%. This is consistent with our clinical experience that the bladder is typically located far enough from the prostate, so that the dose delivered to it is low.

Considering the urethra, the average probability for RTOG grade 2 and grade 3 complications was 38.95% and 9.81%, respectively [28]. There are several dosimetric quantifiers reported in the literature to predict urethral complications, two of which are presented in Table 4 [29-31]. Wallner et al. have shown that maximum urethral doses of 447 Gy and 592 Gy are related to RTOG grade 0-1 and grade 2-3 complications, respectively [32, 33]. However, it is seen in Table 4 that for our group of patients these dose quantifiers provide conflicting results.

Since the radiobiological evaluation in this system is calculated using spatial dose information, is it possible to visualize the response probabilities as one might do with physical dose. More simplistic radiobiological analyses that are based on DVH information do not have this capability. Figures 1 and 2 show that the response probability distribution is different from the dose distribution. This illustrates the advantage of incorporating radiobiological considerations to the implant evaluation, since it provides better estimation of the clinical outcome compared to physical dose alone.

Fig. 1

Images of the physical dose, BED and response probability for patient 3. Organ contours are overlaid on each image

Fig. 2

Images of the physical dose, BED and response probability for patient 11. Organ contours are overlaid on each image

The upper images in Figs. 3 and 4 show that patient 3 represents more conformal treatment, while patient 11 represents a less conformal treatment. The DVH of patient 11 shows that the dose distributions to the rectum and urethra are significantly higher than those of patient 3, while the dose distributions to the tumor are comparable. The initial seed activities for patients 3 and 11 were 0.38 and 0.37 mCi, respectively. The middle plot of Fig. 3 indicates that the P I begins to rise just above this activity. The middle plot of Fig. 4 shows that even for low initial seed activities the PB is saturated and P I is significant. This would seem to indicate that too many seeds were used in this implant. The lower plot of Fig. 3 shows that P + is high for actual implanted of 138 Gy. However, this plot also indicates that the most ideal balance of PB and PI is obtained at a equal to 120 Gy. The lower plot of Fig. 4 shows that the which would be delivered with even low strength seeds is significantly higher than the ideal level. Figures 5 and 6 include the response of the urethra to RTOG grade 3 injury. The higher dose limit for grade 3 injury has the effect of pushing the P I curve to higher doses, thus broadening the P + curve.

Fig. 3

The plots shown above are for patient 3, calculated with urethra RTOG grade 2 injury. The upper image is the DVH. The middle image shows the response probabilities of each organ for an identical implant, but with different initial seed activity. The lower figure shows the P +, P B and P I vs. curves for identical implants with different initial seed activities. The vertical line indicates the clinically prescribed seed activity (middle) or (lower)

Fig. 4

The plots presented above are for patient 11, calculated with urethra RTOG grade 2 injury. The upper image is the DVH. The middle image shows the response probabilities of each organ for an identical implant but with different initial seed activity. The lower figure shows the P+, PB and PI vs. = curves for identical implants with different initial seed activities. The vertical line indicates the clinically prescribed seed activity (middle) or = (lower)

Fig. 5

The plots shown above are for patient 3, calculated with urethra RTOG grade 3 injury. The upper image shows the spatial response probability distribution. The middle image shows the response probabilities of each organ for an identical implant but with different initial seed activity. The lower figure shows the P+, PB and PI vs. = curves for identical implants with different initial seed activities. The vertical line indicates the clinically prescribed seed activity (middle) or = (lower)

Fig. 6

The plots presented above are for patient 11, calculated with urethra RTOG grade 3 injury. The upper image shows the spatial response probability distribution. The middle image shows the response probabilities of each organ for an identical implant but with different initial seed activity. The lower figure shows the P+, PB and PI vs. curves for identical implants with different initial seed activities. The vertical line indicates the clinically prescribed seed activity (middle) and (lower)

Discussion

The evaluation of prostate implants is commonly done using dose quantifiers such as D 100, D 90 and V 150 to compare the implant to known TCP and NTCP. However, such simple quantifiers demonstrate only a small piece of the picture. Prostate implants with similar D90 can have significantly different D 50 and D max. To glean a better idea of the TCP and NTCP one must consider many dose quantifiers. At some point, the increased information that should provide the basis for an excellent evaluation becomes overwhelming. With the proposed methodology based on the spatial distribution of the probability of response, the effectiveness of a given dose distribution is much more clear, easy to assess and clinically relevant.

Figures 2 and 3 presents that the response probability distribution eliminates much of the visual noise due to low doses that may not have a clinical effect. These images may be used to indicate regions that might need further irradiation. Radiobiological evaluation of a treatment provides additional information about the fitness of the implant and a closer association to the clinical outcome. This evaluation takes into account the dose-response characteristics of the tumor and normal tissues that are important in the treatment. The concurrent display of both the radiobiological evaluation and the dosimetric data shows their complementary relation in analyzing the treatment.

The TCP for all of the patients was very high. This result is consistent with the dosimetric data, which shows that all the patients had prostate D 90 well above the acceptable criteria of 145 Gy. Sparing the urethra represents a technical difficulty since the urethra passes directly through the prostate. Thus, one must balance the need for adequate prostate dose coverage while minimizing the urethral dose. Even using dose limiting techniques, the risks of urethra complications was high for several patients. Radiobiological parameters for the urethra are sparsely reported in the literature. For this study, the parameters were chosen such that the calculated response probability was consistent with dose-response information reported by others. Takeda et al. has reported that for a patient cohort with average urethral D 90 = 150 Gy, 90% experienced RTOG grade 2 symptoms and 6.5% experienced grade 3 symptoms [34]. Additionally, Thomas et al. has reported that for a patient cohort with average urethral D 90 = 105 Gy, 36% and 7% experience grade 2 and 3 symptoms, respectively [35]. Thus, the urethra complication rates predicted for the patients in this study are comparable with these results.

Comparing the radiobiological results with the dosimetric parameters, patients that had higher rectum complication probabilities tended to have larger maximum rectal doses. However, the patient with the highest maximum rectal dose did not have the highest rectum complication probability. Authors have reported that rectal complications are most frequent eight months post-implant, with 8-18% of patients reporting complications [27, 36]. The predicted rectal complication rate calculated for our patient cohort was similar, at 11%. Combining all the tissue-response information, the success of the treatment may be given as a single value, the complication-free tumor control probability (P +). The average P + for the patient cohort was 58.9%. However, the value varied amongst the patients, from 10.4% to 90.3%.

Two specific patients were selected in order to present the results of a conformal implant and a less conformal implant. Comparing Figs. 1 and 2, the physical dose distribution for patient 3 seems uniform compared to patient 11 which is characterized by higher inhomogeneities in the dose distribution. The response probability distribution shows that the tumor response for both patients is high. It also shows that while the physical dose distributions look quite different, the tumor responses are similar. However, the response probability distribution for patient 11 shows that the urethra response is high and that a greater area of normal tissue surrounding the tumor is responding compared with patient 3. Figure 3 shows that for patient 3 the implant was ideal with respect to the initial seed strength and overall to the tumor. Figure 4 indicates that an ideal treatment would require much lower activity seeds for patient 11. This would seem to indicate that too many seeds were implanted in this patient for the given source strength. Figures 5 and 6 represent the evaluation where RTOG grade 3 injury is considered as the urethra endpoint of interest. Since the urethra is the most significant component of P I, the higher dose-response limit pushes the P I curve to higher doses. This has the effect of broadening the P + curve toward higher doses.

The clinical outcomes for the patient cohort were consistent with our predicted outcomes. All patients had PSA scores and digital rectal exams indicating satisfactory improvement. Most patients with low normal tissue response probabilities, including patient 3, reported minimal symptoms such as increased nocturia and slight rectal irritation. Patients with higher response probabilities reported more severe symptoms such as erectile dysfunction, bloody stools, increased nocturia, urinary hesitancy and pain with urination. Patient 11 reported the most severe symptoms of the patient cohort.

The main advantage of this evaluation method is that it is based on the response probability calculated from each voxel. This allows for response probability maps to be generated. The major limitation of this system stems from the use of specific radiobiological parameters, which are known to be the subject of various types of uncertainties. Thus, the probabilities calculated should only be considered as estimates of the response for a typical patient. However, this limitation is also shared by all the common radiobiological and dosimetric quantifiers, because the degree of radiation induced effect for a specified dose level is also variable within a patient population.

Conclusions

The goal of this study was to improve the evaluation of prostate seed implants by developing a radiobiological evaluation system based on dose-response probabilities. The system proposed here uses the spatial dose information from the treatment planning system combined with tissue-specific radiobiological factors to estimate tumor control probability and normal tissue complication probabilities for each implant. The evaluation of a small patient cohort with both the proposed methodology and common dosimetric measures shows that the proposed system yields results consistent with those of reported patient outcomes. In summary, this system uses a potentially more comprehensive radiobiological analysis to estimate patient outcomes and the results are reported in a simple and intuitive manner. The authors envision that such a system may be used along with traditional dose quantifiers to provide more inclusive and clinically relevant evaluation of prostate seed implants.