Evaluation of the cell survival curve under radiation exposure based on the kinetics of lesions in relation to dose-delivery time

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 the biological effects of boron neutron capture therapy through cellular DNA damage repair model

BACKGROUND

Boron neutron capture therapy (BNCT) is a targeted radiotherapy that relies on the 10B (n, α) 7Li reaction, which produces secondary particles with high linear energy transfer (LET), leading to a high relative biological effectiveness (RBE) in tumors. The biological effectiveness of BNCT is influenced by factors such as boron distribution and concentration, necessitating improved methods for assessing its radiobiological effects and clarifying the sensitivity of the differences in different factors to the biological effects.

PURPOSE

This paper introduces a method to evaluate the biological effects of BNCT using the cellular repair model. This method aims to overcome some of the limitations of current evaluation approaches. The primary goal is to provide guidance for clinical treatments and the development of boron drugs, as well as to investigate the impact of the synergistic effects of mixed radiation fields in BNCT on treatment outcomes.

METHODS

The approach involves three key steps: first, extending the radial energy deposition distribution of BNCT secondary particles using Geant4-DNA. This allows for the calculation of initial DNA double-strand breaks (DSBs) distributions for a given absorbed dose. Next, the obtained initial DSB distributions are used for DNA damage repair simulations to generate cell survival curves, then thereby quantifying RBE and compound biological effectiveness (CBE). The study also explores the synergistic effects of the mixed radiation fields in BNCT on assessing biological effects were also explored in depth.

RESULTS

The results showed that the RBE of boronophenylalanine (BPA) and sodium borocaptate (BSH) drugs at cell survival fraction 0.01 was 2.50 and 2.15, respectively. The CBE of the boron dose component was 3.60 and 0.73, respectively, and the RBE of the proton component was 3.21, demonstrating that BPA has a significantly higher biological impact than BSH due to superior cellular permeability. The proton dose significance in BNCT treatment is also underscored, necessitating consideration in both experimental and clinical contexts. The study demonstrates that synergistic effects between disparate radiation fields lead to increased misrepairs and enhanced biological impact. Additionally, the biological effect diminishes with rising boron concentration, emphasizing the need to account for intercellular DNA damage heterogeneity.

CONCLUSIONS

This methodology offers valuable insights for the development of new boron compounds and precise assessment of bio-weighted doses in clinical settings and can be adapted to other therapeutic modalities.

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.

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.

Low dose rate radiation impairs early follicles in young mice

Low-dose radiation is generally considered less harmful than high-dose radiation. However, its impact on ovaries remains debated. Since previous reports predominantly employed low-dose radiation delivered at a high dose rate on the ovary, the effect of low-dose radiation at a low dose rate on the ovary remains unknown. We investigated the effect of low-dose ionizing radiation delivered at a low dose rate on murine ovaries. Three- and ten-week-old mice were exposed to 0.1 and 0.5 Gy of radiation at a rate of 6 mGy/h and monitored after 3 and 30 days. While neither body weight nor ovarian area showed significant changes, ovarian cells were damaged, showing apoptosis and a decrease in cell proliferation after exposure to 0.1 and 0.5 Gy radiation. Follicle numbers decreased over time in both age groups proportionally to the radiation dose. Younger mice were more susceptible to radiation damage, as evidenced by decreased follicles in 3-week-old mice after 30 days of 0.1 Gy exposure, while 10-week-old mice showed reduced follicles only following 0.5 Gy exposure. Primordial or primary follicles were the most vulnerable to radiation. These findings suggest that even low-dose radiation, delivered at a low dose rate, can adversely affect ovarian function, particularly in the early follicles of younger mice.

Implications of "flash" radiotherapy for biodosimetry

Extremely high dose rate radiation delivery (FLASH) for cancer treatment has been shown to produce less damage to normal tissues while having the same radiotoxic effect on tumor tissue (referred to as the FLASH effect). Research on the FLASH effect has two very pertinent implications for the field of biodosimetry: (1) FLASH is a good model to simulate delivery of prompt radiation from the initial moments after detonating a nuclear weapon and (2) the FLASH effect elucidates how dose rate impacts the biological mechanisms that underlie most types of biological biodosimetry. The impact of dose rate will likely differ for different types of biodosimetry, depending on the specific underlying mechanisms. The greatest impact of FLASH effects is likely to occur for assays based on biological responses to radiation damage, but the consequences of differential effects of dose rates on the accuracy of dose estimates has not been taken into account.

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.

Differential Effects of Low and High Radiation Dose Rates on Mouse Spermatogenesis

