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

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.

The population percentile allowance method for determining systematic spatial error tolerances for temporary intensity modulated brachytherapy

BACKGROUND

Multiple approaches are under development for delivering temporary intensity modulated brachytherapy (IMBT) using partially shielded applicators wherein the delivered dose distributions are sensitive to spatial uncertainties in both the applicator position and shield orientation, rather than only applicator position as with conventional high-dose-rate brachytherapy (HDR-BT). Sensitivity analyses to spatial uncertainties have been reported as components of publications on these emerging technologies, however, a generalized framework for the rigorous determination of the spatial uncertainty tolerances of dose-volume parameters is needed.

PURPOSE

To derive and present the population percentile allowance (PPA) method, a generalized mathematical and statistical framework to evaluate the tolerance of temporary IMBT approaches to spatial uncertainties in applicator position and shield orientation.

METHODS

A mathematical formalism describing geometric applicator position and shield orientation shifts was derived that supports straight and curved applicators and applies to serial and helical rotating shield brachytherapy (RSBT) and direction modulated brachytherapy (DMBT). The PPA method entails defining the percentage of a patient population receiving a given therapy that is, allowed to receive dose-volume errors in the target volume and specified organs at risk of a defined percentage or less, then determining what combinations of applicator position and shield orientation systematic errors would be expected to produce that outcome in the population. The PPA method was applied to the use case of multi-shield helical 169 Yb-based RSBT for cervical cancer, with 45° and 180° shield emission angles. A total of 37 cervical cancer patients were considered in the population, with average (± 1 standard deviation) HR-CTV volumes of 79 cm3  ± 37 cm3 and optimized baseline treatment plans (no spatial uncertainties applied) created for each patient to meet dose-volume requirements of 85 GyEQD2 (equivalent uniform dose in 2 Gy fraction), with D2cc tolerance doses of 90 GyEQD2 , 75 GyEQD2 , and 75 GyEQD2 for bladder, rectum, and sigmoid colon, respectively.

RESULTS

For the PPA requirement that 90% of cervical cancer patients receiving multi-shield helical RSBT could have a maximum dose-volume uncertainty of 10% for high-risk clinical target volume (HR-CTV) D90 (minimum dose to hottest 90%) and bladder, rectum, and sigmoid colon D2cc (minimum dose to hottest 2 cm3 ), the tolerance systematic applicator position and shield orientation uncertainties were approximately ± 1.0 mm and ± 4.25°, respectively. For ± 1.5 mm and ± 5° systematic applicator position and shield orientation tolerances, 90% of the patients considered would have a maximum dose-volume uncertainty of 12.8% or less.

CONCLUSION

The PPA method was formalized to determine the temporary IMBT spatial uncertainty tolerances that would be expected to result in an allowed percentage of a population of patients receiving relative dose-volume errors above a defined percentage. Multi-shield, helical 169 Yb-based RSBT for cervical cancer was evaluated and tolerances determined, which, if applied on each treatment fraction, would represent an extreme situation. The PPA method is applicable to a variety of temporary IMBT approaches and can be used to rigorously determine the design parameters for the delivery systems such as mechanical driver motor accuracy, shield angle backlash, applicator rotation, and applicator fixation stability.

Learning from the past: a century of accuracy, aspirations, and aspersions in brachytherapy

The oldest form of radiation therapy, brachytherapy, has been investigated and reported in the scientific and medical literature for well over a century. Known by many names over the years, radium-based, empirical practices evolved over decades to contemporary practice. This includes treatment at various dose rates using multiple radionuclides or even electrically generated photon sources. Predictions or prognostications of what may happen in the future enjoy a history that spans centuries, e.g. those by Nostradamus in the 1500s. In this review article, publications from several eras of past practice between the early 1900s and the late 2010s where the authors address the "future of brachytherapy" are presented, and for many of these publications, one can use the benefit of the intervening years to comment on the accuracy or the inaccuracies inherent in those publications. Finally, recently published papers are reviewed to examine current expectations for the future practice of brachytherapy.

Emerging technologies in brachytherapy

