Visual Abstract
Abstract
In radiopharmaceutical therapy, intratumoral uptake of radioactivity usually leads to heterogeneous absorbed dose distribution. The likelihood of treatment success can be estimated with the tumor control probability (TCP), which requires accurate dosimetry, estimating the absorbed dose rate per unit activity to individual tumor cells. Methods: Xenograft cryosections of the prostate cancer cell line LNCaP treated with [177Lu]Lu-PSMA-617 were evaluated with digital autoradiography and stained with hematoxylin and eosin. The digital autoradiography images were used to define the source in a Monte Carlo simulation of the absorbed dose, and the stained sections were used to detect the position of cell nuclei to relate the intratumoral absorbed dose heterogeneity to the cell density. Simulations were performed for 225Ac, 177Lu, and 90Y. TCP was calculated to estimate the mean necessary injected activity for a high TCP. A hypothetical case of activity mainly taken up on the tumor borders was generated and used to simulate the absorbed dose. Results: The absorbed dose per decay to tumor cells was calculated from the staining and simulation results to avoid underestimating the tumor response from low absorbed doses in tumor regions with low cell density. The mean of necessary injected activity to reach a 90% TCP for 225Ac, 177Lu, and 90Y was found to be 18.3 kBq (range, 18–22 kBq), 24.3 MBq (range, 20–29 MBq), and 5.6 MBq (range, 5–6 MBq), respectively. Conclusion: To account for the heterogeneous absorbed dose generated from nonuniform intratumoral activity uptake, dosimetry models can estimate the mean necessary activity to reach a sufficient TCP for treatment response. This approach is necessary to accurately evaluate the efficacy of suggested radiopharmaceuticals for therapy.
- Monte Carlo dosimetry simulation
- radiopharmaceutical therapy
- digital autoradiography
- tumor control probability
- heterogeneity
Radiopharmaceutical therapy has become a promising approach to treating metastatic cancers, such as metastatic castrate-resistant prostate cancer (1,2). Radiopharmaceuticals that target cell-specific epitopes emitting short-range radiation with high linear energy transfer can deliver high absorbed doses to tumors while sparing healthy tissues. Nonuniform uptake of radiopharmaceuticals in a targeted volume can cause heterogeneous energy depositions, leading to large variations in the absorbed dose experienced by the cells.
The MIRD formalism assumes uniform activity in the source and calculates an average absorbed dose to the target volumes (3,4). Although the formalism can be applied on any scale, macroscopic or microscopic, it is commonly used with data from γ-cameras or SPECT imaging, where spatial resolution and sensitivity are limited and mainly organs can be delineated. Instead, digital autoradiography (DAR) can detect the intratumoral distribution of radioactivity (5–8). Chouin et al. used an α-camera to estimate the absorbed dose to cells in micrometastases after treatment with a radioimmunoconjugate labeled with 211At (9). Similar to the study presented here, they correlated the detected activity in cryosections to cells detected in adjacent sections stained with hematoxylin and eosin (HE).
Detailed dosimetry models in preclinical trials can help identify the most promising tracers for radiopharmaceutical therapy. Unexplained failure of tracers could be resolved by improving the tumor penetration of the tracer, optimizing for more homogeneous intratumoral distribution of the tracer, or changing the labeled radionuclide to one with longer-range emission. Our group previously improved the tumor uptake uniformity of 111In-DOTA-hu5A10 by increasing the chelate-to-antibody molar ratio in the labeling process, thereby improving the therapeutic effect in treated xenografted mice (10). Similarly, Howe et al. showed that a combination of carriers labeled with 225Ac with complementary intratumoral distributions generated improved radioactivity distribution and significantly reduced tumor growth compared with the same activity delivered by either of the 2 carriers alone (11).
The tumor control probability (TCP) estimates the probability of killing all cells in a lesion from absorbed dose and cell radiosensitivity data (12–14). The intratumoral radioactivity distribution affects TCP for short-range radiation (15). This paper aims to calculate TCP from dosimetry simulations of heterogeneous activity distributions measured with DAR.
MATERIALS AND METHODS
Sections from xenografts treated with [177Lu]Lu-PSMA-617 from BALB/cAnNRj mice were used to build a dosimetry model. The activity detected in the DAR image pixels was used to define voxels in the source volume in a Monte Carlo simulation of the absorbed dose. By matching DAR and HE images, the cells segmented from the HE stain within an aligned DAR image pixel were assumed to receive the absorbed dose simulated to the corresponding target voxel. TCP, as a function of injected activity, was calculated considering the cell’s simulated absorbed dose per decay (from now on called dose values or dose image). Detailed descriptions of radiolabeling, cell culturing, animal work, and autoradiography can be found in supplemental materials (supplemental materials are available at http://jnm.snmjournals.org).
Monte Carlo dose simulations of 177Lu, 90Y, and 225Ac were performed in GATE version 8.1 (OpenGATE), which in turn uses Geant4 version 10.3.3 (European Organization for Nuclear Research) (16). All image and data processing were performed in MATLAB R2020b (MathWorks). Each tumor’s 8 DAR images were decay-corrected to the same time point and aligned (coregistered), an average of the images was calculated, and grayscale values were normalized. The tumor borders were detected by thresholding the grayscale values. The resulting mask was used to coregister the DAR and HE images. Cell nuclei were segmented from the HE images (details in supplemental materials), and a cell density map was generated in which the number of cells within the corresponding pixel in the coregistered DAR image was calculated. A hypothetical case in which the activity is primarily taken up at the edge of the xenograft because of low tumor penetration was generated by filtering the mean DAR images (description and filter in Supplemental Fig. 1).
GATE Monte Carlo Simulations
The β-decay of both 177Lu and 90Y has a yield of 100%, and the primary β-emissions were included in the simulations. 177Lu emits β-particles with maximum (endpoint) β energy (Eβ(max)) of 497 keV (79.4%), 384 keV (8.9%), 247 keV (0.016%), and 176 keV (11.7%) (17). 90Y emits β-particles with Eβ(max) of 2,279 keV (99.9%), 518 keV (0.01%), and 92 keV (0.0000014%) (17). The mean β-energy spectrums of these emissions considered for 177Lu and 90Y, shown in Supplemental Figure 2A, were collected from the International Atomic Energy Agency’s Livechart of Nuclides (17)—originally computed by Betashape (18,19)—and defined in GATE as histogram sources. The histogram had a step width of 0.5 keV for 177Lu and 0.3 keV for 90Y. No other electron or γ-emissions were considered because it was assumed their contribution to the final absorbed dose was negligible.
For 225Ac simulations, the highest-energy α-emission of 5.8 MeV was considered, together with the highest-energy α-emissions from the daughters 221Fr (6.4 MeV), 217At (7.1 MeV), and 213Po (8.4 MeV), thereby simplifying and speeding up the simulation. The decay chain of 225Ac is shown in Supplemental Figure 2B. The daughters were assumed to remain immobile relative to the mother nuclide. To generate a realistic distribution of these emissions, a discrete spectrum was defined in the GATE simulation, where the 4 α-energies all have the same probability of being emitted as primary particles. Again, no γ-particles or electrons were included in this simulation because they were assumed to be negligible.
In the simulation, the world was defined as a water cube. To build a source volume, each pixel of the mean DAR image was defined as a box-shaped, general-purpose source placed in the water volume as a grid of voxels. The energy spectrum of the simulated radionuclide was ascribed to the general-purpose sources. The intensity of each source voxel was set to the normalized grayscale pixel values of the mean DAR image, thereby giving the relative probability for the emission of a simulated particle. Their x- and y-dimensions remained 50 × 50 μm, like the DAR image pixel size. The z-dimension was set to twice the continuous slowing-down approximation range in liquid water of the maximum energy of the emission of the radionuclide, except for 90Y. Instead, the range was limited to 10 mm, because no axis length of the included tumors was longer. For 225Ac, the daughter 213Po had the highest α-energy considered in the simulation, with an estimated continuous slowing-down approximation range of 0.085 mm (20). For 177Lu, the Eβ(max) 496.8 MeV equated to a continuous slowing-down approximation range of 1.8 mm (20). The electromagnetic interactions of the primary and secondary emissions were simulated by the Geant4 physics list emstandard_opt3 (16).
The energy of simulated primary particles was sampled from the energy spectrum, and its origin was sampled from the voxel intensities. Simulations were run for 109 primary particle emissions. A GATE dose actor was defined that overlapped the source voxel volume with the same resolution as the source in the x- and y-direction (50 × 50 μm2) and 10 μm in the z-direction, mimicking the tumor cryosection thickness. To convert the absorbed dose recorded in the target volume to the absorbed dose per decay (Gy Bq−1 s−1), in the ith voxel was calculated as the total number of nuclear transformations (decays), and the simulated absorbed dose was multiplied with the yield Y of the emissions and the number of primary particle emissions P in the simulation (Equation 1): Eq. 1
For 177Lu and 90Y, the yield Y equals 100%. For 225Ac, the α-decay has a 100% yield. However, because of the included daughter α-emissions, 1 in 4 of the simulated primary particles originate from 225Ac; therefore, the 225Ac yield in the simulation can be described as 400%.
The cell density map was used to overlap the resulting simulated dose image. Cells detected in the corresponding target voxel were assumed to receive its absorbed dose per decay. Histograms of equal binning of target voxels and cell dose values were used to calculate the cumulated dose–volume histograms.
TCP, as a function of injected activity, was calculated from the cell dose values. For each dose value interval j, the absorbed dose was calculated. It is defined by the MIRD formalism (3) as the product of the cumulated activity and the dose value: Eq. 2
where T1/2 is the physical half-life.
The initial activity in the source volume was calculated as the product of the uptake U (percentage injected activity per gram), assuming instantaneous uptake; the injected activity (Bq); and the tumor source volume considered in the simulation (g): Eq. 3
To exemplify a realistic case, average activity uptake was assumed to be 3.6% per injected activity per gram of xenograft tissue, as previously measured by our group (22), with a tissue density of 1.0 g/cm3.
The cumulated activity is the total number of decays from the fraction of the injected activity taken up by the tumor. To exemplify the use of DAR images for TCP calculations, we simplified the model. We assumed no biologic clearance and no redistribution within the tumor. All activity measured 3 d after injection was therefore assumed to have been taken up instantaneously without redistribution during the dose integration period. For a more realistic model, uptake should be measured at several time points to estimate the cumulated activity better. For each radionuclide, a range of injected activities was evaluated within relevant intervals; 0–100 MBq for 177Lu, 0–30 MBq for 90Y, and 0–100 kBq for 225Ac.
The survival probability and TCP were calculated as previously described by Nahum (13) and Bernhardt et al. (12) and later summarized by Uusijärvi et al. (14). For each dose value interval j, the survival probability for the cells in that interval that were assumed to be identical was calculated as in Equation 1. A range of values for the radiosensitivity α in a relevant interval based on radiosensitivity measurements performed by Elgqvist et al. (21) was applied in the calculations to evaluate its effect on the resulting TCP: Eq. 4
Equation 4 is a simplification of the linear-quadratic model (13) and is applicable for radiation with high linear energy transfer, such as α-particles or β-particles delivering an absorbed dose of more than approximately 2 Gy.
TCP, defined as the probability to kill all cells, for all intervals j—normalized for the number of cells in each interval receiving the absorbed dose —is then given by the following equation: Eq. 5
RESULTS
The results for tumor 3 are presented here. Results for tumors 1 and 2 are found in Supplemental Figures 3–25.
The mean of 8 aligned DAR images from tumor 3 is shown in Figure 1A. The modified DAR image, shown in Figure 1B, represents a case of reduced tumor penetration. Activity has been moved from the central parts of the tumor and concentrated on the edges. The cell density map of tumor 3 in Figure 1C found between 0 and 30 cells per pixel. Segmentation results are presented in Supplemental Table 1.
The resulting dose values in the dose actor voxels in tumor 3 are presented in Figure 2. The differing particle range of the 3 radionuclide emissions can be seen in the gradually smoother distribution as the range increases. The long range of the 90Y β-particles generated a smooth dose image in which local variations in activity uptake were indistinguishable. In contrast, the α-particle emissions of 225Ac and its daughters caused a more heterogeneous distribution with hot spots induced by local activity clusters. Effectively, cells residing in a voxel with low activity uptake could still receive a relatively high dose value when simulations were run for 90Y but depended more on local uptake when 225Ac was simulated. Similar results were seen for tumors 1 and 2 (Supplemental Fig. 4).
For simulations performed with the modified activity distributions, the maximum dose values for all radionuclides increased, as did the focus of higher dose values to voxels with higher source intensity. Although 90Y generated the most homogeneous dose value distribution, differences between edges and central parts increased (Supplemental Fig. 5).
The cell dose value distribution, calculated by matching the dose value image to the cell density map, for 177Lu in tumor 3 is presented in the histogram in Figure 3A. Although the mean absorbed dose rate per unit activity in cells was 2.2 × 10−10 Gy Bq−1 s−1, it ranged from close to 0 to more than 4 × 10−10 Gy Bq−1 s−1. This can be compared with the target voxel dose value distribution shown in Figure 3B, where the mean was 1.2 × 10−10 Gy Bq−1 s−1. The shape of the histograms differs greatly, because the cell density varies over the tumor section. Many voxels experienced a low dose value, but these contained few cells, whereas the cell dose value distribution was approximately centered on its mean. The cumulated absorbed dose rate histograms for cells and voxels are plotted in Figure 3C.
The resulting cell dose value histograms of all tumors and radionuclides are shown in Supplemental Figure 15, and the modified activity distributions appear in Supplemental Figure 16. Comparing 90Y and 177Lu, the mean absorbed dose rate per unit activity was higher for 90Y. However, between the two, the highest dose value was received by 177Lu in tumor 2.
From the cell dose value histograms, TCP was calculated for all 3 tumors as a function of injected activity for ranges of activity realistic to inject in a mouse model. Elgqvist et al. investigated the radiosensitivity of several prostate cancer cell lines, including LNCaP (21). By fitting their data to Equation 4, LNCaP should have a radiosensitivity of 1.33 Gy−1 (95% CI, 0.96–1.70 Gy−1) for α-particles and 0.21 Gy−1 (95% CI, 0.16–0.25 Gy−1) for β-particles. Based on this, TCP was calculated as in Equation 5 for intervals of radiosensitivity of 0.5, 1.0, 1.3, and 2.0 Gy−1 for 225Ac and 0.1, 0.2, 0.3, and 0.4 Gy−1 for 177Lu and 90Y. The resulting TCPs for tumor 3 are shown in Figure 4. Corresponding figures for tumors 1 and 2 can be found in Supplemental Figures 17–25. The necessary activities to be injected to reach a 90% TCP are summarized in Table 1, and the modified activity distributions are summarized in Table 2.
When treating tumor 3 with 225Ac, as shown in Figure 4A, an injection of about 20–30 kBq should be expected to give a good tumor response. However, when the activity remained at the border of the tumor, as shown in Figure 1B, the injected activity needed to reach a 90% TCP greatly increased, as seen in Table 2. This TCP level was never reached for a radiosensitivity of 0.5 Gy−1; however, this was a low estimate for radiosensitivity. The necessary activities to be injected might be so high that they cause damage to healthy tissues.
For treatment of tumor 3 with 177Lu, as shown in Figure 4B, an injection of less than 30 MBq would unlikely be curative. The necessary activity for a high TCP increased when the activity distribution in the tumor was modified, as summarized in Table 2. For a radiosensitivity of 0.1 Gy−1, no activity below 100 MBq was sufficient to reach a 90% TCP.
Finally, TCP for tumor 3 treated with 90Y is shown in Figure 4C. The relative difference in activity needed for a 90% TCP when comparing the original and the modified activity distribution is smaller for 90Y than for 177Lu and 225Ac because of its longer β-particle range.
For comparison, if TCP was instead calculated from the voxel dose values, none of the investigated injected activities or radiosensitivities reached a TCP of at least 90% for 225Ac or 177Lu in tumor 3 (results in Supplemental Figs. 17, 20, and 23).
DISCUSSION
A common approach for tumor dosimetry in radiopharmaceutical therapy is to assume a sphere or ellipsoid of evenly distributed tumor cells with homogeneous radioactivity uptake. For short-range radiation, this oversimplification is inappropriate (23), because it risks miscalculating the absorbed dose to cells if activity uptake is heterogeneous. This is accounted for by performing voxel dosimetry, which still ignores the cellular distribution.
In preclinical trials, dosimetry models are necessary to evaluate which radiotracers have the potential to generate good treatment responses. However, overly simplified dosimetry models might mislead researchers instead of guiding their decision-making. An improvement considers intratumoral activity uptake and its relation to tumor cell distribution. Then, based on TCP calculations, realistic activities to be injected for optimal treatment effect can be estimated.
This study simulates heterogeneous absorbed dose distributions within xenografts treated with PSMA-617–ligated radioactivity. To improve calculations of TCP, we connect dose values to the number of cells experiencing them. This generates a dose value distribution that is different from the distribution generated when only voxels are considered, as seen in Figures 3A and 3B. This way, the treatment response to a wasted dose—that is, energy deposited in volumes where few cells reside—and the response of cells in volumes receiving less than necessary for tumor control will not be overestimated.
We estimate the minimum injected activities necessary to reach a 90% TCP for varying radiosensitivities. For 225Ac, assuming a radiosensitivity of 1.3 Gy−1, the mean injected activity would be 18.3 kBq (range, 18–22 kBq). For 177Lu and 90Y, assuming a radiosensitivity of 0.2 Gy−1, a mean of 24.3 MBq (range, 20–29 MBq) and 5.6 MBq (range, 5–6 MBq), respectively, would be necessary. In the case of 177Lu, these numbers agree with the injected activities that resulted in good tumor response when previously investigated by our group (22). The sample consists of only 3 tumors, so the resulting numbers are uncertain.
Along the slicing axis of the sectioned xenograft, the activity distribution and cell density change only slightly. This is an argument for the validity of our approach to evaluating the TCP of the tumor, because we assume the TCP of a section of the tumor is representative of the whole tumor.
We present here a general model; more parameters should be included to increase the accuracy. In our model, uptake is based on a single time point, and there is a lack of pharmacokinetics because we have included only the physical clearance of the radionuclides. Some biologic clearance is expected, so the results likely overestimate TCP. By measuring at several time points, one might improve the estimate of cumulated activity. In addition, no consideration of DNA damage repair, tumor repopulation, or differences in dose rate, which are relevant (24), is made.
Regions of varying cell density and cell type, necrotic areas, vascular structures, etc., are seen throughout the tumor volumes. In addition, the targeted epitopes of the tumor cells might not be equally available because of restricted tumor penetration or varying expression. Bordes et al. (25) used the measured uptake of rituximab by fluorescence microscopy in a multicellular aggregate of lymphoma cells, distinctly limited to the edges, to represent activity uptake in Monte Carlo simulations of the absorbed dose rate per activity unit in multicellular volumes (25). This method is favored by the greater resolution of fluorescence microscopy. However, it is not a direct measurement of actual activity uptake, because uptake might differ when labeling a tracer molecule to a radionuclide rather than a fluorophore.
For short-range radiation with high linear energy transfer, the microscopic energy deposition distribution can affect the absorbed dose to the cell nucleus. Cellular internalization can shorten the distance between decay and nucleus, thereby increasing energy deposited where it is most effective. The DAR images’ pixel size limits the spatial resolution. Similar to the range of α-particles emitted in the 225Ac decay chain, no microscale heterogeneity will be considered. However, Mínguez Gabiña et al. simulated TCP in a cluster of cells with varying 225Ac-PSMA internal uptake and only found a small difference between activity on the cell surface and activity inside the cytoplasm (26).
CONCLUSION
We have shown how to improve preclinical dosimetry in radiopharmaceutical therapy by considering intratumoral activity uptake and its relation to tumor cell distributions. Realistic activities to be injected for optimal treatment effect can be determined with Monte Carlo simulations and TCP calculations. Examples are given for LNCaP xenografts treated with radiolabeled PSMA-617. Our approach, which considers the intratumoral distribution of cells rather than only the voxel volume, avoids underestimating the mean experienced absorbed dose rate per unit activity, because the influence of the wasted dose on the dose value calculations is reduced.
DISCLOSURE
This study was performed with support from the Swedish Cancer Society and Mrs. Berta Kamprad’s Foundation. No other potential conflict of interest relevant to this article was reported.
KEY POINTS
QUESTION: Can combining DAR images with HE-stained xenograft sections in a Monte Carlo dosimetry model improve calculations of TCP?
PERTINENT FINDINGS: The model finds dose values experienced by cells in the tumor. We calculated TCP and estimated the necessary injected activity for LNCaP xenografts treated with PSMA-617 radiolabeled to 225Ac, 177Lu, and 90Y.
IMPLICATIONS FOR PATIENT CARE: Improved dosimetry models are vital to evaluate radiotracers’ potential in a preclinical phase.
ACKNOWLEDGMENTS
We thank Wahed Zedan for cell culturing and animal handling and Anders Örbom for help with DAR imaging.
Footnotes
Published online Jul. 27, 2023.
- © 2023 by the Society of Nuclear Medicine and Molecular Imaging.
Immediate Open Access: Creative Commons Attribution 4.0 International License (CC BY) allows users to share and adapt with attribution, excluding materials credited to previous publications. License: https://creativecommons.org/licenses/by/4.0/. Details: http://jnm.snmjournals.org/site/misc/permission.xhtml.
REFERENCES
- Received for publication January 30, 2023.
- Revision received June 12, 2023.