Liquid water radiolysis induced by secondary electrons generated from MeV-energy carbon ions.

Fast ion beams induce damage to deoxyribonucleic acid (DNA) by chemical products, including secondary electrons, produced from interaction with liquid water in living cells. However, the production process of these chemical products in the Bragg peak region used in particle therapy is not fully understood. To investigate this process, we conducted experiments to evaluate the radiolytic yields produced when a liquid water jet in vacuum is irradiated with MeV-energy carbon beams. We used secondary ion mass spectrometry to measure the products, such as hydronium cations (H3O+) and hydroxyl anions (OH-), produced along with ·OH radicals, which are significant inducers of DNA damage formation. In addition, we simulated the ionization process in liquid water by incident ions and secondary electrons using a Monte Carlo code for radiation transport. Our results showed that secondary electrons, rather than incident ions, are the primary cause of ionization in water. We found that the production yield of H3O+ or OH- was linked to the frequency of ionization by secondary electrons in water, with these electrons having energies between 10.9 and 550 eV. These electrons are responsible for ionizing the outer-shell electrons of water molecules. Finally, we present that the elementary processes contribute to advancing radiation biophysics and biochemistry, which study the formation mechanism of DNA damage.

An Analytical Method for Quantifying the Yields of DNA Double-Strand Breaks Coupled with Strand Breaks by γ-H2AX Focus Formation Assay Based on Track-Structure Simulation.

Complex DNA double-strand break (DSB), which is defined as a DSB coupled with additional strand breaks within 10 bp in this study, induced after ionizing radiation or X-rays, is recognized as fatal damage which can induce cell death with a certain probability. In general, a DSB site inside the nucleus of live cells can be experimentally detected using the γ-H2AX focus formation assay. DSB complexity is believed to be detected by analyzing the focus size using such an assay. However, the relationship between focus size and DSB complexity remains uncertain. In this study, using Monte Carlo (MC) track-structure simulation codes, i.e., an in-house WLTrack code and a Particle and Heavy Ion Transport code System (PHITS), we developed an analytical method for qualifying the DSB complexity induced by photon irradiation from the microscopic image of γ-H2AX foci. First, assuming that events (i.e., ionization and excitation) potentially induce DNA strand breaks, we scored the number of events in a water cube (5.03 × 5.03 × 5.03 nm3) along electron tracks. Second, we obtained the relationship between the number of events and the foci size experimentally measured by the γ-H2AX focus formation assay. Third, using this relationship, we evaluated the degree of DSB complexity induced after photon irradiation for various X-ray spectra using the foci size, and the experimental DSB complexity was compared to the results estimated by the well-verified DNA damage estimation model in the PHITS code. The number of events in a water cube was found to be proportional to foci size, suggesting that the number of events intrinsically related to DSB complexity at the DNA scale. The developed method was applicable to focus data measured for various X-ray spectral situations (i.e., diagnostic kV X-rays and therapeutic MV X-rays). This method would contribute to a precise understanding of the early biological impacts of photon irradiation by means of the γ-H2AX focus formation assay.

A Simplified Cluster Analysis of Electron Track Structure for Estimating Complex DNA Damage Yields.

Complex DNA damage, defined as at least two vicinal lesions within 10-20 base pairs (bp), induced after exposure to ionizing radiation, is recognized as fatal damage to human tissue. Due to the difficulty of directly measuring the aggregation of DNA damage at the nano-meter scale, many cluster analyses of inelastic interactions based on Monte Carlo simulation for radiation track structure in liquid water have been conducted to evaluate DNA damage. Meanwhile, the experimental technique to detect complex DNA damage has evolved in recent decades, so both approaches with simulation and experiment get used for investigating complex DNA damage. During this study, we propose a simplified cluster analysis of ionization and electronic excitation events within 10 bp based on track structure for estimating complex DNA damage yields for electron and X-ray irradiations. We then compare the computational results with the experimental complex DNA damage coupled with base damage (BD) measured by enzymatic cleavage and atomic force microscopy (AFM). The computational results agree well with experimental fractions of complex damage yields, i.e., single and double strand breaks (SSBs, DSBs) and complex BD, when the yield ratio of BD/SSB is assumed to be 1.3. Considering the comparison of complex DSB yields, i.e., DSB + BD and DSB + 2BD, between simulation and experimental data, we find that the aggregation degree of the events along electron tracks reflects the complexity of induced DNA damage, showing 43.5% of DSB induced after 70 kVp X-ray irradiation can be classified as a complex form coupled with BD. The present simulation enables us to quantify the type of complex damage which cannot be measured through in vitro experiments and helps us to interpret the experimental detection efficiency for complex BD measured by AFM. This simple model for estimating complex DNA damage yields contributes to the precise understanding of the DNA damage complexity induced after X-ray and electron irradiations.

