Abstract
123I-Labeled iodoazomycin arabinoside (IAZA) is a marker of hypoxia in vivo. It has been used clinically to image hypoxic tissue in solid tumors, peripheral vascular disease of diabetic origin, blunt brain trauma, and rheumatoid joints and in an animal model of cerebrovascular disease. The radiation dose biodistribution for 123I-IAZA was studied to assess and characterize its suitability as a clinical radiopharmaceutical. Methods: Six healthy volunteers each received a nominal 185-MBq (5 mCi) dose of 123I-IAZA administered as a slow (1–3 min) intravenous injection in the arm. Anterior and posterior whole-body planar images were acquired for each volunteer beginning immediately after injection and at 1–2, 3–4, 6–8, and 20–24 h after injection. Venous blood samples (0 h predose through 28 h after dosing) and 28-h cumulative urine samples were taken from each volunteer for pharmacokinetic analysis. Radiation dose estimates were performed for all volunteers, with “reference adult” (for men) and “adult female” (for women) phantoms, and both the International Commission on Radiological Protection 30 gastrointestinal tract model and the dynamic bladder model, using the MIRDOSE3 program. Two sets of estimates, 1 using a pharmacokinetic analysis of total serum radioactivity and 1 based on scintigraphic image data, were obtained for each volunteer after 123I-IAZA administration. Results: Two compartments were discernible by pharmacokinetic analysis, and 4 compartments were discernible by image analysis. The urinary bladder wall received the greatest radiation dose (6.3E-02 ± 8.7E-03 mGy/MBq), followed by the upper large intestinal wall (5.6E-02 ± 1.2E-02 mGy/MBq), the lower large intestinal wall (5.0E-02 ± 1.2E-02 mGy/MBq), and the thyroid (4.4E-02 ± 1.4E-02 mGy/MBq). Approximately 90% of physiologically eliminated radioactivity was excreted through the kidneys. Radioactivity entering the intestinal tract from the gallbladder constituted <10% of biologically eliminated activity. Conclusion: The dosimetric analysis of 123I-IAZA in 6 healthy volunteers indicated that both disposition kinetics and radiation dosimetry support its clinical use for imaging tissue hypoxia.
Although originally developed for imaging regional hypoxia in patients with cancer (1–5), 123I-iodoazomycin arabinoside (IAZA) has also been used clinically to image hypoxic tissue in peripheral vascular disease of diabetic origin (6), blunt brain trauma (7), and rheumatoid joints (8) and in an animal model of cerebrovascular disease (9,10).
Radiation dosimetry calculations used to predict the radiation dose for patients receiving the radiopharmaceutical take into account radiation type and half-life, the amount of radioactivity administered, and tracer kinetics, including source organ identification and their respective biologic half-lives. Target organ mass, shape, and location and the intraorgan absorbed fraction of energy emitted by all source organs are usually estimated from scintigraphic images. Organ-specific parameters (S values) for a variety of human phantoms and radioisotopes have been published (11) and are conveniently available in the MIRDOSE computer program (Oak Ridge Institute of Science and Education, Oak Ridge, TN). Stypinski et al. (12) reported the clinical pharmacokinetics of IAZA, the radiopharmacokinetics of 123I-IAZA, and total radioactivity kinetics for 6 healthy volunteers. The 123I-IAZA radiation dosimetry estimates for those 6 volunteers are reported in this article, and this dosimetric information is compared with the previously reported pharmacokinetic data.
Materials and Methods
Chemicals and Reagents
IAZA was synthesized and radiolabeled using published procedures established in our laboratory (1,13).
Clinical and Imaging Protocols
The research was conducted in accordance with the tenets of the Declaration of Helsinki (1964) and was approved by the Alberta Cancer Board Research Ethics Committee and the University of Alberta Radiation Safety Committee. Details of the clinical protocol have been described elsewhere (12). Six healthy volunteers (4 men, 2 women; age range, 26–54 y; mean age, 37 ± 13 y) participated with informed consent. Lugol’s solution (0.6 mL, U.S. Pharmacopeia) in orange juice was given orally to block uptake of radioiodide by the thyroid, and shortly thereafter, each subject received a nominal 185-MBq (5 mCi) dose of 123I-IAZA administered as a slow (1–3 min) intravenous injection in the arm.
Anterior and posterior whole-body planar images were acquired at 5 different time periods beginning immediately after injection and at 1–2, 3–4, 6–8, and 20–24 h after injection. All image times were 30 min. The images were acquired using a dual-head, large- field-of-view gamma camera (Odyssey 2000; Picker International Canada, Inc., Brampton, Ontario, Canada) equipped with a low-energy, all-purpose collimator and interfaced to an Odyssey computer (Picker International Canada). The 20% analysis window was set symmetrically over the 159-keV 123I photopeak.
Fourteen venous blood samples, from 0 h (predose) through 28 h after dosing, and 28-h cumulative urine samples were taken for pharmacokinetic analysis from each volunteer (12).
Planar scintigraphic images were exported to a personal computer equipped with MEDisplay (version 1.32s; C-Soft, Edmonton, Alberta, Canada), a medical image analysis package. The number of counts and pixels in selected regions of interest (ROIs) were obtained for each volunteer: liver, left kidney, large intestine, urinary bladder, thyroid, gallbladder, and the whole body (WB). WB counts were determined by drawing a contour around the entire body image. An off-body region drawn to the right of the volunteer’s head was used as instrumental background. To partially compensate for source organ overlap, several overlapping ROIs were included: a region of overlap between the right kidney and liver, a region of left kidney with large intestine overlap, and a region of left kidney where there was no large intestine overlap. The right kidney and liver overlap ROI was established by drawing 2 parallel lines through the patient’s body: 1 line passing through the uppermost point of the left kidney and the other passing through the lowest point of the liver. The left kidney region between the 2 lines was drawn, flipped around its horizontal axis, and dragged to the anatomic position of the right kidney, between the 2 lines. The count contribution from the blood pool to the organ activity was included as part of the organ activity for all ROIs.
Dosimetry Protocols
Because SPECT or CT images were not available for the volunteers, assigning activity to organs with overlapping regions as viewed from a planar 2-dimensional image was based on assumptions that were dependent on image time. The analyses initially were performed manually, and once the protocol was established all calculations were reanalyzed using Excel (version 7.0; Microsoft, Redmond, WA) for Windows 95 (Microsoft). The detailed protocol used to estimate source region activity and overlapping activity distributions is available.
Once the posterior and anterior counts in each organ were determined, the counts were further corrected by subtracting instrumental background activity (as counts per pixel) from the counts-per-pixel organ activity. The geometric mean counts for each organ were then calculated. No linear attenuation correction was made for the anterior and posterior tissues (i.e., bone, muscle, fat, and skin) above or below the organ of interest.
To determine the fractions of radioactivity entering the gut and the bladder, the dosimetry equations include WB, total body (TB), and all body (AB). TB and AB counts were determined from the anterior and posterior images at each acquisition time: Eq. 1
Gastrointestinal tract counts consisted of the sum of the geometric mean of counts in the large intestine and the gallbladder. AB counts were determined according to the formula: Eq. 2
It is important to state at the outset that the pharmacokinetic mean residence time (MRT) and the radiation dosimetry residence time (RT) have somewhat different meanings (14). The RT is the area under a source organ’s time–activity curve divided by the administered activity, whereas the pharmacokinetic MRT may be defined as the area under the concentration–time moment curve divided by the area under the concentration–time curve.
Once the geometric mean counts in an organ were determined for all 5 imaging times, the MRT of the radioactivity in the organ was determined using LAGRAN (version 1.0D; C. Ediss, University of Alberta, Edmonton, Alberta, Canada) (15), by area under the curve (AUC) integration of the ROI count versus time plot. The RT was calculated for the large intestine, liver, kidney, thyroid, WB, AB, and TB. Because each image acquisition took 30 min, it was assumed that the mean count in an organ corresponded to the midtime of 15 min after the start of image acquisition. The only exception was the thyroid, for which it was assumed that the mean counts were acquired at 5 min from the image start.
The 3 regions of interest, WB, TB, and AB, all include the counts from the liver, kidneys, and thyroid. To account for the activity in the body that is not in these organs and has not been excreted into the bladder or the gastrointestinal tract, the remainder of the body (RB) was included as an additional source organ. The RTfr of this source organ was calculated as: Eq. 3
Renal and Gastrointestinal Tract Excretion Fractions
The fraction of the radioactive dose eliminated through the hepatobiliary system into the intestinal tract was determined from the gamma camera images by first calculating the rate of elimination of radioactivity into the gut (Kgut), according to the formula: Eq. 4
Thus, the fraction of radioactive dose eliminated by the gut was determined as: Eq. 5
Because the only other biologic route of excretion from the body available for a radiopharmaceutical is the kidneys, the fraction entering the bladder was calculated as: Eq. 6
Radiation dose estimates were performed for all volunteers using the MIRDOSE3.1 program. The phantoms used were “reference adult” for the male volunteers and “adult female” for the female volunteers. The model used for analysis included a combination of both the International Commission on Radiological Protection 30 gastrointestinal tract model and the dynamic bladder model (16). Two sets of estimates were performed for each volunteer, one using a 2-h and one using a 4.8-h bladder void time.
Results
Immediate and 22-h anterior and posterior gamma camera images of volunteer 3 are shown in Figure 1. Typically, on the earliest images, the organs with the highest radioactivity accumulation were the bladder, liver, and kidneys. The heart was visible in early images, but not at longer times; therefore, the heart activity was attributed to blood pool. The intravenous injection site also was visible on the immediate images of this volunteer, but the percentage of radioactive dose in that area, as determined by image analysis, was negligible and therefore not expected to alter the distribution kinetics of the radiopharmaceutical. The injection site was also visible, to a lesser extent, on the immediate images of volunteer 5 (not shown).
There was a striking absence of radioactivity in the brains of all the volunteers in early- and intermediate-time images, indicating effective exclusion of 123I-IAZA by the blood–brain barrier. Later images, however, showed some redistribution of radioactivity into the brain. Immediate, 1- to 2-h, and 3- to 4-h images also did not show thyroid or large intestine uptake; in fact, these organs did not become readily visible until the 20- to 24-h image acquisition time (Fig. 1).
During the determination of gastrointestinal tract counts, only the large intestine and gallbladder counts were included because the small intestine was never clearly distinguishable in any of the images. Because the radioactivity entering the intestinal tract from the gallbladder constituted a small fraction of the dose (<10% of biologically eliminated activity), it was only visible when either concentrated in a small area (i.e., the gallbladder) or 20 h after dose administration, after a substantial portion of the dose was eliminated from the rest of the body. Some concentration of intestinal activity was evident in the 6- to 8-h image, but at that time it was either localized in the junction of the ileum and the ascending colon or already in the initial part of the ascending colon. In each case, this activity was included as part of the large intestine ROI.
A summary of radiation dose estimates for all volunteers for the 2-h bladder void time is shown in Table 1. The source organs were liver, kidneys, thyroid, and RB, and the activities in these organs were assumed to be uniformly distributed. Radiation dose caused by activity observed in the urinary bladder and the gastrointestinal tract was calculated separately. Although the reference adult phantom used for dose estimates in men does predict radiation dose to the ovaries as well as testis, the adult female phantom does not predict the dose to the testis. Volunteers 4 and 6 were women.
The organs receiving the greatest radiation dose were the thyroid, upper and lower large intestinal wall, urinary bladder wall, and small intestine. The bladder was an important source organ, because >90% of the physiologically eliminated radioactivity was excreted by the kidneys. Encouraging the patients to urinate frequently, especially during the first 8 h after dose administration, would substantially reduce the radiation dose to the urinary bladder wall.
Discussion
The primary source organ for the majority of target organs listed in Table 1 was the RB. Exceptions were the lower large intestinal wall, upper large intestinal wall, small intestine, kidneys, liver, thyroid, and urinary bladder wall, all of which were their own primary sources and targets. For the urinary bladder and upper large intestinal and lower large intestinal walls, it was assumed that the source of activity was the organ content and not the organ wall. The primary contributors for the uterus and TB varied with the bladder void time, with the RB and the urinary bladder acting as the primary contributors for a 2-h void time. Similarly, for the effective dose (16), the primary source organ was the colon or the urinary bladder.
The main sources of error involved in the dose estimates presented in Table 1 include the lack of attenuation correction and removal of overlying background from the ROI source organ activity determinations. The failure to remove overlying background activity from the organ ROI will overestimate the radiation dose; however, this will be partially compensated for by the lack of attenuation correction, which will result in underestimating the dose. The assumption of uniform distribution of radioactivity in an organ is another important factor when using the MIRD schema. This is especially important for organs such as the large intestine, where gastrointestinal tract transit time determines radioactivity localization at the time of image acquisition, and for the kidneys, where the highest activity concentration after 123I-IAZA administration was visible in the renal pelvis area. Not only will different parts of these organs receive variable radiation doses, but also any overlap ROI calculation involving these organs will introduce errors in the ROI counts determination of all the overlapping organs. A final important source of error was the accuracy of ROI drawing, which was estimated at ±10% based on 5 repeated drawings of the same ROI.
The radioactivity in the thyroid gland was thought to be caused by metabolic deiodination of 123I-IAZA, with subsequent active accumulation of 123I-iodide. For the purpose of this study, the radioactivity in the thyroid gland was considered irreversibly bound and was eliminated with the 13.2-h decay half-life of the isotope. The volunteers in this study received a single oral dose of Lugol’s solution at the time of 123I-IAZA injection. Because iodine incorporation into the thyroid gland takes place with a 6- to 8-h half-life, it may be possible to further decrease the radiation dose to the thyroid by administering Lugol’s solution earlier in the protocol (17,18). Although a single dose did not completely block the thyroid, in volunteer 6 the thyroid gland was virtually undetectable even at the time of the 20- to 24-h image acquisition.
Table 2 compares 123I-IAZA and several radiopharmaceuticals that are either commercially available or are being investigated as potential clinically useful imaging agents. Dose estimates are for adults with a 2-h bladder void time. These data emphasize the low doses from 123I-IAZA, because the 123I-IAZA radiation dosimetry is comparable with some 99mTc-labeled agents. 99mTc-sestamibi, for example, is commonly given as a 1.1-GBq injection, resulting in a TB dose of 5 mGy, whereas 123I-IAZA, clinically given as a 370-MBq dose, will produce a TB dose of approximately 3 mGy.
The TB RT calculated using gamma camera ROI information and the pharmacokinetic MRT, which was calculated as the area under the concentration–time moment curve divided by the area under the concentration–time curve, are not identical in concept. Individual organ RTs represent the distribution of activity in various organs in the body as identified by ROI analysis, and the RB represents the tissues and organs with similar distribution characteristics; the whole-body RT of the tracer is the sum of the RTs of each of these components. In multicompartmental pharmacokinetic analysis of serum radioactivity, the organs representing the RB would all constitute a single compartment, whereas organs with different MRTs would each constitute a separate compartment. To compare these approaches, it is important to use pharmacokinetic data representing total radioactivity rather than the radiopharmaceutical alone because no differentiation is possible with gamma camera analysis. After the serum-based pharmacokinetic analysis of total radioactivity after 123I-IAZA administration in these volunteers, only 2 compartments were discernible, as can be seen from the plot of total radioactivity in Figure 2 (12). However, gamma camera scintigraphic images of the WB appear to show at least 4 regions with elevated uptake: the liver, kidney, thyroid, and RB. An attempt was made to discern whether these regions were separate compartments of a multicompartmental pharmacokinetic model.
An explanation for this apparent discrepancy between the serum-based and scintigraphic-based models can be obtained through inspection of the different sources of data sampled by each technique. We observed that the scintigraphic-based RTs calculated for the liver and kidney were 7.35 ± 1.05 h and 7.53 ± 1.47 h, respectively. A 1-way ANOVA (power of significance = 0.05) did not detect statistically significant differences in the organ RTs and, hence, both the liver and the kidney could be ascribed to the same pharmacokinetic compartment (Table 3). Because the thyroid represents <0.5% of the AB cumulated activity, it would not be discernable as a separate compartment by serum analysis. TB imaging data represent activity mainly caused by blood and are therefore expected to be quantitatively similar to the serum-based model.
Classic pharmacokinetic sample acquisition methods generally rely on venous blood samples. Consequently, only quantitatively important compartments with different MRTs can be identified. In the case of 123I-IAZA, the very small fraction of dose accumulated in the thyroid would not be detected by serum analysis, despite the long MRT of activity in this organ. Furthermore, serum analysis would group all organs with similar MRTs into a single compartment, regardless of their respective concentrations of radioactivity. This is in keeping with the tenet that if the rate of change of activity of different organs (as represented by ROIs) is similar, then they can be attributed to the same compartment (22).
The distribution phase of total radioactivity, as determined from the pharmacokinetic analysis of blood samples, was found to represent only 5% of the overall AUC (12). With a 30-min image acquisition time, this distribution phase was too short to be discernable with gamma camera scintigraphic techniques. It was not possible to decrease the whole-body acquisition time caused by gamma camera sensitivity and the amount of 123I-IAZA injected. Thus, the source organs were the sources of radioactivity not by virtue of longer RTs, but because they had relatively higher affinity for the radioactive components, which resulted in greater volumes of distribution of total radioactivity.
In addition to the estimate of the TB RT and MRT, the 2 methods also can be used independently to determine the biologic elimination routes of total radioactivity from the body. In pharmacokinetic analysis, the fraction eliminated in urine was determined directly by analyzing the pooled 0- to 28-h urine samples. The remainder of the dose was assumed to be excreted through the gut (12). Scintigraphic analysis proceeds in the opposite order: The fraction of radioactivity entering the gut is determined from the images, and then remaining loss of activity is considered to be eliminated in the urine.
Gamma camera ROI analysis was used to predict the percentage of total radioactivity excreted in urine for volunteer 3, who had failed to collect a final urine sample. A renal clearance of 111 mL/min was determined for this individual, which was comparable with those of the other 5 volunteers and showed an advantage of gamma camera scintigraphic analysis to compensate for the subject’s noncompliance.
Table 4 shows estimates of TB RT and MRT and the percentage dose eliminated in urine using the 2 methods for each of the volunteers in this study. Percentage dose in urine was based on biologic elimination routes only and did not include clearance caused by physical decay of the radionuclide. Although the 2 methods used completely different approaches, included different arrays of assumptions, and had different errors associated with them, the resulting values were strikingly similar. This is expected when considering that the main sources of error for the ROI dataset would primarily affect the number of counts in a given ROI. Serial images on the same patient would include similar ROI error terms but produce time–count-rate curves that were parallel to, but vertically displaced from, the true distribution.
Although pharmacokinetic analysis was expected to be more accurate than the gamma camera analysis, both methods produced similar estimates of the TB MRT and the renal elimination fraction for total radioactivity. Gamma camera ROI estimation of these parameters required correct identification of the TB ROI, which, according to Equation 1, directly relies on the activity estimates in 3 other ROIs: the WB, the gut, and the urinary bladder. However, gut activity determination depends on numerous other ROIs, including large intestine overlap regions with other organs and the organs that overlap with the large intestine (i.e., the liver and the kidneys). Consequently, the excellent correlation between pharmacokinetic and gamma camera estimates in these volunteers also confirmed the methodology of the ROI drawing and overlap resolution procedures described in the experimental section.
Conclusion
The dosimetric analysis of 123I-IAZA in 6 healthy volunteers indicated that both disposition kinetics and favorable dosimetry support its clinical use for imaging tissue hypoxia. 123I-IAZA has been shown to have comparable dosimetry with several currently marketed diagnostic radiopharmaceuticals. The excellent correlation between gamma camera–based and serum-based kinetic analyses suggested that detailed pharmacokinetic analysis during dosimetry studies can serve as an important aid in validating some of the assumptions made in estimation of location and transit times of radioactive molecules from scintigraphic images.
Acknowledgments
IAZA was synthesized by Dr. Elena Atrazheva and radiolabeled by Ronald Schmidt. This research was supported in part through Alberta Cancer Board grant RI-14 and predoctoral scholarships from the Medical Research Council of Canada and the Alberta Heritage Foundation for Medical Research.
Footnotes
Received Nov. 14, 2000; revision accepted May 2, 2001.
For correspondence or reprints contact: Steve A. McQuarrie, PhD, 3118 Dentistry/Pharmacy Center, University of Alberta, Edmonton, Alberta, Canada T6G 2N8.