Skip to main content
EJNMMI Physics
Download PDF
  • Original research
  • Open access
  • Published:

Distribution modeling of nanoparticles for brachytherapy of human eye tumor

EJNMMI Physics volume 7, Article number: 53 (2020) Cite this article

Abstract

Background

Due to their unique properties, gold nanoparticles (GNPs) have been proposed to be used for a wide range of applications, especially for photon radiation therapy. In addition to experimental works, there are worthwhile simulation-based studies focused on the investigation of the effect of parameters governing the dose enhancement due to the presence of GNPs in tissue. In a recently published study, we found that the distribution of GNPs in a single cell plays an important role in nucleus dose enhancement.

Methods

The present work investigates the sensitivity of dose enhancement of a macroscopic phantom to the modeling of GNPs at the cellular level by using the MCNPX Monte Carlo code. A human eye phantom containing the realistic structures and materials was simulated, with a typical tumor located in its corner filled with three different patterns of distribution of GNPs around the nuclei of the cells. The primary photons emit from a COMS eye plaque brachytherapy containing thirteen 131Cs seeds in the vicinity of the tumor.

Results

The study was extended to estimate dose enhancement for various concentration, size, and density of the GNPs accumulated around the nuclei of the tumor. Moreover, the dose delivered to the healthy eye structures for different models has been investigated and discussed. The results show obvious differences between the dose enhancements in the tumor depending on the modeling of GNPs.

Conclusion

The results emphasized that an appropriate small-scale model for the distribution of GNPs in the cell would be of high importance to estimate the degree of dose enhancement in a macroscopic phantom to provide a trustworthy prediction to move towards clinical application.

Background

Increasing the efficiency of the brachytherapy [1] is possible to selectively increase the dose delivered to the target of interest during the radiation. Thanks to the promotion of nanotechnology in medicine, gold nanoparticles (GNPs) have become appropriate candidates to be used in cancer diagnosis and its treatment [2,3,4,5,6,7]. Due to the interaction between X-ray and high-Z constituent atoms of GNPs, the produced secondary electrons cause the increment in the energy deposited to the target. The effects of using GNPs as radiosensitizers on dose enhancement have extensively been studied in both experimental environments and Monte Carlo simulations. The commonly used term, dose enhancement factor (DEF) quantifies the magnitude of dose enhancement. DEF is defined as the ratio of the absorbed dose to the desired volume in the presence of GNPs to the absorbed dose in the same volume in the absence of nanoparticles.

There are excellent works, accomplished using both experimental [2, 8, 9] and simulation-based [10,11,12] approaches, dealing with the investigation of the effect of various parameters governing the DEF value. The simulation-based studies published so far have considered nanoparticles as gold-water (or gold-tissue) mixture or as the homogeneously distributed spheres in the medium [13,14,15,16]. However, the images reported in experimental documents show that the injected nanoparticles tend to accumulate around the cell nucleus [17,18,19,20,21,22]. Accordingly, the fundamental question is whether the exact simulation of the accumulation of nanoparticles in the medium plays a considerable role in the values of DEF, or gold-water mixture and the homogeneous distribution of GNPs in the medium lead to sufficient results.

There are pioneering studies in the literature which have examined the effect of the approach used for simulating nanoparticles on dose enhancement. For example, by simulating both gold-water mixture and homogeneously distributed GNPs, Zhang et al. [23] investigated the difference between the results corresponding to the mentioned models. They found about 36% dose enhancement for the gold-water mixture and about 28% dose enhancement for a homogenous distribution of gold nanospheres in the tumor. A similar study by Martinov and Thomson [24] showed that the homogeneous distribution of GNPs decreases the values of DEF by about 20% compared to the gold-tissue mixture. Moreover, several studies have inquired into the effect of utilizing different models for the distribution of GNPs at the cellular level. Sung et al. [25] investigated the effect of the cell geometry, and radiosensitization and the biological effectiveness of GNPs distributed in the extracellular media on the nucleus DEF. The study accomplished by Cai et al. [26] evaluates the effects of cell model, subcellular location, size, and the number of GNPs per cell, as well as the photon energy on the nuclear dose enhancement factor (NDEF). The results showed that the localization of GNPs in the nucleus and the increment of their number are two important factors for the growth of NDEF. In the study presented by Xie et al. [27], using a cell model with detailed DNA structure in the central nucleus, the physical and biological radiosensitization effect within the nucleus was set up and nanoparticles were distributed in different regions of the cell. They found that the enhancement in energy deposition increases in the case of GNPs with a diameter of 100 nm distributed on the nucleus surface. Considering the widely documented tendency of GNPs to localize around the nucleus, a recently published study has been devoted to investigate the importance of modeling GNPs distribution in a single cell on cellular dose enhancement and to examine the differences in the behavior of the parameters governing GNPs photon radiosensitization of cell nucleus arisen from using different models [28].

