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.

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.

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.

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.

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.

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.

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.

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.

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.

Application of a simple DNA damage model developed for electrons to proton irradiation.

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keVµm-1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Impact of the Lorentz force on electron track structure and early DNA damage yields in magnetic resonance-guided radiotherapy.

Magnetic resonance-guided radiotherapy (MRgRT) has been developed and installed in recent decades for external radiotherapy in several clinical facilities. Lorentz forces modulate dose distribution by charged particles in MRgRT; however, the impact of Lorentz forces on low-energy electron track structure and early DNA damage induction remain unclear. In this study, we estimated features of electron track structure and biological effects in a static magnetic field (SMF) using a general-purpose Monte Carlo code, particle and heavy ion transport code system (PHITS) that enables us to simulate low-energy electrons down to 1 meV by track-structure mode. The macroscopic dose distributions by electrons above approximately 300 keV initial energy in liquid water are changed by both perpendicular and parallel SMFs against the incident direction, indicating that the Lorentz force plays an important role in calculating dose within tumours. Meanwhile, DNA damage estimation based on the spatial patterns of atomic interactions indicates that the initial yield of DNA double-strand breaks (DSBs) is independent of the SMF intensity. The DSB induction is predominantly attributed to the secondary electrons below a few tens of eV, of which energy deposition patterns are not considerably affected by the Lorentz force. Our simulation study suggests that treatment planning for MRgRT can be made with consideration of only changed dose distribution.

Microdosimetric Modeling of Relative Biological Effectiveness for Skin Reactions: Possible Linkage Between In Vitro and In Vivo Data.

PURPOSE: Precise evaluation of the relative biological effectiveness (RBE) for skin reactions is indispensable for the treatment planning of particle therapy and boron-neutron capture therapy. The evaluation is also needed for the future radiation protection system for mixed radiation fields. Such RBE is often evaluated based on in vitro cell survival data, but its validity remains incompletely understood. This study aimed to develop a model for estimating RBE for skin reactions and dermal cell survival in the same framework and quantitatively discuss a possible link between them. METHODS AND MATERIALS: The microdosimetric kinetic model, which was originally developed for estimating cell surviving fractions for various radiations, was improved to be capable of estimating the mean and uncertainty of RBE for skin reactions. The parameter used in the model was independently determined from in vitro measurements of dermal cell survival and in vivo measurements of skin reactions taken from 8 and 23 articles, respectively. In the parameter determination, the characteristics of the radiation fields employed in each measurement were reproduced in detail by the Particle and Heavy Ion Transport Code System. RESULTS: Our model quantitatively revealed that RBE for skin reactions tend to be higher than those for dermal cell survival. RBE of various monoenergetic radiations calculated from this model confirmed that the past evaluations made by the International Commission on Radiological Protection and the National Council on Radiation Protection and Measurements a few decades ago are still supported by recent experimental data. CONCLUSIONS: Our model can play important roles not only in medical physics for avoiding unnecessary skin reactions in particle therapy and boron-neutron capture therapy but also in radiation protection for future decision-making of the recommended RBE values.

Inflammatory Signaling and DNA Damage Responses after Local Exposure to an Insoluble Radioactive Microparticle.

Cesium-bearing microparticles (Cs-BMPs) can reach the human respiratory system after inhalation, resulting in chronic local internal exposure. We previously investigated the spatial distribution of DNA damage induced in areas around a Cs-BMP; however, the biological impacts have not been fully clarified due to the limited amount of data. Here, we investigated the inflammatory signaling and DNA damage responses after local exposure to a Cs-BMP in vitro. We used two normal human lung cell lines, i.e., lung fibroblast cells (WI-38) and bronchial epithelial cells (HBEC3-KT). After 24 h exposure to a Cs-BMP, inflammation was evaluated by immunofluorescent staining for nuclear factor κB (NF-κB) p65 and cyclooxygenase 2 (COX-2). The number of DNA double-strand breaks (DSBs) was also detected by means of phospholylated histone H2AX (γ-H2AX) focus formation assay. Cs-BMP exposure significantly increased NF-κB p65 and COX-2 expressions, which were related to the number of γ-H2AX foci in the cell nuclei. Compared to the uniform (external) exposure to 137Cs γ-rays, NF-κB tended to be more activated in the cells proximal to the Cs-BMP, while both NF-κB p65 and COX-2 were significantly activated in the distal cells. Experiments with chemical inhibitors for NF-κB p65 and COX-2 suggested the involvement of such inflammatory responses both in the reduced radiosensitivity of the cells proximal to Cs-BMP and the enhanced radiosensitivity of the cells distal from Cs-BMP. The data show that local exposure to Cs-BMP leads to biological effects modified by the NF-κB pathway, suggesting that the radiation risk for Cs-BMP exposure can differ from that estimated based on conventional uniform exposure to normal tissues.

Tumor radioresistance caused by radiation-induced changes of stem-like cell content and sub-lethal damage repair capability.

Cancer stem-like cells (CSCs) within solid tumors exhibit radioresistance, leading to recurrence and distant metastasis after radiotherapy. To experimentally study the characteristics of CSCs, radioresistant cell lines were successfully established using fractionated X-ray irradiation. The fundamental characteristics of CSCs in vitro have been previously reported; however, the relationship between CSC and acquired radioresistance remains uncertain. To efficiently study this relationship, we performed both in vitro experiments and theoretical analysis using a cell-killing model. Four types of human oral squamous carcinoma cell lines, non-radioresistant cell lines (SAS and HSC2), and radioresistant cell lines (SAS-R and HSC2-R), were used to measure the surviving fraction after single-dose irradiation, split-dose irradiation, and multi-fractionated irradiation. The SAS-R and HSC2-R cell lines were more positive for one of the CSC marker aldehyde dehydrogenase activity than the corresponding non-radioresistant cell lines. The theoretical model analysis showed that changes in both the experimental-based ALDH (+) fractions and DNA repair efficiency of ALDH (-) fractions (i.e., sub-lethal damage repair) are required to reproduce the measured cell survival data of non-radioresistant and radioresistant cell lines. These results suggest that the enhanced cell recovery in SAS-R and HSC2-R is important when predicting tumor control probability in radiotherapy to require a long dose-delivery time; in other words, intensity-modulated radiation therapy is ideal. This work provides a precise understanding of the mechanism of radioresistance, which is induced after irradiation of cancer cells.

Track-structure modes in particle and heavy ion transport code system (PHITS): application to radiobiological research.

PURPOSE: In radiation physics, Monte Carlo radiation transport simulations are powerful tools to evaluate the cellular responses after irradiation. When investigating such radiation-induced biological effects, it is essential to perform track structure simulations by explicitly considering each atomic interaction in liquid water at the sub-cellular and DNA scales. The Particle and Heavy-Ion Transport code System (PHITS) is a Monte Carlo code which enables to calculate track structure at DNA scale by employing the track-structure modes for electrons, protons and carbon ions. In this paper, we review the recent development status and future prospects of the track-structure modes in the PHITS code. CONCLUSIONS: To date, the physical features of these modes have been verified using the available experimental data and Monte Carlo simulation results reported in literature. These track-structure modes can be used for calculating microdosimetric distributions to estimate cell survival and for estimating initial DNA damage yields. The use of PHITS track-structure mode is expected not only to clarify the underlying mechanisms of radiation effects but also to predict curative effects in radiation therapy. The results of PHITS simulations coupled with biophysical models will contribute to the radiobiological studies by precisely predicting radiation-induced biological effects based on the Monte Carlo approach.