Electron track structure simulations in a gold nanoparticle using Geant4-DNA

The role of coating layers in gold nanorods' radioenhancement: a Monte Carlo analysis

Gold nanoparticles are promising radiosensitizing agents for nanoparticle-enhanced radiotherapy (NPRT). The coating layer on these nanoparticles can significantly influence their physicochemical characteristics and biological behavior. This study investigates the influence of various coating layers on the radioenhancement efficiency of gold nanorods by modeling the physical interactions and chemical reactions involved. We conducted Monte Carlo simulations using the TOPAS code to model secondary electron generation in gold nanorods exposed to 100 kVp X-rays. Through a multiscale approach, the dose contribution, electron spectrum, and G-values of radiolytic species were determined. Four conventional coating materials were examined and compared to a non-coated nanorod. The simulation results indicate that the addition of coating layers decreases the additional dose due to the gold nanorods by up to 7% across all materials. The assessment of electron spectra revealed that 1% to 8% of electrons with energies below 3.5 keV were absorbed within the various coating layers. In contrast, higher-energy electrons were mainly unaffected. The total G-values for all radiolytic species remained generally unchanged with the addition of the coating layer, regardless of the material used. However, increasing the coating thickness slightly increased the relative yield of chemical species at times beyond 10 ns post-irradiation. While the addition of a coating layer generally resulted in a decrease in electron fluence and dose contribution, the reduction was not as substantial as expected from results previously reported in the literature. This suggests that, from the physics perspective, the influence of coating layers on radioenhancement may be less pronounced than previously believed. Additionally, the observed increase in total G-values with thicker coatings emphasizes the need for further investigation into the effects of coatings on radiolytic species.

Nanoscale dosimetry for a radioisotope-labeled metal nanoparticle using MCNP6.2 and Geant4

BACKGROUND

Metal nanoparticles (MNPs) labeled with radioisotopes (RIs) are utilized as radio-enhancers due to their ability to amplify the radiation dose in their immediate vicinity. A thorough understanding of nanoscale dosimetry around MNPs enables their effective application in radiotherapy. However, nanoscale dosimetry around MNPs still requires further investigation.

PURPOSE

This study aims to provide insight into the radio-enhancement effects of MNPs by elucidating nanoscale dosimetry surrounding MNPs labeled with Auger-emitting RIs. We particularly focus on distinguishing the respective dose contributions of photons and electrons emitted by Auger-emitting RIs in the context of dose enhancement.

METHODS

A 50 nm diameter NP of silver (Ag) core and gold (Au) shell (Ag@Au NP) was assumed to emit mono-energetic electrons and photons (3, 5, 10, 20, and 30 keV), or the energy spectrum corresponding to one of three Auger-emitting RIs (103Pd, 125I, and 131Cs) from the Ag core. Nanoscale radial dose distributions around a single radioactive Ag@Au NP were evaluated in spherical shells of water. Monte Carlo simulations were conducted using single-event and track structure transport methods implemented in MCNP6.2 and Geant4-DNA-Au physics, respectively. To evaluate the extent of radio-enhancement by the Ag@Au NP, two scenarios were considered: Ag@Au NPs (Au shell included) and Ag@water NPs (Au shell replaced by water).

RESULTS

The radial doses of 10, 20, and 30 keV electrons estimated by both codes were comparable. However, the radial doses of 3 and 5 keV electrons by MCNP6.2 were much larger near the NP surface than those by Geant4. There was a dose enhancement of a few % to tens % by the Au shell in the region of the NP surface to 10 µm, depending on the electron energy. The radial doses of photons with the Au shell were higher up to their secondary electron ranges than those without the Au shell. The maximum dose enhancement factor of photons occurred at 20 keV and was 63.4 by MCNP6.2 and 50.5 by Geant4. The overall radial doses of electrons were 1-2 orders of magnitude larger than those of photons. As a result, in cases of RIs emitting both electrons and photons, the radial doses up to electron ranges were dominantly governed by electrons. The dose enhancement estimated by both codes for the RIs ranged from a few % except in the immediate vicinity of the NP surface.

CONCLUSION

Given the dominant contribution of electrons to radial doses of MNP labeled with Auger-emitting RIs, physical dose enhancement expected by interactions with photons was hindered. Since there are no available RIs emitting exclusively photons, achieving enhanced physical doses within a cell through a combination of MNPs and RIs appears currently unattainable. The radial doses of photons near the NP surface exhibited considerable discrepancies between the codes, primarily attributed to low-energy electrons. The difference may arise from higher cross-sections of Au inelastic scattering in Geant4-DNA-Au compared to MCNP6.2.

Radiosensitization with metallic nanoparticles under MeV proton beams: local dose enhancement

In addition to specific dosimetric properties of protons, their higher biological effectiveness makes them superior to X-rays and gamma radiation, in radiation therapy. In recent years, enrichment of tumours with metallic nanoparticles as radiosensitizer agents has generated high interest, with several studies attempting to confirm the efficacy of nanoparticles in proton therapy. In the present study Geant4 Monte Carlo (MC) code was used to quantify the increased nanoscopic dose deposition of 50 nm metallic nanoparticles including gold, bismuth, iridium, and gadolinium in water upon exposure to 5, 25, and 50 MeV protons. Dose enhancement factors, radial dose distributions in nano-scale, as well as secondary electron and photon energy spectra were calculated for the studied nanoparticles and proton beams. The obtained results demonstrated that in the presence of metallic nanoparticles an increase in proton energy leads to a decrease in secondary electron and photon production yield. Additionally, an increase in the radial dose enhancement factor from 1.4 to 16 was calculated for the studied nanoparticles when the proton energy was increased from 5 to 50 MeV. It is concluded that the dosimetric advantages of proton beams could be improved significantly in the presence of metallic nanoparticles.