Although the mentioned studies try to evaluate the importance of using appropriate small-scale models of the distribution of GNPs in the cell on accurate estimation of cellular dose enhancement, however, the effect of distribution model of GNPs in the cell on macroscopic dose enhancement needs to be investigated in detail and to be well quantified. Inspired by the previous researches, the present work aims to investigate the sensitivity of the simulation results to the cellular modeling of nanoparticles in a typical macroscopic phantom. While the GNPs’ dose-enhanced radiotherapy of the eye has been the subject of several studies, the importance of the location of nanoparticles has not been investigated. To address this issue, a cell model containing nucleus and cytoplasm is designed, and three different patterns for distribution of GNPs are proposed: the gold-tissue mixture, homogeneously distributed GNPs in the medium, and the heterogeneous model with GNPs localized in the cytoplasm around the nucleus. These cells are considered to be located in a typical tumor of a simulated adult eye phantom including anatomical structures filled with realistic compositions. Several radioactive seeds of 131Cs that sit within the plaque with a specific geometry are also simulated in the vicinity of the tumor with the aim of brachytherapy. To show the gold concentration per gram of tissue, the notation of mg/g is introduced. Simulations and particle tracking in the medium have been carried out using MCNPX (Version 2.6.0) [29] Monte Carlo code.

Methods

Monte Carlo simulation

In the present work, we benefit from the Monte Carlo method, known as an appropriate way to simulate complex systems with many coupled degrees of freedom which involve the transport of particles through matter, to find the effect of the distribution of GNPs in the medium on the dose delivered to the tumor. Among the several codes commonly employed in nuclear physics problems, the MCNPX (Version 2.6.0) code has been chosen mainly due to its flexibility and possibility to use for solving different problems. Owing to that the Monte Carlo method is based on repeated random sampling and statistical analysis, the adequate number of histories, between 108 and 109, have been used to obtain the results with acceptable relative uncertainties. The uncertainties corresponding to the results have been calculated by applying the error propagation rules on the Monte Carlo relative errors. Variance reduction methods have not been used in the simulations, and the cut-off energies are the default values of the MCNPX code (1 keV for electrons). Detailed descriptions of the geometries, energies, and models simulated are presented in the following subsections.

131Cs seed (model CS-1)

Brachytherapy is a form of internal radiotherapy in which sealed radioactive sources are placed inside or near to the target volume to kill cancer cells and shrink the tumors. This method is extensively used for the treatment of uveal melanoma, one of the most frequent malignant intraocular tumors in adults. The procedure includes presurgical testing, determining the size of the plaque and the duration of time that the patient will need to wear the plaque by the radiation oncologist, stop taking or change in the dose of some of the medications, plaque implantation, treatment, and plaque removal. The American Brachytherapy Society (ABS) recommends a prescription dose of 85 Gy to the apex of the tumor [30]. It is proposed that a dose rate of 0.60 to 1.05 Gy/h delivers the total dose in 3 to 10 consecutive days [31]. To avoid the damage to the critical structures such as optic nerve, the dose ranging between 30 and 60 Gy has been suggested in the literature [32, 33].

While 125I and 103Pd seeds are commonly used in eye plaque brachytherapy implants, the advantages of 131Cs have made it as an appropriate candidate to be used for this purpose [34, 35]. The half-life of 9.7 days, and the average energy of 30.4 keV, which are respectively shorter than and higher than those of 125I and 103Pd seeds (see Table 1), can be listed as the advantages of this source. According to the literature [28, 37], the primary energy of about 30 keV leads to the highest value of DEF. In the present study, we have simulated the new 131Cs brachytherapy seed model CS-1 developed by IsoRay medical Inc [36]. A single CS-1 seed consists of cylindrical gold wire with a diameter of 0.25 mm and a length of 4.1 mm that is surrounded by a glass and ceramic material with a diameter of 0.65 mm has been coated with 131Cs. The seed is encapsulated inside a titanium tube with a diameter of 0.8 mm and a length of 4.5 mm [38]. According to the data published by National Nuclear Data Center (NNDC), the photons’ spectrum of 131Cs with their intensity is as follows: 4.11 keV (8.6%), 29.461 keV (21.1%), 29.782 keV (38.9%), 33.562 keV (3.63%), 33.624 keV (7.02%), 34.419 keV (2.13%) [39]. The activity of each seed has been considered to be 4 mCi.

