A Technique for Stereotactic Radiosurgery Treatment Planning with Helical Tomotherapy
Article Outline
- Abstract
- Introduction
- Materials and Methods
- Results
- Discussion
- Conclusions
- Acknowledgment
- References
- Copyright
Abstract
The purpose of this study was to develop an efficient and effective planning technique for stereotactic radiosurgery using helical tomotherapy. Planning CTs and contours of 20 patients, previously treated in our clinic for brain metastases with linac-based radiosurgery using circular collimators, were used to develop a robust TomoTherapy planning technique. Plan calculation times as well as delivery times were recorded for all patients to allow for an efficiency evaluation. In addition, conformation and homogeneity indices were calculated as metrics to compare plan quality with that which is achieved with conventional radiosurgery delivery systems. A robust and efficient planning technique was identified to produce plans of radiosurgical quality using the TomoTherapy treatment planning system. Dose calculation did not exceed a few hours and resulting delivery times were less than 1 hour, which allows the process to fit into a single day radiosurgery workflow. Plan conformity compared favorably with published results for gamma knife radiosurgery. In addition, plan homogeneity was similar to linac-based approaches. The TomoTherapy planning software can be used to create plans of acceptable quality for stereotactic radiosurgery in a time that is appropriate for a radiosurgery workflow that requires that planning and delivery occur within 1 treatment day.
Key Words: Stereotactic radiosurgery, Tomotherapy, Brain metastases, Treatment planning
Introduction
Tomotherapy has shown potential as a precision stereotactic radiosurgery (SRS) delivery system.1, 2, 3 Recent work carried out by our group has shown that the on-board megavoltage computed tomography (MVCT) can be used for stereotactic localization and the system is capable of submillimeter delivery accuracy if used in conjunction with a precise intracranial stereotactic patient positioning system.1 In addition, planning studies have shown that TomoTherapy produces conformal plans for small targets.4, 5 Because the delivery is nonisocentric, it could potentially provide an efficiency advantage in the case of multiple intracranial targets or large irregularly-shaped targets that would require multiple isocenters when using conventional SRS. Although TomoTherapy lends itself well to fractionated stereotactic procedures, such as stereotactic radiotherapy (SRT) and stereotactic body radiotherapy (SBRT), the unique challenges involved with single-fraction delivery and the use of an invasive fixation system are the focus of this particular work.
Stereotactic localization can be achieved with TomoTherapy using image guidance in combination with rigid head fixation. We therefore propose to use an invasive head frame for fixation of the patient, while on-board MVCT is used for stereotactic localization and an optical tracking system is used for shift verification. Moreover, to avoid difference in table sag between computed tomography (CT) and TomoTherapy, we propose that a thin-slice pretreatment MVCT is used in the SRS treatment planning process. We therefore have designed the following workflow: After head ring placement, the patient is set up on the TomoTherapy couch using a TomoTherapy-specific tabletop frame docking device (InterFix, Integra Radionics, Burlington, MA) for acquisition of the thin-slice reference/planning MVCT image dataset. An array of passive infrared markers is placed on the frame for optical tracking (Dynatrac, 3D Line Medical Systems, Milan, Italy). Camera and couch coordinates are recorded at the time of initial imaging to aid in repeat setup of the patient. After scanning, the thin-slice MVCT image is sent to third-party treatment planning software, where it can be fused with magnetic resonance imaging (MRI) for target delineation. A treatment plan is generated on the TomoTherapy Treatment Planning Station (TPS) using the MVCT for dose calculation. The patient is then repositioned on the couch, using the couch and camera coordinates recorded previously, and a pretreatment verification MVCT is performed. An MVCT-MVCT fusion of the planning and pretreatment verification images is completed to localize the target for delivery. If a shift is required, the shift is applied to the couch, whereas secondary verification of the magnitude and direction of the shift is performed via optical tracking.
From a treatment planning perspective, generating a deliverable radiosurgery plan with the time requirements introduced by this workflow is challenging. Conventional SRS planning can be performed quickly using straightforward planning approaches, in which known beam arrangements are used to create spherical and ellipsoidal shaped dose distributions. TomoTherapy involves a sophisticated inverse planning technique with a relatively lengthy dose calculation. Emphasis must therefore be placed on making good initial choices in planning parameters to avoid a time-consuming trial-and-error process. This work aims to identify an SRS planning technique for TomoTherapy that can be performed quickly yielding high-quality plans for the treatment of brain metastases, which compare well with other forms of radiosurgery.
Materials and Methods
Test patients
Twenty patients with brain metastases, previously treated with cone-based linear accelerator SRS, were used to test for the efficacy and robustness of the proposed TomoTherapy planning technique. Patient details are shown in Table 1. Patients selected had between 1−6 intracranial metastases, with target volumes ranging from 0.02−10.78 mL, with a median volume of 0.38 mL and an average volume of 1.06 mL. Patients' original planning CT and MR-defined target volumes were used for TomoTherapy comparison planning.
Table 1. Information regarding targets and prescription for test patients
| Patient # | Target # | Target Location | Dose (Gy) | Volume (mL) | Total Volume (mL) |
|---|---|---|---|---|---|
| 1 | 1 | Rt | 21 | 0.52 | 0.52 |
| 2 | 4 | Lt | 18 | 0.32 | 6.64 |
| Lt | 15 | 5.81 | |||
| Rt | 18 | 0.16 | |||
| Rt | 18 | 0.35 | |||
| 3 | 3 | Rt | 21 | 0.61 | 3.77 |
| Lt | 21 | 1.24 | |||
| Rt | 21 | 1.92 | |||
| 4 | 3 | Cerebellar | 21 | 0.45 | 1.86 |
| Lt | 21 | 0.34 | |||
| Lt | 21 | 1.07 | |||
| 5 | 3 | Lt | 18 | 0.82 | 2.77 |
| Rt | 14 | 0.65 | |||
| Rt | 18 | 1.30 | |||
| 6 | 6 | Lt | 20 | 0.93 | 3.23 |
| Rt | 20 | 0.60 | |||
| Rt | 20 | 1.45 | |||
| Lt | 20 | 0.08 | |||
| Lt | 20 | 0.12 | |||
| Rt | 20 | 0.04 | |||
| 7 | 3 | Lt | 15 | 2.48 | 5.78 |
| Rt | 18 | 0.04 | |||
| Rt | 15 | 3.26 | |||
| 8 | 4 | Lt | 15 | 0.07 | 0.59 |
| Lt | 15 | 0.40 | |||
| Lt | 17 | 0.04 | |||
| Rt | 21 | 0.07 | |||
| 9 | 2 | Rt | 20 | 0.84 | 7.17 |
| Lt | 20 | 6.33 | |||
| 10 | 6 | Rt | 21 | 1.79 | 2.93 |
| Lt | 21 | 0.43 | |||
| Rt | 21 | 0.16 | |||
| Rt | 21 | 0.03 | |||
| Rt | 21 | 0.07 | |||
| Lt | 17 | 0.46 | |||
| 11 | 3 | Rt | 21 | 0.15 | 0.32 |
| Lt | 15 | 0.10 | |||
| Rt | 21 | 0.08 | |||
| 12 | 4 | Paracentral | 15 | 10.79 | 11.26 |
| Rt | 21 | 0.12 | |||
| Rt | 21 | 0.21 | |||
| Lt. | 21 | 0.14 | |||
| 13 | 3 | Lt. | 18 | 1.92 | 2.17 |
| Rt | 23 | 0.10 | |||
| Rt | 23 | 0.15 | |||
| 14 | 4 | Lt | 21 | 0.06 | 0.22 |
| Lt | 21 | 0.04 | |||
| Rt | 18 | 0.10 | |||
| Lt | 21 | 0.02 | |||
| 15 | 3 | Rt | 21 | 0.14 | 0.24 |
| Lt | 24 | 0.05 | |||
| Lt | 24 | 0.05 | |||
| 16 | 2 | Rt | 18 | 2.61 | 3.02 |
| Lt | 21 | 0.41 | |||
| 17 | 2 | Rt | 21 | 1.62 | 3.54 |
| Lt | 21 | 1.92 | |||
| 18 | 1 | Lt | 15 | 5.13 | 5.13 |
| 19 | 2 | Lt | 24 | 1.03 | 1.23 |
| Rt | 24 | 0.20 | |||
| 20 | 1 | Medial | 21 | 1.19 | 1.19 |
Planning goals
The goal of TomoTherapy planning was to generate clinically acceptable plans within a reasonable planning time. A clinically acceptable plan is one in which dose conforms tightly to the target, and dose fall-off into extra-target normal tissues is as steep as possible given the delivery technique. In other forms of radiosurgery, doses are typically prescribed to 50–80% of target maximum dose depending on the SRS modality to ensure that the tumor boundary is on the steepest part of the dose gradient. The work by Lax6 shows that allowing for heterogeneous target coverage reduces the extra-target dose, which is one of the guiding principles of SRS and has been shown to hold true for TomoTherapy as well.7, 8 Hence, any TomoTherapy SRS planning technique should generate plans that are sufficiently inhomogeneous to allow for a steep dose gradient beyond the target despite the strong tendency of the TomoTherapy optimization algorithm to produce homogeneous dose distributions that are typically desirable for fractionated radiotherapy. The resulting planning techniques should be simple and result in treatment delivery times that are comparable to or shorter than those of conventional SRS techniques.
Treatment planning software
All plans were run using version 3.1.4 of the TomoTherapy Hi-Art treatment planning system (TPS). This software release has 2 primary advantages over previous versions. First, the gradient of the dose fall-off in the longitudinal direction is increased because the first and last delivered beamlets start closer to the tumor boundaries. Second, version 3.1.4 corrects an anomaly in volume element summation present in earlier versions, which impacts dose-volume histogram (DVH) generation for small (<2 mL) target volumes and regions at risk in steep dose gradient regions.9
Planning technique
ContouringBased on our own experience and previously published studies,8, 10 it was determined that additional nonanatomic planning structures would be needed to help create a dose distribution with the desired characteristics. Here the goal was to define useful planning structures that would improve the speed of plan optimization while keeping contouring time to a minimum. For the final technique, two different sets of added planning structures were identified: one for small targets (volume ≤2 mL) and another for large targets (volume >2 mL). In addition to the physician-defined target, the following contours were added:
For larger targets having a volume larger than 2 mL, the following structures may also be needed:
The planning volumes for the two cases have been illustrated in Figure 1.

