Modeling cell survival and change in amount of DNA during protracted irradiation

Effect of Ultrahigh Dose Rate on Biomolecular Radiation Damage

Dose rate is one of the important parameters in radiation-induced biomolecular damage. The effects of dose rate have been known to modify radiation toxicity in biological systems. The rate and extent of sublethal DNA damage (e.g., base damage and single-strand breaks) repair and those of cell proliferation have been manifested by dose rate. However, the recent preclinical application of ultrahigh dose rate [(UHDR) ca. 40 Gy/s and higher] radiation modalities have been shown to lower the type and extent of radiation damage to biological systems. At these UHDR, radiation-induced physicochemical and chemical processes are expected to differ from those observed after irradiation at conventional dose rates (CONV). It is unclear whether these UHDR conditions can affect the quality (type) and quantity (extent) of biomolecular damage such as DNA lesions. Here, we comparatively study the influence of indirect effects of CONV and UHDR on the formation of DNA strand breaks and clustered damage including densely accumulated lesions in an aerated and an anoxic dilute aqueous solution of a plasmid DNA model under low and high hydroxyl radical (•OH) scavenging conditions. Aqueous solutions of purified supercoiled plasmid DNA (pUC19) were prepared in either air- or nitrogen-saturated conditions, with Tris buffer added as the radiation-produced •OH scavenger at low and high scavenging capacities. These DNA samples were irradiated using kV X-ray systems at CONV (0.1 Gy/s) and high dose rate (HDR, 25 Gy/s) as well as UHDR (55 and 125 Gy/s) under different scavenging and environmental conditions. DNA lesions including strand breaks and clustered damage including densely accumulated lesions were quantified by gel electrophoresis and the yields of these lesions were calculated from the dose-response curve. Non-DSB clustered damage including densely accumulated lesions were evaluated by treating DNAs using bacterial endonuclease enzymes (Fpg and Nth) prior to gel electrophoresis. UHDR of 55 and 125 Gy/s induced lower amounts of both isolated strand breaks and clustered DNA damage including densely accumulated lesions at doses >40 Gy in the presence of oxygen, compared to the abundance of these lesions induced by 0.1 and 25 Gy/s irradiation under the same dose conditions. Overall, the strand break and clustered damage including densely accumulated lesions yields decreased by factors of 1.3-3.5 after UHDR. We did not observe these differences either via •OH scavenging or by removing oxygen from the solution. In addition, our results point out that the inter-track recombination reactions did not contribute to the observed dose-rate effects on DNA damage. The effects of dose rate on DNA damage are highly dependent on the total dose, as expected, but also on the •OH scavenging capacity that is employed in the aqueous DNA solutions. These important variables may be relevant in biological systems as well. On a practical level, our in vitro plasmid DNA model, which permits to precisely vary the •OH scavenging capacity and gassing conditions (air saturated vs. N2 saturated) can help to differentiate dose-rate effects on biomolecular damage. Our results indicate that the radical-radical reactions are important in understanding the dose-rate effect on DNA damage.

Assessing tumor volumetric reduction with consideration for setup errors based on mathematical tumor model and microdosimetric kinetic model in single-isocenter VMAT for brain metastases

