BACKGROUND: Radiobiological effectiveness of radiation in cancer treatment can be studied at different scales (molecular till organ scale) and different time post irradiation. The production of free radicals and reactive oxygen species during water radiolysis is particularly relevant to understand the fundamental mechanisms playing a role in observed biological outcomes. The development and validation of Monte Carlo tools integrating the simulation of physical, physico-chemical and chemical stages after radiation is very important to maintain with experiments.
OBJECTIVE: Therefore, in this study, we propose to validate a new Geant4-DNA chemistry module through the simulation of water radiolysis and Fricke dosimetry experiments on a proton preclinical beam line.
METHODS: In this study, we used the GATE Monte Carlo simulation platform (version 9.3) to simulate a 67.5 MeV proton beam produced with the ARRONAX isochronous cyclotron (IBA Cyclone 70XP) at conventional dose rate (0.2 Gy/s) to simulate the irradiation of ultra-pure liquid water samples and Fricke dosimeter. We compared the depth dose profile with measurements performed with a plane parallel Advanced PTW 34045 Markus ionization chamber. Then, a new Geant4-DNA chemistry application proposed from Geant4 version 11.2 has been used to assess the evolution of HO • ${\\mathrm{HO}}^ \\bullet $ , e aq - ${\\mathrm{e}}_{{\\mathrm{aq}}}^ - $ , H 3 O + ${{\\mathrm{H}}}_3{{\\mathrm{O}}}^ + $ , H 2 O 2 ${{\\mathrm{H}}}_2{{\\mathrm{O}}}_2$ , H 2 ${{\\mathrm{H}}}_2$ , HO 2 • ${\\mathrm{HO}}_2^ \\bullet $ , HO 2 - , O 2 • - ${\\mathrm{HO}}_2^ - ,{\\mathrm{\\ O}}_2^{ \\bullet - }$ and HO - ${\\mathrm{HO}}^ - $ reactive species along time until 1-h post-irradiation. In particular, the effect of oxygen and pH has been investigated through comparisons with experimental measurements of radiolytic yields for H 2 O 2 ${{\\mathrm{H}}}_2{{\\mathrm{O}}}_2$ and Fe3+.
RESULTS: GATE simulations reproduced, within 4%, the depth dose profile in liquid water. With Geant4-DNA, we were able to reproduce experimental H 2 O 2 ${{\\mathrm{H}}}_2{{\\mathrm{O}}}_2$ radiolytic yields 1-h post-irradiation in aerated and deaerated conditions, showing the impact of small changes in oxygen concentrations on species evolution along time. For the Fricke dosimeter, simulated G(Fe3+) is 15.97 ± 0.2 molecules/100 eV which is 11% higher than the measured value (14.4 ± 04 molecules/100 eV).
CONCLUSIONS: These results aim to be consolidated by new comparisons involving other radiolytic species, such as e aq - ${\\mathrm{e}}_{{\\mathrm{aq}}}^ - $ or , O 2 • - $,{\\mathrm{\\ O}}_2^{ \\bullet - }$ to further study the mechanisms underlying the FLASH effect observed at ultra-high dose rates (UHDR).

Targeted alpha therapy is an emerging alternative for palliative therapy of a wide range of tumor types. Data from preclinicaland clinical research demonstrates a high potential for the selective killing of tumor cells and minimal toxicity to surroundinghealthy tissues. This article summarizes the developmental stages of alpha-targeted therapy from benchtop to commercialization.It discusses fundamental properties, production pathways, microdosimetry, and possible targeting vectors. Proper coverage hasalso been given to comparing it with other standard treatment procedures while exploring clinical applications of alpha emitters.In the end, like other therapies, the challenges it faces and its future impact on personalized medicine are also illustrated.

In the next decade, the International Commission on Radiological Protection (ICRP) will issue the next set of general recommendations, for which evaluation of relative biological effectiveness (RBE) for various types of tissue reactions would be needed. ICRP has recently classified diseases of the circulatory system (DCS) as a tissue reaction, but has not recommended RBE for DCS. We therefore evaluated the mean and uncertainty of RBE for DCS by applying a microdosimetric kinetic model specialized for RBE estimation of tissue reactions. For this purpose, we analyzed several RBE data for DCS determined by past animal experiments and evaluated the radius of the subnuclear domain best fit to each experiment as a single free parameter included in the model. Our analysis suggested that RBE for DCS tends to be lower than that for skin reactions, and their difference was borderline significant due to large variances of the evaluated parameters. We also found that RBE for DCS following mono-energetic neutron irradiation of the human body is much lower than that for skin reactions, particularly at the thermal energy and around 1 MeV. This tendency is considered attributable not only to the intrinsic difference of neutron RBE between skin reactions and DCS but also to the difference in the contributions of secondary γ-rays to the total absorbed doses between their target organs. These findings will help determine RBE by ICRP for preventing tissue reactions.

