Prediction of therapy tumor-absorbed dose estimates in I-131 radioimmunotherapy using tracer data via a mixed-model fit to time activity

Accuracy and uncertainty analysis of reduced time point imaging effect on time-integrated activity for 177Lu-DOTATATE PRRT in patients and clinically realistic simulations

BACKGROUND

Dosimetry promises many advantages for radiopharmaceutical therapies but repeat post-therapy imaging for dosimetry can burden both patients and clinics. Recent applications of reduced time point imaging for time-integrated activity (TIA) determination for internal dosimetry following ¹⁷⁷Lu-DOTATATE peptide receptor radionuclide therapy have shown promising results that allow for the simplification of patient-specific dosimetry. However, factors such as scheduling can lead to sub-optimal imaging time points, but the resulting impact on dosimetry accuracy is still under investigation. We use four-time point ¹⁷⁷Lu SPECT/CT data for a cohort of patients treated at our clinic to perform a comprehensive analysis of the error and variability in time-integrated activity when reduced time point methods with various combinations of sampling points are employed.

METHODS

The study includes 28 patients with gastroenteropancreatic neuroendocrine tumors who underwent post-therapy SPECT/CT imaging at approximately 4, 24, 96, and 168 h post-therapy (p.t.) following the first cycle of ¹⁷⁷Lu-DOTATATE. The healthy liver, left/right kidney, spleen and up to 5 index tumors were delineated for each patient. Time-activity curves were fit with either monoexponential or biexponential functions for each structure, based on the Akaike information criterion. This fitting was performed using all 4 time points as a reference and various combinations of 2 and 3 time points to determine optimal imaging schedules and associated errors. 2 commonly used methods of single time point (STP) TIA estimation are also evaluated. A simulation study was also performed with data generated by sampling curve fit parameters from log-normal distributions derived from the clinical data and adding realistic measurement noise to sampled activities. For both clinical and simulation studies, error and variability in TIA estimates were estimated with various sampling schedules.

RESULTS

The optimal post-therapy imaging time period for STP estimates of TIA was found to be 3-5 days (71-126 h) p.t. for tumor and organs, with one exception of 6-8 days (144-194 h) p.t. for spleen with one STP approach. At the optimal time point, STP estimates give mean percent errors (MPE) within ± 5% and SD < 9% across all structures with largest magnitude error for kidney TIA (MPE = - 4.1%) and highest variability also for kidney TIA (SD = 8.4%). The optimal sampling schedule for 2TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. for kidney, tumor, and spleen. Using the optimal sampling schedule, the largest magnitude MPE for 2TP estimates is 1.2% for spleen and highest variability is in tumor with SD = 5.8%. The optimal sampling schedule for 3TP estimates of TIA is 1-2 days (21-52 h) p.t. followed by 3-5 days (71-126 h) p.t. and 6-8 days (144-194 h) p.t. for all structures. Using the optimal sampling schedule, the largest magnitude MPE for 3TP estimates is 2.5% for spleen and highest variability is in tumor with SD = 2.1%. Simulated patient results corroborate these findings with similar optimal sampling schedules and errors. Many sub-optimal reduced time point sampling schedules also exhibit low error and variability.

CONCLUSIONS

We show that reduced time point methods can be used to achieve acceptable average TIA errors over a wide range of imaging time points and sampling schedules while maintaining low uncertainty. This information can improve the feasibility of dosimetry for ¹⁷⁷Lu-DOTATATE and elucidate the uncertainty associated with non-ideal conditions.

A Novel Time-Activity Information-Sharing Approach Using Nonlinear Mixed Models for Patient-Specific Dosimetry with Reduced Imaging Time Points: Application in SPECT/CT After 177Lu-DOTATATE