The volumetric reduction rate (VRR) was evaluated with consideration for six degrees-of-freedom (6DoF) patient setup errors based on a mathematical tumor model in single-isocenter volumetric modulated arc therapy (SI-VMAT) for brain metastases. Simulated gross tumor volumes (GTV) of 1.0 cm and dose distribution were created (27 Gy/3 fractions). The distance between the GTV center and isocenter (d) was set at 0-10 cm. The GTV was translated within 0-1.0 mm (Trans) and rotated within 0-1.0° (Rot) in the three axis directions using affine transformation. The tumor growth volume was calculated using a multicomponent mathematical model (MCTM), and lethal effects of irradiation and repair from damage during irradiation were calculated by a microdosimetric kinetic model (MKM) for non-small cell lung cancer (NSCLC) A549 and NCI-H460 (H460) cells. The VRRs were calculated 5 days after the end of irradiation using the physical dose to the GTV for varying d and 6DoF setup errors. The tolerance value of VRR, the GTV volume reduction rate, was set at 5%, based on the pre-irradiation GTV volume. With the exception of the only one A549 condition where (Trans, Rot) = (1.0 mm, 1.0°) was repeated for 3 fractions, all conditions met all the tolerance VRR values for A549 and H460 cells with varying d from 0 to 10 cm. Evaluation based on the mathematical tumor model suggested that if the 6DoF setup errors at each irradiation could be kept within 1.0 mm and 1.0°, there would be little effect on tumor volume regardless of the distance from the isocenter in SI-VMAT.

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.

Dose Rate Effect on Cell Survival in BNCT

The output constancy of the accelerator used for boron neutron capture therapy (BNCT) is essential to ensuring anti-tumor efficacy and safety. BNCT as currently practiced requires a wide variety of beam quality assessments to ensure that RBE dose errors are maintained within 5%. However, the necessity of maintaining a constant beam dose rate has not been fully discussed. We therefore clarified the effect of different physical dose rates of the accelerator BNCT on biological effects. SAS and A172 cells exposed to 10B-boronophenylalanine were irradiated using a neutron beam (normal operating current, 100 μA) at the Aomori Quantum Science Center. Thermal neutron flux was attenuated to 50.0 ± 0.96% under 50 μA irradiation compared to that under 100 μA irradiation. Cells were given physical doses of 1.67 and 3.36 Gy at 30 and 60 mC, respectively, and survival was significantly increased after 50 μA irradiation for both cell types (p = 0.0052 for SAS; p = 0.046 for A172, for 60 mC). Differences in accelerator BNCT beam dose rates have non-negligible effects on biological effects. Dose rate fluctuations and differences should not be easily permitted to obtain consistent biological effects.

Mathematical model combined with microdosimetric kinetic model for tumor volume calculation in stereotactic body radiation therapy

We proposed a new mathematical model that combines an ordinary differential equation (ODE) and microdosimetric kinetic model (MKM) to predict the tumor-cell lethal effect of Stereotactic body radiation therapy (SBRT) applied to non-small cell lung cancer (NSCLC). The tumor growth volume was calculated by the ODE in the multi-component mathematical model (MCM) for the cell lines NSCLC A549 and NCI-H460 (H460). The prescription doses 48 Gy/4 fr and 54 Gy/3 fr were used in the SBRT, and the effect of the SBRT on tumor cells was evaluated by the MKM. We also evaluated the effects of (1) linear quadratic model (LQM) and the MKM, (2) varying the ratio of active and quiescent tumors for the total tumor volume, and (3) the length of the dose-delivery time per fractionated dose (tinter) on the initial tumor volume. We used the ratio of the tumor volume at 1 day after the end of irradiation to the tumor volume before irradiation to define the radiation effectiveness value (REV). The combination of MKM and MCM significantly reduced REV at 48 Gy/4 fr compared to the combination of LQM and MCM. The ratio of active tumors and the prolonging of tinter affected the decrease in the REV for A549 and H460 cells. We evaluated the tumor volume considering a large fractionated dose and the dose-delivery time by combining the MKM with a mathematical model of tumor growth using an ODE in lung SBRT for NSCLC A549 and H460 cells.

Radiobiological evaluation considering the treatment time with stereotactic radiosurgery for brain metastases

Objective

We evaluated the radiobiological effect of the irradiation time with the interruption time of stereotactic radiosurgery (SRS) using CyberKnife® (CK) systemfor brain metastases.

Methods