Inorganic polyphosphate enhances radio-sensitivity in a human non-small cell lung cancer cell line, H1299.

Inorganic polyphosphate is a linear polymer containing tens to hundreds of orthophosphate residues linked by high-energy phosphoanhydride bonds. Polyphosphate has been recognized as a potent anti-metastasis reagent. However, the molecular mechanism underlying polyphosphate action on cancer cells is poorly understood. In this study, we investigated the involvement of polyphosphate in radio-sensitivity using a human non-small cell lung cancer cell line, H1299. We found that polyphosphate treatment decreases cellular adenosine triphosphate levels, suggesting a disruption of energy metabolism. We also found that the induction of DNA double-strand breaks was enhanced in polyphosphate-treated cells after X-ray irradiation and colony formation assay revealed that cell survival decreased compared with that of the control groups. These findings suggest that polyphosphate is a promising radio-sensitizer for cancer cells. Therefore, we hypothesized that polyphosphate treatment disrupts adenosine triphosphate-mediated energy transfer for cellular survival and DNA repair, thereby reducing the cellular capability to resist X-ray irradiation.

Evaluation of relative biological effectiveness for diseases of the circulatory system based on microdosimetry.

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.

Possible mechanisms and simulation modeling of FLASH radiotherapy.

FLASH radiotherapy (FLASH-RT) has great potential to improve patient outcomes. It delivers radiation doses at an ultra-high dose rate (UHDR: ≥ 40 Gy/s) in a single instant or a few pulses. Much higher irradiation doses can be administered to tumors with FLASH-RT than with conventional dose rate (0.01-0.40 Gy/s) radiotherapy. UHDR irradiation can suppress toxicity in normal tissues while sustaining antitumor efficiency, which is referred to as the FLASH effect. However, the mechanisms underlying the effects of the FLASH remain unclear. To clarify these mechanisms, the development of simulation models that can contribute to treatment planning for FLASH-RT is still underway. Previous studies indicated that transient oxygen depletion or augmented reactions between secondary reactive species produced by irradiation may be involved in this process. To discuss the possible mechanisms of the FLASH effect and its clinical potential, we summarized the physicochemical, chemical, and biological perspectives as well as the development of simulation modeling for FLASH-RT.

Reply to comment on 'Modeling for predicting survival fraction of cells after ultra-high dose rate irradiation'.

Liew and Mairani commented on our paper 'Modeling for predicting survival fraction of cells after ultra-high dose rate irradiation' (Shiraishiet al2024aPhys. Med. Biol.69015017), which proposed a biophysical model to predict the dose-response curve of surviving cell fractions after ultra-high dose rate irradiation following conventional dose rate irradiation by considering DNA damage yields. They suggested the need to consider oxygen concentration in our prediction model and possible issues related to the data selection process used for the benchmarking test in our paper. In this reply, we discuss the limitations of both the present model and the available experimental data for determining the model's parameters. We also demonstrate that our proposed model can reproduce the experimental survival data even when using only the experimental DNA damage data measured reliably under normoxic conditions.

The impact of dose rate on responses of human lens epithelial cells to ionizing irradiation.