Fig. 1.
Examples of added contours for treatment planning. Small target volumes require only a ring and a very small central subvolume (CSV), whereas larger tumors require 2 subvolumes (the CSV and a larger SV), a ring inside the periphery of the volume, and another ring external to the target.
Initial plan parameters
All plans were calculated using the “Fine” dose grid, resulting in a dose voxel size of 1.41 × 1.41 × 1.25 mm (where image pixels are directly used as dose voxels). A field width of 10 mm was used for treatment planning to maximize the rate of dose fall-off at the superior and inferior boundaries of the target.
The modulation factor (MF) determines the range of allowable intensity values for the leaves of the multileaf collimator (MLC) and is defined as the ratio of maximum to mean MLC leaf opening times for the leaves that are open. As shown in Figure 2, increasing modulation can be used to increase dose inhomogeneity within the target. Figure 3 shows the corresponding increase in the steepness of the peripheral dose gradient. However, increasing modulation also proportionally increases the gantry period (when it is above its minimum value) and therefore leads to increased delivery times. Plan quality and efficiency must be considered together to determine an ideal MF. A MF of 1.7 was chosen for all plans to obtain a balance between dose inhomogeneity and treatment time although a range between about 1.5 and 2.0 would be acceptable for most cases.

Fig. 2.
Dose inhomogeneity, as measured with a homogeneity index, as a function of modulation factor for a single metastasis.