The adverse effects of radiation are proportional to the total dose and dose rate. We aimed to investigate the effects of radiation dose rate on different organs in mice. The mice were subjected to low dose rate (LDR, ~3.4 mGy/h) and high dose rate (HDR, ~51 Gy/h) radiation. LDR radiation caused severe tissue toxicity, as observed in the histological analysis of testis. It adversely influenced sperm production, including sperm count and motility, and induced greater sperm abnormalities. The expression of markers of early stage spermatogonial stem cells, such as Plzf, c-Kit, and Oct4, decreased significantly after LDR irradiation, compared to that following exposure of HDR radiation, in qPCR analysis. The compositional ratios of all stages of spermatogonia and meiotic cells, except round spermatid, were considerably reduced by LDR in FACS analysis. Therefore, LDR radiation caused more adverse testicular damage than that by HDR radiation, contrary to the response observed in other organs. Therefore, the dose rate of radiation may have differential effects, depending on the organ; it is necessary to evaluate the effect of radiation in terms of radiation dose, dose rate, organ type, and other conditions.

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.

A theoretical cell-killing model to evaluate oxygen enhancement ratios at DNA damage and cell survival endpoints in radiation therapy

Radio-resistance induced under low oxygen pressure plays an important role in malignant progression in fractionated radiotherapy. For the general approach to predict cell killing under hypoxia, cell-killing models (e.g. the Linear-Quadratic model) have to be fitted to in vitro experimental survival data for both normoxia and hypoxia to obtain the oxygen enhancement ratio (OER). In such a case, model parameters for every oxygen condition needs to be considered by model-fitting approaches. This is inefficient for fractionated irradiation planning. Here, we present an efficient model for fractionated radiotherapy the integrated microdosimetric-kinetic model including cell-cycle distribution and the OER at DNA double-strand break endpoint (OER_DSB). The cell survival curves described by this model can reproduce the in vitro experimental survival data for both acute and chronic low oxygen concentrations. The OER_DSB used for calculating cell survival agrees well with experimental DSB ratio of normoxia to hypoxia. The important parameters of the model are oxygen pressure and cell-cycle distribution, which enables us to predict cell survival probabilities under chronic hypoxia and chronic anoxia. This work provides biological effective dose (BED) under various oxygen conditions including its uncertainty, which can contribute to creating fractionated regimens for multi-fractionated radiotherapy. If the oxygen concentration in a tumor can be quantified by medical imaging, the present model will make it possible to estimate the cell-killing and BED under hypoxia in more realistic intravital situations.

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.

Radiobiological effects of flattening filter-free photon beams on A549 non-small-cell lung cancer cells

Flattening filter-free (FFF) photon beams minimize the intrafraction motion of tumors, and this feature is useful in pulmonary malignancies, such as non-small-cell lung cancer (NSCLC). However, the radiobiological effects of such beams on NSCLC cells, which are often treated with stereotactic body radiotherapy (SBRT), have not been investigated sufficiently. Although cell motility may be promoted by photon beams with a low dose, the relationship between cell motility and the dose rate of photon beams has not been evaluated. The purpose of this study was to evaluate the radiobiological effects of FFF photon beams on cell survival and motility in NSCLC. A human lung cancer cell line (A549) was irradiated with conventional flattening filter (FF) and FFF photon beams at dose rates of 300 (FF), 500 and 2000 MU/min (FFF). While cell survival was estimated using the colony formation assay, cell motility was evaluated using the Boyden chamber and Matrigel invasion assays. FFF photon beams with a high dose rate neither affected the survival of A549 cells nor caused any significant difference in their motility. On the other hand, high-dose irradiation reduced cell survival and motility regardless of the dose rate. Photon beams with a high dose rate used for radiation therapy are suitable for SBRT from the standpoint of both cell survival and motility, in addition to their physical characteristics.

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.

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.

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

Hyper-radiosensitivity (HRS) is a well-known bioresponse under low-dose or low-dose-rate exposures. Although disorder of the DNA repair function, non-targeted effects and accumulation of cells in G2 have been experimentally observed, the mechanism for inducing HRS by long-term irradiation is still unclear. On the basis of biological experiments and a theoretical study, we have shown that change in the amount of DNA associated with accumulation of cells in G2 enhances radiosensitivity. To demonstrate continuous irradiation with 250 kVp X-rays, we adopted a fractionated regimen of 0.186 or 1.00 Gy per fraction at intervals of 1 h (i.e. 0.186 Gy/h, 1.00 Gy/h on average) to Chinese Hamster Ovary (CHO)-K1 cells. The change in the amount of DNA during irradiation was quantified by flow cytometric analysis with propidium iodide (PI). Concurrently, we attempted a theoretical evaluation of the DNA damage by using a microdosimetric-kinetic (MK) model that was modified to incorporate the change in the amount of DNA. Our experimental results showed that the fraction of the cells in G2/M phase increased by 6.7% with 0.186 Gy/h and by 22.1% with 1.00 Gy/h after the 12th irradiation. The MK model considering the change in amount of DNA during the irradiation exhibited a higher radiosensitivity at a high dose range, which could account for the experimental clonogenic survival. The theoretical results suggest that HRS in the high dose range is associated with an increase in the total amount of DNA during irradiation.