The knowledge on responses of human lens epithelial cells (HLECs) to ionizing radiation exposure is important to understand mechanisms of radiation cataracts that are of concern in the field of radiation protection and radiation therapy. However, biological effects in HLECs following protracted exposure have not yet fully been explored. Here, we investigated the temporal kinetics of γ-H2AX foci as a marker for DNA double-strand breaks (DSBs) and cell survival in HLECs after exposure to photon beams at various dose rates (i.e., 150 kVp X-rays at 1.82, 0.1, and 0.033 Gy/min, and 137Cs γ-rays at 0.00461 Gy/min (27.7 cGy/h) and 0.00081 Gy/min (4.9 cGy/h)), compared to those in human lung fibroblasts (WI-38). In parallel, we quantified the recovery for DSBs and cell survival using a biophysical model. The study revealed that HLECs have a lower DSB repair rate than WI-38 cells. There is no significant impact of dose rate on cell survival in both cell lines in the dose-rate range of 0.033-1.82 Gy/min. In contrast, the experimental residual γ-H2AX foci showed inverse dose rate effects (IDREs) compared to the model prediction, highlighting the importance of the IDREs in evaluating radiation effects on the ocular lens.

In vitro and in silico study of biological effects on cancer cells in the presence of metallic materials during radiotherapy.

X-ray therapy aims to eliminate tumours while minimizing side effects. Intense mucositis is sometimes induced when irradiating the oral cavity with a dental metal crown (DMC). However, the underlying mechanisms of such inducing radiosensitization by DMC remain uncertain. This study explored the radiosensitizing mechanisms around DMCs in an interdisciplinary approach with cell experiments and Monte Carlo simulation with the PHITS code. Clonogenic survival and nuclear 53BP1 foci of a cell line derived from cervical cancer cells (HeLa cells) were measured post-irradiation with therapeutic X-rays near high-Z materials such as Pb or Au plates, and the experimental sensitizer enhancement ratio (SER) was obtained. Meanwhile, the dose enhancement ratio (DER) and relative biological effectiveness for DNA damage yields were calculated using the PHITS code, by considering the corresponding experimental condition. The experiments show the experimental SER values for cell survival and 53BP1 foci near metals are 1.2-1.4, which agrees well with the calculated DER values. These suggest that the radiosensitizing effects near metal are predominantly attributed to the dose increase. In addition, as a preclinical evaluation, the spatial distributions of DER near DMC are calculated using Computed Tomography Digital Imaging and Communications in Medicine (CT-DICOM) data and a simple tooth model. As a result, the DER values evaluated using the CT-DICOM data were lower than those from a simple tooth model. These findings highlight the challenge of evaluating radiosensitizing effects near DMCs using Digital Imaging and Communications in Medicine (DICOM) images due to volume-averaging effects and emphasize the need for a high-resolution (<1 mm) dose assessment method unaffected by these effects.

A step-by-step simulation code for estimating yields of water radiolysis species based on electron track-structure mode in the PHITS code.

Objective. Time-dependent yields of chemical products resulting from water radiolysis play a great role in evaluating DNA damage response after exposure to ionizing radiation. Particle and Heavy Ion Transport code System (PHITS) is a general-purpose Monte Carlo simulation code for radiation transport, which simulates atomic interactions originating from discrete energy levels of ionizations and electronic excitations as well as molecular excitations as physical stages. However, no chemical code for simulating water radiolysis products exists in the PHITS package.Approach.Here, we developed a chemical simulation code dedicated to the PHITS code, hereafter calledPHITS-Chemcode, which enables the calculation of theGvalues of water radiolysis species (•OH, eaq-, H2, H2O2etc) by electron beams.Main results.The estimatedGvalues during 1 µs are in agreement with the experimental ones and other simulations. ThisPHITS-Chemcode also simulates the radiolysis in the presence of OH radical scavengers, such as tris(hydroxymethyl)aminomethane and dimethyl sulfoxide. Thank to this feature, the contributions of direct and indirect effects on DNA damage induction under various scavenging capacities can be analyzed.Significance.This chemical code coupled with PHITS could contribute to elucidating the mechanism of radiation effects by connecting physical, physicochemical, and chemical processes.

DNA damage response in a 2D-culture model by diffusing alpha-emitters radiation therapy (Alpha-DaRT).

