Technical Note: Monte Carlo calculations of the AAPM TG-43 brachytherapy dosimetry parameters for a new titanium-encapsulated Yb-169 source

Brachytherapy at the nanoscale with protein functionalized and intrinsically radiolabeled [169Yb]Yb2O3 nanoseeds

PURPOSE

Classical brachytherapy of solid malignant tumors is an invasive procedure which often results in an uneven dose distribution, while requiring surgical removal of sealed radioactive seed sources after a certain period of time. To circumvent these issues, we report the synthesis of intrinsically radiolabeled and gum Arabic glycoprotein functionalized [169Yb]Yb2O3 nanoseeds as a novel nanoscale brachytherapy agent, which could directly be administered via intratumoral injection for tumor therapy.

METHODS

169Yb (T½ = 32 days) was produced by neutron irradiation of enriched (15.2% in 168Yb) Yb2O3 target in a nuclear reactor, radiochemically converted to [169Yb]YbCl3 and used for nanoparticle (NP) synthesis. Intrinsically radiolabeled NP were synthesized by controlled hydrolysis of Yb3+ ions in gum Arabic glycoprotein medium. In vivo SPECT/CT imaging, autoradiography, and biodistribution studies were performed after intratumoral injection of radiolabeled NP in B16F10 tumor bearing C57BL/6 mice. Systematic tumor regression studies and histopathological analyses were performed to demonstrate therapeutic efficacy in the same mice model.

RESULTS

The nanoformulation was a clear solution having high colloidal and radiochemical stability. Uniform distribution and retention of the radiolabeled nanoformulation in the tumor mass were observed via SPECT/CT imaging and autoradiography studies. In a tumor regression study, tumor growth was significantly arrested with different doses of radiolabeled NP compared to the control and the best treatment effect was observed with ~ 27.8 MBq dose. In histopathological analysis, loss of mitotic cells was apparent in tumor tissue of treated groups, whereas no significant damage in kidney, lungs, and liver tissue morphology was observed.

CONCLUSIONS

These results hold promise for nanoscale brachytherapy to become a clinically practical treatment modality for unresectable solid cancers.

A re-activation model for 169Yb intensity modulated brachytherapy sources accounting for spatiotemporal isotopic composition

BACKGROUND

Intensity modulated brachytherapy based on partially shielded intracavitary and interstitial applicators is possible with a cost-effective 169Yb production method. 169Yb is a traditionally expensive isotope suitable for this purpose, with an average γ-ray energy of 93 keV. Re-activating a single 169Yb source multiple times in a nuclear reactor between clinical uses was shown to theoretically reduce cost by approximately 75% relative to conventional single-activation sources. With re-activation, substantial spatiotemporal variation in isotopic source composition is expected between activations via 168Yb burnup and 169Yb decay, resulting in time dependent neutron transmission, precursor usage, and reactor time needed per re-activation.

PURPOSE

To introduce a generalized model of radioactive source production that accounts for spatiotemporal variation in isotopic source composition to improve the efficiency estimate of the 169Yb production process, with and without re-activation.

METHODS AND MATERIALS

A time-dependent thermal neutron transport, isotope transmutation, and decay model was developed. Thermal neutron flux within partitioned sub-volumes of a cylindrical active source was calculated by raytracing through the spatiotemporal dependent isotopic composition throughout the source, accounting for thermal neutron attenuation along each ray. The model was benchmarked, generalized, and applied to a variety of active source dimensions with radii ranging from 0.4 to 1.0 mm, lengths from 2.5 to 10.5 mm, and volumes from 0.31 to 7.85 mm3, at thermal neutron fluxes from 1 × 1014 to 1 × 1015 n cm-2 s-1. The 168Yb-Yb2O3 density was 8.5 g cm-3 with 82% 168Yb-enrichment. As an example, a reference re-activatable 169Yb active source (RRS) constructed of 82%-enriched 168Yb-Yb2O3 precursor was modeled, with 0.6 mm diameter, 10.5 mm length, 3 mm3 volume, 8.5 g cm-3 density, and a thermal neutron activation flux of 4 × 1014 neutrons cm-2 s-1.

RESULTS

