Quantitative estimation of DNA damage by photon irradiation based on the microdosimetric-kinetic model

Towards an adequate description of the dose-response relationship in BNCT of glioblastoma multiforme

BACKGROUND

Boron Neutron Capture Therapy (BNCT) is a binary radiotherapy based on the intravenous administration of a borated drug to the patient and the subsequent irradiation with a low-energy neutron beam. The borated formulation accumulates in the tumor cells, and when neutrons interact with boron, a nuclear capture reaction occurs, releasing high-linear energy transfer, short-range particles that cause lethal damage to the cancer cells. Due to its selectivity, BNCT has the potential to treat aggressive brain tumors such as glioblastoma multiforme (GBM), minimizing the side effects. GBM is a brain neoplasia that poses significant treatment challenges due to its invasiveness and resistance to conventional treatments.

PURPOSE

This work aims to find a suitable model for calculating the photon isoeffective dose for GBM, producing ad hoc radiobiological data to feed the model.

METHODS

To describe adequately the dose-effect relation of BNCT for GBM, the following strategy has been applied 1.We studied the impact of choosing two different photon radiation types (x- or gamma- rays) 2.We assumed that the correct description of the photon-equivalent dose is obtained with the photon isoeffective dose model. This model calculates the photon dose that equals the cell survival obtained with BNCT, taking into account synergism and sub-lethal damage (SLD). 3.Survival curves as a function of the dose for the human GBM U87 cell line were constructed using the clonogenic assays for irradiation with photons (reference), neutron beam, and BNCT. 4.Survival curves were fitted with the modified linear quadratic model, using SLD repair times derived for U87. The radiobiological parameters were determined for the photon isoeffective dose model. 5.The model was applied to a clinical case that received BNCT in Taiwan. Treatment planning has been simulated using an accelerator-based designed neutron beam following the real treatment process and parameters. The results were discussed and compared to the current method, which employs relative biological effectiveness (RBE) factors to obtain BNCT dosimetry in photon-equivalent units.

RESULTS

The dose-survival curves have been obtained with two different photon radiation sources as the reference with a thermal neutron beam and neutrons in the presence of boron. The fitted parameters have been obtained as the input for the photon isoeffective dose and the traditional RBE model. For the first time, the radiobiological parameters of a photon isoeffective dose model were produced for BNCT of GBM. Photon isoeffective dose value can differ up to 32% using gamma photons and low-energy x-rays. Photon isoeffective dose values are lower (17%) than the RBE model currently employed in clinical trials.

CONCLUSION

The results highlight the impact of the reference radiation chosen for the isoeffective dose calculation and the importance of feeding the model with the appropriate radiobiological parameters.The dosimetry obtained with the new radiobiological data is consistent with the dose delivered in modern stereotactic radiotherapy, enabling tumor control predictions.

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.

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.

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.

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.

Chromosomal Location of Genes Differentially Expressed in Tumor Cells Surviving High-Dose X-ray Irradiation: A Preliminary Study on Radio-Fragile Sites

Altered gene expression is a common feature of tumor cells after irradiation. Our previous study showed that this phenomenon is not only an acute response to cytotoxic stress, instead, it was persistently detected in tumor cells that survived 10 Gy irradiation (IR cells). The current understanding is that DNA double-strand breaks (DSBs) are recognized by the phosphorylation of histone H2AX (H2AX) and triggers the ataxia-telangiectasia mutated (ATM) protein or the ATM- and Rad3-related (ATR) pathway, which activate or inactivate the DNA repair or apoptotic or senescence related molecules and causes the expression of genes in many instances. However, because changes in gene expression persist after passaging in IR cells, it may be due to the different pathways from these transient intracellular signaling pathways caused by DSBs. We performed microarray analysis of 30,000 genes in radiation-surviving cells (H1299-IR and MCF7-IR) and found an interesting relation between altered genes and their chromosomal loci. These loci formed a cluster on the chromosome, especially on 1q21 and 6p21-p22 in both irradiated cell lines. These chromosome sites might be regarded as “radio-fragile” sites.