In silico estimation of polyethylene glycol coating effect on metallic NPs radio-sensitization in kilovoltage energy beams

PURPOSE

Nanoparticles (NPs) as radiosensitizers present a promising strategy for enhancing radiotherapy effectiveness, but their potential is significantly influenced by the properties of their surface coating, which can impact treatment outcomes. Most Monte Carlo studies have focused on metallic NPs without considering the impact of coating layers on radiosensitization. In this study, we aim to assess both the physical and radiobiological effects of nanoparticle coatings in nanoparticle-based radiation therapy.

MATERIALS AND METHODS

In this simulation study, we used Geant4 Monte Carlo (MC) toolkit (v10.07.p02) and simulated the bismuth, gold, iridium and gadolinium NPs coated with polyethylene glycol (PEG-400: Density: 1.13 g/cm³, Molar mass: 380-420 g/mol) as radiosensitizer for photon beams of 30, 60 and 100 keV. Secondary electron number and reactive oxygen species enhancement factor were estimated. Also, dose enhancement factor (DEF) was determined in spherical shells with logarithmic scale thickness from the nanoparticle surface to 4 mm.

RESULTS

Secondary electron emission was highest at 30 keV for gold, bismuth, and iridium NPs, while gadolinium NPs peaked at 60 keV. Coating reduced electron emissions across all energies, with thicker coatings leading to a more significant decrease. DEF values declined with increasing radial distance from the NP surface and were lower with thicker coatings. For gadolinium NPs, DEF behavior differed due to K-edge energy effects. Reactive species generation varied, showing maximum production at 30 keV for gold, bismuth, and iridium NPs, while gadolinium NPs showed peak activity at 60 keV. PEG coatings enhanced reactive species formation at 100 keV.

CONCLUSION

The findings indicate that the coating layer thickness and material not only influence the emission of secondary particles and DEF but also affect the generation of reactive species from water radiolysis. Specifically, thicker coatings reduce secondary particle emission and DEF, while PEG coatings demonstrate a dual behavior, offering both protective and enhancing effects depending on photon energy. These insights underscore the importance of optimizing NP design and coating in future studies to maximize therapeutic efficacy in nanoparticle-based radiation therapy.

Radial dependence of ionization clustering around a gold nanoparticle irradiated by X-rays under charged particle equilibrium

Objective.This work explores the enhancement of ionization clustering and its radial dependence around a gold nanoparticle (NP), indicative of the induction of DNA lesions, a potential trigger for cell-death.Approach.Monte Carlo track structure simulations were performed to determine (a) the spectral fluence of incident photons and electrons in water around a gold NP under charged particle equilibrium conditions and (b) the density of ionization clusters produced on average as well as conditional on the occurrence of at least one interaction in the NP using Associated Volume Clustering. Absorbed dose was determined for comparison with a recent benchmark intercomparison. Reported quantities are normalized to primary fluence, allowing to establish a connection to macroscopic dosimetric quantities.Main results.The modification of the electron spectral fluence by the gold NP is minor and mainly occurs at low energies. The net fluence of electrons emitted from the NP is dominated by electrons resulting from photon interactions. Similar to the known dose enhancement, increased ionization clustering is limited to a distance from the NP surface of up to200nm. The number of clusters per energy imparted is increased at distances of up to150nm, and accordingly the enhancement in clustering notably surpasses that of dose enhancement. Smaller NPs cause noticeable peaks in the conditional frequency of clusters between50nm-100nmfrom the NP surface.Significance.This work shows that low energy electrons emitted by NPs lead to an increase of ionization clustering in their vicinity exceeding that of energy imparted. While the electron component of the radiation field plays an important role in determining the background contribution to ionization clustering and energy imparted, the dosimetric effects of NPs are governed by the interplay of secondary electron production by photon interaction and their ability to leave the NP.

Microdosimetric Simulation of Gold-Nanoparticle-Enhanced Radiotherapy

Conventional X-ray therapy (XRT) is commonly applied to suppress cancerous tumors; however, it often inflicts collateral damage to nearby healthy tissue. In order to provide a better conformity of the dose distribution in the irradiated tumor, proton therapy (PT) is increasingly being used to treat solid tumors. Furthermore, radiosensitization with gold nanoparticles (GNPs) has been extensively studied to increase the therapeutic ratio. The mechanism of radiosensitization is assumed to be connected to an enhancement of the absorbed dose due to huge photoelectric cross-sections with gold. Nevertheless, numerous theoretical studies, mostly based on Monte Carlo (MC) simulations, did not provide a consistent and thorough picture of dose enhancement and, therefore, the radiosensitization effect. Radiosensitization by nanoparticles in PT is even less studied than in XRT. Therefore, we investigate the physics picture of GNP-enhanced RT using an MC simulation with Geant4 equipped with the most recent physics models, taking into account a wide range of physics processes relevant for realistic PT and XRT. Namely, we measured dose enhancement factors in the vicinity of GNP, with diameters ranging from 10 nm to 80 nm. The dose enhancement in the vicinity of GNP reaches high values for XRT, while it is very modest for PT. The macroscopic dose enhancement factors for realistic therapeutic GNP concentrations are rather low for all RT scenarios; therefore, other physico-chemical and biological mechanisms should be additionally invoked for an explanation of the radiosensitization effect observed in many experiments.