The average clinical 169Yb activity for a 0.99 versus 0.31 mm3 source dropped from 20.1 to 7.5 Ci for a 4 × 1014 n cm-2 s-1 activation flux and from 20.9 to 8.7 Ci for a 1 × 1015 n cm-2 s-1 activation flux. For thermal neutron fluxes ≥2 × 1014 n cm-2 s-1, total precursor and reactor time per clinic-year were maximized at a source volume of 0.99 mm3 and reached a near minimum at 3 mm3. When the spatiotemporal isotopic composition effect was accounted for, average thermal neutron transmission increased over RRS lifetime from 23.6% to 55.9%. A 28% reduction (42.5 days to 30.6 days) in the reactor time needed per clinic-year for the RRS is predicted relative to a model that does not account for spatiotemporal isotopic composition effects.

CONCLUSIONS

Accounting for spatiotemporal isotopic composition effects within the RRS results in a 28% reduction in the reactor time per clinic-year relative to the case in which such changes are not accounted for. Smaller volume sources had a disadvantage in that average clinical 169Yb activity decreased substantially below 20 Ci for source volumes under 1 mm3. Increasing source volume above 3 mm3 adds little value in precursor and reactor time savings and has a geometric disadvantage.

Dosimetric impact of a robust optimization approach to mitigate effects from rotational uncertainty in prostate intensity-modulated brachytherapy

BACKGROUND

Intensity-modulated brachytherapy (IMBT) is an emerging technology for cancer treatment, in which radiation sources are shielded to shape the dose distribution. The rotatable shields provide an additional degree of freedom, but also introduce an additional, directional, type of uncertainty, compared to conventional high-dose-rate brachytherapy (HDR BT).

PURPOSE

We propose and evaluate a robust optimization approach to mitigate the effects of rotational uncertainty in the shields with respect to planning criteria.

METHODS

A previously suggested prototype for platinum-shielded prostate ¹⁶⁹ Yb-based dynamic IMBT is considered. We study a retrospective patient data set (anatomical contours and catheter placement) from two clinics, consisting of six patients that had previously undergone conventional ¹⁹² Ir HDR BT treatment. The Monte Carlo-based treatment planning software RapidBrachyMCTPS is used for dose calculations. In our computational experiments, we investigate systematic rotational shield errors of ±10° and ±20°, and the same systematic error is applied to all dwell positions in each scenario. This gives us three scenarios, one nominal and two with errors. The robust optimization approach finds a compromise between the average and worst-case scenario outcomes.

RESULTS

We compare dose plans obtained from standard models and their robust counterparts. With dwell times obtained from a linear penalty model (LPM), for 10° errors, the dose to urethra ( D 0.1 c c $D_{0.1cc}$ ) and rectum ( D 0.1 c c $D_{0.1cc}$ and D 1 c c $D_{1cc}$ ) increase with up to 5% and 7%, respectively, in the worst-case scenario, while with the robust counterpart, the corresponding increases were 3% and 3%. For all patients and all evaluated criteria, the worst-case scenario outcome with the robust approach had lower deviation compared to the standard model, without compromising target coverage. We also evaluated shield errors up to 20° and while the deviations increased to a large extent with the standard models, the robust models were capable of handling even such large errors.

CONCLUSIONS

We conclude that robust optimization can be used to mitigate the effects from rotational uncertainty and to ensure the treatment plan quality of IMBT.

169 Yb-based rotating shield brachytherapy for prostate cancer

PURPOSE

To present a system for the treatment of prostate cancer in a single-fraction regimen using ¹⁶⁹ Yb-based rotating shield brachytherapy (RSBT) with a single-catheter robotic delivery system. The proposed system is innovative because it can deliver RSBT through multiple implanted needles independently, in serial, using flexible catheters, with no inter-needle shielding effects and without the need to rotate multiple shielded catheters inside the needles simultaneously, resulting in a simple, mechanically robust, delivery approach. RSBT was compared to conventional ¹⁹² Ir-based high-dose-rate brachytherapy (HDR-BT) in a treatment planning study with dose escalation and urethral sparing goals, representing single-fraction brachytherapy monotherapy and brachytherapy as a boost to external beam radiotherapy, respectively. A prototype mechanical delivery system was constructed and quantitatively evaluated as a proof of concept.

