Modeling gold nanoparticle radiosensitization using a clustering algorithm to quantitate DNA double-strand breaks with mixed-physics Monte Carlo simulation.

Abstract

The radiosensitization properties of gold nanoparticles (GNPs) are investigated using a simple Geant4 cell model considering a realistic cell geometry and a clustering algorithm to characterize the number of DNA double-strand breaks (DSBs). A mixed-physics approach is taken for accurate modeling of low-energy photon interactions in the different regions of the model using Geant4-DNA physics within the cell, and Livermore physics within gold. Density-based spatial clustering of applications with noise (DBSCAN), a clustering algorithm, is used to directly quantitate DNA DSBs after irradiation. The simulation was run using different sizes of GNPs, different distances of GNPs from the cell nucleus, and several combinations of these two conditions. Four types of radiation were simulated in the work: 80-keV monoenergetic photons, 100-keV monoenergetic photons, a 250-kVp photon spectrum, and a 6-MV flattening filter free (FFF) photon spectrum. A variable enhancement in DSB yield, nucleus dose, and cell dose was observed when there are GNPs in the cell cytoplasm, and increases with larger GNPs and proximity to the nucleus. The distance of the GNPs from the nucleus has a large impact on the DSB yield and nucleus dose, but little to no effect on the cell dose. The cell dose enhancement factor of 80keV photons varies from 1.037-1.125 at 0.2µm for 30-100nm GNPs to 1.040-1.127 at 4μm. The DSB enhancement factor varies from 1.050 to 1.174 at 0.2µm to a marginal effect of <1.01 at 4μm. For 100keV, the dose enhancement factor is from 1.142-1.470 at 0.2µm to 1.106-1.371 at 4μm. The DSB enhancement factor varies from 1.249-1.813 at 0.2µm to almost no effect at 4μm. For 250kVp, the dose enhancement factor is from 1.117-1.393 at 0.2μm to 1.110-1.342 at 4μm. The DSB enhancement factor varies from 1.183-1.600 at 0.2μm to a marginal effect of ~1.03 at 4μm. A 6-MV FFF shows a dose enhancement factor of 1.061-1.103 at 0.2μm and 1.053-1.107 at 4μm. The DSB yield varies from 1.070-1.143 at 0.2μm to a marginal effect at 4μm. The stark difference in behavior for DSB yield when compared to cell dose highlights the importance of evaluating more complex radiobiological quantities rather than dose alone when evaluating the radiosensitization properties from metallic nanomaterials. The nucleus dose showed similar characteristics to the DSB yield demonstrating the ability of the method to predict DNA damage and its relationship with nuclear dose. The proposed method provides a way to explore the radiobiological mechanisms of radiation-induced DNA damages, and it aids to evaluate the physical radiosensitization properties of GNP-aided radiotherapy, which can be easily combined with radiochemical DSB quantitation in order to better understand the intricate DNA damage induction mechanisms that are involved in GNP-aided radiotherapy.