Multiple-time-point SPECT/CT imaging for dosimetry is burdensome for patients and lacks statistical efficiency. A novel method for joint kidney time-activity estimation based on a statistical mixed model, a prior cohort of patients with complete time-activity data, and only 1 or 2 imaging points for new patients was compared with previously proposed single-time-point methods in virtual and clinical patient data. Methods: Data were available for 10 patients with neuroendocrine tumors treated with 177Lu-DOTATATE and imaged up to 4 times between days 0 and 7 using SPECT/CT. Mixed models using 1 or 2 time points were evaluated retrospectively in the clinical cohort, using the multiple-time-point fit as the reference. Time-activity data for 250 virtual patients were generated using parameter values from the clinical cohort. Mixed models were fit using 1 (∼96 h) and 2 (4 h, ∼96 h) time points for each virtual patient combined with complete data for the other patients in each dataset. Time-integrated activities (TIAs) calculated from mixed model fits and other reduced-time-point methods were compared with known values. Results: All mixed models and single-time-point methods performed well overall, achieving mean bias < 7% in the virtual cohort. Mixed models exhibited lower bias, greater precision, and substantially fewer outliers than did single-time-point methods. For clinical patients, 1- and 2-time-point mixed models resulted in more accurate TIA estimates for 94% (17/18) and 72% (13/18) of kidneys, respectively. In virtual patients, mixed models resulted in more than a 2-fold reduction in the proportion of kidneys with |bias| > 10% (6% vs. 15%). Conclusion: Mixed models based on a historical cohort of patients with complete time-activity data and new patients with only 1 or 2 SPECT/CT scans demonstrate less bias on average and significantly fewer outliers when estimating kidney TIA, compared with popular reduced-time-point methods. Use of mixed models allows for reduction of the imaging burden while maintaining accuracy, which is crucial for clinical implementation of dosimetry-based treatment.

Mapping nephron mass in vivo using positron emission tomography

Nephron number varies widely in humans. A low nephron endowment at birth or a loss of functioning nephrons is strongly linked to increased susceptibility to chronic kidney disease. In this work, we developed a contrast agent, radiolabeled cationic ferritin (RadioCF), to map functioning glomeruli in vivo in the kidney using positron emission tomography (PET). PET radiotracers can be detected in trace doses (<30 nmol), making them useful for rapid clinical translation. RadioCF is formed from cationic ferritin (CF) and with a radioisotope, Cu-64, incorporated into the ferritin core. We showed that RadioCF binds specifically to kidney glomeruli after intravenous injection in mice, whereas radiolabeled noncationic ferritin (RadioNF) and free Cu-64 do not. We then showed that RadioCF-PET can distinguish kidneys in healthy wild-type (WT) mice from kidneys in mice with oligosyndactylism (Os/+), a model of congenital hypoplasia and low nephron mass. The average standardized uptake value (SUV) measured by PET 90 min after injection was 21% higher in WT mice than in Os/+ mice, consistent with the higher glomerular density in WT mice. The difference in peak SUV from SUV at 90 min correlated with glomerular density in male mice from both WT and Os/+ cohorts (R2 = 0.98). Finally, we used RadioCF-PET to map functioning glomeruli in a donated human kidney. SUV within the kidney correlated with glomerular number (R2= 0.78) measured by CF-enhanced magnetic resonance imaging in the same locations. This work suggests that RadioCF-PET appears to accurately detect nephron mass and has the potential for clinical translation.

Standardised quantitative radioiodine SPECT/CT Imaging for multicentre dosimetry trials in molecular radiotherapy

The SEL-I-METRY trial (EudraCT No 2015-002269-47) is the first multicentre trial to investigate the role of 123I and 131I SPECT/CT-based tumour dosimetry to predict response to radioiodine therapy. Standardised dosimetry methodology is essential to provide a robust evidence-base for absorbed dose-response thresholds for molecular radiotherapy (MRT). In this paper a practical standardised protocol is used to establish the first network of centres with consistent methods of radioiodine activity quantification. Nine SPECT/CT systems at eight centres were set-up for quantitative radioiodine imaging. The dead-time of the systems was characterised for up to 2.8 GBq 131I. Volume dependent calibration factors were measured on centrally reconstructed images of 123I and 131I in six (0.8-196 ml) cylinders. Validation of image quantification using these calibration factors was performed on three systems, by imaging a 3D-printed phantom mimicking a patient's activity distribution. The percentage differences between the activities measured in the SPECT/CT image and those measured by the radionuclide calibrator were calculated. Additionally uncertainties on the SPECT/CT-based activities were calculated to indicate the limit on the quantitative accuracy of this method. For systems set-up to image high 131I count rates, the count rate versus activity did not peak below 2.8 GBq and fit a non-paralysable model. The dead-times and volume-dependent calibration factors were comparable between systems of the same model and crystal thickness. Therefore a global calibration curve could be fitted to each. The errors on the validation phantom activities' were comparable to the measurement uncertainties derived from uncertainty analysis, at 10% and 16% on average for 123I and 131I respectively in a 5 cm sphere. In conclusion, the dead-time and calibration factors varied between centres, with different models of system. However, global calibration factors may be applied to the same system model with the same crystal thickness, to simplify set-up of future multi-centre MRT studies.