Table 1. The data corresponding to the 125I, 103Pd, and 131Cs radioactive sources as commonly used seeds for eye plaque brachytherapy [36]
Full size table

Models for distribution of GNPs

To investigate the effect of the distribution of GNPs on the nucleus dose enhancement, three various models are proposed:

  1. i)

    Gold-tissue mixture. In this model, the material around the nucleus is considered as a homogeneous mixture of tissue and gold with a specific concentration. Owing to that the GNPs are injected in the tumor volume, the compositions of the tissue are considered to be those of eye melanoma [40].

  2. ii)

    Homogeneously distributed GNPs in the medium. In this model, the spheres of GNP have homogeneously been distributed in the whole medium and in the cytoplasm. The remainder volume, including the nucleus and space between the GNPs, was defined as tumor tissue.

  3. iii)

    The heterogeneous model. In this model, the GNPs are localized inside the cytoplasm and around the nucleus.

Schematic diagrams of a single cell, modeled with each of the foregoing patterns, have been shown in Fig. 1. According to the specified dimensions, which have been estimated based on the in vitro images reported by Rezaei Kanavi et al. [41], the tumor of our simulated phantom contains ~ 26 × 106 cells. These cells have been simulated using the availability of simulating the repeated-structures by using the FILL, UNIVERSE, and LAT cards in MCNPX code. Owing that the damage to the DNA in the cell nucleus results in cell death, the physical doses due to the irradiation of photons are scored in the nucleus of each cell using *F8 tally. It is worth to mention that this work does not deal with the calculation of either the dose or damage delivered to the DNA and subcellular structures of a single nucleus. The values of the dose reported in this work are the mean dose delivered to the nuclei of the simulated macroscopic tumor.

Fig. 1.
figure 1

A schematic cross-sectional view of the models considered for distribution of GNPs in the cell. a Gold-tissue mixture. b Homogeneously distributed GNPs in the medium. c Heterogeneous model, with GNPs distributed inside the cytoplasm surrounding the nucleus. In the mixture model, the medium is filled with a homogeneous mixture of eye melanoma and gold. In all three models, the medium is a cubic box of 50 × 50 × 50 μm3, and the nucleus is filled with water

Full size image

The eye phantom

A geometrically and compositionally realistic three-dimensional model of the human eye developed by Lesperance et al. [40] is used as phantom. The dimensions, the shapes of ocular structures, and the elemental composition of the eye are also taken from the mentioned study. This phantom is representative of an adult eye and includes two concentric spheres as the outer and inner sclera, a hollow cylinder as iris, two spherical shells limited by sclera as the cornea, and a volume bounded by two spheres as the lens. In this model, a typical tumor, simulated as a sphere of 12 mm in diameter, is located in the corner of the phantom so that its apex lies 6 mm inward from the inner edge of the sclera. More details can be found in Ref. [42], and a schematic representation of the simulated phantom is sown in Fig. 2a. As the figure shows, the tumor volume is filled with the cells containing GNPs with the models explained in the previous subsection.

Fig. 2
figure 2

a Schematic cross-sectional view of the three-dimensional eye phantom including COMS plaque and tumor. The voxels along the Y-axis are designed for depth-dose calculations. b The simulated 14 mm (inner diameter) eye plaque containing thirteen 131Cs seeds, with various positions and orientations. c The geometry, materials, and dimensions of a single 131Cs seed

Full size image

COMS eye plaque brachytherapy

The radioactive seeds of 131Cs that sit within the plaque with a specific geometry are simulated in the vicinity of the tumor. This method, known as the collaborative ocular melanoma study (COMS) eye plaques, is common in treating eye tumors. The size of COMS eye plaques ranges between 10 and 22 mm, depending on the size of the ocular tumor. In the present study, a typical 14 mm inner diameter COMS eye plaque containing thirteen 131Cs seeds, with various positions and orientations, is used. The seeds are fixed by the silastic carrier with a density of 1.12 g cm−3 and are covered by gold alloy (77% gold, 14% silver, 8% copper, and 1% palladium) with the thickness of 0.5 mm and the density of 17.4 g cm−3 [43]. The simulated plaque is shown in Fig. 2b, and its location near the tumor in our human eye phantom are presented in Fig. 2a. The geometry and dimensions of a single 131Cs seed are also shown in Fig. 2c.

Results and discussion

Figure 3 shows the values of DEF for the tumor filled with our three models for the distribution of 100 nm GNPs with a concentration of 30 mg/g. These values have been calculated in the voxels located in the depth of the phantom (see Fig. 2a). For the same concentration and dimension of GNPs, Fig. 4 shows the depth-dose rate (in Gy/h) in the simulated phantom.