We used the DICOM data and irradiation log file of the 10 patients with brain metastases from non-small-cell lung cancer (NSCLC) who underwent brain SRS. We defined the treatment time as the sum of the dose-delivery time and the interruption time during irradiations, and we used a microdosimetric kinetic model (MKM) to evaluate the radiobiological effects of the treatment time. The biological parameters, i.e. α0, β0, and the DNA repair constant rate (a + c), were acquired from NCI-H460 cell for the MKM. We calculated the radiobiological dose for the gross tumor volume (GTVbio) to evaluate the treatment time's effect compared with no treatment time as a reference. The D95 (%) and the Radiation Therapy Oncology Group conformity index (RCI) and Paddick conformity index (PCI) were calculated as dosimetric indices. We used several DNA repair constant rates (a + c) (0.46, 1.0, and 2.0) to assess the radiobiological effect by varying the DNA repair date (a + c) values.

Results

The mean values of D95 (%), RCI, and PCI for GTVbio were 98.8%, 0.90, and 0.80, respectively, and decreased with increasing treatment time. The mean values of D95 (%), RCI, and PCI of GTVbio at 2.0 (a+c) value were 94.9%, 0.71, and 0.49, respectively.

Conclusion

The radiobiological effect of the treatment time on tumors was accurately evaluated with brain SRS using CK.

Advances in knowledge

There has been no published investigation of the radiobiological impact of the longer treatment time with multiple interruptions of SRS using a CK on the target dose distribution in a comparison with the use of a linac. Radiobiological dose assessment that takes into account treatment time in the physical dose in this study may allow more accurate dose assessment in SRS for metastatic brain tumors using CK.

Radiation dose rate effects: what is new and what is needed?

Despite decades of research to understand the biological effects of ionising radiation, there is still much uncertainty over the role of dose rate. Motivated by a virtual workshop on the "Effects of spatial and temporal variation in dose delivery" organised in November 2020 by the Multidisciplinary Low Dose Initiative (MELODI), here, we review studies to date exploring dose rate effects, highlighting significant findings, recent advances and to provide perspective and recommendations for requirements and direction of future work. A comprehensive range of studies is considered, including molecular, cellular, animal, and human studies, with a focus on low linear-energy-transfer radiation exposure. Limits and advantages of each type of study are discussed, and a focus is made on future research needs.

Radiobiological effects of the interruption time with Monte Carlo Simulation on multiple fields in photon beams

PURPOSE

The interruption time is the irradiation interruption that occurs at sites and operations such as the gantry, collimator, couch rotation, and patient setup within the field in radiotherapy. However, the radiobiological effect of prolonging the treatment time by the interruption time for tumor cells is little evaluated. We investigated the effect of the interruption time on the radiobiological effectiveness with photon beams based on a modified microdosimetric kinetic (mMK) model.

METHODS

The dose-mean lineal energy y_D (keV/µm) of 6-MV photon beams was calculated by the particle and heavy ion transport system (PHITS). We set the absorbed dose to 2 or 8 Gy, and the interruption time (τ) was set to 1, 3, 5, 10, 30, and 60 min. The biological parameters such as α_0, β_0, and DNA repair constant rate (a + c) values were acquired from a human non-small-cell lung cancer cell line (NCI-H460) for the mMK model. We used two-field and four-field irradiation with a constant dose rate (3 Gy/min); the photon beams were paused for interruption time τ. We calculated the relative biological effectiveness (RBE) to evaluate the interruption time's effect compared with no interrupted as a reference.

RESULTS

The y_D of 6-MV photon beams was 2.32 (keV/µm), and there was little effect by changing the water depth (standard deviation was 0.01). The RBE with four-field irradiation for 8 Gy was decreased to 0.997, 0.975, 0.900, and 0.836 τ = 1, 10, 30, 60 min, respectively. In addition, the RBE was affected by the repair constant rate (a + c) value, the greater the decrease in RBE with the longer the interruption time when the (a + c) value was large.

CONCLUSION