Brachytherapy is a mature treatment modality. The literature is abundant in terms of review articles and comprehensive books on the latest established as well as evolving clinical practices. The intent of this article is to part ways and look beyond the current state-of-the-art and review emerging technologies that are noteworthy and perhaps may drive the future innovations in the field. There are plenty of candidate topics that deserve a deeper look, of course, but with practical limits in this communicative platform, we explore four topics that perhaps is worthwhile to review in detail at this time. First, intensity modulated brachytherapy (IMBT) is reviewed. The IMBT takes advantage ofanisotropicradiation profile generated through intelligent high-density shielding designs incorporated onto sources and applicators such to achieve high quality plans. Second, emerging applications of 3D printing (i.e. additive manufacturing) in brachytherapy are reviewed. With the advent of 3D printing, interest in this technology in brachytherapy has been immense and translation swift due to their potential to tailor applicators and treatments customizable to each individual patient. This is followed by, in third, innovations in treatment planning concerning catheter placement and dwell times where new modelling approaches, solution algorithms, and technological advances are reviewed. And, fourth and lastly, applications of a new machine learning technique, called deep learning, which has the potential to improve and automate all aspects of brachytherapy workflow, are reviewed. We do not expect that all ideas and innovations reviewed in this article will ultimately reach clinic but, nonetheless, this review provides a decent glimpse of what is to come. It would be exciting to monitor as IMBT, 3D printing, novel optimization algorithms, and deep learning technologies evolve over time and translate into pilot testing and sensibly phased clinical trials, and ultimately make a difference for cancer patients. Today's fancy is tomorrow's reality. The future is bright for brachytherapy.

Optimization in treatment planning of high dose-rate brachytherapy - Review and analysis of mathematical models

Treatment planning in high dose-rate brachytherapy has traditionally been conducted with manual forward planning, but inverse planning is today increasingly used in clinical practice. There is a large variety of proposed optimization models and algorithms to model and solve the treatment planning problem. Two major parts of inverse treatment planning for which mathematical optimization can be used are the decisions about catheter placement and dwell time distributions. Both these problems as well as integrated approaches are included in this review. The proposed models include linear penalty models, dose-volume models, mean-tail dose models, quadratic penalty models, radiobiological models, and multiobjective models. The aim of this survey is twofold: (i) to give a broad overview over mathematical optimization models used for treatment planning of brachytherapy and (ii) to provide mathematical analyses and comparisons between models. New technologies for brachytherapy treatments and methods for treatment planning are also discussed. Of particular interest for future research is a thorough comparison between optimization models and algorithms on the same dataset, and clinical validation of proposed optimization approaches with respect to patient outcome.

On the impact of absorbed dose specification, tissue heterogeneities, and applicator heterogeneities on Monte Carlo-based dosimetry of Ir-192, Se-75, and Yb-169 in conventional and intensity-modulated brachytherapy for the treatment of cervical cancer

PURPOSE

The purpose of this study was to evaluate the impact of dose reporting schemes and tissue/applicator heterogeneities for ¹⁹² Ir-, ⁷⁵ Se-, and ¹⁶⁹ Yb-based MRI-guided conventional and intensity-modulated brachytherapy.

METHODS AND MATERIALS

Treatment plans using a variety of dose reporting and tissue/applicator segmentation schemes were generated for a cohort (n = 10) of cervical cancer patients treated with ¹⁹² Ir-based Venezia brachytherapy. Dose calculations were performed using RapidBrachyMCTPS, a Geant4-based research Monte Carlo treatment planning system. Ultimately, five dose calculation scenarios were evaluated: (a) dose to water in water (D_w,w ); (b) D_w,w taking the applicator material into consideration (D_w,wApp ); (c) dose to water in medium (D_w,m ); (d and e) dose to medium in medium with mass densities assigned either nominally per structure (D_{m,m (Nom)} ) or voxel-by-voxel (D_m,m ).

RESULTS

