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

Linear Boltzmann equation solver for voxel-level dosimetry in radiopharmaceutical therapy: Comparison with Monte Carlo and kernel convolution

BACKGROUND

With recent interest in patient-specific dosimetry for radiopharmaceutical therapy (RPT) and selective internal radiation therapy (SIRT), an increasing number of voxel-based algorithms are being evaluated. Monte Carlo (MC) radiation transport, generally considered to be the most accurate among different methods for voxel-level absorbed dose estimation, can be computationally inefficient for routine clinical use.

PURPOSE

This work demonstrates a recently implemented grid-based linear Boltzmann transport equation (LBTE) solver for fast and accurate voxel-based dosimetry in RPT and SIRT and benchmarks it against MC.

METHODS

A deterministic LBTE solver (Acuros MRT) was implemented within a commercial RPT dosimetry package (Velocity 4.1). The LBTE is directly discretized using an adaptive mesh refined grid and then the coupled photon-electron radiation transport is iteratively solved inside specified volumes to estimate radiation doses from both photons and charged particles in heterogeneous media. To evaluate the performance of the LBTE solver for RPT and SIRT applications, 177Lu SPECT/CT, 90Y PET/CT, and 131I SPECT/CT images of phantoms and patients were used. Multiple lesions (2-1052 mL) and normal organs were delineated for each study. Voxel dosimetry was performed with the LBTE solver, dose voxel kernel (DVK) convolution with density correction, and a validated in-house MC code using the same time-integrated activity and density maps as input to the different dose engines. The resulting dose maps, difference maps, and dose-volume-histogram (DVH) metrics were compared, to assess the voxel-level agreement. Evaluation of mean absorbed dose included comparison with structure-level estimates from OLINDA.

RESULTS

In the phantom inserts/compartments, the LBTE solver versus MC and DVK convolution demonstrated good agreement with mean absorbed dose and DVH metrics agreeing to within 5% except for the D90 and D70 metrics of a very low activity concentration insert of 90Y where the agreement was within 15%. In the patient studies (five patients imaged after 177Lu DOTATATE RPT, five after 90Y SIRT, and two after 131I radioimmunotherapy), in general, there was better agreement between the LBTE solver and MC than between LBTE solver and DVK convolution for mean absorbed dose and voxel-level evaluations. Across all patients for all three radionuclides, for soft tissue structures (kidney, liver, lesions), the mean absorbed dose estimates from the LBTE solver were in good agreement with those from MC (median difference < 1%, maximum 9%) and those from DVK (median difference < 5%, maximum 9%). The LBTE and OLINDA estimates for mean absorbed dose in kidneys and liver agreed to within 10%, but differences for lesions were larger with a maximum 14% for 177Lu, 23% for 90Y, and 26% for 131I. For bone regions, the agreement in mean absorbed doses between LBTE and both MC and DVK were similar (median < 11%, max 11%) while for lung the agreement between LBTE and MC (median < 1%, max 8%) was substantially better than between LBTE and DVK (median < 16%, max 33%). Voxel level estimates for soft tissue structures also showed good agreement between the LBTE solver and both MC and DVK with a median difference < 5% (maximum < 13%) for the DVH metrics with all three radionuclides. The largest difference in DVH metrics was for the D90 and D70 metric in lung and bone where the uptake was low. Here, the difference between LBTE and MC had a median value < 14% (maximum 23%) for bone and < 4% (maximum 37%) for lung, while the corresponding differences between LBTE and DVK were < 23% (maximum 31%) and < 67% (maximum 313%), respectively. For a typical patient with a matrix size of 166 × 166 × 129 (voxel size 3 × 3 × 3 mm3), voxel dosimetry using the LBTE solver was as fast as ∼2 min on a desktop computer.

CONCLUSION

Having established good agreement between the LBTE solver and MC for RPT and SIRT applications, the LBTE solver is a viable option for voxel dosimetry that can be faster than MC. Further analysis is being performed to encompass the broad range of radionuclides and conditions encountered clinically.

Fit of biokinetic data in molecular radiotherapy: a machine learning approach

BACKGROUND

In literature are reported different analytical methods (AM) to choose the proper fit model and to fit data of the time-activity curve (TAC). On the other hand, Machine Learning algorithms (ML) are increasingly used for both classification and regression tasks. The aim of this work was to investigate the possibility of employing ML both to classify the most appropriate fit model and to predict the area under the curve (τ).

METHODS

Two different ML systems have been developed for classifying the fit model and to predict the biokinetic parameters. The two systems were trained and tested with synthetic TACs simulating a whole-body Fraction Injected Activity for patients affected by metastatic Differentiated Thyroid Carcinoma, administered with [131I]I-NaI. Test performances, defined as classification accuracy (CA) and percentage difference between the actual and the estimated area under the curve (Δτ), were compared with those obtained using AM varying the number of points (N) of the TACs. A comparison between AM and ML were performed using data of 20 real patients.

RESULTS

As N varies, CA remains constant for ML (about 98%), while it improves for F-test (from 62 to 92%) and AICc (from 50 to 92%), as N increases. With AM, [Formula: see text] can reach down to - 67%, while using ML [Formula: see text] ranges within ± 25%. Using real TACs, there is a good agreement between τ obtained with ML system and AM.

CONCLUSIONS

The employing of ML systems may be feasible, having both a better classification and a better estimation of biokinetic parameters.

Predictive Dosimetry and Outcomes of Hepatocellular Carcinoma Treated by Yttrium-90 Resin Microsphere Radioembolization: A Retrospective Analysis Using Technetium-99m Macroaggregated Albumin Single Photon Emission CT/CT and Planning Software

PURPOSE

To characterize estimated mean absorbed tumor dose (ADT), objective response (OR), and estimated target dose of hepatocellular carcinoma (HCC) after resin microsphere yttrium-90 (90Y) radioembolization using partition dosimetry.

MATERIALS AND METHODS

In this retrospective, single-center study, multicompartment dosimetry of index tumors receiving 90Y radioembolization between October 2015 and June 2022 was performed using a commercial software package and pretreatment technetium-99m macroaggregated albumin single photon emission computed tomography (SPECT)/computed tomography (CT). In total, 101 patients with HCC underwent 102 treatments of 127 index tumors. Patients underwent imaging every 2-3 months after treatment to determine best response per modified Response Evaluation Criteria in Solid Tumors (mRECIST). Best response was defined as the greatest response category per mRECIST and categorized as OR or nonresponse (NR). A Cox proportional hazards model evaluated the probability of tumor OR and progression-free survival using ADT.

RESULTS

The median follow-up period was 148 days (interquartile range [IQR], 92-273 days). The median ADT of OR was 141.9 Gy (IQR, 89.4-215.8 Gy) compared with the median ADT of NR treatments of 70.8 Gy (IQR, 42.0-135.3 Gy; P < .001). Only ADT was predictive of response (hazard ratio = 2.79 [95% confidence interval {CI}: 1.44-5.40]; P = .003). At 6 months, an ADT of 157 Gy predicted 90.0% (95% CI: 41.3%-98.3%) probability of OR. At 1 year, an ADT of 157 Gy predicted 91.6% (95% CI: 78.3%-100%) probability of progression-free survival. Partition modeling and delivered activity were predictive of progression (P = .021 and P = .003, respectively).

CONCLUSIONS

For HCC treated with resin microspheres, tumors receiving higher ADT exhibited higher rates of OR. An ADT of 157 Gy predicted 90.0% OR at 6 months.

