Nanodosimetric Simulation of Direct Ion-Induced DNA Damage Using Different Chromatin Geometry Models

Monte Carlo based simulation has proven useful in investigating the effect of proton-induced DNA damage and the processes through which this damage occurs. Clustering of ionizations within a small volume can be related to DNA damage through the principles of nanodosimetry. For simulation, it is standard to construct a small volume of water and determine spatial clusters. More recently, realistic DNA geometries have been used, tracking energy depositions within DNA backbone volumes. Traditionally a chromatin fiber is built within the simulation and identically replicated throughout a cell nucleus, representing the cell in interphase. However, the in vivo geometry of the chromatin fiber is still unknown within the literature, with many proposed models. In this work, the Geant4-DNA toolkit was used to build three chromatin models: the solenoid, zig-zag and cross-linked geometries. All fibers were built to the same chromatin density of 4.2 nucleosomes/11 nm. The fibers were then irradiated with protons (LET 5-80 keV/μm) or alpha particles (LET 63-226 keV/μm). Nanodosimetric parameters were scored for each fiber after each LET and used as a comparator among the models. Statistically significant differences were observed in the double-strand break backbone size distributions among the models, although nonsignificant differences were noted among the nanodosimetric parameters. From the data presented in this article, we conclude that selection of the solenoid, zig-zag or cross-linked chromatin model does not significantly affect the calculated nanodosimetric parameters. This allows for a simulation-based cell model to make use of any of these chromatin models for the scoring of direct ion-induced DNA damage.