METHODS

Treatment plans for twenty-six patients with single fraction prescriptions of 20.5 and 15 Gy, were created for dose escalation and urethral sparing, respectively. The RSBT and HDR-BT delivery systems were modeled with one partially shielded 999 GBq (27 Ci) ¹⁶⁹ Yb source and one 370 GBq (10 Ci) ¹⁹² Ir source, respectively. A prototype angular drive system for helical source delivery was constructed. Mechanical accuracy measurements of source translational position and angular orientation in a simulated treatment delivery setup were obtained using the prototype system.

RESULTS

For dose escalation, with equivalent urethra D_10% , PTV D_90% for RSBT vs HDR-BT increased from 22.6 ± 0.0 Gy (average ± standard deviation) to 29.3 ± 0.9 Gy, or 29.9 % ± 3.0%, with treatment times of 51.4 ± 6.1 min for RSBT and 15.8 ± 2.3 min for 10 Ci ¹⁹² Ir-based HDR-BT. For urethra sparing, with equivalent PTV D_{90 %} , urethra D_10% for RSBT vs HDR-BT decreased for RSBT vs HDR-BT from 15.6 ± 0.4 Gy to 12.0 ± 0.4 Gy, or 23.1% ± 3.5%, with treatment times of 30.0 ± 3.7 min for RSBT and 12.3 ± 1.8 min for HDR-BT. Differences between measured vs predicted rotating catheter positions (corresponding to source position) were within 0.18 mm ± 0.12 mm longitudinally and 0.07° ± 0.78°.

CONCLUSION

¹⁶⁹ Yb-based RSBT can increase PTV D_90% or decrease urethral D_10% relative to HDR-BT with treatment times of less than 1 h using a single-source robotic delivery system with treatment delivered in a single fraction. The prototype helical delivery system was able to demonstrate adequate mechanical accuracy.

Monte Carlo dosimetric characterization of a new high dose rate 169 Yb brachytherapy source and independent verification using a multipoint plastic scintillator detector

PURPOSE

A prototype 169 Yb source was developed in combination with a dynamic rotating platinum shield system (AIM-Brachy) to deliver intensity modulated brachytherapy (IMBT). The purpose of this study was to evaluate the dosimetric characteristics of the bare/shielded 169 Yb source using Monte Carlo (MC) simulations and perform an independent dose verification using a dosimetry platform based on a multipoint plastic scintillator detector (mPSD).

METHODS

The TG-43U1 dosimetric parameters were calculated for the source model using RapidBrachyMCTPS. Real-time dose rate measurements were performed in a water tank for both the bare/shielded source using a custom remote afterloader. For each dwell position, the dose rate was independently measured by the three scintillators (BCF-10, BCF-12, and BCF-60). For the bare source, dose rate was measured at distances up to 3 cm away from the source over a range of 7 cm along the catheter. For the shielded source, measurements were performed with the mPSD placed at 1 cm from the source at four different azimuthal angles ( 0 ∘ , 9 0 ∘ , 18 0 ∘ , and 27 0 ∘ ).

RESULTS

The dosimetric parameters were tabulated for the source model. For the bare source, differences between measured and calculated along-away dose rates were generally below 5-10%. Along the transverse axis, deviations were, on average (range), 3.3% (0.6-6.2%) for BCF-10, 1.7% (0.9-2.9%) for BCF-12, and 2.2% (0.3-4.4%) for BCF-60. The maximum dose rate reduction due to shielding at a radial distance of 1 cm was 88.8 ± 1.2%, compared to 83.5 ± 0.5% as calculated by MC.

CONCLUSIONS

The dose distribution for the bare/shielded 169 Yb source was independently verified using mPSD with good agreement in regions close to the source. The 169 Yb source coupled with the partial-shielding system is an effective technique to deliver IMBT.

Needle-free cervical cancer treatment using helical multishield intracavitary rotating shield brachytherapy with the 169 Yb Isotope

PURPOSE

To assess the capability of an intracavitary ¹⁶⁹ Yb-based helical multishield rotating shield brachytherapy (RSBT) delivery system to treat cervical cancer. The proposed RSBT delivery system contains a pair of 1.25 mm thick platinum partial shields with 45° and 180° emission angles, which travel in a helical pattern within the applicator.

