Dose and DNA damage modelling of diffusing alpha-emitters radiation therapy using Geant4.

PURPOSE: Diffusing alpha-emitters radiation therapy (DaRT) is a brachytherapy technique using α-particles to treat solid tumours. The high linear energy transfer (LET) and short range of α-particles make them good candidates for the targeted treatment of cancer. Treatment planning of DaRT requires a good understanding of the dose from α-particles and the other particles released in the 224Ra decay chain. METHODS: The Geant4 Monte Carlo toolkit has been used to simulate a DaRT seed to better understand the dose contribution from all particles and simulate the DNA damage due to this treatment. RESULTS: Close to the seed α-particles deliver the majority of dose, however at radial distances greater than 4 mm, the contribution of ß-particles is greater. The RBE has been estimated as a function of number of double strand breaks (DSBs) and complex DSBs. A maximum seed spacing of 5.5 mm and 6.5 mm was found to deliver at least 20 Gy RBE weighted dose between the seeds for RBEDSB and RBEcDSB respectively. CONCLUSIONS: The DNA damage changes with radial distance from the seed and has been found to become less complex with distance, which is potentially easier for the cell to repair. Close to the seed α-particles contribute the majority of dose, however the contribution from other particles cannot be neglected and may influence the choice of seed spacing.

First in vitro measurement of VHEE relative biological effectiveness (RBE) in lung and prostate cancer cells using the ARES linac at DESY.

Very high energy electrons (VHEE) are a potential candidate for radiotherapy applications. This includes tumours in inhomogeneous regions such as lung and prostate cancers, due to the insensitivity of VHEE to inhomogeneities. This study explores how electrons in the VHEE range can be used to perform successful in vitro radiobiological studies. The ARES (accelerator research experiment at SINBAD) facility at DESY, Hamburg, Germany was used to deliver 154 MeV electrons to both prostate (PC3) and lung (A549) cancer cells in suspension. Dose was delivered to samples with repeatability and uniformity, quantified with Gafchromic film. Cell survival in response to VHEE was measured using the clonogenic assay to determine the biological effectiveness of VHEE in cancer cells for the first time using this method. Equivalent experiments were performed using 300 kVp X-rays, to enable VHEE irradiated cells to be compared with conventional photons. VHEE irradiated cancer cell survival was fitted to the linear quadratic (LQ) model (R2 = 0.96-0.97). The damage from VHEE and X-ray irradiated cells at doses between 1.41 and 6.33 Gy are comparable, suggesting similar relative biological effectiveness (RBE) between the two modalities. This suggests VHEE is as damaging as photon radiotherapy and therefore could be used to successfully damage cancer cells during radiotherapy. The RBE of VHEE was quantified as the relative doses required for 50% (D0.5) and 10% (D0.1) cell survival. Using these values, VHEE RBE was measured as 0.93 (D0.5) and 0.99 (D0.1) for A549 and 0.74 (D0.5) and 0.93 (D0.1) for PC3 cell lines respectively. For the first time, this study has shown that 154 MeV electrons can be used to effectively kill lung and prostate cancer cells, suggesting that VHEE would be a viable radiotherapy modality. Several studies have shown that VHEE has characteristics that would offer significant improvements over conventional photon radiotherapy for example, electrons are relatively easy to steer and can be used to deliver dose rapidly and with high efficiency. Studies have shown improved dose distribution with VHEE in treatment plans, in comparison to VMAT, indicating that VHEE can offer improved and safer treatment plans with reduced side effects. The biological response of cancer cells to VHEE has not been sufficiently studied as of yet, however this initial study provides some initial insights into cell damage. VHEE offers significant benefits over photon radiotherapy and therefore more studies are required to fully understand the biological effectiveness of VHEE.

Estimation of the relative biological effectiveness for double strand break induction of clinical kilovoltage beams using Monte Carlo simulations.