The biological mechanisms triggered by low-dose exposure still need to be explored in depth. In this study, the potential mechanisms of low-dose radiation when irradiating the BEAS-2B cell lines with a Cs-137 gamma-ray source were investigated through simulations and experiments. Monolayer cell population models were constructed for simulating and analyzing distributions of nucleus-specific energy within cell populations combined with the Monte Carlo method and microdosimetric analysis. Furthermore, the 10 × Genomics single-cell sequencing technology was employed to capture the heterogeneity of individual cell responses to low-dose radiation in the same irradiated sample. The numerical uncertainties can be found both in the specific energy distribution in microdosimetry and in differential gene expressions in radiation cytogenetics. Subsequently, the distribution of nucleus-specific energy was compared with the distribution of differential gene expressions to guide the selection of differential genes bioinformatics analysis. Dose inhomogeneity is pronounced at low doses, where an increase in dose corresponds to a decrease in the dispersion of cellular-specific energy distribution. Multiple screening of differential genes by microdosimetric features and statistical analysis indicate a number of potential pathways induced by low-dose exposure. It also provides a novel perspective on the selection of sensitive biomarkers that respond to low-dose radiation.

We present a numerical method for studying reversible electroporation on normal and cancerous cervical cells. This microdosimetry analysis builds on a unique approach for extracting contours of free and overlapping cervical cells in the cluster from the Extended Depth of Field (EDF) images. The algorithm used for extracting the contours is a joint optimization of multiple-level set function along with the Gaussian mixture model and Maximally Stable Extremal Regions. These contours are then exported to a multi-physics domain solver, where a variable frequency pulsed electric field is applied. The trans-Membrane voltage (TMV) developed across the cell membrane is computed using the Maxwell equation coupled with a statistical approach, employing the asymptotic Smoluchowski equation. The numerical model was validated by successful replication of existing experimental configurations that employed low-frequency uni-polar pulses on the overlapping cells to obtain reversible electroporation, wherein, several overlapping clumps of cervical cells were targeted. For high-frequency calculation, a combination of normal and cancerous cells is introduced to the computational domain. The cells are assumed to be dispersive and the Debye dispersion equation is used for further calculations. We also present the resulting strength-duration relationship for achieving the threshold value of electroporation between the normal and cancerous cervical cells due to their size and conductivity differences. The dye uptake modulation during the high-frequency electric field electroporation is further advocated by a mathematical model.

Objective.This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50-250 MeV at different depths in water using TOPAS Monte Carlo code.Approach.The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions.Main results.The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04-1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15-1.33, 1.13-1.21, 1.11-1.17, 1.13-1.18 and 1.17-1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively.Significance.The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

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.

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.

OBJECTIVE: 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.

BACKGROUND: Although the benefits of breast screening and early diagnosis are known for reducing breast cancer mortality rates, the effects and risks of low radiation doses to the cells in the breast are still ongoing topics of study.
OBJECTIVE: To study specific energy distributions ( f ( z , D g ) $f(z,D_{g})$ ) in cytoplasm and nuclei of cells corresponding to glandular tissue for different x-ray breast imaging modalities.
METHODS: A cubic lattice (500 μm length side) containing 4064 spherical cells was irradiated with photons loaded from phase space files with varying glandular voxel doses ( D g $D_{g}$ ). Specific energy distributions were scored for nucleus and cytoplasm compartments using the PENELOPE (v. 2018) + penEasy (v. 2020) Monte Carlo (MC) code. The phase space files, generated in part I of this work, were obtained from MC simulations in a voxelized anthropomorphic phantom corresponding to glandular voxels for different breast imaging modalities, including digital mammography (DM), digital breast tomosynthesis (DBT), contrast enhanced digital mammography (CEDM) and breast CT (BCT).
RESULTS: In general, the average specific energy in nuclei is higher than the respective glandular dose scored in the same region, by up to 10%. The specific energy distributions for nucleus and cytoplasm are directly related to the magnitude of the glandular dose in the voxel ( D g $D_{g}$ ), with little dependence on the spatial location. For similar D g $D_{g}$ values, f ( z , D g ) $f(z,D_{g})$ for nuclei is different between DM/DBT and CEDM/BCT, indicating that distinct x-ray spectra play significant roles in f ( z , D g ) $f(z,D_{g})$ . In addition, this behavior is also present when the specific energy distribution ( F g ( z ) $F_{g}(z)$ ) is considered taking into account the GDD in the breast.
CONCLUSIONS: Microdosimetry studies are complementary to the traditional macroscopic breast dosimetry based on the mean glandular dose (MGD). For the same MGD, the specific energy distribution in glandular tissue varies between breast imaging modalities, indicating that this effect could be considered for studying the risks of exposing the breast to ionizing radiation.