Nature of the physicochemical process in water photolysis uncovered by a computer simulation.

In this work, we investigate the physicochemical process of water photolysis to bridge physical and chemical processes by a newly developed first-principles calculation code. The deceleration, thermalization, delocalization, and initial hydration of the extremely low-energy electrons ejected by water photolysis are sequentially tracked in the condensed phase. We show herein the calculated results for these sequential phenomena during 300 fs. Our results indicate that the mechanisms heavily depend on the intermolecular vibration and rotation modes peculiar to water and the momentum transfer between the electrons and the water medium. We suggest that using our results for the delocalized electron distribution will reproduce successive chemical reactions measured by photolysis experiments using a chemical reaction code. We expect our approach to become a powerful technique for various scientific fields related to water photolysis and radiolysis.

No Intercellular Regulation of the Cell Cycle among Human Cervical Carcinoma HeLa Cells Expressing Fluorescent Ubiquitination-Based Cell-Cycle Indicators in Modulated Radiation Fields.

The non-targeted effects of radiation have been known to induce significant alternations in cell survival. Although the effects might govern the progression of tumor sites following advanced radiotherapy, the impacts on the intercellular control of the cell cycle following radiation exposure with a modified field, remain to be determined. Recently, a fluorescent ubiquitination-based cell-cycle indicator (FUCCI), which can visualize the cell-cycle phases with fluorescence microscopy in real time, was developed for biological cell research. In this study, we investigated the non-targeted effects on the regulation of the cell cycle of human cervical carcinoma (HeLa) cells with imperfect p53 function that express the FUCCI (HeLa-FUCCI cells). The possible effects on the cell-cycle phases via soluble factors were analyzed following exposure to different field configurations, which were delivered using a 150 kVp X-ray irradiator. In addition, using synchrotron-generated, 5.35 keV monochromatic X-ray microbeams, high-precision 200 µm-slit microbeam irradiation was performed to investigate the possible impacts on the cell-cycle phases via cell-cell contacts. Collectively, we could not detect the intercellular regulation of the cell cycle in HeLa-FUCCI cells, which suggested that the unregulated cell growth was a malignant tumor. Our findings indicated that there was no significant intercellular control system of the cell cycle in malignant tumors during or after radiotherapy, highlighting the differences between normal tissue and tumor characteristics.

Spirulina platensis and its ingredient biopterin glucoside improved insulin sensitivity in non-alcoholic steatohepatitis model.

Non-alcoholic steatohepatitis is the chronic liver disease leading to cirrhosis and cancer and its prevalence is increasing. Some agents are under clinical trials for non-alcoholic steatohepatitis treatment. We previously reported Spirulina (Arthrospira) platensis effectively prevented non-alcoholic steatohepatitis progression in our model rats. The contribution of phycocyanin, an ingredient of Spirulina (Arthrospira) platensis, was limited. We, therefore, have looked for more active components of Spirulina (Arthrospira) platensis. In this study, we pursued the effect of biopterin glucoside, another bioactive ingredient of Spirulina (Arthrospira) platensis. We found Spirulina (Arthrospira) platensis and biopterin glucoside oral administrations effectively alleviated oxidative stress, inflammation and insulin signal failure, and prevented fibroblast growth factor 21 gene overexpression in non-alcoholic steatohepatitis rat livers. We concluded biopterin glucoside is a major component of Spirulina (Arthrospira) platensis action.

Spatially Fractionated Microbeam Analysis of Tissue-sparing Effect for Spermatogenesis.

Spatially fractionated radiation therapy (SFRT) has been based on the delivery of a single high-dose fraction to a large treatment area that has been divided into several smaller fields, reducing the overall toxicity and adverse effects. Complementary microbeam studies have also shown an effective tissue-sparing effect (TSE) in various tissue types and species after spatially fractionated irradiation at the microscale level; however, the underlying biological mechanism remains elusive. In the current study, using the combination of an ex vivo mouse spermatogenesis model and high-precision X-ray microbeams, we revealed the significant TSE for maintaining spermatogenesis after spatially fractionated microbeam irradiation. We used the following ratios of the irradiated to nonirradiated areas: 50:50, 150:50 and 350:50 µm-slit, where approximately 50, 75 and 87.5% of the sample was irradiated (using center-to-center distances of 100, 200 and 400 µm, respectively). We found that the 50 and 75% micro-slit irradiated testicular tissues showed an almost unadulterated TSE for spermatogenesis, whereas the 87.5% micro-slit irradiated tissues showed an incomplete TSE. This suggests that the TSE efficiency for spermatogenesis is dependent on the size of the nonirradiated spermatogonial stem cell pool in the irradiated testicular tissues. In addition, there would be a spatiotemporal limitation of stem cell migration/competition, resulting in the insufficient TSE for 87.5% micro-slit irradiated tissues. These stem cell characteristics are essential for the accurate prediction of tissue-level responses during or after SFRT, indicating the clinical potential for achieving better outcomes while preventing adverse effects.