Review of chemical models and applications in Geant4-DNA: Report from the ESA BioRad III Project

A chemistry module has been implemented in Geant4-DNA since Geant4 version 10.1 to simulate the radiolysis of water after irradiation. It has been used in a number of applications, including the calculation of G-values and early DNA damage, allowing the comparison with experimental data. Since the first version, numerous modifications have been made to the module to improve the computational efficiency and extend the simulation to homogeneous kinetics in bulk solution. With these new developments, new applications have been proposed and released as Geant4 examples, showing how to use chemical processes and models. This work reviews the models implemented and application developments for modeling water radiolysis in Geant4-DNA as reported in the ESA BioRad III Project.

In vitro validation of helium ion irradiations as a function of linear energy transfer in radioresistant human malignant cells

PURPOSE

Based on considerable interest to enlarge the experimental database of radioresistant cells after their irradiation with helium ions, HTB140, MCF-7 and HTB177 human malignant cells are exposed to helium ion beams having different linear energy transfer (LET).

MATERIALS AND METHODS

The cells are irradiated along the widened 62 MeV/u helium ion Bragg peak, providing LET of 4.9, 9.8, 23.4 and 36.8 keV/µm. Numerical simulations with the Geant4 toolkit are used for the experimental design. Cell survival is evaluated and compared with reference γ-rays. DNA double strand breaks are assessed via γ-H2AX foci.

RESULTS

With the increase of LET, surviving fractions at 2 Gy decrease, while RBE (2 Gy, γ) gradually increase. For HTB140 cells, above the dose of 4 Gy, a slight saturation of survival is observed while the increase of RBE (2 Gy, γ) remains unaffected. With the increase of LET the increase of γ-H2AX foci is revealed at 0.5 h after irradiation. There is no significant difference in the number of foci between the cell lines for the same LET. From 0.5 to 24 h, the number of foci drops reaching its residual level. For each time point, there are small differences in DNA DSB among the three cell lines.

CONCLUSION

Analyses of data acquired for the three cell lines irradiated by helium ions, having different LET, reveal high elimination capacity and creation of a large number of DNA DSB with respect to γ-rays, and are between those reported for protons and carbon ions.

XMEA: A New Hybrid Diamond Multielectrode Array for the In Situ Assessment of the Radiation Dose Enhancement by Nanoparticles

This work presents a novel multielectrode array (MEA) to quantitatively assess the dose enhancement factor (DEF) produced in a medium by embedded nanoparticles. The MEA has 16 nanocrystalline diamond electrodes (in a cell-culture well), and a single-crystal diamond divided into four quadrants for X-ray dosimetry. DEF was assessed in water solutions with up to a 1000 µg/mL concentration of silver, platinum, and gold nanoparticles. The X-ray detectors showed a linear response to radiation dose (r2 ≥ 0.9999). Overall, platinum and gold nanoparticles produced a dose enhancement in the medium (maximum of 1.9 and 3.1, respectively), while silver nanoparticles produced a shielding effect (maximum of 37%), lowering the dose in the medium. This work shows that the novel MEA can be a useful tool in the quantitative assessment of radiation dose enhancement due to nanoparticles. Together with its suitability for cells' exocytosis studies, it proves to be a highly versatile device for several applications.

Clonogenic assay and computational modeling using real cell images to study physical enhancement and cellular sensitization induced by metal nanoparticles under MV and kV X-ray irradiation

This study was initiated due to the physically unexplainable tumor controls resulting from metal nanoparticle (MNP) experiments even under MV X-ray irradiation. A more accurate explanation of the mechanism of radiosensitization induced by MNP is warranted, considering both its physical dose enhancement and biological sensitization, as related research is lacking. Thus, we aimed to examine the intricate dynamics involved in MNP-induced radiosensitization. We conducted specifically designed clonogenic assays for the A549 lung cancer cell line with MNP irradiated by 6 MV and 300 kVp X-rays. Two types of MNP were employed: one based on iron oxide, promoting ferroptosis, and the other on gold nanoparticles known for inducing a significant dose enhancement, particularly at low-energy X-rays. We introduced the lethality enhancement factor (LEF) as the fraction in the cell killing attributed to biological sensitization. Subsequently, Monte Carlo simulations were conducted to evaluate the radial dose profiles for each MNP, corresponding to the physical enhancement. Finally, the local effect model was applied to the clonogenic assay results on real cell images. The LEF and the dose enhancement in the cytoplasm were incorporated to increase the accuracy in the average lethal events and, consequently, in the survival fraction. The results reveal an increased cell killing for both of the MNP under MV and kV X-ray irradiation. In both types of MNP, the LEF reveals a biological sensitization evident. The sensitizer enhancement ratio, derived from the calculations, exhibited only 3% and 1% relative differences compared to the conventional linear-quadratic model for gold and ferroptosis inducer nanoparticles, respectively. These findings indicate that MNPs sensitize cells via radiation through mechanisms akin to ferroptosis inducers, not exclusively relying on a physical dose enhancement. Their own contributions to survival fractions were successfully integrated into computational modeling.

Development and implementation of a metaphase DNA model for ionizing radiation induced DNA damage calculation