Microdosimetric-Kinetic Model for Radio-enhancement of Gold Nanoparticles: Comparison with LEM

Numerous studies have strongly supported the application of gold nanoparticles (GNPs) as radio-enhanced agents. In our previous study, the local effect model (LEM I) was adopted to predict the cell survival for MDA-MB-231 cells exposed to 150 kVp X rays after 500 µg/ml GNPs treatment. However, microdosimetric quantities could not be obtained, which were correlated with biological effects on cells. Thus, we developed microdosimetric kinetic model (MKM) for GNP radio-enhancement (GNP-MKM), which uses the microdosimetric quantities such as dose-mean lineal energy with subcellular domain size. Using the Monte Carlo simulation tool Geant4, we estimated the dose-mean lineal energy with secondary radiations from GNPs and absorbed dose in the nucleus. The variations in MKM parameters for different domain sizes, and GNP concentrations, were calculated to compare the survival fractions predicted by both models. With a domain radius of 500 nm and a threshold dose of 20 Gy, the sensitizer enhancement ratio predicted by GNP-MKM and GNP-LEM was 1.41 and 1.29, respectively. The GNP-MKM predictions were much more strongly dependent on the domain size than were the GNP-LEM on the threshold dose. These findings provide another method to predict survival fraction for the GNP radio-enhancement.

DNA damage in lens epithelial cells exposed to occupationally-relevant X-ray doses and role in cataract formation

The current framework of radiological protection of occupational exposed medical workers reduced the eye-lens equivalent dose limit from 150 to 20 mSv per year requiring an accurate dosimetric evaluation and an increase understanding of radiation induced effects on Lens cells considering the typical scenario of occupational exposed medical operators. Indeed, it is widely accepted that genomic damage of Lens epithelial cells (LEC) is a key mechanism of cataractogenesis. However, the relationship between apoptosis and cataractogenesis is still controversial. In this study biological and physical data are combined to improve the understanding of radiation induced effects on LEC. To characterize the occupational exposure of medical workers during angiographic procedures an INNOVA 4100 (General Electric Healthcare) equipment was used (scenario A). Additional experiments were conducted using a research tube (scenario B). For both scenarios, the frequencies of binucleated cells, micronuclei, p21-positive cells were assessed with different doses and dose rates. A Monte-Carlo study was conducted using a model for the photon generation with the X-ray tubes and with the Petri dishes considering the two different scenarios (A and B) to reproduce the experimental conditions and validate the irradiation setups to the cells. The simulation results have been tallied using the Monte Carlo code MCNP6. The spectral characteristics of the different X-ray beams have been estimated. All irradiated samples showed frequencies of micronuclei and p21-positive cells higher than the unirradiated controls. Differences in frequencies increased with the delivered dose measured with Gafchromic films XR-RV3. The spectrum incident on eye lens and Petri, as estimated with MCNP6, was in good agreement in the scenario A (confirming the experimental setup), while the mean energy spectrum was higher in the scenario B. Nevertheless, the response of LEC seemed mainly related to the measured absorbed dose. No effects on viability were detected. Our results support the hypothesis that apoptosis is not responsible for cataract induced by low doses of X-ray (i.e. 25 mGy) while the induction of transient p21 may interfere with the disassembly of the nuclear envelop in differentiating LEC, leading to cataract formation. Further studies are needed to better clarify the relationship we suggested between DNA damage, transient p21 induction and the inability of LEC enucleation.

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.

Dose compensation based on biological effectiveness due to interruption time for photon radiation therapy

OBJECTIVE

To evaluate the biological effectiveness of dose associated with interruption time; and propose the dose compensation method based on biological effectiveness when an interruption occurs during photon radiation therapy.

METHODS

