Investigation of effect of variations in bone fraction and red marrow cellularity on bone marrow dosimetry in radio-immunotherapy.

A method is described for computing patient-specific absorbed dose rates to active marrow which accounts for spatial variation in bone volume fraction and marrow cellularity. A module has been added to the 3D Monte Carlo dosimetry program DPM to treat energy deposition in the components of bone spongiosa distinctly. Homogeneous voxels in regions containing bone spongiosa (as defined on CT images) are assumed to be comprised only of bone, active (red) marrow and inactive (yellow) marrow. Cellularities are determined from biopsy, and bone volume fractions are computed from cellularities and CT-derived voxel densities. Electrons are assumed to deposit energy locally in the three constituent components in proportions determined by electron energy absorption fractions which depend on energy, cellularity, and bone volume fraction, and which are either taken from the literature or are derived from Monte Carlo simulations using EGS5. Separate algorithms are used to model primary ß particles and secondary electrons generated after photon interactions. Treating energy deposition distinctly in bone spongiosa constituents leads to marrow dosimetry results which differ from homogeneous spongiosa dosimetry by up to 20%. Dose rates in active marrow regions with cellularities of 20, 50, and 80% can vary by up to 20%, and can differ by up to 10% as a function of bone volume fraction. Dose to bone marrow exhibits a strong dependence on marrow cellularity and a potentially significant dependence on bone volume fraction.

Dosimetry in patients with B-cell lymphoma treated with [(90)Y]ibritumomab tiuxetan or [(131)I]tositumomab.

Radioimmunotherapy involves the use of radiolabeled monoclonal antibodies (MAbs) to treat malignancy. The therapeutic effect is determined by the radiopharmaceutical, the radiation absorbed dose and previous treatments. There are currently two approved radiopharmaceuticals for the treatment of B-cell lymphoma - the (90)Y-labeled ibritumomab and the (131)I-labeled tositumomab. Both are directed against CD20, albeit not against the same epitope. This paper summarizes current results of dose-responses for normal tissues and tumours of [(131)I]tositumomab and [(90)Y]ibritumomab tiuxetan, discusses them in the context of dosimetry methods used and highlights the assumptions being made in the different dosimetry methodologies. Moreover, we wish to point at the possibility of performing low-cost therapy bremsstrahlung imaging for [(90)Y]ibritumomab tiuxetan to confirm biodistribution, and possibly also for dosimetric calculations.

Method for Fast CT/SPECT-Based 3D Monte Carlo Absorbed Dose Computations in Internal Emitter Therapy.

The DPM (Dose Planning Method) Monte Carlo electron and photon transport program, designed for fast computation of radiation absorbed dose in external beam radiotherapy, has been adapted to the calculation of absorbed dose in patient-specific internal emitter therapy. Because both its photon and electron transport mechanics algorithms have been optimized for fast computation in 3D voxelized geometries (in particular, those derived from CT scans), DPM is perfectly suited for performing patient-specific absorbed dose calculations in internal emitter therapy. In the updated version of DPM developed for the current work, the necessary inputs are a patient CT image, a registered SPECT image, and any number of registered masks defining regions of interest. DPM has been benchmarked for internal emitter therapy applications by comparing computed absorption fractions for a variety of organs using a Zubal phantom with reference results from the Medical Internal Radionuclide Dose (MIRD) Committee standards. In addition, the ß decay source algorithm and the photon tracking algorithm of DPM have been further benchmarked by comparison to experimental data. This paper presents a description of the program, the results of the benchmark studies, and some sample computations using patient data from radioimmunotherapy studies using (131)I.

Monte Carlo evaluation of object shape effects in iodine-131 SPET tumor activity quantification.

