Search

One of the main appeals of using gold nanoparticles (GNPs) as radiosensitizers is that their surface coatings can be altered to manipulate their pharmacokinetic properties. However, Monte Carlo studies of GNP dosimetry tend to neglect these coatings, potentially changing the dosimetric results. This study quantifies the dosimetric effects of including a polyethylene glycol (PEG) surface coating on GNPs over both nanoscopic and microscopic ranges. Two dosimetric scales were explored using PENELOPE Monte Carlo simulations. In microscopic simulations, 500-1000 GNPs, with and without coatings, were placed in cavities of side lengths 0.8-4 µm, and the reduction of dose deposited to surrounding medium within these volumes due to the coating was quantified. Including PEG surface coatings of up to 20 nm thickness resulted in reductions of up to 7.5%, 4.0%, and 2.0% for GNP diameters of 10, 20, and 50 nm, respectively. Nanoscopic simulations observed the dose falloff in the first 500 nm surrounding a single GNP both with and without surface coatings of various thicknesses. Over the first 500 nm surrounding a single GNP, the presence of a PEG surface coating reduced dose by 5-26%, 8-28%, 8-30%, and 8-34% for 2, 10, 20, and 50 nm diameter GNPs, respectively, for various energies and coating thicknesses. Reductions in dose enhancement due to the inclusion of a GNP surface coating are non-negligible and should be taken into consideration when investigating GNP dose enhancement. Further studies should be carried out to investigate the biological effects of these coatings.

To account for particle interactions in the extracellular (physical) environment, information from the cell-level Monte Carlo damage simulation (MCDS) for DNA double strand break (DSB) induction has been integrated into the general purpose Monte Carlo N-particle (MCNP) radiation transport code system. The effort to integrate these models is motivated by the need for a computationally efficient model to accurately predict particle relative biological effectiveness (RBE) in cell cultures and in vivo. To illustrate the approach and highlight the impact of the larger scale physical environment (e.g. establishing charged particle equilibrium), we examined the RBE for DSB induction (RBEDSB) of x-rays, (137)Cs γ-rays, neutrons and light ions relative to γ-rays from (60)Co in monolayer cell cultures at various depths in water. Under normoxic conditions, we found that (137)Cs γ-rays are about 1.7% more effective at creating DSB than γ-rays from (60)Co (RBEDSB  =  1.017) whereas 60-250 kV x-rays are 1.1 to 1.25 times more efficient at creating DSB than (60)Co. Under anoxic conditions, kV x-rays may have an RBEDSB up to 1.51 times as large as (60)Co γ-rays. Fission neutrons passing through monolayer cell cultures have an RBEDSB that ranges from 2.6 to 3.0 in normoxic cells, but may be as large as 9.93 for anoxic cells. For proton pencil beams, Monte Carlo simulations suggest an RBEDSB of about 1.2 at the tip of the Bragg peak and up to 1.6 a few mm beyond the Bragg peak. Bragg peak RBEDSB increases with decreasing oxygen concentration, which may create opportunities to apply proton dose painting to help address tumor hypoxia. Modeling of the particle RBE for DSB induction across multiple physical and biological scales has the potential to aid in the interpretation of laboratory experiments and provide useful information to advance the safety and effectiveness of hadron therapy in the treatment of cancer.