The lineal energy distribution for human salivary gland tumor was calculated by Monte Carlo simulation using a photon beam. The biological dose (D_bio) was estimated using the microdosimetric kinetic model. The dose compensating factor with the physical dose for the difference of the D_bio with and without interruption (Δ) was derived. The interruption time (τ) was varied to 0.1, 0.2, 0.3, 0.4, 0.5, 1, 2, 3, 4, 5, 10, 20, 30, 40, 50, 75, and 120 min. The dose per fraction and dose rate varied from 2 to 8 Gy and 0.1 to 24 Gy/min, respectively.

RESULTS

The maximum Δ with 1 Gy/min occurred when the interruption occurred at half the dose. The Δ with 1 Gy/min at half of the dose was over 3% for τ >= 20 min for 2 Gy, τ = 10 min for 5 Gy, and τ = 10 min for 8 Gy. The maximum difference of the Δ due to the dose rate was within 3% for 2 and 5 Gy, and achieving values of 4.0% for 8 Gy. The dose compensating factor was larger with a high dose per fraction and high-dose rate beams.

CONCLUSION

A loss of biological effectiveness occurs due to interruption. Our proposal method could correct for the unexpected decrease of the biological effectiveness caused by interruption time.

ADVANCES IN KNOWLEDGE

For photon radiotherapy, the interruption causes the sublethal damage repair. The current study proposed the dose compensation method for the decrease of the biological effect by the interruption.

A Model for Estimating Dose-Rate Effects on Cell-Killing of Human Melanoma after Boron Neutron Capture Therapy

Boron neutron capture therapy (BNCT) is a type of radiation therapy for eradicating tumor cells through a ¹⁰B(n,α)⁷Li reaction in the presence of ¹⁰B in cancer cells. When delivering a high absorbed dose to cancer cells using BNCT, both the timeline of ¹⁰B concentrations and the relative long dose-delivery time compared to photon therapy must be considered. Changes in radiosensitivity during such a long dose-delivery time can reduce the probability of tumor control; however, such changes have not yet been evaluated. Here, we propose an improved integrated microdosimetric-kinetic model that accounts for changes in microdosimetric quantities and dose rates depending on the ¹⁰B concentration and investigate the cell recovery (dose-rate effects) of melanoma during BNCT irradiation. The integrated microdosimetric–kinetic model used in this study considers both sub-lethal damage repair and changes in microdosimetric quantities during irradiation. The model, coupled with the Monte Carlo track structure simulation code of the Particle and Heavy Ion Transport code System, shows good agreement with in vitro experimental data for acute exposure to ⁶⁰Co γ-rays, thermal neutrons, and BNCT with ¹⁰B concentrations of 10 ppm. This indicates that microdosimetric quantities are important parameters for predicting dose-response curves for cell survival under BNCT irradiations. Furthermore, the model estimation at the endpoint of the mean activation dose exhibits a reduced impact of cell recovery during BNCT irradiations with high linear energy transfer (LET) compared to ⁶⁰Co γ-rays irradiation with low LET. Throughout this study, we discuss the advantages of BNCT for enhancing the killing of cancer cells with a reduced dose-rate dependency. If the neutron spectrum and the timelines for drug and dose delivery are provided, the present model will make it possible to predict radiosensitivity for more realistic dose-delivery schemes in BNCT irradiations.

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.

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.

Track Structure Study for Energy Dependency of Electrons and X-rays on DNA Double-Strand Break Induction