Fig. 3.
Change in maximum dose to the target and the minimum dose at 3 mm from the target, with increasing modulation factor for one target. Results show that as the intra-target hot spot increases, peripheral fall off is steeper.
For helical delivery, pitch refers to the fraction of the field width the couch moves into the bore per rotation of the gantry. This parameter thus defines the degree of beam overlap between each rotation and determines the number of rotations over which a target voxel will be within the beam projections. For SRS planning, there are two main considerations in terms of pitch. First, when pitch is reduced, the gantry period is decreased because an increased number of rotations are available to deliver dose to a given voxel and each rotation delivers a smaller dose. Second, as the dose per fraction increases, the gantry period increases. Therefore in SRS, where the dose per fraction is high, it might not be possible to deliver the entire dose while keeping the gantry rotation period within the 60 second limit and multiple “passes” might be necessary to deliver the prescribed dose. If one is intent on keeping the delivery to a single pass, then the pitch should be as small as possible. However, from a workflow perspective a small pitch increases the beamlet computation time. Hence, it will be inefficient to choose a pitch that is too small. For a patient with four metastases, Figure 4 shows the beamlet computation time as a function of pitch for a clinically used 16-node cluster. The goal here was to determine the largest pitch that will allow for delivery of the entire dose in a single pass. Ideal pitches consistently produce a gantry period between 50 and 60 sec.

