Abstract
In peptide receptor radionuclide therapy with 90Y-labeled DOTATATE, the kidney absorbed dose limits the maximum amount of total activity that can be safely administered in many patients. A higher tumor-to-kidney absorbed dose ratio might be achieved by optimizing the amount of injected peptide and activity, as recent studies have shown different degrees of receptor saturation for normal tissue and tumor. The aim of this work was to develop and implement a modeling method for treatment planning to determine the optimal combination of peptide amount and pertaining therapeutic activity for each patient. Methods: A whole-body physiologically based pharmacokinetic (PBPK) model was developed. General physiologic parameters were taken from the literature. Individual model parameters were fitted to a series (n = 12) of planar γ-camera and serum measurements (111In-DOTATATE) of patients with meningioma or neuroendocrine tumors (NETs). Using the PBPK model and the individually estimated parameters, we determined the tumor, liver, spleen, and red marrow biologically effective doses (BEDs) for a maximal kidney BED (20 Gy2.5) for different peptide amounts and activities. The optimal combination of peptide amount and activity for maximal tumor BED, considering the additional constraint of a red marrow BED less than 1 Gy15, was individually quantified. Results: The PBPK model describes the biokinetic data well considering the criteria of visual inspection, the coefficients of determination, the relative standard errors (<50%), and the correlation of the parameters (<0.8). All fitted parameters were in a physiologically reasonable range but varied considerably between patients, especially tumor perfusion (meningioma, 0.1–1 mL·g−1·min−1, and NETs, 0.02–1 mL·g−1·min−1) and receptor density (meningioma, 5–34 nmol·L−1, and NETs, 7–35 nmol·L−1). Using the proposed method, we identified the optimal amount and pertaining activity to be 76 ± 46 nmol (118 ± 71 μg) and 4.2 ± 1.8 GBq for meningioma and 87 ± 50 nmol (135 ± 78 μg) and 5.1 ± 2.8 GBq for NET patients. Conclusion: The presented work suggests that to achieve higher efficacy and safety for 90Y-DOATATE therapy, both the administered amount of peptide and the activity should be optimized in treatment planning using the proposed method. This approach could also be adapted for therapy with other peptides.
Peptide receptor radionuclide therapy (PRRT) is a promising method in the treatment of somatostatin receptor–expressing neuroendocrine tumors (NETs) (1,2). However, specific and unspecific kidney uptake limits a more effective treatment (3). Although, the administration of amino acids has considerably reduced unspecific uptake, the absorbed dose to the kidneys is still the most relevant hindrance to increasing the tumor absorbed dose (3).
A strategy to further improve the ratio of tumor to kidney absorbed dose might be adjusting the amount of injected peptide. The effect of peptide amount was investigated in animal studies (4,5). In humans, a strong effect of substance amount on the therapeutic index has been observed in previous work with monoclonal antibodies (6–8) but also peptides (9,10). Sabet et al. showed that saturation of tumor-binding sites requires considerably higher amounts of peptide than for those in normal tissue (10). Recently, we developed a pharmacokinetic model to demonstrate that different amounts of peptide for dosimetry and therapy could lead to considerably different absorbed doses (9). Until now, the effects of peptide amount on the tumor–to–kidney biologically effective dose (BED) ratio has not been considered in treatment planning.
Therefore, the aim of this work was to develop a treatment planning approach that allows for the estimation of the optimal combination of amount and activity for PRRT for a given maximal kidney BED. For that purpose, a whole-body physiologically based pharmacokinetic (PBPK) model was developed. The PBPK model was fitted to time–activity data derived from 12 series of planar γ-camera images and serum data. On the basis of the individually estimated parameters, the tumor, red marrow, liver, and spleen BEDs were calculated for increasing amounts of peptide using the boundary condition of a fixed BED to the kidneys of 20 Gy2.5.
MATERIALS AND METHODS
Patients/Data
Nine patients with metastasizing NETs (n = 5) or meningioma (n = 4) were included in this study. The institutional review board (or equivalent) approved this retrospective study, and the requirement to obtain informed consent was waived. For patients 1 and 2, pretherapeutic dosimetry was conducted for 3 and 2 cycles of PRRT, respectively. In total, 12 datasets with the same sampling schedule, with an additional CT for volume measurement of the liver, kidneys, spleen, and 1 or 2 tumor lesions and with a glomerular filtration rate measurement (51Cr-ethylenediaminetetraacetic acid [EDTA]), were included (Table 1). For application, 111In-DOTATATE was prepared according to a protocol (11) that was a modified procedure as published earlier (12). 111In-DOTATATE in phosphate-buffered saline showed a 98% or greater purity at the end of synthesis, and thus it was used directly without further purification. However, quality control was performed by radio–high-performance liquid chromatography using the same conditions as reported before for 90Y-DOTA-DUPA-Pep (11). For pretherapeutic measurements, 116 ± 16 μg (75 ± 10 nmol) DOTATATE labeled (12) with (140 ± 14 MBq) 111In were intravenously injected as a 51 ± 8 min infusion. Lysine and arginine (1,000 mL, 2.5% infusion) were coadministered over 2 h, beginning 0.5 h before the administration of 111In-DOTATATE. Blood was drawn at 5 and 15 min; 0.5, 1, 2, and 4 h; and 1, 2, and 3 d after the end of the 111In-DOTATATE infusion. Serum activity was measured using a γ-counter (Auto-γ-5003; Canberra Packard). Planar whole-body scintigraphy (anterior and posterior) with a double-head γ-camera (ECAM; Siemens) was performed at 2 and 4 h and 1, 2, and 3 d after infusion (Fig. 1). For the first measurement series (N1), the reference time was the start of the infusion. To obtain the time–activity data, ULMDOS (13) and NUKDOS (14), which contain background correction, self-attenuation, and scatter corrections according to MIRD pamphlet 16 (15) for conjugate view counting, were used. Regions of interest were drawn for 1 (for patients 2, 4, and 5) or 2 tumor sites, liver, left kidney, spleen, and total body. For red marrow, no direct time–activity data could be obtained. However, serum measurements and whole-body data were available (supplemental data; supplemental materials are available at http://jnm.snmjournals.org).
PBPK Model Structure
A whole-body PBPK model was developed. It describes all major physiologic and physical mechanisms—that is, distribution via blood flow, extravasation, specific binding, internalization, degradation and release, physical decay, and clearance (supplemental data). The model consists of 2 systems, 1 for labeled and 1 for unlabeled peptide (6). The systems are coupled by the competition for binding to free receptors and by physical decay (6). All physiologic parameters are assumed equal for the labeled and unlabeled substance. The organs are connected via blood flow.
For those somatostatin receptor subtype 2 (sst2)–positive organs that were not delineated, receptor densities relative to the kidney were calculated using values for regional sst2 expression reported in Boy et al. (supplemental data) (16). Nonlinear peptide-receptor binding was implemented. Binding to sst1 and sst3–sst5 was neglected (17). The internalization rates, λi,int, and the total numbers of sst2 receptors, Ri,0, were assumed to be constant over time. 111In-DOTATATE was metabolized and subsequently released (λi,release) from the cell. In the model, for reasons of parsimony, the released fragments or free 111In were directly excreted from the body (18). Equal rates were assumed for 90Y-DOTATATE (19).
In contrast to the previously developed model (9), the binding of peptide to serum proteins was explicitly modeled to adequately describe serum activity for time points greater than 24 h.
Kidney uptake was modeled sst2-specific and -unspecific. All unspecific uptake mechanisms were described with the fractional flow of peptide into and out of the cells back to the serum compartment (9).
PBPK Model Parameters
All adjustable parameters are presented in Table 2. In the supplemental data, all equations and parameters are described. The dissociation constant and rate, KD and koff, for 111In-DOTATATE (for fitting) and 90Y-DOTATATE (for simulating) were set to 0.4 nmol/L (9) and 0.04 min−1 (supplemental data), respectively. It was assumed that only 2% of the peptide in the kidneys was filtered back to the serum due to the coadministration of amino acids (9). The individually measured glomerular filtration rates were used (0.66) and scaled to DOTATATE molecule size (20) because the macromolecular sieving coefficient of 51Cr-EDTA is smaller than for DOTATATE. Cell release rates, λi,release, were estimated and assumed to be different for tumor and normal tissue. In Antunes et al. (21), internalization of 111In-DOTATATE for tumor cells but not for kidney cells was directly measured. Thus, the average ratio (1.7) of internalization rates of kidney and tumor, estimated for other 111In-labeled peptides (21), was used. Preliminary fits yielded tumor internalization in the order of 10−3 min−1. Therefore, the internalization rates were fixed to λTU, int = 0.001 min−1 and λK, int = 0.0017 min−1, respectively. The internalization rates of all other sst2-positive tissues were assumed equal to that of the kidneys. For 2 patients, a sensitivity analysis for these parameters was conducted (supplemental data).
Blood flow to the tumor is an important quantity because tissue uptake is blood-flow-limited for DOTATATE. However, in our case, the earliest organ measurement was 2 h after injection and therefore it was difficult to estimate from the data for some patients. Therefore, in models 1 and 2 the blood flow rates were set to 0.9 mL·g−1·min−1 for meningioma (22) and 1 mL·g−1·min−1 for NETs (23) and to a 10-fold-lower value, respectively. In model 3, the blood flow rates were fitted. The structure and parameters of models 1–3 were identical except for the assumption regarding the tumor perfusion.
PBPK Model Fitting and Selection
SAAMII (24) and SAAMII Population kinetics (versions 2.2 and 2.0; The ε Group) were used for model development with the settings as described earlier (9). The adjustable parameters of models 1, 2, and 3 were fitted to the time–activity data of each patient. The ratio of data points to adjustable parameters was greater than 3 for all patients (supplemental data). The corrected Akaike information criterion (25) was used for model selection and parameter averaging. Models with relative standard errors greater than 50% for any estimated parameter, with elements of the correlation matrix greater than 0.8 (26) and weights less than 10% (27), were not used for model averaging.
Biodistribution Simulations and BED Calculations
To automate the simulations, the model was implemented in MATLAB/Simulink (The MathWorks Inc.). The BEDs for tumor, kidneys, spleen, liver (all delineated organs), and red marrow were calculated for 2–213 nmol DOTATATE in steps of factor 2 (injection duration, 60 min) for each patient. Figure 2 describes the automated workflow for the calculation of the BEDs for a fixed BED to the kidneys of 20 Gy2.5 and a red marrow BED less than 1 Gy15, as 2 cycles and a maximal cumulative BED for kidney 40 Gy2.5 (28) and red marrow 2 Gy15 (29) were assumed. The simulations were conducted retrospectively. For the actual treatment of the patients, the peptide amount was not optimal, but similar amounts for pretherapeutic measurements and therapy were applied. The actual amounts and activities were used to determine the actual tumor and kidney absorbed doses. In addition, the BED of the kidneys and red marrow was calculated using the herein estimated parameters for NET patients and a standard amount (150 μg) and activity (5 GBq) from Vinjamuri et al. (30) for NET patients.
S values were taken from OLINDA/EXM (31) and scaled using CT-based liver, spleen, and kidney masses. For tumor, spheric shapes were assumed. Radiobiologic parameters α/β and μ for all organs were taken from the literature (supplemental data) (28,32,33).
RESULTS
Fitting and Validation
Visual inspection showed good fits. The coefficient of determination R2 was greater than 0.9 for all curves except the spleen for N4 (0.84), N5 (0.89), N6 (0.81), N7 (0.67), and N11 (0.51); the tumor for N11 (0.80); and the kidneys for N7 (0.78). A typical fit is depicted in Figure 3. Fitted curves and estimated parameters for all measurement series are provided in the supplemental data. Models 1, 2, and 3 were most supported in 3, 1, and 8 cases, respectively. Model parameter averaging was conducted only for N8, as in all other cases only 1 model fit passed the quality criteria or had an Akaike weight ω greater than 10%. The mean values of the estimated parameters, which compare favorably to literature values, are presented in Table 2.
Simulations
Figure 4 shows the simulations for patient P2 (N5). The graphs of all other patients are presented in the supplemental data. The optimal amount and pertaining activity was 76 ± 46 nmol with 4.2 ± 1.8 GBq and 87 ± 50 nmol with 5.1 ± 2.8 GBq, for meningioma and NET patients, respectively. The results for all patients are provided in Table 3. The peptide amount affected the ratio of tumor to kidney BED strongly for measurement N2, N5, N7, and N8 and moderately for N1 and N6. For N3, N4, and N9, the effect was negligible. For the treated patients (patients 1, 2, 4, and 6), the actually used amounts were close to the optimum. However, in some patients the activity was not well selected. For the NET patients, when therapy was simulated with 5 GBq and 150 μg the median kidney BED was 11 Gy2.5 (range, 6–19 Gy2.5).
DISCUSSION
We developed a treatment planning approach based on a whole-body PBPK model that allowed for the estimation of the individually optimal combination of amount and activity in PRRT for a given maximal kidney BED.
The developed PBPK model described the data well, and the fitted parameter values were all in a physiologically reasonable range. The simulations showed that for patients with a high ratio of tumor to sst2 receptor density in combination with a low ratio of tumor to kidney perfusion, the peptide amount strongly influenced the BED ratio. Thus, the peptide amount is important for treatment planning. For the investigated patients, the mean optimal amount was 81 ± 47 nmol. The pertaining activities to these optimal amounts for the kidney constraint of a BED of 20 Gy2.5 and a red marrow BED less than 1 Gy15 could be determined with the presented approach. The therapeutic window was wide for the actually treated patients (patients 1, 2, 4, and 6) in contrast to patients 5, 7, and 8. Figure 4 shows that along the line of allowed combinations (with kidney BED = 20 Gy2.5 and red marrow BED < 1 Gy15) from 50 nmol labeled with 4 GBq to 500 nmol with 20 GBq (not recommended due to high red marrow irradiation), the tumor BED is almost equal. Therefore, the simulations showed that the injected amounts for pretherapeutic measurements (74 ± 9 nmol) and therapy (63 ± 10 nmol) for the actually treated patients led to an almost optimal tumor-to-kidney ratio, although the effect of peptide amount was not considered. However, in some patients the optimal activity was not well selected. Thus, if the amount of peptide is not optimized for therapy, at least exactly the same amount of peptide should be used for pretherapeutic measurements (on which dosimetry is based) and therapy. For our NET patients, a standard amount and activity as reported by Vinjamuri et al. (30) (5 GBq with 100–200 μg) would have led to considerable undertreatment because the median kidney BED was 11 Gy2.5 (range, 6–19 Gy2.5), assuming maximal 3 cycles.
Compared with the recently published model (9), the structure of the PBPK model presented here has been refined (supplemental data), and the simulations have been extended to calculations of the absorbed doses and BEDs for all relevant organs. The herein observed effects are in agreement with our previous results.
For accurate prediction of the therapeutic biodistribution, it is necessary that the estimated physiologic parameters using individual pretherapeutic measurements and assumptions (minor fraction of unspecific kidney uptake) must be true for higher or lower peptide amounts. That necessitates the selection of proper measurement time points. The up-to-now used sampling protocol, which was adopted from radioimmunotherapy, is not optimal. For the determination of the perfusion rate, early measurements (in the range of minutes) would improve the identifiability. This work represents a general approach to quantify and include the effect of peptide amount in treatment planning. Although this effect could be adequately demonstrated and implemented using the corrected 2-dimensional data, 3-dimensional measurement techniques would reduce bias and uncertainty in the pharmacokinetic data and would thus allow a more accurate and precise parameter estimation.
The assumptions of a constant internalization rate (i.e., an average value for the measured and simulated time course) and constant total number of receptors (sst2 is immediately recycled back to the cell surface) had to be included as the design of the up-to-now used sampling protocol did not allow more explicit modeling. The model using these assumptions could describe the data well. The effect of these assumptions is, however, not entirely clear and needs to be further investigated.
Our presented method to improve the BED to the tumor could be applied to other therapies with radiolabeled peptides or antibodies, for example, to 177Lu-labeled peptides, for which less kidney toxicity was observed.
To further validate the presented approach, different amounts of peptide and pertaining activities could be administered in the same patient for various cycles or in similar patients.
CONCLUSION
The results of this work based on PBPK modeling suggest that in PRRT with 90Y-DOTATATE the effect of peptide amount on the tumor-to-kidney BED ratio is substantial in most patients and therefore needs to be considered in PRRT. For the investigated patients, the optimal administered peptide amount and the pertaining activity could be identified. In particular, the absorbed doses to poorly perfused tumors (or tumor regions) with average to high receptor concentration might be considerably improved. The proposed treatment planning method might also be important for other therapies using radiolabeled ligands. The PBPK model can be adapted to other peptides such as prostate-specific membrane antigen ligands while using the same algorithm to find the individual optimal combination of activity and amount.
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 work was supported by the Deutsche Forschungsgemeinschaft (DFG KL 2742/2-1, BE 4393/1-1 and GL 236/11-1). No other potential conflict of interest relevant to this article was reported.
Footnotes
Published online Dec. 17, 2015.
- © 2016 by the Society of Nuclear Medicine and Molecular Imaging, Inc.
REFERENCES
- Received for publication August 3, 2015.
- Accepted for publication November 23, 2015.