Radiation weighting factor wR for photons and electrons has been defined as unity independently of the energy of the particles. However, the biological effects depend on the incident energies according to in vitro experimental data. In this study, we have quantified the energy concentration along electron tracks in terms of dose-mean lineal energy (yD) on chromosome (micro-meter) and DNA (nano-meter) order scales by Monte Carlo simulations, and evaluated the impact of photon energies on DNA double-strand break (DNA-DSB) induction from an experimental study of irradiated cells. Our simulation result shows that the yD values for diagnostic X-rays (60-250 kVp) are higher than that for therapeutic X-rays (linac 6 MV), which agrees well with the tissue equivalent proportional counter (TEPC) measurements. The relation between the yD values and the numbers of γ-H2AX foci for various photon energy spectra suggests that low energy X-rays induce DNA-DSB more efficiently than higher energy X-rays even at the same absorbed dose (e.g., 1.0 Gy). The relative biological effectiveness based on DNA-DSBs number (RBEDSB) is proportionally enhanced as the yD value increases, demonstrating that the biological impact of the photon irradiation depends on energy concentration along radiation tracks of electrons produced in the bio-tissues. Ultimately, our study implies that the value of wR for photons varies depending on their energies.

DNA damage induction during localized chronic exposure to an insoluble radioactive microparticle

Insoluble radioactive microparticles emitted by the incident at the Fukushima nuclear power plant have drawn keen interests from the viewpoint of radiation protection. Cs-bearing particles have been assumed to adhere in the long term to trachea after aspirated into respiratory system, leading to heterogeneous dose distribution within healthy tissue around the particles. However, the biological effects posed by an insoluble radioactive particle remain unclear. Here, we show cumulative DNA damage in normal human lung cells proximal and distal to the particle (β-ray and γ-ray-dominant areas, respectively) under localized chronic exposure in comparison with uniform exposure. We put a Cs-bearing particle into a microcapillary tip and placed it onto a glass-base dish containing fibroblast or epithelial cells cultured in vitro. A Monte Carlo simulation with PHITS code provides the radial distribution of absorbed dose-rate around the particle, and subsequently we observed a significant change in nuclear γ-H2AX foci after 24 h or 48 h exposure to the particle. The nuclear foci in the cells distal to the particle increased even under low-dose-rate exposure compared with uniform exposure to 137Cs γ-rays, which was suppressed by a treatment with a scavenger of reactive oxygen species. In contrast, such focus formation was less manifested in the exposed cells proximal to the particle compared with uniform exposure. These data suggest that the localized exposure to a Cs-bearing particle leads to not only disadvantage to distal cells but also advantage to proximal cells. This study is the first to provide quantitative evaluation for the spatial distribution of DNA double strand breaks after the heterogeneous chronic exposure to a Cs-bearing particle in comparison with uniform Cs exposure.

Rationale for Combining Radiotherapy and Immune Checkpoint Inhibition for Patients With Hypoxic Tumors

In order to compensate for the increased oxygen consumption in growing tumors, tumors need angiogenesis and vasculogenesis to increase the supply. Insufficiency in this process or in the microcirculation leads to hypoxic tumor areas with a significantly reduced pO2, which in turn leads to alterations in the biology of cancer cells as well as in the tumor microenvironment. Cancer cells develop more aggressive phenotypes, stem cell features and are more prone to metastasis formation and migration. In addition, intratumoral hypoxia confers therapy resistance, specifically radioresistance. Reactive oxygen species are crucial in fixing DNA breaks after ionizing radiation. Thus, hypoxic tumor cells show a two- to threefold increase in radioresistance. The microenvironment is enriched with chemokines (e.g., SDF-1) and growth factors (e.g., TGFβ) additionally reducing radiosensitivity. During recent years hypoxia has also been identified as a major factor for immune suppression in the tumor microenvironment. Hypoxic tumors show increased numbers of myeloid derived suppressor cells (MDSCs) as well as regulatory T cells (Tregs) and decreased infiltration and activation of cytotoxic T cells. The combination of radiotherapy with immune checkpoint inhibition is on the rise in the treatment of metastatic cancer patients, but is also tested in multiple curative treatment settings. There is a strong rationale for synergistic effects, such as increased T cell infiltration in irradiated tumors and mitigation of radiation-induced immunosuppressive mechanisms such as PD-L1 upregulation by immune checkpoint inhibition. Given the worse prognosis of patients with hypoxic tumors due to local therapy resistance but also increased rate of distant metastases and the strong immune suppression induced by hypoxia, we hypothesize that the subgroup of patients with hypoxic tumors might be of special interest for combining immune checkpoint inhibition with radiotherapy.