BACKGROUND: The Relative Biological Effectiveness (RBE) of kilovoltage photon beams has been previously investigated in vitro and in silico using analytical methods. The estimated values range from 1.03 to 1.82 depending on the methodology and beam energies examined. PURPOSE: The focus of this work was to independently estimate RBE values for a range of clinically used kilovoltage beams (70-200 kVp) while investigating the suitability of using TOPAS-nBio for this task. METHODS: Previously validated spectra of clinical beams were used to generate secondary electron spectra at several depths in a water tank phantom via TOPAS Monte Carlo (MC) simulations. Cell geometry was irradiated with the secondary electrons in TOPAS-nBio MC simulations. The deposited dose and the calculated number of DNA strand breaks were used to estimate RBE values. RESULTS: Monoenergetic secondary electron simulations revealed the highest direct and indirect double strand break yield at approximately 20 keV. The average RBE value for the kilovoltage beams was calculated to be 1.14. CONCLUSIONS: TOPAS-nBio was successfully used to estimate the RBE values for a range of clinical radiotherapy beams. The calculated value was in agreement with previous estimates, providing confidence in its clinical use in the future.

A method to predict space radiation biological effectiveness for non-cancer effects following intense Solar Particle Events.

In addition to the continuous exposure to cosmic rays, astronauts in space are occasionally exposed to Solar Particle Events (SPE), which involve less energetic particles but can deliver much higher doses. The latter can exceed several Gy in a few hours for the most intense SPEs, for which non-stochastic effects are thus a major concern. To identify adequate shielding conditions that would allow respecting the dose limits established by the various space agencies, the absorbed dose in the considered organ/tissue must be multiplied by the corresponding Relative Biological Effectiveness (RBE), which is a complex quantity depending on several factors including particle type and energy, considered biological effect, level of effect (and thus absorbed dose), etc. While in several studies only the particle-type dependence of RBE is taken into account, in this work we developed and applied a new approach where, thanks to an interface between the FLUKA Monte Carlo transport code and the BIANCA biophysical model, the RBE dependence on particle energy and absorbed dose was also considered. Furthermore, we included in the considered SPE spectra primary particles heavier than protons, which in many studies are neglected. This approach was then applied to the October 2003 SPE (the most intense SPE of solar cycle 23, also known as "Halloween event") and the January 2005 event, which was characterized by a lower fluence but a harder spectrum, i.e., with higher-energy particles. The calculation outcomes were then discussed and compared with the current dose limits established for skin and blood forming organs in case of 30-days missions. This work showed that the BIANCA model, if interfaced to a radiation transport code, can be used to calculate the RBE values associated to Solar Particle Events. More generally, this work emphasizes the importance of taking into account the RBE dependence on particle energy and dose when calculating equivalent doses.

Evaluation of specific RBE in different cells of hippocampus under high-dose proton irradiation in rats.

The study aimed to determine the specific relative biological effectiveness (RBE) of various cells in the hippocampus following proton irradiation. Sixty Sprague-Dawley rats were randomly allocated to 5 groups receiving 20 or 30 Gy of proton or photon irradiation. Pathomorphological neuronal damage in the hippocampus was assessed using Hematoxylin-eosin (HE) staining. The expression level of NeuN, Nestin, Caspase-3, Olig2, CD68 and CD45 were determined by immunohistochemistry (IHC). The RBE range established by comparing the effects of proton and photon irradiation at equivalent biological outcomes. Proton20Gy induced more severe damage to neurons than photon20Gy, but showed no difference compared to photon30Gy. The RBE of neuron was determined to be 1.65. Similarly, both proton20Gy and proton30Gy resulted in more inhibition of oligodendrocytes and activation of microglia in the hippocampal regions than photon20Gy and photon30Gy. However, the expression of Olig2 was higher and CD68 was lower in the proton20Gy group than in the photon30Gy group. The RBE of oligodendrocyte and microglia was estimated to be between 1.1 to 1.65. For neural stem cells (NSCs) and immune cells, there were no significant difference in the expression of Nestin and CD45 between proton and photon irradiation (both 20 and 30 Gy). Therefore, the RBE for NSCs and immune cell was determined to be 1.1. These findings highlight the varying RBE values of different cells in the hippocampus in vivo. Moreover, the actual RBE of the hippocampus may be higher than 1.1, suggesting that using as RBE value of 1.1 in clinical practice may underestimate the toxicities induced by proton radiation.

Dose and dose rate dependence of the tissue sparing effect at ultra-high dose rate studied for proton and electron beams using the zebrafish embryo model.

