Sato, Tatsuhiko

Microdosimetry plays a critical role in particle therapy by quantifying energy deposition within microscopic domains to assess biological effects. This study evaluates the influence of different microdosimetric functions (MFs) and domain geometries (DGs) on relative biological effectiveness (RBE) predictions in carbon ion radiotherapy. Specifically, we compare the analytical microdosimetric function (AMF), calculated for spherical domains and implemented in PHITS, with the Kiefer–Chatterjee (KC) track structure model, which is conventionally applied to cylindrical geometries. To enable a direct comparison, we also introduce a novel implementation of the KC model for spherical domains. Using both models, specific energy distributions were calculated across a range of domain sizes and geometries. These distributions were input into the modified microdosimetric kinetic model (mMKM) to calculate RBE for the HSG cell line and compared against published in vitro data. The results show that both microdosimetric function and domain geometry significantly affect microdosimetric spectra and the resulting RBE, with deviations exceeding 10% when fixed mMKM parameters are used. Parameter optimization within the mMKM enables alignment across models. Our findings emphasize that microdosimetric function and domain geometry selection must be explicitly accounted for in microdosimetry-based RBE modeling, and that model parameters must be tuned accordingly to ensure consistent and biologically accurate predictions.

AbstractBackgroundPhenomenological relative biological effectiveness (RBE) models for proton therapy, based on the dose‐averaged linear energy transfer (LET), have been developed to address the apparent RBE increase towards the end of the proton range. The results of these phenomenological models substantially differ due to varying empirical assumptions and fitting functions. In contrast, more theory‐based approaches are used in carbon ion radiotherapy, such as the microdosimetric kinetic model (MKM). However, implementing microdosimetry‐based models in LET‐based proton therapy treatment planning systems poses challenges.PurposeThis work presents a LET‐based version of the MKM that is practical for clinical use in proton radiotherapy.MethodsAt first, we derived an approximation of the Mayo Clinic Florida (MCF) MKM for relatively‐sparsely ionizing radiation such as protons. The mathematical formalism of the proposed model is equivalent to the original MKM, but it maintains some key features of the MCF MKM, such as the determination of model parameters from measurable cell characteristics. Subsequently, we carried out Monte Carlo calculations with PHITS in different simulated scenarios to establish a heuristic correlation between microdosimetric quantities and the dose averaged LET of protons.ResultsA simple allometric function was found able to describe the relationship between the dose‐averaged LET of protons and the dose‐mean lineal energy, which includes the contributions of secondary particles. The LET‐based MKM was used to model the in vitro clonogenic survival RBE of five human and rodent cell lines (A549, AG01522, CHO, T98G, and U87) exposed to pristine and spread‐out Bragg peak (SOBP) proton beams. The results of the LET‐based MKM agree well with the biological data in a comparable or better way with respect to the other models included in the study. A sensitivity analysis on the model results was also performed.ConclusionsThe LET‐based MKM integrates the predictive theoretical framework of the MCF MKM with a straightforward mathematical description of the RBE based on the dose‐averaged LET, a physical quantity readily available in modern treatment planning systems for proton therapy.

The analytical microdosimetric function (AMF) implemented in the Monte Carlo code PHITS is a unique tool that bridges the gap between macro- and microscopic scales of radiation interactions, enabling accurate microdosimetric calculations over macroscopic bodies. The original AMF was published in 2006, based on the results of track structure calculations. Recently, a newer version of the AMF was proposed, incorporating an improved description of the energy loss at the microscopic scale. This study compares the older and the newer AMFs in computing microdosimetric probability distributions, mean values, and the relative biological effectiveness (RBE). To this end, 16000 microdosimetric lineal energy probability density distributions were simulated with PHITS for ions from 1H to 238U over a broad energy range (1–1000 MeV/n). The newer AMF was found to offer superior performance, particularly for very heavy ions, producing results that align more closely with published in vitro clonogenic survival experiments. These findings suggest that the updated AMF provides a more reliable tool for microdosimetric calculations and RBE modeling, essential for ion radiation therapy and space radiation protection.

Abstract Objective . Estimation of the probability density of the microdosimetric quantities in macroscopic matter is indispensable for applying the concept of microdosimetry to medical physics and radiological protection. The Particle and Heavy Ion Transport code System (PHITS) enables estimating the microdosimetric probability densities due to its unique hybrid modality between the Monte Carlo and analytical approaches called the microdosimetric function. It can convert the deposition energies calculated by the macroscopic Monte Carlo radiation transport simulation to microdosimetric probability densities in water using an analytical function based on the track-structure simulations. Approach . In this study, we improved this function using the latest track-structure simulation codes implemented in PHITS. The improved function is capable of calculating the probability densities of not only the conventional microdosimetric quantities such as lineal energy but also the number of ionization events occurring in a target site, the so-called ionization cluster size distribution, for arbitrary site diameters from 3 nm to 1 μ m. Main results . The accuracy of the improved function was well verified by comparing the microdosimetric probability densities measured by tissue-equivalent proportional counters with the corresponding data calculated in this study. Test calculations for clonogenic cell survival using the improved function coupled with the modified microdosimetric kinetic model suggested a slight increase of its relative biological effectiveness compared with our previous estimations. As a new application of the improved function, we calculated the relative biological effectiveness of the single-strand break and double-strand break yields for proton irradiations using the updated PHITS coupled with the simplified DNA damage estimation model, and confirmed its equivalence in accuracy and its superiority in computational time compared to our previously proposed method based on the track-structure simulation. Significance . From these features, we concluded that the improved function could expand the application fields of PHITS by bridging the gap between microdosimetry and macrodosimetry.

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