Ignoring the plastic Venezia applicator (D_w,wApp ) overestimates D_m,m by up to 1% (average) with high energy source (¹⁹² Ir and ⁷⁵ Se) and up to 2% with ¹⁶⁹ Yb. Scoring dose to water (D_w,wApp or D_w,m ) generally overestimates dose and this effect increases with decreasing photon energy. Reporting dose other than D_m,m (or D_{m,m Nom} ) for ¹⁶⁹ Yb-based conventional and intensity-modulated brachytherapy leads to a simultaneous overestimation (up to 4%) of CTV_HR D₉₀ and underestimation (up to 2%) of bladder D_2cc due to a significant dip in the mass-energy absorption ratios at the depths of nearby targets and OARs. Using a nominal mass-density assignment per structure, rather than a CT-derived voxel-by-voxel assignment for MRI-guided brachytherapy, amounts to a dose error up to 1% for all radionuclides considered.

CONCLUSIONS

The effects of the considered dose reporting schemes trend correspondingly between conventional and intensity-modulated brachytherapy. In the absence of CT-derived mass densities, MRI-only-based dosimetry can adequately approximate D_m,m by assigning nominal mass densities to structures. Tissue and applicator heterogeneities do not significantly impact dosimetry for ¹⁹² Ir and ⁷⁵ Se, but do for ¹⁶⁹ Yb; dose reporting must be explicitly defined since D_w,m and D_w,w may overstate the dosimetric benefits.

A novel minimally invasive dynamic-shield, intensity-modulated brachytherapy system for the treatment of cervical cancer

PURPOSE

To present a novel, MRI-compatible dynamicshield intensity modulated brachytherapy (IMBT) applicator and delivery system using ¹⁹² Ir, ⁷⁵ Se, and ¹⁶⁹ Yb radioisotopes for the treatment of locally advanced cervical cancer. Needle-free IMBT is a promising technique for improving target coverage and organs at risk (OAR) sparing.

METHODS AND MATERIALS

The IMBT delivery system dynamically controls the rotation of a novel tungsten shield placed inside an MRI-compatible, 6-mm wide intrauterine tandem. Using 36 cervical cancer cases, conventional intracavitary brachytherapy (IC-BT) and intracavitary/interstitial brachytherapy (IC/IS-BT) (10Ci ¹⁹² Ir) plans were compared to IMBT (10Ci ¹⁹² Ir; 11.5Ci ⁷⁵ Se; 44Ci ¹⁶⁹ Yb). All plans were generated using the Geant4-based Monte Carlo dose calculation engine, RapidBrachyMC. Treatment plans were optimized then normalized to the same high-risk clinical target volume (HR-CTV) D₉₀ and the D_2cc for bladder, rectum, and sigmoid in the research brachytherapy planning system, RapidBrachyMCTPS. Plans were renormalized until either of the three OAR reached dose limits to calculate the maximum achievable HR-CTV D₉₀ and D₉₈ .

RESULTS

Compared to IC-BT, IMBT with either of the three radionuclides significantly improves the HR-CTV D₉₀ and D₉₈ by up to 5.2% ± 0.3% (P < 0.001) and 6.7% ± 0.5% (P < 0.001), respectively, with the largest dosimetric enhancement when using ¹⁶⁹ Yb followed by ⁷⁵ Se and then ¹⁹² Ir. Similarly, D_2cc for all OAR improved with IMBT by up to 7.7% ± 0.6% (P < 0.001). For IC/IS-BT cases, needle-free IMBT achieved clinically acceptable plans with ¹⁶⁹ Yb-based IMBT further improving HR-CTV D₉₈ by 1.5% ± 0.2% (P = 0.034) and decreasing sigmoid D_2cc by 1.9% ± 0.4% (P = 0.048). Delivery times for IMBT are increased by a factor of 1.7, 3.3, and 2.3 for ¹⁹² Ir, ⁷⁵ Se, and ¹⁶⁹ Yb, respectively, relative to conventional ¹⁹² Ir BT.

CONCLUSIONS

Dynamic shield IMBT provides a promising alternative to conventional IC- and IC/IS-BT techniques with significant dosimetric enhancements and even greater improvements with intermediate energy radionuclides. The ability to deliver a highly conformal, OAR-sparing dose without IS needles provides a simplified method for improving the therapeutic ratio less invasively and in a less resource intensive manner.

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.