Diffusing alpha-emitters radiation therapy (Alpha-DaRT) is a unique method, in which interstitial sources carrying 224Ra release a chain of short-lived daughter atoms from their surface. Although DNA damage response (DDR) is crucial to inducing cell death after irradiation, how the DDR occurs during Alpha-DaRT treatment has not yet been explored. In this study, we temporo-spatially characterized DDR such as kinetics of DNA double-strand breaks (DSBs) and cell cycle, in two-dimensional (2D) culture conditions qualitatively mimicking Alpha-DaRT treatments, by employing HeLa cells expressing the Fucci cell cycle-visualizing system. The distribution of the alpha-particle pits detected by a plastic nuclear track detector, CR-39, strongly correlated with γH2AX staining, a marker of DSBs, around the 224Ra source, but the area of G2 arrested cells was more widely spread 24 h from the start of the exposure. Thereafter, close time-lapse observation revealed varying cell cycle kinetics, depending on the distance from the source. A medium containing daughter nuclides prepared from 224Ra sources allowed us to estimate the radiation dose after 24 h of exposure, and determine surviving fractions. The present experimental model revealed for the first time temporo-spatial information of DDR occurring around the source in its early stages.

Generality Assessment of a Model Considering Heterogeneous Cancer Cells for Predicting Tumor Control Probability for Stereotactic Body Radiation Therapy Against Non-Small Cell Lung Cancer.

The generality of a model for predicting tumor control probability from in vitro clonogenic survival considering of cancer stem-like cells, the so-called integrated microdosimetric-kinetic model, is presented by comparing the model to public data on stereotactic body radiation therapy for non-small cell lung cancer cells.

Cell-cycle dependence on the biological effects of boron neutron capture therapy and its modification by polyvinyl alcohol.

Boron neutron capture therapy (BNCT) is a unique radiotherapy of selectively eradicating tumor cells using boron compounds (e.g., 4-borono-L-phenylalanine [BPA]) that are heterogeneously taken up at the cellular level. Such heterogenicity potentially reduces the curative efficiency. However, the effects of temporospatial heterogenicity on cell killing remain unclear. With the technical combination of radiation track detector and biophysical simulations, this study revealed the cell cycle-dependent heterogenicity of BPA uptake and subsequent biological effects of BNCT on HeLa cells expressing fluorescent ubiquitination-based cell cycle indicators, as well as the modification effects of polyvinyl alcohol (PVA). The results showed that the BPA concentration in the S/G2/M phase was higher than that in the G1/S phase and that PVA enhances the biological effects both by improving the uptake and by canceling the heterogenicity. These findings might contribute to a maximization of therapeutic efficacy when BNCT is combined with PVA and/or cell cycle-specific anticancer agents.

Radiation research trends by young scientists and the future tasks in Northern Japan: report on 'the 10th educational symposium on radiation and health (ESRAH) by young scientists in 2023'.

PURPOSE: Since 2014, an educational activity on radiation and health in northern Japan has been carried out by young scientists, the so-called 'Educational Symposium on Radiation and Health (ESRAH)'. Close cooperation has been continued in preparing for any possible emergency response to radiation accidents because several facilities, e.g., the Tomari Nuclear Power Plant in Hokkaido and the Low-Level Radioactive Waste Disposal Facility in Aomori prefecture. The ESRAH meeting has provided informational exchange and discussion forum on a broad range of subjects in various. In 2023, the 10th Memorial ESRAH meeting took place to boost scientific understanding and multidisciplinary collaborations for young scientists. Herein, we report on the ESRAH2023 symposium and analyze the research categories of young scientists from the past 10-year presentations. CONCLUSIONS: To date, the ESRAH meeting has successfully provided a chance for multi-disciplinary research, which accounted for 27% of the total despite the COVID-19 pandemic. We found that the fraction of multi-disciplinary research in 2023 was the highest during 10-year ESRAH meetings. Meanwhile, amongst the research categories, physics, chemistry, and pharmacological studies were indicated to be less for young scientists. It is desired that further collaboration between physics, chemistry, and pharmacology in addition to the current fields would not only clarify radiation effects on the human body but also promote emergency medical care for radiation exposure in the future.

