Recent Developments on gMicroMC: Transport Simulations of Proton and Heavy Ions and Concurrent Transport of Radicals and DNA

High anhydrous proton conductivity and a smart proton transportation approach of a sulfate coordination polymer

We successfully synthesized a one-dimensional cobalt sulfate coordinating polymer, whose simple hydrogen bond web structure facilitated the analysis of the proton transfer process. At 175 °C, without humidity, the conductivity is 0.0311 S cm-1, which exceeds those of most of the organic inorganic hybrid materials under anhydrous conditions (world record rank 8). Based on its crystal structure and theoretical calculations, the subversive proton conduction pathway was inferred clearly. We, for the first time, found that the proton smartly chose the path with a lower energy barrier but not the one with short distance to transport avoiding short circuit.

Mechanistic modelling of relative biological effectiveness of carbon ion beams and comparison with experiments

Objective.Relative biological effectiveness (RBE) plays a vital role in carbon ion radiotherapy, which is a promising treatment method for reducing toxic effects on normal tissues and improving treatment efficacy. It is important to have an effective and precise way of obtaining RBE values to support clinical decisions. A method of calculating RBE from a mechanistic perspective is reported.Approach.Ratio of dose to obtain the same number of double strand breaks (DSBs) between different radiation types was used to evaluate RBE. Package gMicroMC was used to simulate DSB yields. The DSB inductions were then analyzed to calculate RBE. The RBE values were compared with experimental results.Main results.Furusawa's experiment yielded RBE values of 1.27, 2.22, 3.00 and 3.37 for carbon ion beam with dose-averaged LET of 30.3 keVμm-1, 54.5 keVμm-1, 88 keVμm-1and 137 keVμm-1, respectively. RBE values computed from gMicroMC simulations were 1.75, 2.22, 2.87 and 2.97. When it came to a more sophisticated carbon ion beam with 6 cm spread-out Bragg peak, RBE values were 1.61, 1.63, 2.19 and 2.36 for proximal, middle, distal and distal end part, respectively. Values simulated by gMicroMC were 1.50, 1.87, 2.19 and 2.34. The simulated results were in reasonable agreement with the experimental data.Significance.As a mechanistic way for the evaluation of RBE for carbon ion radiotherapy by combining the macroscopic simulation of energy spectrum and microscopic simulation of DNA damages, this work provides a promising tool for RBE calculation supporting clinical applications such as treatment planning.

Development and implementation of a metaphase DNA model for ionizing radiation induced DNA damage calculation

Objective. To develop a metaphase chromosome model representing the complete genome of a human lymphocyte cell to support microscopic Monte Carlo (MMC) simulation-based radiation-induced DNA damage studies.Approach. We first employed coarse-grained polymer physics simulation to obtain a rod-shaped chromatid segment of 730 nm in diameter and 460 nm in height to match Hi-C data. We then voxelized the segment with a voxel size of 11 nm per side and connected the chromatid with 30 types of pre-constructed nucleosomes and 6 types of linker DNAs in base pair (bp) resolutions. Afterward, we piled different numbers of voxelized chromatid segments to create 23 pairs of chromosomes of 1-5μm long. Finally, we arranged the chromosomes at the cell metaphase plate of 5.5μm in radius to create the complete set of metaphase chromosomes. We implemented the model in gMicroMC simulation by denoting the DNA structure in a four-level hierarchical tree: nucleotide pairs, nucleosomes and linker DNAs, chromatid segments, and chromosomes. We applied the model to compute DNA damage under different radiation conditions and compared the results to those obtained with G0/G1 model and experimental measurements. We also performed uncertainty analysis for relevant simulation parameters.Main results. The chromatid segment was successfully voxelized and connected in bps resolution, containing 26.8 mega bps (Mbps) of DNA. With 466 segments, we obtained the metaphase chromosome containing 12.5 Gbps of DNA. Applying it to compute the radiation-induced DNA damage, the obtained results were self-consistent and agreed with experimental measurements. Through the parameter uncertainty study, we found that the DNA damage ratio between metaphase and G0/G1 phase models was not sensitive to the chemical simulation time. The damage was also not sensitive to the specific parameter settings in the polymer physics simulation, as long as the produced metaphase model followed a similar contact map distribution.Significance. Experimental data reveal that ionizing radiation induced DNA damage is cell cycle dependent. Yet, DNA chromosome models, except for the G0/G1 phase, are not available in the state-of-the-art MMC simulation. For the first time, we successfully built a metaphase chromosome model and implemented it into MMC simulation for radiation-induced DNA damage computation.