PURPOSE: A better characterization of the dependence of the tissue sparing effect at ultra-high dose rate (UHDR) on physical beam parameters (dose, dose rate, radiation quality) would be helpful towards a mechanistic understanding of the FLASH effect and for its broader clinical translation. To address this, a comprehensive study on the normal tissue sparing at UHDR using the zebrafish embryo (ZFE) model was conducted. METHODS: One-day-old ZFE were irradiated over a wide dose range (15-95 Gy) in three different beams (proton entrance channel, proton spread out Bragg peak and 30 MeV electrons) at UHDR and reference dose rate. After irradiation the ZFE were incubated for 4 days and then analyzed for four different biological endpoints (pericardial edema, curved spine, embryo length and eye diameter). RESULTS: Dose-effect curves were obtained and a sparing effect at UHDR was observed for all three beams. It was demonstrated that proton relative biological effectiveness and UHDR sparing are both relevant to predict the resulting dose response. Dose dependent FLASH modifying factors (FMF) for ZFE were found to be compatible with rodent data from the literature. It was found that the UHDR sparing effect saturates at doses above âˆ¼ 50 Gy with an FMF of âˆ¼ 0.7-0.8. A strong dose rate dependence of the tissue sparing effect in ZFE was observed. The magnitude of the maximum sparing effect was comparable for all studied biological endpoints. CONCLUSION: The ZFE model was shown to be a suitable pre-clinical high-throughput model for radiobiological studies on FLASH radiotherapy, providing results comparable to rodent models. This underlines the relevance of ZFE studies for FLASH radiotherapy research.

Benchmarking proton RBE models.

Objective.To biologically optimise proton therapy, models which can accurately predict variations in proton relative biological effectiveness (RBE) are essential. Current phenomenological models show large disagreements in RBE predictions, due to different model assumptions and differences in the data to which they were fit. In this work, thirteen RBE models were benchmarked against a comprehensive proton RBE dataset to evaluate predictions when all models are fit using the same data and fitting techniques, and to assess the statistical robustness of the models.Approach.Model performance was initially evaluated by fitting to the full dataset, and then a cross-validation approach was applied to assess model generalisability and robustness. The impact of weighting the fit and the choice of biological endpoint (either single or multiple survival levels) was also evaluated.Main results.Fitting the models to a common dataset reduced differences between their predictions, however significant disagreements remained due to different underlying assumptions. All models performed poorly under cross-validation in the weighted fits, suggesting that some uncertainties on the experimental data were significantly underestimated, resulting in over-fitting and poor performance on unseen data. The simplest model, which depends linearly on the LET but has no tissue or dose dependence, performed best for a single survival level. However, when fitting to multiple survival levels simultaneously, more complex models with tissue dependence performed better. All models had significant residual uncertainty in their predictions compared to experimental data.Significance.This analysis highlights that poor quality of error estimation on the dose response parameters introduces substantial uncertainty in model fitting. The significant residual error present in all approaches illustrates the challenges inherent in fitting to large, heterogeneous datasets and the importance of robust statistical validation of RBE models.

Hypoxia-informed RBE-weighted beam orientation optimization for intensity modulated proton therapy.

BACKGROUND: Variable relative biological effectiveness (RBE) models in treatment planning have been proposed to optimize the therapeutic ratio of proton therapy. It has been reported that proton RBE decreases with increasing tumor oxygen level, offering an opportunity to address hypoxia-related radioresistance with RBE-weighted optimization. PURPOSE: Here, we obtain a voxel-level estimation of partial oxygen pressure to weigh RBE values in a single biologically informed beam orientation optimization (BOO) algorithm. METHODS: Three glioblastoma patients with [18 F]-fluoromisonidazole (FMISO)-PET/CT images were selected from the institutional database. Oxygen values were derived from tracer uptake using a nonlinear least squares curve fitting. McNamara RBE, calculated from proton dose, was then weighed using oxygen enhancement ratios (OER) for each voxel and incorporated into the dose fidelity term of the BOO algorithm. The nonlinear optimization problem was solved using a split-Bregman approach, with FISTA as the solver. The proposed hypoxia informed RBE-weighted method (HypRBE) was compared to dose fidelity terms using the constant RBE of 1.1 (cRBE) and the normoxic McNamara RBE model (RegRBE). Tumor homogeneity index (HI), maximum biological dose (Dmax), and D95%, as well as OAR therapeutic index (TI = gEUDCTV /gEUDOAR ) were evaluated along with worst-case statistics after normalization to normal tissue isotoxicity. RESULTS: Compared to [cRBE, RegRBE], HypRBE increased tumor HI, Dmax, and D95% across all plans by on average [31.3%, 31.8%], [48.6%, 27.1%], and [50.4%, 23.8%], respectively. In the worst-case scenario, the parameters increase on average by [12.5%, 14.7%], [7.3%,-8.9%], and [22.3%, 2.1%]. Despite increased OAR Dmean and Dmax by [8.0%, 3.0%] and [13.1%, -0.1%], HypRBE increased average TI by [22.0%, 21.1%]. Worst-case OAR Dmean, Dmax, and TI worsened by [17.9%, 4.3%], [24.5%, -1.2%], and [9.6%, 10.5%], but in the best cases, HypRBE escalates tumor coverage significantly without compromising OAR dose, increasing the therapeutic ratio. CONCLUSIONS: We have developed an optimization algorithm whose dose fidelity term accounts for hypoxia-informed RBE values. We have shown that HypRBE selects bE:\Alok\aaeams better suited to deliver high physical dose to low RBE, hypoxic tumor regions while sparing the radiosensitive normal tissue.