Quantitative Estimation of the Equivalent Radiation Dose Escalation using Radiofrequency Hyperthermia in Mouse Xenograft Models of Human Lung Cancer

Hyperthermia is a potent radiosensitizer, and its effect varies according to the different types of cancer cells. In the present study, the radiosensitizing effect of hyperthermia on lung cancer cell lines A549 and NCI-H1299 was determined based on the equivalent radiation dose escalation. In vitro cell experiments were conducted using lung cancer cell lines A549 and NCI-H1299 to determine thermal radiosensitivity. In vivo experiments were conducted using mouse heterotopic xenograft models to determine the treatment response and increase in the temperature of tumors using a 13.56 MHz radiofrequency (RF) hyperthermia device. Using the α and β values of the linear–quadratic equations of cell survival curves, numerical simulations were performed to calculate the equivalent radiation dose escalations. The dielectric properties of tumors were measured, and their effect on the calculated equivalent radiation dose was analyzed. Hyperthermia increased the equivalent radiation dose of lung cancer xenografts and a higher escalation was found in NCI-H1299 cells compared with that observed in A549 cells. An underestimation of the calculated equivalent radiation dose was observed when the dielectric property of the tumor was varied. This study may contribute to the effective planning of thermoradiotherapy in clinics.

A comparison of mechanism-inspired models for particle relative biological effectiveness (RBE)

BACKGROUND AND SIGNIFICANCE

The application of heavy ion beams in cancer therapy must account for the increasing relative biological effectiveness (RBE) with increasing penetration depth when determining dose prescriptions and organ at risk (OAR) constraints in treatment planning. Because RBE depends in a complex manner on factors such as the ion type, energy, cell and tissue radiosensitivity, physical dose, biological endpoint, and position within and outside treatment fields, biophysical models reflecting these dependencies are required for the personalization and optimization of treatment plans.

AIM

To review and compare three mechanism-inspired models which predict the complexities of particle RBE for various ion types, energies, linear energy transfer (LET) values and tissue radiation sensitivities.

METHODS

The review of models and mechanisms focuses on the Local Effect Model (LEM), the Microdosimetric-Kinetic (MK) model, and the Repair-Misrepair-Fixation (RMF) model in combination with the Monte Carlo Damage Simulation (MCDS). These models relate the induction of potentially lethal double strand breaks (DSBs) to the subsequent interactions and biological processing of DSB into more lethal forms of damage. A key element to explain the increased biological effectiveness of high LET ions compared to MV x rays is the characterization of the number and local complexity (clustering) of the initial DSB produced within a cell. For high LET ions, the spatial density of DSB induction along an ion's trajectory is much greater than along the path of a low LET electron, such as the secondary electrons produced by the megavoltage (MV) x rays used in conventional radiation therapy. The main aspects of the three models are introduced and the conceptual similarities and differences are critiqued and highlighted. Model predictions are compared in terms of the RBE for DSB induction and for reproductive cell survival.

RESULTS AND CONCLUSIONS

Comparisons of the RBE for DSB induction and for cell survival are presented for proton (¹ H), helium (⁴ He), and carbon (¹² C) ions for the therapeutically most relevant range of ion beam energies. The reviewed models embody mechanisms of action acting over the spatial scales underlying the biological processing of potentially lethal DSB into more lethal forms of damage. Differences among the number and types of input parameters, relevant biological targets, and the computational approaches among the LEM, MK and RMF models are summarized and critiqued. Potential experiments to test some of the seemingly contradictory aspects of the models are discussed.

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.

Relative biological effectiveness study of Lipiodol based on microdosimetric-kinetic model

OBJECTIVES