METHODS

A helically threaded tandem applicator with a 45° tandem curvature containing a helically threaded catheter was designed. A 0.6 mm diameter ¹⁶⁹ Yb source with a length of 10.5 mm was simulated. A 37-patient treatment planning study, based on Monte Carlo dose calculations using MCNP5, was conducted with high-risk clinical target volumes (HR-CTVs) of 41.2-192.8 cm³ (average ± standard deviation of 79.9 ± 35.8 cm³ ). All patients were assumed to receive 25 fractions of 1.8 Gy of external beam radiation therapy (EBRT) before receiving 5 fractions of high-dose-rate brachytherapy (HDR-BT). For each patient, ¹⁹² Ir-based intracavitary (IC) HDR-BT, ¹⁹² Ir-based intracavitary/interstitial (IC/IS) HDR-BT using a hybrid applicator with eight IS needles, and ¹⁶⁹ Yb-based RSBT plans were generated.

RESULTS

For the IC, IC/IS, and RSBT treatment plans, 38%, 84%, and 86% of the plans, respectively, met the planning goal of an HR-CTV D₉₀ (minimum dose to hottest 90%) of 85 Gy_EQD2 (α/β = 10 Gy). Median (25th percentile, 75th percentile) treatment times for IC, IC/IS, and RSBT were 11.71 (6.62, 15.40) min, 68.00 (45.02, 80.02) min, and 25.30 (13.87, 35.39) min, respectively. ¹⁹² Ir activities ranging from 159.1-370 GBq (4.3-10 Ci) and ¹⁶⁹ Yb activities ranging from 429.2-999 GBq (11.6-27 Ci) were used, which correspond to the same clinical ranges of dose rates at 1 cm off-source-axis in water. Extra needle insertion and planning time beyond that needed for intracavitary-only approaches was accounted for in the IC/IS treatment time calculations.

CONCLUSION

¹⁶⁹ Yb-based RSBT for cervical cancer met the HR-CTV D₉₀ goal of 85 Gy in a greater percentage of the patients considered than IC/IS (86% vs 84%, respectively) and can reduce overall treatment time relative to IC/IS. ¹⁶⁹ Yb-based RSBT could be used to replace IC/IS in instances where IC/IS treatment is not available, especially in instances when HR-CTV volumes are ≥30 cm³ .

Maximum RBE change in 192Ir, 125I, and 169Yb brachytherapy and the corresponding effect on treatment planning

PURPOSE

The purpose of this study was to examine RBE variation as a function of distance from the radioactive source, and the potential impact of this variation on a realistic prostate brachytherapy treatment plan.

METHODS

Three brachytherapy sources (¹²⁵I, ¹⁹²Ir, and ¹⁶⁹Yb) were modelled in Geant4 Monte Carlo code, and the resulting electron energy spectrum in water in 3D space around these sources was scored (voxel size of 2 mm³). With this energy spectrum, microdosimetric techniques were used to calculate the maximum RBE, RBE_M, as a function of distance from the source. RBE_M of ¹²⁵I relative to ¹⁹²Ir was calculated in order to validate simulations against literature; all other RBE_M calculations were done by normalizing electron fluence at various distances to the source position. In order to examine the impact of RBE_M variation in treatment planning, a realistic ¹⁹²Ir prostate plan was re-evaluated in terms of RBE instead of absorbed dose.

RESULTS

The RBE_M of ¹²⁵I, ¹⁹²Ir, and ¹⁶⁹Yb at 8 cm away from the source was 0.994 (+/-0.002), 1.030 (+/-0.003), and 1.066 (+/-0.008), respectively. RBE_M in the HDR prostate treatment plan exhibited several hot (+3.6% in RBE_M) spots.

CONCLUSIONS

The large increase RBE_M observed in ¹⁶⁹Yb has not yet been described in the literature. Despite the presence of radiobiological hotspots in the HDR treatment, these variations are likely nominal and clinically insignificant.

Efficient 169 Yb high-dose-rate brachytherapy source production using reactivation

PURPOSE