The ~10-min interruption of 6-MV photon beams did not significantly impact the radiobiological effectiveness, since the RBE decrease was <3%. Nevertheless, the RBE's effect on tumor cells was decreased about 30% by increasing the 60 min interruption time at 8 Gy with four-field irradiation. It is thus necessary to make the interruption time as short as possible.

Intensity Modulated Radiation Fields Induce Protective Effects and Reduce Importance of Dose-Rate Effects

In advanced radiotherapy, intensity modulated radiation fields and complex dose-delivery are utilized to prescribe higher doses to tumours. Here, we investigated the impact of modulated radiation fields on radio-sensitivity and cell recovery during dose delivery. We generated experimental survival data after single-dose, split-dose and fractionated irradiation in normal human skin fibroblast cells (AGO1522) and human prostate cancer cells (DU145). The dose was delivered to either 50% of the area of a T25 flask containing the cells (half-field) or 100% of the flask (uniform-field). We also modelled the impact of dose-rate effects and intercellular signalling on cell-killing. Applying the model to the survival data, it is found that (i) in-field cell survival under half-field exposure is higher than uniform-field exposure for the same delivered dose; (ii) the importance of sub-lethal damage repair (SLDR) in AGO1522 cells is reduced under half-field exposure; (iii) the yield of initial DNA lesions measured with half-field exposure is smaller than that with uniform-field exposure. These results suggest that increased cell survival under half-field exposure is predominantly attributed not to rescue effects (increased SLDR) but protective effects (reduced induction of initial DNA lesions). In support of these protective effects, the reduced DNA damage leads to modulation of cell-cycle dynamics, i.e., less G₁ arrest 6 h after irradiation. These findings provide a new understanding of the impact of dose-rate effects and protective effects measured after modulated field irradiation.

Effect of dose-delivery time for flattened and flattening filter-free photon beams based on microdosimetric kinetic model

The effect of dose-delivery time with flattening filter (FF) and flattening filter-free (FFF) photon beams based on microdosimetric kinetic model (MKM) was investigated in this study. Monte Carlo simulation with the particle and heavy ion transport code system (PHITS) was performed to calculate the dose-mean lineal energy y_D (keV/μm) of FF and FFF 6 MV photon beams using the IAEA phase-space files of Varian TrueBeam linear accelerator. Human non-small cell lung cancer NCI-H460 cells were used to determine the MKM parameters under the condition that dose-delivery times with continuous irradiation were 1, 5, 10, 30, and 60 min, and the adsorbed dose was 2, 4, and 8 Gy in this study. In addition, the relative biological effectiveness (RBE) of FF and FFF photon beams were calculated for evaluating the effect of dose delivery time. The RBE of FF decreased to 99.8% and 97.5% with 5 and 60 min for 2 Gy in comparison to 99.6% and 95.1% for 8 Gy, respectively. Meanwhile, that of FFF decreased to 99.5% and 94.9% with 5 and 60 min for 2 Gy in comparison to 99.5% and 94.9% for 8 Gy, respectively. Dose-delivery time has an effect on the RBE with photon beams. In other words, the dose-delivery time should be considered during radiation therapy. Furthermore, FFF photon beams were an effective irradiation method compared to FF in dose-delivery time on account of improving clinic throughput.

Investigation of dose-rate effects and cell-cycle distribution under protracted exposure to ionizing radiation for various dose-rates

During exposure to ionizing radiation, sub-lethal damage repair (SLDR) competes with DNA damage induction in cultured cells. By virtue of SLDR, cell survival increases with decrease of dose-rate, so-called dose-rate effects (DREs). Here, we focused on a wide dose-rate range and investigated the change of cell-cycle distribution during X-ray protracted exposure and dose-response curves via hybrid analysis with a combination of in vitro experiments and mathematical modelling. In the course of flow-cytometric cell-cycle analysis and clonogenic assays, we found the following responses in CHO-K1 cells: (1) The fraction of cells in S phase gradually increases during 6 h exposure at 3.0 Gy/h, which leads to radio-resistance. (2) Slight cell accumulation in S and G2/M phases is observed after exposure at 6.0 Gy/h for more than 10 hours. This suggests that an increase of SLDR rate for cells in S phase during irradiation may be a reproducible factor to describe changes in the dose-response curve at dose-rates of 3.0 and 6.0 Gy/h. By re-evaluating cell survival for various dose-rates of 0.186-60.0 Gy/h considering experimental-based DNA content and SLDR, it is suggested that the change of S phase fraction during irradiation modulates the dose-response curve and is possibly responsible for some inverse DREs.