Objective. To develop a metaphase chromosome model representing the complete genome of a human lymphocyte cell to support microscopic Monte Carlo (MMC) simulation-based radiation-induced DNA damage studies.Approach. We first employed coarse-grained polymer physics simulation to obtain a rod-shaped chromatid segment of 730 nm in diameter and 460 nm in height to match Hi-C data. We then voxelized the segment with a voxel size of 11 nm per side and connected the chromatid with 30 types of pre-constructed nucleosomes and 6 types of linker DNAs in base pair (bp) resolutions. Afterward, we piled different numbers of voxelized chromatid segments to create 23 pairs of chromosomes of 1-5μm long. Finally, we arranged the chromosomes at the cell metaphase plate of 5.5μm in radius to create the complete set of metaphase chromosomes. We implemented the model in gMicroMC simulation by denoting the DNA structure in a four-level hierarchical tree: nucleotide pairs, nucleosomes and linker DNAs, chromatid segments, and chromosomes. We applied the model to compute DNA damage under different radiation conditions and compared the results to those obtained with G0/G1 model and experimental measurements. We also performed uncertainty analysis for relevant simulation parameters.Main results. The chromatid segment was successfully voxelized and connected in bps resolution, containing 26.8 mega bps (Mbps) of DNA. With 466 segments, we obtained the metaphase chromosome containing 12.5 Gbps of DNA. Applying it to compute the radiation-induced DNA damage, the obtained results were self-consistent and agreed with experimental measurements. Through the parameter uncertainty study, we found that the DNA damage ratio between metaphase and G0/G1 phase models was not sensitive to the chemical simulation time. The damage was also not sensitive to the specific parameter settings in the polymer physics simulation, as long as the produced metaphase model followed a similar contact map distribution.Significance. Experimental data reveal that ionizing radiation induced DNA damage is cell cycle dependent. Yet, DNA chromosome models, except for the G0/G1 phase, are not available in the state-of-the-art MMC simulation. For the first time, we successfully built a metaphase chromosome model and implemented it into MMC simulation for radiation-induced DNA damage computation.

Evaluation of Dosimetric Effect of Bone Scatter on Nanoparticle-Enhanced Orthovoltage Radiotherapy: A Monte Carlo Phantom Study

In nanoparticle (NP)-enhanced orthovoltage radiotherapy, bone scatter affected dose enhancement at the skin lesion in areas such as the forehead, chest wall, and knee. Since each of these treatment sites have a bone, such as the frontal bone, rib, or patella, underneath the skin lesion and this bone is not considered in dose delivery calculations, uncertainty arises in the evaluation of dose enhancement with the addition of NPs in radiotherapy. To investigate the impact of neglecting the effect of bone scatter, Monte Carlo simulations based on heterogeneous phantoms were carried out to determine and compare the dose enhancement ratio (DER), when a bone was and was not present underneath the skin lesion. For skin lesions with added NPs, Monte Carlo simulations were used to calculate the DER values using different elemental NPs (gold, platinum, silver, iodine, as well as iron oxide), in varying NP concentrations (3−40 mg/mL), at two different photon beam energies (105 and 220 kVp). It was found that DER values at the skin lesion increased with the presence of bone when there was a higher atomic number of NPs, a higher NP concentration, and a lower photon beam energy. When comparing DER values with and without bone, using the same NP elements, NP concentration, and beam energy, differences were found in the range 0.04−3.55%, and a higher difference was found when the NP concentration increased. By considering the uncertainty in the DER calculation, the effect of bone scatter became significant to the dose enhancement (>2%) when the NP concentration was higher than 18 mg/mL. This resulted in an underestimation of dose enhancement at the skin lesion, when the bone underneath the tumour was neglected during orthovoltage radiotherapy.

Quantification of Nanoscale Dose Enhancement in Gold Nanoparticle-Aided External Photon Beam Radiotherapy

The recent progress in Nanotechnology has introduced Gold Nanoparticles (AuNPs) as promising radiosensitizing agents in radiation oncology. This work aims to estimate dose enhancement due to the presence of AuNPs inside an irradiated water region through Monte Carlo calculations. The GATE platform was used to simulate 6 MV photon histories generated from a TrueBeam^® linear accelerator with and without a Flattening Filter (FF) and model AuNPs clusters. The AuNPs size, concentration and distribution pattern were examined. To investigate different clinical irradiation conditions, the effect of field size, presence of FF and placement of AuNPs in water were evaluated. The range of Dose Enhancement Factors (DEF = Dose_Au/Dose_Water) calculated in this study is 0.99 ± 0.01-1.26 ± 0.02 depending on photon beam quality, distance from AuNPs surface, AuNPs size and concentration and pattern of distribution. The highest DEF is reported for irradiation using un-flattened photon beams and at close distances from AuNPs. The obtained findings suggest that dose deposition could be increased in regions that represent whole cells or subcellular targets (mitochondria, cell nucleus, etc.). Nevertheless, further and consistent research is needed in order to make a step toward AuNP-aided radiotherapy in clinical practice.

Monte Carlo simulations of nanodosimetry and radiolytic species production for monoenergetic proton and electron beams: Benchmarking of GEANT4-DNA and LPCHEM codes

PURPOSE

In hadrontherapy, biophysical models can be used to predict the biological effect received by cancerous tissues and organs at risk. The input data of these models generally consist of information on nano/micro dosimetric quantities and, concerning some models, reactive species produced in water radiolysis. In order to fully account for the radiation stochastic effects, these input data have to be provided by Monte Carlo track structure (MCTS) codes allowing to estimate physical, physico-chemical, and chemical effects of radiation at the molecular scale. The objective of this study is to benchmark two MCTS codes, Geant4-DNA and LPCHEM, that are useful codes for estimating the biological effects of ions during radiation therapy treatments.