The Tissue-Sparing Effect of Spatially Fractionated X-rays for Maintaining Spermatogenesis: A Radiobiological Approach for the Preservation of Male Fertility after Radiotherapy.

Radiotherapy can result in temporary or permanent gonadal toxicity in male cancer patients despite the high precision and accuracy of modern radiation treatment techniques. Previous radiobiological studies have shown an effective tissue-sparing response in various tissue types and species following exposure to spatially fractionated radiation. In the present study, we used an ex vivo mouse testicular tissue culture model and a conventional X-ray irradiation device to evaluate the tissue-sparing effect (TSE) of spatially fractionated X-rays for the protection of male fertility from radiotherapy-related adverse effects. We revealed a significant TSE for maintaining spermatogenesis in the ex vivo testes model following spatially fractionated X-ray irradiation. Moreover, we experimentally propose a possible mechanism by which the migration of spermatogonial cells, from the non-irradiated areas to the irradiated ones, in irradiated testicular tissue, is essential for the TSE and maintaining spermatogenesis. Therefore, our findings demonstrate that the control of TSE following spatially fractionated X-rays in the testes has a considerable potential for clinical application. Interdisciplinary research will be essential for further expanding the applicability of this method as an approach for the preservation of male fertility during or after radiotherapy.

High-precision microbeam radiotherapy reveals testicular tissue-sparing effects for male fertility preservation.

Microbeam radiotherapy (MRT) is based on a spatial fractionation of synchrotron X-ray microbeams at the microscale level. Although the tissue-sparing effect (TSE) in response to non-uniform radiation fields was recognized more than one century ago, the TSE of MRT in the testes and its clinical importance for preventing male fertility remain to be determined. In this study, using the combination of MRT techniques and a unique ex vivo testes organ culture, we show, for the first time, the MRT-mediated TSE for the preservation of spermatogenesis. Furthermore, our high-precision microbeam analysis revealed that the survival and potential migration steps of the non-irradiated germ stem cells in the irradiated testes tissue would be needed for the effective TSE for spermatogenesis. Our findings indicated the distribution of dose irradiated in the testes at the microscale level is of clinical importance for delivering high doses of radiation to the tumor, while still preserving male fertility.

Application of an Ex Vivo Tissue Model to Investigate Radiobiological Effects on Spermatogenesis.

The formation of sperm by the testes through the process of spermatogenesis is highly radiosensitive and can be affected by environmental, occupational and therapeutic radiation exposures. In this study, we applied an ex vivo mouse testis organ culture as an experimental model of spermatogenesis to investigate the radiobiological effects and to demonstrate its feasibility as a tool to determine response to complex, modulated radiation fields. This model uses Acr-GFP transgenic mice, which express the marker green fluorescent proteins specific for meiosis to allow observation of functional changes in real-time that can be used to analyze radiation-induced changes in the process of spermatogenesis. Our results showed that the model can accurately reproduce radiation-induced male germ cell toxicity, such as temporary infertility and permanent sterility. Furthermore, using a monochromatic X-ray microbeam, we applied this model to investigate the effects of heterogeneous radiation fields on testis tissue ex vivo. Our model represents a unique application in the field, which offers significant potential for gaining further mechanistic insight into radiation effects on the process of spermatogenesis.

A significant role of non-thermal equilibrated electrons in the formation of deleterious complex DNA damage.