Biodistribution and Dosimetry of Indigenously Produced 131I-Rituximab in B-Cell Lymphoma: Pilot Study Estimating Patient-Specific Dose Comparing 2 Different Dosimetric Methods

Cost containment through indigenous production of radioimmunotherapy agents for non-Hodgkin lymphoma (NHL) would be a pivotal step toward wider clinical availability, especially in developing countries. We examined the biodistribution and dosimetry of indigenously developed and radiolabeled ¹³¹I-rituximab, using the monoclonal antibody of chimeric origin, in patients with B-cell lymphoma for potential use in radioimmunotherapy. Methods: This prospective study included 13 patients with B-cell NHL who underwent low-dose diagnostic scanning for dosimetric and biodistribution studies. Soon after rituximab infusion, a diagnostic dose of radioiodinated rituximab was administered. Serial planar whole-body γ-camera images were taken soon afterward and on days 1, 2, 4, and 6. A source of ¹³¹I with known activity was used as a reference standard for dosimetry calculations. Results: The patient-specific administered dose that would give a whole-body absorbed radiation dose of 75 cGy, calculated by the MIRD schema, ranged from 3,095.42 to 6,330.33 MBq (83.66-171.09 mCi), with a mean of 3,986.01 ± 863.95 MBq (107.73 ± 23.35 mCi) and a median of 3,697.41 MBq (99.93 mCi). The mean residence time was 69.54 h. Within the first 48 h at least 50% of the injected activity was cleared, and by 144 h at least 80% was cleared. The patient-specific administered dose that would give a whole-body absorbed radiation dose of 75 cGy, calculated by mean residence time and activity-hours, ranged from 2,654.75 to 6,210.45 MBq (71.75-167.85 mCi), with a mean of 3,576.42 ± 927.59 MBq (96.66 ± 25.07 mCi) and a median of 3,421.02 MBq (92.46 mCi). With respect to organ-specific dosimetry, the mean absorbed doses to organs (apart from blood pool [3.77 Gy] and spleen [4.02 Gy]) were 0.97 Gy to the lungs, 0.69 Gy to the liver, and 0.7 Gy to the kidneys. Conclusion: The indigenous product had kinetics similar to commercial radiopharmaceuticals, with the advantage of a lower human antimouse antibody response because of the pharmaceutical's being a chimeric antibody rather than a murine antibody. Hence, clinical administration was safe. In none of the organs did dose-limiting radiation exposure occur at the proposed therapeutic dose.

Pretreatment CLR 124 Positron Emission Tomography Accurately Predicts CLR 131 Three-Dimensional Dosimetry in a Triple-Negative Breast Cancer Patient

INTRODUCTION

CLR1404 is a theranostic molecular agent that can be radiolabeled with ¹²⁴I (CLR 124) for positron emission tomography (PET) imaging, or ¹³¹I (CLR 131) for single-photon emission computed tomography (SPECT) imaging and targeted radionuclide therapy. This pilot study evaluated a pretreatment dosimetry methodology in a triple-negative breast cancer patient who was uniquely enrolled in both a CLR 124 PET imaging clinical trial and a CLR 131 therapeutic dose escalation clinical trial.

MATERIALS AND METHODS

Three-dimensional PET/CT images were acquired at 1, 3, 24, 48, and 120 h postinjection of 178 MBq CLR 124. One month later, pretherapy 2D whole-body planar images were acquired at 0.25, 5, 24, 48, and 144 h postinjection of 370 MBq CLR 131. Following the therapeutic administration of 1990 MBq CLR 131, 3D SPECT/CT images were acquired at 74, 147, 334, and 505 h postinjection. The therapeutic CLR 131 voxel-level absorbed dose was estimated from PET (RAPID PET) and SPECT (RAPID SPECT) images using a Geant4-based Monte Carlo dosimetry platform called RAPID (Radiopharmaceutical Assessment Platform for Internal Dosimetry), and region of interest (ROI) mean doses were also estimated using the OLINDA/EXM software based on PET (OLINDA PET), SPECT (OLINDA SPECT), and planar (OLINDA planar) images.

RESULTS

The RAPID PET and OLINDA PET tracer-predicted ROI mean doses correlated well (m ≥ 0.631, R² ≥ 0.694, p ≤ 0.01) with both the RAPID SPECT and OLINDA SPECT therapeutic mean doses. The 2D planar images did not have any significant correlations. The ROI mean doses differed by -4% to -43% between RAPID and OLINDA/EXM, and by -19% to 29% between PET and SPECT. The 3D dose distributions and dose volume histograms calculated with RAPID were similar for the PET/CT and SPECT/CT.