Purpose: Radiation damage to mitochondria has been shown to alter cellular processes and even lead to apoptosis. Gold nanoparticles (AuNPs) may be used to enhance these effects in scenarios where they collect on the outer membranes of mitochondria. A Monte Carlo (MC) approach is used to estimate mitochondrial dose enhancement under a variety of conditions., Methods: The penelope MC code was used to generate dose distributions resulting from photons striking a 13 nm diameter AuNP with various thicknesses of water-equivalent coatings. Similar dose distributions were generated with the AuNP replaced by water so as to estimate the gain in dose on a microscopic scale due to the presence of AuNPs within an irradiated volume. Models of mitochondria with AuNPs affixed to their outer membrane were then generated-considering variation in mitochondrial size and shape, number of affixed AuNPs, and AuNP coating thickness-and exposed (in a dose calculation sense) to source spectra ranging from 6 MV to 90 kVp. Subsequently dose enhancement ratios (DERs), or the dose with the AuNPs present to that for no AuNPs, for the entire mitochondrion and its components were tallied under these scenarios., Results: For a representative case of a 1000 nm diameter mitochondrion affixed with 565 AuNPs, each with a 13 nm thick coating, the mean DER over the whole organelle ranged from roughly 1.1 to 1.6 for the kilovoltage sources, but was generally less than 1.01 for the megavoltage sources. The outer membrane DERs remained less than 1.01 for the megavoltage sources, but rose to 2.3 for 90 kVp. The voxel maximum DER values were as high as 8.2 for the 90 kVp source and increased further when the particles clustered together. The DER exhibited dependence on the mitochondrion dimensions, number of AuNPs, and the AuNP coating thickness., Conclusions: Substantial dose enhancement directly to the mitochondria can be achieved under the conditions modeled. If the mitochondrion dose can be directly enhanced, as these simulations show, this work suggests the potential for both a tool to study the role of mitochondria in cellular response to radiation and a novel avenue for radiation therapy in that the mitochondria may be targeted, rather than the nuclear DNA.

Purpose: In radiation biology, rats are often irradiated, but the precise dose distributions are often lacking, particularly in areas that receive scatter radiation. We used a non-dedicated set of resources to calculate detailed dose distributions, including doses to peripheral organs well outside of the primary field, in common rat exposure settings., Materials and Methods: We conducted a detailed dose reconstruction in a rat through an analog to the conventional human treatment planning process. The process consisted of: (i) Characterizing source properties of an X-ray irradiator system, (ii) acquiring a computed tomography (CT) scan of a rat model, and (iii) using a Monte Carlo (MC) dose calculation engine to generate the dose distribution within the rat model. We considered cranial and liver irradiation scenarios where the rest of the body was protected by a lead shield. Organs of interest were the brain, liver and gonads. The study also included paired scenarios where the dose to adjacent, shielded rats was determined as a potential control for analysis of bystander effects., Results: We established the precise doses and dose distributions delivered to the peripheral organs in single and paired rats. Mean doses to non-targeted organs in irradiated rats ranged from 0.03-0.1% of the reference platform dose. Mean doses to the adjacent rat peripheral organs were consistent to within 10% those of the directly irradiated rat., Conclusions: This work provided details of dose distributions in rat models under common irradiation conditions and established an effective scenario for delivering only scattered radiation consistent with that in a directly irradiated rat.

Due to the higher LET of kilovoltage (kV) radiation, there is potential for an increase in relative biological effectiveness (RBE) of absorbed doses of radiation from kV cone beam computed tomography (CBCT) sources in reference to megavoltage or Co-60 doses. This work develops a method for accurately coupling a Monte Carlo (MC) radiation transport code (PENELOPE) with the damage simulation (MCDS) to predict relative numbers of DNA double strand breaks (DSBs). The MCDS accounts for slowing down of electrons and delta ray production within the cell nucleus; however, determining the spectrum of electrons incident on the cell nucleus from photons interacting in a larger region of tissue is not trivial. PENELOPE simulations were conducted with a novel tally algorithm invoked where electrons incident on a detection material were tracked and both the incident energy and the final deposited dose were recorded. The DSB yield predicted by a set of MCDS runs of monoenergetic electrons was then looked up in a table and weighted by the specific energy of the incident electron. Our results indicate that the RBE for DSB induction is 1.1 for diagnostic x-rays with energies from 80 to 125 kVp. We found no significant change in RBE with depth or filtration. The predicted absolute DSB yields are about three times lower for cells irradiated under anoxic conditions than the yield in cells irradiated under normoxic (5%) or fully aerobic (100%) conditions. However, oxygen concentration has a negligible (± 0.02) effect on the RBE of kV CBCT x-rays.