Although most of the radiation damage to genomic DNA could be rendered harmless using repair enzymes in a living cell, a certain fraction of the damage is persistent resulting in serious genetic effects, such as mutation induction. In order to understand the mechanisms of the deleterious DNA damage formation in terms of its earliest physical stage at the radiation track end, dynamics of low energy electrons and their thermalization processes around DNA molecules were investigated using a dynamic Monte Carlo code. The primary incident (1 keV) electrons multiply collide within 1 nm (equivalent to three DNA-base-pairs, 3bp) and generate secondary electrons which show non-Gaussian and non-thermal equilibrium distributions within 300 fs. On the other hand, the secondary electrons are mainly distributed within approximately 10 nm from their parent cations although approximately 5% of the electrons are localized within 1 nm of the cations owing to the interaction of their Coulombic fields. The mean electron energy is 0.7 eV; however, more than 10% of the electrons fall into a much lower-energy region than 0.1 eV at 300 fs. These results indicate that pre-hydrated electrons are formed from the extremely decelerated electrons over a few nm from the cations. DNA damage sites comprising multiple nucleobase lesions or single strand breaks can therefore be formed by multiple collisions of these electrons within 3bp. This multiple damage site is hardly processed by base excision repair enzymes. However, pre-hydrated electrons can also be produced resulting in an additional base lesion (or a strand break) more than 3bp away from the multi-damage site. These damage sites may be finally converted into a double strand break (DSB) when base excision enzymes process the additional base lesions. This DSB includes another base lesion(s) at their termini, and may introduce miss-rejoining by DSB repair enzymes, and hence may result in biological effects such as mutation in surviving cells.

Dynamic Behavior of Secondary Electrons in Liquid Water at the Earliest Stage upon Irradiation: Implications for DNA Damage Localization Mechanism.

To clarify the formation of radiation damage in DNA, the dynamic behavior of low-energy secondary electrons produced by ionizing radiation in water was studied by using a dynamic Monte Carlo code that considers the Coulombic force between electrons and their parent cations. The calculated time evolution of the mean energy, total track length, and mean traveling distance of the electrons indicated that the prehydration of the electrons occurs competitively with thermalization on a time scale of hundreds of femtoseconds. The decelerating electrons are gradually attracted to their parent cations by Coulombic force within hundreds of femtoseconds, and finally about 12.6% electrons are distributed within 2 nm of the cations. The collision fraction for ionization and electronic excitation within 1 nm of the cation was estimated to be about 40%. If these electrons are decelerated in a living cell, they may cause highly localized lesions around a cation in a DNA molecule through additional dissociative electron transfer (DET) as well as ionization and electronic excitation (EXC), possibly resulting in cell death or mutation.

Deceleration processes of secondary electrons produced by a high-energy Auger electron in a biological context.

PURPOSE: To simulate the deceleration processes of secondary electrons produced by a high-energy Auger electron in water, and particularly to focus on the spatial and temporal distributions of the secondary electron and the collision events (e.g. ionization, electronic excitation, and dissociative electron attachment) that are involved in the multiplication of lesions at sites of DNA damage. MATERIALS AND METHODS: We developed a dynamic Monte Carlo code that considers the Coulombic force between an ejected electron and its parent cation produced by the Auger electron in water. Thus our code can simulate some return electrons to the parent cations. Using the code, we calculated to within the order of femtoseconds the temporal evolution of collision events, the mean energy, and the mean traveling distance (including its spatial probability distribution) of the electron at an ejected energy of 20 eV. RESULTS: Some of the decelerating electrons in water in the Coulombic field were attracted to the ionized atoms (cations) by the Coulombic force within hundreds of femtoseconds, although the force did not significantly enhance the number of ionization, electronic excitation, and dissociative electron attachment collision events leading to water radiolysis. CONCLUSIONS: The secondary electrons are decelerated in water by the Coulombic force and recombined to the ionized atoms (cations). Furthermore, the some return electrons might be prehydrated in water layer near the parent cation in DNA if the electrons might be emitted from the DNA. The prehydrated electron originated from the return electron might play a significant role in inducing DNA damage.

Enhancement of DNA double-strand break induction and cell killing by K-shell absorption of phosphorus in human cell lines.