Beam quality correction factors for ionization chambers in a 0.35 T magnetic resonance (MR)-linac - A Monte Carlo study.

PURPOSE: The purpose of this study was to directly calculate [Formula: see text] correction factors for four cylindrical ICs for a 0.35 T MR-linac using the Monte Carlo (MC) method. METHODS: A previously-validated TOPAS/GEANT4 MC head model of the 0.35 T MR-linac was employed. The MR-compatible Exradin A12, A1SL, A26, and A28 cylindrical ICs were modeled considering the dead volume in the air cavity. The [Formula: see text] correction factor was determined for initial electron energies of 5-7 MeV. The correction factor was calculated for all four angular orientations in the lateral plane. The impact of the 0.35 T magnetic field on the IC response was also investigated. RESULTS: The maximum beam quality dependence in the [Formula: see text] exhibited by the A12, A1SL, A26, and A28 ICs was 1.10 %, 2.17 %, 0.81 %, and 1.75 %, respectively, considering all angular orientations. The magnetic field dependence was < 1 % and the maximum [Formula: see text] correction was < 2 % when the detector was aligned along the direction of the magnetic field at 0° and 180° angles. The A12 IC over-responded up to 5.40 % for the orthogonal orientation. An asymmetry in the response of up to 8.30 % was noted for the A28 IC aligned at 90° and 270° angles. CONCLUSIONS: A parallel orientation for the IC, with respect to the magnetic field, is recommended for reference dosimetry in MRgRT. Both over and under-response in the IC signal was noted for the orthogonal orientations, which is highly dependent on the cavity diameter, cavity length, and the dead volume.

Equivalent uniform RBE-weighted dose in eye plaque brachytherapy.

BACKGROUND: Brachytherapy for ocular melanoma is based on the application of eye plaques with different spatial dose nonuniformity, time-dependent dose rates and relative biological effectiveness (RBE). PURPOSE: We propose a parameter called the equivalent uniform RBE-weighted dose (EUDRBE) that can be used for quantitative characterization of integrated cell survival in radiotherapy modalities with the variable RBE, dose nonuniformity and dose rate. The EUDRBE is applied to brachytherapy with 125I eye plaques designed by the Collaborative Ocular Melanoma Study (COMS). METHODS: The EUDRBE is defined as the uniform dose distribution with RBE = 1 that causes equal cell survival for a given nonuniform dose distribution with the variable RBE > 1. The EUDRBE can be used for comparison of cell survival for nonuniform dose distributions with different RBE, because they are compared to the reference dose with RBE = 1. The EUDRBE is applied to brachytherapy with 125I COMS eye plaques that are characterized by a steep dose gradient in tumor base-apex direction, protracted irradiation during time intervals of 3-8 days, and variable dose-rate dependent RBE with a maximum of about 1.4. The simulations are based on dose of 85 Gy prescribed to the farthest intraocular extent of the tumor (tumor apex). To compute the EUDRBE in eye plaque brachytherapy and correct for protracted irradiation, the distributions of physical dose have been converted to non-uniform distributions of biologically effective dose (BED) to include the biological effects of sublethal cellular repair, Our radiobiological analysis considers the combined effects of different time-dependent dose rates, spatial dose non-uniformity, dose fractionation and different RBE and can be used to derive optimized dose regimens brachytherapy. RESULTS: Our simulations show that the EUDRBE increases with the prescription depths and the maximum increase may achieve 6% for the tumor height of 12 mm. This effect stems from a steep dose gradient within the tumor that increases with the prescription depth. The simulations also show that the EUDRBE increase may achieve 12% with increasing the dose rate when implant duration decreases. The combined effect of dose nonuniformity and dose rate may change the EUDRBE up to 18% for the same dose prescription of 85 Gy to tumor apex. The absolute dose range of 48-61 Gy (RBE) for the EUDRBE computed using 4 or 5 fractions is comparable to the dose prescriptions used in stereotactic body radiation therapy (SBRT) with megavoltage X-rays (RBE = 1) for different cancers. The tumor control probabilities in SBRT and eye plaque brachytherapy are very similar at the level of 80% or higher that support the hypothesis that the selected approximations for the EUDRBE are valid. CONCLUSIONS: The computed range of the EUDRBE in 125I COMS eye plaque brachytherapy suggests that the selected models and hypotheses are acceptable. The EUDRBE can be useful for analysis of treatment outcomes and comparison of different dose regimens in eye plaque brachytherapy.

