Monte Carlo track-structure for the radionuclide Copper-64: characterization of S-values, nanodosimetry and quantification of direct damage to DNA

TOPAS-nBio was used to simulate, collision-to-collision, the complete trajectories of electrons in water generated during the explicit simulation of ⁶⁴Cu decay. S-values and direct damage to the DNA were calculated representing the cell (C) and the cell nucleus (N) with concentric spheres of 5 μm and 4 μm in radius, respectively. The considered 'target'←'source' configurations, including the cell surface (Cs) and cytoplasm (Cy), were: C←C, C←Cs, N←N, N←Cy and N←Cs. Ionization cluster size distributions were also calculated in a cylinder immersed in water corresponding to a DNA segment of 10 base-pairs in length (diameter 2.3 nm, length 3.4 nm), modeling a radioactive point source moving from the central axis to the edge of the cylinder. For that, the first moment (M₁) and cumulative probability of having a cluster size of 2 or more ionizations in the cylindrical volume (F₂) were obtained. Finally, the direct damage to the DNA was estimated by quantifying double-strand breaks (DSBs) using the clustering algorithm DBSCAN.

The S-values obtained with TOPAS-nBio for ⁶⁴Cu were 7.879 × 10⁻⁴ ± 5 × 10⁻⁷, 4.351 × 10⁻⁴ ± 6 × 10⁻⁷, 1.442 × 10⁻³ ± 1 × 10⁻⁶, 2.596 × 10⁻⁴ ± 8 × 10⁻⁷, 1.127 × 10⁻⁴ ± 4 × 10⁻⁷ Gy Bq-s⁻¹ for the configurations C←C, C←Cs, N←N, N←Cy and N←Cs, respectively. The difference of these values, compared with previously reported S-values for ⁶⁴Cu with the code MNCP and software MIRDCell, ranged from −4% to −25% for the configurations N←N and N←Cs, respectively. On the other hand, F₂ was maximum with the source at the center of the cylinder 0.373 ± 0.001, and monotonically decreased until reaching a value of 0.058 ± 0.001 at 2.3 nm. The same behavior was observed for M₁ with values ranging from 2.188 ± 0.004 to 0.242 ± 0.002. Finally, the DBSCAN algorithm showed that the mean number of DNA DSBs per decay were 0.187 ± 0.001, 0.0317 ± 0.0005, and 0.0125 ± 0.0002 DSB-(Bq-s)⁻¹ for the configurations N←N, N←Cs, and N←Cy, respectively.

In conclusion, the results of the S-values show that the absorbed dose strongly depends on the distribution of the radionuclide in the cell, the dose being higher when ⁶⁴Cu is internalized in the cell nucleus, which is reinforced by the nanodosimetric study by the presence of DNA DSBs attributable to the Auger electrons emitted during the decay of ⁶⁴Cu. (paper)