Abstract
Several radionuclides used in medical imaging emit Auger electrons, which, depending on the targeting strategy, either may be exploited for therapeutic purposes or may contribute to an unintentional mean absorbed dose burden. In this study, the virtues of 12 Auger electron–emitting radionuclides were evaluated in terms of cellular S values in concentric and eccentric cell–nucleus arrangements and by comparing their dose-point kernels. Methods: The Monte Carlo code PENELOPE was used to transport the full particulate spectrum of 67Ga, 80mBr, 89Zr, 90Nb, 99mTc, 111In, 117mSn, 119Sb, 123I, 125I, 195mPt, and 201Tl by means of event-by-event simulations. Cellular S values were calculated for varying cell and nucleus radii, and the effects of cell eccentricity on S values were evaluated. Dose-point kernels were determined up to 30 μm. Energy deposition at DNA scales was also compared with an α emitter, 223Ra. Results: PENELOPE-determined S values were generally within 10% of MIRD values when the source and target regions strongly overlapped, that is, S(nucleus←nucleus) configurations, but greater differences were noted for S(nucleus←cytoplasm) and S(nucleus←cell surface) configurations. Cell eccentricity had the greatest effect when the nucleus was small, compared with the cell size, and when the radiation sources were on the cell surface. Dose-point kernels taken together with the energy spectra of the radionuclides can account for some of the differences in energy deposition patterns between the radionuclides. The energy deposition of most Auger electron emitters at DNA scales of 2 nm or less exceeded that of a monoenergetic 5.77-MeV α particle, but not for 223Ra. Conclusion: A single-cell dosimetric approach is required to evaluate the efficacy of individual radionuclides for theranostic purposes, taking cell geometry into account, with internalizing and noninternalizing targeting strategies.
- Auger electron emitters
- cellular dosimetry
- targeted radiotherapy
- theranostics
- Monte Carlo simulation
- S values
- dose point kernels
The therapeutic rationale for molecularly targeted radiotherapy is the selective delivery of a radionuclide to tumor cells via a targeting moiety, thereby enhancing the therapeutic index of the agent. Several Auger electron–emitting radionuclides have been proposed for molecularly targeted radiotherapy of small metastases and disseminated cancer cells, with some promising clinical results (1–3). These radionuclides are well suited to serving as molecularly targeted radiotherapy agents because of the extremely short range in matter (nanometers to a few micrometers) of the low-energy, intermediate–linear-energy-transfer (LET) Auger and Coster–Kronig electrons they emit (4). These electrons account for high energy deposition in the immediate vicinity of the decay site, and because of their short range, irradiation of normal neighboring cells is limited, thus reducing nonspecific radiotoxicity. In addition, radiation emitted during the nuclear decay can be exploited for imaging purposes either with SPECT (in the case of γ rays in the energy range of 70–360 keV or Bremsstrahlung imaging for pure β− emitters) or PET (in the case of annihilation photons), thus making Auger electron–emitting radionuclides ideal as theranostic agents (5,6).
When evaluating Auger electron–emitting radionuclides as potential theranostic agents, one should address the following aspects. The therapeutic efficacy of a radionuclide depends on the total number of electrons emitted per decay (including Auger electrons, Coster–Kronig electrons, and internal conversion electrons, as well as β particles), along with the total energy released. When the total energy released is carried by a small number of relatively high-energy long-range electrons, then targeting of the nucleus from the cytoplasm or cell surface is possible (7). Another key factor is the physical half-life. Because the maximum theoretic specific activity of a radionuclide is inversely proportional to physical half-life, a prolonged physical half-life may result in redistribution of the therapeutic agent before sufficient decays have occurred to cause lethal damage. It is also important to consider the ratio of penetrating (x rays and γ rays) to nonpenetrating (electrons or β particles; i.e., p/e) forms of ionizing radiation because the moderate-to-high-energy but low-LET γ emissions from some radionuclides may irradiate and potentially kill nontargeted normal cells. It has therefore been proposed that an ideal therapeutic radionuclide should have a p.e. ratio of 2 (8). Conversely, for imaging purposes a high proportion of γ emissions is required, which poses a tradeoff between sparing healthy tissue surrounding the target region and providing an adequate mean absorbed dose to the target region.
The most critical point to consider in the use of Auger electron–emitting radionuclides for molecularly targeted radiotherapy is the short range of Auger and Coster–Kronig electrons, as this necessitates intranuclear accumulation if therapeutic effect is to be maximized. Since the dimensions of the different DNA condensation states (e.g., chromatin fibers, nucleosomes, and double-stranded DNA) are within the range of typical Auger and Coster–Kronig electrons, nuclear incorporation leads to extreme radiotoxicity, resembling high-LET radiation with a relative biological effectiveness similar to that achieved by α-emitting radionuclides (9). Several strategies have been proposed to optimally localize the radionuclides with respect to the sensitive targets in cells (5,10,11). However, recent observations suggest that nuclear accumulation may not be required in order for an Auger electron-emitter to produce the high-LET type of radiotoxicity (12). In contrast, radionuclides bound outside the cell nucleus (e.g., in the cytoplasm), on the cellular membrane, or extracellularly do not produce severe lethal effects and have a relative biological effectiveness comparable to that observed for low-LET radiation (13).
Regardless of the targeting strategy adopted, cellular geometry can influence the mean absorbed dose to the nucleus and thus the biological effect of an Auger electron–emitting theranostic agent. The S value estimates provided by the MIRD Committee (14) assume a spheric cell geometry, but it is conceded that cellular geometry may affect these values. This was demonstrated by Nettleton et al. (15), who noted differences between S value calculations in spheric and ellipsoid cell geometries, especially toward the edge of the cell. Considering that many cells exhibit irregular geometries and eccentric cell–nucleus arrangements, the dose to the nucleus may be over- or underestimated when MIRD-tabulated S values are used.
The aim of this study was to evaluate the properties of 12 Auger electron–emitting radionuclides that have been proposed as theranostic agents, namely 67Ga, 80mBr, 89Zr1, 90Nb1, 99mTc, 111In, 117mSn, 119Sb1, 123I, 125I, 195mPt, and 201Tl (Supplemental Table 1; supplemental materials are available at http://jnm.snmjournals.org) (16–18). (89Zr, 90Nb, and 119Sb are not included in the MIRD monograph (19).)
For these radionuclides, we evaluated S values from Monte Carlo simulations, the effect of cellular geometry and eccentric cell or nucleus arrangements on S values, dose-point kernels (DPKs), and energy deposition on a DNA scale in comparison with an α emitter, 223Ra.
MATERIALS AND METHODS
Monte Carlo Simulations: The PENELOPE Code
The S values and DPKs were calculated with the general-purpose Monte Carlo code PENELOPE (20). PENELOPE simulates the coupled electron–photon transport in arbitrary materials from 50 eV to 1 GeV. The simulation is controlled by 7 user-defined parameters: Eabs(1), Eabs(2), Eabs(3), C1, C2, Wcc, and Wcr. The first 3 parameters fix the absorption energy for electrons, photons, and positrons, defining the cutoff energy below which simulation is discontinued (50 eV) and the residual energy of the particle deposited locally. The remaining parameters control the mixed simulation algorithm for the transport of electrons and positrons. To force detailed (event-by-event) simulation, the latter parameters were set to zero (21).
The cell model consists of 2 homogeneous spheres of liquid water (mass density ρ = 1 gcm−3), representing the cell and its nucleus, immersed in an infinite water medium (Fig. 1). Cell (RC) and nucleus (RN) radii combinations as tabulated by the MIRD Committee (14) were considered, and the list of RC and RN values was expanded to include larger cell geometries (up to RC = 12 μm and RN = 11 μm) (19). Typically, 2 × 109 primary particles were simulated in each run.
S Value Calculations
Monte Carlo transport of the complete radiation spectra based on the unabridged nuclear decay data (RADTABS software, version 2.2) (22) was used for all radionuclides except 90Nb and 117mSn, for which the condensed Auger electron + Coster–Kronig electron + internal conversion electron (AE+CK+IE) spectrum was used because the unabridged spectrum was not provided in the tabulations. In the case of 89Zr, 90Nb, and 99mTc, the full β spectra were included. The AE+CK+IE and β particle contributions to S values were determined separately. Taking the nucleus as the target, simulations were run assuming uniformly distributed activity in the nucleus (N←N), in the cytoplasm (N←Cy), or on the cell surface (N←CS). Cellular S values are the mean absorbed dose to the target region (T) per unit cumulated activity in the source region (S); S(T←S) according to the MIRD formalism isEq. 1
where is the number of electrons (yield) emitted per transition j with energy , and is the fraction of the source energy deposited in the target region T (mass ) from activity in the source region S. S values included in Supplemental Tables 2–13 were calculated from an event-by-event Monte Carlo simulation implementing Equation 1. The resulting S values except for the nuclides 89Zr, 90Nb, and 119Sb, which were not included in the MIRD monograph, were compared with those in the MIRD tabulations (14) derived from Equation 2. The S value in the continuous-slowing-down, straight-trajectory approximation is given byEq. 2
where is the geometric reduction factor (14), is a semiempiric electronic (collision) stopping power of an electron with initial energy after passing a distance r through the medium, and rCSDA is the range in the continuous-slowing-down approximation.
DPKs
For DPKs, a point isotropic radiation source was placed in an infinite liquid water medium and the mean absorbed dose from the emitted electrons (AE+CK+IEs and β particles) was scored in 1-nm-thick spheric shells around the decay site. Mean absorbed doses were tallied up to a radial distance of 30 μm from the point source. This corresponded to the radius of a sphere in which 100% of all emitted energy from the Auger electron and Coster–Kronig electron spectra was absorbed for 67Ga, 80mBr, 89Zr, 90Nb, 99mTc, 111In, 117mSn, 119Sb, 123I, 125I, and 90% was absorbed for 195mPt and 201Tl.
To compare the dose deposition of the radionuclides with that of 223Ra (22), the ratio of the DPKs calculated for spheres representing different DNA condensation states was determined. The energy deposited by the 5.77-MeV α particles was approximated by multiplying the mass electronic stopping power (800 MeVcm2/g) over the path length (i.e., the radius of the sphere), whereas energy deposition of the AE+CK+IE spectrum was determined by event-by-event simulation.
RESULTS
Monte Carlo–Calculated S Values
Cellular S values calculated from Monte Carlo simulations are provided in Supplemental Tables 2–13. The statistical uncertainty of the S values was less than 1.3% with 95% probability (type A uncertainty). The contribution of photons is much smaller than that of electrons (<2% of total S values) and was disregarded in the present S value tabulations; this was also done to facilitate comparison with MIRD S values, for which the contribution from photon radiations was not included (Supplemental Fig. 1). For each radionuclide, the contributions to S values of the Auger electrons and Coster–Kronig electrons, as well as, when appropriate, the β spectra, is provided. In the case of 99mTc, the contribution of the β spectra to the cellular S values is negligible (on average <0.003% to S(N←N), as expected from the β− yield, and <0.03% to S(N←Cy) and S(N←CS)), and therefore the contribution of S(β spectra) to the total was omitted. S values were verified using the simulated DPKs, with differences of 2% or less (data not shown).
In summary, the S values calculated from PENELOPE Monte Carlo simulations were in excellent agreement with the MIRD data when the source and target regions strongly overlapped, that is, for N←N configurations. The differences between PENELOPE and MIRD for S(N←N) values were generally smaller than 10%, with MIRD values consistently smaller than PENELOPE values. The differences tended to decrease with increasing RN; however, discrepancies increase when the source is far from the target region, as is the case for N←Cy (≤30%) and N←CS (≤60%). This increased discrepancy can be ascribed to the respective calculation approaches; that is, the MIRD method propagates electrons in straight trajectories, not taking energy straggling into account. PENELOPE as a rule overestimates the dose contributions for all 3 source–target configurations except for 67Ga, 99mTc, and 201Tl, for which PENELOPE underestimates the N←Cy and N←CS contributions compared with MIRD S values.
Effect of Cellular Geometry on S Values
The effect of cellular geometry on the S values, taking into account only contributions from the AE+CK+IE spectrum, is summarized in Supplemental Tables 2–13. The self-dose to the nucleus is obviously not influenced by the position of the nucleus relative to the cell (concentric vs. eccentric nucleus arrangements) or the shape of the cell. However, the position of the nucleus, especially for eccentric arrangements, can contribute significantly to crossfire of neighboring cells, and this would specifically be seen for radionuclides with longer-range Auger electrons and internal conversion electrons. Generally, the contribution to S(N←Cy) is less for eccentric than concentric cell arrangements, with the greatest differences (≤30%) noted in smaller RC/RN configurations. These differences become less pronounced (<10%) when RC and RN increase. For eccentric versus concentric nuclear arrangements in which the activity is uniformly distributed on the cell surface, the greatest contribution to the mean absorbed dose delivered to the nucleus is again seen for small RN relative to RC (Supplemental Fig. 2, RC = 5 μm and RN = 2 μm). It follows that S(N←CS) contributions in eccentric nucleus configurations increase as the size of the nucleus decreases and thus the distance of the nucleus from the cell surface decreases in relation to the range of the particles. This is particularly apparent for 89Zr (Supplemental Fig. 4), for which up to a 60-fold increase in nuclear dose is observed for RC = 10 μm and RN = 5 μm. However, this effect decreases when the contribution from the β spectra is included (Fig. 2A).
To assess the effect of the nucleus position (in relation to the cell surface) on S(N←Cy) and S(N←CS), the eccentricity was evaluated for a cell with dimensions RC = 10 μm and RN = 5 μm (Fig. 2). Eccentricity is expressed as the ratio of S values for eccentric compared with concentric cell–nucleus arrangements. The nucleus position is shifted by 0.5-μm increments from the center until it is contiguous with the cell surface. In the case of 89Zr and 90Nb, the contribution of the β spectra is included. Nucleus eccentricity had the least effect on 80mBr, 117mSn, 123I, and 125I, with the contribution from cytoplasmic mean absorbed dose S(N←Cy) decreasing by less than 25% and the contribution from cell surface mean absorbed dose S(N←CS) increasing by less than 30%. For 89Zr, the β spectra contributes approximately 25% to the overall mean absorbed dose for S(N←Cy) configurations, whereas this contribution decreases from approximately 95% to 30% for S(N←CS) configurations as the nucleus shifts toward the cell surface. Conversely, the 90Nb β spectra contribution to S(N←Cy) and S(N←CS) configurations varied by approximately 30% and 30%–45%, respectively. The effect of eccentricity on dose calculations is discussed in Supplemental Figure 3.
Contribution of Auger Electron to S Values
Self-dose to the nucleus for most of the radionuclides studied is almost exclusively due to Auger electrons, with internal conversion electrons (and, when appropriate, β particles) responsible for the remaining mean absorbed doses (Supplemental Tables 2–13). Cytoplasmically distributed Auger electrons contribute less when RC increases than when RN increases. The contribution of Auger electrons to S(N←CS) can be roughly divided into 3 categories: the contribution increases when the distance between RC and RN decreases (67Ga, 80mBr, 99mTc, 195mPt, and 201Tl); the contribution increases when the distance between RC and RN increases (89Zr, 90Nb, and 119Sb); and the contribution of surface-bound Auger electrons to the nucleus is influenced by the end of the range of some electrons in the spectra (111In, 117mSn, 123I, and 125I). In the last category, 111In, 123I, and 125I are largely unaffected by RC and RN combinations, which makes these radionuclides more attractive for targeting strategies not relying on internalization of the construct.
DPKs
The DPKs for the radionuclides studied are presented in Figure 3. The DPKs calculated for 89Zr and 90Nb include the contributions of their β spectra (Supplemental Fig. 4). For all radionuclides considered, the energy deposition in the first 1-nm shell ranged from a minimum of 0.14 keV (67Ga) to a maximum of 1.21 keV (195mPt), with corresponding DPKs of 5.23 and 46.4 MGy, respectively. If a DNA targeting strategy is adopted, the radionuclides that will result in the highest DPK over a radial distance of 11 nm, chosen to represent chromatin structure (i.e., DNA double-helix wrapped around histones) are, in decreasing order of efficacy, 195mPt, 201Tl, 125I, 119Sb, 111In, 123I, 117mSn, 80mBr, 90Nb, 89Zr, 99mTc, and 67Ga. DPKs decrease dramatically by 8–9 orders of magnitude over the 30-μm distance considered, falling to mean absorbed doses of a few milligrays.
Relative Dose of Auger Electrons Compared with α Particles
A comparison of the energy deposition from α (223Ra) and Auger electron emitters (AE+CK+IE and β particles) in spheres of diameters representing different DNA condensation states is shown in Figure 4. The energy deposition in a sphere 2 nm in diameter, representing the DNA double-helix, exceeds that of a monoenergetic 5.77-MeV α particle for all the radionuclides apart from 99mTc and 67Ga (Fig. 4A). When the Auger electron emitters are compared with the complete particulate spectrum of 223Ra (5.77-MeV α, AE+CK+IE, and daughters), the effect is vastly reduced, with only 195mPt depositing a comparable amount of energy (Fig. 4B). In general, the relative dose of the Auger electron–emitting radionuclides diminishes dramatically as diameters representing DNA condensation states increase.
DISCUSSION
To evaluate the virtues of Auger electron–emitting radionuclides that have previously been suggested for therapeutic or imaging purposes, a single-cell dosimetric approach was undertaken by considering all energy deposition events and related probabilities. First, cellular S values were determined by Monte Carlo transport of all particulate radiation following the MIRD formalism, and second, DPKs were calculated for all the radionuclides considered in 1-nm radial bins up to 30 μm. In this regard, it is crucial to use a Monte Carlo code that provides an event-by-event simulation, as a loss of spatial resolution during particle transport from condensed simulation algorithms (i.e., grouping elastic, inelastic, and radiative events), and underestimation of secondary electrons, have a large effect on energy deposition (23). PENELOPE can perform event-by-event coupled photon–electron transport simulations, thus providing a more accurate estimation of the energy deposition than is possible with other general-purpose Monte Carlo codes (24,25).
PENELOPE cellular S values for overlapping source and target agree well with MIRD and previously published S values (26), with differences of less than 10% and 5%, respectively. However, larger discrepancies are seen when the source and target volumes are farther apart. The greatest difference (≤60%) between PENELOPE and MIRD S values was when activity was assumed to be distributed on the cell surface in accordance with the observations of Uusijärvi et al. (25). These differences can be ascribed to the fact that MIRD S values are based on approximate DPKs calculated from the collision stopping power (i.e., assuming straight electron trajectories and neglecting energy straggling). Differences might also be attributed to the different energy spectra used. Whereas the unabridged radiation spectrum from the MIRD monograph (22) was adopted in the Monte Carlo simulations, MIRD S values were generated from the energy spectra provided by Eckerman et al. (27). A limitation of the PENELOPE code when Auger electrons are transported is the set cutoff energy of 50 eV. Many of the radionuclides evaluated here have an abundance of low-energy Auger electrons, less than 50 eV, which thus are not transported. This could potentially lead to an overestimation of energy deposited within the first few nanometers of the DPKs calculated.
It was recently demonstrated that a small variation in mean absorbed dose can have a significant impact on tumor control probability (28). It is therefore crucial that the selection of a radionuclide-delivery vehicle addresses not only the distribution (particularly internalization) of the radionuclide but also the geometry of the targeted cells. As is shown, eccentric cell–nucleus arrangements can lead to increased S values for some of the radionuclides studied here. The greatest effect of nucleus eccentricity was noted for S(N←CS) configurations. For example, the dose to the nucleus from cell surface source arrangements for 201Tl increased by more than 200%. As is the case with the other radionuclides (67Ga, 89Zr, 90Nb, 99mTc, 111In, 119Sb, and 195mPt), this huge increase in S(N←CS) is due to an abundance of very low-energy electrons in their spectra. This property could be advantageous for cell types with eccentric nucleus arrangements for which targeting is achieved via surface-bound receptors, but it could also have a deleterious effect on surrounding normal cells in close contact with the targeted cell. Another point of consideration with eccentric cell–nucleus arrangements is cell size; for smaller cells and nuclei the effect of eccentricity is rather marked, with up to a 2-fold difference in dose.
The spatial dose gradients in the respective DPK spectra are the result of the ends of ranges of major low-energy electron groups. The sharp drop in dose after the first few nanometers illustrates the highly localized energy deposition caused by the Auger electrons, with local mean absorbed doses over this range in excess of 10 MGy. Radial dose distributions diverge by less than 25% compared with those previously reported for 111In and 125I (29,30) and might be ascribed to the different radiation spectra used.
Although β emitters are clinically widely used because of their long range, the recent success of the α-emitting radionuclide 223Ra (31) has focused attention on the use of shorter-range radionuclides for molecularly targeted radiotherapy. From the viewpoint of their cell-killing potential, the advantage of Auger electron–emitting radionuclides is their extremely short range and localized dose deposition. New strategies are being developed for delivery of Auger electron–emitting radionuclides to the nucleus, such as carbon nanotubes (33), gold nanoparticles (34), antibody-based peptides (35), and cell-penetrating peptides (36). Intranuclear delivery of Auger electron–emitting constructs results in relative biological effectiveness similar to that of α emitters but with a reduced crossfire effect compared with α emitters, making them more suitable for single-cell irradiation (9). The DPKs of 223Ra and the Auger electron emitters show that with respect to the energy deposited in spheres of DNA dimensions, only the higher-mass-number Auger electron emitters deposit comparable amounts of energy. However, comparison with a monoenergetic 5.77-MeV α particle notes the major advantage of Auger electron emitters at distances less than 11 nm, which is in agreement with data presented by Charlton (36).
CONCLUSION
Many Auger electron-emitters are suitable for theranostic applications, enabling simultaneous detection (PET and SPECT imaging) and treatment. The data presented here suggest that single-cell characteristics should be considered when one is designing these molecularly targeted agents. The challenge will be to match the delivery strategy with the physical properties of a particular radionuclide.
DISCLOSURE
The costs of publication of this article were defrayed in part by the payment of page charges. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734. This research was supported by Cancer Research U.K. (CR-UK), the Medical Research Council (MRC), the CR-UK Oxford Centre, and the Generalitat de Catalunya (project 2014 SGR 846). No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Jul. 23, 2015.
- © 2015 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication December 23, 2014.
- Accepted for publication July 8, 2015.