Primary study of the relative and compound biological effectiveness model for boron neutron capture therapy based on nanodosimetry.

BACKGROUND: The current radiobiological model employed for boron neutron capture therapy (BNCT) treatment planning, which relies on microdosimetry, fails to provide an accurate representation the biological effects of BNCT. The precision in calculating the relative biological effectiveness (RBE) and compound biological effectiveness (CBE) plays a pivotal role in determining the therapeutic efficacy of BNCT. Therefore, this study focuses on how to improve the accuracy of the biological effects of BNCT. PURPOSE: The purpose of this study is to propose new radiation biology models based on nanodosimetry to accurately assess RBE and CBE for BNCT. METHODS: Nanodosimetry, rooted in ionization cluster size distributions (ICSD), introduces a novel approach to characterize radiation quality by effectively delineating RBE through the ion track structure at the nanoscale. In the context of prior research, this study presents a computational model for the nanoscale assessment of RBE and CBE. We establish a simplified model of DNA chromatin fiber using the Monte Carlo code TOPAS-nBio to evaluate the applicability of ICSD to BNCT and compute nanodosimetric parameters. RESULTS: Our investigation reveals that both homogeneous and heterogeneous nanodosimetric parameters, as well as the corresponding biological model coefficients α and ß, along with RBE values, exhibit variations in response to varying intracellular 10B concentrations. Notably, the nanodosimetric parameter M 1 C 2 $M_1^{{{\mathrm{C}}}_2}$ effectively captures the fluctuations in model coefficients α and RBE. CONCLUSION: Our model facilitates a nanoscale analysis of BNCT, enabling predictions of nanodosimetric quantities for secondary ions as well as RBE, CBE, and other essential biological metrics related to the distribution of boron. This contribution significantly enhances the precision of RBE calculations and holds substantial promise for future applications in treatment planning.

Proton relative biological effectiveness for the induction of DNA double strand breaks based on Geant4.

Uncertainties in the relative biological effectiveness (RBE) of proton remains a major barrier to the biological optimization of proton therapy. A large amount of experimental data suggest that proton RBE is variable. As an evolving Monte Carlo code toolkit, Geant4-DNA is able to simulate the initial DNA damage caused by particle beams through physical and chemical interactions at the nanometer scale over a short period of time. This contributes to evaluating the radiobiological effects induced by ionizing radiation. Based on the Geant4-DNA toolkit, this study constructed a DNA geometric model containing 6.32Gbp, simulated the relationship between radiochemical yields (G-values) and their corresponding chemical constructors, and calculated a detailed calculation of the sources of damage and the complexity of damage in DNA strand breaks. The damage model constructed in this study can simulate the relative biological effectiveness (RBE) in the proton Bragg peak region. The results indicate that: (1) When the electron energy is below 400 keV, the yield of OH·account for 18.1% to 25.3% of the total water radiolysis yields. (2) Under the influence of histone clearance function, the yield of indirect damage account for over 72.93% of the yield of DNA strand breaks (SBs). When linear energy transfer (LET) increased from 29.79 (keV/µm) to 64.29 (keV/µm), the yield of double strand breaks (DSB) increased from 17.27% to 32.65%. (3) By investigating the effect of proton Bragg peak depth on the yield of direct DSB (DSBdirect) and total DSB (DSBtotal), theRBEDSBtotandRBEDSBdirlevels of cells show that the RBE value of protons reaches 2.2 in the Bragg peak region.

