Implementation of simplified stochastic microdosimetric kinetic models into PHITS for application to radiation treatment planning

Radiation research trends by young scientists and the future tasks in Northern Japan: report on 'the 10th educational symposium on radiation and health (ESRAH) by young scientists in 2023'

PURPOSE

Since 2014, an educational activity on radiation and health in northern Japan has been carried out by young scientists, the so-called 'Educational Symposium on Radiation and Health (ESRAH)'. Close cooperation has been continued in preparing for any possible emergency response to radiation accidents because several facilities, e.g., the Tomari Nuclear Power Plant in Hokkaido and the Low-Level Radioactive Waste Disposal Facility in Aomori prefecture. The ESRAH meeting has provided informational exchange and discussion forum on a broad range of subjects in various. In 2023, the 10th Memorial ESRAH meeting took place to boost scientific understanding and multidisciplinary collaborations for young scientists. Herein, we report on the ESRAH2023 symposium and analyze the research categories of young scientists from the past 10-year presentations.

CONCLUSIONS

To date, the ESRAH meeting has successfully provided a chance for multi-disciplinary research, which accounted for 27% of the total despite the COVID-19 pandemic. We found that the fraction of multi-disciplinary research in 2023 was the highest during 10-year ESRAH meetings. Meanwhile, amongst the research categories, physics, chemistry, and pharmacological studies were indicated to be less for young scientists. It is desired that further collaboration between physics, chemistry, and pharmacology in addition to the current fields would not only clarify radiation effects on the human body but also promote emergency medical care for radiation exposure in the future.

An automation of Monte Carlo workflow for dosimetry study of an Elekta LINAC delivery system in radiotherapy

Purpose

This study aims to automate the Monte Carlo (MC) workflow utilized for radiotherapy dosimetry, focusing on an Elekta LINAC delivery system. It addresses the challenge of integrating MC simulations into routine clinical practice, making this accurate yet complex method more accessible and efficient for radiotherapy dosimetry.

Methods and Materials

We developed a user-friendly software featuring a graphical user interface (GUI) that integrates EGSnrc for MC simulations. The software streamlines the process from retrieving Digital Imaging and Communications in Medicine (DICOM) data to executing dose calculations and comparing dose distributions. To validate our proposed tool, we compared its computed doses for IMRT and VMAT plans from the Pinnacle TPS for an Elekta Versa HD linear accelerator against MC simulation results. This comparison utilized our in-house software and GUI as the tool, covering various treatment sites and prescriptions.

Results

The automated MC workflow demonstrated high accuracy in dose calculations and streamlined integration with clinical workflows. The comparison between the MC-simulated and TPS-calculated doses revealed excellent agreement, highlighting the reliability of MC for independent dose verification in complex treatment scenarios.

Conclusions

The automated MC workflow developed represents a substantial improvement in the practicality and efficiency of MC simulations in radiotherapy. This advancement not only simplifies the dosimetry process but also ensures high accuracy, establishing it as a valuable tool for routine patient-specific quality assurance and the development of specialized treatment procedures.

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

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

Application of a simple DNA damage model developed for electrons to proton irradiation

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.

The Mayo Clinic Florida microdosimetric kinetic model of clonogenic survival: formalism and first benchmark against in vitro and in silico data

Objective. To develop a new model (Mayo Clinic Florida microdosimetric kinetic model, MCF MKM) capable of accurately describing the in vitro clonogenic survival at low and high linear energy transfer (LET) using single-event microdosimetric spectra in a single target. Methodology. The MCF MKM is based on the ‘post-processing average’ implementation of the non-Poisson microdosimetric kinetic model and includes a novel expression to compute the particle-specific quadratic-dependence of the cell survival with respect to dose (β of the linear-quadratic model). A new methodology to a priori calculate the mean radius of the MCF MKM subnuclear domains is also introduced. Lineal energy spectra were simulated with the Particle and Heavy Ion Transport code System (PHITS) for 1H, 4He, 12C, 20Ne, 40Ar, 56Fe, and 132Xe ions and used in combination with the MCF MKM to calculate the ion-specific LET-dependence of the relative biological effectiveness (RBE) for Chinese hamster lung fibroblasts (V79 cell line) and human salivary gland tumor cells (HSG cell line). The results were compared with in vitro data from the Particle Irradiation Data Ensemble (PIDE) and in silico results of different models. The possibility of performing experiment-specific predictions to explain the scatter in the in vitro RBE data was also investigated. Finally, a sensitivity analysis on the model parameters is also included. Main results. The RBE values predicted with the MCF MKM were found to be in good agreement with the in vitro data for all tested conditions. Though all MCF MKM model parameters were determined a priori, the accuracy of the MCF MKM was found to be comparable or superior to that of other models. The model parameters determined a priori were in good agreement with the ones obtained by fitting all available in vitro data. Significance. The MCF MKM will be considered for implementation in cancer radiotherapy treatment planning with accelerated ions.