Fig. 3
figure 3

The values of DEF in the depth of human eye phantom (see Fig. 2a) containing a tumor filled with 100 nm GNPs with a concentration of 30 mg/g of different models. Error bars indicate the relative uncertainties

Full size image
Fig. 4
figure 4

The values of the dose rate (Gy/h) calculated in the depth of human eye phantom containing the tumor filled with GNPs of three models for the concentration of 30 mg/g with 100 nm diameter. The depth-dose rates in the absence of GNPs have also been reported. Error bars indicate the relative uncertainties

Full size image

As can be seen, DEFs and dose rates corresponding to the depth of the tumor for the mixture model are larger than those of the homogeneous model. The results agree with those reported in the published works covering the overestimation of dose enhancement caused by gold-tissue mixture compared with homogeneous distribution [23, 24]. Also, DEF values for the heterogeneous model are considerably higher than those of two other models, which can be justified by the accumulation of GNPs in the volume close to the target (nucleus). In other words, in this model, the energy deposition in the nucleus due to the secondary electrons with ranges of several micrometers, generated because of the interaction of irradiated photons with GNPs, is more than two other models. For example, in the heterogeneous model, the dose rate evaluated in the first voxel of the tumor deviates from that of the sclera by about 591.8%. This deviation is about 113.4% and 185.4% for homogeneous and gold-tissue mixture models, respectively. It should be noted that the values calculated for DEF and dose rates depend on details of the models used in the problem, and our results can quantitatively change by changing the geometries and the dimensions. However, the behavior is the same.

The other interesting point in these figures is the reduction of the dose delivered to the surrounding healthy tissues in the presence of GNPs compared with the common eye plaque brachytherapy. This effect is the result of the reduction in the time that the implant needs to be used to reach the therapeutic prescribed radiation dose to the tumor (85 Gy). For example, according to Fig. 4, the apex of the tumor reaches this value in about 25 h in the heterogeneous model. This duration time is about 165 h in the absence of GNPs.

Figure 5 examines the sensitivity of the depth-dose rates to the concentration of GNPs in the tumor for homogeneous and heterogeneous distributions. In these calculations, the tested concentrations are typical values of 7, 15, and 30 mg/g, and the size of GNPs is considered to be the constant value of 100 nm in diameter. The results show that compared with the homogeneous model, the concentration is a more important factor in the dose values for the heterogeneous model. According to the results, the dose rate in the initial voxel of the tumor for the concentration of 30 mg/g deviates from that of 7 mg/g by about 174% and 69% for heterogeneous and homogeneous models, respectively. The more the concentration of GNPs increase, the more the values of dose achieves. In a constant volume, the higher density involves a larger number of GNPs. This causes increasing the interactions between photons and gold which obviously leads to growing the number of secondary electrons effective on dose enhancement. However, owing that the presence of heavy elements with a relatively high concentration into the tissue leads to cellular toxicity, increasing the concentrations to a desirable value is not biologically allowed [44,45,46].

Fig. 5
figure 5

The dose rates calculated in the depth of the eye phantom containing a tumor filled with 100 nm GNPs with concentrations of 7, 15, 30 mg/g for a the homogeneous model. b The heterogeneous model

Full size image

To assess the effect of GNP size on the value of DEF, three diameters of 20, 50, and 100 nm have been investigated. Table 2 reports the DEFs corresponding to the homogeneous and heterogeneous models for concentrations of 7, 15, and 30 mg/g. As the results show, for a given concentration of GNPS in the homogeneous model, the values of DEF increase by increasing the size of the nanoparticle. However, this is not the only key factor for governing the dose enhancement, and the total mass of nanoparticles and the density of GNPs around the nucleus are also important. The results show that in the heterogeneous model in which the GNPs are accumulated in the small volume of cytoplasm around the target volume, the effect of GNPs’ size lessens and regardless of the size chosen for the synthesis of nanoparticles, one approximately will achieve the same dose enhancement in this model. These results, which are in agreement with those of Xie et al. [27], can be justified by increasing the number of low-energy electrons which do not contribute to the dose delivered to the nuclei because of being trapped inside the GNPs and difficulties to escape toward the target of interest.

Table 2 The values of DEF calculated in the nuclei of the human eye tumor for various concentrations and diameters in homogeneous and heterogeneous distributions
Full size table