CONCLUSIONS

This pilot study demonstrated that CLR 124 pretreatment PET images can be used to predict CLR 131 3D therapeutic dosimetry better than CLR 131 2D planar images. In addition, unlike OLINDA/EXM, Monte Carlo dosimetry methods were capable of accurately predicting dose heterogeneity, which is important for predicting dose-response relationships and clinical outcomes.

Beyond Dose: Using Pretherapy Biomarkers to Improve Dose Prediction of Outcomes for Radioimmunotherapy of Non-Hodgkin Lymphoma

INTRODUCTION

Non-Hodgkin Lymphoma patients respond differently to therapy according to inherent biological variations. Pretherapy biomarkers may improve dose-response prediction.

MATERIALS AND METHODS

Hybrid single-photon emission computed tomography (SPECT)/computed tomography (CT) three-dimensional imaging at multiple time points plus follow-up positron emission tomography (PET)/CT or CT at 2 and 6 months post therapy were used to fit tumor response to combined biological effect and cell clearance models from which three biological effect response parameters (radiosensitivity, cold effect sensitivity, and proliferation potential) were determined per patient. A correlation of biological effect parameters and pretherapy biomarker data (ki67, p53, and phospho-histone H3) allowed a dose-based equivalent biological effect (EBE) to be calculated for each patient.

RESULTS

Significant correlations were found between biological effect parameters and pretherapy biomarkers. Optimum correlations were found by splitting the patient data according to p53 status. Response correlation of progression free survival (PFS) and EBE were significantly improved compared with PFS and absorbed dose alone.

CONCLUSIONS

It is possible and desirable to use pretherapy biomarkers to enhance the predictive potential of dose calculations for patient-specific treatment planning.

Cellular Response to Exponentially Increasing and Decreasing Dose Rates: Implications for Treatment Planning in Targeted Radionuclide Therapy

The treatment of cancer using targeted radionuclide therapy is of interest to nuclear medicine and radiation oncology because of its potential for killing tumor cells while minimizing dose-limiting toxicities to normal tissue. The ionizing radiations emitted by radiopharmaceuticals deliver radiation absorbed doses over protracted periods of time with continuously varying dose rates. As targeted radionuclide therapy becomes a more prominent part of cancer therapy, accurate models for estimating the biologically effective dose (BED) or equieffective dose (EQD2α/β) will become essential for treatment planning. This study examines the radiobiological impact of the dose rate increase half-time during the uptake phase of the radiopharmaceutical. MDA-MB-231 human breast cancer cells and V79 Chinese hamster lung fibroblasts were irradiated chronically with 662 keV γ rays delivered with time-varying dose rates that are clinically relevant. The temporal dose-rate patterns were: 1. acute, 2. exponential decrease with a half-time of 64 h (Td = 64 h), 3. initial exponential increase to a maximum (half time Ti = 2, 8 or 24 h) followed by exponential decrease (Td = 64 h). Cell survival assays were conducted and surviving fractions were determined. There was a marked reduction in biological effect when Ti was increased. Cell survival data were tested against existing dose-response models to assess their capacity to predict response. Currently accepted models that are used in radiation oncology overestimated BED and EQD2α/β at low-dose rates and underestimated them at high-dose rates. This appears to be caused by an adaptive response arising as a consequence of the initial low-dose-rate phase of exposure. An adaptive response function was derived that yields more accurate BED and EQD2α/β values over the spectrum of dose rates and absorbed doses delivered. Our experimental data demonstrate a marked increase in cell survival when the dose-rate-increase half-time is increased, thereby suggesting an adaptive response arising as a consequence of this phase of exposure. We have modified conventional radiobiological models used in the clinic for brachytherapy and external beams of radiation to account for this phenomenon and facilitate their use for treatment planning in targeted radionuclide therapy.

VIDA: a voxel-based dosimetry method for targeted radionuclide therapy using Geant4