Fig. 4.
Beamlet computation time as a function of pitch for one patient with 4 metastases using a clinically configured planning station and cluster.
Table 2 shows the range of pitch values that allow for single-pass delivery as a function of the minimal peripheral dose and a MF of around 1.7. However, if the pitch results in a gantry period of slightly longer than 60 sec, a small decrease in MF can be used to increase the gantry speed as described by Woch et al.11
Table 2. Pitch values as a function of prescribed dose
| Maximum Prescribed Dose (Gy) | Pitch |
|---|---|
| 15 | 0.18 |
| 18 | 0.17 |
| 20 | 0.14 |
| 21 | 0.13–0.14 |
| 24 | 0.13 |
Plan objective and optimization
Through iteration, it was determined that the planning objectives shown in Table 3a and Table 3b reliably yield inhomogeneous dose distributions for small (<2 mL) and large (>2 mL) targets, respectively. The aim is to “weight” the minimum dose to the target and the maximum dose to the surrounding 3-mm ring structure higher than all other target goals so that the prescription isodose coincides with the target boundary. The maximum dose to the target is set to the maximum allowed (120 Gy) to ensure that it is removed from the optimization. The maximum dose in the target is then set by the CSV, which is set to 125% of the prescribed dose to mimic prescribing the 80% line in conventional SRS. For larger targets, instead of assigning optimization goals to the target structure itself, one uses a larger SV at the center of the target (SV) and the 2-mm outer shell as “Tumor” structures to better shape the dose, whereas the target structure is given a “Sensitive Structure” type for target DVH evaluation. Two additional elevated weights are given to the maximum dose to the SV and the uniform (max and min) dose to the 2-mm inner shell to ensure the shell is at the prescribed dose while the SV is at 125% of the prescribed dose. The CSV is still included and prescribed to 125% of the prescription dose to aid in plan scaling. The cutoff where the small-volume technique breaks down is somewhat dependent on target shape and location but any target greater than 2 mL may require the large volume technique to achieve the desired dose distribution.
Table 3a. Initial optimization constraints for TomoTherapy radiosurgery planning for small (<2 mL) targets
| Tumor Settings | |||||||
|---|---|---|---|---|---|---|---|
| Type | I | Max Dose | Max Dose Penalty | DVH % | DVH Dose | Min Dose | Min Dose Penalty |
| Target | 1 | 120 Gy⁎ | 1 | 100 | Prescription | Prescription | 100 |
| CSV† | 1 | Prescription/0.8 | 1 | 100 | Prescription/0.8 | Prescription/0.8 | 1 |
| Sensitive Structure Settings | ||||||
|---|---|---|---|---|---|---|
| Type | I | Max Dose | Max Dose Penalty | DVH % | DVH Dose | DVH Penalty |
| 3-mm ring | 1 | Prescription | 100 | 50 | Prescription/2 | 1 |
⁎Maximal allowed dose to be prescribed. |
†Indicates one of these structure types is used as the prescription. |
Table 3b. Initial optimization constraints for TomoTherapy radiosurgery planning for large (>2 mL) targets
| Tumor Settings | |||||||
|---|---|---|---|---|---|---|---|
| Type | I | Max Dose | Max Dose Penalty | DVH % | DVH Dose | Min Dose | Min Dose Penalty |
| CSV⁎ | 1 | Prescription/0.8 | 1 | 100 | Prescription/0.8 | Prescription/0.8 | 1 |
| 2-mm RS | 1 | Prescription | 1 | 100 | Prescription | Prescription | 100 |
| SV | 1 | Prescription/0.8 | 1 | 100 | Presciption/0.8 | 1 | Same as CSV |
| Sensitive Structure Settings | ||||||
|---|---|---|---|---|---|---|
| Type | I | Max Dose | Max Dose Penalty | DVH % | DVH Dose | DVH Penalty |
| Target | 1 | 120Gy† | 1 | 100 | Prescription | 100 |
| 3 mm ring | 1 | Prescription | 100 | 50 | Prescription/2 | 1 |
⁎Indicates that one of these structure types is used as the prescription. |
†Maximal allowed dose to be prescribed. |
Controlling the dose near the outer edge of the target with the 2-mm shell helps prevent this peripheral dose from becoming elevated, as shown in Figure 5 for a large target planned with both the small volume and large target techniques. This elevated dose is caused by simultaneously trying to keep the target dose above the prescription and trying to create a fast dose fall-off outside the target—a phenomenon known as “Gibbs ringing.”12 The frequency of this ringing is such that small targets do not suffer from it and the 2-mm shell is unnecessary.