MATERIAL AND METHODS

In this study we considered the simulation of specific energy spectra for monoenergetic proton beams (10 MeV) as well as radiolysis species production for both electron (1 MeV) and proton (10 MeV) beams with Geant4-DNA and LPCHEM codes. Options 2, 4, and 6 of the Geant4-DNA physics lists have been benchmarked against LPCHEM. We compared probability distributions of energy transfer points in cylindrical nanometric targets (10 nm) positioned in a liquid water box. Then, radiochemical species (· OH, e aq - ${\rm{e}}_{{\rm{aq}}}^ - $ , H 3 O + , H 2 O 2 ${{\rm{H}}_3}{{\rm{O}}^ + },{\rm{\;}}{{\rm{H}}_2}{{\rm{O}}_2}$ , H₂ , and O H - ) ${\rm{O}}{{\rm{H}}^ - }){\rm{\;}}$ yields simulated between 10^-12 and 10^-6 s after irradiation are compared.

RESULTS

Overall, the specific energy spectra and the chemical yields obtained by the two codes are in good agreement considering the uncertainties on experimental data used to calibrate the parameters of the MCTS codes. For 10 MeV proton beams, ionization and excitation processes are the major contributors to the specific energy deposition (larger than 90%) while attachment, solvation, and vibration processes are minor contributors. LPCHEM simulates tracks with slightly more concentrated energy depositions than Geant4-DNA which translates into slightly faster recombination than Geant4-DNA. Relative deviations (C_EV ) with respect to the average of evolution rates of the radical yields between 10^-12 and 10^-6 s remain below 10%. When comparing execution times between the codes, we showed that LPCHEM is faster than Geant4-DNA by a factor of about four for 1000 primary particles in all simulation stages (physical, physico-chemical, and chemical). In multi-thread mode (four threads), Geant4-DNA computing times are reduced but remain slower than LPCHEM by ∼20% up to ∼50%.

CONCLUSIONS

For the first time, the entire physical, physico-chemical, and chemical models of two track structure Monte Carlo codes have been benchmarked along with an extensive analysis on the effects on the water radiolysis simulation. This study opens up new perspectives in using specific energy distributions and radiolytic species yields from monoenergetic ions in biophysical models integrated to Monte Carlo software.

Comparing Geant4 physics models for proton-induced dose deposition and radiolysis enhancement from a gold nanoparticle

Gold nanoparticles (GNPs) are materials that make the tumor cells more radiosensitive when irradiated with ionizing radiation. The present study aimed to evaluate the impact of different physical interaction models on the dose calculations and radiochemical results around the GNP. By applying the Geant4 Monte Carlo (MC) toolkit, a single 50-nm GNP was simulated, which was immersed in a water phantom and irradiated with 5, 50, and 150 MeV proton beams. The present work assessed various parameters including the secondary electron spectra, secondary photon spectra, radial dose distribution (RDD), dose enhancement factor (DEF), and radiochemical yields around the GNP. The results with an acceptable statistical uncertainty of less than 1% indicated that low-energy electrons deriving from the ionization process formed a significant part of the total number of secondary particles generated in the presence of GNP; the Penelope model produced a larger number of these electrons by a factor of about 30%. Discrepancies of the secondary electron spectrum between Livermore and Penelope were more obvious at energies of less than 1 keV and reached the factor of about 30% at energies between 250 eV and 1 keV. The RDDs for Livermore and Penelope models were very similar with small variations within the first 6 nm from NP surface by a factor of 10%. In addition, neither the G-value nor the REF was affected by the choice of physical interaction models with the same energy cut-off. This work illustrated the similarity of the Livermore and Penelope models (within 15%) available in Geant4 for future simulation studies of GNP enhanced proton therapy with physical, physicochemical, and chemical mechanisms.

CPOP: An open source C++ cell POPulation modeler for radiation biology applications

PURPOSE

Multicellular tumor spheroids are realistic in-vitro systems used in radiation biology research to study the effect of anticancer drugs or to evaluate the resistance of cancer cells under specific conditions. When combining the modeling of spheroids together with the simulation of radiation using Monte Carlo methods, one could estimate cell and DNA damage to be compared with experimental data. We developed a Cell Population (CPOP) modeler combined to Geant4 simulations in order to tackle how energy depositions are allocated to cells, especially when enhancing radiation outcomes using high-Z nanoparticles. CPOP manages to model large three-dimensional cell populations with independent deformable cells described with their nucleus, cytoplasm and membranes together with force law systems to manage cell-cell interactions.

METHODS

CPOP is an opensource platform written in C++. It is divided into two main libraries: a "Modeler" library, for cell geometry modeling using meshes, and a Multi Agent System (MAS) library, simulating all agent (cell) interactions among the population. CPOP is fully interfaced with the Geant4 Monte Carlo toolkit and is able to directly launch Geant4 simulations after compilation. We modeled a full and realistic 3D cell population from SK-MEL28 melanoma cell population cultured experimentally. The spheroid diameter of 550 ± 40 µm corresponds to a population of approximately 1000 cells having a diameter of 17.2 ± 2.5 µm and a nucleus diameter of 11.2 ± 2.0 µm. We decided to reproduce cell irradiations performed with a X-RAD 320 Biological Irradiator (Precision XRay Inc., North Branford, CT).

RESULTS