We examine the contrast agent Lipiodol effect on the relative biological effectiveness (RBE) values for flattening filter free (FFF) and flattening filter (FF) beams of 6 MV-Xray (6 MVX) and 10 MVX.

METHODS

Lipiodol was placed at 5 cm depth in water. According to the microdosimetric kinetic model, the RBE values for killing the human liver hepatocellular cells were calculated from dose and lineal energy (yd(y)) from Monte Carlo simulations. RBE_200kVX and RBE_Co were defined as the ratios of dose using reference radiation (200 kVX, Co-ɤ) to the dose of test radiation (FFF and FF beams for 6 MV and 10 MV) to produce the same biological effects. The dose enhancement RBE (RBE_DE) was defined as the ratios of a dose without Lipiodol to with Lipiodol using to produce the same biological effects. The dose needed to achieve 10% (D_10%) and 1% cell survival (D_1%) was evaluated by cell surviving fraction (SF) formula.

RESULTS

The deviation of mean y‾_D values with and without Lipiodol were 3.9-4.8% for 6 MVX and 3.5-3.6% for 10 MVX. The RBE_200kVX and RBE_Co with Lipiodol were larger than that without Lipiodol. The RBE_DE was larger for FFF beam than for FF beam. The deviation of RBE_DE for FFF and FF beams of 6 MVX was larger than that of 10 MVX.

CONCLUSION

The presence of Lipiodol seemed to locally increase the absorbed dose and to also cause an enhancement of the relative biological effectiveness.

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.

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.

Comparison of dicentric dose response curves of 6MV LINAC X-rays and 60Co γ-rays for biodosimetry application

In vitro Dicentric (DC) dose response curves of 6MV X-rays (3Gy/min) and ⁶⁰Co γ-rays (0.63Gy/min) were generated and compared (intra-laboratory) to understand their similarities and differences. Human peripheral blood samples exposed to ten different doses (0, 0.05, 0.1, 0.25, 0.5, 1, 2, 3, 4 and 5Gy) of 6MV X-rays and ⁶⁰Co γ-rays were subjected to Dicentric Chromosome Assay (DCA) and dicentrics, excess acentric fragments (AF) and chromatid breaks (CB) were scored. Totally about 31,553 metaphase spreads were scored for the purpose. Dose response curves of both radiation qualities were almost same except for a 13.8% higher β value for 6MV X-rays. However, blind tests results revealed that both these curves are biologically equivalent and exhibited good dose prediction accuracy for the entire dose range. This demonstrated the feasibility of interchangeable use of these curves in biodosimetry. Consequently it has been suggested that LINAC facilities worldwide can be roped in for biodosimetry capacity augmentation towards managing nuclear emergency situations involving γ-radiation exposures.

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.

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.

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

We have investigated the dose rate effects on cell damage caused by photon-beam irradiation. During a relatively long dose-delivery time with a low dose rate, lesions created in cells may undergo some reactions, such as DNA repair. In order to investigate these reactions quantitatively, we adopted the microdosimetric-kinetic (MK) model and deduced a cell surviving fraction (SF) formula for continuous irradiation. This model enabled us to estimate the SF from dose and dose rate. The parameters in the MK model were determined so as to generate the SF, and we attempted to evaluate the dose rate effects on the SF. To deduce the cell-specific parameters in the SF formula, including the dose rate, we performed a split-dose experiment and a single-dose experiment with a constant dose-delivery time (10 min) (to retain the condition for equivalent behavior of cell lesions) by means of a clonogenic assay. Then, using the MK model parameters, the SFs were reproduced for a variety of dose rates (1.0, 0.31, 0.18, 0.025 and 0.0031 Gy/min) and were compared with reported experimental data. The SF curves predicted by the MK model agreed well with the experimental data, suggesting that the dose rate effects appear in the kinetics of cell lesions during the dose-delivery time. From fitting the analysis of the model formula to the experimental data, it was shown that the MK model could illustrate the characteristics of log-SF in a rectilinear form at a high dose range with a relatively low dose rate.