PURPOSE: To investigate an enhancement of DNA double-strand break (DSB) induction and cell killing effect by K-shell ionization of phosphorus atoms and Auger electrons on human cell lines. MATERIALS AND METHODS: Induction of DSB, DNA damage responses, cell cycle distributions, and cell killing effects were investigated after exposures of the cells with monochromatic synchrotron radiation soft X-rays of 2153 and 2147 eV, which were the resonance peak and off peak, respectively, of the K-shell photoabsorption of phosphorus. RESULTS: Higher biological effects in the cells irradiated with soft X-rays at 2153 eV than at 2147 eV were observed in (i) the efficiency of 53BP1/γ-H2AX co-localized foci formation per dose and residual number of foci, (ii) prolonged phosphorylation levels of DSB repair and/or cell cycle checkpoint related proteins and G2 arrest, (iii) the cell killing effects at the 10% survival level of normal human fibroblasts, HeLa cells, and human glioblastoma M059K cells (1.2-1.5 times higher) and that of human ataxia telangiectasia mutated (ATM)-defective cells and glioblastoma DNA-dependent protein kinase catalytic subunit (DNA-PKcs)-defective cells (1.2 times). CONCLUSION: The yield of DSB and partly less-reparable complex DNA damage induction in human cells was enhanced by K-shell photoabsorption of phosphorus and low-energy Auger electrons.

Evaluation of DNA damage induced by Auger electrons from 137Cs.

PURPOSE: To understand the biological effect of external and internal exposure from 137Cs, DNA damage spectrum induced by directly emitted electrons (γ-rays, internal conversion electrons, Auger electrons) from 137Cs was compared with that induced by 137Cs γ-rays. METHODS: Monte Carlo track simulation method was used to calculate the microscopic energy deposition pattern in liquid water. Simulation was performed for the two simple target systems in microscale. Radiation sources were placed inside for one system and outside for another system. To simulate the energy deposition by directly emitted electrons from 137Cs placed inside the system, the multiple ejections of electrons after internal conversion were considered. In the target systems, induction process of DNA damage was modeled and simulated for both direct energy deposition and the water radical reaction on the DNA. The yield and spatial distribution of simple and complex DNA damage including strand breaks and base lesions were calculated for irradiation by electrons and γ-rays from 137Cs. RESULTS: The simulation showed that the significant difference in DNA damage spectrum was not caused by directly ejected electrons and γ-rays from 137Cs. CONCLUSIONS: The result supports the existing perception that the biological effects by internal and external exposure by 137Cs are equivalent.

Cellular automaton-based model for radiation-induced bystander effects.

BACKGROUND: The radiation-induced bystander effect is a biological response observed in non-irradiated cells surrounding an irradiated cell. The bystander effect is known to be induced by two intercellular signaling pathways, the medium-mediated pathway (MDP) and the gap junctional pathway (GJP). To investigate the relative contribution of each signaling pathway, we have developed a mathematical model of the cellular response through these two pathways, with a particular focus on cell-cycle modification. METHODS: The model is based on a cellular automaton and consists of four components: (1) irradiation, (2) generation and diffusion of intercellular signals, (3) induction of DNA double-strand breaks (DSBs), and (4) cell-cycle modification or cell death. The intercellular signals are generated in and released from irradiated cells. The signals through the MDP and the GJP are modeled independently based on diffusion equations. The irradiation and both signals raise the number of DSBs, which determines transitions of cellular states, such as cell-cycle arrest or cell death. RESULTS: Our model reproduced fairly well previously reported experimental data on the number of DSBs and cell survival curves. We examined how radiation dose and intercellular signaling dynamically affect the cell cycle. The analysis of model dynamics for the bystander cells revealed that the number of arrested cells did not increase linearly with dose. Arrested cells were more efficiently accumulated by the GJP than by the MDP. CONCLUSIONS: We present here a mathematical model that integrates various bystander responses, such as MDP and GJP signaling, DSB induction, cell-cycle arrest, and cell death. Because it simulates spatial and temporal conditions of irradiation and cellular characteristics, our model will be a powerful tool to predict dynamical radiobiological responses of a cellular population in which irradiated and non-irradiated cells co-exist.

Spectrum of Radiation-Induced Clustered Non-DSB Damage - A Monte Carlo Track Structure Modeling and Calculations.