Purpose: The magnetic fields of linac-MR systems modify the path of contaminant electrons in photon beams, which alters patient skin dose. To accurately quantify the magnitude of changes in skin dose, the authors use Monte Carlo calculations that incorporate realistic 3D magnetic field models of longitudinal and transverse linac-MR systems., Methods: Finite element method (FEM) is used to generate complete 3D magnetic field maps for 0.56 T longitudinal and transverse linac-MR magnet assemblies, as well as for representative 0.5 and 1.0 T Helmholtz MRI systems. EGSnrc simulations implementing these 3D magnetic fields are performed. The geometry for the BEAMnrc simulations incorporates the Varian 600C 6 MV linac, magnet poles, the yoke, and the magnetic shields of the linac-MRIs. Resulting phase-space files are used to calculate the central axis percent depth-doses in a water phantom and 2D skin dose distributions for 70 μm entrance and exit layers using DOSXYZnrc. For comparison, skin doses are also calculated in the absence of magnetic field, and using a 1D magnetic field with an unrealistically large fringe field. The effects of photon field size, air gap (longitudinal configuration), and angle of obliquity (transverse configuration) are also investigated., Results: Realistic modeling of the 3D magnetic fields shows that fringe fields decay rapidly and have a very small magnitude at the linac head. As a result, longitudinal linac-MR systems mostly confine contaminant electrons that are generated in the air gap and have an insignificant effect on electrons produced further upstream. The increase in the skin dose for the longitudinal configuration compared to the zero B-field case varies from ∼1% to ∼14% for air gaps of 5-31 cm, respectively. (All dose changes are reported as a % of D(max).) The increase is also field-size dependent, ranging from ∼3% at 20 × 20 cm(2) to ∼11% at 5 × 5 cm(2). The small changes in skin dose are in contrast to significant increases that are calculated for the unrealistic 1D magnetic field. For the transverse configuration, the entrance skin dose is equal or smaller than that of the zero B-field case for perpendicular beams. For a 10 × 10 cm(2) oblique beam the transverse magnetic field decreases the entry skin dose for oblique angles less than ±20° and increases it by no more than 10% for larger angles up to ±45°. The exit skin dose is increased by 42% for a 10 × 10 cm(2) perpendicular beam, but appreciably drops and approaches the zero B-field case for large oblique angles of incidence., Conclusions: For longitudinal linac-MR systems only a small increase in the entrance skin dose is predicted, due to the rapid decay of the realistic magnetic fringe fields. For transverse linac-MR systems, changes to the entrance skin dose are small for most scenarios. For the same geometry, on the exit side a fairly large increase is observed for perpendicular beams, but significantly drops for large oblique angles of incidence. The observed effects on skin dose are not expected to limit the application of linac-MR systems in either the longitudinal or transverse configuration.