Modeling for predicting survival fraction of cells after ultra-high dose rate irradiation.

Objective. FLASH radiotherapy (FLASH-RT) with ultra-high dose rate (UHDR) irradiation (i.e. > 40 Gy s-1) spares the function of normal tissues while preserving antitumor efficacy, known as the FLASH effect. The biological effects after conventional dose rate-radiotherapy (CONV-RT) with ≤0.1 Gy s-1have been well modeled by considering microdosimetry and DNA repair processes, meanwhile modeling of radiosensitivities under UHDR irradiation is insufficient. Here, we developed anintegrated microdosimetric-kinetic(IMK)model for UHDR-irradiationenabling the prediction of surviving fraction after UHDR irradiation.Approach.TheIMK model for UHDR-irradiationconsiders the initial DNA damage yields by the modification of indirect effects under UHDR compared to CONV dose rate. The developed model is based on the linear-quadratic (LQ) nature with the dose and dose square coefficients, considering the reduction of DNA damage yields as a function of dose rate.Main results.The estimate by the developed model could successfully reproduce thein vitroexperimental dose-response curve for various cell line types and dose rates.Significance.The developed model would be useful for predicting the biological effects under the UHDR irradiation.

First-principles simulation of an ejected electron produced by monochromatic deposition energy to water at the femtosecond order.

This study uses a time-dependent first-principles simulation code to investigate the transient dynamics of an ejected electron produced in the monochromatic deposition energy from 11 to 19 eV in water. The energy deposition forms a three-body single spur comprising a hydroxyl radical (OH˙), hydronium ion (H3O+), and hydrated electron (eaq-). The earliest formation involves electron thermalization and delocalization dominated by the molecular excitation of water. Our simulation results show that the transient electron dynamics primarily depends on the amount of deposition energy to water; the thermalization time varies from 200 to 500 fs, and the delocalization varies from 3 to 10 nm in this energy range. These features are crucial for determining the earliest single-spur formation and facilitating a sequential simulation from an energy deposition to a chemical reaction in water photolysis or radiolysis. The spur radius obtained from the simulation correlates reasonably with the experimental-based estimations. Our results should provide universalistic insights for analysing ultrafast phenomena dominated by the molecular excitation of water in the femtosecond order.

Improvement of the hybrid approach between Monte Carlo simulation and analytical function for calculating microdosimetric probability densities in macroscopic matter.

Objective. Estimation of the probability density of the microdosimetric quantities in macroscopic matter is indispensable for applying the concept of microdosimetry to medical physics and radiological protection. The Particle and Heavy Ion Transport code System (PHITS) enables estimating the microdosimetric probability densities due to its unique hybrid modality between the Monte Carlo and analytical approaches called the microdosimetric function. It can convert the deposition energies calculated by the macroscopic Monte Carlo radiation transport simulation to microdosimetric probability densities in water using an analytical function based on the track-structure simulations.Approach. In this study, we improved this function using the latest track-structure simulation codes implemented in PHITS. The improved function is capable of calculating the probability densities of not only the conventional microdosimetric quantities such as lineal energy but also the number of ionization events occurring in a target site, the so-called ionization cluster size distribution, for arbitrary site diameters from 3 nm to 1µm.Main results. The accuracy of the improved function was well verified by comparing the microdosimetric probability densities measured by tissue-equivalent proportional counters with the corresponding data calculated in this study. Test calculations for clonogenic cell survival using the improved function coupled with the modified microdosimetric kinetic model suggested a slight increase of its relative biological effectiveness compared with our previous estimations. As a new application of the improved function, we calculated the relative biological effectiveness of the single-strand break and double-strand break yields for proton irradiations using the updated PHITS coupled with the simplified DNA damage estimation model, and confirmed its equivalence in accuracy and its superiority in computational time compared to our previously proposed method based on the track-structure simulation.Significance. From these features, we concluded that the improved function could expand the application fields of PHITS by bridging the gap between microdosimetry and macrodosimetry.

Initial yield of hydrated electron production from water radiolysis based on first-principles calculation.