Sensitivity of a mini-TEPC to radiation quality variations in clinical proton beams.

PURPOSE: This work aims at studying the sensitivity of a miniaturized Tissue-Equivalent Proportional Counter to variations of beam quality in clinical radiation fields, to further investigate its performances as radiation quality monitor. METHODS: Measurements were taken at the CATANA facility (INFN-LNS, Catania, Italy), in a monoenergetic and an energy-modulated proton beam with the same initial energy of 62 MeV. PMMA layers were placed in front of the detector to measure at different depths along the depth-dose profile. The frequency- and dose-mean lineal energy were compared to the track- and dose-averaged LET calculated by Monte Carlo simulations. A microdosimetric evaluation of the Relative Biological Effectiveness (RBE) was performed and compared with cell survival experiments. RESULTS: Microdosimetric distributions measured at identical depths in the two beams show spectral differences reflecting their different radiation quality. Discrepancies are most evident at depths corresponding to the Spread-Out Bragg Peak, while spectra at the entrance and in the dose fall-off regions are similar. This can be explained by the different energy components that compose the pristine and spread-out peaks at each depth. The trend of microdosimetric mean values matches that of calculated LET averages along the entire penetration depth, and the microdosimetric estimation of RBE is consistent with radiobiological data not only at 2 Gy but also at lower dose levels, such as those absorbed by healthy tissues. CONCLUSIONS: The mini-TEPC is sensitive to differences in radiation quality resulting from different modulations of the proton beam, confirming its potential for beam quality monitoring in proton therapy.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool.

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Boron-Containing MOF Nanoparticles with Stable Metabolism in U87-MG Cells Combining Microdosimetry To Evaluate Relative Biological Effectiveness of Boron Neutron Capture Therapy.

Accurate prediction of the relative biological effectiveness (RBE) of boron neutron capture therapy (BNCT) is challenging. The therapy is different from other radiotherapy; the dynamic distribution of boron-containing compounds in tumor cells affects the therapeutic outcome considerably and hampers accurate measurement of the neutron-absorbed dose. Herein, we used boron-containing metal-organic framework nanoparticles (BMOFs) with high boron content to target U87-MG cells and maintain the concentration of the 10B isotope in cells. The content of boron in the cells could maintain 90% (60 ppm) within 20 min compared with that at the beginning; therefore, the accurate RBE of BNCT can be acquired. The effects of BNCT upon cells after neutron irradiation were observed, and the neutron-absorbed dose was obtained by Monte Carlo simulations. The RBE of BMOFs was 6.78, which was 4.1-fold higher than that of a small-molecule boron-containing agent (boric acid). The energy spectrum of various particles was analyzed by Monte Carlo simulations, and the RBE was verified theoretically. Our results suggested that the use of nanoparticle-based boron carriers in BNCT may have many advantages and that maintaining a stable boron distribution within cells may significantly improve the efficiency of BNCT.

Cell Survival Computation via the Generalized Stochastic Microdosimetric Model (GSM2); Part II: Numerical Results.

In the present paper we numerically investigate, using Monte Carlo simulation, the theoretical results predicted by the Generalized Stochastic Microdosimetric Model (GSM2), as shown in the published companion paper. Taking advantage of the particle irradiation data ensemble (PIDE) dataset, we calculated GSM2 biological parameters of human salivary gland (HSG) and V79 cell lines. Further, exploiting the TOPAS-microdosimetric extension, we simulated the microdosimetric spectra of different radiation fields of therapeutic interest generated by four different ions (protons, helium-4, carbon-12 and oxygen-16) each at three different residual ranges. We investigated the properties of the initial damage distributions as well as the cell survival curve predicted by GSM2, focusing especially on the non-Poissonian effects naturally included in the model. GSM2 successfully computed cell survival curves, accurately describing experimental behavior even under challenging LET and dose conditions.

Enhanced RBE of Particle Radiation Depends on Beam Size in the Micrometer Range.

