Abstract
The aim of this work was to simulate the effect of prostate-specific membrane antigen (PSMA)–positive total tumor volume (TTV) on the biologically effective doses (BEDs) to tumors and organs at risk in patients with metastatic castration-resistant prostate cancer who are undergoing 177Lu-PSMA radioligand therapy. Methods: A physiologically based pharmacokinetic model was fitted to the data of 13 patients treated with 177Lu-PSMA I&T (a PSMA inhibitor for imaging and therapy). The tumor, kidney, and salivary gland BEDs were simulated for TTVs of 0.1–10 L. The activity and peptide amounts leading to an optimal tumor-to-kidneys BED ratio were also investigated. Results: When the TTV was increased from 0.3 to 3 L, the simulated BEDs to tumors, kidneys, parotid glands, and submandibular glands decreased from 22 ± 15 to 11.0 ± 6.0 Gy1.49, 6.5 ± 2.3 to 3.7 ± 1.4 Gy2.5, 11.0 ± 2.7 to 6.4 ± 1.9 Gy4.5, and 10.9 ± 2.7 to 6.3 ± 1.9 Gy4.5, respectively (where the subscripts denote that an α/β of 1.49, 2.5, or 4.5 Gy was used to calculate the BED). The BED to the red marrow increased from 0.17 ± 0.05 to 0.32 ± 0.11 Gy15. For patients with a TTV of more than 0.3 L, the optimal amount of peptide was 273 ± 136 nmol and the optimal activity was 10.4 ± 4.4 GBq. Conclusion: This simulation study suggests that in patients with large PSMA-positive tumor volumes, higher activities and peptide amounts can be safely administered to maximize tumor BEDs without exceeding the tolerable BED to the organs at risk.
- prostate-specific membrane antigen (PSMA)
- biologically effective dose (BED)
- total tumor volume (TTV)
- physiologically based pharmacokinetic (PBPK) modeling
Radioligand therapy (RLT) using radiolabeled prostate-specific membrane antigen (PSMA)–specific peptides has been proven to be a safe, effective, and thus promising method for the treatment of metastatic castration-resistant prostate cancer (1–4). A one-activity-fits-all approach is widely used. However, RLT might be improved by considering the high interpatient pharmacokinetic variability reported in several dosimetry studies (5–8). In RLT with somatostatin receptor type 2–specific peptides, it has been shown that tumor burden, body size, and renal function significantly affect the activity concentration in organs at risk (OARs) and tumor tissue (9). However, so far, in patients with metastatic castration-resistant prostate cancer, a standard activity and peptide amount is usually injected, although the tumor burden varies considerably (from milliliters to several liters). In a recent study, it was shown that large tumor loads lead to a considerably lower SUVmean in dose-limiting organs (10). Therefore, the biologically effective doses (BEDs) to normal tissues and tumor tissue might also vary considerably (9). For example, large tumor burdens most presumably lead to a less effective treatment when a standard activity and amount of peptide are used (9). Physiologically based pharmacokinetic (PBPK) models are most suitable to systematically and quantitatively investigate important pharmacokinetic parameters, because they take into account the individual patient physiology of the OARs and the pathophysiology of the tumor (8,11,12).
The aims of this study were to quantitatively investigate, for 177Lu-labeled PSMA-specific peptides, the effects of different total tumor volumes (TTVs) and varying activities and peptide amounts on the BEDs to tumors and OARs. Therefore, the parameters of a recently developed PBPK model (8) were simultaneously fitted to the PET/CT and planar imaging data of 13 patients treated with RLT with 177Lu-PSMA I&T (1). On the basis of the estimated parameters, simulations were conducted for a range of tumor volumes (0.1–10 L) and different peptide amounts and activities.
MATERIALS AND METHODS
Patients and Therapeutic Measurements
The data of 13 patients with metastatic castration-resistant prostate cancer treated with 177Lu-PSMA I&T RLT were included (Table 1 and Supplemental Table B1; supplemental materials are available at http://jnm.snmjournals.org). These patients were also part of a recently published study (1). The inclusion criteria were at least 3 whole-body planar scintigraphs, including 1 late scintigraph (6–8 d after injection) (Supplemental Table B2). Only the first cycle of RLT was used. The patients received an average (±SD) of 91 ± 5.0 nmol of 177Lu-PSMA I&T labeled with 7.3 ± 0.34 GBq of 177Lu (1). Amino acids were coinfused as applied in peptide receptor radionuclide therapy. Quantification was performed as previously described (1). The time–activity data of the kidneys, 2 tumor lesions showing no significant overlap with other PSMA-ligand–positive tissues, and the total body were used, and the activities in underlying tissue (muscle and adipose) were corrected within the PBPK model (Supplemental Eq. B1). The institutional review board of the Technische Universität München approved the compassionate use of 177Lu-PSMA I&T in metastatic castration-resistant prostate cancer patients who had no other therapeutic options. All subjects signed a written informed consent form.
Pretherapy Measurements
The 68Ga-PSMA HBED-CC (1.6 ± 0.3 nmol labeled with 115 ± 16 MBq) PET/CT images before therapy (1–9 wk) were used to determine the volumes for tumors, kidneys, salivary glands, and lacrimal glands (1) with SyngoVia (Siemens Healthcare). TrueD (Siemens Healthcare) was used to estimate the activity concentration. Correction for background activity in the tumor lesion region of interest was performed (Supplemental Eq. B2). To estimate the volume of a single tumor lesion, the volume of interest was defined with an isocontour of 20%–50% of SUVmax, with adjustment of the volume of interest to the anatomic configuration of the lesion based on the CT scan (1). The cumulated activity and volume of all tumor lesions other than lesions 1 and 2 were obtained by adding all lesions slice by slice using isocontours of 15%–20%. This cumulated volume, which consisted of hundreds of lesions, is given the term tumor REST here. For this tumor REST mass, CT was not used. In PBPK models, measured data can be linked to several compartments (Supplemental Eqs. B1–B4). The PET/CT data of the kidney and the 2 tumor lesions were corrected for background using 2 different volumes of interest (Supplemental Eqs. B2 and B4). The PET/CT data of the tumor REST were assigned to the pertaining compartments, with a correction factor considering potential over- or underestimation (Supplemental Eq. B3; Supplemental Table C1).
Therefore, to correct for the overestimation or underestimation of the tumor volume, a correction factor was introduced and estimated in the fitting process (Supplemental Eq. B3).
For the kidneys, the activity concentration was derived using the postprocessing software Syngo TrueD (Siemens Healthcare) with 15% or 20% of SUVmax and correction of background activity (Supplemental Eq. B4).
PBPK Modeling
To investigate the influence of the TTV on the BED, the application of PBPK modeling is required, as the tumor volume and all parameters that determine the biokinetics are explicitly modeled. A whole-body PBPK model for PSMA-specific ligands was recently developed (8). To account for tumor growth between the PET/CT imaging and the therapy, the PBPK model was modified to include exponential tumor growth (Supplemental Table A1). The assignment of data to the specific compartments was also improved (part B of the supplemental materials). The model includes all physiologically and physically relevant mechanisms such as blood flow, PSMA-specific binding, internalization and release of 177Lu from the cells, excretion, and physical decay. The model and parameters are described in detail in part A of the supplemental materials.
The simulation analysis and modeling software (SAAM II and Popkinetics, version 2.2; The Epsilon Group) was used for fitting (8). For the 2 investigated tumor lesions, the receptor density, RTU,0, and the blood flow density, fTU, were estimated for each patient. The cell release rates of all tumor lesions, λTU,release, were assumed to be equal within one patient. For the kidneys, the blood flow density, fK; receptor density, [RK,0]; and release rate, λK,release, were also fitted. Because differences in salivary and lacrimal gland uptake do not considerably alter the tumor and kidney biokinetics, and planar imaging is inaccurate for these small organs, we assumed the literature values for the receptor densities and flows (8). The fitting was conducted in a Bayesian framework. Goodness of fit was evaluated using visual inspection, the SE of the fitted parameters, the correlation matrix, and the coefficient of determination (13).
BEDs
The BEDs to each lesion were calculated (Supplemental Eqs. A14–A18) on the basis of the determined parameters and the radiobiologic parameters (for tumors, ɑ/β = 1.49 Gy (14) and a repair half-life μ of 1.9 h (14)) from radiotherapy and brachytherapy. For calculation of the absorbed dose to the kidneys, salivary glands, and tumors, only the self-dose was considered. Determination of the S values for each tumor lesion and the salivary glands was based on the data of OLINDA/EXM for 177Lu for spheres (Supplemental Table A2). The BED model for red marrow was implemented according to Hindorf et al. (15) and calculated as described in part C of the supplemental materials.
Simulations
For the simulations and BED calculations, MATLAB/Simulink, version 2015b (The MathWorks, Inc.), was used. On the basis of the fitted parameters, the pharmacokinetics were simulated for, first, a TTV (i.e., the volume of the tumor REST plus the 2 tumor lesions) of 0.1, 0.3, 1, 3, and 10 L and, second, varying combinations of peptide amount (2–210 nmol) and activity leading to a kidney dose of 5 Gy2.5 (where the subscript denotes that an α/β of 2.5 Gy was used to calculate the BED). The limit of 5 Gy2.5 was chosen because it is a commonly assumed tolerable dose to the kidneys allowing multiple cycles—for example, 8 when assuming a maximal cumulative dose of 40 Gy2.5 and no additional risk factors. With these simulations, the BEDs were obtained for 2 tumor lesions, the tumor REST, the salivary glands, and the kidneys. The ratio of tumor BED to kidney BED was also determined. Mean values were derived for all patients and for just patients with a TTV of 0.3–3 L.
Dosimetry models for the salivary glands and the bone marrow (especially for highly metastasized bones and diffuse red marrow infiltration) include large uncertainties. However, to see the potential influence of the TTV (despite the known uncertainties) on the red marrow BED and salivary gland BED, the simulations were also conducted for these organs.
RESULTS
The goodness of fits was excellent for all patients as ascertained by visual inspection and the coefficient of determination (R2 > 0.9) (Supplemental Fig. C1). All elements of the correlation matrix were within −0.8 and 0.8 in all patients. The relative SD for all fitted parameters in all patients was less than 28%. All estimated parameter values (Supplemental Table C2) were in a physiologically reasonable range.
The dependence of BED on TTV is shown in Figure 1 for a typical patient (patient 4) with high tumor receptor densities (∼80 nmol/L). The results for all patients are provided in Supplemental Table C3A for BED and in Supplemental Table C3B for absorbed dose. The decreases in BED within a TTV of 0.3–3 L were considered most relevant, because the median of the patients’ TTV was 0.22 L; 10 of 13 patients had a TTV of more than 0.1 L and only one a TTV of more than 3 L. The simulated BEDs to the tumor lesions and kidneys decreased by 58.2% and 58.4%, respectively, from a TTV of 0.3 L to a TTV of 3 L (Table 2). For all patients, the changes in the tumor-to-kidney BED ratio for standard amounts of peptide and activities were maximally 12%. The BEDs to the submandibular glands, parotid glands, and red marrow changed by 59%, 59%, and 185%, respectively, from a TTV of 0.3 L to a TTV of 3 L (Table 2). For a TTV of 10 L, the average BED to the bone marrow increased to 0.39 Gy15.
Figure 2 shows the optimal amount of peptide and pertaining activity for patient 4 with the maximal BED to the kidneys fixed to 5 Gy2.5. The individual results for all patients are provided in Supplemental Table C4. Patients with a high tumor burden (>0.3 L) could be administered higher activities (10.4 ± 4.4 GBq) and peptide amounts (273 ± 136 nmol) without exceeding the tolerable BED to the kidneys.
DISCUSSION
The influence of larger tumor masses on pharmacokinetics is known as the tumor sink effect (9). This may lead to decreased activity concentrations in healthy tissue (10) but also in the tumor, with the risk of underdosing the tumor. In this work, we have quantitatively investigated the influence of TTV on the pharmacokinetics and thus on the BED to the OARs and tumors for PSMA-specific peptides. We used a recently developed PBPK model (8) that was simultaneously fitted to PET/CT and planar images to determine the pharmacokinetic parameters for each patient. The PET/CT information and the PBPK model were used to correct for nontumor activity in the tumor regions of interest on the planar images. However, SPECT/CT imaging during therapy, assuming a well-calibrated measurement system, would improve the accuracy of the individual pharmacokinetic parameters. An approximately 50-fold higher amount of peptide was used for therapeutic applications than for PET imaging. The parameters therefore describe the pharmacokinetics of smaller amounts and larger amounts, thus improving the reliability of the simulations. Within the PBPK model, the molecular weight, affinity, activity, physical decay, and peptide amount of the 68Ga-PSMA HBED-CC and 177Lu-PSMA I&T injections were considered. The fitted parameters were in a physiologically reasonable range. On the basis of these parameters (e.g., receptor densities), simulations with TTVs of between 0.1 and 10 L were conducted. The BED to the kidneys and tumors decreased considerably. Our results are remarkably similar to the SUVmean decline in kidneys (55%), submandibular glands (64%), and parotid glands (57%) between patients with a low tumor load and patients who were classified as having a “high tumor load” in a study by Gaertner et al. (10). The increase in red marrow BED with increasing tumor burden occurs because of higher total-body retention.
Garske et al., in a case report on 177Lu DOTA-octreotate therapy, noted that one lesson learned was the relevance of the amount of peptide relative to the number of receptors (16). With increasing tumor load (using the same activity and peptide amount), the tumor-to-kidney BED ratio slightly decreased as the saturation effects of the kidneys were reduced with less available peptide. However, the ratio can be improved by applying higher amounts of peptide with adapted activity. For antibodies (17–19) and peptides (8,11,20), it has been shown that an optimal tumor-to-OAR ratio exists for each patient. Using an algorithm developed by Kletting et al. (11), we determined, for each individual patient, the optimal activity and peptide amount for a maximal kidney BED of 5 Gy2.5 for 1 cycle. The tumor radiobiologic parameters from Fowler et al. (14) were used, although there is ongoing controversy about the actual range (21) and its applicability to RLT. Therefore, the absorbed doses and the BED for an α/β of 1.8, reported by Dearnaley et al. (21), are also provided in Supplemental Tables C3B and C3C. The maximal deviation for the tumor BED was 5% over the range of α/β values. The results indicate that through individualized dosimetry based on our PBPK model, the BED to the tumors could be safely maximized by optimizing the tumor-to-kidney BED ratio for a given limiting BED. This holds especially true for patients with a high tumor burden. However, because the number of possible cycles could also be limited by the bone marrow BED (Supplemental Table C4) (22) or the salivary gland BED, an overall optimization approach simultaneously taking into account all OARs and treatment effects is desirable. Future studies with quantitative 3-dimensional imaging (during therapy) are required to accurately estimate the parameters for the salivary (and lacrimal) glands. In addition, improved models for red marrow dosimetry are needed to accurately estimate the red marrow dose. The Simulink model and the algorithm can easily be modified to consider more than one OAR or limiting dose.
The idea of using the diagnostic PET/CT scan for treatment planning (23) is intriguing but requires several elements: a PPBK model that considers different affinities, peptide amounts, and half-lives of the PET and therapeutic agents; population parameters with an acceptably low SD for the release rates and the receptor densities; modeling of tumor growth and potential changes in salivary gland function due to cooling or candies; and an algorithm that objectively determines the PSMA-positive tumor volume without information on the therapeutic planar total-body time–activity data. Then, future clinical management could be changed by adapting the activity and peptide amount individually before the first therapeutic cycle. Individual adaptation is also relevant for early phases of the disease with a low tumor burden, as the optimal combination of activity and peptide amount might be considerably different from that in the patient population with a high tumor burden.
CONCLUSION
This modeling and simulation study indicated that the tumor burden considerably affects the BEDs to tumors and OARs. Our data suggest that patients with a large PSMA-ligand–positive tumor mass might be safely administered higher activities and peptide amounts to maximize the tumor BED without exceeding the critical dose to the OARs.
DISCLOSURE
This work was supported by LPDP-Lembaga Pengelola Dana Pendidikan (the Indonesia Endowment Fund for Education, Ministry of Finance, Indonesia) (s-2512/LPDP.3/2016) and the Deutsche Forschungsgemeinschaft (the German Research Foundation) (KL2742/2-1, BE4393/1-1, GL236/11-1, and SFB824 project B11). No other potential conflict of interest relevant to this article was reported.
Footnotes
- © 2018 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication October 10, 2017.
- Accepted for publication January 20, 2018.