We simulated the energy spectrum of secondary particles generated in the vicinity of the spheroid and plotted the different energy spectra recovered internally to the spheroid. We evaluated also the impact of AGuIX (Gadolinium) nanoparticles modeled into the spheroid with their corresponding secondary energy spectra.

CONCLUSIONS

We succeeded into modeling cell populations and combined them with Geant4 simulations. The next step will be to integrate DNA geometrical models into cell nuclei and to use the Geant4-DNA physics and radiolysis modeling capabilities in order to evaluate early strand breaks induced on DNA.

Impact of Hypoxia on Relative Biological Effectiveness and Oxygen Enhancement Ratio for a 62-MeV Therapeutic Proton Beam

This study uses the yields of double-strand breaks (DSBs) to determine the relative biological effectiveness (RBE) of proton beams, using cell survival as a biological endpoint. DSB induction is determined when cells locate at different depths (6 positions) along the track of 62 MeV proton beams. The DNA damage yields are estimated using Monte Carlo Damage Simulation (MCDS) software. The repair outcomes are estimated using Monte Carlo excision repair (MCER) simulations. The RBE for cell survival at different oxygen concentrations is calculated using the repair-misrepair-fixation (RMF) model. Using 60Co γ-rays (linear energy transfer (LET) = 2.4 keV/μm) as the reference radiation, the RBE for DSB induction and enzymatic DSB under aerobic condition (21% O2) are in the range 1.0-1.5 and 1.0-1.6 along the track depth, respectively. In accord with RBE obtained from experimental data, RMF model-derived RBE values for cell survival are in the range of 1.0-3.0. The oxygen enhancement ratio (OER) for cell survival (10%) decreases from 3.0 to 2.5 as LET increases from 1.1 to 22.6 keV/μm. The RBE values for severe hypoxia (0.1% O2) are in the range of 1.1-4.4 as LET increases, indicating greater contributions of direct effects for protons. Compared with photon therapy, the overall effect of 62 MeV proton beams results in greater cell death and is further intensified under hypoxic conditions.

Monte Carlo-based calculation of nano-scale dose enhancement factor and relative biological effectiveness in using different nanoparticles as a radiosensitizer

INTRODUCTION

Nowadays, some nanoparticles (NPs) are known and used as radiosensitizers in radiotherapy and radiobiology, due to their desired biological, physical, and chemical effects on cells. This study aimed to evaluate and compare the dose enhancement factor (DEF) and the biological effectiveness of some common NPs through EGSnrc and MCDS Monte Carlo (MC) simulation codes.

MATERIALS AND METHODS

To evaluate considered NPs' DEF, a single NP with 50 nm diameter was simulated at the center of concentric spheres. NP irradiations were done with 30, 60, and 100 keV photon energies. The secondary electron spectra were scored at the surface of considered NPs, and the dose values were scored at surrounding water-filled spherical shells which were distributed up to 4000 nm from the NP surface. The electron spectra were used in the MCDS code to obtain different initial DNA damages for the calculation of enhanced relative biological effectiveness (eRBE).

RESULTS

By decreasing the photon energy, an increment of DEF was seen for all studied NPs. The maximum DEF at 30, 60, and 100 keV photon energies were respectively related to silver (Ag), gadolinium (Gd), and bismuth (Bi) NPs. The maximum double-strand break (DSB) related (eRBEDSB) values for the 30 keV photon belonged to Ag, while BiNPs showed the maximum values at other photon energies. The minimum eRBEDSB values were also related to iron (Fe) NPs at the entire range of studied photon energies.

CONCLUSIONS

The compared nanoscale physical and biological results of our study can be helpful in the selection of optimum NP as a radiosensitizer in future radiobiological studies. Bi, gold (Au), Ag, and platinum (Pt) NPs had great potential, respectively, as radiosensitizers relative to the other studied NPs.

IDDRRA: A novel platform, based on Geant4-DNA to quantify DNA damage by ionizing radiation

PURPOSE

This study proposes a novel computational platform that we refer to as IDDRRA (DNA Damage Response to Ionizing RAdiation), which uses Monte Carlo (MC) simulations to score radiation induced DNA damage. MC simulations provide results of high accuracy on the interaction of radiation with matter while scoring the energy deposition based on state-of-the-art physics and chemistry models and probabilistic methods.

METHODS

The IDDRRA software is based on the Geant4-DNA toolkit together with new tools that were developed for the purpose of this study, including a new algorithm that was developed in Python for the design of the DNA molecules. New classes were developed in C++ to integrate the GUI and produce the simulation's output in text format. An algorithm was also developed to analyze the simulation's output in terms of energy deposition, Single Strand Breaks (SSB), Double Strand Breaks (DSB) and Cluster Damage Sites (CDS). Finally, a new tool was developed to implement probabilistic SSB and DSB repair models using MC techniques.

RESULTS

This article provides the first benchmarks that the user of the IDDRRA tool can use to validate the functionality of the software as well as to provide a starting point to produce different types of DNA simulations. These benchmarks incorporate different kind of particles (e-, e+, protons, electron spectrum) and DNA molecules.

CONCLUSION

We have developed the IDDRRA tool and demonstrated its use to study various aspects of the modeling and simulation of a DNA irradiation experiment. The tool is expandable and can be expanded by other users with new benchmarks and applications based on the user's needs and experience. New functionality will be added over time, including the quantification of the indirect damage.

Combined cell and nanoparticle models for TOPAS to study radiation dose enhancement in cell organelles

