We have developed a time-dependent computational framework, hematological dose (HEDOS), to estimate dose to circulating blood cells from radiation therapy treatment fields for any treatment site. Two independent dynamic models were implemented in HEDOS: one describing the spatiotemporal distribution of blood particles (BPs) in organs and the second describing the time-dependent radiation field delivery. A whole-body blood flow network based on blood volumes and flow rates from ICRP Publication 89 was simulated to produce the spatiotemporal distribution of BPs in organs across the entire body using a discrete-time Markov process. Constant or time-varying transition probabilities were applied and their impact on transition time was investigated. The impact of treatment time and anatomical site were investigated using imaging data and dose distributions from a liver cancer and a brain cancer patient. The simulations revealed different dose levels to the circulating blood for brain irradiation compared to liver irradiation even for similar field sizes due to the different blood flow properties of the two organs. The volume of blood receiving any dose (V >0 Gy) after a single radiation fraction increases from 1.2% for a 1 s delivery time to 20.9% for 120 s delivery time for the brain cancer treatment, and from 10% (1 s) to 48.7% (120 s) for a liver cancer treatment. Other measures of the low-dose bath to the circulating blood such as the dose to small volumes of blood (D 2%) decreases with longer delivery time. Furthermore, we demonstrate that the blood dose-volume histogram is highly sensitive to changes in the treatment time, indicating that dynamic modeling of blood flow and radiation fields is necessary to evaluate dose to circulating blood cells for the assessment of radiation-induced lymphopenia. HEDOS is publicly available and allows for the estimation of patient-specific dose to circulating blood cells based on organ DVHs, thus enabling the study of the impact of different treatment plans, dose rates, and fractionation schemes.

We have developed a novel 4D dynamic liver blood flow model, capable of accurate dose estimation to circulating blood cells during liver-directed external beam radiotherapy, accounting for blood recirculation and radiation delivery time structure. Adult male and adult female liver computational phantoms with detailed vascular trees were developed to include the hepatic arterial, hepatic portal venous, and hepatic venous trees. A discrete time Markov Chain approach was applied to determine the spatiotemporal distribution of 105 blood particles (BP) in the human body based on reference values for cardiac output and organ blood volumes. For BPs entering the liver, an explicit Monte Carlo simulation was implemented to track their propagation along ~2000 distinct vascular pathways through the liver. The model tracks accumulated absorbed dose from time-dependent radiation fields with a 0.1s time resolution. The computational model was then evaluated for 3 male and 3 female patients receiving photon (VMAT and IMRT) and proton (passive SOBP and active PBS) treatments. The dosimetric impact of treatment modality, delivery time, and fractionation on circulating blood cells was investigated and quantified using the mean dose (μdose,b), V>0Gy, V>0.125Gy, and D2%. Average reductions in μdose,b, V>0Gy, V>0.125Gy, and D2% of 45%, 6%, 53%, 19% respectively, were observed for proton treatments as compared to photon treatments. Our simulation also showed that V>0Gy, V>0.125Gy, and D2% were highly sensitive to the beam-on time. Both V>0Gy and V>0.125Gy increased with beam-on time, whereas D2% decreased with increasing beam-on time, demonstrating the tradeoff between low dose to a large fraction of blood cells and high dose to a small fraction of blood cells. Consequently, proton treatments are not necessarily advantageous in terms of dose to the blood simply based on integral dose considerations. Instead, both integral dose and beam-on time can substantially impact relevant dosimetric indices.

Purpose To develop a model of the internal vasculature of the adult liver and demonstrate its application to the differentiation of radiopharmaceutical decay sites within liver parenchyma from those within organ blood. Method Computer-generated models of hepatic arterial (HA), hepatic venous (HV), and hepatic portal venous (HPV) vascular trees were algorithmically created within individual lobes of the ICRP adult female and male livers (AFL/AML). For each iteration of the algorithm, pressure, blood flow, and vessel radii within each tree were updated as each new vessel was created and connected to a viable bifurcation site. The vascular networks created inside the AFL/AML were then tetrahedralized for coupling to the PHITS radiation transport code. Specific absorbed fractions (SAF) were computed for monoenergetic alpha particles, electrons, positrons, and photons. Dual-region liver models of the AFL/AML were proposed, and particle-specific SAF values were computed assuming radionuclide decays in blood within two locations: (1) sites within explicitly modeled hepatic vessels, and (2) sites within the hepatic blood pool residing outside these vessels to include the capillaries and blood sinuses. S values for 22 and 10 radionuclides commonly used in radiopharmaceutical therapy and imaging, respectively, were computed using the dual-region liver models and compared to those obtained in the existing single-region liver model. Results Liver models with virtual vasculatures of ~ 6000 non-intersecting straight cylinders representing the HA, HPV, and HV circulations were created for the ICRP reference. For alpha emitters and for beta and auger-electron emitters, S values using the single-region models were approximately 11% (AML) to 14% (AFL) and 11% (AML) to 13% (AFL) higher than the S values obtained using the dual-region models, respectively. Conclusions The methodology employed in this study has shown improvements in organ parenchymal dosimetry through explicit consideration of blood self-dose for alpha particles (all energies) and for electrons at energies below ~ 100 keV.

No abstract

Objective: Phantoms of the International Commission on Radiological Protection provide a framework for standardized dosimetry. The modeling of internal blood vessels – essential to tracking circulating blood cells exposed during external beam radiotherapy and to account for radiopharmaceutical decays while still in blood circulation – is, however, limited to the major inter-organ arteries and veins. Intra-organ blood is accounted for only through the assignment of a homogeneous mixture of parenchyma and blood [single-region (SR) organs]. Our goal was to develop explicit dual-region (DR) models of intra-organ blood vasculature of the adult male brain (AMB) and adult female brain (AFB). Approach: A total of 4000 vessels were created amongst 26 vascular trees. The AMB and AFB models were then tetrahedralized for coupling to the PHITS radiation transport code. Absorbed fractions were computed for monoenergetic alpha particles, electrons, positrons, and photons for both decay sites within the blood vessels and for tissues outside these vessels. Radionuclide S-values were computed for 22 and 10 radionuclides commonly employed in radiopharmaceutical therapy and nuclear medicine diagnostic imaging, respectively. Main Results: For radionuclide decays, values of S(brain tissue ← brain blood) assessed in the traditional manner (SR) were higher than those computed using our DR models by factors of 1.92, 1.49, and 1.57 for therapeutic alpha-emitters, beta-emitters, and auger electron-emitters, respectively in the AFB and by factors of 1.65, 1.37, and 1.42 for these same radionuclide categories in the AMB. Corresponding ratios of SR and DR values of S(brain tissue ← brain blood) were 1.34 (AFB) and 1.26 (AMB) for four SPECT radionuclides, and were 1.32 (AFB) and 1.24 (AMB) for six common PET radionuclides. Significance: The methodology employed in this study can be explored in other organs of the body for proper accounting of blood self-dose for that fraction of the radiopharmaceutical still in general circulation.