The depth-dose and depth-DEF curves reported show comparisons between the beam performances in the depth of realistic eye phantom in a specific direction for the designed models. To evaluate planar variations in absorbed dose (or DEF), the results are presented in the form of isodose curves. An example of these curves, including the results corresponding to 100 nm GNPs with a concentration of 30 mg/g, is shown in Fig. 6. As the curves show, in the case of the heterogeneous model, a given location in the tumor receives higher deposited energy compared with the same location in the homogeneous model. However, the more interesting result is that this effect will be inversed in the healthy tissue beyond the tumor, leading to the more preservation of healthy tissues. In other words, the localization of the GNPs in the small volume of cytoplasm concentrates the secondary particles in a smaller region and limits the scattering of the radiations to the non-target tissues.

Fig. 6.
figure 6

The isodose curves for the homogeneous model (solid lines) and heterogeneous model (dashed lines) in the depth of eye phantom including a tumor loaded with 100 nm GNPs of 30 mg/g

Full size image

It is also interesting to examine the dose behavior in the surface perpendicular to the axis along with the depth of the phantom (Fig. 2a). Obviously, the dose curves in the X-Z plane are depth-dependent. As an example, Fig. 7 shows the color-filled contour plots on this surface in the central depth of the tumor. As expected, the doses in both models are approximately symmetric respect to X and Z axes. Moreover, it can be found that the mean dose in the inner layers of the tumor in the heterogeneous model is larger than that of the homogeneous model, which is in agreement with the results of Fig. 5.

Fig. 7.
figure 7

The color-filled contour plots in the X-Z plane in the central depth of the tumor (y = 0.45 cm) for a homogeneous and b heterogeneous models in the eye phantom including a tumor (the radius of the tumor area in the central depth is 0.5 cm) loaded with 100 nm GNPs of 30 mg/g. The color scale (in percent) has also been shown

Full size image

The main goal in radiotherapy is the achievement of high dose deposition in the tumor while sparing the surrounding healthy tissue as much as possible. To assess the performance of gold radiosensitizers in getting close to this goal, the values of DEF in various structures of the eye phantom for two distribution models are calculated. Table 3 reports these values for concentrations of 7, 15, and 30 mg/g of GNPs with a typical diameter of 50 nm.

Table 3 The values of DEF calculated in different structures of the eye including a tumor loaded with 50 nm GNPs with concentrations of 7, 15, and 30 mg/g for homogeneous and heterogeneous models. The last row presents the average DEFs calculated in the nuclei of the tumor
Full size table

The results exhibit that the DEFs calculated in different structures of the eye containing a tumor loaded with GNPs are dependent on the distribution model. As is expected, the values of DEF in the nuclei of the tumor of the heterogeneous model are higher than those of the homogeneous model. For our tested concentrations, this difference ranges between 65.6% and 464.9%. However, the heterogeneous model leads to the decrement of the dose delivered to the sensitive structures such as lens and iris. This decrement increases with increasing the GNPs’ concentration in the tumor. For example, for 50 nm nanoparticles with the concentration of 30 mg/g, the values of DEF corresponding to the homogeneous model in the lens, iris, and vitreous deviates from those of heterogeneous model by about 18.3%, 16.6%, and 16.2%, respectively. A more detailed data on these comparisons can be found in Fig. 8, which presents the deviation from DEF in the heterogeneous model of DEF in the homogeneous model for various eye structures. These data confirm that the healthy eye structures in the homogeneous model are more at risk of receiving high doses, and consequently damage, compared with the same structures in the heterogeneous model.

Fig. 8.
figure 8

The deviation from DEF in the heterogeneous model (DEFhetro) of DEF in the homogeneous model (DEFhomo) for various eye structures. These data correspond to the presence of 50 nm GNPs with concentrations of 7, 15, and 30 mg/g in the tumor

Full size image

Conclusions

The importance of modeling the distribution of GNPs in the cellular level on macroscopic dose enhancement has been investigated. The results show that for a given concentration, the dose delivered to the tumor in the heterogeneous concentration is considerably larger than that of other models. For example for 100 nm nanoparticles with a concentration of 30 mg/g, the DEF in the voxel at the beginning of the tumor of the heterogeneous model deviates from those of homogeneous and gold-tissue mixture models by about 256.8% and 142%, respectively. It was found that the homogeneous model is sensitive to the size of GNPs so that the DEFs increase by increasing of the GNPs size. However, for the heterogeneous model, change in GNPs’ size will not significantly affect the energy deposition in the nucleus. Also, the results indicate that the presence of GNP within the tumor in the heterogeneous model leads to the reduction of the absorbed dose in different structures of the eye, especially in sensitive organs such as lens and iris, compared with the homogeneous model. Moreover, using GNPs reduces the time that the plaque needs to be used to reach the therapeutic prescribed radiation dose to the tumor, which decreases the radiation dose to the surrounding healthy tissues.

