Abstract
Radioimmunotherapy is a method to selectively deliver radioactivity to cancer cells via specific antibodies. A strategy to enhance the efficacy of radioimmunotherapy is the prior application of unlabeled antibody, resulting in an increase in the dose to the target tissue and a decrease in the burden to other organs. It was suggested that optimizing this approach might considerably improve radioimmunotherapy with anti-CD45 antibody. The present work develops a physiologically based pharmacokinetic model to individually determine the optimal preload for radioimmunotherapy with the YAML568 anti-CD45 antibody for each patient. Methods: A physiologically based pharmacokinetic model was developed to describe the biodistribution of anti-CD45 antibody. The transport of antibody to the organs of interest via blood flow, competitive binding of unlabeled and labeled antibody, degradation and excretion of antibody, and physical decay were included in the model. The model was fitted to the biokinetics data of 5 patients with acute myeloid leukemia. On the basis of the estimated parameters, simulations for a 0- to 534-nmol preload of unlabeled antibody were conducted and the organ residence times were calculated. Results: The measured data could be adequately described by the constructed model. The estimated numbers of accessible antigens in the respective organ, in nanomoles, were 97 ± 33 for red marrow, 49 ± 24 for liver, 34 ± 18 for spleen, 38 ± 31 for lymph nodes, and 0.9 ± 0.4 for blood. These ranges indicate high interpatient variability. The optimal amount of unlabeled antibody identified by simulations would improve the ratio of residence time in red marrow to residence time in liver by a factor of 1.6–2.4. Conclusion: The efficacy of radioimmunotherapy using anti-CD45 antibody can be considerably increased with the presented model. A more selective delivery of radioactivity to the target organ and a reduction in the toxicity to normal tissue are achieved by determining the optimal preload. Furthermore, the adverse effects of radioimmunotherapy might be drastically reduced while saving antibody expenses. The validation of the model is ongoing. The model is easily extendible and therefore most probably applicable to radioimmunotherapy of other hematologic malignancies, such as antibodies targeted to CD20, CD33, or CD66.
In radioimmunotherapy, a specific antigen is targeted to selectively deliver radioactivity to cancer cells. Radioimmunotherapy of blood-borne cancers has been the most successful because of high accessibility of the target antigen expressed on inherently radiosensitive cells (1). For treatment of acute myeloid leukemia, monoclonal antibodies have been directed to CD33 (2,3), CD66 (4–6), and CD45 (7–9), and studies using anti-CD45 have shown the most promising results (10).
In general, various methods or concepts are applied to improve the biodistribution of radiolabeled monoclonal antibodies (11). Pretargeting and dose fractionation have proved to be efficient in increasing the delivery of radiation to the target tumor (1,12). Another strategy to achieve a more favorable biodistribution for radioimmunotherapy is the application of unlabeled antibodies, that is, preloading (9,13–16). For radioimmunotherapy with anti-CD45 antibody, in some studies different amounts of unlabeled antibodies have been injected before the application of labeled antibody (9,15,16). These studies showed the positive effect of unlabeled antibodies on biodistribution and high interpatient variability. Therefore, the optimal amount of unlabeled antibody seems to depend on the individual patient (9,12). However, selective delivery has not been optimized, and an individualized approach is needed to avoid supersaturation or insufficient saturation, as both could strongly decrease the efficacy of the therapy.
An individualized method to quantify the optimal amount of administered unlabeled antibody must include all relevant information about the patient and must account for biologic variability. The use of an empiric nonparametric modeling approach (17) is not appropriate for this purpose, since it is not possible to model the saturation effect, which is based on different blood flows and numbers of antigens in the organs. For modeling the biokinetic behavior of antibodies, physiologically based pharmacokinetic (PBPK) models have the advantage of being able to determine parameters with physiologic meaning (18,19). If parameters represent a physiologic quantity, the described process is better understood and other inputs (such as the varying amounts of unlabeled antibody) can be simulated. PBPK models have been developed previously to describe the biokinetics of antibodies that are directed to solid tumors (20–22). However, to our knowledge, no model for hematologic malignancies in this form and no model with competitive binding of labeled and unlabeled antibody have been presented and evaluated for clinical relevance.
Here, a PBPK model is developed for the anti-CD45 antibody YAML568, a non−2FcR-binding rat IgG2a (9). On the basis of estimated model parameters, various preloading conditions—that is, different amounts of unlabeled antibody—were investigated to find the most favorable biodistribution. The potential improvement applying this individualized approach is estimated.
MATERIALS AND METHODS
Experimental Protocol
Patient characteristics (n = 5) are presented in Table 1. All patients signed an informed consent form for this study, which was approved by the local ethics committee and the national radiation protection authorities.
The measurements were performed as described elsewhere (9). In brief, the patients received an intravenous injection of 30–47 mg (0.5 mg/kg) of unlabeled antibody 16–34 min before the intravenous injection of 121 ± 21 MBq of 111In-labeled anti-CD45 antibodies as a slow intravenous push. Blood samples were obtained at 5, 10, and 30 min; at 1 and 2 h; and at 1 d, 2 d, 3 or 5 d, and 6 d after injection. The 111In 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 at 1 d, 2 d, 3 or 5 d, and 6 d after injection to evaluate the distribution and elimination of the radiotracer. Percentage administered doses were calculated using the ULMDOS program (23).
Model Development
A PBPK model was constructed to describe the biodistribution of labeled and unlabeled anti-CD45 antibodies using SAAM2 (University of Washington) (24). Based on the studies of Matthews et al. (25,26), the model includes blood, liver, spleen, red marrow, and the lymph nodes as sites of antigens (Figs. 1 and 2), since the CD45 antigen is expressed by most hematopoietic cells, lymphoblasts, myeloblasts, and myeloid leukemic cells but not by nonhematopoietic cells (9,10). Specific antibody uptake was incorporated from the plasma compartment of each organ to the antigen compartment directly (27), since in liver, spleen, red marrow, and lymph nodes blood flows through discontinuous capillaries (openings are up to 300-fold greater than the diameter of the antibody) (28,29). The accessibility of cancer cells is also high, since in more than 90%−95% of all cases acute myeloid leukemia is a diffuse disease with no focal or solid lesions. Thus, prompt uptake of antibody to cancer cells and to bystander cells expressing CD45 can be assumed (27). This is in contrast to PBPK models including solid tumors, where the transport of the antibody through the capillary wall to the interstitial spaces has to be explicitly considered, such as in the 2-pore model of Rippe and Haraldson (20–22).
For the blood flow, Fi, to the organs and the corresponding vascular volumes, the values of the blood circulation model for reference man were used (30). However, only vascular spaces, Vi, of the organs of interest and of the gastrointestinal tract are explicitly implemented; other vascular spaces (e.g., veins, arteries, and heart) are merged into a single compartment. To ensure parsimony, the distribution of antibody to the interstitial spaces was not explicitly incorporated into the model because ready access (due to sinusoid capillaries) of antibodies to antigens in liver, spleen, red marrow, and lymph nodes has been reported (25). Therefore, the exchange across the capillary wall of other tissue is slow, compared with the fast transport and binding to antigen sites in such organs (28).
The distribution of labeled and unlabeled antibody was modeled in 2 separate circulation systems (Figs. 1 and 2). Competitive binding between the 2 circulation systems was implemented by the constraint of a stationary number of antigens, Agi, as shown in the supplemental equations and parameters (supplemental materials are available online only at http://jnm.snmjournals.org). In addition, it is assumed that labeled antibody moves into the unlabeled antibody compartment when the radiolabel decays (rate λphys).
Other important processes are the degradation of bound (rate λdb) and unbound (rate λdu) antibody and the clearance (rate λcl) of degraded antibody. We assume that bound antibody is degraded at the same rate from all organs (27). Nonbound antibody degradation and clearance of the fraction of 111In-diethylenetriaminepentaacetic acid and 111In-low-molecular-weight biodistribution are described according to the empiric submodel of Eger et al. (Fig. 1) (27). Reported values of degradation in liver and other organs lie within a wide range, indicating the relevance of the specific chemical properties of the antibody and the applied dose. For most antibodies, Fc-receptor–mediated elimination plays an important role (21). However, the used antibody YAML568 has a modified Fc-receptor (9). To investigate the implications of that modification, we compared the model with degradation according to Eger et al. (model A, Fig. 1) to the model without degradation of unbound antibody (model B, Fig. 1). The corrected Akaike information criterion was used for model selection (31). Using this model selection approach, we calculated the probability (Akaike weights) that model A or B would be the best model.
The measured radiochemical purity was modeled by 2 injections into the respective compartments: the free antibody compartment and the MetaP compartment (Fig. 1).
Data Fitting and Parameter Estimation
The estimated parameters of the models are the numbers of antigens in all organs, Agi; the total blood serum volume, Vtotal; the degradation constant, λdb, of bound antibody; the association rate, kon; and the effective dissociation constant, KD. KD is determined individually, since it varies greatly depending on the experimental conditions (32). One factor that in our case could strongly affect the value of KD is the ratio of applied antibody and the number of antigens, which is previously not known and has to be estimated. kon is also estimated since the values are likely to be different from the values calculated in vitro (Table 2) (33). However, kon has to be determined globally, that is, fitted with the data of all patients, to keep the fitting process stable. This approximation will have a minor impact on the results, since KD is the relevant parameter, whereas kon needs to be estimated but hardly affects the results (33). Studies have shown that blood flow is among the most powerful parameter of a PBPK model (19) and that acute leukemia can alter normal blood supply (34); therefore, blood flow was fitted with a Bayesian term deriving the SD from values published by Leggett and Williams (30) as a population mean. Radioactive decay was included in the differential equations (35) provided in the supplemental data.
The model was fitted to the data using SAAM2 (24). A relative data-weighting scheme with a fractional SD of 0.1 was assigned to all (n = 160) data. The Rosenbrock least squares algorithm was used, and the convergence criterion was set to 10−5 (24).
Simulations
To find the optimal preloading for the individual patients, we conducted simulations with 14 different preloading amounts of unlabeled antibodies (0–534 nmol, corresponding to 0–80 mg) using the model with its individually estimated parameters. The area under the curve, which is the relative number of decays (and therefore proportional to the residence time), was computed using SAAM2. The ratios of the residence times in liver and red marrow were used as a measure of favorable biodistribution. In addition, the ratios of residence times, τrm/τliver, for the optimal and actual preloading amount of unlabeled antibody (0.5 mg/kg) were compared, with analysis of the potential for improving radioimmunotherapy by optimizing the preload.
The simulations and calculations were also performed for a KD 10-fold smaller because KD most presumably depends on the ratio of antibody and antigens, which is expected to decrease with a decreasing amount of preloaded antibody (32).
RESULTS
Model A, in which nonbound antibody was degraded according to the submodel of Eger et al., was found to be best, with a probability of 100% as determined by Akaike weights. Visual inspection of the graphs showed good fits (Fig. 3).
The values of the estimated parameters of model A are presented in Table 3. They indicate high interpatient variability. Especially the estimated number of antigens and the blood flows differ from patient to patient considerably. The calculated values (0.35–4 nmol) of the number of antigens in the plasma using the counted leukocyte number and the antigen expression (2 × 105 per cell) were of the same magnitude as the estimated values (Table 3). The blood volume calculated according to a standard formula (given in the supplemental data) was in good agreement with the estimated volume (mean value of the ratio of estimated and calculated value was 1.01, and the range was 0.90–1.24). The fit yielded an affinity rate, kon, of 1.62 × 107 L/mol/min, which is in the typical range (33).
The results of all simulations of each patient are depicted in Figure 4. Lack of preloading leads to the worst biodistribution, and administering too much antibody obviously saturates the target organ. The maxima of the simulated ratios of residence times, τrm/τliver, depend on the individual patients and are presented in Table 4. On the basis of the presented model, the optimal amount of unlabeled antibody would yield a residence time ratio, τrm/τliver, of 2.2–7.7, which would be a 1.6- to 2.4-fold improvement over the actually used preload (Table 4). At the same time, unlabeled antibody preload could be reduced to 24%−89% of the actually administered amount. Changing the dissociation constant, KD, to a 10-fold smaller value led to the same optimal preload and an even higher (3.2- to 7.1-fold) improvement.
DISCUSSION
To find the optimal dose of unlabeled antibody, which could further improve radioimmunotherapy using anti-CD45 antibody, we developed a PBPK model that describes measured time–activity data. PBPK models are advantageous since parameters have a physiologic meaning and thus are measurable (18). Models that are physiologically based also allow the simulation of other conditions, such as the varying amounts of unlabeled antibody. In our model, which includes blood flow to the organs and competitive binding of labeled and unlabeled antibody, some parameter values were taken from the literature and others were estimated from the individual patient data. An individualized approach for each patient was developed since it has been suggested that the optimal preloading dose is patient-dependent (9).
In vitro measurements of KD do not represent in vivo values (33). For example, a lower concentration might lead to a double binding that is up to 105-fold stronger or higher than the values determined in vitro (36). Therefore, the KD values in our model have to be considered effective values, which were estimated to be 63- to 146-fold lower than the in vitro values. Patient 2 had a 10-fold lower KD (Table 3), reflecting the higher Ak/Ag ratio of 2.5. All results for KD agree well with the measurements of other studies (33). The estimated values for KD vary with the ratio of administered antibody and antigens. This result confirms the finding of Ong and Mattes (32) that KD strongly depends on the experimental conditions.
The total plasma volume had to be fitted in our model, since standard methods to estimate the total plasma value (given in the supplemental data) do not include all patient-specific deviations, especially those caused by diseases. However, the fitted parameters agreed well with values calculated using such standard methods, demonstrating the reliability of the model. A further indication that the fitted parameters actually represent physiologic quantities is that for patients 2 and 4 hypocellularity of the red marrow was confirmed by a lower estimated number of antigens. The predosimetric calculation of antigens (other than in the blood) is not reasonable since the number of “readily accessible” (25) antigens depends on the number of leukocytes and their precursors, the actual expression of antigens per cell, and the accessibility of the cells. This makes an individualized approach indispensable. Noteworthy, the degradation rate of bound antibody is on the same order of magnitude as the value determined by Eger et al. (27).
The residence times for different preloading conditions for each organ and all patients were calculated on the basis of the individually estimated PBPK model parameters. The ratio of residence times, τrm/τliver, could be improved by a factor of up to 2.4 (Table 4). Especially for patients in whom high blood activity was observed at 10 h after injection (patients 2 and 3), a considerably improved biodistribution was obtained by the simulations with lower doses of unlabeled antibody.
For radioimmunotherapy, targeting CD33 (37) and CD66 (4) antigens leads to ratios of approximately 3 and 5, respectively, for residence time in red marrow to residence time in liver. Optimizing the amount of unlabeled antibody, we could reach ratios of 2.2–7.7 for residence time in red marrow to residence time in liver (τliver/τrm) using anti-CD45 in this study—ratios that compare favorably with the ratios of anti-CD33 and anti-CD66 antibodies.
Notably, the measured residence times obtained for our patients with 10 mg or no unlabeled antibody (data of these patients were not adequate to be fitted) (9) yield values of the same magnitude as predicted for the patients investigated here. This can be considered a first confirmation of the predictive power of the model. Nevertheless, our findings need further validation using additional measurements: the model has been established and validated in a small population (5 cases), and for the investigated patients the optimal amount of unlabeled antibody was found to be considerably lower than the actual administered amount (Table 4). Especially, the extrapolation from the injected amount of antibody to the optimal amount needs further experimental proof. A straightforward investigation of the predictive power of this model can be investigated by obtaining 2 consecutive measurements with 2 different preloading conditions in the same patient. Further simulation work has to be done to find the best combination of measurements for this validation of the model.
The validated model will be integrated into clinical routine as follows: first, standard dosimetry data are used to estimate the patient-specific model parameters, which are the basis for subsequent simulations of the therapy. Second, simulations with various amounts of unlabeled antibody (preloads) are conducted to identify the most favorable biodistribution for therapy. Third, the calculated optimal therapy is performed (including additional quality control measurements).
The presented PBPK model is flexible and extendible for radioimmunotherapy of other hematologic malignancies for which blood flow and competitive binding of antibody are relevant, such as those in which CD20, CD33, or CD66 is targeted. The fact that CD33 is also expressed on tissue macrophages in the liver (10) indicates that radioimmunotherapy with anti-CD33 might be improved with a preload optimized by the presented approach. Furthermore, our presented model is also applicable for describing the biodistribution of antibody–chemotherapy conjugates in which the application of a preload of unlabeled antibody would reduce chemical toxicity, especially of the liver (38).
CONCLUSION
The presented PBPK model adequately describes the biodistribution of anti-CD45 antibody. The need for an individualized approach has been confirmed, since the fits yielded considerable differences in the number of antigens and the blood flow to the red marrow and spleen for the investigated patients. The model allows simulating different preloading conditions to determine the optimal unlabeled amount for each individual. A 1.6- to 2.4-fold improvement of the ratio of residence time in red marrow to residence time in liver shows the great potential of this approach. The presented method is a major step toward the goal of individually determining pharmacokinetics-based optimal preloading for 90Y-labeled anti-CD45 antibody therapy of hematologic malignancies. The method potentially can be extended to radioimmunotherapy with other antibodies for which unwanted binding of labeled antibody in the liver or spleen occurs. The validation of the prediction power of the model is ongoing.
Acknowledgments
We gratefully acknowledge a grant (GL 236/7-1) from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation).
Footnotes
-
COPYRIGHT © 2009 by the Society of Nuclear Medicine, Inc.
References
- Received for publication May 8, 2008.
- Accepted for publication October 23, 2008.