To present and quantify the effectiveness of a method for the efficient production of ¹⁶⁹ Yb high-dose-rate brachytherapy sources with 27 Ci activity upon clinical delivery, which have about the same dose rate in water at 1 cm from the source center as 10 Ci ¹⁹² Ir sources.

MATERIALS

A theoretical framework for ¹⁶⁹ Yb source activation and reactivation using thermal neutrons in a research reactor and ¹⁶⁸ Yb-Yb₂ O₃ precursor is derived and benchmarked against published data. The model is dependent primarily on precursor ¹⁶⁸ Yb enrichment percentage, active source volume of the active element, and average thermal neutron flux within the active source.

RESULTS

Efficiency gains in ¹⁶⁹ Yb source production are achievable through reactivation, and the gains increase with active source volume. For an average thermal neutron flux within the active source of 1 × 10¹⁴  n cm^-2  s^-1 , increasing the active source volume from 1 to 3 mm³ decreased reactor-days needed to generate one clinic-year of ¹⁶⁹ Yb from 256 days yr^-1 to 59 days yr^-1 , and 82%-enriched precursor dropped from 80 mg yr^-1 to 21 mg yr^-1 . A resource reduction of 74%-77% is predicted for an active source volume increase from 1 to 3 mm³ .

CONCLUSIONS

Dramatic cost savings are achievable in ¹⁶⁹ Yb source production costs through reactivation if active sources larger than 1 mm³ are used.

Evaluation of dose point kernel rescaling methods for nanoscale dose estimation around gold nanoparticles using Geant4 Monte Carlo simulations

The absence of proper nanoscale experimental techniques to investigate the dose-enhancing properties of gold nanoparticles (GNPs) interacting with radiation has prompted the development of various Monte Carlo (MC)-based nanodosimetry techniques that generally require considerable computational knowledge, time and specific tools/platforms. Thus, this study investigated a hybrid computational framework, based on the electron dose point kernel (DPK) method, by combining Geant4 MC simulations with an analytical approach. This hybrid framework was applied to estimate the dose distributions around GNPs due to the secondary electrons emitted from GNPs irradiated by various photon sources. Specifically, the equivalent path length approximation was used to rescale the homogeneous DPKs for heterogeneous GNPs embedded in water/tissue. Compared with Geant4 simulations, the hybrid framework halved calculation time while utilizing fewer computer resources, and also resulted in mean discrepancies less than 20 and 5% for Yb-169 and 6 MV photon irradiation, respectively. Its appropriateness and computational efficiency in handling more complex cases were also demonstrated using an example derived from a transmission electron microscopy image of a cancer cell containing internalized GNPs. Overall, the currently proposed hybrid computational framework can be a practical alternative to full-fledged MC simulations, benefiting a wide range of GNP- and radiation-related applications.

Radiosensitization of Prostate Cancers In Vitro and In Vivo to Erbium-filtered Orthovoltage X-rays Using Actively Targeted Gold Nanoparticles

Theoretical investigations suggest that gold nanoparticle (GNP)-mediated radiation dose enhancement and radiosensitization can be maximized when photons interact with gold, predominantly via photoelectric absorption. This makes ytterbium (Yb)-169, which emits photons with an average energy of 93 keV (just above the K-edge of gold), an ideal radioisotope for such purposes. This investigation tests the feasibility of tumor-specific prostate brachytherapy achievable with Yb-169 and actively targeted GNPs, using an external beam surrogate of Yb-169 created from an exotic filter material - erbium (Er) and a standard copper-filtered 250 kVp beam. The current in vitro study shows that treatment of prostate cancer cells with goserelin-conjugated gold nanorods (gGNRs) promotes gonadotropin releasing hormone receptor-mediated internalization and enhances radiosensitivity to both Er-filtered and standard 250 kVp beams, 14 and 10%, respectively. While the degree of GNP-mediated radiosensitization as seen from the in vitro study may be considered moderate, the current in vivo study shows that gGNR treatment plus Er-filtered x-ray irradiation is considerably more effective than radiation treatment alone (p < 0.0005), resulting in a striking reduction in tumor volume (50% smaller) 2 months following treatment. Overall, the current results provide strong evidence for the feasibility of tumor-specific prostate brachytherapy with Yb-169 and gGNRs.