Many scientific insights into water radiolysis have been applied for developing life science, including radiation-induced phenomena, such as DNA damage and mutation induction or carcinogenesis. However, the generation mechanism of free radicals due to radiolysis remains to be fully understood. Consequently, we have encountered a crucial problem in that the initial yields connecting radiation physics to chemistry must be parameterized. We have been challenged in the development of a simulation tool that can unravel the initial free radical yields, from physical interaction by radiation. The presented code enables the first-principles calculation of low energy secondary electrons resulting from the ionization, in which the secondary electron dynamics are simulated while considering dominant collision and polarization effects in water. In this study, using this code, we predicted the yield ratio between ionization and electronic excitation from a delocalization distribution of secondary electrons. The simulation result presented a theoretical initial yield of hydrated electrons. In radiation physics, the initial yield predicted from parameter analysis of radiolysis experiments in radiation chemistry was successfully reproduced. Our simulation code helps realize a reasonable spatiotemporal connection from radiation physics to chemistry, which would contribute to providing new scientific insights for precise understanding of underlying mechanisms of DNA damage induction.

Modelling oxygen effects on the in- and out-of-field radiosensitivity of cells exposed to intensity-modulated radiation fields.

Objective. The delivery of intensity-modulated radiation fields has improved the conformity of dose to tumour targets during radiotherapy (RT). Previously, it has been shown that intercellular communication between cells positioned in- and outside of the radiation field impacts cellular radiosensitivity under hypoxic and normoxic conditions. However, the mechanism of intercellular communication in hypoxia remains to be fully understood. In this study, the cell-killing effects of intercellular communication in hypoxia were modelled in an effort to better understand the underlying mechanisms of response.Approach. By irradiating a 50% area of the culture dish (half-field exposure), experimental dose-response curves for cell survival and residual DNA double-strand breaks (DSBs) were generated in prostate (DU145) and non-small cell lung cancer (H1299) cells. The oxygen enhancement ratio (OER) was determined from early DSB yields (corresponding to relative direct damage) and used to model the in- and out-of-field radiosensitivity.Main results. The developed integrated microdosimetric-kinetic (IMK) model successfully predicted the experimental dose responses for survival and lethal lesions, and provides a mechanistic interpretation that the probability of hits for releasing cell-killing signals is dependent on oxygen. This experimental and modelling study also suggests that residual DSBs correspond to logarithmic survival fraction (meaning lethal lesions) for in- and out-of-field cells. Our data suggest that the OER value determined using uniform-field exposure can be applied to predict the in- and out-of-field radiosensitivity of cells following exposure to intensity modulated beams.Significance. The developed IMK model facilitates a more precise understanding of intercellular signalling following exposure to intensity-modulated radiation fields.

Translational study for stereotactic body radiotherapy against non-small cell lung cancer, including oligometastases, considering cancer stem-like cells enable predicting clinical outcome from in vitro data.

BACKGROUND: Curative effects of stereotactic body radiotherapy (SBRT) for non-small cell lung cancer (NSCLC) have been evaluated using various biophysical models. Because such model parameters are empirically determined based on clinical experience, there is a large gap between in vitro and clinical studies. In this study, considering the heterogeneous cell population, we performed a translational study to realize the possible linkage based on a modeling approach. METHODS: We modeled cell-killing and tumor control probability (TCP) considering two populations: progeny and cancer stem-like cells. The model parameters were determined from in vitro survival data of A549 and EBC-1 cells. Based on the cellular parameters, we predicted TCP and compared it with the corresponding clinical data from 553 patients collected at Hirosaki University Hospital. RESULTS: Using an all-in-one developed model, the so-called integrated microdosimetric-kinetic (IMK) model, we successfully reproduced both in vitro survival after acute irradiation and the 3-year TCP with various fractionation schemes (6-10 Gy per fraction). From the conventional prediction without considering cancer stem cells (CSCs), this study revealed that radioresistant CSCs play a key role in the linkage between in vitro and clinical outcomes. CONCLUSIONS: This modeling study provides a possible generalized biophysical model that enables precise estimation of SBRT worldwide.