A paired-image radiation transport model for skeletal dosimetry

Toxicity of the hematopoietically active bone marrow continues to be a primary limitation in radionuclide therapies of cancer. Improved techniques for patient-specific skeletal dosimetry are thus crucial to the development of dose-response relationships needed to optimize these therapies (i.e., avoid both marrow toxicity and tumor underdosing). Current clinical methods of skeletal dose assessment rely heavily on a single set of bone and marrow cavity chord-length distributions in which particle energy deposition is tracked within an infinite extent of trabecular spongiosa, with no allowance for particle escape to cortical bone. In the present study, we introduce a paired-image radiation transport (PIRT) model that can provide a more realistic 3-dimensional geometry for particle transport of the skeletal site at both microscopic and macroscopic levels of its histology.

Ex vivo CT scans were acquired of the lumbar vertebra and right proximal femur excised from a 66-y male cadaver (body mass index, 22.7 kg m(-2)). For both skeletal sites, regions of trabecular spongiosa and cortical bone were identified and segmented. Physical sections of interior spongiosa were then taken and subjected to nuclear magnetic resonance (NMR) microscopy. Voxels within the resulting NMR microimages were segmented and labeled into regions of bone trabeculae, endosteum, active marrow, and inactive marrow. The PIRT methodology was then implemented within the EGSnrc radiation transport code, whereby electrons of various initial energies are simultaneously tracked within both the ex vivo CT macroimage and the NMR microimage of the skeletal site.

At electron initial energies greater than 50-200 keV, a divergence in absorbed fractions to active marrow is noted between PIRT model simulations and those estimated under infinite spongiosa transport techniques. Calculations of radionuclide S values under both methodologies imply that current chord-based models used in clinical skeletal dosimetry can overestimate dose to active bone marrow in these 2 skeletal sites by approximately 4%-23% for low-energy beta-emitters ((33)P, (169)Er, and (177)Lu), by approximately 4%-25% for intermediate-energy beta-emitters ((153)Sm, (186)Re, and (89)Sr), and by approximately 11%-30% for high-energy beta-emitters ((32)P, (188)Re, and (90)Y).

The PIRT methodology allows for detailed modeling of the 3D macrostructure of individual marrow-containing bones within the skeleton, thus permitting improved estimates of absorbed fractions and radionuclide S values for intermediate-to-high beta-emitters.