TY - JOUR
T1 - Mechanistic modelling of oxygen enhancement ratio of radiation via Monte Carlo simulation-based DNA damage calculation
AU - Lai, Youfang
AU - Chi, Yujie
AU - Jia, Xun
N1 - Funding Information:
This study is supported in part by the Cancer Prevention and Research Institute of Texas (Grant Number RP160661) and by the National Institutes of Health (Grant Number R15CA256668, R37CA214639).
Publisher Copyright:
© 2022 Institute of Physics and Engineering in Medicine.
PY - 2022/9/7
Y1 - 2022/9/7
N2 - 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 of t 0. Finally, we defined two parameters P 0 and P 1, 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 (t 0, P 0 and P 1) 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 as t 0 = 3.8 ms, P 0 = 0.08, and P 1 = 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 of P 0 and P 1, while the curve was subject to a horizontal shift when adjusting t 0. Significance. This study developed a mechanistic model that fully relies on microscopic MC simulations to compute OER.
AB - 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 of t 0. Finally, we defined two parameters P 0 and P 1, 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 (t 0, P 0 and P 1) 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 as t 0 = 3.8 ms, P 0 = 0.08, and P 1 = 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 of P 0 and P 1, while the curve was subject to a horizontal shift when adjusting t 0. Significance. This study developed a mechanistic model that fully relies on microscopic MC simulations to compute OER.
KW - DNA damage
KW - GPU
KW - Monte Carlo
KW - oxygen enhancement ratio
UR - http://www.scopus.com/inward/record.url?scp=85137267885&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85137267885&partnerID=8YFLogxK
U2 - 10.1088/1361-6560/ac8853
DO - 10.1088/1361-6560/ac8853
M3 - Article
C2 - 35944522
AN - SCOPUS:85137267885
SN - 0031-9155
VL - 67
JO - Physics in medicine and biology
JF - Physics in medicine and biology
IS - 17
M1 - 175009
ER -