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

Abstract

Background: The introduction of Gold NanoParticles (GNPs) in radiotherapy treatments necessitates considerations such as GNP size, location, and quantity, as well as patient geometry and beam quality. Physics considerations span length scales across many orders of magnitude (nanometer-to-centimeter), presenting challenges that often limit the scope of dosimetric studies to either micro-or macroscopic scales. Purpose: To investigate GNP dose-enhanced radiation Therapy (GNPT) through Monte Carlo (MC) simulations that bridge micro-to-macroscopic scales. The work is presented in two parts, with Part I (this work) investigating accurate and efficient MC modeling at the single cell level to calculate nucleus and cytoplasm Dose Enhancement Factors (n,cDEFs), considering a broad parameter space including GNP concentration, GNP intracellular distribution, cell size, and incident photon energy. Part II then evaluates cell dose enhancement factors across macroscopic (tumor) length scales. Methods: Different methods of modeling gold within cells are compared, from a contiguous volume of either pure gold or gold-tissue mixture to discrete GNPs in a hexagonal close-packed lattice. MC simulations with EGSnrc are performed to calculate n,cDEF for a cell with radius r cell = 7.35 µm and nucleus r nuc = 5 µm considering 10 to 370 keV incident photons, gold concentrations from 4 to 24 mg Au /g tissue , and three different GNP configurations within the cell: GNPs distributed around the surface of the nucleus (perinuclear) or GNPs packed into one (or four) endosome(s). Select simulations are extended to cells with different cell (and nucleus) sizes: 5 µm (2, 3, and 4 µm), 7.35 µm (4 and 6 µm), and 10 µm (7, 8, and 9 µm). Results: n,cDEFs are sensitive to the method of modeling gold in the cell, with differences of up to 17% observed; the hexagonal lattice of GNPs is chosen (as the most realistic model) for all subsequent simulations. Across cell/nucleus radii, source energies, and gold concentrations, both nDEF and cDEF are highest for GNPs in the perinuclear configuration, compared with GNPs in one (or four) endosome(s). Across all simulations of the (r cell , r nuc ) = (7.35, 5) µm cell, nDEFs and cDEFs range from unity to 6.83 and 3.87, respectively. Including different cell sizes, nDEFs and cDEFs as