In silico Neutron Relative Biological Effectiveness Estimations For Pre-DNA Repair And Post-DNA Repair Endpoints

A comprehensive understanding of the energy-dependent stochastic risks associated with neutron exposure is crucial to develop robust radioprotection systems. However, the scarcity of experimental data presents significant challenges in this domain. Track-structure Monte Carlo simulations with DNA models have demonstrated their potential to further our fundamental understanding of neutron-induced stochastic risks. To date, most track-structure Monte Carlo studies on the relative biological effectiveness (RBE) of neutrons have focused on various types of DNA damage clusters defined using base pair distances. In this study, we extend these methodologies by incorporating the simulation of non-homologous end joining (NHEJ) DNA repair in order to evaluate the RBE of neutrons for misrepairs. To achieve this, we adapted our previously published Monte Carlo DNA damage simulation pipeline, which combines condensed-history and track-structure Monte Carlo methods, to support the standard DNA damage (SDD) data format. This adaptation enabled seamless integration of neutron-induced DNA damage results with the DNA Mechanistic Repair Simulator (DaMaRiS) toolkit. Additionally, we developed a clustering algorithm that reproduces pre-repair endpoints studied in prior works, as well as novel damage clusters based on Euclidean distances. The neutron RBE for misrepairs obtained in this study exhibits a qualitatively similar shape as the RBE obtained for previously reported pre-repair endpoints. However, it peaks higher, reaching a maximum RBE value of 23(1) at a neutron energy of 0.5 MeV. Furthermore, double-strand break (DSB) pairs with a separation distance of 11 nm better matched the RBE for misrepairs than did the DSB clusters defined with base pair distances that were used in previous in silico neutron RBE studies. These results suggest that there are some features to the spatial distribution of neutron-induced DNA damage that are best characterized either by explicit repair mechanism simulations or by cluster analysis based on Euclidean distances.