The aim of this report is to present the spectrum of initial radiation-induced cellular DNA damage [with particular focus on non-double-strand break (DSB) damage] generated by computer simulations. The radiation types modeled in this study were monoenergetic electrons (100 eV-1.5 keV), ultrasoft X-ray photons Ck, AlK and TiK, as well as some selected ions including 3.2 MeV/u proton; 0.74 and 2.4 MeV/u helium ions; 29 MeV/u nitrogen ions and 950 MeV/u iron ions. Monte Carlo track structure methods were used to simulate damage induction by these radiation types in a cell-mimetic condition from a single-track action. The simulations took into account the action of direct energy deposition events and the reaction of hydroxyl radicals on atomistic linear B-DNA segments of a few helical turns including the water of hydration. Our results permitted the following conclusions: a. The absolute levels of different types of damage [base damage, simple and complex single-strand breaks (SSBs) and DSBs] vary depending on the radiation type; b. Within each damage class, the relative proportions of simple and complex damage vary with radiation type, the latter being higher with high-LET radiations; c. Overall, for both low- and high-LET radiations, the ratios of the yields of base damage to SSBs are similar, being about 3.0 ± 0.2; d. Base damage contributes more to the complexity of both SSBs and DSBs, than additional SSB damage and this is true for both low- and high-LET radiations; and e. The average SSB/DSB ratio for low-LET radiations is about 18, which is about 5 times higher than that for high-LET radiations. The hypothesis that clustered DNA damage is more difficult for cells to repair has gained currency among radiobiologists. However, as yet, there is no direct in vivo experimental method to validate the dependence of kinetics of DNA repair on DNA damage complexity (both DSB and non-DSB types). The data on the detailed spectrum of DNA damage presented here, in particular the non-DSB type, provide a good basis for testing mechanistic models of DNA repair kinetics such as base excision repair.

A mathematical model of radiation-induced responses in a cellular population including cell-to-cell communications.

Cell-to-cell communication is an important factor for understanding the mechanisms of radiation-induced responses such as bystander effects. In this study, a new mathematical model of intercellular signalling between individual cells in a cellular population is proposed. The authors considered two types of transmission of signals: via culture medium and via gap junction. They focus on the effects that radiation and intercellular signalling have on cell-cycle modification. The cell cycle is represented as a virtual clock that includes several checkpoint pathways within a cyclic process. They also develop a grid model and set up diffusion equations to model the propagation of signals to and from spatially located cells. The authors have also considered the role that DNA damage plays in the cycle of cells which can progress through the cell cycle or stop at the G1, S, G2 or M-phase checkpoints. Results of testing show that the proposed model can simulate intercellular signalling and cell-cycle progression in individual cells during and after irradiation.

Radial dependence of lineal energy distribution of 290-MeV/u carbon and 500-MeV/u iron ion beams using a wall-less tissue-equivalent proportional counter.

Using a wall-less tissue-equivalent proportional counter for a 0.72-µm site in tissue, we measured the radial dependence of the lineal energy distribution, yf(y), of 290-MeV/u carbon ions and 500-MeV/u iron ion beams. The measured yf(y) distributions and the dose-mean of y, [Formula: see text], were compared with calculations performed with the track structure simulation code TRACION and the microdosimetric function of the Particle and Heavy Ion Transport code System (PHITS). The values of the measured [Formula: see text] were consistent with calculated results within an error of 2%, but differences in the shape of yf(y) were observed for iron ion irradiation. This result indicates that further improvement of the calculation model for yf(y) distribution in PHITS is needed for the analytical function that describes energy deposition by delta rays, particularly for primary ions having linear energy transfer in excess of a few hundred keV µm(-1).

Visualisation of γH2AX foci caused by heavy ion particle traversal; distinction between core track versus non-track damage.

Heavy particle irradiation produces complex DNA double strand breaks (DSBs) which can arise from primary ionisation events within the particle trajectory. Additionally, secondary electrons, termed delta-electrons, which have a range of distributions can create low linear energy transfer (LET) damage within but also distant from the track. DNA damage by delta-electrons distant from the track has not previously been carefully characterised. Using imaging with deconvolution, we show that at 8 hours after exposure to Fe (∼200 keV/µm) ions, γH2AX foci forming at DSBs within the particle track are large and encompass multiple smaller and closely localised foci, which we designate as clustered γH2AX foci. These foci are repaired with slow kinetics by DNA non-homologous end-joining (NHEJ) in G1 phase with the magnitude of complexity diminishing with time. These clustered foci (containing 10 or more individual foci) represent a signature of DSBs caused by high LET heavy particle radiation. We also identified simple γH2AX foci distant from the track, which resemble those arising after X-ray exposure, which we attribute to low LET delta-electron induced DSBs. They are rapidly repaired by NHEJ. Clustered γH2AX foci induced by heavy particle radiation cause prolonged checkpoint arrest compared to simple γH2AX foci following X-irradiation. However, mitotic entry was observed when ∼10 clustered foci remain. Thus, cells can progress into mitosis with multiple clusters of DSBs following the traversal of a heavy particle.