Helical tomotherapy delivers intensity modulated radiation therapy using a binary multileaf collimator (MLC) to modulate a fan beam of radiation. This delivery occurs while the linac gantry and treatment couch are both in constant motion, so the beam describes, from a patient/phantom perspective, a spiral or helix of dose. The planning system models this continuous delivery as a large number (51) of discrete gantry positions per rotation, and given the small jaw/fan width setting typically used (1 or 2.5 cm) and the number of overlapping rotations used to cover the target (pitch often <0.5), the treatment planning system (TPS) potentially employs a very large number of static beam directions and leaf opening configurations to model the modulated fields. All dose calculations performed by the system employ a convolution/superposition model. In this work the authors perform a full Monte Carlo (MC) dose calculation of tomotherapy deliveries to phantom computed tomography (CT) data sets to verify the TPS calculations. All MC calculations are performed with the EGSnrc-based MC simulation codes, BEAMnrc and DOSXYZnrc. Simulations are performed by taking the sinogram (leaf opening versus time) of the treatment plan and decomposing it into 51 different projections per rotation, as does the TPS, each of which is segmented further into multiple MLC opening configurations, each with different weights that correspond to leaf opening times. Then the projection is simulated by the summing of all of the opening configurations, and the overall rotational treatment is simulated by the summing of all of the projection simulations. Commissioning of the source model was verified by comparing measured and simulated values for the percent depth dose and beam profiles shapes for various jaw settings. The accuracy of the MLC leaf width and tongue and groove spacing were verified by comparing measured and simulated values for the MLC leakage and a picket fence pattern. The validated source and MLC configuration were then used to simulate a complex modulated delivery from fixed gantry angle. Further, a preliminary rotational treatment plan to a delivery quality assurance phantom (the "cheese" phantom) CT data set was simulated. Simulations were compared with measured results taken with an A1SL ionization chamber or EDR2 film measurements in a water tank or in a solid water phantom, respectively. The source and MLC MC simulations agree with the film measurements, with an acceptable number of pixels passing the 2%/1 mm gamma criterion. 99.8% of voxels of the MC calculation in the planning target volume (PTV) of the preliminary plan passed the 2%/2 mm gamma value test. 87.0% and 66.2% of the voxels in two organs at risk (OARs) passed the 2%/2 mm tests. For a 3%/3 mm criterion, the PTV and OARs show 100%, 93.2%, and 86.6% agreement, respectively. All voxels passed the gamma value test with a criterion of 5%/3 mm. The Tomo-Therapy TPS showed comparable results.

Existing studies have suggested some debate on whether the quality of radiation that delivers dose outside of the primary field of a radiotherapy photon beam can be considered the same as that inside the primary field. We used a Monte Carlo approach to simulate the electron fluence differential in energy inside a water phantom in response to irradiation by a 6 MV photon beam. The goal was to quantify how significantly the electron fluence changes when moving from a volume exposed to the primary field to one outside of the primary field, and understand any potential biological implications. We scored the electron fluence outwards in annular volumes in response to a 5 cm radius 6 MV beam and at the central axis in response to a rectangular 6 MV beam partially blocked by an MLC. The resulting fluence spectra were compared to different low-LET sources for which biological response in the form of chromosomal aberrations has been published. Our results show a significant increase in the low energy component of the fluence spectra outside of the primary field, which increases the mean LET to values similar to that seen in response to a 137Cs photon source. In turn, it is shown that this has the potential to increase the RBE.

Helical tomotherapy is an increasingly common form of intensity modulated radiation therapy that allows for image guided adaptive radiotherapy. Its treatment planning system (TPS) uses a convolution superposition algorithm for dose distribution calculations. The accuracy of this algorithm in the presence of heterogeneities was evaluated against Monte Carlo (MC) calculations and measurements. This work performed BEAMnrc-and DOSXYZnrc-based MC dose calculations of tomotherapy deliveries to a CIRS anthropomorphic heterogeneous phantom with typical clinical inverse planning and delivery settings. Point measurements with A1SL ion chambers and relative measurements with Kodak EDR2 film were carried out in the phantom. The experimental results were used to evaluate both the TPS and MC dose calculations. Furthermore, the dose distribution for a clinical head-and-neck cancer plan was calculated on the TPS and MC systems. The results support this MC system as a viable option for the accurate simulation of the tomotherapy process in the presence of heterogeneities where direct measurement may not be practical. Ion chamber measurements in the CIRS phantom suggested the TPS has an average relative difference of 2.3%, with the largest difference being -4.1% in one of the organs at risk. The MC system accurately predicted the dose to these measurement points within statistical uncertainty. The film measurements in the CIRS phantom demonstrated 90.7% (of pixels) agreed with the MC system using a ±3%/3 mm acceptance criteria, where only 50.3% agreed with the TPS. In the clinical head-and-neck cancer plan evaluation where MC served as a reference against which to compare the TPS result, an average of 92.7% of the voxels within volumes of interest passed a 3%/3 mm criteria. The PTV54 region showed the worst agreement with 85.4% of the volume passing the 3%/3 mm criteria. In general, the ±3%/3 mm criterion was found to be a challenge for the TPS in the presence of lung inhomogeneity. [ABSTRACT FROM AUTHOR]