The contest between internal and external-beam dosimetry: The Zeno's paradox of Achilles and the tortoise

Radionuclide therapy, also called molecular radiotherapy (MRT), has come of age, with several novel radiopharmaceuticals being approved for clinical use or under development in the last decade. External beam radiotherapy (EBRT) is a well-established treatment modality, with about half of all oncologic patients expected to receive at least one external radiation treatment over their disease course. The efficacy and the toxicity of both types of treatment rely on the interaction of radiation with biological tissues. Dosimetry played a fundamental role in the scientific and technological evolution of EBRT, and absorbed doses to the target and to the organs at risk are calculated on a routine basis. In contrast, in MRT the usefulness of internal dosimetry has long been questioned, and a structured path to include absorbed dose calculation is missing. However, following a similar route of development as EBRT, MRT treatments could probably be optimized in a significant proportion of patients, likely based on dosimetry and radiobiology. In the present paper we describe the differences and the similarities between internal and external-beam dosimetry in the context of radiation treatments, and we retrace the main stages of their development over the last decades.

Development and Utility of an Imaging System for Internal Dosimetry of Astatine-211 in Mice

In targeted radionuclide therapy, determining the absorbed dose of the ligand distributed to the whole body is vital due to its direct influence on therapeutic and adverse effects. However, many targeted alpha therapy drugs present challenges for in vivo quantitative imaging. To address this issue, we developed a planar imaging system equipped with a cadmium telluride semiconductor detector that offers high energy resolution. This system also comprised a 3D-printed tungsten collimator optimized for high sensitivity to astatine-211, an alpha-emitting radionuclide, and adequate spatial resolution for mouse imaging. The imager revealed a spectrum with a distinct peak for X-rays from astatine-211 owing to the high energy resolution, clearly distinguishing these X-rays from the fluorescent X-rays of tungsten. High collimator efficiency (4.5 × 10-4) was achieved, with the maintenance of the spatial resolution required for discerning mouse tissues. Using this system, the activity of astatine-211 in thyroid cancer tumors with and without the expression of the sodium iodide symporter (K1-NIS/K1, respectively) was evaluated through in vivo imaging. The K1-NIS tumors had significantly higher astatine-211 activity (sign test, p = 0.031, n = 6) and significantly decreased post-treatment tumor volume (Student's t-test, p = 0.005, n = 6). The concurrent examination of intratumor drug distribution and treatment outcome could be performed with the same mice.

EANM enabling guide: how to improve the accessibility of clinical dosimetry

Dosimetry can be a useful tool for personalization of molecular radiotherapy (MRT) procedures, enabling the continuous development of theranostic concepts. However, the additional resource requirements are often seen as a barrier to implementation. This guide discusses the requirements for dosimetry and demonstrates how a dosimetry regimen can be tailored to the available facilities of a centre. The aim is to help centres wishing to initiate a dosimetry service but may not have the experience or resources of some of the more established therapy and dosimetry centres. The multidisciplinary approach and different personnel requirements are discussed and key equipment reviewed example protocols demonstrating these factors are given in the supplementary material for the main therapies carried out in nuclear medicine, including [131I]-NaI for benign thyroid disorders, [177Lu]-DOTATATE and 131I-mIBG for neuroendocrine tumours and [90Y]-microspheres for unresectable hepatic carcinoma.

Radioimmunotherapy of Non-Hodgkin B-cell Lymphoma: An update

Systemic radioimmunotherapy (RIT) is arguably the most effective and least toxic anticancer treatment for non-Hodgkin lymphoma (NHL). In treatment-naïve patients with indolent NHL, the efficacy of a single injection of RIT compares with that of multiple cycles of combination chemotherapy. However, 20 years following the approval of the first CD20-targeting radioimmunoconjugates ⁹⁰Y-Ibritumomab-tiuxetan (Zevalin) and ¹³¹I-tositumomab (Bexxar), the number of patients referred for RIT in western countries has dramatically decreased. Notwithstanding this, the development of RIT has continued. Therapeutic targets other than CD20 have been identified, new vector molecules have been produced allowing for faster delivery of RIT to the target, and innovative radionuclides with favorable physical characteristics such as alpha emitters have been more widely available. In this article, we reviewed the current status of RIT in NHL, with particular focus on recent clinical and preclinical developments.

Experimental validation of Monte Carlo dosimetry for therapeutic beta emitters with radiochromic film in a 3D-printed phantom

PURPOSE

Validation of dosimetry software, such as Monte Carlo (MC) radiation transport codes used for patient-specific absorbed dose estimation, is critical prior to their use in clinical decision making. However, direct experimental validation in the clinic is generally not performed for low/medium-energy beta emitters used in radiopharmaceutical therapy (RPT) due to the challenges of measuring energy deposited by short-range particles. Our objective was to design a practical phantom geometry for radiochromic film (RF)-based absorbed dose measurements of beta-emitting radionuclides and perform experiments to directly validate our in-house developed Dose Planning Method (DPM) MC code dedicated to internal dosimetry.

METHODS

The experimental setup was designed for measuring absorbed dose from beta emitters that have a range sufficiently penetrating to ∼200 μm in water as well as to capture any photon contributions to absorbed dose. Assayed ¹⁷⁷ Lu and ⁹⁰ Y liquid sources, 13-450 MBq estimated to deliver 0.5-10 Gy to the sensitive layer of the RF, were injected into the cavity of two 3D-printed half-cylinders that had been sealed with 12.7 μm or 25.4 μm thick Kapton Tape. A 3.8 × 6 cm strip of GafChromic EBT3 RF was sandwiched between the two taped half-cylinders. After 2-48 h exposures, films were retrieved and wipe tested for contamination. Absorbed dose to the RF was measured using a commercial triple-channel dosimetry optimization method and a calibration generated via 6 MV photon beam. Profiles were analyzed across the central 1 cm² area of the RF for validation. Eleven experiments were completed with ¹⁷⁷ Lu and nine with ⁹⁰ Y both in saline and a bone equivalent solution. Depth dose curves were generated for ¹⁷⁷ Lu and ⁹⁰ Y stacking multiple RF strips between a single filled half-cylinder and an acrylic backing. All experiments were modeled in DPM to generate voxelized MC absorbed dose estimates. We extended our study to benchmark general purpose MC codes MCNP6 and EGSnrc against the experimental results as well.

RESULTS

A total of 20 experiments showed that both the 3D-printed phantoms and the final absorbed dose values were reproducible. The agreement between the absorbed dose estimates from the RF measurements and DPM was on average -4.0% (range -10.9% to 3.2%) for all single film ¹⁷⁷ Lu experiments and was on average -1.0% (range -2.7% to 0.7%) for all single film ⁹⁰ Y experiments. Absorbed depth dose estimates by DPM agreed with RF on average 1.2% (range -8.0% to 15.2%) across all depths for ¹⁷⁷ Lu and on average 4.0% (range -5.0% to 9.3%) across all depths for ⁹⁰ Y. DPM absorbed dose estimates agreed with estimates from EGSnrc and MCNP across the board, within 4.7% and within 3.4% for ¹⁷⁷ Lu and ⁹⁰ Y respectively, for all geometries and across all depths. MC showed that absorbed dose to RF from betas was greater than 92% of the total (betas + other radiations) for ¹⁷⁷ Lu, indicating measurement of dominant beta contribution with our design.

CONCLUSIONS

The reproducible results with a RF insert in a simple phantom designed for liquid sources demonstrate that this is a reliable setup for experimentally validating dosimetry algorithms used in therapies with beta-emitting unsealed sources. Absorbed doses estimated with the DPM MC code showed close agreement with RF measurement and with results from two general purpose MC codes, thereby validating the use of this algorithms for clinical RPT dosimetry.

Radioembolization Dosimetry with Total-Body 90Y PET

Transarterial radioembolization (TARE) is a locoregional radiopharmaceutical therapy based on the delivery of radioactive ⁹⁰Y microspheres to liver tumors. The importance of personalized dosimetry to make TARE safer and more effective has been demonstrated in recent clinical studies, stressing the need for quantification of the dose-response relationship to ultimately optimize the administered activity before treatment and image it after treatment. ⁹⁰Y dosimetric studies are challenging because of the lack of accurate and precise methods but are best realized with PET combined with Monte Carlo simulations and other image modalities to calculate a segmental dose distribution. The aim of this study was to assess the suitability of imaging ⁹⁰Y PET patients with the total-body PET/CT uEXPLORER and to investigate possible improvements in TARE ⁹⁰Y PET-based dosimetry. The uEXPLORER is the first commercially available ultra-high-resolution (171 cps/kBq) total-body digital PET/CT device with a 194-cm axial PET field of view that enables the whole body to be scanned at a single bed position. Methods: Two PET/CT scanners were evaluated in this study: the Biograph mCT and the total-body uEXPLORER. Images of a National Electrical Manufacturers Association (NEMA) image-quality phantom and 2 patients were reconstructed using our standard clinical oncology protocol. A late portal phase contrast-enhanced CT scan was used to contour the liver segments and create corresponding volumes of interest. To calculate the absorbed dose, Monte Carlo simulations were performed using Geant4 Application for Tomographic Emission (GATE). The absorbed dose and dose-volume histograms were calculated for all 6 spheres (diameters ranging from 10 to 37 mm) of the NEMA phantom, the liver segments, and the entire liver. Differences between the phantom doses and an analytic ground truth were quantified through the root mean squared error. Results: The uEXPLORER showed a higher signal-to-noise ratio at 10- and 13-mm diameters, consistent with its high spatial resolution and system sensitivity. The total liver-absorbed dose showed excellent agreement between the uEXPLORER and the mCT for both patients, with differences lower than 0.2%. Larger differences of up to 60% were observed when comparing the liver segment doses. All dose-volume histograms were in good agreement, with narrower tails for the uEXPLORER in all segments, indicating lower image noise. Conclusion: This patient study is compelling for the use of total-body ⁹⁰Y PET for liver dosimetry. The uEXPLORER scanner showed a better signal-to-noise ratio than mCT, especially in lower-count regions of interest, which is expected to improve dose quantification and tumor dosimetry.

Optimized cocktail of 90Y/177Lu for radionuclide therapy of neuroendocrine tumors of various sizes: a simulation study

BACKGROUND AND OBJECTIVES

There is significant interest and potential in the treatment of neuroendocrine tumors via peptide receptor radionuclide therapy (PRRT) using one or both of 90Y and 177Lu-labeled peptides. Given the presence of different tumor sizes in patients and differing radionuclide dose delivery properties, the present study aims to use Monte Carlo simulations to estimate S-values to spherical tumors of various sizes with 90Y and 177Lu separately and in combination. The goal is to determine ratios of 90Y to 177Lu that result in the largest absorbed doses per decay of the radionuclides and the most suitable dose profiles to treat tumors of specific sizes.

MATERIAL AND METHODS

Particle transfer calculations and simulations were performed using the Monte Carlo GATE simulation software. Spherical tumors of different sizes, ranging from 0.5 to 20 mm in radius, were designed. Activities of 177Lu and 90Y, individually and in combination, were homogeneously placed within the total volume of the tumors. We determined the S-values to the tumors, and to the external volume outside of the tumors (cross-dose) which was used to approximate background tissue. The dose profiles were obtained for each of the different tumor sizes, and the uniformity of dose within each tumor was calculated.

RESULTS

For all tumor sizes, the self-dose and cross-dose per decay from 90Y were higher than that from 177Lu. We observed that 177Lu had the most uniform dose distribution within tumors with radii less than 5 mm. For tumors greater than 5 mm in radius, a ratio of 25% 90Y to 75% 177Lu resulted in the most uniform doses. When the ratio of 177Lu to 90Y was smaller, the uniformity improved more with increasing tumor size. The cross-dose stayed approximately constant for tumors larger than 15 mm for all ratios of 177Lu to 90Y. Finally, as the size of the tumor increased, differences in the S-values between different ratios of 177Lu to 90Y decreased.

CONCLUSION

Our work showed that to achieve a more uniform dose distribution within the tumor, 177Lu alone is more effective for small tumors. For medium and large tumors, a ratio of 90Y to 177Lu with more or less 177Lu, respectively, is recommended.

FDG PET/CT and Dosimetric Studies of 177Lu-Lilotomab Satetraxetan in a First-in-Human Trial for Relapsed Indolent non-Hodgkin Lymphoma-Are We Hitting the Target?

Purpose

[¹⁷⁷Lu]Lu-lilotomab satetraxetan, a novel CD37 directed radioimmunotherapy (RIT), has been investigated in a first-in-human phase 1/2a study for relapsed indolent non-Hodgkin lymphoma. In this study, new methods were assessed to calculate the mean absorbed dose to the total tumor volume, with the aim of establishing potential dose–response relationships based on 2-deoxy-2-[18F]fluoro-d-glucose (FDG) positron emission tomography (PET) parameters and clinical response. Our second aim was to study if higher total tumor burden induces reduction in the ¹⁷⁷Lu-lilotomab satetraxetan accumulation in tumor.

Procedures

Fifteen patients with different pre-dosing (non-radioactive lilotomab) regimens were included and the cohort was divided into low and high non-radioactive lilotomab pre-dosing groups for some of the analyses. ¹⁷⁷Lu-lilotomab satetraxetan was administered at dosage levels of 10, 15, or 20 MBq/kg. Mean absorbed doses to the total tumor volume (tTAD) were calculated from posttreatment single-photon emission tomography (SPECT)/computed tomography (CT) acquisitions. Total values of metabolic tumor volume (tMTV), total lesion glycolysis (tTLG) and the percent change in these parameters were calculated from FDG PET/CT performed at baseline, and at 3 and 6 months after RIT. Clinical responses were evaluated at 6 months as complete remission (CR), partial remission (PR), stable disease (SD), or progressive disease (PD).

Results

Significant decreases in tMTV and tTLG were observed at 3 months for patients receiving tTAD ≥ 200 cGy compared to patients receiving tTAD < 200 cGy (p = .03 for both). All non-responders had tTAD < 200 cGy. Large variations in tTAD were observed in responders. Reduction in ¹⁷⁷Lu-lilotomab satetraxetan uptake in tumor volume was not observed in patients with higher baseline tumor burden (tTMV).

Conclusion

tTAD of ≥ 200 cGy may prove valuable to ensure clinical response, but further studies are needed to confirm this in a larger patient population. Furthermore, this work indicates that higher baseline tumor burden (up to 585 cm³) did not induce reduction in radioimmunoconjugate accumulation in tumor.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11307-022-01731-3.

Comparison of absorbed dose extrapolation methods for mouse-to-human translation of radiolabelled macromolecules

BACKGROUND

Extrapolation of human absorbed doses (ADs) from biodistribution experiments on laboratory animals is used to predict the efficacy and toxicity profiles of new radiopharmaceuticals. Comparative studies between available animal-to-human dosimetry extrapolation methods are missing. We compared five computational methods for mice-to-human AD extrapolations, using two different radiopharmaceuticals, namely [111In]CHX-DTPA-scFv78-Fc and [68Ga]NODAGA-RGDyK. Human organ-specific time-integrated activity coefficients (TIACs) were derived from biodistribution studies previously conducted in our centre. The five computational methods adopted are based on simple direct application of mice TIACs to human organs (M1), relative mass scaling (M2), metabolic time scaling (M3), combined mass and time scaling (M4), and organ-specific allometric scaling (M5), respectively. For [68Ga]NODAGA-RGDyK, these methods for mice-to-human extrapolations were tested against the ADs obtained on patients, previously published by our group. Lastly, an average [68Ga]NODAGA-RGDyK-specific allometric parameter αnew was calculated from the organ-specific biological half-lives in mouse and humans and retrospectively applied to M3 and M4 to assess differences in human AD predictions with the α = 0.25 recommended by previous studies.

RESULTS

For both radiopharmaceuticals, the five extrapolation methods showed significantly different AD results (p < 0.0001). In general, organ ADs obtained with M3 were higher than those obtained with the other methods. For [68Ga]NODAGA-RGDyK, no significant differences were found between ADs calculated with M3 and those obtained directly on human subjects (H) (p = 0.99; average M3/H AD ratio = 1.03). All other methods for dose extrapolations resulted in ADs significantly different from those calculated directly on humans (all p ≤ 0.0001). Organ-specific allometric parameters calculated using combined experimental [68Ga]NODAGA-RGDyK mice and human biodistribution data varied significantly. ADs calculated with M3 and M4 after the application of αnew = 0.17 were significantly different from those obtained by the application of α = 0.25 (both p < 0.001).

CONCLUSIONS

Available methods for mouse-to-human dosimetry extrapolations provided significantly different results in two different experimental models. For [68Ga]NODAGA-RGDyK, the best approximation of human dosimetry was shown by M3, applying a metabolic scaling to the mouse organ TIACs. The accuracy of more refined extrapolation algorithms adopting model-specific metabolic scaling parameters should be further investigated.

Personalized Dosimetry in Targeted Radiation Therapy: A Look to Methods, Tools and Critical Aspects

Targeted radiation therapy (TRT) is a strategy increasingly adopted for the treatment of different types of cancer. The urge for optimization, as stated by the European Council Directive (2013/59/EURATOM), requires the implementation of a personalized dosimetric approach, similar to what already happens in external beam radiation therapy (EBRT). The purpose of this paper is to provide a thorough introduction to the field of personalized dosimetry in TRT, explaining its rationale in the context of optimization and describing the currently available methodologies. After listing the main therapies currently employed, the clinical workflow for the absorbed dose calculation is described, based on works of the most experienced authors in the literature and recent guidelines. Moreover, the widespread software packages for internal dosimetry are presented and critical aspects discussed. Overall, a selection of the most important and recent articles about this topic is provided.

Implications of physics, chemistry and biology for dosimetry calculations using theranostic pairs

Theranostics is an emerging paradigm that combines imaging and therapy in order to personalize patient treatment. In nuclear medicine, this is achieved by using radiopharmaceuticals that target identical molecular targets for both imaging (using emitted gamma rays) and radiopharmaceutical therapy (using emitted beta, alpha or Auger-electron particles) for the treatment of various diseases, such as cancer. If the therapeutic radiopharmaceutical cannot be imaged quantitatively, a “theranostic pair” imaging surrogate can be used to predict the absorbed radiation doses from the therapeutic radiopharmaceutical. However, theranostic dosimetry assumes that the pharmacokinetics and biodistributions of both radiopharmaceuticals in the pair are identical or very similar, an assumption that still requires further validation for many theranostic pairs. In this review, we consider both same-element and different-element theranostic pairs and attempt to determine if factors exist which may cause inaccurate dose extrapolations in theranostic dosimetry, either intrinsic (e.g. chemical differences) or extrinsic (e.g. injecting different amounts of each radiopharmaceutical) to the radiopharmaceuticals. We discuss the basis behind theranostic dosimetry and present common theranostic pairs and their therapeutic applications in oncology. We investigate general factors that could create alterations in the behavior of the radiopharmaceuticals or the quantitative accuracy of imaging them. Finally, we attempt to determine if there is evidence showing some specific pairs as suitable for theranostic dosimetry. We show that there are a variety of intrinsic and extrinsic factors which can significantly alter the behavior among pairs of radiopharmaceuticals, even if they belong to the same chemical element. More research is needed to determine the impact of these factors on theranostic dosimetry estimates and on patient outcomes, and how to correctly account for them.

Dosimetry of nuclear medicine therapies: current controversies and impact on treatment optimization

Nuclear medicine therapeutic procedures have considerably expanded over the last few years, and their number is expected to grow exponentially in the future. Internal dosimetry has significantly developed as well, but has not yet been uniformly accepted as a valuable tool for prediction of therapeutic efficacy and toxicity. In this paper, we briefly summarize some of the arguments about the implementation of internal dosimetry in clinical practice. In addition, we provide a few examples of radionuclide anticancer therapies for which internal dosimetry demonstrated a significant impact on treatment optimization and patient outcome.

Imaging and dosimetry for alpha-particle emitter radiopharmaceutical therapy: improving radiopharmaceutical therapy by looking into the black box

Radiopharmaceutical therapy using α-particle emitting radionuclides (αRPT) is a novel treatment modality that delivers highly potent alpha-particles to cancer cells or their environment. We review the advantages and challenges of imaging and dosimetry in implementing αRPT for cancer patients.

Dosimetry in Clinical Radiopharmaceutical Therapy of Cancer: Practicality Versus Perfection in Current Practice

The use of radiopharmaceutical therapies (RPTs) in the treatment of cancers is growing rapidly, with more agents becoming available for clinical use in last few years and many new RPTs being in development. Dosimetry assessment is critical for personalized RPT, insofar as administered activity should be assessed and optimized in order to maximize tumor-absorbed dose while keeping normal organs within defined safe dosages. However, many current clinical RPTs do not require patient-specific dosimetry based on current Food and Drug Administration-labeled approvals, and overall, dosimetry for RPT in clinical practice and trials is highly varied and underutilized. Several factors impede rigorous use of dosimetry, as compared with the more convenient and less resource-intensive practice of empiric dosing. We review various approaches to applying dosimetry for the assessment of activity in RPT and key clinical trials, the extent of dosimetry use, the relative pros and cons of dosimetry-based versus fixed activity, and practical limiting factors pertaining to current clinical practice.

Radioimmunotherapy in Oncology: Overview of the Last Decade Clinical Trials

Simple Summary

Monoclonal antibody-bearing radionuclides have been under clinical investigation over the last two decades for their use in theranostic (diagnostic and therapeutic) applications in cancer. However, despite the numerous trials that have been conducted, only two radioimmunotherapies (RIT) have been approved by the FDA for the targeted therapy of hematologic tumors expressing CD20 antigens. Moreover, RIT applications for solid cancers faced major issues—such as radiotoxicity due to low antibodies penetrance requiring substantial curative dose—where new discoveries concerning antibody engineering or radionuclides are trying to overcome. Here, we performed an overview of the last 11-year clinical trials involving RIT for solid and non-solid cancers conducted either with full antibodies or antibody fragments. We discussed the low-to-moderate efficiency of RIT compared to conventional therapies and described the last advances in clinic for antibodies carriers (F(ab′)₂, Fab′, ScFv). Finally, we discussed about the complexity of RIT as a therapy and depicted both the issues and the prospects of such a strategy.

Abstract

The specific irradiation of tumors with selective radiolabeled antibodies constitutes an attractive therapeutic approach. Consequent preclinical research has been conducted by both biologists to identify pertinent targets and to select corresponding antibodies (mAb) and by radiochemists to radiolabel mAbs. These numerous preclinical investigations have ascertained the therapeutic interest of radioimmunotherapy (RIT) protocols in mice models. Here, we summarize the clinical studies that have been performed the last decade, including clinical trials (phases I, II, and III), prospective and retrospective studies, and cases series. We thereby reported 92 clinical studies. Among them, 62 concern the treatment of hematological malignancies, and 30 concern solid tumors. For hematologic diseases, the analysis was complex due to the high discrepancy of therapeutic strategies (first-line therapy, consolidation, stem cell transplantation conditioning) as well as the high variety of malignancies that were treated. The clinical studies from the last decade failed to expand anti-CD20 RIT indications but confirmed that RIT using radiolabeled anti-CD20 remains a pertinent choice for patients with relapse follicular lymphomas. For solid tumors, the positive benefit of RIT is more mitigated, apart for few malignancies that can be treated locally. Clinical trials also demonstrated the potential of some antibody formats, such as F(ab′)₂, which has already been approved by the China State FDA under the trend name Licartin®. Despite disparate results, mAb fragments are an interesting prospect for the improvement of RIT efficiency as well as for pretargeted strategies that delay the injection of radioactive treatments from the mAb ones.

Overview of the First NRG Oncology-National Cancer Institute Workshop on Dosimetry of Systemic Radiopharmaceutical Therapy

In 2018, the National Cancer Institute and NRG Oncology partnered for the first time to host a joint workshop on systemic radiopharmaceutical therapy (RPT) to specifically address dosimetry issues and strategies for future clinical trials. The workshop focused on current dosimetric approaches for clinical trials, strategies under development that would optimize dose reporting, and future desired or optimized approaches for novel emerging radionuclides and carriers in development. In this article, we review the main approaches that are applied clinically to calculate the absorbed dose. These include absorbed doses calculated over a variety of spatial scales, including whole body, organ, suborgan, and voxel, the last 3 of which are achievable within the MIRD schema (S value) and can be calculated with analytic methods or Monte Carlo methods, the latter in most circumstances. This article will also contrast currently available methods and tools with those used in the past, to propose a pathway whereby dosimetry helps the field by optimizing the biologic effect of the treatment and trial design in the drug approval process to reduce financial and logistical costs. We also briefly discuss the dosimetric equivalent of biomarkers to help bring a precision medicine approach to RPT implementation when merited by evidence collected during early-phase trial investigations. Advances in the methodology and related tools have made dosimetry the optimum biomarker for RPT.

Prospective SPECT-CT Organ Dosimetry-Driven Radiation-Absorbed Dose Escalation Using the In-111 (111In)/Yttrium 90 (90Y) Ibritumomab Tiuxetan (Zevalin®) Theranostic Pair in Patients with Lymphoma at Myeloablative Dose Levels

Simple Summary

We prospectively evaluated the feasibility of SPECT-CT/planar organ dosimetry-based radiation dose escalation radioimmunotherapy in patients with recurrent non-Hodgkin’s lymphoma using the theranostic pair of ¹¹¹In and ⁹⁰Y anti-CD20 ibritumomab tiuxetan (Zevalin^®) at myeloablative radiation-absorbed doses with autologous stem cell support. Unlike most routine dose escalation approaches, our approach used patient-individualized measurements of organ radiation absorbed dose from the tracer study, with patient-specific adjustments of the injected therapy dose to deliver a pre-specified radiation absorbed dose to the liver. Our approach was feasible, stem cell engraftment was swift, resulted in an 89% tumor response rate in treated patients, demonstrated over 3 fold variability in liver dosimetry/injected activity among patients, allowed us to exceed the FDA approved administered activity by over 5 fold and demonstrated the normal liver maximum tolerated dose to exceed 28 Gy. Dose escalation was not continued due to lack of drug availability. With modern dosimetry approaches, patient specific dosimetry-driven radiation dose escalation is feasible, allows adjustment of administered activity for heterogeneous pharmacokinetics and allows marked dose escalation vs. non-dosimetry driven approaches.

Abstract

Purpose: We prospectively evaluated the feasibility of SPECT-CT/planar organ dosimetry-based radiation dose escalation radioimmunotherapy in patients with recurrent non-Hodgkin’s lymphoma using the theranostic pair of ¹¹¹In and ⁹⁰Y anti-CD20 ibritumomab tiuxetan (Zevalin^®) at myeloablative radiation-absorbed doses with autologous stem cell support. We also assessed acute non-hematopoietic toxicity and early tumor response in this two-center outpatient study. Methods: 24 patients with CD20-positive relapsed or refractory rituximab-sensitive, low-grade, mantle cell, or diffuse large-cell NHL, with normal organ function, platelet counts > 75,000/mm³, and <35% tumor involvement in the marrow were treated with Rituximab (375 mg/m²) weekly for 4 consecutive weeks, then one dose of cyclophosphamide 2.5 g/m² with filgrastim 10 mcg/kg/day until stem cell collection. Of these, 18 patients with successful stem cell collection (at least 2 × 10⁶ CD34 cells/kg) proceeded to RIT. A dosimetric administration of ¹¹¹In ibritumomab tiuxetan (185 MBq) followed by five sequential quantitative planar and one SPECT/CT scan was used to determine predicted organ radiation-absorbed dose. Two weeks later, ⁹⁰Y ibritumomab tiuxetan was administered in an outpatient setting at a cohort- and patient-specific predicted organ radiation-absorbed dose guided by a Continuous Response Assessment (CRM) methodology with the following cohorts for dose escalation: 14.8 MBq/kg, and targeted 18, 24, 28, and 30.5 Gy to the liver. Autologous stem cell infusion occurred when the estimated marrow radiation-absorbed dose rate was predicted to be <1 cGy/h. Feasibility, short-term toxicities, and tumor response were assessed. Results: Patient-specific hybrid SPECT/CT + planar organ dosimetry was feasible in all 18 cases and used to determine the patient-specific therapeutic dose and guide dose escalation (26.8 ± 7.3 MBq/kg (mean), 26.3 MBq/kg (median) of ⁹⁰Y (range: 12.1–41.4 MBq/kg)) of ibritumomab tiuxetan that was required to deliver 10 Gy to the liver. Infused stem cells engrafted rapidly. The most common treatment-related toxicities were hematological and were reversible following stem cell infusion. No significant hepatotoxicity was seen. One patient died from probable treatment-related causes—pneumonia at day 27 post-transplant. One patient at dose level 18 Gy developed myelodysplastic syndrome (MDS), 4 patients required admission post-⁹⁰Y RIT for febrile neutropenia, 16/18 patients receiving ⁹⁰Y ibritumomab tiuxetan (89%) responded to the therapy, with 13 CR (72%) and 3/18 PR (17%), at 60 days post-treatment. Two patients had progressive disease at sixty days. One patient was lost to follow-up. Median time to progression was estimated to be at least 13 months. MTD to the liver is greater than 28 Gy, but the MTD was not reached as the study was terminated due to unexpected discontinuation of availability of the therapeutic agent. Conclusions: Patient-specific outpatient ⁹⁰Y ibritumomab tiuxetan RIT with myeloablative doses of RIT up to a targeted 30.5 Gy to the liver is feasible, guided by prospective SPECT/CT + planar imaging with the theranostic pair of ¹¹¹In and ⁹⁰Y anti-CD20, with outpatient autologous stem cell transplant support. Administered activity over 5 times the standard FDA-approved activity was well-tolerated. The non-hematopoietic MTD in this study exceeds 28 Gy to the liver. Initial tumor responses were common at all dose levels. This study supports the feasibility of organ dosimetry-driven patient-specific dose escalation in the treatment of NHL with stem cell transplant and provides additional information on the radiation tolerance of the normal liver to radiopharmaceutical therapy.

Recommendations for Multicentre Clinical Trials Involving Dosimetry for Molecular Radiotherapy

Multicentre clinical trials involving a dosimetry component are becoming more prevalent in molecular radiotherapy and are essential to generate the evidence to support individualised approaches to treatment planning and to ensure that sufficient patients are recruited to achieve the statistical significance required. Quality assurance programmes should be considered to support the standardisation required to achieve meaningful results. Trials should be designed to ensure that dosimetry results from image acquisition systems across centres are comparable by incorporating steps to standardise the methodologies used for the quantification of images and dosimetry. Furthermore, it is essential to assess the expertise and resources available at each participating site prior to trial commencement. A quality assurance plan should be drawn up and training provided if necessary. Standardisation of quantification and dosimetry methodologies used in a trial are essential to ensure that results from different centres may be collated. In addition, appropriate uncertainty analysis should be carried out to correct for differences in methodologies between centres. Recommendations are provided to support dosimetry studies based on the experience of several previous and ongoing multicentre trials.

FDG PET/CT parameters and correlations with tumor-absorbed doses in a phase 1 trial of 177Lu-lilotomab satetraxetan for treatment of relapsed non-Hodgkin lymphoma

Purpose

¹⁷⁷Lu-lilotomab satetraxetan targets the CD37 antigen and has been investigated in a first-in-human phase 1/2a study for relapsed non-Hodgkin lymphoma (NHL). Tumor dosimetry and response evaluation can be challenging after radioimmunotherapy (RIT). Changes in FDG PET/CT parameters after RIT and correlations with tumor-absorbed doses has not been examined previously in patients with lymphoma. Treatment-induced changes were measured at FDG PET/CT and ceCT to evaluate response at the lesion level after treatment, and correlations with tumor-absorbed doses were investigated.

Methods

Forty-five tumors in 16 patients, with different pre-treatment and pre-dosing regimens, were included. Dosimetry was performed based on multiple SPECT/CT images. FDG PET/CT was performed at baseline and at 3 and 6 months. SUV_max, MTV, TLG, and changes in these parameters were calculated for each tumor. Lesion response was evaluated at 3 and 6 months (PET_3months and PET_6months) based on Deauville criteria. Anatomical changes based on ceCT at baseline and at 6 and 12 months were investigated by the sum of perpendiculars (SPD).

Results

Tumor-absorbed doses ranged from 35 to 859 cGy. Intra- and interpatient variations were observed. Mean decreases in PET parameters from baseline to 3 months were ΔSUV_max-3months 61%, ΔMTV_3months 80%, and ΔTLG_3months 77%. There was no overall correlation between tumor-absorbed dose and change in FDG PET or ceCT parameters at the lesion level or significant difference in tumor-absorbed doses between metabolic responders and non-responders after treatment.

Conclusion

Our analysis does not show any correlation between tumor-absorbed doses and changes in FDG PET or ceCT parameters for the included lesions. The combination regimen, including cold antibodies, may be one of the factors precluding such a correlation. Increased intra-patient response with increased tumor-absorbed doses was observed for most patients, implying individual variations in radiation sensitivity or biology.

Trial registration

ClinicalTrials.gov Identifier (NCT01796171). Registered December 2012

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-020-05098-x.

Current Status of Radiopharmaceutical Therapy

In radiopharmaceutical therapy (RPT), a radionuclide is systemically or locally delivered with the goal of targeting and delivering radiation to cancer cells while minimizing radiation exposure to untargeted cells. Examples of current RPTs include thyroid ablation with the administration of ¹³¹I, treatment of liver cancer with ⁹⁰Y microspheres, the treatment of bony metastases with ²²³Ra, and the treatment of neuroendocrine tumors with ¹⁷⁷Lu-DOTATATE. New RPTs are being developed where radionuclides are incorporated into systemic targeted therapies. To assure that RPT is appropriately implemented, advances in targeting need to be matched with advances in quantitative imaging and dosimetry methods. Currently, radiopharmaceutical therapy is administered by intravenous or locoregional injection, and the treatment planning has typically been implemented like chemotherapy, where the activity administered is either fixed or based on a patient's body weight or body surface area. RPT pharmacokinetics are measurable by quantitative imaging and are known to vary across patients, both in tumors and normal tissues. Therefore, fixed or weight-based activity prescriptions are not currently optimized to deliver a cytotoxic dose to targets while remaining within the tolerance dose of organs at risk. Methods that provide dose estimates to individual patients rather than to reference geometries are needed to assess and adjust the injected RPT dose. Accurate doses to targets and organs at risk will benefit the individual patients and decrease uncertainties in clinical trials. Imaging can be used to measure activity distribution in vivo, and this information can be used to determine patient-specific treatment plans where the dose to the targets and organs at risk can be calculated. The development and adoption of imaging-based dosimetry methods is particularly beneficial in early clinical trials. In this work we discuss dosimetric accuracy needs in modern radiation oncology, uncertainties in the dosimetry in RPT, and best approaches for imaging and dosimetry of internal radionuclide therapy.

Radiopharmaceutical therapy in cancer: clinical advances and challenges

Radiopharmaceutical therapy (RPT) is emerging as a safe and effective targeted approach to treating many types of cancer. In RPT, radiation is systemically or locally delivered using pharmaceuticals that either bind preferentially to cancer cells or accumulate by physiological mechanisms. Almost all radionuclides used in RPT emit photons that can be imaged, enabling non-invasive visualization of the biodistribution of the therapeutic agent. Compared with almost all other systemic cancer treatment options, RPT has shown efficacy with minimal toxicity. With the recent FDA approval of several RPT agents, the remarkable potential of this treatment is now being recognized. This Review covers the fundamental properties, clinical development and associated challenges of RPT.

Radiopharmaceutical therapy is emerging as a safe and effective approach for the treatment of cancer, offering several advantages over existing therapeutic strategies. Here, Sgouros and colleagues provide an overview of the fundamental properties of radiopharmaceutical therapy, discuss agents in use and in clinical development and highlight the associated translational challenges.

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.

Development and Validation of RAPID: A Patient-Specific Monte Carlo Three-Dimensional Internal Dosimetry Platform

This work describes the development and validation of a patient-specific Monte Carlo internal dosimetry platform called RAPID (Radiopharmaceutical Assessment Platform for Internal Dosimetry). RAPID utilizes serial PET/CT or SPECT/CT images to calculate voxelized three-dimensional (3D) internal dose distributions with the Monte Carlo code Geant4. RAPID's dosimetry calculations were benchmarked against previously published S-values and specific absorbed fractions (SAFs) calculated for monoenergetic photon and electron sources within the Zubal phantom and for S-values calculated for a variety of radionuclides within spherical tumor phantoms with sizes ranging from 1 to 1000 g. The majority of the S-values and SAFs calculated in the Zubal Phantom were within 5% of the previously published values with the exception of a few 10 keV photon SAFs that agreed within 10%, and one value within 16%. The S-values calculated in the spherical tumor phantoms agreed within 2% for ¹⁷⁷Lu, ¹³¹I, ¹²⁵I, ¹⁸F, and ⁶⁴Cu, within 3.5% for ²¹¹At and ²¹³Bi, within 6.5% for ¹⁵³Sm, ¹¹¹In, ⁸⁹Zr, and ²²³Ra, and within 9% for ⁹⁰Y, ⁶⁸Ga, and ¹²⁴I. In conclusion, RAPID is capable of calculating accurate internal dosimetry at the voxel-level for a wide variety of radionuclides and could be a useful tool for calculating patient-specific 3D dose distributions.

Variations in the practice of molecular radiotherapy and implementation of dosimetry: results from a European survey

BACKGROUND

Currently, the implementation of dosimetry in molecular radiotherapy (MRT) is not well investigated, and in view of the Council Directive (2013/59/Euratom), there is a need to understand the current availability of dosimetry-based MRT in clinical practice and research studies. The aim of this study was to assess the current practice of MRT and dosimetry across European countries.

METHODS

An electronic questionnaire was distributed to European countries. This addressed 18 explicitly considered therapies, and for each therapy, a similar set of questions were included. Questions covered the number of patients and treatments during 2015, involvement of medical specialties and medical physicists, implementation of absorbed dose planning, post-therapy imaging and dosimetry, and the basis of therapy prescription.

RESULTS

Responses were obtained from 26 countries and 208 hospitals, administering in total 42,853 treatments. The most common therapies were 131I-NaI for benign thyroid diseases and thyroid ablation of adults. The involvement of a medical physicist (mean over all 18 therapies) was reported to be either minority or never by 32% of the responders. The percentage of responders that reported that dosimetry was included on an always/majority basis differed between the therapies and showed a median value of 36%. The highest percentages were obtained for 177Lu-PSMA therapy (100%), 90Y microspheres of glass (84%) and resin (82%), 131I-mIBG for neuroblastoma (59%), and 131I-NaI for benign thyroid diseases (54%). The majority of therapies were prescribed based on fixed-activity protocols. The highest number of absorbed-dose based prescriptions were reported for 90Y microsphere treatments in the liver (64% and 96% of responses for resin and glass, respectively), 131I-NaI treatment of benign thyroid diseases (38% of responses), and for 131I-mIBG treatment of neuroblastoma (18% of responses).

CONCLUSIONS

There is a wide variation in MRT practice across Europe and for different therapies, including the extent of medical-physicist involvement and the implementation of dosimetry-guided treatments.

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.

Voxel-based multimodel fitting method for modeling time activity curves in SPECT images

PURPOSE

Estimating the biodistribution and the pharmacokinetics from time-sequence SPECT images on a per-voxel basis is useful for studying activity nonuniformity or computing absorbed dose distributions by convolution of voxel kernels or Monte-Carlo radiation transport. Current approaches are either region-based, thus assuming uniform activity within the region, or voxel-based but using the same fitting model for all voxels.

METHODS

We propose a voxel-based multimodel fitting method (VoMM) that estimates a fitting function for each voxel by automatically selecting the most appropriate model among a predetermined set with Akaike criteria. This approach can be used to compute the time integrated activity (TIA) for all voxels in the image. To control fitting optimization that may fail due to excessive image noise, an approximated version based on trapezoid integration, named restricted method, is also studied. From this comparison, the number of failed fittings within images was estimated and analyzed. Numerical experiments were used to quantify uncertainties and feasibility was demonstrated with real patient data.

RESULTS

Regarding numerical experiments, root mean square errors of TIA obtained with VoMM were similar to those obtained with bi-exponential fitting functions, and were lower (< 5% vs.  > 10%) than with single model approaches that consider the same fitting function for all voxels. Failure rates were lower with VoMM and restricted approaches than with single-model methods. On real clinical data, VoMM was able to fit 90% of the voxels and led to less failed fits than single-model approaches. On regions of interest (ROI) analysis, the difference between ROI-based and voxel-based TIA estimations was low, less than 4%. However, the computation of the mean residence time exhibited larger differences, up to 25%.

CONCLUSIONS

The proposed voxel-based multimodel fitting method, VoMM, is feasible on patient data. VoMM leads organ-based TIA estimations similar to conventional ROI-based method. However, for pharmacokinetics analysis, studies of spatial heterogeneity or voxel-based absorbed dose assessment, VoMM could be used preferentially as it prevents model overfitting.

Impact of SPECT corrections on 3D-dosimetry for liver transarterial radioembolization using the patient relative calibration methodology

PURPOSE

Many centers aim to plan liver transarterial radioembolization (TARE) with dosimetry, even without CT-based attenuation correction (AC), or with unoptimized scatter correction (SC) methods. This work investigates the impact of presence vs absence of such corrections, and limited spatial resolution, on 3D dosimetry for TARE.

METHODS

Three voxelized phantoms were derived from CT images of real patients with different body sizes. Simulations of (99m)Tc-SPECT projections were performed with the SIMIND code, assuming three activity distributions in the liver: uniform, inside a "liver's segment," or distributing multiple uptaking nodules ("nonuniform liver"), with a tumoral liver/healthy parenchyma ratio of 5:1. Projection data were reconstructed by a commercial workstation, with OSEM protocol not specifically optimized for dosimetry (spatial resolution of 12.6 mm), with/without SC (optimized, or with parameters predefined by the manufacturer; dual energy window), and with/without AC. Activity in voxels was calculated by a relative calibration, assuming identical microspheres and (99m)Tc-SPECT counts spatial distribution. 3D dose distributions were calculated by convolution with (90)Y voxel S-values, assuming permanent trapping of microspheres. Cumulative dose-volume histograms in lesions and healthy parenchyma from different reconstructions were compared with those obtained from the reference biodistribution (the "gold standard," GS), assessing differences for D95%, D70%, and D50% (i.e., minimum value of the absorbed dose to a percentage of the irradiated volume). γ tool analysis with tolerance of 3%/13 mm was used to evaluate the agreement between GS and simulated cases. The influence of deep-breathing was studied, blurring the reference biodistributions with a 3D anisotropic gaussian kernel, and performing the simulations once again.

RESULTS

Differences of the dosimetric indicators were noticeable in some cases, always negative for lesions and distributed around zero for parenchyma. Application of AC and SC reduced systematically the differences for lesions by 5%-14% for a liver segment, and by 7%-12% for a nonuniform liver. For parenchyma, the data trend was less clear, but the overall range of variability passed from -10%/40% for a liver segment, and -10%/20% for a nonuniform liver, to -13%/6% in both cases. Applying AC, SC with preset parameters gave similar results to optimized SC, as confirmed by γ tool analysis. Moreover, γ analysis confirmed that solely AC and SC are not sufficient to obtain accurate 3D dose distribution. With breathing, the accuracy worsened severely for all dosimetric indicators, above all for lesions: with AC and optimized SC, -38%/-13% in liver's segment, -61%/-40% in the nonuniform liver. For parenchyma, D50% resulted always less sensitive to breathing and sub-optimal correction methods (difference overall range: -7%/13%).

CONCLUSIONS

Reconstruction protocol optimization, AC, SC, PVE and respiratory motion corrections should be implemented to obtain the best possible dosimetric accuracy. On the other side, thanks to the relative calibration, D50% inaccuracy for the healthy parenchyma from absence of AC was less than expected, while the optimization of SC was scarcely influent. The relative calibration therefore allows to perform TARE planning, basing on D50% for the healthy parenchyma, even without AC or with suboptimal corrections, rather than rely on nondosimetric methods.

Tumor-Absorbed Dose for Non-Hodgkin Lymphoma Patients Treated with the Anti-CD37 Antibody Radionuclide Conjugate 177Lu-Lilotomab Satetraxetan

¹⁷⁷Lu-lilotomab satetraxetan is a novel antibody radionuclide conjugate currently tested in a phase 1/2a first-in-human dosage escalation trial for patients with relapsed CD37+ indolent non-Hodgkin lymphoma. The aim of this work was to develop dosimetric methods and calculate tumor-absorbed radiation doses for patients treated with ¹⁷⁷Lu-lilotomab satetraxetan.

METHODS

Patients were treated at escalating injected activities (10, 15 and 20 MBq/kg) of ¹⁷⁷Lu-lilotomab satetraxetan and with different predosing, with or without 40 mg of unlabeled lilotomab. Eight patients were included for the tumor dosimetry study. Tumor radioactivity concentrations were calculated from SPECT acquisitions at multiple time points, and tumor masses were delineated from corresponding CT scans. Tumor-absorbed doses were then calculated using the OLINDA sphere model. To perform voxel dosimetry, the SPECT/CT data and an in-house-developed MATLAB program were combined to investigate the dose rate homogeneity.

RESULTS

Twenty-six tumors in 8 patients were ascribed a mean tumor-absorbed dose. Absorbed doses ranged from 75 to 794 cGy, with a median of 264 cGy across different dosage levels and different predosing. A significant correlation between the dosage level and tumor-absorbed dose was found. Twenty-one tumors were included for voxel dosimetry and parameters describing dose-volume coverage calculated. The investigation of intratumor voxel doses indicates that mean tumor dose is correlated to these parameters.

CONCLUSION

Tumor-absorbed doses for patients treated with ¹⁷⁷Lu-lilotomab satetraxetan are comparable to doses reported for other radioimmunotherapy compounds. Although the intertumor variability was considerable, a correlation between tumor dose and patient dosage level was found. Our results indicate that mean dose may be used as the sole dosimetric parameter on the lesion level.

A no-gold-standard technique for objective assessment of quantitative nuclear-medicine imaging methods

The objective optimization and evaluation of nuclear-medicine quantitative imaging methods using patient data is highly desirable but often hindered by the lack of a gold standard. Previously, a regression-without-truth (RWT) approach has been proposed for evaluating quantitative imaging methods in the absence of a gold standard, but this approach implicitly assumes that bounds on the distribution of true values are known. Several quantitative imaging methods in nuclear-medicine imaging measure parameters where these bounds are not known, such as the activity concentration in an organ or the volume of a tumor. We extended upon the RWT approach to develop a no-gold-standard (NGS) technique for objectively evaluating such quantitative nuclear-medicine imaging methods with patient data in the absence of any ground truth. Using the parameters estimated with the NGS technique, a figure of merit, the noise-to-slope ratio (NSR), can be computed, which can rank the methods on the basis of precision. An issue with NGS evaluation techniques is the requirement of a large number of patient studies. To reduce this requirement, the proposed method explored the use of multiple quantitative measurements from the same patient, such as the activity concentration values from different organs in the same patient. The proposed technique was evaluated using rigorous numerical experiments and using data from realistic simulation studies. The numerical experiments demonstrated that the NSR was estimated accurately using the proposed NGS technique when the bounds on the distribution of true values were not precisely known, thus serving as a very reliable metric for ranking the methods on the basis of precision. In the realistic simulation study, the NGS technique was used to rank reconstruction methods for quantitative single-photon emission computed tomography (SPECT) based on their performance on the task of estimating the mean activity concentration within a known volume of interest. Results showed that the proposed technique provided accurate ranking of the reconstruction methods for 97.5% of the 50 noise realizations. Further, the technique was robust to the choice of evaluated reconstruction methods. The simulation study pointed to possible violations of the assumptions made in the NGS technique under clinical scenarios. However, numerical experiments indicated that the NGS technique was robust in ranking methods even when there was some degree of such violation.

Activity estimation in radioimmunotherapy using magnetic nanoparticles

OBJECTIVE

Estimation of activity accumulated in tumor and organs is very important in predicting the response of radiopharmaceuticals treatment. In this study, we synthesized (177)Lutetium ((177)Lu)-trastuzumab-iron oxide nanoparticles as a double radiopharmaceutical agent for treatment and better estimation of organ activity in a new way by magnetic resonance imaging (MRI).

METHODS

(177)Lu-trastuzumab-iron oxide nanoparticles were synthesized and all the quality control tests such as labeling yield, nanoparticle size determination, stability in buffer and blood serum up to 4 d, immunoreactivity and biodistribution in normal mice were determined. In mice bearing breast tumor, liver and tumor activities were calculated with three methods: single photon emission computed tomography (SPECT), MRI and organ extraction, which were compared with each other.

RESULTS

The good results of quality control tests (labeling yield: 61%±2%, mean nanoparticle hydrodynamic size: 41±15 nm, stability in buffer: 86%±5%, stability in blood serum: 80%±3%, immunoreactivity: 80%±2%) indicated that (177)Lu-trastuzumab-iron oxide nanoparticles could be used as a double radiopharmaceutical agent in mice bearing tumor. Results showed that (177)Lu-trastuzumab-iron oxide nanoparticles with MRI had the ability to measure organ activities more accurate than SPECT.

CONCLUSIONS

Co-conjugating radiopharmaceutical to MRI contrast agents such as iron oxide nanoparticles may be a good way for better dosimetry in nuclear medicine treatment.

Introduction to radiobiology of targeted radionuclide therapy

During the last decades, new radionuclide-based targeted therapies have emerged as efficient tools for cancer treatment. Targeted radionuclide therapies (TRTs) are based on a multidisciplinary approach that involves the cooperation of specialists in several research fields. Among them, radiobiologists investigate the biological effects of ionizing radiation, specifically the molecular and cellular mechanisms involved in the radiation response. Most of the knowledge about radiation effects concerns external beam radiation therapy (EBRT) and radiobiology has then strongly contributed to the development of this therapeutic approach. Similarly, radiobiology and dosimetry are also assumed to be ways for improving TRT, in particular in the therapy of solid tumors, which are radioresistant. However, extrapolation of EBRT radiobiology to TRT is not straightforward. Indeed, the specific physical characteristics of TRT (heterogeneous and mixed irradiation, protracted exposure, and low absorbed dose rate) differ from those of conventional EBRT (homogeneous irradiation, short exposure, and high absorbed dose rate), and consequently the response of irradiated tissues might be different. Therefore, specific TRT radiobiology needs to be explored. Determining dose-effect correlation is also a prerequisite for rigorous preclinical radiobiology studies because dosimetry provides the necessary referential to all TRT situations. It is required too for developing patient-tailored TRT in the clinic in order to estimate the best dose for tumor control, while protecting the healthy tissues, thereby improving therapeutic efficacy. Finally, it will allow to determine the relative contribution of targeted effects (assumed to be dose-related) and non-targeted effects (assumed to be non-dose-related) of ionizing radiation. However, conversely to EBRT where it is routinely used, dosimetry is still challenging in TRT. Therefore, it constitutes with radiobiology, one of the main challenges of TRT in the future.