Background: The steep radiation dose gradients in cervical cancer brachytherapy (BT) necessitate a thorough understanding of the behavior of afterloader source cables or needles in the curved channels of (patient-tailored) applicators. Purpose: The purpose of this study is to develop and validate computer models to simulate: (1) BT source positions, and (2) insertion forces of needles in curved applicator channels. The methodology presented can be used to improve the knowledge of instrument behavior in current applicators and aid the development of novel (3D-printed) BT applicators. Methods: For the computer models, BT instruments were discretized in finite elements. Simulations were performed in SPACAR by formulating nodal contact force and motion input models and specifying the instruments’ kinematic and dynamic properties. To evaluate the source cable model, simulated source paths in ring applicators were compared with manufacturer-measured source paths. The impact of discrepancies on the dosimetry was estimated for standard plans. To validate needle models, simulated needle insertion forces in curved channels with varying curvature, torsion, and clearance, were compared with force measurements in dedicated 3D-printed templates. Results: Comparison of simulated with manufacturer-measured source positions showed 0.5–1.2 mm median and <2.0 mm maximum differences, in all but one applicator geometry. The resulting maximum relative dose differences at the lateral surface and at 5 mm depth were 5.5% and 4.7%, respectively. Simulated insertion forces for BT needles in curved channels accurately resembled the forces experimentally obtained by including experimental uncertainties in the simulation. Conclusion: The models developed can accurately predict source positions and insertion forces in BT applicators. Insights from these models can aid novel applicator design with improved motion and force transmission of BT instruments, and contribute to the estimation of overall treatment precision. The methodology presented can be extended to study other applicator geometries, flexible instruments, and afterloading systems.

Prostate cancer patients with an enlarged prostate and/or excessive pubic arch interference (PAI) are generally considered non-eligible for high-dose-rate (HDR) brachytherapy (BT). Steerable needles have been developed to make these patients eligible again. This study aims to validate the dosimetric impact and performance of steerable needles within the conventional clinical setting. HDR BT treatment plans were generated, needle implantations were performed in a prostate phantom, with prostate volume > 55 cm³ and excessive PAI of 10 mm, and pre- and post-implant dosimetry were compared considering the dosimetric constraints: prostate V₁₀₀ > 95 % (13.50 Gy), urethra D_0.1cm³ < 115 % (15.53 Gy) and rectum D_1cm³ < 75 % (10.13 Gy). The inclusion of steerable needles resulted in a notable enhancement of the dose distribution and prostate V₁₀₀ compared to treatment plans exclusively employing rigid needles to address PAI. Furthermore, the steerable needle plan demonstrated better agreement between pre- and post-implant dosimetry (prostate V₁₀₀: 96.24 % vs. 93.74 %) compared to the rigid needle plans (79.13 % vs. 72.86 % and 87.70 % vs. 81.76 %), with no major changes in the clinical workflow and no changes in the clinical set-up. The steerable needle approach allows for more flexibility in needle positioning, ensuring a highly conformal dose distribution, and hence, HDR BT is a feasible treatment option again for prostate cancer patients with an enlarged prostate and/or excessive PAI.

Purpose. Patient-tailored intracavitary/interstitial (IC/IS) brachytherapy (BT) applicators may increase dose conformity in cervical cancer patients. Current configuration planning methods in these custom applicators rely on manual specification or a small set of (straight) needles. This work introduces and validates a two-stage approach for establishing channel configurations in the 3D printed patient-tailored ARCHITECT applicator. Methods. For each patient, the patient-tailored applicator shape was based on the first BT application with a commercial applicator and integrated connectors to a commercial (Geneva) intrauterine tube and two lunar ring channels. First, a large candidate set was generated of channels that steer the needle to desired poses in the target region and are contained in the applicator. The channels’ centrelines were represented by Bézier curves. Channels running between straight target segments and entry points were optimised and refined to ensure (dynamic) feasibility. Second, channel configurations were selected using geometric coverage optimisation. This workflow was applied to establish patient-tailored geometries for twenty-two patients previously treated using the Venezia applicator. Treatment plans were automatically generated using the in-house developed algorithm BiCycle. Plans for the clinically used configuration, T P clin , and patient-tailored configuration, T P arch , were compared. Results. Channel configurations could be generated in clinically feasible time (median: 2651 s, range 1826-3812 s). All T P arch and T P clin plans were acceptable, but planning aims were more frequently attained with patient-tailored configurations (115/132 versus 100/132 instances). Median CTV_IR D 98 and bladder D 2 c m 3 doses significantly improved ( p < 0.001 and p < 0.01 respectively) in T P arch plans in comparison with T P clin plans, and in approximately half of the patients dosimetric indices improved. Conclusion. Automated patient-tailored BT channel configuration planning for 3D printed applicators is clinically feasible. A treatment planning study showed that all plans met planning limits for the patient-tailored configurations, and in selected cases improved the plan quality in comparison with commercial applicator configurations.