Fig. 5.
Two dose distributions for a 10.5 mL target. Although both dose distributions are inhomogeneous, in (a) the highest dose is located peripherally, and in (b) it is located centrally. The result on the right is achieved by setting harder constraints on the 2-mm outer ring structure (2 mm RS), shown in purple, whose outer edge is formed by the periphery of the treatment target, and using the additional larger subvolume (SV), shown in red, at the center of the target to force a higher dose at the center.
This planning strategy works for reasons directly related to the TomoTherapy planning system's objective function. The TomoTherapy TPS uses an objective function that favors homogeneous dose in structures designated as targets. The optimizer objective function is given by13
(1)
.is the prescribed dose to structure x, always 0 for sensitive structures
.is the dose deposited in voxel i
is the importance value for structure x
.is the number of voxels in structure x
is target structure x's maximum dose penalty if voxel i is greater than x's maximum dose, otherwise 1
is target structure x's minimum dose penalty if voxel i is less than x's minimum dose, otherwise 1
is sensitive structure x's maximum dose penalty if voxel i is greater than x's maximum dose, otherwise 0
is sensitive structure x's DVH point penalty if voxel i violates x's DVH constraint, otherwise 0This objective function favors uniformity in targets because of the leftmost summation of the 2 summations. For a target, even if no constraint is violated, all of its voxels still contribute to the sum because in that case
, and their individual contribution to the overall sum is equal to their difference from the intended prescription dose, which is the same for all voxels in the target. For sensitive structures, which are governed by the rightmost summation, only voxels that violate a constraint make a nonzero contribution to the objective function.
To achieve nonuniform dose distributions, it is therefore necessary to counter the tendency of the optimizer to generate a homogeneous dose distribution within the target as well as to control how the nonuniform dose is distributed within it. The extra planning structures considered here thus accomplish three goals:
For small targets (see Table 3a), the original target structure is not used as the prescription structure, so the plan is not normalized to its DVH dose during optimization. Instead, the CSV is used as the prescription structure and is prescribed to more than the intended dose. Because the CSV is small, uniformity in it is trivially achieved and its constraints are easily satisfied. The CSV thus serves only to normalize the plan to a known scale. The question then becomes one of adequately covering the treatment target, and controlling the shape of the inhomogeneous dose distribution.
If only the CSV were used as a target, not enough beamlets would be enabled to cover the treatment target well. At least some beamlets that do not directly strike the CSV are needed. These “grazing beamlets” strike the target without striking the CSV and provide coverage of voxels toward the periphery of the treatment target. For this reason, the treatment target itself is also included as a target structure so that sufficient beamlets are enabled at the start of optimization. However, if the “grazing” beamlets are not constrained, the optimizer will adjust them such that uniformity within the target is achieved; therefore, a 3-mm ring is added around the treatment target as a sensitive structure to ensure a steep dose fall-off beyond the target.
The technique for larger targets (see Table 3b for constraints) is to build a pair of structures inside the treatment target to control the dose distribution within the target. The SV is the size of a small target (≤2 mL) to achieve a high dose in the central portion of the target, whereas the smaller CSV inside it allows one to scale the plan. The target becomes a “Sensitive Structure” and in its place a 2-mm-thick ring structure whose outer edge is formed by the periphery of the treatment target is used as a “tumor” structure in the optimization.
As indicated in Table 3a, Table 3b, the CSV structure (or one of the CSV structures in the case of multiple targets) is used for the prescription. In general, for multiple targets the CSV for the largest treatment target should be used for prescription but one can use any other CSV if desired to aid in plan scaling. The ability to achieve target coverage at the 80% prescription isodose level depends on the target volume. If the target is too small to achieve the desired gradient, the maximum achievable target inhomogeneity becomes apparent in the first few iterations of plan optimization and the dose distribution for that target can be scaled to ensure the target is not overdosed.
For plan scaling, the easiest is to use the CSV as a hard constraint in the optimization to scale the entire planned dose to the treatment target up or down as needed. However, in the case of multiple treatment targets, the CSV technique alone works best for the treatment target to which the prescription CSV belongs. For the other treatment targets, it is usually necessary to use a combination of 3 parameters: increasing the importance of the treatment targets, increasing the maximum dose penalty to the ring structure (to scale the dose down), or increasing the minimum dose penalty to the target (to scale the dose up). To maintain control of the plan, it is best to move in small increments when increasing importance and penalty values. If the plan starts moving in an undesirable direction, it was found best to cancel and start the optimization again as opposed to continuing to use the same techniques.
All plans were run until the plan showed that 100% of each target was getting close to the prescribed dose with a reasonable gradient; typically, a few hundred iterations (with each iteration requiring 5−7 sec), but no restrictions were placed on this except that the planning time was kept to a time that would reasonably fit into our proposed radiosurgery workflow. In addition, in many cases increased iterations led to a more homogenous tumor dose so plans were not left to run unattended.
Multiple plans
In the case of multiple treatment targets, it is sometimes beneficial to plan targets separated by a large longitudinal distance as separate plans. Not only does this improve delivery efficiency in the case that there are several unused rotations between tumors during delivery, but it also simplifies the planning. For example, a hard constraint on one target will then not compromise the ability to deliver a desired dose to other targets. If targets in different areas of the brain are prescribed to different doses, allowing for increased pitch for the lower dose treatment targets can further reduce delivery and computation time. The effects of splitting the plan are shown in Table 4 for a subset of patients.
Table 4. Table shows a comparison in delivery time for a single plan and for multiple plans for patients with a large (>5 cm) longitudinal gap between targets
| Patient | No. Plans | Active Rotations | Time (min) | Time Saved (min) |
|---|---|---|---|---|
| 4 | 1 | 71.9 | 72.1 | 40.3 |
| 2 | 33.2 | 31.8 | ||
| 8 | 1 | 70.1 | 58.6 | 28.9 |
| 2 | 29.4 | 29.7 | ||
| 11 | 1 | 58.8 | 54.1 | 16.2 |
| 2 | 38.9 | 37.9 | ||
| 16 | 1 | 85.3 | 69.8 | 39.2 |
| 2 | 32.3 | 30.6 | ||
| 19 | 1 | 59.7 | 59.9 | 31.8 |
| 2 | 28.8 | 28.1 |
Plan quality comparison
Although the detailed plan comparison between the TomoTherapy SRS plans and linear accelerator−based SRS plans is outside the scope of this work, for the purpose of comparing this technique to more established radiosurgery techniques, the prescription isodose volume to tumor volume ratio (PITV),13 conformation number (CN), and a homogeneity index (HI) were calculated for all 59 targets. In addition, the volume of normal tissue receiving 12 Gy was reported as an indication of dose fall-off. The CN, originally described by van't Reit et al.,14 takes into account the relative locations of the prescription isodose volume and the target volume by incorporating target coverage. Currently, both metrics are regularly reported in the literature as objective measures of conformity so both will be considered here to allow for comparison with other work. The PITV and CN are defined below:
(2)
(3)
(4)Because gamma knife radiosurgery is often perceived of as the gold standard of treatment delivery, a comparison is made to the data collected by Nakamura et al.15 for 1338 lesions treated with the gamma knife unit. Nakamura et al.15 used what is referred to in their paper as the “New Conformity Index,” which is the inverse of the conformity index originally proposed by van't Reit et al.14 and Paddick.16 For the purpose of comparing with this study, the inverse CN is considered. To compare the plan homogeneity achieved with this technique to preexisting radiosurgery planning techniques, the average, range, and median HI were also compared with the original linear accelerator–based SRS plans. In addition, 12-Gy volumes were compared with the original linear accelerator–based SRS plans to determine any changes in normal tissues dose between the techniques. The 12-Gy volume has been associated with the risk of symptomatic radiation necrosis (S-NEC) by several authors17, 18, 19, 20 and can be an indication of dose fall-off into normal brain.
Results
Planning technique robustness and delivery parameters
A reliable and efficient planning technique was identified for SRS treatments on TomoTherapy. Final plan details for the 20 patients studied are shown in Table 5. For many patients, delivery time is improved with TomoTherapy over conventional linear accelerator–based SRS treatment plans in which multiple isocenters are used. Duration of linac-SRS delivery can be estimated by multiplying the number of isocenters by the time it takes to deliver each isocenter, which is typically between 15–20 min in our clinic. Thus, for example, a 4-lesion delivery requiring one isocenter per lesion requires a total delivery time of 60–80 min. However, many lesions themselves require more than 1 isocenter, and therefore these estimates represents a lower bound on the delivery time of linac-based SRS. Times are typically more than 120 min for multiple lesions requiring multiple isocenters per lesion. As can be seen from Table 5, the maximum delivery time for a TomoTherapy-based SRS was 66 min, with a median time of 31.2 min for the group of patients studied, which compares favorably with the estimated lower bounds for the delivery times for linac-SRS. Furthermore, it was always possible to deliver the entire prescribed dose in a single treatment fraction with the plan parameters described here.
Table 5. Table shows pitch, modulation factor, number of rotations, total length, delivery time, and beamlet computation time for 10 patients
| Patient | No. of Plans | No. Targets | Max Dose | Pitch | Active Rots. | Mod. Factor | Gantry Period | Delivery Time (min) | DC2 Beamlet Comp Time |
|---|---|---|---|---|---|---|---|---|---|
| 1 | 1 | 1 | 21 | 0.14 | 12.8 | 1.7 | 57 | 12.3 | 0:33:31 |
| 2 | 1 | 4 | 18 | 0.17 | 44.1 | 1.6 | 56 | 41.4 | 2:46:30 |
| 3 | 1 | 3 | 21 | 0.14 | 53.4 | 1.7 | 55 | 49.2 | 2:06:33 |
| 4 | 2 | 3 | 21 | 0.14 | 33.2 | 1.7 | 54, 60 | 31.8 | 1:28:04 |
| 5 | 1 | 3 | 18 | 0.16 | 61.1 | 1.7 | 56 | 57.2 | 1:42:43 |
| 6 | 1 | 6 | 20 | 0.14 | 74.5 | 1.7 | 53 | 66.0 | 2:44:00 |
| 7 | 1 | 3 | 18 | 0.17 | 46.5 | 1.7 | 54 | 42.1 | 1:52:54 |
| 8 | 2 | 4 | 21, 17 | 0.14, 0.17 | 29.3 | 1.5, 1.7 | 60 | 29.7 | 1:08:20 |
| 9 | 1 | 2 | 20 | 0.15 | 30.7 | 1.7 | 53 | 27.3 | 2:06:18 |
| 10 | 1 | 6 | 21 | 0.14 | 60.0 | 1.7 | 60 | 60.0 | 2:37:09 |
| 11 | 2 | 3 | 21 | 0.14 | 38.9 | 1.6 | 58 | 37.9 | 0:51:28 |
| 12 | 1 | 4 | 21 | 0.14 | 50.0 | 1.7 | 50 | 41.9 | 3:02:09 |
| 13 | 1 | 3 | 23 | 0.13 | 29.0 | 1.7 | 60 | 29.2 | 1:46:47 |
| 14 | 1 | 4 | 21 | 0.14 | 47.1 | 1.7 | 57 | 44.9 | 1:03:00 |
| 15 | 1 | 3 | 24 | 0.13 | 46.3 | 1.7 | 60 | 46.5 | 0:59:03 |
| 16 | 2 | 2 | 21, 18 | 0.13 | 32.3 | 1.7 | 56 | 30.6 | 1:57:19 |
| 17 | 1 | 2 | 21 | 0.13 | 19.2 | 1.7 | 53 | 17.1 | 2:08:36 |
| 18 | 1 | 1 | 15 | 0.18 | 16.7 | 1.7 | 55 | 15.5 | 1:22:36 |
| 19 | 2 | 2 | 24 | 0.13 | 28.7 | 1.7 | 56, 60 | 28.1 | 1:04:57 |
| 20 | 1 | 1 | 21 | 0.14 | 13.9 | 1.7 | 55 | 13 | 0:47:15 |
TomoTherapy delivery time is quite predictable, given that this is simply the number of gantry rotations multiplied by the gantry period, which, as already mentioned can be manipulated via adjusting the pitch to be in the range of 50–60 sec. Table 2 indicates that pitches between 0.13 and 0.18 give a gantry period in this range for delivered doses of 15–24 Gy (with a pitch of 0.13, corresponding to the high end of the dose range). It would appear that even for metastases scattered throughout the longitudinal extent of the brain with no gaps to allow for multiple plans, the treatment time would be unlikely to exceed 90 min, rivaling alternate delivery methods.
Finally, beamlet computation time with the 16-node cluster was not more than a few hours for any patient and should fit into a radiosurgery workflow. It is anticipated that this time will only decrease as computational capabilities improve over time.
Comparison with conventional radiosurgery
To determine the credibility of any planning technique, it is important to ensure that it reliably produces plans of radiosurgical quality. We aimed to achieve with the TomoTherapy system a conformity that is comparable to preexisting radiosurgery systems with dose fall-off that is as steep as possible. To measure conformity, data for the 59 lesions studied here is compared with data reported by Nakamura et al.15 for 1338 lesions treated with gamma knife radiosurgery. They report median target coverage, median PITV, and median inverse CNs of 0.97, 1.67, and 1.78, respectively. For the TomoTherapy dataset, these median values were 0.99, 1.59, and 1.66, respectively. In addition Nakamura et al.15 split their data into 4-volume quartiles and reported median PITV values for each. The first quartile included volumes ranging from 0.01–0.25 mL, the second includes 0.26–1.5 mL, the third includes 1.6–5.5 mL, and the fourth includes 5.6–55.9 mL. They report median PITV ratios of 3.43,1.85,1.53, and 1.35 for the 4 quartiles, starting with the first quartile. Using the same volume cutoffs, we report median ratios of 2.26, 1.51, 1.44, and 1.18 for the same 4 quartiles, which compare favorably to the gamma knife values. Moreover, these results also compare favorably with a study by Gutierrez et al.,21 who made the same comparison for brain metastases treated in combination with whole-brain radiation therapy. They report PITV values of 2.22, 1.34, and 1.17 for the second, third, and fourth quartiles, respectively. Using the CN, we also looked at average values for 3 additional subgroups: <0.1 mL, 0.1–1 mL, and >1 mL. Average CNs for the 3 groups were 3.84, 1.88, and 1.46, respectively. Both results show an improvement of conformity with volume. Although comparatively, worse conformity is observed for very small metastases, Nakamura et al.15 note that because no complications were observed for patients with smaller-volume lesions (<1 mL), conformity is not as critical as it is for larger lesions, where a much larger volume of normal brain is irradiated with decreasing conformity.
Homogeneity was compared with the original linear accelerator–based plans using circular collimators to show that the planning technique described here does allow for adequate inhomogeneity to increase peripheral dose fall-off. Results are shown in Table 6. Achieved homogeneity indices compared favorably with linear accelerator–based SRS. Median homogeneity for the linac-based SRS and TomoTherapy plans, respectively, were 1.25 and 1.24, with average values of 1.30 and 1.24.
Table 6. Comparison of the homogeneity index (maximum dose/prescription dose) achieved for a linear accelerator–based planning approach using circular collimators and helical TomoTherapy
| Linac | Tomotherapy | |
|---|---|---|
| Maximum HI | 2.02 | 1.52 |
| Minimum HI | 1.06 | 1.07 |
| Median HI | 1.25 | 1.24 |
| Average HI (± SD) | 1.30 | 1.24 |
The 12-Gy volume has been postulated to correlate with the risk of S-NEC.17, 18, 19, 20 Because of the small size, and predominantly spherical shape of the targets in this patient cohort, the changes in low-dose volumes were not dramatic. The average increase for the 60 targets was 1.1 mL (range –1.9 mL to +5.2 mL), with an average total increase of 3.4 mL (range –1.9 mL to +12.1 mL) for the 20 patients. Clearly, the increase in the 12-Gy volume is dependent on number and size of the targets. The effect of these changes on relative risk of S-NEC will be explored in future work.
Discussion
The work presented above shows that our proposed TomoTherapy SRS planning technique allows one to develop plans for multiple brain metastases that can be delivered efficiently in a time frame that is comparable to or shorter than what is currently achievable using linac-based SRS or gamma knife SRS. Multiple sites can be treated within a single plan and there is no need to have each target positioned at the isocenter, as with conventional SRS techniques.
This work only considers the case of brain metastases, which are often well circumscribed spherical targets. It is possible that the planning technique described here might not be applicable for other diseases treated with SRS. For instance, for lesions where dose homogeneity is important (e.g., those infiltrating normal brain) the special techniques described here to achieve a greater degree of inhomogeneity may not be useful. In addition, targets with a complex shape might require higher modulation, and hence pitch will have to be reduced to allow for delivery in a single fraction.
It is important to note that this technique applies to the TomoTherapy optimization algorithm as currently implemented. If there are changes to the algorithm described in Eq (1), the described technique will have to be scrutinized and a new planning technique may have to be developed. All beamlet calculation and optimization time estimates only pertain to calculations made with the 16-node cluster currently in operation at our clinic. With a newer cluster design, it would be expected that dose calculation times would be shorter. It may even be possible to calculate on 512 × 512 resolution images, reducing the voxel size and improving plan accuracy. We found that we could not increase resolution and calculate beamlets in a reasonable time frame with our current cluster configuration. The addition of dynamic jaws as described by Chen et al.22 could improve delivery efficiency by allowing for increased couch speed between targets in the case of multiple targets. Future work will aim to look at the impact of dynamic jaws on delivery efficiency.
Because all of the plans here were preformed using CT datasets for patients treated with linear accelerator–based SRS, detailed comparison of these plans to the clinically used linac-based SRS plans will be considered in future work and is beyond the scope of this technical paper. Here, we report on the detailed development of the planning technique for this nonstandard application of the TomoTherapy system.
Conclusions
This work describes the development and evaluation of a new planning technique for SRS for brain metastases using TomoTherapy. An efficient planning technique has been described that will in turn allow for efficient treatment delivery. We have presented here one possible solution, of a multitude of possible solutions, which yields reliable results that are consistent with SRS planning goals. This work also establishes the feasibility of a one-day SRS workflow on TomoTherapy because planning time can be kept to 2 hours or less for a number of targets, which compares well with planning times necessary for multiple targets in conventional linear accelerator–based SRS.
Acknowledgment
Supported in part by NIH R01 109656 and R01 118365.
References
- Quality assurance of an image guided intracranial stereotactic positioning system. Technol. Cancer. Res. Treat. 2009;8:39–49
- TomoTherapy for cranial radiosurgery/radiotherapy. IFMBE Proc. 2007;14:1852–1855
- Stereotactic image-guided intensity modulated radiotherapy using the Hi-ART II helical tomotherapy system. Med. Dosim. 2008;33:135–148
- Tomotherapy planning of small brain tumors. Radiother. Oncol. 2005;74:49–52
- Comparison of advanced irradiation techniques with photons for benign intracranial tumours. Radiother. Oncol. 2006;80:268–273
- . Target dose versus extratarget dose in stereotactic radiosurgery. Acta. Oncol. 1993;32:453–457
- Dosimetric comparisons of helical tomotherapy treatment plans and step-and-shoot intensity-modulated radiosurgery treatment plans in intracranial stereotactic radiosurgery. Int. J. Radiat. Oncol. Biol. Phys. 2006;65:608–616
- Integrated subvolume boosting technique to reduce normal Tissue irradiation during SBRT treatment planning using helical tomotherapy. Med. Phys. 2008;25:2957
- TomoTherapy Incorporated (Field Safety Notice 3449). http://tomoexchange.tomotherapy.com/documentation/category/product_advisory_notices2008;Accessed December 16, 2009
- Acceleration of tomotherapy treatment delivery by increasing pitch and decreasing modulation. Med. Phys. 2008;35:2636–2637
- On the making of sharp longitudinal dose profiles with helical tomotherapy. Phys. Med. Biol. 2007;52:6497–6510
- Tomo® Planning Guide Version 3.x. Madison, WI: TomoTherapy Inc; 2007;
- Radiation Therapy Oncology Group: Radiosurgery Quality Assurance Guidelines. Int. J. Radiat. Oncol. Biol. Phys. 1993;27:1231–1239
- A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. Int. J. Radiat. Oncol. Biol. Phys. 1997;37:731–736
- Dose conformity of gamma knife radiosurgery and risk factors for complications. Int. J. Radiat. Oncol. Biol. Phys. 2001;51:1313–1319
- . A simple scoring ratio to index the conformity of radiosurgical treatment plans (Technical note). J. Neurosurg. 2000;93:219–222
- . Estimation of complications for linear accelerator radiosurgery with the integrated logistic formula. Int. J. Radiat. Oncol. Biol. Phys. 1990;19:143–148
- Development of a model to predict permanent symptomatic postradiosurgery injury for arteriovenous malformation patients (Arteriovenous Malformation Radiosurgery Study Group). Int. J. Radiat. Oncol. Biol. Phys. 2000;46:1143–1148
- . Intracranial lesions made by the Bragg peak of a proton beam. In: Response of the nervous system to ionizing irradiation. Boston: Little, Brown, and Co; 1964;
- 12 Gy gamma knife radiosurgical volume is a predictor for radiation necrosis in non-AVM intracranial tumors. Int. J. Radiat. Oncol. Biol. Phys. 2006;64:419–424
- Whole brain radiotherapy with hippocampal avoidance and simultaneously integrated brain metastases boost: a planning study. Int. J. Radiat. Oncol. Biol. Phys. 2007;69:589–597
- Dynamic jaw for helical tomotherapy delivery. Med. Phys. 2008;35:2744–2745
PII: S0958-3947(09)00129-0
doi:10.1016/j.meddos.2009.11.003
© 2011 American Association of Medical Dosimetrists. Published by Elsevier Inc. All rights reserved.