High-linear energy transfer (LET) radiation, such as heavy ions is associated with a higher relative biological effectiveness (RBE) than low-LET radiation, such as photons. Irradiation with low- and high-LET particles differ in the interaction with the cellular matter and therefore in the spatial dose distribution. When a single high-LET particle interacts with matter, it results in doses of up to thousands of gray (Gy) locally concentrated around the ion trajectory, whereas the mean dose averaged over the target, such as a cell nucleus is only in the range of a Gy. DNA damage therefore accumulates in this small volume. In contrast, up to hundreds of low-LET particle hits are required to achieve the same mean dose, resulting in a quasi-homogeneous damage distribution throughout the cell nucleus. In this study, we investigated the dependence of RBE from different spatial dose depositions using different focused beam spot sizes of proton radiation with respect to the induction of chromosome aberrations and clonogenic cell survival. Human-hamster hybrid (AL) as well as Chinese hamster ovary cells (CHO-K1) were irradiated with focused low LET protons of 20 MeV (LET = 2.6 keV/µm) beam energy with a mean dose of 1.7 Gy in a quadratic matrix pattern with point spacing of 5.4 × 5.4 µm2 and 117 protons per matrix point at the ion microbeam SNAKE using different beam spot sizes between 0.8 µm and 2.8 µm (full width at half maximum). The dose-response curves of X-ray reference radiation were used to determine the RBE after a 1.7 Gy dose of radiation. The RBE for the induction of dicentric chromosomes and cell inactivation was increased after irradiation with the smallest beam spot diameter (0.8 µm for chromosome aberration experiments and 1.0 µm for cell survival experiments) compared to homogeneous proton radiation but was still below the RBE of a corresponding high LET single ion hit. By increasing the spot size to 1.6-1.8 µm, the RBE decreased but was still higher than for homogeneously distributed protons. By further increasing the spot size to 2.7-2.8 µm, the RBE was no longer different from the homogeneous radiation. Our experiments demonstrate that varying spot size of low-LET radiation gradually modifies the RBE. This underlines that a substantial fraction of enhanced RBE originates from inhomogeneous energy concentrations on the µm scale (mean intertrack distances of low-LET particles below 0.1 µm) and quantifies the link between such energy concentration and RBE. The missing fraction of RBE enhancement when comparing with high-LET ions is attributed to the high inner track energy deposition on the nanometer scale. The results are compared with model results of PARTRAC and LEM for chromosomal aberration and cell survival, respectively, which suggest mechanistic interpretations of the observed radiation effects.

Effects of cellular radioresponse on therapeutic helium-, carbon-, oxygen-, and neon-ion beams: a simulation study.

Objective. Helium, oxygen, and neon ions in addition to carbon ions will be used for hypofractionated multi-ion therapy to maximize the therapeutic effectiveness of charged-particle therapy. To use new ions in cancer treatments based on the dose-fractionation protocols established in carbon-ion therapy, this study examined the cell-line-specific radioresponse to therapeutic helium-, oxygen-, and neon-ion beams within wide dose ranges.Approach. Response of cells to ions was described by the stochastic microdosimetric kinetic model. First, simulations were made for the irradiation of one-field spread-out Bragg peak beams in water with helium, carbon, oxygen, and neon ions to achieve uniform survival fractions at 37%, 10%, and 1% for human salivary gland tumor (HSG) cells, the reference cell line for the Japanese relative biological effectiveness weighted dose system, within the target region defined at depths from 90 to 150 mm. The HSG cells were then replaced by other cell lines with different radioresponses to evaluate differences in the biological dose distributions of each ion beam with respect to those of carbon-ion beams.Main results. For oxygen- and neon-ion beams, the biological dose distributions within the target region were almost equivalent to those of carbon-ion beams, differing by less than 5% in most cases. In contrast, for helium-ion beams, the biological dose distributions within the target region were largely different from those of carbon-ion beams, more than 10% in several cases.Significance.From the standpoint of tumor control evaluated by the clonogenic cell survival, this study suggests that the dose-fractionation protocols established in carbon-ion therapy could be reasonably applied to oxygen- and neon-ion beams while some modifications in dose prescription would be needed when the protocols are applied to helium-ion beams. This study bridges the gap between carbon-ion therapy and hypofractionated multi-ion therapy.

Deep learning-based synthetic dose-weighted LET map generation for intensity modulated proton therapy.

