Corresponding author: Ivan A. Konobeev ( beo0@mail.ru ) Academic editor: Yury Korovin
© 2019 Ivan A. Konobeev, Yurij A. Kurachenko, Igor’ N. Sheino.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Konobeev IA, Kurachenko YA, Sheino IN (2019) Impact of secondary particles on microdistribution of deposited dose in biological tissue in the presence of gold and gadolinium nanoparticles under photon beam irradiation. Nuclear Energy and Technology 5(2): 109-116. https://doi.org/10.3897/nucet.5.35798
|
It is experimentally proven that nanoparticles of high-Z materials can be used as radiosensitizers for photon beam therapy. In the authors’ opinion, data available as of today on the impact of secondary particles (electrons, photons and positrons generated in biological tissue by penetrating beam of primary photons) on the distribution of deposited dose during photon beam therapy in the presence of nanoparticles, are insufficient. Investigation of this impact constituted the main goal of this work.
Two-stage simulation was performed using Geant4 platform. During the first stage a layer of biological tissue (water) was irradiated by monoenergetic photon sources with energies ranging from 10 keV to 6 MeV. As the result of this modeling spectra of electrons, photons and positrons were obtained at the depth of 5 cm. During the second stage the obtained photon spectra were used to irradiate gold, gadolinium and water nanoparticles. Radial distributions of energy deposited around nanoparticles were obtained as the result of this modeling.
Radial DEF (Dose Enhancement Factor) values around nanoparticles of gold and gadolinium positioned in water at the depth of 5 cm were obtained after processing the collected data. Contributions from primary photons and secondary particles (electrons, photons and positrons generated in the layer of water with 5-cm thickness by the penetrating beam of primary photons) in the additional dose deposited around the nanoparticles were calculated as well.
It was demonstrated that layer of biological tissue placed between the source of photons and nanoparticles considerably changes the initial spectrum of photons and this change is significant in the analysis of mechanism of radiosensitization of biological tissues by nanoparticles for all energies of photon sources (up to 6 MeV).
It was established that interaction of electrons and positrons with nanoparticles does not lead to significant increase of additional dose in the vicinity of their surfaces and can be most likely excluded from consideration in the analysis of radiosensitization mechanism of nanoparticles.
Nanoparticles of gold and gadolinium, deposited dose, photon beam therapy, Geant4, Monte Carlo simulation
Methods for physical aiming of radiation on the tumor, such as three-dimensional conformal radiation therapy (3D CRT), intensity-modulated radiotherapy (IMRT), image guided radiation therapy (IGRT), stereotactic radiation surgery (SRT) ensuring maximum precision of radiation dose rendering to tumor target have reached the limits of their further refinement (
One of the options of radiosensitization is the introduction in the biological medium of elements with significantly higher radiation absorption cross-section than that of the biological tissue per se. Emerging secondary short-range radiation localizes energy absorption in the vicinity of these elements and affects only adjacent biological structures.
Main principles of binary radiation therapy technologies such as the neutron capture therapy (NCT) or photon capture therapy (PCT) refer to the methods of radiation sensitization. In PCT technology additional energy yield under irradiation of biological tissue with photons is caused by the emission of photoelectrons and accompanying Auger-cascade from atoms of “heavy” elements (high-Z elements): 53I, 64Gd, 78Pt, 79Au and other elements included in the compositions of a number of preparations (
The most widely spread approach to the explanation of radiation sensitization of biological tissue under photon radiation is the local increase of deposited dose in the biological medium due to the presence in the medium of elements with high atomic number (
In vitro and in vivo experiments demonstrated significant effects of radiation sensitization by nanoparticles of gold for different cellular lines with application of X-ray beams within kilovolt range (
The reason of discrepancy between calculated values of increase of deposited dose and observable biological effect of necrocytosis is the significant heterogeneity of the dose field emerging in the case of manifested localization of energy absorption in the direct vicinity from “heavy” elements which may reach extremely high values (
Such significant heterogeneities of dose can become the cause of a whole series of biochemical processes, including high generation of deleterious and hydroxyl radicals and other ensuing processes of complex biological damage (
The heterogeneity in question appears under the effects on heavy elements of both photons emitted by the source, and of secondary particles generated in the biological medium (photons, electrons and positrons).
As of today, available data on the effects of secondary particles on dose heterogeneity in the vicinity of heavy elements are insufficient. Effects of electrons are often disregarded without indicating the uncertainty introduced by such approximation.
The objective of the present study is the investigation of distributions of deposited dose in biological tissue in the presence of nanoparticles of gold and gadolinium depending on the photon source energy and effects of secondary particles generated in the biological medium on these distributions.
All photons are divided in the paper into two groups: “primary” and “secondary” photons. Photons are considered to be “primary” if their energy is equal to the energy of photons emitted from monoenergetic source. In particular, photons undergoing Raleigh scattering remain to be “primary”. Photons are considered to be “secondary” if their energy is less than the energy of photons emitted from the source. In particular, photons that underwent Compton scattering or photons emitted as the result of annihilation are attributed to the group of “secondary” photons.
One nanoparticle positioned at the depth of 5 cm from flat monoenergetic photon source was simulated for investigating distribution of additional dose associated with presence of nanoparticles in the water (as a tissue-equivalent material). The geometry in question is a simplified model of small-size tumor located at the depth of 5 cm in human body and containing nanoparticles. The problem had to be solved in two stages because of the disproportion between the dimensions of the photon beam and the nanoparticle size. Each stage represented separate calculation using Geant4 platform of version 10.1 (
The main objective of the first stage was obtaining spectra of electrons, photons and positrons generated in water after irradiation by monoenergetic photon source with energy within the range from 10°keV to 6°MeV.
The following problem was simulated. Flat monoenergetic homogenous photon source with 10-cm diameter is located in the center of water cylinder base with 20-cm height and 20-cm diameter. Photons escape the source at right angle to its surface. The cylinder contains detector positioned at the distance of 5 cm from its base and registering the spectra of particles passing through it. The detector is a cylinder of water with thickness equal to 1 μm and diameter equal to 20 cm (Fig.
Livermore (in a larger cylinder) and DNA (in the detector) physical models (
Spectra of electrons, photons and positrons were stored with pitch equal to 1°keV. Therefore, spectra obtained for source with maximum energy of 6 MeV contained 6000 points each. Spectrum of electrons with energies up to 1 keV was stored separately with pitch equal to 1 eV and contained 1000 points. Spectra of only those particles which were flying at angles [0;π/2) against the direction of flight of photons escaping the source were stored during simulation.
The main objective of this stage was the determination of radial distributions of absorbed energy around nanoparticles irradiated with equivalent source of photons, electrons and positrons with spectra obtained during the first stage of simulation.
The following problem was simulated. Single nanoparticle (homogenous spherical ball with diameter equal to 2 nm and 30 nm consisting of gold, gadolinium and water) was irradiated by circular source with diameter equal to the diameter of the nanoparticle positioned directly adjacent to it (see Fig.
Fluences of particles per primary photon generated at the depth of 5 cm from the monoenergetic photon source inside and outside of the original beam: 1 – electrons inside the beam; 2 – electrons outside the beam; 3 – primary photons inside the beam; 4 – primary photons outside the beam; 5 – secondary photons inside the beam; 6 – secondary photons outside the beam; 7 – positrons inside the beam; 8 – positrons outside the beam. Thin lines indicate the corridors of statistical uncertainties of calculations; upper straight line corresponds to the initial fluence of the source (photon fluence at the depth of 0 cm).
Energy absorbed in spherical layers with 1 nm thickness surrounding nanoparticle up to the distance of 10 μm from the center of the nanoparticle was stored during simulation. Absorbed energy was stored separately for the following four types of particles emitted from the source: electrons, “primary” photons (photons with energy equal to the energy of primary photons during the first stage of simulation), “secondary” photons (photons with energies smaller than the energy of primary photons during the first stage of simulation) and positrons.
Livermore (inside nanoparticle) and DNA (inside the spherical water ball) physical models were used for describing propagation of electrons. Threshold of secondary particle generation in the nanoparticle was equal to 10 eV. Propagation of electrons in the water ball was simulated down to the energy equal to 5.1 eV. Propagation of ions and radicals was not taken into account.
Radial DEF (Dose Enhancement Factor, i.e. the ratio of dose in the presence of nanoparticle of gold or gadolinium to the dose value in the absence of nanoparticle) values and contributions from electrons, primary photons, secondary photons and positrons in the additional deposited dose within the area from the surface to 100 nm from the center of the nanoparticle (area of peak radial dose, see Fig.
DEFXNP = 1 + (DXNP – DWNP)/D*,
where DXNP is the calculated average dose in the spherical layer surrounding the nanoparticle consisting of material X (gold or gadolinium); DWNP is the calculated average dose in the same spherical layer surrounding the nanoparticle consisting of water and irradiated with the same particle beam; D* is the average dose in the water in the absence of nanoparticles consisting of material X calculated as the average calculated dose inside nanoparticle consisting of water. Each dose value was calculated as the ratio of energy absorbed inside the volume to its mass.
Thus obtained DEF values are close to those which are expected to be observed in experiment.
Contribution of each type of particles in the additional dose was calculated as the difference between deposited doses from the given type of particles when they irradiate nanoparticle made of agent material and nanoparticle consisting of water.
Dependences of fluences of electrons, photons and positrons emitted at the depth of 5 cm on the energy of primary photons for monoenergetic source are presented in Figure
At the energy of the source equal to 10 keV fluences for all particles at the depth of 5 cm are equal to zero. Corridors of uncertainties for photons inside the beam are very narrow and, therefore, they are indistinguishable in the plots.
Spectra of photons and electrons emitted in the water at the depth of 5 cm from the monoenergetic photon sources with energies equal to 300 keV, 1000 keV and 6000 keV are presented in Figure
As it has been described above in the section “Materials and Methods” spectra were stored with pitch along the energy axis equal to 1 keV. Therefore, fluence of particles with energy [E; E + 1) keV, where E is the positive integer value corresponds to the point with energy value equal to (E + 0.5)°keV.
A large number of electrons “inside the beam” with energies within 1 keV (about 4. 1 and 0.6% of the total number of electrons “inside the beam” for 300, 1000 and 6000 keV, respectively) are also present in the spectra. Spectra in the region have pronounced peak within the range of 8 – 20 eV. In the process of irradiation of nanoparticle these electrons reduce additional deposited dose associated with presence of the nanoparticle by being efficiently absorbed inside the nanoparticle.
As it is clear from Fig. 4а, low energy component appears in the photon spectrum during penetration of monoenergetic photon beam through the layer of water. At the energies of primary photons exceeding two electron rest energies (1022 keV) large number of photons with energy of 511 keV (annihilation peak) are formed in the spectrum.
Radial distributions of DEF in the layers surrounding nanoparticles positioned in the water at the depth of 5 cm and irradiated with monoenergetic photon sources (from the first stage) are presented in Figure
Similar to the cross-section of photon interaction with agent material, additional dose decreases with increased photon energy and has a jump corresponding to the penetration through K-shell of the agent material (energy of K-shell amounts to approximately 80.7 keV for gold and to 50.2°keV for gadolinium). Within the range from 600 keV to 6000°keV “inside-the-beam” DEF curves are very close to each other and are not discernible within the limits of calculation uncertainty.
It has to be noted that while radial DEF in Fig.
Radial DEF for nanoparticles: a) – Au ϕ 30 nm, inside the beam; b) – Gd ϕ 30 nm, inside the beam; c) – Au ϕ 30 nm, outside the beam; d) – Au ϕ 2 nm, inside the beam. Figures near the curves indicate energies of primary photons in keV. Curves for 81 keV and 100°keV in Fig. c) practically coincide.
Histograms of contributions from each type of particles escaping from the equivalent source (electrons, primary photons, secondary photons and positrons) in the additional deposited dose within the area from the surface to 100 nm from center of the nanoparticles with 30-nm diameter positioned inside the beam, are presented in Figure
Thus, photon spectrum changes significantly when photon beam penetrates through the water (Fig. 4а). Photons with lower energies which ionize nanoparticles more intensively appear in the spectrum. Contribution of secondary photons in the additional dose in the vicinity of nanoparticle surface appears to be significant for all examined energies of the primary photon source (see Fig.
Contributions in the additional dose within the area from the surface to 100 nm from the center of nanoparticles of gold and gadolinium with ϕ 30 nm positioned inside the beam from electrons, primary photons, secondary photons and positrons: 1 – electrons; 2 – primary photons; 3 – secondary photons; 4 – positrons. Statistical uncertainties are indicated with bars.
Conclusion can be drawn by comparing Figures
Let us note that while water layer increases efficiency of the photon spectrum per se in its presence DEF around nanoparticles will be significantly lower than in its absence because of considerably higher flux of electrons which cannot efficiently ionize nanoparticle and, therefore, they will be leveling the ratio of values of doses in the presence of nanoparticle and in its.
Sharp jumps of deposited dose in the vicinity of surfaces of nanoparticles emerge under irradiation of biological tissue with photons with energies up to 6 MeV. This is explained, in the first place, by photoelectric absorption of photons in nanoparticles with subsequent emission of short-range electrons.
Irradiation with photons of tissue equivalent phantom results in the emergence of low-energy components in the photon spectrum inside the tissue and in the appearance of fluxes of electrons and positrons. It was established that such variation of photon spectrum is significant in the analysis of the mechanism of radiation sensitization of biological tissue by nanoparticles for any energies of the photon source (up to 6°MeV).
It was demonstrated that interaction of electrons and positrons with nanoparticles does not result in significant growth of additional dose in the vicinity of surfaces of nanoparticles and, most likely, can be excluded from further examination in the analysis of the mechanism of radiation sensitization of biological tissue by nanoparticles.
For quantitative description of radiation sensitization of tissues by nanoparticles it is necessary to take into account a number of biochemical processes, including high rate of generation of deleterious hydroxyl radicals and other subsequent processes of complex biological damages.