In our clinical iodine-131 single-photon emission tomography (SPET) quantification for radioimmunotherapy, calibration and partial volume correction are based on measurements with phantoms containing spheres to simulate patient tumors even though real tumors are frequently nonspherical. In this study, Monte Carlo simulation was used to evaluate how object shape influences "spill-out" and "spill-in", which are major sources of quantification error associated with the poor spatial resolution of 131I SPET. Objects that varied in shape (spheres, cylinders, and an irregular structure) but were identical in activity and volume were simulated. Iterative reconstruction employed both attenuation and triple-energy-window scatter compensation. VOIs were defined in the reconstructed images both using physical boundaries and using expanded boundaries to allow for the limited resolution. When physical boundaries were used, both spill-out and spill-in were more significant for nonspherical structures than for spherical structures. Over the range of object volumes (50-200 ml) and at all background levels, VOI counts in cylinders were lower than VOI counts in spheres. This underestimation increased with decrease in object size (for the cold background -18% at 200 ml and -39% at 50 ml). It also decreased with increase in background activity because spill-in partially compensated for spill-out. It was shown that with a VOI larger than physical size, the results are independent of object shape and size only in the case of cold background. Activity quantification was carried out using a procedure similar to that used in our clinic. Quantification of nonspherical objects was improved by simple sphere-based partial volume correction, but the error was still large in some cases (for example, -39% for a 50-ml cylinder in a cold background and -35% for a 200-ml irregular structure defined on the basis of a typical tumor outlined on an X-ray computed tomography scan of a patient with non-Hodgkin's lymphoma). Partial volume correction by patient-specific Monte Carlo simulation may provide better quantification accuracy.

Accuracy of 131I tumor quantification in radioimmunotherapy using SPECT imaging with an ultra-high-energy collimator: Monte Carlo study.

UNLABELLED: Accuracy of 131I tumor quantification after radioimmunotherapy (RIT) was investigated for SPECT imaging with an ultra-high-energy (UHE) collimator designed for imaging 511-keV photons. METHODS: First, measurements and Monte Carlo simulations were carried out to compare the UHE collimator with a conventionally used, high-energy collimator. On the basis of this comparison, the UHE collimator was selected for this investigation, which was carried out by simulation of spherical tumors in a phantom. Reconstruction was by an expectation-maximization algorithm that included scatter and attenuation correction. Keeping the tumor activity constant, simulations were carried out to assess how volume-of-interest (VOI) counts vary with background activity, radius of rotation (ROR), tumor location, and size. The constant calibration factor for quantification was determined from VOI counts corresponding to a 3.63-cm-radius sphere of known activity. Tight VOIs corresponding to the physical size of the spheres or tumors were used. RESULTS: Use of the UHE collimator resulted in a large reduction in 131I penetration, which is especially significant in RIT where background uptake is high. With the UHE collimator, typical patient images showed an improvement in contrast. Considering the desired geometric events, sensitivity was reduced, but only by a factor of 1.6. Simulation results for a 3.63-cm-radius tumor showed that VOI counts vary with background, location, and ROR by less than 3.2%, 3%, and 5.3%, respectively. The variation with tumor size was more significant and was a function of the background. Good quantification accuracy (<6.5% error) was achieved when tumor size was the same as the sphere size used in the calibration, irrespective of the other parameters. For smaller tumors, activities were underestimated by up to -15% for the 2.88-cm-radius sphere, -23% for the 2.29-cm-radius sphere, and -47% for the 1.68-cm-radius sphere. CONCLUSION: Reasonable accuracy can be achieved for VOI quantification of 131I using SPECT with an UHE collimator and a constant calibration factor. Difference in tumor size relative to the size of the calibration sphere had the biggest effect on accuracy, and recovery coefficients are needed to improve quantification of small tumors.

Characterization of scatter and penetration using Monte Carlo simulation in 131I imaging.

UNLABELLED: In 131I SPECT, image quality and quantification accuracy are degraded by object scatter as well as scatter and penetration in the collimator. The characterization of energy and spatial distributions of scatter and penetration performed in this study by Monte Carlo simulation will be useful for the development and evaluation of techniques that compensate for such events in 131I imaging. METHODS: First, to test the accuracy of the Monte Carlo model, simulated and measured data were compared for both a point source and a phantom. Next, simulations to investigate scatter and penetration were performed for four geometries: point source in air, point source in a water-filled cylinder, hot sphere in a cylinder filled with nonradioactive water, and hot sphere in a cylinder filled with radioactive water. Energy spectra were separated according to order of scatter, type of interaction, and gamma-ray emission energy. A preliminary evaluation of the triple-energy window (TEW) scatter correction method was performed. RESULTS: The accuracy of the Monte Carlo model was verified by the good agreement between measured and simulated energy spectra and radial point spread functions. For a point source in air, simulations show that 73% of events in the photopeak window had either scattered in or penetrated the collimator, indicating the significance of collimator interactions. For a point source in a water-filled phantom, the separated energy spectra showed that a 20% photopeak window can be used to eliminate events that scatter more than two times in the phantom. For the hot sphere phantoms, it was shown that in the photopeak region the spectrum shape of penetration events is very similar to that of primary (no scatter and no penetration) events. For the hot sphere regions of interest, the percentage difference between true scatter counts and the TEW estimate of scatter counts was <12%. CONCLUSION: In 131I SPECT, object scatter as well as collimator scatter and penetration are significant. The TEW method provides a reasonable correction for scatter, but the similarity between the 364-keV primary and penetration energy spectra makes it difficult to compensate for these penetration events using techniques that are based on spectral analysis.