The advantage of proton therapy as compared to photon therapy stems from the Bragg peak effect, which allows protons to deposit most of their energy directly at the tumor while sparing healthy tissue. However, even with such benefits, proton therapy does present certain challenges. The biological effectiveness differences between protons and photons are not fully incorporated into clinical treatment planning processes. In current clinical practice, the relative biological effectiveness (RBE) between protons and photons is set as constant 1.1. Numerous studies have suggested that the RBE of protons can exhibit significant variability. Given these findings, there is a substantial interest in refining proton therapy treatment planning to better account for the variable RBE. Dose-average linear energy transfer (LETd) is a key physical parameter for evaluating the RBE of proton therapy and aids in optimizing proton treatment plans. Calculating precise LETddistributions necessitates the use of intricate physical models and the execution of specialized Monte-Carlo simulation software, which is a computationally intensive and time-consuming progress. In response to these challenges, we propose a deep learning based framework designed to predict the LETddistribution map using the dose distribution map. This approach aims to simplify the process and increase the speed of LETdmap generation in clinical settings. The proposed CycleGAN model has demonstrated superior performance over other GAN-based models. The mean absolute error (MAE), peak signal-to-noise ratio and normalized cross correlation of the LETdmaps generated by the proposed method are 0.096 ± 0.019 keVµm-1, 24.203 ± 2.683 dB, and 0.997 ± 0.002, respectively. The MAE of the proposed method in the clinical target volume, bladder, and rectum are 0.193 ± 0.103, 0.277 ± 0.112, and 0.211 ± 0.086 keVµm-1, respectively. The proposed framework has demonstrated the feasibility of generating synthetic LETdmaps from dose maps and has the potential to improve proton therapy planning by providing accurate LETdinformation.

Experimental determination of magnetic field quality conversion factors for eleven ionization chambers in 1.5 T and 0.35 T MR-linac systems.

BACKGROUND: The static magnetic field present in magnetic resonance (MR)-guided radiotherapy systems can influence dose deposition and charged particle collection in air-filled ionization chambers. Thus, accurately quantifying the effect of the magnetic field on ionization chamber response is critical for output calibration. Formalisms for reference dosimetry in a magnetic field have been proposed, whereby a magnetic field quality conversion factor kB,Q is defined to account for the combined effects of the magnetic field on the radiation detector. Determination of kB,Q in the literature has focused on Monte Carlo simulation studies, with experimental validation limited to only a few ionization chamber models. PURPOSE: The purpose of this study is to experimentally measure kB,Q for 11 ionization chamber models in two commercially available MR-guided radiotherapy systems: Elekta Unity and ViewRay MRIdian. METHODS: Eleven ionization chamber models were characterized in this study: Exradin A12, A12S, A28, and A26, PTW T31010, T31021, and T31022, and IBA FC23-C, CC25, CC13, and CC08. The experimental method to measure kB,Q utilized cross-calibration against a reference Exradin A1SL chamber. Absorbed dose to water was measured for the reference A1SL chamber positioned parallel to the magnetic field with its centroid placed at the machine isocenter at a depth of 10 cm in water for a 10 × 10 cm2 field size at that depth. Output was subsequently measured with the test chamber at the same point of measurement. kB,Q for the test chamber was computed as the ratio of reference dose to test chamber output, with this procedure repeated for each chamber in each MR-guided radiotherapy system. For the high-field 1.5 T Elekta Unity system, the dependence of kB,Q on the chamber orientation relative to the magnetic field was quantified by rotating the chamber about the machine isocenter. RESULTS: Measured kB,Q values for our test dataset of ionization chamber models ranged from 0.991 to 1.002, and 0.995 to 1.004 for the Elekta Unity and ViewRay MRIdian, respectively, with kB,Q tending to increase as the chamber sensitive volume increased. Measured kB,Q values largely agreed within uncertainty to published Monte Carlo simulation data and available experimental data. kB,Q deviation from unity was minimized for ionization chamber orientation parallel or antiparallel to the magnetic field, with increased deviations observed at perpendicular orientations. Overall (k = 1) uncertainty in the experimental determination of the magnetic field quality conversion factor, kB,Q was 0.71% and 0.72% for the Elekta Unity and ViewRay MRIdian systems, respectively. CONCLUSIONS: For a high-field MR-linac, the characterization of ionization chamber performance as angular orientation varied relative to the magnetic field confirmed that the ideal orientation for output calibration is parallel. For most of these chamber models, this study represents the first experimental characterization of chamber performance in clinical MR-linac beams. This is a critical step toward accurate output calibration for MR-guided radiotherapy systems and the measured kB,Q values will be an important reference data source for forthcoming MR-linac reference dosimetry protocols.