Induction of DNA damage, including abasic sites, in plasmid DNA by carbon ion and X-ray irradiation.

DNA from plasmid pUC18 was irradiated with low-LET (13 keV/µm) or high-LET (60 keV/µm) carbon ions or X-rays (4 keV/µm) in solutions containing several concentrations of Tris (0.66-200 mM) to determine the yield of abasic (AP) sites and the effect of scavenging capacity. The yield of AP sites, detected as single-strand breaks (SSB) after digestion with E. coli endonuclease IV (Nfo), was compared with that of SSB and base lesions. At higher concentrations of Tris, the yields of single or clustered AP sites were significantly lower than those of single or clustered base lesions. The relative yields of single AP sites and AP clusters were less than 10 and 7 %, respectively, of the total damage produced at a scavenger capacity mimicking that in cells. The dependence of the yield of AP sites on scavenging capacity was similar to that of prompt strand breaks. The ratios of the yield of isolated AP sites to that of SSB induced by carbon ion or X-ray irradiation were relatively constant at 0.45 ± 0.15 over the tested range of scavenger capacity, although the ratio of SSB to double-strand breaks (DSB) showed the characteristic dependence on both scavenging capacity and radiation quality. These results indicate that the reaction of water radiolysis products, presumably OH radicals, with the sugar-phosphate moieties in the DNA backbone induces both AP sites and SSB with similar efficiency. Direct ionization of DNA is notably more involved in the production of DSB and base lesion clusters than in the production of AP site clusters.

Development of a simulation method for dynamics of electrons ejected from DNA molecules irradiated with X-rays.

PURPOSE: To develop a method for simulating the dynamics of the photoelectrons and Auger electrons ejected from DNA molecules irradiated with pulsed monochromatic X-rays. MATERIALS AND METHODS: A 30-base-pair (bp) DNA molecule was used as the target model, and the X-rays were assumed to have a Gaussian-shaped time distribution. Photoionization and Auger decay were considered as the atomic processes. The atoms from which the photoelectrons or Auger electrons were emitted were specified in the DNA molecule (or DNA ion) using the Monte Carlo method, and the trajectory of each electron in the electric field formed around the positively charged DNA molecule was calculated with a Newtonian equation. The kinetics of the electrons produced by irradiation with X-rays at an intensity ranging from 1 × 10(12) to 1 × 10(16) photons/mm(2) and energies of 380 eV (below the carbon K-edge), 435 eV (above the nitrogen K-edge), and 560 eV (above the oxygen K-edge) were evaluated. RESULTS: It was found that at an X-ray intensity of 1 × 10(14) photons/mm(2) or less, all the produced electrons escaped from the target. However, above an X-ray intensity of 1 × 10(15) photons/mm(2) and an energy of 560 eV, some photoelectrons that were ejected from the oxygen atoms were trapped near the target DNA. CONCLUSIONS: A simulation method for studying the trajectories of electrons ejected from a 30-bp DNA molecule irradiated with pulsed monochromatic X-rays has been developed. The present results show that electron dynamics are strongly dependent on the charged density induced in DNA by pulsed X-ray irradiation.

Applications of the microdosimetric function implemented in the macroscopic particle transport simulation code PHITS.

PURPOSE: Microdosimetric quantities such as lineal energy are generally considered to be better indices than linear energy transfer (LET) for expressing the relative biological effectiveness (RBE) of high charge and energy particles. To calculate their probability densities (PD) in macroscopic matter, it is necessary to integrate microdosimetric tools such as track-structure simulation codes with macroscopic particle transport simulation codes. METHODS: As an integration approach, the mathematical model for calculating the PD of microdosimetric quantities developed based on track-structure simulations was incorporated into the macroscopic particle transport simulation code PHITS (Particle and Heavy Ion Transport code System). The improved PHITS enables the PD in macroscopic matter to be calculated within a reasonable computation time, while taking their stochastic nature into account. APPLICATIONS: The microdosimetric function of PHITS was applied to biological dose estimation for charged-particle therapy and risk estimation for astronauts. The former application was performed in combination with the microdosimetric kinetic model, while the latter employed the radiation quality factor expressed as a function of lineal energy. CONCLUSION: Owing to the unique features of the microdosimetric function, the improved PHITS has the potential to establish more sophisticated systems for radiological protection in space as well as for the treatment planning of charged-particle therapy.