We have developed the Voxel-Based Internal Dosimetry Application (VIDA) to provide patient-specific dosimetry in targeted radionuclide therapy performing Monte Carlo simulations of radiation transport with the Geant4 toolkit. The code generates voxel-level dose rate maps using anatomical and physiological data taken from individual patients. Voxel level dose rate curves are then fit and integrated to yield a spatial map of radiation absorbed dose. In this article, we present validation studies using established dosimetry results, including self-dose factors (DFs) from the OLINDA/EXM program for uniform activity in unit density spheres and organ self- and cross-organ DFs in the Radiation Dose Assessment Resource (RADAR) reference adult phantom. The comparison with reference data demonstrated agreement within 5% for self-DFs to spheres and reference phantom source organs for four common radionuclides used in targeted therapy ((131)I, (90)Y, (111)In, (177)Lu). Agreement within 9% was achieved for cross-organ DFs. We also present dose estimates to normal tissues and tumors from studies of two non-Hodgkin Lymphoma patients treated by (131)I radioimmunotherapy, with comparison to results generated independently with another dosimetry code. A relative difference of 12% or less was found between methods for mean absorbed tumor doses accounting for tumor regression.

Spatial resolution and image qualities of Zr-89 on Siemens Biograph TruePoint PET/CT

PURPOSE

Zirconium-89 (t(1/2)=78.41 hours) is an ideal metallic radioisotope for immuno-positron emission tomography (PET), given that its physical half-life closely matches the biological half-life of monoclonal antibodies. In this study, the authors measured the spatial resolution and image quality of Zr-89 PET and compared the results against those obtained using F-18 PET, which is widely regarded as the gold standard for comparison of imaging characteristics.

MATERIALS AND METHODS

The spatial resolution and image qualities of Zr-89 were measured on the Siemens Biograph Truepoint TrueV PET/CT scanner, partly according to NEMA NU2-2007 standards. For spatial resolution measurement, the Zr-89 point source was located at the center of the axial field of view (FOV) and offset 1/4 axial FOV from the center. For image quality measurements, an NEMA IEC Phantom was used. The NEMA IEC Phantom consists of six hot spheres that were filled with Zr-89 solution. Spatial resolution and image quality (%contrast, %background variability [BV], and source to background ratio [SBR]) were assessed to compare the imaging characteristics of F-18 with those of Siemens Biograph Truepoint TrueV.

RESULTS

The transverse and axial spatial resolutions at 1 cm were 4.5 and 4.7 mm for Zr-89, respectively. The %contrast of Zr-89 was 25.5% for the smallest 10 mm sized sphere and 89.8% for the largest 37 mm sized sphere, and for F-18, it was 32.5% for the smallest 10 mm sized sphere and 103.9% for the largest 37 mm sized sphere using the ordered subset expectation maximization (OSEM) reconstruction method. The %BV of F-18 PET was 6.4% for the smallest 10 mm sized sphere and 3.5% for the largest 37 mm sized sphere using the OSEM reconstruction. The SBR of Zr-89 was 1.8 for the smallest 10 mm sized sphere and 3.7 for the largest 37 mm sized sphere, and for F-18, it was 2.0 for the smallest 10 mm sized sphere and 4.1 for the largest 37 mm sized sphere using the OSEM reconstruction method.

CONCLUSIONS

This study assessed Zr-89 imaging characteristics using a Siemens Biograph Truepoint TrueV PET/CT scanner and compared the results with those obtained for F-18 PET. Although spatial resolution and image quality of Zr-89 PET were lower compared with F-18 PET, due to longer positron range and low positron branching ratio, Zr-89 is advantageous for immuno-PET due to well-matched half-life with monoclonal antibodies.

Local Retention and Combination Effects of Biocompatible Doxorubicin-Loaded and Radioiodine-Labeled Microhydrogels in Cancer Therapy

I-131-labeled chitosan microhydrogels (I-131-CMH) that are retained at an injection site without leaking free I-131 into normal tissue can provide opportunities to improve cancer therapy. This study focuses on the development of doxorubicin-loaded I-131-CMH (Dox-I-131-CMH) for use in radiochemotherapy against cancer. The radiolabeling of I-131-CMH was found to be stable over a period of 2 weeks with no disassociation of free I-131, and Dox showed a sustained release from the CMH. When I-131-CMH were injected into the thigh muscle or tumor tissue, in vivo gamma imaging showed a retention at the injection site with no significant leakage of I-131 into other areas of normal tissue, and after an intrahepatic arterial injection, I-131-CMH were selectively retained in the liver. Dox-I-131-CMH had significant synergistic therapeutic effects of radiation and chemotherapy on mouse breast cancer models. In this regard, Dox-I-131-CMH may be a new alternative agent for cancer therapy.

