Multiscale Monte Carlo simulations of gold nanoparticle dose-enhanced radiotherapy I. Cellular dose enhancement in microscopic models

Martin P. Martinov,Elizabeth M. Fletcher,Rowan M. Thomson
2023-05-03
Abstract:Purpose: The introduction of Gold NanoParticles (GNPs) in radiotherapy requires consideration of GNP size, location, and quantity, and more, as well as beam quality. The work is presented in two parts, with Part I (this work) investigating accurate and efficient MC modeling to calculate nucleus and cytoplasm Dose Enhancement Factors (n,cDEFs). Part II then evaluates cell dose enhancement factors across macroscopic length scales. Methods: Different methods of modeling gold within cells are compared: a contiguous volume of either pure gold or gold-tissue mixture or discrete GNPs. MC simulations with EGSnrc are performed to calculate n,cDEF for a cell considering 10 to 370 keV incident photons, gold concentrations from 4 to 24~mg of gold per gram, and three different GNP configurations: GNPs distributed around the surface of the nucleus or packed into one or four endosomes. Results: n,cDEFs are sensitive to the model of gold in the cell, with differences of up to 17% observed between models. Across cell/nucleus radii, source energies, and gold concentrations, n,cDEF are highest for GNPs in the perinuclear configuration, compared with GNPs in endosomes. Across all simulations of the default cell, nDEFs and cDEFs range from one to 6.83 and 3.87, respectively. If the cell size is changed, much higher n,cDEFs are observed. n,cDEF are maximized at photon energies at or just above the K- or L-edges of gold. Conclusions: This work demonstrates many physics trends on DEFs at the cellular level, including that cellular DEFs are sensitive to gold modeling approach, GNP configuration, cell/nucleus size, gold concentration, and incident source energy. These data will allow one to optimize or estimate DEF using not only GNP uptake, but account for incident photon energy and intracellular configuration of GNPs. Part II will expand the investigation, applying the cell model in cm-scale phantoms.
Medical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to accurately evaluate the dose - enhancement effect (DEF) when using gold nanoparticles (GNPs) as sensitizers in radiotherapy. Specifically, the paper focuses on studying the dose - enhancement effect of gold nanoparticles at the cellular level through multi - scale Monte Carlo (MC) simulations, and how this effect is influenced by factors such as the concentration of gold nanoparticles, intracellular distribution, cell size, and incident photon energy. In addition, the paper also explores the impact of different modeling methods on the results, aiming to find an accurate and efficient simulation method to guide treatment planning in clinical applications. ### Background and Purpose of the Paper **Background**: - The application of gold nanoparticles (GNPs) in radiotherapy needs to consider multiple factors, such as the size, position, and quantity of GNPs, as well as the patient's geometric structure and beam quality. - Physical considerations span multiple length scales (from nanometers to centimeters), which pose challenges to dosimetry research and usually limit the research scope to only the micro - or macro - scale. **Purpose**: - Through Monte Carlo (MC) simulations, study the effect of gold nanoparticle - enhanced radiotherapy (GNPT) between the micro - and macro - scales. - The first part (this paper) focuses on accurate and efficient MC modeling at the single - cell level, calculating the nuclear and cytoplasmic dose - enhancement factors (nDEF and cDEF), considering a wide range of parameter spaces, including the concentration of GNPs, intracellular distribution, cell size, and incident photon energy. - The second part evaluates the cell - dose - enhancement factor at the macro (tumor) scale. ### Methods - **Different gold - modeling methods**: Compared three methods: continuous pure - gold volume, gold - tissue mixture, and hexagonal close - packed lattice. - **MC simulations**: Used EGSnrc to conduct MC simulations, calculating nDEF and cDEF for cells with a radius of 7.35 µm and nuclei with a radius of 5 µm under incident photons from 10 to 370 keV, considering gold concentrations from 4 to 24 mg Au/g tissue, and three different GNP configurations: perinuclear, 1 - endo, and 4 - endo. - **Extended simulations**: Selected some cells with different cell sizes (5 µm, 7.35 µm, and 10 µm) for simulation to further explore the impact of cell size on the dose - enhancement factor. ### Results - **Sensitivity of modeling methods**: Differences in DEF caused by different gold - modeling methods can reach 17%, among which the hexagonal close - packed lattice is considered the best model. - **Impact of intracellular GNP configuration**: For the perinuclear configuration, nDEF and cDEF are the highest, followed by the 4 - endo configuration, and the lowest is the 1 - endo configuration. - **Impact of cell size**: The sizes of cells and nuclei have a significant impact on DEF, especially for high gold concentrations and specific photon energies (10 - 20 keV). - **Impact of photon energy**: DEF values are the highest at photon energies (10 - 20 keV) near the K - edge and L - edge of gold. ### Conclusions - Through 5,000 unique simulation scenarios, comprehensively demonstrated multiple physical trends of DEF at the cellular level, including the impact of modeling methods, intracellular GNP configuration, cell / nuclear size, gold concentration, and incident source energy on DEF. - These data are very useful for research and treatment planning, and can optimize or estimate DEF, considering not only GNP uptake, but also the average tumor cell size, incident photon energy, and intracellular GNP configuration. - The second part will further expand the research by applying the cell model in the first part to centimeter - scale simulations. ### Keywords - Gold nanoparticles - Monte Carlo - Microdosimetry