The combination of interstitial hyperthermia treatment (IHT) with high dose rate brachytherapy (HDR‐BT) can improve clinical outcomes since it highly enhances the efficiency of cell kill, especially when applied simultaneously. Therefore, we have developed the ThermoBrachy applicators. To effectively apply optimal targeted IHT, treatment planning is considered essential. However, treatment planning in IHT is rarely applied as it is regarded as difficult to accurately calculate the deposited energy in the tissue in a short enough time for clinical practice. In this study, we investigated various time‐efficient methods for fast computation of the electromagnetic (EM) energy deposition resulting from the ThermoBrachy applicators. Initially, we investigated the use of an electro‐quasistatic solver. Next, we extended our investigation to the application of geometric simplifications. Furthermore, we investigated the validity of the superpositioning principle, which can enable adaptive treatment plan optimization without the need for continuous recomputation of the EM field. Finally, we evaluated the accuracy of the methods by comparing them to the golden standard Finite‐Difference Time‐Domain calculation method using gamma‐index analysis. The simplifications considerably reduced the computation time needed, improving from >12 h to a few seconds. All investigated methods showed excellent agreement with the golden standard by showing a >99% passing rate with 1%/0.5 mm Dose Difference and Distance‐to‐Agreement criteria. These results allow the proposed electromagnetic simulation method to be used for fast and accurate adaptive treatment planning.

In high-dose-rate brachytherapy (HDR-BT) for prostate cancer treatment, interstitial hyperthermia (IHT) is applied to sensitize the tumor to the radiation (RT) dose, aiming at a more efficient treatment. Simultaneous application of HDR-BT and IHT is anticipated to provide maximum radiosensitization of the tumor. With this rationale, the ThermoBrachyTherapy applicators have been designed and developed, enabling simultaneous irradiation and heating. In this research, we present a method to optimize the three-dimensional temperature distribution for simultaneous HDR-BT and IHT based on the resulting equivalent physical dose (EQDphys) of the combined treatment. First, the temperature resulting from each electrode is precomputed. Then, for a given set of electrode settings and a precomputed radiation dose, the EQDphys is calculated based on the temperature-dependent linear-quadratic model. Finally, the optimum set of electrode settings is found through an optimization algorithm. The method is applied on implant geometries and anatomical data of 10 previously irradiated patients, using reported thermoradiobiological parameters and physical doses. We found that an equal equivalent dose coverage of the target can be achieved with a physical RT dose reduction of 20% together with a significantly lower EQDphys to the organs at risk (p-value < 0.001), even in the least favorable scenarios. As a result, simultaneous ThermoBrachy-Therapy could lead to a relevant therapeutic benefit for patients with prostate cancer.