Investigation of dose-rate effects and cell-cycle distribution under protracted exposure to ionizing radiation for various dose-rates

During exposure to ionizing radiation, sub-lethal damage repair (SLDR) competes with DNA damage induction in cultured cells. By virtue of SLDR, cell survival increases with decrease of dose-rate, so-called dose-rate effects (DREs). Here, we focused on a wide dose-rate range and investigated the change of cell-cycle distribution during X-ray protracted exposure and dose-response curves via hybrid analysis with a combination of in vitro experiments and mathematical modelling. In the course of flow-cytometric cell-cycle analysis and clonogenic assays, we found the following responses in CHO-K1 cells: (1) The fraction of cells in S phase gradually increases during 6 h exposure at 3.0 Gy/h, which leads to radio-resistance. (2) Slight cell accumulation in S and G₂/M phases is observed after exposure at 6.0 Gy/h for more than 10 hours. This suggests that an increase of SLDR rate for cells in S phase during irradiation may be a reproducible factor to describe changes in the dose-response curve at dose-rates of 3.0 and 6.0 Gy/h. By re-evaluating cell survival for various dose-rates of 0.186-60.0 Gy/h considering experimental-based DNA content and SLDR, it is suggested that the change of S phase fraction during irradiation modulates the dose-response curve and is possibly responsible for some inverse DREs.

Estimation of the radiation-induced DNA double-strand breaks number by considering cell cycle and absorbed dose per cell nucleus

DNA double-strand breaks (DSBs) are thought to be the main cause of cell death after irradiation. In this study, we estimated the probability distribution of the number of DSBs per cell nucleus by considering the DNA amount in a cell nucleus (which depends on the cell cycle) and the statistical variation in the energy imparted to the cell nucleus by X-ray irradiation. The probability estimation of DSB induction was made following these procedures: (i) making use of the Chinese Hamster Ovary (CHO)-K1 cell line as the target example, the amounts of DNA per nucleus in the logarithmic and the plateau phases of the growth curve were measured by flow cytometry with propidium iodide (PI) dyeing; (ii) the probability distribution of the DSB number per cell nucleus for each phase after irradiation with 1.0 Gy of 200 kVp X-rays was measured by means of γ-H2AX immunofluorescent staining; (iii) the distribution of the cell-specific energy deposition via secondary electrons produced by the incident X-rays was calculated by WLTrack (in-house Monte Carlo code); (iv) according to a mathematical model for estimating the DSB number per nucleus, we deduced the induction probability density of DSBs based on the measured DNA amount (depending on the cell cycle) and the calculated dose per nucleus. The model exhibited DSB induction probabilities in good agreement with the experimental results for the two phases, suggesting that the DNA amount (depending on the cell cycle) and the statistical variation in the local energy deposition are essential for estimating the DSB induction probability after X-ray exposure.

Integrated Modelling of Cell Responses after Irradiation for DNA-Targeted Effects and Non-Targeted Effects