Dose enhancement by gold nanoparticles (AuNP) increases the biological effectiveness of radiation damage in biomolecules and tissue. To apply them effectively during cancer therapy their influence on the locally delivered dose has to be determined. Hereby, the AuNP locations strongly influence the energy deposit in the nucleus, mitochondria, membrane and the cytosol of the targeted cells. To estimate these effects, particle scattering simulations are applied. In general, different approaches for modeling the AuNP and their distribution within the cell are possible. In this work, two newly developed continuous and discrete-geometric models for simulations of AuNP in cells are presented. These models are applicable to simulations of internal emitters and external radiation sources. Most of the current studies on AuNP focus on external beam therapy. In contrast, we apply the presented models in Monte-Carlo particle scattering simulations to characterize the energy deposit in cell organelles by radioactive ¹⁹⁸AuNP. They emit beta and gamma rays and are therefore considered for applications with solid tumors. Differences in local dose enhancement between randomly distributed and nucleus targeted nanoparticles are compared. Hereby nucleus targeted nanoparticels showed a strong local dose enhancement in the radio sensitive nucleus. These results are the foundation for future experimental work which aims to obtain a mechanistic understanding of cell death induced by radioactive ¹⁹⁸Au.

Intercomparison of Monte Carlo calculated dose enhancement ratios for gold nanoparticles irradiated by X-rays: Assessing the uncertainty and correct methodology for extended beams

Results of a Monte Carlo code intercomparison exercise for simulations of the dose enhancement from a gold nanoparticle (GNP) irradiated by X-rays have been recently reported. To highlight potential differences between codes, the dose enhancement ratios (DERs) were shown for the narrow-beam geometry used in the simulations, which leads to values significantly higher than unity over distances in the order of several tens of micrometers from the GNP surface. As it has come to our attention that the figures in our paper have given rise to misinterpretation as showing 'the' DERs of GNPs under diagnostic X-ray irradiation, this article presents estimates of the DERs that would have been obtained with realistic radiation field extensions and presence of secondary particle equilibrium (SPE). These DER values are much smaller than those for a narrow-beam irradiation shown in our paper, and significant dose enhancement is only found within a few hundred nanometers around the GNP. The approach used to obtain these estimates required the development of a methodology to identify and, where possible, correct results from simulations whose implementation deviated from the initial exercise definition. Based on this methodology, literature on Monte Carlo simulated DERs has been critically assessed.

Influence of gold nanoparticles embedded in water on nanodosimetry for keV photon irradiation

PURPOSE

For the past two decades, high-Z nanoparticles have been of high interest to improve the therapeutic outcomes of radiation therapy, especially for low-energy x-rays. Monte Carlo (MC) simulations have been used to evaluate the boost of dose deposition induced by Auger electrons near the nanoparticle surface, by calculating average energy deposition at the nanoscale. In this study, we propose to go beyond average quantities and quantify the stochastic nature of energy deposition at such a scale. We present results of probability density of the specific energy (restricted to ionization, excitation and electron attachment events) in cylindrical nanotargets of height and radius set at 10 nm. This quantity was evaluated for nanotargets located within 200 nm around 5-50 nm gold nanoparticles (GNPs), for 20-90 keV photon irradiation.

METHODS

This nanodosimetry study was based on the MC simulation MDM that allows tracking of electrons down to thermalization energy. We introduced a new quantity, namely the probability enhancement ratio (PER), by estimating the probability of imparting to nanotargets a restricted specific energy larger than a threshold z 0 (1, 10, and 20 kGy), normalized to the probability for pure water. The PER was calculated as a function of the distance between the nanotarget and the GNP surface. The threshold values were chosen in light of the biophysical model NanOx that predicts cell survival by calculating local lethal events based on the restricted specific energy and an effective local lethal function. z 0 then represents a threshold above which the nanotarget damages induce efficiently cell death.

RESULTS

Our calculations showed that the PER varied a lot with the GNP radius, the photon energy, z 0 and the distance of the GNP to the nanotarget. The highest PER was 95 when the nanotarget was located at 5 nm from the GNP surface, for a photon energy of 20 keV, a threshold of 20 kGy, and a GNP radius of 50 nm. This enhancement dramatically decreased with increasing GNP-nanotarget distances as it went below 1.5 for distances larger than 200 nm.

CONCLUSIONS

The PER seems better adapted than the mean dose deposition to describe the formation of biological damages. The significant increase of the PER within 200 nm around the GNP suggests that severe damages could occur for biological nanotargets located near the GNP. These calculations will be used as an input of the biophysical model NanOx to convert PER into estimation of radiation-induced cell death enhanced by GNPs.

A detailed experimental and Monte Carlo analysis of gold nanoparticle dose enhancement using 6 MV and 18 MV external beam energies in a macroscopic scale

Dose enhancement due to gold nanoparticles (GNPs) has been quantified experimentally and through Monte Carlo simulations for external beam radiation therapy energies of 6 and 18 MV. The highest enhancement was observed for the 18 MV beam at the highest GNP concentration tested, amounting to a DEF of 1.02. DEF is shown to increase with increasing concentration of gold and increasing energy in the megavoltage energy range. The largest difference in measured vs. simulated DEF across all data sets was 0.3%, showing good agreement.

An in silico study on the effect of host tissue at brachytherapy dose enhancement by gold nanoparticles

PURPOSE

Iridium-192 brachytherapy dose enhancement by gold nanoparticles was investigated in five different tumor tissues to observe the tissue-related differences as an effective environmental factor in the applications of nanoparticles as radio-enhancer agents.