Tumor-Absorbed Dose Predicts Progression-Free Survival Following (131)I-Tositumomab Radioimmunotherapy

The study aimed at identifying patient-specific dosimetric and nondosimetric factors predicting outcome of non-Hodgkin lymphoma patients after (131)I-tositumomab radioimmunotherapy for potential use in treatment planning.

METHODS

Tumor-absorbed dose measures were estimated for 130 tumors in 39 relapsed or refractory non-Hodgkin lymphoma patients by coupling SPECT/CT imaging with the Dose Planning Method (DPM) Monte Carlo code. Equivalent biologic effect was calculated to assess the biologic effects of nonuniform absorbed dose including the effects of the unlabeled antibody. Evaluated nondosimetric covariates included histology, presence of bulky disease, and prior treatment history. Tumor level outcome was based on volume shrinkage assessed on follow-up CT. Patient level outcome measures were overall response (OR), complete response (CR), and progression-free survival (PFS), determined from clinical assessments that included PET/CT.

RESULTS

The estimated mean tumor-absorbed dose had a median value of 275 cGy (range, 94-711 cGy). A high correlation was observed between tracer-predicted and therapy-delivered mean tumor-absorbed doses (P < 0.001; r = 0.85). In univariate tumor-level analysis, tumor shrinkage correlated significantly with almost all of the evaluated dosimetric factors, including equivalent biologic effect. Regression analysis showed that OR, CR, and PFS were associated with the dosimetric factors and equivalent biologic effect. Both mean tumor-absorbed dose (P = 0.025) and equivalent biologic effect (P = 0.035) were significant predictors of PFS whereas none of the nondosimetric covariates were found to be statistically significant factors affecting PFS. The most important finding of the study was that in Kaplan-Meier curves stratified by mean dose, longer PFS was observed in patients receiving mean tumor-absorbed doses greater than 200 cGy than in those receiving 200 cGy or less (median PFS, 13.6 vs. 1.9 mo for the 2 dose groups; log-rank P < 0.0001).

CONCLUSION

A higher mean tumor-absorbed dose was significantly predictive of improved PFS after (131)I-tositumomab radioimmunotherapy. Hence tumor-absorbed dose, which can be estimated before therapy, can potentially be used to design radioimmunotherapy protocols to improve efficacy.

MIRD pamphlet No. 24: Guidelines for quantitative 131I SPECT in dosimetry applications

The reliability of radiation dose estimates in internal radionuclide therapy is directly related to the accuracy of activity estimates obtained at each imaging time point. The recently published MIRD pamphlet no. 23 provided a general overview of quantitative SPECT imaging for dosimetry. The present document is the first in a series of isotope-specific guidelines that will follow MIRD 23 and focuses on one of the most commonly used therapeutic radionuclides, (131)I. The purpose of this document is to provide guidance on the development of protocols for quantitative (131)I SPECT in radionuclide therapy applications that require regional (normal organs, lesions) and 3-dimensional dosimetry.

Biological-effect modeling of radioimmunotherapy for non-hodgkins lymphoma: determination of model parameters

Treatment with Tositumomab and 131I tositumomab anti-CD20 radioimmunotherapy (Bexxar) yields a nonradioactive antibody antitumor response (the so-called cold effect) and a radiation response. Numerical parameter determination by least-squares (LS) fitting was implemented for more accurate parameter estimates in equivalent biological-effect calculations.

METHODS

One hundred thirty-two tumors in 37 patients were followed using five or six SPECT/CT studies per patient, three each (typical) post-tracer (0.2 GBq) and post-therapy (∼3 GBq) injections. The SPECT/CT data were used to calculate position- and time-dependent dose rates and antibody concentrations for each tumor. CT-defined tumor volumes were used to track tumor volume changes. Combined biological-effect and cell-clearance models were fit to tumor volume changes. Optimized parameter values determined using LS fitting were compared to previous fitted values that were determined by matching calculated to measured tumor volume changes using visual assessment. Absorbed dose sensitivity (α) and cold-effect sensitivity (λp) parameters were the primary fitted parameters, yielding equivalent biological-effect (E) values.

RESULTS

Individual parameter uncertainties were approximately 10% and 30% for α and λp, respectively. LS versus previously fit parameter values were highly correlated, although the averaged α value decreased and the averaged λp value increased for the LS fits compared to the previous fits. Correlation of E with 2-month tumor shrinkage data was similar for the two fitting techniques. The LS fitting yielded improved fit quality and likely improved parameter estimation.