The results emphasize the importance of the GNPs’ distribution modeling in the cell on the overall dose enhancement of a realistic human eye phantom, as an example of macroscopic volumes extensively used for studies on gold radiosensitization. However, it is worth mentioning that while the distribution of the GNPs around the nucleus is not completely uniform, it has not been considered in our heterogeneous model. Taking into account the non-uniform distribution of GNPs in the cytoplasm offers new ways to accomplish more researches in this field.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

References

  1. Skowronek J. Current status of brachytherapy in cancer treatment–short overview. J Contemp Brachytherapy. 2017;9:581–9.

    PubMed  PubMed Central  Google Scholar 

  2. Kuncic Z, Sandrine L. Nanoparticle radio-enhancement: principles, progress and application to cancer treatment. Phys Med Bio. 2018;63:02TR01.

    Google Scholar 

  3. Hainfeld JF, Avraham Dilmanian F, Slatkin DN, Smilowitz HM. Radiotherapy enhancement with gold nanoparticles. J Pharm Pharmacol. 2008;60:977–85.

    CAS  PubMed  Google Scholar 

  4. Kwatra D, Venugopal A, Anant S. Nanoparticles in radiation therapy: a summary of various approaches to enhance radiosensitization in cancer. Transl Cancer Res. 2013;2(4):330–42.

    CAS  Google Scholar 

  5. Jaffray DA. Image-guided radiotherapy: from current concept to future perspectives. Nat Rev Clin Oncol. 2012 Dec;9(12):688.

    CAS  PubMed  Google Scholar 

  6. Lim ZZ, Li JE, Ng CT, Yung LY, Bay BH. Gold nanoparticles in cancer therapy. Acta Pharm Sinic. 2011;32(8):983–90.

    CAS  Google Scholar 

  7. Jain S, Hirst DG, O'sullivan JM. Gold nanoparticles as novel agents for cancer therapy. Brit J Radiol. 2012;85(1010):101–13.

    CAS  PubMed  Google Scholar 

  8. Ngwa W, Korideck H, Kassis AI, Kumar R, Sridhar S, Makrigiorgos GM, Cormack RA. In vitro radiosensitization by gold nanoparticles during continuous low-dose-rate gamma irradiation with I-125 brachytherapy seeds. Nanomedicine: Nanotechnology, Biology and Medicine. 2013;9(1):25–7.

    CAS  Google Scholar 

  9. Huang X. Gold nanoparticles used in cancer cell diagnostics, selective photothermal therapy and catalysis of NADH oxidation reaction (Doctoral dissertation, Georgia Institute of Technology).

  10. Jones BL, Krishnan S, Cho SH. Estimation of microscopic dose enhancement factor around gold nanoparticles by Monte Carlo calculations. Med Phys. 2010;37(7Part1):3809–16.

    CAS  PubMed  Google Scholar 

  11. Cho SH. Estimation of tumour dose enhancement due to gold nanoparticles during typical radiation treatments: a preliminary Monte Carlo study. Phys Med Biol. 2005;50(15):N163.

    PubMed  Google Scholar 

  12. Cho SH, Jones BL, Krishnan S. The dosimetric feasibility of gold nanoparticle-aided radiation therapy (GNRT) via brachytherapy using low-energy gamma-/x-ray sources. Phys Med Biol. 2009;54(16):4889.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Corde S, Joubert A, Adam JF, Charvet AM, Le Bas JF, Esteve F, Elleaume H, Balosso J. Synchrotron radiation-based experimental determination of the optimal energy for cell radiotoxicity enhancement following photoelectric effect on stable iodinated compounds. Brit J Cancer. 2004;91(3):544–51.

    CAS  PubMed  Google Scholar 

  14. Al-Musywel HA, Laref A. Effect of gold nanoparticles on radiation doses in tumor treatment: A monte carlo study. Laser Med Sci. 2017;32(9):2073–80.

    CAS  Google Scholar 

  15. Koger B, Kirkby C. Optimization of photon beam energies in gold nanoparticle enhanced arc radiation therapy using Monte Carlo methods. Phys Med Biol. 2016;61(24):8839.

    CAS  PubMed  Google Scholar 

  16. Hwang C, Kim JM, Kim J. Influence of concentration, nanoparticle size, beam energy, and material on dose enhancement in radiation therapy. J Radiat Res. 2017;58(4):405–11.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Navas-Moreno M, Mehrpouyan M, Chernenko T, Candas D, Fan M, Li JJ, Yan M, Chan JW. Nanoparticles for live cell microscopy: A surface-enhanced Raman scattering perspective. Sci Rep. 2017;7(1):1-0.

  18. Zhang X, Shastry S, Bradforth SE, Nadeau JL. Nuclear uptake of ultrasmall gold-doxorubicin conjugates imaged by fluorescence lifetime imaging microscopy (FLIM) and electron microscopy. Nanoscale. 2015;7(1):240–51.

    CAS  PubMed  Google Scholar 

  19. Liu F, Lou J, Hristov D. X-Ray responsive nanoparticles with triggered release of nitrite, a precursor of reactive nitrogen species, for enhanced cancer radiosensitization. Nanoscale. 2017;9(38):14627–34.

    CAS  PubMed  Google Scholar 

  20. Coulter JA, Jain S, Butterworth KT, Taggart LE, Dickson GR, McMahon SJ, Hyland WB, Muir MF, Trainor C, Hounsell AR, O’Sullivan JM. Cell type-dependent uptake, localization, and cytotoxicity of 1.9 nm gold nanoparticles. Int J Nanomed. 2012;7:2673.

    CAS  Google Scholar 

  21. Kong T, Zeng J, Wang X, Yang X, Yang J, McQuarrie S, McEwan A, Roa W, Chen J, Xing JZ. Enhancement of radiation cytotoxicity in breast-cancer cells by localized attachment of gold nanoparticles. Small. 2008;4(9):1537–43.

    CAS  PubMed  Google Scholar 

  22. Lechtman E, Mashouf S, Chattopadhyay N, Keller BM, Lai P, Cai Z, Reilly RM, Pignol JP. A Monte Carlo-based model of gold nanoparticle radiosensitization accounting for increased radiobiological effectiveness. Phys Med Biol. 2013;58(10):3075.

    CAS  PubMed  Google Scholar 

  23. Zhang SX, Gao J, Buchholz TA, Wang Z, Salehpour MR, Drezek RA, Yu TK. Quantifying tumor-selective radiation dose enhancements using gold nanoparticles: a monte carlo simulation study. Biomed Microdevices. 2009;11(4):925.

    CAS  PubMed  Google Scholar 

  24. Martinov MP, Thomson RM. Heterogeneous multiscale Monte Carlo simulations for gold nanoparticle radiosensitization. Med Phys. 2017;44(2):644–53.

    CAS  PubMed  Google Scholar 

  25. Sung W, Ye SJ, McNamara AL, McMahon SJ, Hainfeld J, Shin J, Smilowitz HM, Paganetti H, Schuemann J. Dependence of gold nanoparticle radiosensitization on cell geometry. Nanoscale. 2017;9(18):5843–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Cai Z, Pignol JP, Chattopadhyay N, Kwon YL, Lechtman E, Reilly RM. Investigation of the effects of cell model and subcellular location of gold nanoparticles on nuclear dose enhancement factors using Monte Carlo simulation. Med Phys. 2013;40(11):114101.

    PubMed  Google Scholar 

  27. Xie WZ, Friedland W, Li WB, Li CY, Oeh U, Qiu R, Li JL, Hoeschen C. Simulation on the molecular radiosensitization effect of gold nanoparticles in cells irradiated by x-rays. Phys Med Biol. 2015;60(16):6195.

    CAS  PubMed  Google Scholar 

  28. Rasouli FS, Masoudi SF. Monte Carlo investigation of the effect of gold nanoparticles’ distribution on cellular dose enhancement. Radiat Phys Chem. 2019;158:6–12.

    CAS  Google Scholar 

  29. Pelowitz, B., 2008. MCNPX User's Manual Version 2.6.0. Los Alamos national laboratory (Report LA-CP-07-1473).

  30. Nag S, Quivey JM, Earle JD, Followill D, Fontanesi J, Finger PT, Society AB. The American Brachytherapy Society recommendations for brachytherapy of uveal melanomas. Int J Radiat Oncol. 2003;56(2):544–55.

    Google Scholar 

  31. Rivard MJ, Coursey BM, DeWerd LA, Hanson WF, Saiful Huq M, Ibbott GS, Mitch MG, Nath R, Williamson JF. Update of AAPM Task Group No. 43 Report: A revised AAPM protocol for brachytherapy dose calculations. Med Phys. 2004;31(3):633–74.

    PubMed  Google Scholar 

  32. Miguel D, de Frutos JM, López-Lara F, Alonso P, Diezhandino P, Saornil MA, García C. Treatment planning considerations for 125I eye plaque brachytherapy. J Contem Brachyther. 2019;11(3):280.

    Google Scholar 

  33. Gragoudas ES, Li W, Lane AM, Munzenrider J, Egan KM. Risk factors for radiation maculopathy and papillopathy after intraocular irradiation. Ophthalmology. 1999;106(8):1571–8.

    CAS  PubMed  Google Scholar 

  34. Wittman RS, Fisher DR. Multiple-estimate Monte Carlo calculation of the dose rate constant for a cesium-131 interstitial brachytherapy seed. Med Phys. 2007 Jan;34(1):49–54.

    CAS  PubMed  Google Scholar 

  35. Murphy MK, Piper RK, Greenwood LR, Mitch MG, Lamperti PJ, Seltzer SM, Bales MJ, Phillips MH. Evaluation of the new cesium-131 seed for use in low-energy x-ray brachytherapy. Med Phys. 2004 Jun;31(6):1529–38.

    CAS  PubMed  Google Scholar 

  36. IsoRay Medical Inc., CS-1 131Cs product information sheet and company website http://www.isoray.com/.

  37. Rahman WN, Corde S, Yagi N, Aziz SA, Annabell N, Geso M. Optimal energy for cell radiosensitivity enhancement by gold nanoparticles using synchrotron-based monoenergetic photon beams. Int J Nanomed. 2014;9:2459.

    Google Scholar 

  38. Chen Z, Bongiorni P, Nath R. Dose rate constant of a cesium-131 interstitial brachytherapy seed measured by thermoluminescent dosimetry and gamma-ray spectrometry. Medical physics. 2005;32(11):3279–85.

    CAS  PubMed  Google Scholar 

  39. NUDAT 2.0, National Nuclear Data Center, Brookhaven National Laboratory, http://www.nndc.bnl.gov/nudat2/index.jsp (last accessed 27 January 2007).

  40. Lesperance M, Inglis-Whalen M, Thomson RM. Model-based dose calculations for COMS eye plaque brachytherapy using an anatomically realistic eye phantom. Med Phys. 2014;41(2):021717.

    PubMed  Google Scholar 

  41. Kanavi MR, Asadi S, Balagholi S, Alikarami F, Nosrati H, Ahmadieh H. Gamma irradiation of ocular melanoma and lymphoma cells in the presence of gold nanoparticles: in vitro study. J Appl Clin Med Phys. 2018;19(3):268–75.

    PubMed  PubMed Central  Google Scholar 

  42. Rasouli FS, Masoudi SF, Keshazare S, Jette D. Effect of elemental compositions on Monte Carlo dose calculations in proton therapy of eye tumors. Radiat Phys Chem. 2015;117:112–9.

    CAS  Google Scholar 

  43. Zhang H, Martin D, Chiu-Tsao ST, Meigooni A, Thomadsen BR. A comprehensive dosimetric comparison between 131Cs and 125I brachytherapy sources for COMS eye plaque implant. Brachytherapy. 2010;9(4):362–72.

    PubMed  Google Scholar 

  44. Zhang XD, Guo ML, Wu HY, Sun YM, Ding YQ, Feng X, Zhang LA. Irradiation stability and cytotoxicity of gold nanoparticles for radiotherapy. Int J Nanomed. 2009;4:165.

    CAS  Google Scholar 

  45. Coradeghini R, Gioria S, García CP, Nativo P, Franchini F, Gilliland D, Ponti J, Rossi F. Size-dependent toxicity and cell interaction mechanisms of gold nanoparticles on mouse fibroblasts. Toxicol Lett. 2013;217(3):205–16.

    CAS  PubMed  Google Scholar 

  46. Haume K, Rosa S, Grellet S, Śmiałek MA, Butterworth KT, Solov’yov AV, Prise KM, Golding J, Mason NJ. Gold nanoparticles for cancer radiotherapy: a review. Cancer Nanotechnol. 2016;7(1):8.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The present research work has not been founded by any organization.

Author information

Authors and Affiliations

  1. Department of Physics, K.N. Toosi University of Technology, P.O. Box 15875-4416, Tehran, Iran

    S. Farhad Masoudi, Fahimeh S. Daryabari & Fatemeh S. Rasouli

Authors
  1. S. Farhad Masoudi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Fahimeh S. Daryabari
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fatemeh S. Rasouli
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

All authors conceived the project, prepared the data and the figures, and analyzed the results. They have collaborated in writing the main text of the manuscript. All authors read and approved the final article.

Corresponding author

Correspondence to S. Farhad Masoudi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Masoudi, S.F., Daryabari, F.S. & Rasouli, F.S. Distribution modeling of nanoparticles for brachytherapy of human eye tumor. EJNMMI Phys 7, 53 (2020). https://doi.org/10.1186/s40658-020-00321-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40658-020-00321-y

Keywords