METHODS AND MATERIALS

The brachytherapy high-dose-rate source of BEBIG Ir-192, a tumor volume with five different tissues including water, Plexiglas, soft tissue, adipose, and bone with and without a uniform distribution of gold nanoparticles were mimicked by MCNPX Monte Carlo simulation code using lattice feature. Dose enhancement factors in the tumor volume were measured separately regarding the types of tissue, and a previous study using GEometry ANd Tracking 4 simulation was used for result validation.

RESULTS

The results demonstrated that various types of tissue, as the host of gold nanoparticles, lead to different dose enhancement level, so that the bone and adipose have the lowest and the highest amount of dose enhancement factor with values 20.8% and 39.75%, respectively. The maximum difference of 4.8% was achieved from data benchmarking.

CONCLUSIONS

The results of this study indicate that the MCNPX code can be used as a valid tool for dose measurement in the presence of nanoparticles. Moreover, tissue types of tumor as an environmental feature, alongside with the nanoparticle's size and concentration as well as the conditions of radiotherapy, should be considered in the dose calculation.

Advances in modelling gold nanoparticle radiosensitization using new Geant4-DNA physics models

Gold nanoparticles have demonstrated significant radiosensitization of cancer treatment with x-ray radiotherapy. To understand the mechanisms at the basis of nanoparticle radiosensitization, Monte Carlo simulations are used to investigate the dose enhancement, given a certain nanoparticle concentration and distribution in the biological medium. Earlier studies have ordinarily used condensed history physics models to predict nanoscale dose enhancement with nanoparticles. This study uses Geant4-DNA complemented with novel track structure physics models to accurately describe electron interactions in gold and to calculate the dose surrounding gold nanoparticle structures at nanoscale level. The computed dose in silico due to a clinical kilovoltage beam and the presence of gold nanoparticles was related to in vitro brain cancer cell survival using the local effect model. The comparison of the simulation results with radiobiological experimental measurements shows that Geant4-DNA and local effect model can be used to predict cell survival in silico in the case of x-ray kilovoltage beams.

Advancements in the use of Auger electrons in science and medicine during the period 2015-2019

Auger electrons can be highly radiotoxic when they are used to irradiate specific molecular sites. This has spurred basic science investigations of their radiobiological effects and clinical investigations of their potential for therapy. Focused symposia on the biophysical aspects of Auger processes have been held quadrennially. This 9th International Symposium on Physical, Molecular, Cellular, and Medical Aspects of Auger Processes at Oxford University brought together scientists from many different fields to review past findings, discuss the latest studies, and plot the future work to be done. This review article examines the research in this field that was published during the years 2015-2019 which corresponds to the period since the last meeting in Japan. In addition, this article points to future work yet to be done. There have been a plethora of advancements in our understanding of Auger processes. These advancements range from basic atomic and molecular physics to new ways to implement Auger electron emitters in radiopharmaceutical therapy. The highly localized doses of radiation that are deposited within a 10 nm of the decay site make them precision tools for discovery across the physical, chemical, biological, and medical sciences.

Ionizing Radiation and Complex DNA Damage: Quantifying the Radiobiological Damage Using Monte Carlo Simulations

Ionizing radiation is a common tool in medical procedures. Monte Carlo (MC) techniques are widely used when dosimetry is the matter of investigation. The scientific community has invested, over the last 20 years, a lot of effort into improving the knowledge of radiation biology. The present article aims to summarize the understanding of the field of DNA damage response (DDR) to ionizing radiation by providing an overview on MC simulation studies that try to explain several aspects of radiation biology. The need for accurate techniques for the quantification of DNA damage is crucial, as it becomes a clinical need to evaluate the outcome of various applications including both low- and high-energy radiation medical procedures. Understanding DNA repair processes would improve radiation therapy procedures. Monte Carlo simulations are a promising tool in radiobiology studies, as there are clear prospects for more advanced tools that could be used in multidisciplinary studies, in the fields of physics, medicine, biology and chemistry. Still, lot of effort is needed to evolve MC simulation tools and apply them in multiscale studies starting from small DNA segments and reaching a population of cells.

Radio-enhancement effects by radiolabeled nanoparticles

In cancer radiation therapy, dose enhancement by nanoparticles has to date been investigated only for external beam radiotherapy (EBRT). Here, we report on an in silico study of nanoparticle-enhanced radiation damage in the context of internal radionuclide therapy. We demonstrate the proof-of-principle that clinically relevant radiotherapeutic isotopes (i.e. 213Bi, 223Ra, 90Y, 177Lu, 67Cu, 64Cu and 89Zr) labeled to clinically relevant superparamagnetic iron oxide nanoparticles results in enhanced radiation damage effects localized to sub-micron scales. We find that radiation dose can be enhanced by up to 20%, vastly outperforming nanoparticle dose enhancement in conventional EBRT. Our results demonstrate that in addition to the favorable spectral characteristics of the isotopes and their proximity to the nanoparticles, clustering of the nanoparticles results in a nonlinear collective effect that amplifies nanoscale radiation damage effects by electron-mediated inter-nanoparticle interactions. In this way, optimal radio-enhancement is achieved when the inter-nanoparticle distance is less than the mean range of the secondary electrons. For the radioisotopes studied here, this corresponds to inter-nanoparticle distances <50 nm, with the strongest effects within 20 nm. The results of this study suggest that radiolabeled nanoparticles offer a novel and potentially highly effective platform for developing next-generation theranostic strategies for cancer medicine.