A potential revolution in cancer treatment: A topical review of FLASH radiotherapy

FLASH radiotherapy (RT) is a novel technique in which the ultrahigh dose rate (UHDR) (≥40 Gy/s) is delivered to the entire treatment volume. Recent outcomes of in vivo studies show that the UHDR RT has the potential to spare normal tissue without sacrificing tumor control. There is a growing interest in the application of FLASH RT, and the ultrahigh dose irradiation delivery has been achieved by a few experimental and modified linear accelerators. The underlying mechanism of FLASH effect is yet to be fully understood, but the oxygen depletion in normal tissue providing extra protection during FLASH irradiation is a hypothesis that attracts most attention currently. Monte Carlo simulation is playing an important role in FLASH, enabling the understanding of its dosimetry calculations and hardware design. More advanced Monte Carlo simulation tools are under development to fulfill the challenge of reproducing the radiolysis and radiobiology processes in FLASH irradiation. FLASH RT may become one of standard treatment modalities for tumor treatment in the future. This paper presents the history and status of FLASH RT studies with a focus on FLASH irradiation delivery modalities, underlying mechanism of FLASH effect, in vivo and vitro experiments, and simulation studies. Existing challenges and prospects of this novel technique are discussed in this manuscript.

Mechanistic modelling of oxygen enhancement ratio of radiation via Monte Carlo simulation-based DNA damage calculation

Objective.Oxygen plays an important role in affecting the cellular radio-sensitivity to ionizing radiation. The objective of this study is to build a mechanistic model to compute oxygen enhancement ratio (OER) using a GPU-based Monte Carlo (MC) simulation package gMicroMC for microscopic radiation transport simulation and DNA damage calculation.Approach.We first simulated the water radiolysis process in the presence of DNA and oxygen for 1 ns and recorded the produced DNA damages. In this process, chemical reactions among oxygen, water radiolysis free radicals and DNA molecules were considered. We then applied a probabilistic approach to model the reactions between oxygen and indirect DNA damages for a maximal reaction time oft0. Finally, we defined two parametersP0andP1, representing probabilities for DNA damages without and with oxygen fixation effect not being restored in the repair process, to compute the final DNA double strand breaks (DSBs). As cell survival fraction is mainly determined by the number of DSBs, we assumed that the same numbers of DSBs resulted in the same cell survival rates, which enabled us to compute the OER as the ratio of doses producing the same number of DSBs without and with oxygen. We determined the three parameters (t0,P0andP1) by fitting the OERs obtained in our computation to a set of published experimental data under x-ray irradiation. We then validated the model by performing OER studies under proton irradiation and studied model sensitivity to parameter values.Main results.We obtained the model parameters ast0= 3.8 ms,P0= 0.08, andP1= 0.28 with a mean difference of 3.8% between the OERs computed by our model and that obtained from experimental measurements under x-ray irradiation. Applying the established model to proton irradiation, we obtained OERs as functions of oxygen concentration, LET, and dose values, which generally agreed with published experimental data. The parameter sensitivity analysis revealed that the absolute magnitude of the OER curve relied on the values ofP0andP1, while the curve was subject to a horizontal shift when adjustingt0.Significance.This study developed a mechanistic model that fully relies on microscopic MC simulations to compute OER.