Intercellular communication after ionizing radiation exposure, so-called non-targeted effects (NTEs), reduces cell survival. Here we describe an integrated cell-killing model considering NTEs and DNA damage along radiation particle tracks, known as DNA-targeted effects (TEs) based on repair kinetics of DNA damage. The proposed model was applied to a series of experimental data, i.e., signal concentration, DNA damage kinetics, cell survival curve and medium transfer bystander effects (MTBEs). To reproduce the experimental data, the model considers the following assumptions: (i) the linear-quadratic (LQ) function as absorbed dose to express the hit probability to emit cell-killing signals, (ii) the potentially repair of DNA lesions induced by NTEs, and (iii) lower efficiency of repair for the damage in NTEs than that in TEs. By comparing the model results with experimental data, we found that signal-induced DNA damage and lower repair efficiency in non-hit cells are responsible for NTE-related repair kinetics of DNA damage, cell survival curve with low-dose hyper-radiosensitivity (HRS) and MTBEs. From the standpoint of modelling, the integrated cell-killing model with the LQ relation and a different repair function for NTEs provide a reasonable signal-emission probability and a new estimation of low-dose HRS linked to DNA repair efficiency.

Markov chain Monte Carlo analysis for the selection of a cell-killing model under high-dose-rate irradiation

PURPOSE

High-dose-rate irradiation with 6 MV linac x rays is a wide-spread means to treat cancer tissue in radiotherapy. The treatment planning relies on a mathematical description of surviving fraction (SF), such as the linear-quadratic model (LQM) formula. However, even in the case of high-dose-rate treatment, the repair kinetics of DNA damage during dose-delivery time plays a function in predicting the dose-SF relation. This may call the SF model selection into question when considering the dose-delivery time or dose-rate effects (DREs) in radiotherapy and in vitro cell experiments. In this study, we demonstrate the importance of dose-delivery time at high-dose-rate irradiations used in radiotherapy by means of Bayesian estimation.

METHODS

To evaluate the model selection for SF, three types of models, the LQM and two microdosimetric-kinetic models with and without DREs (MKM_DR and MKM) were applied to describe in vitroSF data (our work and references). The parameters in each model were evaluated by a Markov chain Monte Carlo (MCMC) simulation.

RESULTS

The MCMC analysis shows that the cell survival curve by the MKM_DR fits the experimental data the best in terms of the deviance information criterion (DIC). In the fractionated regimen with 30 fractions to a total dose of 60 Gy, the final cell survival estimated by the MKM_DR was higher than that by the LQM. This suggests that additional fractions are required for attaining the total dose equivalent to yield the same effect as the conventional regimen using the LQM in fractionated radiotherapy.

CONCLUSIONS

Damage repair during dose-delivery time plays a key role in precisely estimating cell survival even at a high dose rate in radiotherapy. Consequently, it was suggested that the cell-killing model without repair factor during a short dose-delivery time may overestimate actual cell killing in fractionated radiotherapy.

Educational Activity for the Radiation Emergency System in the Northern Part of Japan: Meeting Report on "The 3rd Educational Symposium on Radiation and Health (ESRAH) by Young Scientists in 2016"

In the northern part of Japan, close cooperation is essential in preparing for any possible emergency response to radiation accidents because several facilities, such as the Low-Level Radioactive Waste Disposal Facility, the MOX Fuel Fabrication Plant and the Vitrified Waste Storage Center, exist in Rokkasho Village (Aomori Prefecture). After the accident at Fukushima Daiichi Nuclear Power Plant in 2011, special attention should be given to the relationship between radiation and human health, as well as establishing a system for managing with a radiation emergency. In the area of Hokkaido and Aomori prefectures in Japan, since 2008 an exchange meeting between Hokkaido University and Hirosaki University has been held every year to have opportunities to discuss radiation effects on human health and to collect the latest news on monitoring environmental radiation. This meeting was elevated to an international meeting in 2014 titled "Educational Symposium on Radiation and Health (ESRAH) by Young Scientists". The 3rd ESRAH meeting took place in 2016, with the theme "Investigating Radiation Impact on the Environmental and Health". Here we report the meeting findings on the continuing educational efforts after the Fukushima incident, what was accomplished in terms of building a community educational approaches, and future goals.