Brachytherapy is a common treatment in cervical, uterine and vaginal cancer management. The technique is characterised by rapid developments in the fields of medical imaging, dosimetry planning and personalised medical device design. To reduce unnecessary burden on patients, assessments and training of these technologies should preferable be done using high-fidelity physical phantoms. In this study, anthropomorphic deformable phantoms of the vaginal wall and cavity were developed for image-guided adaptive brachytherapy, in which vaginal wall biomechanics were mimicked. Phantoms were produced from both silicone and polyvinyl alcohol materials. Material characterisations were performed with uniaxial tensile tests, via which Young’s moduli and toughness were quantified. In addition, the contrast between adjacent phantom layers was quantified in magnetic resonance images. The results showed that stress-strain curves of the silicone phantoms were within the range of those found in healthy human vaginal wall tissues. Sample preconditioning had a large effect on Young’s moduli, which ranged between 2.13 and 6.94 MPa in silicone. Toughness was a more robust and accurate metric for biomechanical matching, and ranged between 0.23 and 0.28 ·10⁶ J·m^-3 as a result of preconditioning. The polyvinyl alcohol phantoms were not stiff or tough enough, with a Young’s modulus of 0.16 MPa and toughness of 0.02 ·10⁶ J·m^-3. All materials used could be clearly delineated in magnetic resonance images, although the MRI sequence did affect layer contrast. In conclusion, we developed anthropomorphic deformable phantoms that mimic vaginal wall tissue and are well visible in magnetic resonance images. These phantoms will be used to evaluate the properties and to optimise the development and use of personalised brachytherapy applicators.

Purpose: To develop and evaluate a fast, automated multi-criterial treatment planning approach for adaptive high-dose-rate (HDR) intracavitary + interstitial brachytherapy (BT) for locally advanced cervical cancer. Methods and materials: Twenty-two previously delivered single fraction MRI-based HDR treatment plans (SF_clin) were used to guide training of our in-house system for multi-criterial autoplanning, aiming for an autoplan quality superior to the training plans, while respecting the clinically desired “pear-shaped” dose distribution. Next, the configured algorithm was used to automatically generate treatment plans for 63 other fractions (SF_auto). The SF_auto plans were compared to the corresponding SF_clin plans in blind pairwise comparisons by an expert clinician. Then, the effect of adaptive autoplanning on total treatment (TT) plans (external beam + 3 BT fractions) was evaluated for 16 patients by simulating the clinically applied adaptive strategy to generate TT_auto plans and compare them with the corresponding clinical treatments (TT_clin). Results: In the blind comparisons, all SF_auto plans were considered clinically acceptable. In 62/63 comparisons, SF_auto plans were considered at least as good as, or better than the corresponding SF_clin. The average optimization time for autoplanning was 20.5 ± 19.2 s (range 4.4–106.4 s) per plan. In 14 of 16 TT_auto plans, the desired total dose of 90 Gy (EQD₂) was obtained, compared to only 9 in the corresponding TT_clin, while autoplanning also decreased bladder and rectum doses. Conclusions: Fast, fully-automated multi-criterial treatment planning for adaptive HDR-BT for locally advanced cervical cancer is feasible. Autoplans were superior to corresponding clinical plans.

We developed a fast and fully-automated, multi-criteria treatment planning workflow for high dose rate brachytherapy (HDR-BT). In this workflow, the patient-CT with catheter reconstructions and dwell positions are imported from the clinical TPS into a novel system for automated dwell time optimisation. The optimised dwell times are then imported into the clinical TPS. The aims of automation were (1) planner-independent, enhanced plan quality, (2) short optimisation times. Our in-house developed system for fully automated, multi-criteria external beam radiotherapy (EBRT) treatment planning (Erasmus-iCycle) was adapted for optimisation of HDR-BT dose distributions. The investigations were performed with planning CT scans with catheter reconstructions and delineations of twenty-five low- A nd intermediate-risk prostate cancer patients who were previously treated in our center with 4 × 9.5 Gy HDR-BT. Automatically generated plans (autoplans) were compared to the corresponding clinical plans. All evaluations were performed in the clinical TPS. The requested 95% tumour coverage was obtained for all autoplans, while this was only observed in 23/25 clinical plans. All autoplans showed a consistent reduction of the D1% for the highest prioritised OAR, the urethra. The average and maximum reductions were 6.3%-point and 12.1%-point of the prescribed dose, respectively. In addition, conformality of the autoplans was higher. The autoplans had slightly smaller delivery times. Autoplanning took on average 4.6 s, including computation of the dose kernels.