- 1Department of Medical Physics and Radiobiology, Nuclear Medicine Research, Infrastructure (NuMeRI) NPC, Pretoria, South Africa
- 2Department of Nuclear Medicine, Steve Biko Academic Hospital, Pretoria, South Africa
- 3Department of Nuclear Medicine, Faculty of Health Sciences, University of Pretoria, Pretoria, South Africa
Radiopharmaceutical therapy has been widely adopted owing primarily to the development of novel radiopharmaceuticals. To fully utilize the potential of these RPTs in the era of precision medicine, therapy must be optimized to the patient's tumor characteristics. The vastly disparate dosimetry methodologies need to be harmonized as the first step towards this. Multiple factors play a crucial role in the shift from empirical activity administration to patient-specific dosimetry-based administrations from RPT. Factors such as variable responses seen in patients with presumably similar clinical characteristics underscore the need to standardize and validate dosimetry calculations. These efforts combined with ongoing initiatives to streamline the dosimetry process facilitate the implementation of radiomolecular precision oncology. However, various challenges hinder the widespread adoption of personalized dosimetry-based activity administration, particularly when compared to the more convenient and resource-efficient approach of empiric activity administration. This review outlines the fundamental principles, procedures, and methodologies related to image activity quantification and dosimetry with a specific focus on 177Lutetium-based radiopharmaceuticals.
1 Introduction
The medical community including nuclear medicine (NM) is strongly in favor of personalized treatment. Within the NM field, the recent advancements in molecular medicine have led to a surge in the development of radiopharmaceutical therapy (RPT). Both diagnostic and therapeutic radiopharmaceuticals necessitate dosimetry estimations before clinical use to assess radiation toxicity and overall effectiveness for disease diagnosis or treatment. Dosimetry plays a crucial role in diagnostic nuclear medicine, primarily to evaluate the cancer risk associated with imaging procedures. It entails determining the mean absorbed dose by organs based on representative anatomical models for assessing risk data related to exposed populations (1). In therapeutic applications, dosimetry focuses on organ toxicity and tumor control for individual patients. Ultimately, the medical decision to treat a patient depends on the evaluation of tumors and organs at risk (OAR) (2). Dosimetry data may contribute to adjusting the administered activity for treatment purposes. Based on the dosimetry findings, the administered activity may be increased or decreased for personalized medicine tailored to each patient's needs while minimizing the radiation-induced toxicity to the critical organs and adhering to their threshold radiation doses. The European Council (EC) Directive 2013/59/Euratom establishes fundamental safety standards for protection against the risks associated with exposure to ionising radiation (3). The new EC Directive 2013/59 mandates the justification of medical exposure to ionizing radiation and reaffirms the practice of optimization which may be achieved in close collaboration with medical physicists. According to EC Directive 2013/59/Euratom, radiation exposure for therapeutic purposes should be individually planned and verified to minimize doses to non-target volumes and tissues. The directive also specifies that patients undergoing treatment or diagnosis with radiopharmaceuticals should receive information on the risks of ionizing radiation and appropriate instructions to limit doses to persons in contact with the patient as much as possible.
Patient-specific dosimetry represents a significant opportunity to improve therapeutic outcomes by revealing critical insights into the trade-off between treatment efficacy and patient-specific toxicity. However, to date, patient-specific dosimetry has not undergone a thorough investigation to the full extent of its potential in the optimization of RPT. In this era of personalized therapy, there is a potential to combine RPT and external beam radiotherapy (EBRT) to treat cancers (4). The potential of using combinational therapy has been demonstrated in cancers such as prostate cancer and hepatocellular carcinoma (4–7). Despite both EBRT and RPT employing ionizing radiation for cancer treatment, they exhibit distinct differences, which makes it necessary to distinguish systematic radiation from RPT and EBRT. RPT represents a targeted therapy, where a radionuclide is precisely directed toward its intended target using a pharmaceutical agent specifically designed to bind to the target tissue. The three essential components of RPT involve a radioisotope with the associated chelater, a targeting agent, and a linker that connects the two (8), as illustrated in Figure 1. These RPTs are typically introduced into the bloodstream via intravenous infusion. The targeting agent actively seeks out the intended target i.e., the cancerous cells, and facilitates radiological damage from the energy released by the radionuclide to the cell deoxyribonucleic acid (DNA). This radiation-induced DNA damage may be irreversible, leading to cell death. Cancer cells are notably sensitive to such damage. The cytotoxic effect of the RPT and its extent on surrounding cells is dependent on the physical properties of the radionuclide (9, 10). In the case of beta and alpha particles, the extent to which the emitted radiation can penetrate cells bound to the radiopharmaceutical and surrounding cells varies depending on the specific radionuclide employed and the energy it releases. The radiopharmaceutical must remain within the cell long enough to induce radiation damage.
Figure 1. Illustration of the fundamental components in radiopharmaceutical therapy employing a linker to connect a target agent to the radioactive molecule (11).
The phenomenon of RPT is distinguished by the biochemical and physical processes occurring during radiological decay. Crucial patient-specific bio-kinetic data, such as cellular uptake and release, non-uniform distribution within the target, and metabolic behavior, play a pivotal role in determining patient-specific dosimetry (12). RPT employs a continuous, gradually decreasing radiation dose rate. This approach results in a more intricate temporal and spatial distribution of radiation dose when compared to EBRT (12, 13).
In contrast, EBRT uses fractionated high dose rates to induce lethal damage to cancer cells, combined with non-treatment intervals for cellular repair (12). Another key difference between these modalities, as depicted in Figure 2, is that EBRT administers a consistent absorbed dose per cell, regardless of the number of cells exposed to irradiation. Depending on the size of the tumor, the absorbed dose per cell might be changed at the position where the Bragg-peak occurs. On the other hand, in RPT, the absorbed dose per cell due to internal emissions from neighboring cells depends on the emission ranges of particles and the quantity of targeted cells (14).
Figure 2. Depiction of tumor cell irradiation with (A) external beam radiation therapy and (B) radiopharmaceutical therapy (14).
The primary objective of RPT is to selectively deliver the highest possible absorbed dose to tumors while minimizing potential toxicity to OAR. Achieving this goal ensures that the absorbed doses in OAR remain below respective tolerance levels to minimize adverse side effects (15). The effectiveness of this treatment approach depends on the targeted accumulation and prolonged retention of the radiopharmaceutical in the tumor cells, thus sparing normal tissue. The degree of cell destruction, in the cancerous and normal cells, is described by the absorbed radiation dose, representing the amount of energy transferred to the target tissue per unit mass. Dosimetry defines this energy deposition within the patient's body. The absorbed dose in a specific tissue, measured in Gray's (joules per kilogram), is a key determinant of the biological effect, such as radiation-induced cell death (16). Consequently, estimating the absorbed radiation dose i.e., dosimetry, is an essential tool for the efficacy of radiation-based treatments.
In EBRT, dosimetry is a well-established technique for routine treatment planning. Parameters in EBRT that are reasonably managed include the treatment site, the target volume to be irradiated, as well as the duration and intensity of the radiation. In contrast, the intricate pharmacokinetics and physical interactions of radiopharmaceuticals on a whole-body (WB) level are part of the internal dosimetry process needed for RPT treatments. These fundamental differences in RPT dosimetry as opposed to EBRT increase the challenge of translating research measurements into standard treatments, which usually depend on the total delivered radioactivity rather than radiation doses absorbed at individual sites (17). In contrast to the focal EBRT, radionuclide therapies are administered systemically and, as a result, can target more widespread diseases (18). Consequently, dosimetry should be given equal importance, as is the case with EBRT, as a driving force for the advancement of radionuclide therapies.
RPT is currently undergoing significant development, driven by the introduction of novel imaging and therapeutic radiopharmaceuticals. Dosimetry plays a critical role in achieving a balance between effectively delivering the maximum dose to tumor cells and minimizing damage to healthy tissues. The unique strength of RPT lies in its ability to image and quantitatively assess the likely biological outcomes of treatment through dosimetry and treatment planning. Clinical trials in their early stages should encompass both imaging and dosimetry to utilise the value of these distinctive RPT features. Through this approach, RPT clinical trials can undergo a comprehensive evaluation and be compared to established therapeutic methods (14). Accurate dosimetry will not only benefit the outcomes of individual patients but also reduce the uncertainties in clinical trials and practice (19). Implementing patient-specific treatment guided by dosimetry in RPT can enhance patient outcomes and survival (20). This approach offers multiple benefits, which include establishing optimal absorbed dose tolerance levels, understanding the dose-response relationship between tumors and normal tissue, as well as comparing dose-response patterns between different radiopharmaceuticals and patients. Personalized treatment is standard in EBRT for improved tumor control and reduced normal tissue toxicity, and this principle should be extended to RPT (12). This aligns with the EC Directive 2013/59/Euratom, which recommends implementing individual dose assessment in RPT (3). Patient-specific dosimetry in RPT has the potential to prevent under- or over-dosing associated with fixed administered activity regimes.
1.1 Therapeutic radionuclides with emphasis on 177Lutetium
An understanding of the RPT phenomenon requires awareness of three distinct types of radiation: photons (penetrating radiation), electrons, and alpha particles (non-penetrating radiation). Photons are primarily used to image the bio-distribution of RPT but are not employed for delivering localized cytotoxic therapeutic radiation. Typically, therapeutic radionuclides emit short-range charged particles such as beta, alpha particles, and Auger electrons, and possess longer physical half-lives compared to diagnostic radionuclides. Toxicity continues to be a concern with alpha-emitting radionuclides due to the nuclear recoil effect, which results in the release of radioactive daughter nuclei from initial radiopharmaceutical preparations, potentially leading to unintended irradiation of healthy tissues (21). Recoil energy experienced by daughters during alpha decay can be over 100 times larger than the binding energy of any chemical compound. As a result, bond rupture always follows alpha decay, suggesting that released daughters, often themselves alpha emitters, could be significantly harmful. Three strategies are suggested in the literature to mitigate recoil issues. These include encapsulation in a nano-carrier, rapid uptake of the alpha emitters in tumor cells, and local administration through intratumoral injection (22). Managing alpha particles presents challenges primarily linked with radioisotope daughters (23).
Alpha particles and Auger electrons, characterized by higher linear energy transfers (LET), deposit their energy over shorter distances, leading to increased cell death and limited repair of DNA damage (24). Smaller tumors respond well to treatments with Auger electrons and alpha particles (25, 26). The efficiency of beta particles for treating small tumors depends on their energy and penetration depth into the tissue. The range of emitted beta particles should correspond to the size of the tumor, ensuring that minimal radiation dose is delivered to the adjacent healthy tissue (10). Various radionuclides emit both short-range, non-penetrating radiation and longer-range, penetrating radiation, serving the dual purpose of therapy and subsequent imaging to support patient-specific dosimetry. A more extended physical half-life is desirable for prolonged retention, facilitating the cumulative irradiation of the tumor and maintaining a high ratio of penetrating to non-penetrating radiation. The physical half-life of the radionuclide needs to align with the biological half-life of the radiopharmaceutical to ensure effective treatment (27).
Numerous radionuclides that possess the above-mentioned therapeutic decay characteristics are shown in Table 1 and have been documented extensively (24). 131Iodine continues to play a significant role as a therapeutic radionuclide, particularly in thyroid gland ablation using [131I]-NaI (28). The advancements in RPT, primarily with the introduction of radiolabeled antibodies and molecules like somatostatin analogs and ligands, have resulted in the widespread use of [177Lu]Lutetium-PSMA for treating castrate-resistant prostate cancer (CRPC) (29), as well as [223Ra]Radium-Chloride (30). Additionally, 90Yttrium microspheres have found application in the treatment of both primary and metastatic liver cancer. Given the high biological effectiveness of alpha emitters, there is a growing interest in RPT utilizing radionuclides such as 227Thorium, 212Lead/212Bismut, 211Astatine, 213Bismuth, and 225Actinium (31). Efforts to enhance the specific activity of radionuclides, such as 153Samarium, have been explored for potential applications in RPT for metastatic bone palliation (32, 33). There have only been two clinical trials thus far that summarize the clinical circumstances surrounding the less commonly utilized 211Astatine with antibodies as one of the carriers (34). These investigations involve, amongst others, colon cancer, glioma leukemic, and neuroblastoma.
[177Lu]Lutetium-DOTA-TATE/NOC/TOC has been used in RPT applications to treat late-stage neuroendocrine tumors (NETs) with extensive metastases while [90Y]Yttrium-DOTATOC is less frequently employed (40–42). This choice is attributed to 177Lutetium's favorable decay characteristics and compatibility with peptide labeling, making it an ideal candidate for the treatment of NETs. RPTs using 177Lutetium-based radioligands have the advantage that they can be applied to primary and metastatic cancers with relatively low toxicity (14).
177Lutetium has been proven as an established RPT radionuclide due to its application in treating patients with NETs and CRPC. Even though several isotopes are accessible for use in radionuclide treatment applications, 177Lutetium persists as the preferred radionuclide used routinely in clinics (43–45). It also has the additional benefit of being readily available commercially. As a result, our emphasis on activity quantification and dosimetry will focus on 177Lutetium.
1.1.1 177Lutetium decay characteristics
With a physical half-life of 6.7 days and three excited levels above the ground state, 177Lutetium decays by beta emission to the stable ground state of 177Hafnium (depicted in Figure 3). The physical half-life is sufficiently long to ensure that transportation, storage, and delivery are manageable, yet it is short enough to spare OAR from excessively high doses (46). The highest energy (Emax) that may be produced through Beta emission is 498.3 keV (79.3%) (47). 12.0% of the emissions involve Beta particles with an Emax of 177.0 keV, which yields an excited state of 177Hafnium at 321.3 keV above the ground state. Furthermore, Beta particles with an Emax of 385.4 keV (9.1%) are encountered, leading to an excited state of 177Hafnium at 249.7 keV above the ground state. The decay process produces trace amounts of Auger electrons and x-rays. 177Hafnium de-excites to the ground state by emission of gamma rays, with energies most prevalent at 112.9 keV (6.2%) and 208.4 keV (10.4%). At energies of 54.6 keV (1.6%) and 55.8 keV (2.8%), characteristic x-rays are obtained (47). Owing to the Beta particles' interaction with the tissue, the Bremsstrahlung yield is extremely low (>11%), with the majority of the photons (∼85%) having energies below 50 keV (48).
Figure 3. A Beta decay scheme depicting the emission process of 177Lutetium (47).
The two above-mentioned gamma rays have been effectively employed for imaging (49) to facilitate dosimetry (50). The Beta particles (Emax of 498.3 keV) are successfully applied in RPT (51) and have a short mean soft tissue penetration depth of 0.7 mm (52). The Beta particles from 177Lutetium decay ensure a uniform absorbed dose at the cellular level due to the high energy while still avoiding harm to the surrounding healthy tissues (53).
1.2 Theranostics with 177Lutetium
177Lutetium is considered a theranostic radionuclide since it can be used to assess tumor uptake and cancer progression and for treatment purposes (54). The theranostic approach with 177Lutetium employed for personalized patient management may be direct or indirect. The direct method takes advantage of 177Lutetium's emission of both gamma and beta radiation as demonstrated by Figure 4. This gamma radiation can be employed for image-based treatment planning before sequential treatment cycles and for RPT. This method employs image-based dosimetry from the preceding treatment to predict potential radiation-induced adverse side effects and to ensure that the doses received by the OAR remain within acceptable tolerance levels. This approach determines the feasibility of administering additional dose cycles following the initial 177Lutetium therapy. The use of image-based indirect theranostic approaches with [68Ga]Gallium-PSMA or [68Ga]Gallium-DOTA-TATE/NOC/TOC is employed to diagnose patients with CRPC and NETs, eligible for [177Lu]Lutetium-PSMA and [177Lu]Lutetium-DOTA-TATE/NOC/TOC treatment respectively.
Figure 4. Illustration of the theranostic direct approach using [177Lu]Lutetium-PSMA (55).
1.2.1 177Lutetium dose-limiting factors and use in neuroendocrine tumors and prostate cancer
1.2.1.1 [177Lu]Lutetium-somatostatin receptors
The use of [177Lu]Lutetium-DOTATATE for the treatment of patients with midgut NETs was approved as a result of the NETTER-1 trial. The dose-limiting OAR for RPT of [177Lu]Lutetium- DOTATATE in patients with NETs are the kidneys and bone marrow (56). The activity and the number of cycles administered for RPT are governed by the tolerance values of the OAR. The decision to administer a subsequent cycle relies on the estimated absorbed dose determined from the dosimetry data of the preceding cycle, to ensure that the cumulated dose does not exceed the tolerance levels of the OAR (57). Furthermore, randomized clinical trials depend on the threshold tolerance values of the dose-limiting organs applicable to RPT. Since acute hematological toxicity is common, the administration of additional cycles during treatment may need to consider bone marrow deficiency. The permissible maximum absorbed dose to the bone marrow is 2 Gy (58). Although bone marrow involvement is taken into account, it is not the primary organ that limits the dose administered (59). The most critical OAR for [177Lu]Lutetium-DOTATATE RPT is the kidneys. This is owing to the proximal tubular reabsorption and retention of the radiopeptide in the interstitium and subsequently renal irradiation. Positively charged amino acid infusion may reduce kidney doses to some extent by lowering high renal retention and limiting proximal re-absorption (60). Since the kidneys are a late-responding tissue and the radiation-induced harm from deterministic effects may not become apparent for one year or more, it is critical to monitor kidney dosage and kidney function in patients. This finding prompted a research cascade in renal dosimetry of [177Lu]Lutetium RPT (61–67).
It has been commonly accepted that kidney doses obtained via RPT should not exceed the maximum tolerance threshold of 23 Gy (68). This value is based on the fractionated EBRT tolerance (69) which results in a 5% probability of developing late kidney damage in 5 years. Garkavij et al. (56) contended that the suggestion might be called into doubt due to the notable differences in the absorbed dose rate between EBRT and RPT. For this reason, the application of dose-response data applied to RPT has been reviewed (70–72). The tolerance was raised to 27 Gy in research by Valkema et al. (73), citing the fact that EBRT provides a dose rate that is significantly higher than RPT. Clinical trials have been carried out where additional treatment cycles are administered until the calculated biological effective dose (BED) for the kidneys reaches either 27 Gy or 40 Gy, with the choice between these levels determined by patient-specific risk factors (74). Despite these debates, the generally accepted upper limit of kidney dose that can be administered with RPT remains 23 Gy. Evidence-based trials would need to demonstrate that patient-individualized RPT with [177Lu]Lutetium-DOTATATE is superior to the 7.4 GBq over the four-cycle standard as indicated in the NETTER trial (75).
1.2.1.2 [177Lu]Lutetium-PSMA-617
The use of [177Lu]Lutetium-PSMA-617 for the treatment of patients with metastatic CRPC has been approved by the FDA as PluvictoTM after the results of the VISION trial (76). The dose-limiting organs for [177Lu]Lutetium-PSMA-617 are similar to those of [177Lu]Lutetium-DOTATATE with a potential concern of xerostomia radiation-induced effects on the salivary glands. Similarly taken from EBRT data, the tolerance threshold of the salivary gland is based on EBRT data as 25 Gy to ensure a less than 25% probability of long-term radiation damage. Salivary gland dysfunction varies amongst patients and dose reduction measures such as hydration can mitigate the radiation-induced effects (77). Radiation-induced salivary gland toxicity has been found to have minor clinical relevance for [177Lu]Lutetium-PSMA (78).
The response variation of patients treated with 7.4 GBq of [177Lu]Lutetium has been ascribed to the comparatively low activity of “one dose fits all” (79). For this reason, RPT finds itself at a crossroads between patient-specific treatment and the rigid “one dose fits all” fixed activity regime. Table 2 outlines the vectors explored with the [177Lu]Lutetium treatment. Among these vectors [177Lu]Lutetium-PSMA-617 and [177Lu]Lutetium-DOTATATE are extensively used.
1.2.2 [225Ac]Actinium-PSMA-617 referencing [177Lu]Lutetium-PSMA-617
Notably, in instances where therapy with the beta-emitting [177Lu]Lutetium-PSMA-617 has failed, RPT targeting PSMA with the alpha-emitting [225Ac]Actinium-PSMA-617 has demonstrated therapeutic success (93). This has the benefit of the higher LET from the alpha particles, that if the radionuclide is coupled to a vector that targets antigens expressed on tumor cells, the radiation from the alpha particles can be efficiently delivered to the malignant cells. In this manner, the malignant cells will receive high energy delivery with minimal radiation harm to the surrounding healthy cells, which has been beneficial in the treatment of prostate cancer (94–96). To achieve reliable dosimetry data, the dose deposition should be accurately described to the range of the dose-depositing particles.
The emission range of the RPT particles varies, from millimeters and micrometers to nanometers for beta, alpha, and Auger electrons respectively. Dosimetry computations may be cumbersome with alpha particles due to their limited tissue range and high LETs, resulting in a high-energy deposition close to the emission site. Since dosimetry is computed from positron emission tomography (PET) and single=photon emission computed tomography (SPECT) images the adaption to RPT with alpha-emitting radionuclides administered with low activities leads to low signal-to-noise ratios from the gamma rays, causing image quality challenges (97, 98). These challenges are also brought on by the decay properties of alpha emitters, such as 225Ac, with limited gamma emission and the competing bremsstrahlung radiation (99).
Previous dosimetry on [225Ac]Actinium-PSMA-617 and [225Ac]Actinium-PSMA-I&T extrapolated the uptake of respective 177Lutetium-labeled analogs on imaging to mitigate the imaging challenges with alpha emitters (99–101). The derivation of time activity curves (TAC) from [177Lu]Lutetium-PSMA-617 data extrapolated to the physical half-life of [225Ac]Actinium-PSMA-617 by assuming instantaneous decay of daughter nuclides results in uncertainty of bio-kinetic data at the cellular level (102). Alternatively, the dosimetry of [213Bi]Bismut-PSMA-617 was previously determined by extrapolating the uptake of [68Ga]Gallium-PSMA-617 from PET images (103). Recent investigations have demonstrated that the degree of uptake with radiolabeled PSMA-617 varies depending on the radionuclide (104, 105). This suggests that an independent assessment of the uptake of [225Ac]Actinium-PSMA-617 would improve its dosimetry.
Because alpha particles have short ranges in tissue, small-scale dosimetry techniques such as microdosimetry and autoradiography have been identified to determine the dose distribution on sub-organ levels (97, 106–110). To date, the integration of small-scale dosimetry into clinical practice has not been achieved (14, 15). This is because the activity distribution needs to be quantified to the cellular level for the range of the alpha particles. Even with voxel-based dosimetry, source-target combinations on this scale are challenging to measure, and the distribution of activity over time needs to be quantified at the subcellular level (1). In clinical practice, the quantitative information used as the input data for dosimetry is conventionally obtained from PET and SPECT scans, which have spatial resolution in the order of millimeters. For these reasons, preclinical validations have benefited more from sub-organ dosimetry (107, 111, 112).
1.3 Clinical intent of dosimetry
The approach to RPT should be tailored to match the intended therapeutic aim, which could be curative or palliative (113). For a treatment intended to increase the chances of a cure, the focus would be on delivering a substantial therapeutic dose within a relatively short timeframe. The biological response to RPT is determined by the absorbed dose and the LET (113, 114). Individualized administered activity in RPT has been referred to as a “quest for the holy gray” because there are no threshold values from RPT and dosimetry is based on EBRT threshold values (12). Although RPT dosimetry has undergone tremendous effort, concerns such as “to what extent is dosimetry needed for RPT applications” arise from the threshold uncertainty (115). In this exciting era of innovative RPT agents with the possibility of precision medicine, it is still up to the NM community to optimize and customize patient treatment. Initiating patient-specific dosimetry for personalized administered activity would pave the ground for this endeavor considering the therapeutic objective of RPT, whether it is for palliative or curative purposes.
Image quantification of activity distributions and absorbed dose modeling are two aspects of dosimetry. The basic requirements of dosimetry are to determine the radiopharmaceutical bio-distribution and calculate the absorbed dose to the organs of interest based on this distribution. There are numerous processes involved in the procedures used to obtain the bio-distribution data, and the assumptions made vary greatly depending on the protocols and radiopharmaceuticals. Some of the factors contributing to the relatively low acceptance of dosimetry include the large variety in absorbed dose techniques and the predicted “tolerated holy gray.” A foundation for standardizing dosimetry methods is provided by an understanding of the importance and constraints of each step in the clinical dosimetry chain.
2 Quantification of image-based activity for dosimetry
The quantities for image-based dosimetry include a) the number of decays in the organs of interest for the various times post-administration, b) the decay particles' temporal pattern which determines how the energy is released from the time-activity data, and c) the absorbed energy in the respective mass or volume of the organ of interest. Since the absorbed dose to the radiation-exposed organs is determined from the image-based activity distributions, it pertinently follows that the accuracy of the activity quantification is a crucial factor for dosimetry. A series of steps is required for the calculation of absorbed doses. These steps form part of a chain, illustrated in Figure 5. The process begins with (i) measuring activity in the dose calibrator and obtaining a gamma camera calibration factor, followed by (ii) acquiring patient data through a designated imaging protocol. In the case of SPECT/CT acquisitions, the activity distribution is quantified from the reconstructed corrected registered images which are (iii) segmented using a volume of interest definition to derive absolute activity quantitative data from organs of interest by applying a previously determined calibration factor. Subsequently, (iv) a time-activity curve (TAC) is generated, and TAC analysis and integration are performed to produce a time-integrated activity (TIA). Finally, (v) absorbed dose computations are conducted.
Every step propagates the dosimetry inaccuracy and precision. The foundation for achieving accurate dosimetry lies in the accurate quantification of activity, which is dependent on how faithfully the images represent the true activity distribution. SPECT/CT has overcome the challenges associated with planar imaging. When SPECT/CT modalities are used, the corrections for image degrading factors are carried out as part of the clinical workflow. Therefore, the discussion on image acquisitions shall be made with reference to SPECT/CT data.
2.1 Gamma camera calibration factor
SPECT images are generally considered to be non-quantitative, in contrast to PET images. For images to be used for quantification and dosimetry purposes, SPECT images must be obtained in units of activity concentration (kBq/cm3) rather than the conventional counts. For certain gamma cameras, to obtain SPECT images in absolute units of activity concentration, a calibration factor (CF) must still be applied to the reconstructed image. The system sensitivity is frequently used as a CF in SPECT images to gain quantitative data (116). It is determined from a planar image of a Petri dish in the air with known activity. The accuracy of the dose calibrator used to quantify the activity for the CF also determines the effective activity administered to the patient, which is another crucial factor to take into account (117). To increase the accuracy of CFs, a variety of geometries have been explored (118).
For activity quantification, the process of obtaining the CF ought to approximate the method used in the quantified clinical studies. The CF source should provide reconstruction and compensation techniques for SPECT applications. This enhances quantitative accuracy and lessens the impact of inaccurate scatter and attenuation corrections (119). The CF measurement is a function of the source geometry and depends on how the volume used to obtain the counts for the CF is defined and incorporated with the recovery coefficient (RC) (120). Techniques for image-based CF that employ SPECT/CT and corrected planar patient data have been presented and compared with conventional phantom CF (121). For the patient data CF, the authors discovered an error rate of less than 2%, resulting in an overall quantification accuracy of 7%.
Establishing a standard for absolute quantitative SPECT by harmonizing the CF in a multicenter and multivendor setting where reconstruction techniques are established is a key step to standardizing dosimetry. Such a multicentre and multivendor study demonstrated that a high patient body mass index (BMI ≥ 47 kg/m2) increased CF variability between systems and made it more difficult to quantify minor lesions (less than 10 mm3) (122). A standard CF, for gamma camera vendors and models, can be used as input for absolute SPECT quantification in radionuclide therapy research; this will aid with complex clinical dosimetry, especially in multi-center research endeavors. To support this, image voxels from reconstructed SPECT images have recently become available in radioactive concentration units in recent gamma cameras such as the Siemens Symbia IntevoTM scanner (Siemens Healthineers). Concerning PET imaging, the PET systems are usually calibrated to measure accurate concentrations of 18Fluorine. The reconstructed standard uptake value (SUV) should be checked for the radionuclide in question (123). For radionuclides other than 18Fluorine the reconstruction should incorporate the physical characteristics of the radionuclides to obtain accurate quantitative images.
2.2 Image acquisition protocol
Patient dosimetry is based on gamma camera images that have been tainted by uncertainties in the imaging procedure and the related protocols. Figure 6 depicts various physical factors that degrade the gamma camera images from a perfect representation of the activity distribution within the body, as shown in Figure 6A. These factors include Photon attenuation (absorption) and Compton scattering (scatter), the effects of which are demonstrated in Figure 6B. Figure 6C illustrates the effects of collimator-detector response (CDR) or collimator blurring, while Figure 6D demonstrates partial volume effects (PVEs) and noise. The influence of these variables on patient dosimetry is contingent upon the gamma camera's chosen imaging protocol and the corrections applied to account for these degradation factors. For this reason, image activity quantification has been a subject of investigation for many years (124–132). Over the years, there has been consistent progress in improving the accuracy of activity quantification and image analysis (120, 133–138).
Figure 6. Illustration of degrading physical factors inherent with gamma camera imaging using Monte Carlo simulations corresponding to (A) image simulated with perfect gamma camera resolution. (B) Image including degradation from photon attenuation and scatter. (C) Image obtained in “b” including normal gamma camera resolution. (D) The image obtained in “C” includes partial volume effects and clinically realistic noise levels (127).
Selecting the appropriate collimators is the first step toward optimizing the imaging protocol. This is based on the energy of the gamma-ray to be imaged and the trade-off between sensitivity and spatial resolution. In the case of 177Lutetium, the imaging protocol, the collimator, the imaged gamma-ray energy, and energy window settings have been thoroughly investigated using Monte Carlo (MC) simulations (127). The 208 keV photopeak should be used with a medium energy (ME) collimator and a 20% energy window setting (49, 56, 67, 118, 127, 139–144). Due to the low gamma-ray emission of 177Lutetium, the effects of dead time are minimal and imaging can commence immediately after therapeutic activity values are administered (145).
2.3 SPECT image activity quantification
Traditionally, 2D planar anterior and posterior WB images have been used for activity image quantification in dosimetry (146–149). Planar WB methods offer a quick and simple method to image the patient's entire body and extract biodistribution information. However, organ overlap, superimposed background activity, and the absence of organ volume information are known drawbacks of planar imaging. When compared to planar images, SPECT images have shown superior image quality and improved quantitative accuracy (150–152). Planar images are often considered to have poorer quantitative accuracy since not all centers consistently undertake additional image acquisitions such as blank and transmission scans. These additional scans are used to compensate for attenuation and further analysis of the images in sub-windows are necessary for scatter corrections. Hybrid planar WB SPECT/CT imaging techniques, also known as hybrid WB/SPECT images, have resolved the trade-off between the faster multiple biokinetic image data collections for dosimetry that planar imaging advocates and the enhanced accuracy provided by SPECT/CT images (153–155). Even though SPECT/CT has overcome the challenges associated with planar imaging, the reconstruction algorithm's integrity and the impact of the inherent image acquisition degradation factors complicate SPECT imaging. Therefore, several steps must be taken to increase the SPECT quantitative accuracy, which is covered in more detail in the following sections.
2.3.1 SPECT image reconstruction
Iterative reconstruction algorithms such as the maximum a posterior (MAP) (156), maximum likelihood expectation maximization (ML-EM) (157), and ordered subset expectation maximization (OS-EM) (158) can include modeling of the physical characteristics of the imaging process. These algorithms mainly consist of compensations for collimator and object scatter, system geometry, and finite detector resolution. These algorithms result in reconstructed images with better image quality and quantitative accuracy and are less prone to artifacts compared to analytical methods such as filtered back projection (159). The OS-EM algorithm has become a standard algorithm with most clinical SPECT processing units. It is highly recommended and commonly used to obtain improved quantitative SPECT data (119, 160). When using the OS-EM reconstruction algorithm, it is important to consider the optimum number of updates, defined as the product of the number of subsets and iterations for a particular SPECT system's reconstruction algorithm. The trade-off is that more updates lead to higher levels of image noise but yield more accurate quantitative activity distributions from the reconstructed images (127, 161, 162). To improve image quality and activity quantification accuracy, the current 3D OS-EM reconstruction algorithms available with most SPECT systems include, in addition to attenuation and scatter corrections, compensation for CDR.
Reconstructing an image from raw projection data is an inverse problem. Although artificial intelligence (AI) technology, particularly deep learning-based solutions, has emerged as a promising solution to the reconstruction of emission images, the inverse problem remains unsolvable with AI. Three distinct systems, namely static scan, dynamic scan, and hybrid fusion, have realized the introduction of AI applications centered around NM image reconstruction (163). AI essentially provides a mapping connection to address specific significant reconstructive challenges, such as completing the transition from the sinogram domain to the image domain or substituting for regularization in traditional algorithms with a data-driven approach. The majority of the image reconstruction work in NM using AI technology is in PET reconstruction (164).
2.3.2 Attenuation correction
Attenuation is the most significant factor that reduces quantitative accuracy, particularly in large patients (165). The routine implementation of attenuation correction has been made easier by the spatially and temporally co-registered PET/CT and SPECT/CT data. Accurate SPECT (or PET) and computed tomography (CT) data registration is necessary to maintain the integrity of the attenuation correction from the CT data (166). Converting the CT images into acceptable attenuation maps using linear attenuation coefficients is the first step in attenuation correction. This entails employing bi-linear models (167), as shown in the example for 177Lutetium depicted in Figure 7, to map the effective CT energy to the radionuclide's primary emission energy. Increased CT numbers in patients with metallic implants and those who have taken contrast agents lead to inaccurate SPECT attenuation coefficients, which in turn causes overestimation of radioactive uptake and false-positive results (133).
Figure 7. An example of a bi-linear model used to map the effective CT energy to the imaging emission energy of the radionuclide of interest.
Advances in AI technology have made it possible to make improvements in attenuation and scatter correction for PET and SPECT scans. AI enhances PET/magnetic resonance imaging (MRI), as well as PET-only and SPECT-only gamma cameras, by facilitating the use of synthetic attenuation maps derived from uncorrected emission images. This eliminates the necessity for a CT scan for attenuation and scatter correction. Accurate attenuation correction maps produced by AI for myocardial perfusion SPECT imaging have been demonstrated (168). Moreover, AI-based attenuation correction prevents misregistration artifacts and simplifies identifying artifacts for physicians (169). Attenuation correction using deep learning-based methods has shown promising results in quantitative SPECT imaging (169–171).
2.3.3 Scatter correction
It is assumed that a photon is eliminated if it is completely absorbed or scattered by the linear attenuation coefficients that are produced using CT data for corrections (172). To increase the quantitative accuracy, scatter correction seeks to eliminate scattered gamma rays that have occurred away from the emission location. Multiple window scatter compensation approaches have practical applications, therefore most gamma camera vendors conduct scatter correction for SPECT data using methods such as the dual-energy window (DEW) (173) and the triple-energy window (TEW) (174). Down-scatter from higher energy gamma rays than those for which the imaging window was set is considered by the TEW, while the DEW only considers the self-scatter from the source. The process of scatter correction involves deducting the predicted scatter projections from the primary photopeak projection image, on a pixel-by-pixel basis. The photopeak and sub-windows record distinct spatial distributions of the scattered gamma rays, and the gamma rays in the sub-windows do not undergo the same scattering as those in the photopeak windows. Negative counts and noise amplification are possible outcomes of the subtraction procedure (175).
In the reconstruction process, scatter may be modeled using the effective source scatter estimation (ESSE) approach (176, 177). To estimate the scatter contribution in the photopeak, the ESSE makes use of pre-calculated scatter kernels of a point source at a specified distance from the collimator face in a uniform water-filled slab phantom based on MC simulations. The 3D OS-EM iterative reconstruction may incorporate the ESSE scatter kernels to account for the scatter. According to reports, the ESSE produces quantification accuracies that are superior to those of the TEW and DEW scatter correction approaches (139). Although practical, multiple energy window methods for scatter correction involve subtracting counts which may lead to reduced image sensitivity. AI has emerged as a promising solution for scatter correction, where the scatter sinogram can be generated from the raw data of emission and attenuation data obtained from PET or SPECT, or it can be created using uncorrected PET images as input data (178–180). AI promises to improve both patient throughput and image reconstruction speed when it is applied to scatter correction (181).
2.3.4 Collimator detector response
Source-to-detector distance determines the CDR, which is the primary factor affecting the spatial picture resolution in SPECT images, as depicted in Figure 8 (119). Gamma rays pass through and interact with the collimator and detector of the gamma camera, resulting in three factors that alter the CDR (182). The collimator resolution, which is based on the geometrical acceptance angle of the collimator holes and decreases linearly with the source-to-collimator distance, is the first factor (172). The second factor is the detector's intrinsic resolution. To precisely estimate gamma ray interaction positions, this component is restricted by the uncertainty present in the crystal and gamma camera positioning electronics. The probability that the gamma rays will penetrate the septa makes up the third factor.
Figure 8. Illustration of the SPECT collimator-detector response determined by the source-detector distance. (A) Narrow point spread function closer to the detector. (B) A broader point spread function extends further from the detector.
When gamma rays scatter in the collimator septa and are detected in the main energy window, this is known as septal scatter. Each of these components exacerbates the blurring of the image and reduces the spatial resolution. The effects of medium- and high-energy gamma rays are more noticeable in terms of dispersion and septal penetration. These factors reduce the quantitative accuracy by affecting the resolution of SPECT-reconstructed images of individual pixels and the capacity to identify microscopic diseased tissue, such as metastasis.
The geometric response, which consists of the first two components mentioned above, is easily compensated for by most commercial systems as part of the CDR during the iterative reconstruction process. The manufacturer sets the specifications of the modeling. Increased quantification accuracy has been reported for 177Lutetium spheres positioned in a phantom, employing MC modeling scatter compensation techniques similar to the convolution-based forced detection scatter correction which incorporates collimator and detector modeling (126).
CDR correction partially mitigates activity spill-out by reducing the poor resolution effects that cause PVEs. It is important to note that the Gibbs ringing artifacts could appear in the vicinity of sharp boundaries of activity distributions while employing CDR (183). Even while CDR compensation helps minimize PVEs, partial volume corrections (PVCs) still need to be applied to the reconstructed images near the gamma camera resolution limit to improve the accuracy of activity quantification (184).
2.3.5 Partial volume effect
The PVE refers to the phenomenon when activity concentration from an emission reconstructed image is not only confined in the respective voxel but also smeared out into the neighboring voxels. This occurs when an imaged object is similar in size to a multiple of the gamma camera's estimated spatial resolution assessed in terms of its full width at half maximum (FWHM), that is, less than 2 × FWHM (173). Image blurring of objects below the aforementioned resolution volumes is a result of the gamma camera's limited spatial resolution and image sampling in an image matrix, causing the PVE.
The measured activity distribution is said to be underestimated by the PVE, illustrated in Figure 9A, which is commonly observed in tumor imaging (185). This is one of the challenges in determining dose-response relationships in RPT for tissues other than OAR, which directly impacts dosimetry (186).
Figure 9. (A) Profiles drawn through reconstructed spheres of varying sizes above and below the gamma camera resolution limit, with the spread functions illustrating the activity underestimation. (B) Illustration of a recovery coefficient used for partial volume correction (120, 187).
The RC is one of the PVC techniques describing the ratio between the actual and measured activity concentration (188). For well-defined shapes, often spheres, an RC curve can be produced as a function of object size (Figure 9B). In phantom investigations where the true activity and object size can be determined, this can be used to readily generate characterization of the gamma camera's PVEs for particular shapes. Nonetheless, this compensatory technique might provide challenges in clinical investigations when the sizes of the organs and tumors are unknown and have irregular shapes. Higher-resolution modalities that include CT or MRI examinations may be employed to determine the size of the tumors and organs. However, since every clinical study would be different, modeling partial volume for all irregular geometries would be cumbersome.
Dosimetry evaluated with a SPECT clinical example ranged from >99% underestimation in the smallest lesion (4 × 5 mm) to more than 60% underestimation in the greatest lesion (28 × 22 mm) (189). Several PVC quantification techniques were examined for 177Lutetium, in a 3D-built kidney phantom featuring cortex and medulla compartments (190). The authors recommended against using a sphere-based RC for characterizing organs like the kidney and instead using a geometry-specific RC. Large variations in PVE magnitude are caused by variances in kidney shape. Moreover, the distribution of intrarenal activity has a significant impact on the severity of PVE.
Consequently, it is highly unlikely that RCs created from simpler phantoms will be adequate to rectify the PVE in patient images. Furthermore, renal RCs have been reported to be well-modeled by the surface area-to-volume ratio; this method may also be used for other geometries (191). Numerous software-based PVCs have been proposed (119) but, due to their complexity, they have not been widely used in clinical settings. This explains why sphere-based RCs for PVC are still commonly used, particularly in the measurement of tumors. When assessing radionuclide uptake in vivo, the PVE is still a significant component, particularly in small volumes. It will be especially challenging to estimate absorbed doses from RPT accurately until robust methods to account for the PVE are developed.
2.3.6 Volume of interest definitions
To achieve accurate activity quantification the VOI definition is another crucial factor to take into account (192, 193). There isn't a widely recognized technique for defining VOI in NM images. A common approach has been to employ anatomical CT data from SPECT/CT images, a method that has been extensively used in clinical and phantom investigations (128, 137, 152, 190, 194). SPECT quantification accuracy is affected by errors such as misdefinition (variability in organ delineation) and misregistration (between emission and transmission data) errors. Notably, even a single voxel misdefinition can lead to a quantification error of up to 8% (133). These errors become more pronounced in small organs with low levels of activity uptake. While a fully automated segmentation method using CT images and convolutional neural networks has shown accelerated organ segmentation and high accuracy in kidney dosimetry for 177Lutetium, RPT, it still necessitates expert supervision and corrections, primarily due to misalignments in the co-registration of SPECT and CT images. Image sampling affects the variability in the VOI definition, and the precision of the quantification accuracy may be improved when using smaller pixels (195). As new radiotracers are consistently introduced into NM dosimetry, the utilization of automated AI-based segmentation is increasingly seen as an advantageous initial step in the dosimetry process (196).
Compensating for the aforementioned factors contributes to an improved estimation of the activity distribution. The selection of appropriate correction techniques depends on various factors, including the specific clinical study under investigation, the accessibility of the correction methods, and the radionuclide under consideration. To assess the accuracy of activity quantification, the percent difference between the actual activity distribution and the quantified activity distribution obtained from the acquired images is estimated. The accuracy of dosimetry is directly linked to the accuracy of the quantified activity.
3 Radiopharmaceutical dosimetry
The physical measure used to assess the effects of ionizing radiation in tissue is the absorbed dose (27), which is often used to characterize the energy delivered in a volume with a specific mass. Gray ((J)⁄(kg)) is the international system unit used to express absorbed dose. From this, it follows that the definition of the mean absorbed dose is given by Equation (1).
Where is the mean energy imparted and is the mass of a specific tissue volume. Equation (1) is extended to account for the different source-target relationships and the radionuclide decay characteristics expressed in Equation (2).
Where is the time-integrated activity (TIA), which depends on the bio-kinetics of the activity distribution and represents the cumulated activity for a specific period, while is the mean absorbed dose deposited in the target organ per TIA unit present in the source organ (Figure 10). Equation (2) forms the foundational basis for RPT dosimetry. This basis was initially developed and presented in the 1960s by the Medical Internal Radiation Dose (MIRD) Committee of the Society of Nuclear Medicine, known as the MIRD formalism. The MIRD formalism gained wide acceptance as the standard method for absorbed dose calculations (197). The absorbed dose computation is the product of two quantities, the imparted energy summations over all emission types (S-coefficient) and the cumulative activity (TIA), which will be elaborated upon in the subsequent discussions.
Figure 10. The absorbed radiation dose directed to target organs, symbolized by healthy tissue, originating from a source organ, which is depicted as cancerous cells in the liver containing radioactivity.
3.1 Time activity data
3.1.1 SPECT/CT images for time activity curves
SPECT/CT images have been widely recognized for their improved capability to estimate activity distributions (124, 129, 198). In the process of activity quantification, the count rate per voxel in the reconstructed SPECT image is proportional to the activity concentration. This proportional relationship is achieved given that previously mentioned corrections are applied (199), encompassing scatter correction using the TEW method (119, 200), CT-based attenuation correction, CDR compensation (119, 201), and corrections for PVEs using RCs (133). The corrected count rate per voxel is then divided by the calibration factor [(counts per second per voxel)/(kilo becquerels per milliliter)] to determine the activity concentration. Even with improved accuracy, the repeated acquisition of multiple SPECT/CT images for dosimetry can be time-consuming for clinical facilities and can impact patient throughput. Consequently, caution should be exercised in the selection of optimal image sampling schedules, as suboptimal choices may lead to significant over- or underestimations of absorbed dose estimates, particularly for organs at risk such as the kidneys (62). It is crucial to understand the purpose of dosimetry and the resources available in a clinical setting to achieve accurate results. While simplifications in dosimetry methodology may affect its accuracy, an acceptable level of uncertainty can be determined based on clinical needs. Personalized dosimetry often requires SPECT/CT studies at multiple time points, but efforts have been made to reduce the burden on patients and clinics. In centers with limited capacity, dosimetry can be performed at alternate cycles or by using less quantified images for subsequent cycles. The use of only the initial cycle, as well as post-therapy imaging, can be developed into a quantitative image for absorbed dose estimations. Gamma camera availability limitations can be addressed by using a hybrid approach of WB/SPECT. Pre-therapy images can predict therapeutic absorbed doses in theranostic applications, allowing for tailored activity prescriptions for optimized therapy success. Hybrid WB/SPECT imaging methods have been introduced to streamline the imaging process for dosimetry (153). Furthermore, the use of single time-point post-treatment imaging for dosimetry, based on SPECT/CT data, has recently been described in the literature; however, its widespread adoption has yet to be achieved (202).
3.2 Time activity curve
The TAC gives temporal information regarding the patient-specific variations in radiopharmaceutical uptake, retention, and excretion. Successive quantitative images obtained at various time points after the administration of radiopharmaceuticals are used to generate TACs. If the TAC is not determined optimally, it could lead to considerable dosimetry inaccuracies. Once the TACs for various source organs of interest are obtained, the TIA is determined by integrating these curves. The accuracy of the TAC is influenced by the number of imaging time points, the frequency of the image sampling schedule, the appropriate imaging span period following administration (integration period), and the model applied for TAC fitting, especially for TIA calculations (203–205). For accurate results, at least three data points should be gathered throughout two to three effective half-lives, with the integration time matching the study's biological endpoint (27). Gleisner et al. (206) observed detectable levels of [177Lu]Lutetium-DOTATATE five to seven weeks after injection due to tumor retention. Tumor dosimetry at a time point beyond the conventional seven days following administration might be useful.
3.2.1 Time-integrated activity
The TIA, also referred to as cumulated activity, is computed by taking the area under the TAC derived from the series of imaging time points. Multi-exponential functions that are integrated analytically can be fitted to the TACs. Uncertainties in the TAC fitting models have been reported (207). The constants determined from the exponential functions describe the bio-kinetic data in the organ of interest. Statistical analyses have been conducted to determine the best-fitting functions for the TIA (203). The total number of decays (S-coefficient) obtained within a specific volume must be obtained to calculate the absorbed dose once the following conditions have been met: (i) The quantified activity has been accurately determined, (ii) the optimal image sampling has been obtained, and (iii) the optimum fitting model to compute the TIA has been established. Performing two-time point imaging can yield TIA estimates with an average error below 5% of the reference TIA for both tumors and kidneys. Similarly, three-time-point imaging offers a comparable level of error but exhibits less variability (208).
Optimal sampling of hybrid WB/SPECT images for kidney dosimetry has revealed that the variability in TIA depends on the number of post-administration activity images, especially in scenarios with two-time points (166). In such cases, prior knowledge of population averages for biokinetic data is necessary, and this approach may not be suitable if patients' biokinetics deviate significantly from the population average, particularly in the context of toxicity detection. The application of a single time point may deviate from the principle of patient-specific dosimetry for treatment planning.
3.3 Absorbed dose computations
The S-coefficient accounts for the mean energy emitted per decay of the radionuclide, the mass of the target organ, and the absorbed fraction of the energy emitted from the source organ. The absorbed fraction depends on the type and energy of the emitted radiation, the size, shape, and composition of the source and target regions, and the distance and type of material separating the source-target regions. The S-coefficients are stored as look-up tables for a variety of radionuclides and multiple source-target organ combinations (209). The MIRD schema is commonly used for organ-level S-coefficient estimates (210). The dosimetry models employed for the computation of S-coefficients have advanced from simple geometric shapes to more intricate voxel-based phantom series using NURBS models, as illustrated in Figure 11 (211, 212). These models rely on standardized published organ masses (213, 214). More accurate patient-specific dosimetry can be achieved if the S-coefficients are scaled by the mass of the patient's organ. The progression of these models, encompassing various source and target organs, has been well-documented (215).
Figure 11. Original stylized adult male model of (212): (A) exterior view, (B) skeleton and internal organs. Anterior views of the Segars NURBS models (216) (C) adult male and (D) adult female.
Three approaches are available to compute the S-coefficients, namely, local energy deposition (LED), convolution using dose point kernels (DPK), and direct (full) Monte Carlo transport. A selection criterion, detailed in reference (102) and illustrated in Figure 12, may serve as a practical guideline for selecting the most suitable algorithm for absorbed dose calculations. To elaborate on the guidelines, if the radiation in question is non-penetrating, LED may be an appropriate choice. However, for penetrating radiation with a high emission yield, where energy deposition from radiation sources outside the original volume could contribute significantly to cross-doses, the other two algorithms should be considered. In instances where radiation propagation within tissue is homogeneous, the behavior of LED is primarily influenced by the distance from the emission point (217). When the propagation medium is homogeneous, but radiation particles penetrate the tissue, convolution with DPK can be a viable approach. When dealing with penetrating radiation in a non-homogeneous medium, direct MC methods should be considered, particularly for accurate absorbed dose estimates.
3.3.1 Local energy deposition
A simple method to calculate absorbed dose is the LED which only considers self-dose within the voxel and assumes that all produced energy is entirely deposited in the originating voxel. Absorbent fractions are set to one using this method (220). The range of the charged particles released during radioactive decay, i.e., the alpha, beta particles, and Auger electrons, which are the primary contributors to energy deposition, determines the validity of this assumption (221). The assumption is valid for the majority of radionuclides used in RPT, where the projected path length of charged particles in tissue is less than voxel dimensions (128). It is also applicable for 177Lutetium, given that the charged particles' range is typically within the dimensions used for clinical SPECT images (153). Even for radionuclides such as 90Yttrium, where the range of emitted charged particles is longer, most of the emitted charged particles remain within a projected path length of 5 mm, which aligns with SPECT voxel dimensions.
The validity of the LED assumption is contingent upon the gamma-ray yield, making it less applicable to gamma rays. Because the 177Lutetium 208.4 keV gamma ray has a relatively low yield (10.4%), there is a low possibility of cross-dose contribution from target organ gamma rays. This makes the method particularly accurate for assessing 177Lutetium toxicity studies (222–224). Radionuclides such as 131Iodine with higher gamma-ray yield a higher probability of cross-dose from the target organ. As a result, different approaches such as convolution with DPK (225, 226) and direct MC transport calculations (128, 227–230) are considered.
3.3.2 Convolution with point voxel kernels
Convolution with DPK is considered when photons have a path length that is longer than the spatial resolution of the reconstructed SPECT image. This approach employs a DPK that describes the radial absorbed dose in a uniform water medium when an isotropic point source is positioned at the center. The total energy released per unit mass (TERMA) in conjunction with the kernel is used to calculate the dose (19, 231). The TIA image is convolved with these functions to estimate the deposited energy distribution, thereby deriving the absorbed dose (232). Instead of calculating a continuous dose point kernel, a discrete DPK is often determined. The TIA image is then convolved with tissue-specific DPKs (233).
Originally, these methods used simulations in homogeneous phantoms (234, 235) that included comprehensive decay data for photons and mono-energetic charged particles. The techniques have since been expanded to density scaling for non-homogeneous tissue 177Lutetium (225, 226, 236). The methods proved less accurate when applied to a non-homogeneous medium. For inhomogeneities, DPKs have been implemented and compared with the direct MC gold standard (233). Due to the inadequate kernel size used in DPK dosimetry, the modified DPK overestimated the mean absorbed dose from the MC technique by 5% to 8%. Neural networks have been trained with DPKs to estimate absorbed dose calculations in kidneys (189, 237).
An AI-based deep learning approach for whole-body organ level dosimetry considering tissue inhomogeneity and patient-specific anatomy has also been described (238). In this novel method using patient-specific anatomy and the S-coefficient kernel, deep learning was used to predict the energy deposition and compared well with the direct MC dosimetry. Additionally, the techniques were marketed as a way around MC dosimetry's computational burden restriction. Even though deep learning-assisted dosimetry is a very promising advancement in absorbed dose estimation, the outcomes rely on how many training kernels must be assembled for a sizable patient cohort, all relevant tissue combinations, and a vast number of voxels to minimize inter-patient variability. The methods rely on simplified computations to estimate absorbed doses, as accounting for the complex effects of tissue heterogeneity in human anatomy on energy propagation and deposition proves challenging.
3.3.3 Direct Monte Carlo-based dosimetry
The complete meaning of patient-specific dosimetry entails a detailed evaluation of the absorbed dose in tumor and normal tissue provided by the MC-based pre-calculated S-coefficients that are tailored to the unique anatomy and heterogeneity of each patient. Direct MC-based absorbed dose calculations are accepted as the gold standard (239). However these methods are still computationally demanding for widespread clinical application and, as a result, have received limited acceptance as a standard for clinical dosimetry (128, 134). However, the concept has potential due to the rapidly increasing computational capacity.
In addition to computing patient-specific S-coefficients, MC absorbed dose calculations are capable of handling non-uniform absorbed dose estimates at a voxel level. This allows them to overcome the assumptions involved in the previously described approaches. Particle transport is simulated by the CT images, while radiation decay characteristics are simulated by the SPECT images. MC-based dosimetry uses predefined patient atomic composition and density information from CT images to model the emission from SPECT images and the transport of photons as well as charged particles across different structures in the body. To sample decay locations, it is assumed from the quantified SPECT images that each voxel represents activity from the relevant volume element in the patient. This is one of the reasons that accurate activity quantification is important. The density image volumes simulate the passage of charged particles and photons. In an image matrix with the same voxel dimensions as the reconstructed SPECT image, the energy deposition at each interaction site is scored to create the absorbed dose rate images (128). The patient's unique geometry can be used to construct source-target combinations rather than relying on reference phantom models. The S-coefficients are computed for each radiation source to the target combination and consider the radiation range for each source-to-target geometry (217). In addition to accounting for tissue heterogeneity and secondary particle emissions, MC dosimetry also takes tumor geometries and tissue type transitions into consideration (226, 240).
Small-scale dosimetry, which involves calculating radiation absorbed dose at sub-organ and sub-tumor levels, may be addressed using DPK and direct MC radiation transport simulations. In addition to accounting for tissue non-homogeneities, another advantage of MC simulation is its applicability to conditions where charged-particle equilibrium is not achieved, such as tissue interfaces (241). To use DPK in heterogeneous media, simple scaling factors are applied to those in water-equivalent media, producing results that closely approximate those of MC but with reduced computational times. Voxel dosimetry involves calculating radiation absorbed dose to tissue regions ranging from a few centimeters to hundreds of micrometers. This method is commonly associated with tomographic imaging such as PET/CT and SPECT/CT or autoradiographic techniques for activity quantification. However, in the context of multicellular, cellular, and subcellular dosimetry, there is a need for quantifying activity at smaller scales, typically ranging from tens to hundreds of micrometers (242). Acquiring activity distribution data directly from clinical tomographic gamma camera images at these small scales remains challenging due to the gamma camera's limited spatial resolution. Nonetheless, autoradiographic techniques may offer a solution by enabling the quantification of activity within groups of cells (multicellular) and even within single cells. In preclinical settings, small-scale or voxel dosimetry has become a more widely utilized approach due to the improved resolution capabilities of preclinical imaging (112). However, it's worth noting that conventional preclinical dose estimations often assume uniform distribution of activity and dose deposition within organs. This assumption may not reflect reality, particularly for β-/α-emitting radiopharmaceuticals, where tissue activity distribution can be heterogeneous.
Tissue non-homogeneities may be represented using BEDs and dose-volume histograms (DVHs) unique to each patient. DVHs demonstrate the relationship between a volume (%) that has received a specific absorbed dose as a function of the absorbed dose (27). In EBRT, DVHs are frequently employed to depict the dose distributions of tumors and OAR. It is unclear how to appropriately present the idea of DVHs as a reporting mechanism in the NM community (243). Applications of DVH and BEDs are still being developed for NM due to the limited spatial resolution of SPECT/CT images that do not support fine sub-region uptake (244). The use of DVH and BED in the NM clinics has been restricted due to the limited extent of their clinical validation.
Ideally, a dosimetry program should encompass all the necessary steps required in the clinical dosimetry workflow, commencing with dose calibrator and gamma camera calibrations and concluding with computation of dose estimates. This facilitates a modular approach, and the workflow becomes more user-friendly with incentives to complete all steps in a clinical context required to obtain accurate dosimetry. Several commercial software programs, listed in Table 3, are available for RPT dosimetry. These tools are based on different scientific methodologies comparable to the ones discussed above and have become valuable tools for clinical applications involving 177Lutetium dosimetry. It's worth noting that not all commercial dosimetry program encompasses every step of the dosimetry workflow. To enhance the precision of dose estimations, it is advisable to employ optimization techniques at each stage of the dosimetry process to minimize errors (207).
4 Practical prospects
RPT is an evolving field that requires continuous refinement. Its clinical significance involves comparing the administered activity tailored to individual patients with generic or fixed-activity regimens, demonstrating its impact on clinical outcomes. This strategy will help collect the crucial data needed for establishing the dose-response relationship in RPT. The focus of these initiatives ought to be on well-planned, multi-center clinical trials that compare the one-size-fits-all approach with dosimetry-based activity regimens. Following the defined tailored treatment in EBRT, patient-specific dosimetry-driven activity administrations ought to be given priority as standard procedures. Technological developments in imaging hardware resolve some of the historical issues that have been barriers to this endeavor. PET gamma cameras with full-ring detector geometries have made whole-body imaging more rapid and practical, while SPECT gamma cameras with solid-state detectors offer better energy resolution. Harmonizing the crucial first steps in the dosimetry chain holds the key to generating consistent dose estimates using patient-specific dosimetry. The calibration of dose calibrators and precise determination of gamma camera calibration parameters should be the first steps in the pursuit of this objective. These are easy but necessary first steps in obtaining reliable and consistent dosimetry readings.
5 Conclusion
There is growing evidence demonstrating inter-patient variability, suggesting that minimizing OAR toxicity and optimizing tumor management must be balanced. Despite these progressions, a fixed administered activity regime continues to be predominant, primarily because of uncertainties associated with dosimetry calculations for RPT. This state stems from several obstacles including the lack of standardized methods in the clinical dosimetry workflow, the arduous effort for successive imaging examinations, and the lack of comprehensive documentation correlating the administered activity to patient outcomes.
While patient-specific dosimetry plays a significant role in ensuring safety, its application lacks empirical support, a notion reinforced by the “tolerated holy gray” of maximum threshold values. Consequently, the more established empirical activity administration is often preferred over patient-specific dosimetry-driven activity administrations, which are perceived as complicated and time-consuming (253). Patient-specific dosimetry provides the evidence-based data necessary for personalized RPT, making it essential to achieving RPT's full potential as a precision-based cancer therapeutic alternative. The pursuit for personalized administered activities in RPT must remain a priority, guided by Einstein's dictum to be “…as simple as possible but no simpler.”
Author contributions
KR: Conceptualization, Methodology, Writing – original draft, Writing – review & editing, Formal Analysis, Visualization. MQ: Conceptualization, Writing – review & editing. HN: Writing – review & editing. CS: Writing – review & editing. LS: Writing – review & editing. SM: Writing – review & editing. AM: Writing – review & editing. MS: Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
KR extends gratitude to Johannes Van Staden and Hanlie du Raan, from the University of the Free State, Bloemfontein, South Africa, for their valuable contributions and insights.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Sgouros G, Hobbs RF. Dosimetry for radiopharmaceutical therapy. Semin Nucl Med. (2014) 44(3):172–8. doi: 10.1053/j.semnuclmed.2014.03.007
2. Gupta SK, Rangarajan V. Radiopharmaceutical biodistribution and dosimetry. In: Badria FA, editor. Radiopharmaceuticals—current Research for Better Diagnosis and Therapy. Eastern Virginia: IntechOpen (2022). p. 6–8.
3. Safety and health at work EU-OSHA. Directive 2013/59/Euratom—protection against ionising radiation. Updated (2019). Available online at: https://osha.europa.eu/en/legislation/directives/directive-2013-59-euratom-protection-against-ionising-radiation
4. Dietrich A, Koi L, Zöphel K, Sihver W, Kotzerke J, Baumann M, et al. Improving external beam radiotherapy by combination with internal irradiation. Br J Radiol. (2015) 88(1051):20150042. doi: 10.1259/bjr.20150042
5. van der Sar ECA, Braat AJAT, van der Voort- van Zyp JRN, van der Veen BS, van Leeuwen PJ, de Vries-Huizing DMV, et al. Tolerability of concurrent external beam radiotherapy and [177Lu]Lu-PSMA-617 for node-positive prostate cancer in treatment naïve patients, phase I study (PROQURE-I trial). BMC Cancer. (2023) 23:268. doi: 10.1186/s12885-023-10725-5
6. Yang Y, Xiong L, Li M, Jiang P, Wang J, Li C. Advances in radiotherapy and immunity in hepatocellular carcinoma. J Transl Med. (2023) 21(1):526. doi: 10.1186/s12967-023-04386-y
7. Zeng ZC, Tang ZY, Yang BH, Liu KD, Wu ZQ, Fan J, et al. Comparison between radioimmunotherapy and external beam radiation therapy for patients with hepatocellular carcinoma. Eur J Nucl Med Mol Imaging. (2002) 29(12):1657–68. doi: 10.1007/s00259-002-0996-x
8. NCI. Radiopharmaceuticals Emerging as New Cancer Therapy—NCI. (2020) (cited Oct 30, 2023). Available online at: https://www.cancer.gov/news-events/cancer-currents-blog/2020/radiopharmaceuticals-cancer-radiation-therapy#:∼:text=Radiopharmaceuticals%20consist%20of%20a%20radioactive,linker%20that%20joins%20the%20two.&text=The%20past%20two%20decades%20have,grow%2C%20divide%2C%20and%20spread
9. Stokke C, Kvassheim M, Blakkisrud J. Radionuclides for targeted therapy: physical properties. Molecules. (2022) 27(17):5429. doi: 10.3390/molecules27175429
10. Morris ZS, Wang AZ, Knox SJ. The radiobiology of radiopharmaceuticals. Semin Radiat Oncol. (2021) 31(1):20–7. doi: 10.1016/j.semradonc.2020.07.002
11. Plichta K, Graves S, Buatti J. Prostate-specific membrane antigen (PSMA) theranostics for treatment of oligometastatic prostate cancer. Int J Mol Sci. (2021) 22(22):2–10. doi: 10.3390/ijms222212095
12. Brans B, Bodei L, Giammarile F, Linden O, Luster M, Oyen WJG, et al. Clinical radionuclide therapy dosimetry: the quest for the ‘Holy Gray’. Eur J Nucl Med Mol Imaging. (2007) 34(5):772–86. doi: 10.1007/s00259-006-0338-5
13. Dale R, Carabe-Fernandez A. The radiobiology of conventional radiotherapy and its application to radionuclide therapy. Cancer Biother Radiopharm. (2005) 20(1):47–51. doi: 10.1089/cbr.2005.20.47
14. Sgouros G, Bodei L, McDevitt MR, Nedrow JR. Radiopharmaceutical therapy in cancer: clinical advances and challenges. Nat Rev Drug Discov. (2020) 19(9):589–608. doi: 10.1038/s41573-020-0073-9
15. Sgouros G, Chiu S, Pentlow KS, Brewster LJ, Kalaigian H, Baldwin B, et al. Three-dimensional dosimetry for radioimmunotherapy treatment planning. J Nucl Med. (1993) 34(9):1595–601.8394886
16. ICRU. ICRU Report 85a (Revised) Fundamental Quantities and Units for Ionizing Radiation—ICRU. Available online at: https://www.icru.org/report/fundamental-quantities-and-units-for-ionizing-radiation-icru-report-85a-revised/ (Accessed February 24, 2024).
17. Sapienza MT, Willegaignon J. Radionuclide therapy: current status and prospects for internal dosimetry in individualized therapeutic planning. Clinics (Sao Paulo). (2019) 74:e835. doi: 10.6061/clinics/2019/e835
18. Bodei L, Herrmann K, Schöder H, Scott AM, Lewis JS. Radiotheranostics in oncology: current challenges and emerging opportunities. Nat Rev Clin Oncol. (2022) 19(8):534–50. doi: 10.1038/s41571-022-00652-y
19. St. James S, Bednarz B, Benedict S, Buchsbaum JC, Dewaraja Y, Frey E, et al. Current status of radiopharmaceutical therapy. Int J Radiat Oncol. (2021) 109(4):891–901. doi: 10.1016/j.ijrobp.2020.08.035
20. Strigari L, Konijnenberg M, Chiesa C, Bardies M, Du Y, Gleisner KS, et al. The evidence base for the use of internal dosimetry in the clinical practice of molecular radiotherapy. Eur J Nucl Med Mol Imaging. (2014) 41(10):1976–88. doi: 10.1007/s00259-014-2824-5
21. Kozempel J, Mokhodoeva O, Vlk M. Progress in targeted alpha-particle therapy. What we learned about recoils release from in vivo generators. Molecules. (2018) 23(3):581. doi: 10.3390/molecules23030581
22. De Kruijff RM, Wolterbeek HT, Denkova AG. A critical review of alpha radionuclide therapy—how to deal with recoiling daughters? Pharmaceuticals. (2015) 8(2):321–36. doi: 10.3390/ph8020321
23. Trujillo-Nolasco M, Morales-Avila E, Cruz-Nova P, Katti KV, Ocampo-García B. Nanoradiopharmaceuticals based on alpha emitters: recent developments for medical applications. Pharmaceutics. (2021) 13(8):1123. doi: 10.3390/pharmaceutics13081123
24. Volkert WA, Goeckeler WF, Ehrhardt GJ, Ketring AR. Therapeutic radionuclides: production and decay property considerations. J Nucl Med. (1991) 32(1):174–85.1988628
25. Ersahin D, Doddamane I, Cheng D. Targeted radionuclide therapy. Cancers (Basel). (2011) 3(4):3838–55. doi: 10.3390/cancers3043838
26. National Research Council (US) and Institute of Medicine (US) Committee on State of the Science of Nuclear Medicine. Advancing Nuclear Medicine Through Innovation. Washington (DC): National Academies Press (US) (2007). (The National Academies Collection: Reports funded by National Institutes of Health).
27. Hindorf C. Internal Dosimetry. Chapter 18. Nuclear Medicine Physics: A Handbook for Teachers and Students Endorsed by: American Association of Physicists in Medicine (AAPM), Asia–Oceania Federation of Organizations for Medical Physics (AFOMP), Australasian College of Physical Scientists and Engineers in Medicine (ACPSEM), European Federation of Organisations for Medical Physics (EFOMP), Federation of African Medical Physics Organisations (FAMPO), World Federation of Nuclear Medicine and Biology (WFNMB). (2014); Available online at: http://inis.iaea.org/Search/search.aspx?orig_q=RN:46034096
28. Sjögreen Gleisner K, Spezi E, Solny P, Gabina PM, Cicone F, Stokke C, et al. Variations in the practice of molecular radiotherapy and implementation of dosimetry: results from a European survey. EJNMMI Phys. (2017) 4(1):28. doi: 10.1186/s40658-017-0193-4
29. Fendler WP, Rahbar K, Herrmann K, Kratochwil C, Eiber M. 177Lu-PSMA radioligand therapy for prostate cancer. J Nucl Med. (2017) 58(8):1196–200. doi: 10.2967/jnumed.117.191023
30. Stokke C, Gabiña PM, Solný P, Cicone F, Sandström M, Gleisner KS, et al. Dosimetry-based treatment planning for molecular radiotherapy: a summary of the 2017 report from the internal dosimetry task force. EJNMMI Phys. (2017) 4(1):27. doi: 10.1186/s40658-017-0194-3
31. Nelson BJB, Andersson JD, Wuest F. Targeted alpha therapy: progress in radionuclide production, radiochemistry, and applications. Pharmaceutics. (2020) 13(1):49. doi: 10.3390/pharmaceutics13010049
32. Van de Voorde M, Duchemin C, Heinke R, Lambert L, Chevallay E, Schneider T, et al. Production of Sm-153 with very high specific activity for targeted radionuclide therapy. Front Med (Lausanne). (2021) 8:675221. doi: 10.3389/fmed.2021.675221
33. Parlak Y, Gumuser G, Sayit E. Samarium-153 therapy and radiation dose for prostate cancer. In: Mohan R, editor. Prostate Cancer—leading-Edge Diagnostic Procedures and Treatments. Eastern Virginia: IntechOpen (2016). p. 81–91. https://www.intechopen.com/chapters/51764
34. Albertsson P, Bäck T, Bergmark K, Hallqvist A, Johansson M, Aneheim E, et al. Astatine-211 based radionuclide therapy: current clinical trial landscape. Front Med (Lausanne). (2023) 9:1076210. doi: 10.3389/fmed.2022.1076210
35. Eckerman K, Endo A. ICRP publication 107. Nuclear decay data for dosimetric calculations. Ann ICRP. (2008) 38(3):7–96. doi: 10.1016/j.icrp.2008.10.001
36. Simpkin DJ, Mackie TR. EGS4 Monte Carlo determination of the beta dose kernel in water. Med Phys. (1990) 17(2):179–86. doi: 10.1118/1.596565
37. Papadimitroulas P, Loudos G, Nikiforidis GC, Kagadis GC. A dose point kernel database using GATE Monte Carlo simulation toolkit for nuclear medicine applications: comparison with other Monte Carlo codes. Med Phys. (2012) 39(8):5238–47. doi: 10.1118/1.4737096
38. Duchemin C, Essayan M, Guertin A, Haddad F, Michel N, Métivier V. How to produce high specific activity tin-117 m using alpha particle beam. Appl Radiat Isot. (2016) 115:113–24. doi: 10.1016/j.apradiso.2016.06.016
39. Lepareur N, Lacœuille F, Bouvry C, Hindré F, Garcion E, Chérel M, et al. Rhenium-188 labeled radiopharmaceuticals: current clinical applications in oncology and promising perspectives. Front Med (Lausanne). (2019) 6:1–19. doi: 10.3389/fmed.2019.00132
40. Bodei L, Cremonesi M, Grana CM, Fazio N, Iodice S, Baio SM, et al. Peptide receptor radionuclide therapy with 177Lu-DOTATATE: the IEO phase I-II study. Eur J Nucl Med Mol Imaging. (2011) 38(12):2125–35. doi: 10.1007/s00259-011-1902-1
41. Kwekkeboom DJ, Bakker WH, Kooij PP, Konijnenberg MW, Srinivasan A, Erion JL, et al. [177Lu-DOTAOTyr3]octreotate: comparison with [111In-DTPAo]octreotide in patients. Eur J Nucl Med. (2001) 28(9):1319–25. doi: 10.1007/s002590100574
42. Pauwels S, Barone R, Walrand S, Borson-Chazot F, Valkema R, Kvols LK, et al. Practical dosimetry of peptide receptor radionuclide therapy with (90)Y-labeled somatostatin analogs. J Nucl Med. (2005) 46(Suppl 1):92S–8S.15653657
43. Naeem Z, Zahra UB, Numair Younis M, Khan IU, Shahid A. Lutetium-177 prostate specific membrane antigen therapy in a patient with double malignancy and single functioning kidney: a case report. Cureus. (2023) 15(3):e36938. doi: 10.7759/cureus.36938
44. Sutherland DEK, Kashyap R, Jackson P, Buteau JP, Murphy DG, Kelly B, et al. Safety of Lutetium-177 prostate-specific membrane antigen-617 (PSMA-617) radioligand therapy in the setting of severe renal impairment: a case report and literature review. Ther Adv Med Oncol. (2023) 15:17588359231177018. doi: 10.1177/1758835923117701
45. Okwundu N, Weil CR, Soares HP, Fine GC, Cannon DM. Case report: efficacy of lutetium-177 oxodotreotide for neuroendocrine tumor with central nervous system metastases. Frontiers in Nuclear Medicine. (2023) 3:1–7. doi: 10.3389/fnume.2023.1074948
46. Gudkov SV, Shilyagina NY, Vodeneev VA, Zvyagin AV. Targeted radionuclide therapy of human tumors. Int J Mol Sci. (2015) 17(1):33. doi: 10.3390/ijms17010033
47. LNHB. (cited Nov 6, 2021). Available online at: http://www.nucleide.org/
48. Uribe CF, Esquinas PL, Gonzalez M, Celler A. Characteristics of bremsstrahlung emissions of (177)Lu, (188)Re, and (90)Y for SPECT/CT quantification in radionuclide therapy. Phys Med. (2016) 32(5):691–700. doi: 10.1016/j.ejmp.2016.04.014
49. Beauregard JM, Hofman MS, Pereira JM, Eu P, Hicks RJ. Quantitative (177)Lu SPECT (QSPECT) imaging using a commercially available SPECT/CT system. Cancer Imaging. (2011) 11:56–66. doi: 10.1102/1470-7330.2011.0012
50. Sandström M, Garske-Román U, Johansson S, Granberg D, Sundin A, Freedman N. Kidney dosimetry during 177Lu-DOTATATE therapy in patients with neuroendocrine tumors: aspects on calculation and tolerance. Acta Oncol (Madr). (2018) 57(4):516–21. doi: 10.1080/0284186X.2017.1378431
51. Zaknun JJ, Bodei L, Mueller-Brand J, Pavel ME, Baum RP, Hörsch D, et al. The joint IAEA, EANM, and SNMMI practical guidance on peptide receptor radionuclide therapy (PRRNT) in neuroendocrine tumours. Eur J Nucl Med Mol Imaging. (2013) 40(5):800–16. doi: 10.1007/s00259-012-2330-6
52. Dash A, Pillai MRA, Knapp FF. Production of (177)Lu for targeted radionuclide therapy: available options. Nucl Med Mol Imaging. (2015) 49(2):85–107. doi: 10.1007/s13139-014-0315-z
53. Reijonen V, Kanninen LK, Hippeläinen E, Lou YR, Salli E, Sofiev A, et al. Multicellular dosimetric chain for molecular radiotherapy exemplified with dose simulations on 3D cell spheroids. Phys Med. (2017) 40:72–8. doi: 10.1016/j.ejmp.2017.07.012
54. Das T, Banerjee S. Theranostic applications of lutetium-177 in radionuclide therapy. Curr Radiopharm. (2016) 9(1):94–101. doi: 10.2174/1874471008666150313114644
55. Becx MN, Minczeles NS, Brabander T, de Herder WW, Nonnekens J, Hofland J. A clinical guide to peptide receptor radionuclide therapy with 177Lu-DOTATATE in neuroendocrine tumor patients. Cancers (Basel). (2022) 14(23):5792. doi: 10.3390/cancers14235792
56. Garkavij M, Nickel M, Sjögreen-Gleisner K, Ljungberg M, Ohlsson T, Wingårdh K, et al. 177Lu-[DOTA0,Tyr3] Octreotate therapy in patients with disseminated neuroendocrine tumors: analysis of dosimetry with impact on future therapeutic strategy. Cancer. (2010) 116(4 Suppl):1084–92. doi: 10.1002/cncr.24796
57. Arveschoug AK, Bekker AC, Iversen P, Bluhme H, Villadsen GE, Staanum PF. Extravasation of [177Lu]lu-DOTATOC: case report and discussion. EJNMMI Res. (2020) 10(1):68. doi: 10.1186/s13550-020-00658-6
58. Kam BLR, Teunissen JJM, Krenning EP, de Herder WW, Khan S, van Vliet EI, et al. Lutetium-labelled peptides for therapy of neuroendocrine tumours. Eur J Nucl Med Mol Imaging. (2012) 39(Suppl 1):103–12. doi: 10.1007/s00259-011-2039-y
59. Oyen WJG, Bodei L, Giammarile F, Maecke HR, Tennvall J, Luster M, et al. Targeted therapy in nuclear medicine—current status and future prospects. Ann Oncol. (2007) 18(11):1782–92. doi: 10.1093/annonc/mdm111
60. Kwekkeboom DJ, Teunissen JJ, Bakker WH, Kooij PP, de Herder WW, Feelders RA, et al. Radiolabeled somatostatin analog [177Lu-DOTA0,Tyr3]octreotate in patients with endocrine gastroenteropancreatic tumors. J Clin Oncol. (2005 Apr 20) 23(12):2754–62. doi: 10.1200/JCO.2005.08.066
61. Dewaraja Y, Wilderman S, Niedbala J, Frey K, Wong KK. Multi SPECT/CT-based patient specific lesion and kidney dosimetry for verification of simpler approaches for treatment planning in Lu-177 DOTATATE PRRT. J Nucl Med. (2019) 60(supplement 1):1626–1626.
62. Sandström M, Freedman N, Fröss-Baron K, Kahn T, Sundin A. Kidney dosimetry in 777 patients during 177Lu-DOTATATE therapy: aspects on extrapolations and measurement time points. EJNMMI Phys. (2020) 7(1):73. doi: 10.1186/s40658-020-00339-2
63. Sundlöv A, Sjögreen-Gleisner K, Svensson J, Ljungberg M, Olsson T, Bernhardt P, et al. Individualised 177Lu-DOTATATE treatment of neuroendocrine tumours based on kidney dosimetry. Eur J Nucl Med Mol Imaging. (2017) 44(9):1480–9. doi: 10.1007/s00259-017-3678-4
64. Willowson KP, Eslick E, Ryu H, Poon A, Bernard EJ, Bailey DL. Feasibility and accuracy of single time point imaging for renal dosimetry following 177Lu-DOTATATE (‘lutate’) therapy. EJNMMI Phys. (2018) 5(1):33. doi: 10.1186/s40658-018-0232-9
65. Kurth J, Heuschkel M, Tonn A, Schildt A, Hakenberg OW, Krause BJ, et al. Streamlined schemes for dosimetry of 177Lu-labeled PSMA targeting radioligands in therapy of prostate cancer. Cancers (Basel). (2021) 13(15):3884. doi: 10.3390/cancers13153884
66. Resch S, Takayama Fouladgar S, Zacherl M, Sheikh GT, Liubchenko G, Rumiantcev M, et al. Investigation of image-based lesion and kidney dosimetry protocols for 177Lu-PSMA-I&T therapy with and without a late SPECT/CT acquisition. EJNMMI Phys. (2023) 10(1):11. doi: 10.1186/s40658-023-00529-8
67. Delker A, Schleske M, Liubchenko G, Berg I, Zacherl MJ, Brendel M, et al. Biodistribution and dosimetry for combined [177Lu]Lu-PSMA-I&T/[225Ac]Ac-PSMA-I&T therapy using multi-isotope quantitative SPECT imaging. Eur J Nucl Med Mol Imaging. (2023) 50(5):1280–90. doi: 10.1007/s00259-022-06092-1
68. Emami B, Lyman J, Brown A, Coia L, Goitein M, Munzenrider JE, et al. Tolerance of normal tissue to therapeutic irradiation. Int J Radiat Oncol Biol Phys. (1991) 21(1):109–22. doi: 10.1016/0360-3016(91)90171-Y
69. Cassady JR. Clinical radiation nephropathy. Int J Radiat Oncol Biol Phys. (1995) 31(5):1249–56. doi: 10.1016/0360-3016(94)00428-N
70. Fowler JF. Radiobiological aspects of low dose rates in radioimmunotherapy. Int J Radiat Oncol Biol Phys. (1990) 18(5):1261–9. doi: 10.1016/0360-3016(90)90467-X
71. International C. ICRU Report 67: absorbed-dose specification in nuclear medicine. International C, editor. J ICRU. (2002) 2(1):5–105.
72. Langmuir VK, Fowler JF, Knox SJ, Wessels BW, Sutherland RM, Wong JY. Radiobiology of radiolabeled antibody therapy as applied to tumor dosimetry. Med Phys. (1993) 20(2 Pt 2):601–10. doi: 10.1118/1.597055
73. Valkema R, Pauwels SA, Kvols LK, Kwekkeboom DJ, Jamar F, de Jong M, et al. Long-term follow-up of renal function after peptide receptor radiation therapy with (90)Y-DOTA(0),tyr(3)-octreotide and (177)Lu-DOTA(0), tyr(3)-octreotate. J Nucl Med. (2005) 46(Suppl 1):83S–91S.15653656
74. Bodei L, Cremonesi M, Ferrari M, Pacifici M, Grana CM, Bartolomei M, et al. Long-term evaluation of renal toxicity after peptide receptor radionuclide therapy with 90Y-DOTATOC and 177Lu-DOTATATE: the role of associated risk factors. Eur J Nucl Med Mol Imaging. (2008) 35(10):1847–56. doi: 10.1007/s00259-008-0778-1
75. Strosberg J, El-Haddad G, Wolin E, Hendifar A, Yao J, Chasen B, et al. Phase 3 trial of 177Lu-dotatate for midgut neuroendocrine tumors. N Engl J Med. (2017) 376(2):125–35. doi: 10.1056/NEJMoa1607427
76. Sartor O, de Bono J, Chi KN, Fizazi K, Herrmann K, Rahbar K, et al. Lutetium-177-PSMA-617 for metastatic castration-resistant prostate cancer. N Engl J Med. (2021) 385(12):1091–103. doi: 10.1056/NEJMoa2107322
77. Mahajan S, Grewal RK, Friedman KP, Schöder H, Pandit-Taskar N. Assessment of salivary gland function after 177Lu-PSMA radioligand therapy: current concepts in imaging and management. Transl Oncol. (2022) 21:101445. doi: 10.1016/j.tranon.2022.101445
78. Langbein T, Kulkarni HR, Schuchardt C, Mueller D, Volk GF, Baum RP. Salivary gland toxicity of PSMA-targeted radioligand therapy with 177Lu-PSMA and combined 225Ac- and 177Lu-labeled PSMA ligands (TANDEM-PRLT) in advanced prostate cancer: a single-center systematic investigation. Diagnostics (Basel). (2022) 12(8):1926. doi: 10.3390/diagnostics12081926
79. Radojewski P, Dumont R, Marincek N, Brunner P, Mäcke HR, Müller-Brand J, et al. Towards tailored radiopeptide therapy. Eur J Nucl Med Mol Imaging. (2015) 42(8):1231–7. doi: 10.1007/s00259-015-3030-9
80. Carollo A, Papi S, Chinol M. Lutetium-177 labeled peptides: the European institute of oncology experience. Curr Radiopharm. (2016) 9(1):19–32. doi: 10.2174/1874471008666150313111633
81. Liu Z, Ma T, Liu H, Jin Z, Sun X, Zhao H, et al. 177Lu-labeled antibodies for EGFR-targeted SPECT/CT imaging and radioimmunotherapy in a preclinical head and neck carcinoma model. Mol Pharmaceutics. (2014) 11(3):800–7. doi: 10.1021/mp4005047
82. Shinto AS, Shibu D, Kamaleshwaran KK, Das T, Chakraborty S, Banerjee S, et al. 177Lu-EDTMP for treatment of bone pain in patients with disseminated skeletal metastases. J Nucl Med Technol. (2014) 42(1):55–61. doi: 10.2967/jnmt.113.132266
83. Almqvist Y, Steffen AC, Tolmachev V, Divgi CR, Sundin A. In vitro and in vivo characterization of 177Lu-huA33: a radioimmunoconjugate against colorectal cancer. Nucl Med Biol. (2006) 33(8):991–8. doi: 10.1016/j.nucmedbio.2006.09.003
84. Kim K, Kim SJ. Lu-177-based peptide receptor radionuclide therapy for advanced neuroendocrine tumors. Nucl Med Mol Imaging. (2018) 52(3):208–15. doi: 10.1007/s13139-017-0505-6
85. Niaz MO, Sun M, Ramirez-Fort MK, Niaz MJ. Review of lutetium-177-labeled anti-prostate-specific membrane antigen monoclonal antibody J591 for the treatment of metastatic castration-resistant prostate cancer. Cureus. (2020) 12(2):e7107. doi: 10.7759/cureus.7107
86. Cunningham AD, Qvit N, Mochly-Rosen D. Peptides and peptidomimetics as regulators of protein-protein interactions. Curr Opin Struct Biol. (2017) 44:59–66. doi: 10.1016/j.sbi.2016.12.009
87. Grob NM, Schibli R, Béhé M, Mindt TL. Improved tumor-targeting with peptidomimetic analogs of minigastrin 177Lu-PP-F11N. Cancers (Basel). (2021) 13(11):2629. doi: 10.3390/cancers13112629
88. Ladrière T, Faudemer J, Levigoureux E, Peyronnet D, Desmonts C, Vigne J. Safety and therapeutic optimization of lutetium-177 based radiopharmaceuticals. Pharmaceutics. (2023) 15(4):1240. doi: 10.3390/pharmaceutics15041240
89. Repetto-Llamazares AHV, Larsen RH, Mollatt C, Lassmann M, Dahle J. Biodistribution and dosimetry of 177Lu-tetulomab, a new radioimmunoconjugate for treatment of non-hodgkin lymphoma. Curr Radiopharm. (2013) 6(1):20. doi: 10.2174/1874471011306010004
90. Chen J, Qi L, Tang Y, Tang G, Gan Y, Cai Y. Current role of prostate-specific membrane antigen-based imaging and radioligand therapy in castration-resistant prostate cancer. Front Cell Dev Biol. (2022) 10:958180. doi: 10.3389/fcell.2022.958180
91. Sjögreen Gleisner K, Chouin N, Gabina PM, Cicone F, Gnesin S, Stokke C, et al. EANM dosimetry committee recommendations for dosimetry of 177Lu-labelled somatostatin-receptor- and PSMA-targeting ligands. Eur J Nucl Med Mol Imaging. (2022) 49(6):1778–809. doi: 10.1007/s00259-022-05727-7
92. Zakaly HMH, Mostafa MYA, Deryabina D, Zhukovsky M. Comparative studies on the potential use of 177Lu-based radiopharmaceuticals for the palliative therapy of bone metastases. Int J Radiat Biol. (2020) 96(6):779–89. doi: 10.1080/09553002.2020.1729441
93. Lee H. Relative efficacy of 225Ac-PSMA-617 and 177Lu-PSMA-617 in prostate cancer based on subcellular dosimetry. Mol Imaging Radionucl Ther. (2022) 31(1):1–6. doi: 10.4274/mirt.galenos.2021.63308
94. Makvandi M, Dupis E, Engle JW, Nortier FM, Fassbender ME, Simon S, et al. Alpha-emitters and targeted alpha therapy in oncology: from basic science to clinical investigations. Targ Oncol. (2018) 13(2):189–203. doi: 10.1007/s11523-018-0550-9
95. Kratochwil C, Bruchertseifer F, Giesel FL, Weis M, Verburg FA, Mottaghy F, et al. 225Ac-PSMA-617 for PSMA-targeted α-radiation therapy of metastatic castration-resistant prostate cancer. J Nucl Med. (2016) 57(12):1941–4. doi: 10.2967/jnumed.116.178673
96. Frantellizzi V, Cosma L, Brunotti G, Pani A, Spanu A, Nuvoli S, et al. Targeted alpha therapy with thorium-227. Cancer Biother Radiopharm. (2020) 35(6):437–45. doi: 10.1089/cbr.2019.3105
97. Seo Y. Quantitative imaging of alpha-emitting therapeutic radiopharmaceuticals. Nucl Med Mol Imaging. (2019) 53(3):182–8. doi: 10.1007/s13139-019-00589-8
98. Hindorf C, Chittenden S, Aksnes AK, Parker C, Flux GD. Quantitative imaging of 223Ra-chloride (alpharadin) for targeted alpha-emitting radionuclide therapy of bone metastases. Nucl Med Commun. (2012) 33(7):726. doi: 10.1097/MNM.0b013e328353bb6e
99. Kratochwil C, Bruchertseifer F, Rathke H, Bronzel M, Apostolidis C, Weichert W, et al. Targeted α-therapy of metastatic castration-resistant prostate cancer with 225Ac-PSMA-617: dosimetry estimate and empiric dose finding. J Nucl Med. (2017) 58(10):1624–31. doi: 10.2967/jnumed.117.191395
100. Belli ML, Sarnelli A, Mezzenga E, Cesarini F, Caroli P, Di Iorio V, et al. Targeted alpha therapy in mCRPC (metastatic castration-resistant prostate cancer) patients: predictive dosimetry and toxicity modeling of 225Ac-PSMA (prostate-specific membrane antigen). Front Oncol. (2020) 10:531660. doi: 10.3389/fonc.2020.531660
101. Gosewisch A, Schleske M, Gildehaus FJ, Berg I, Kaiser L, Brosch J, et al. Image-based dosimetry for 225Ac-PSMA-I&T therapy using quantitative SPECT. Eur J Nucl Med Mol Imaging. (2021) 48(4):1260–1. doi: 10.1007/s00259-020-05024-1
102. Sgouros G, Bolch WE, Chiti A, Dewaraja YK, Emfietzoglou D, Hobbs RF, et al. ICRU REPORT 96, dosimetry-guided radiopharmaceutical therapy. J ICRU. (2021) 21(1):1–212. doi: 10.1177/14736691211060117
103. Kratochwil C, Schmidt K, Afshar-Oromieh A, Bruchertseifer F, Rathke H, Morgenstern A, et al. Targeted alpha therapy of mCRPC: dosimetry estimate of 213Bismuth-PSMA-617. Eur J Nucl Med Mol Imaging. (2018) 45(1):31–7. doi: 10.1007/s00259-017-3817-y
104. Sinnes JP, Bauder-Wüst U, Schäfer M, Moon ES, Kopka K, Rösch F. 68Ga, 44sc and 177Lu-labeled AAZTA5-PSMA-617: synthesis, radiolabeling, stability and cell binding compared to DOTA-PSMA-617 analogues. EJNMMI Radiopharm Chem. (2020) 5(1):28. doi: 10.1186/s41181-020-00107-8
105. Umbricht CA, Benešová M, Schmid RM, Türler A, Schibli R, van der Meulen NP, et al. 44Sc-PSMA-617 For radiotheragnostics in tandem with 177Lu-PSMA-617-preclinical investigations in comparison with 68Ga-PSMA-11 and 68Ga-PSMA-617. EJNMMI Res. (2017) 7(1):9. doi: 10.1186/s13550-017-0257-4
106. Radiations NRC (US) Committee on the Biological Effects of Ionising Radiations. Dosimetry of alpha particles. In: Health Risks of Radon and Other Internally Deposited Alpha-Emitters: Beir IV. Washingtom: National Academies Press (US) (1988). p. 397–414. https://www.ncbi.nlm.nih.gov/books/NBK218131/
107. Peter R, Sandmaier BM, Dion MP, Frost SHL, Santos EB, Kenoyer A, et al. Small-scale (sub-organ and cellular level) alpha-particle dosimetry methods using an iQID digital autoradiography imaging system. Sci Rep. (26) 12(1):17934. doi: 10.1038/s41598-022-22664-5
108. Sgouros G, Roeske JC, McDevitt MR, Palm S, Allen BJ, Fisher DR, et al. MIRD pamphlet no. 22 (abridged): radiobiology and dosimetry of α-particle emitters for targeted radionuclide therapy. J Nucl Med. (2010) 51(2):311–28. doi: 10.2967/jnumed.108.058651
109. Li WB, Hofmann W, Friedland W. Microdosimetry and nanodosimetry for internal emitters. Radiat Meas. (2018) 115:29–42. doi: 10.1016/j.radmeas.2018.05.013
110. Rumiantcev M, Li WB, Lindner S, Liubchenko G, Resch S, Bartenstein P, et al. Estimation of relative biological effectiveness of 225Ac compared to 177Lu during [225Ac]ac-PSMA and [177Lu]lu-PSMA radiopharmaceutical therapy using TOPAS/TOPAS-nBio/MEDRAS. EJNMMI Phys. (2023) 10(1):53. doi: 10.1186/s40658-023-00567-2
111. Bolcaen J, Combrink N, Spoormans K, More S, Vandevoorde C, Fisher R, et al. Biodosimetry, can it find its way to the nuclear medicine clinic? Front Nucl Med. (2023) 3:1–17. Available online at: https://www.frontiersin.org/articles/10.3389/fnume.2023.1209823. doi: 10.3389/fnume.2023.1209823
112. Vargas CS, Struelens L, D’Huyvetter M, Caveliers V, Covens P. A realistic multiregion mouse kidney dosimetry model to support the preclinical evaluation of potential nephrotoxicity of radiopharmaceutical therapy. J Nucl Med. (2023) 64(3):493–9. doi: 10.2967/jnumed.122.264453
113. O’Donoghue J, Zanzonico P, Humm J, Kesner A. Dosimetry in radiopharmaceutical therapy. J Nucl Med. (2022) 63(10):1467–74. doi: 10.2967/jnumed.121.262305
114. Salerno KE, Roy S, Ribaudo C, Fisher T, Patel RB, Mena E, et al. A primer on radiopharmaceutical therapy. Int J Radiat Oncol Biol Phys. (2023) 115(1):48–59. doi: 10.1016/j.ijrobp.2022.08.010
115. Valkema R, Konijnenberg M, Pauwels S, Kvols L, Krenning E. How much dosimetry is needed for therapeutic applications? (2010). Available online at: https://inis.iaea.org/collection/NCLCollectionStore/_Public/42/022/42022126.pdf (Accessed December 14, 2023).
116. Shcherbinin S, Celler A, Belhocine T, Vanderwerf R, Driedger A. Accuracy of quantitative reconstructions in SPECT/CT imaging. Phys Med Biol. (2008) 53(17):4595–604. doi: 10.1088/0031-9155/53/17/009
117. Van Wyk BP, Hasford F, Nyakale NE, Vangu MMDT, Oelofse B, Leboea HM. Critical appraisal of radionuclide calibrators and gamma cameras prior to lutetium-177 internal dosimetry at two South African hospitals. World J Nucl Med. (2022) 21(1):44–51. doi: 10.1055/s-0042-1746173
118. D’Arienzo M, Cazzato M, Cozzella ML, Cox M, D’Andrea M, Fazio A, et al. Gamma camera calibration and validation for quantitative SPECT imaging with (177)Lu. Appl Radiat Isot. (2016) 112:156–64. doi: 10.1016/j.apradiso.2016.03.007
119. Dewaraja YK, Frey EC, Sgouros G, Brill AB, Roberson P, Zanzonico PB, et al. MIRD pamphlet No. 23: quantitative SPECT for patient-specific 3-dimensional dosimetry in internal radionuclide therapy. J Nucl Med. (2012) 53(8):1310–25. doi: 10.2967/jnumed.111.100123
120. Ramonaheng K, van Staden JA, du Raan H. The effect of calibration factors and recovery coefficients on 177Lu SPECT activity quantification accuracy: a Monte Carlo study. EJNMMI Phys. (2021) 8(1):27. doi: 10.1186/s40658-021-00365-8
121. Halty A, Badel JN, Kochebina O, Sarrut D. Image-based SPECT calibration based on the evaluation of the fraction of activity in the field of view. EJNMMI Phys. (2018) 5:11. doi: 10.1186/s40658-018-0209-8
122. Peters SMB, Van Der Werf NR, Segbers M, Van Velden FHP, Wierts R, Blokland KAK, et al. Towards standardization of absolute SPECT/CT quantification: a multi-center and multi-vendor phantom study. EJNMMI Phys. (2019) 6(1):29. doi: 10.1186/s40658-019-0268-5
123. Bailey DL, Hofman MS, Forwood NJ, O’Keefe GJ, Scott AM, Wyngaardt Wv, et al. Accuracy of dose calibrators for 68Ga PET imaging: unexpected findings in a multicenter clinical pretrial assessment. J Nucl Med. (2018) 59(4):636–8. doi: 10.2967/jnumed.117.202861
124. Bailey DL, Hennessy TM, Willowson KP, Henry EC, Chan DLH, Aslani A, et al. In vivo quantification of 177Lu with planar whole-body and SPECT/CT gamma camera imaging. EJNMMI Phys. (2015) 2(1):20. doi: 10.1186/s40658-015-0123-2
125. Bailey DL, Willowson KP. Quantitative SPECT/CT: SPECT joins PET as a quantitative imaging modality. Eur J Nucl Med Mol Imaging. (2014) 41(Suppl 1):S17–25. doi: 10.1007/s00259-013-2542-4
126. Hippeläinen E, Tenhunen M, Mäenpää H, Sohlberg A. Quantitative accuracy of (177)Lu SPECT reconstruction using different compensation methods: phantom and patient studies. EJNMMI Res. (2016) 6(1):16. doi: 10.1186/s13550-016-0172-0
127. Ljungberg M, Celler A, Konijnenberg MW, Eckerman KF, Dewaraja YK, Sjögreen-Gleisner K, et al. MIRD pamphlet no. 26: joint EANM/MIRD guidelines for quantitative 177Lu SPECT applied for dosimetry of radiopharmaceutical therapy. J Nucl Med. (2016) 57(1):151–62. doi: 10.2967/jnumed.115.159012
128. Ljungberg M, Sjögreen Gleisner K. Personalized dosimetry for radionuclide therapy using molecular imaging tools. Biomedicines. (2016) 4(4):25. doi: 10.3390/biomedicines4040025
129. Peters SMB, Meyer Viol SL, van der Werf NR, de Jong N, van Velden FHP, Meeuwis A, et al. Variability in lutetium-177 SPECT quantification between different state-of-the-art SPECT/CT systems. EJNMMI Phys. (2020) 7(1):9. doi: 10.1186/s40658-020-0278-3
130. International Atomic Energy Agency. Quantitative SPECT OR SPECT/CT imaging. In: Buvat I, Frey EC, Green AJ, Ljungberg M, editors. IAEA Quantitative Nuclear Medicine Imaging: Concepts, Requirements and Methods, Human Health Report No.9. Vienna, Austria: IAEA (2014). p. 27–36.
131. Wevrett J, Fenwick A, Scuffham J, Johansson L, Gear J, Schlögl S, et al. Inter-comparison of quantitative imaging of lutetium-177 (177Lu) in European hospitals. EJNMMI Phys. (2018) 5(1):17. doi: 10.1186/s40658-018-0213-z
132. Stabin M, Brill A, Farncombe T, King M. Image quantification for radiation dosimetry analyses. J Nucl Med. (2006) 47(suppl 1):491P–491P.
133. Frey EC, Humm JL, Ljungberg M. Accuracy and precision of radioactivity quantification in nuclear medicine images. Semin Nucl Med. (2012) 42(3):208–18. doi: 10.1053/j.semnuclmed.2011.11.003
134. Ljungberg M, Sjögreen-Gleisner K. The accuracy of absorbed dose estimates in tumours determined by quantitative SPECT: a Monte Carlo study. Acta Oncol. (2011) 50(6):981–9. doi: 10.3109/0284186X.2011.584559
135. Marin G, Vanderlinden B, Karfis I, Guiot T, Wimana Z, Flamen P, et al. Accuracy and precision assessment for activity quantification in individualized dosimetry of 177Lu-DOTATATE therapy. EJNMMI Phys. (2017) 4:1–15. doi: 10.1186/s40658-017-0174-7
136. Sandström M, Garske U, Granberg D, Sundin A, Lundqvist H. Individualized dosimetry in patients undergoing therapy with (177)Lu-DOTA-D-phe (1)-tyr (3)-octreotate. Eur J Nucl Med Mol Imaging. (2010) 37(2):212–25. doi: 10.1007/s00259-009-1216-8
137. Uribe CF, Esquinas PL, Tanguay J, Gonzalez M, Gaudin E, Beauregard JM, et al. Accuracy of 177Lu activity quantification in SPECT imaging: a phantom study. EJNMMI Phys. (2017) 4:5238–47. doi: 10.1186/s40658-016-0170-3
138. Raskin S, Gamliel D, Abookasis D, Ben-Haim S, Chicheportiche A. Towards accurate 177Lu SPECT activity quantification and standardization using lesion-to-background voxel ratio. EJNMMI Phys. (2023) 10(1):5. doi: 10.1186/s40658-023-00526-x
139. de Nijs R, Lagerburg V, Klausen TL, Holm S. Improving quantitative dosimetry in (177)Lu-DOTATATE SPECT by energy window-based scatter corrections. Nucl Med Commun. (2014) 35(5):522–33. doi: 10.1097/MNM.0000000000000079
140. He B, Nikolopoulou A, Osborne J, Vallabhajosula S, Goldsmith S. Quantitative SPECT imaging with Lu-177: a physical phantom evaluation. J Nucl Med. (2012) 53(supplement 1):2407–2407.
141. Shcherbinin S, Piwowarska-Bilska H, Celler A, Birkenfeld B. Quantitative SPECT/CT reconstruction for 177Lu and 177Lu/90Y targeted radionuclide therapies. Phys Med Biol. (2012) 57(18):5733. doi: 10.1088/0031-9155/57/18/5733
142. Huizing DMV, Sinaasappel M, Dekker MC, Stokkel MPM, van der Veen BJ. 177Lutetium SPECT/CT: evaluation of collimator, photopeak and scatter correction. J Appl Clin Med Phys. (2020) 21(9):272. doi: 10.1002/acm2.12991
143. Asmi H, Bentayeb F, Bouzekraoui Y, Bonutti F, Douama S. Energy window and collimator optimization in Lutetium-177 single-photon emission computed tomography imaging using monte carlo simulation. Indian J Nucl Med. (2020) 35(1):36–9. doi: 10.4103/ijnm.IJNM_121_19
144. Karimi Ghodoosi E, D’Alessandria C, Li Y, Bartel A, Köhner M, Höllriegl V, et al. The effect of attenuation map, scatter energy window width, and volume of interest on the calibration factor calculation in quantitative 177Lu SPECT imaging: simulation and phantom study. Phys Med. (2018) 56:74–80. doi: 10.1016/j.ejmp.2018.11.009
145. Delker A. Investigation of Models with Temporal and Spatial Interference in Image Based Dosimetry of 177Lu- Labelled Radioligand Therapies. (82).
146. Miller C, Filipow L, Jackson S. A review of activity quantification by planar imaging methods. J Nucl Med Technol. (1995) 23(1):7.
147. Minarik D, Sjögreen K, Ljungberg M. A new method to obtain transmission images for planar whole-body activity quantification. Cancer Biother Radiopharm. (2005) 20(1):72–6. doi: 10.1089/cbr.2005.20.72
148. Pereira JM, Stabin MG, Lima FRA, Guimarães MICC, Forrester JW. Image quantification for radiation dose calculations–limitations and uncertainties. Health Phys. (2010) 99(5):688–701. doi: 10.1097/HP.0b013e3181e28cdb
149. Sjögreen K, Ljungberg M, Wingårdh K, Minarik D, Strand SE. The LundADose method for planar image activity quantification and absorbed-dose assessment in radionuclide therapy. Cancer Biother Radiopharm. (2005) 20(1):92–7.
150. Cachovan M, Vija AH, Hornegger J, Kuwert T. Quantification of 99mTc-DPD concentration in the lumbar spine with SPECT/CT. EJNMMI Res. (2013) 3:45. doi: 10.1186/2191-219X-3-45
151. Gnesin S, Leite Ferreira P, Malterre J, Laub P, Prior JO, Verdun FR. Phantom validation of Tc-99m absolute quantification in a SPECT/CT commercial device. Comput Math Methods Med. (2016) 2016:1–6. doi: 10.1155/2016/4360371
152. Willowson K, Bailey DL, Baldock C. Quantitative SPECT reconstruction using CT-derived corrections. Phys Med Biol. (2008) 53(12):3099–112. doi: 10.1088/0031-9155/53/12/002
153. Ljungberg M, Sjögreen Gleisner K. Hybrid imaging for patient-specific dosimetry in radionuclide therapy. Diagnostics. (2015) 5(3):296–317. doi: 10.3390/diagnostics5030296
154. Beykan S, Tran-Gia J, Borup Jensen S, Lassmann M. Is a single late SPECT/CT based kidney 177Lu-dosimetry superior to hybrid dosimetry with sequential multiple time-point whole-body planar scans in combination with an early SPECT/CT? Phys Med. (2022) 100:39–50. doi: 10.1016/j.ejmp.2022.06.002
155. Roth D, Gustafsson J, Sundlöv A, Sjögreen Gleisner K. A method for tumor dosimetry based on hybrid planar-SPECT/CT images and semiautomatic segmentation. Med Phys. (2018) 45(11):5004–18. doi: 10.1002/mp.13178
156. Denisova N. Bayesian maximum-A-Posteriori approach with global and local regularization to image reconstruction problem in medical emission tomography. Entropy (Basel). (2019) 21(11):1108. doi: 10.3390/e21111108
157. Shepp LA, Vardi Y. Maximum likelihood reconstruction for emission tomography. IEEE Transactions on Medical Imaging. (1982) 1(2):113–22. doi: 10.1109/TMI.1982.4307558
158. Hudson HM, Larkin RS. Accelerated image reconstruction using ordered subsets of projection data. IEEE Trans Med Imaging. (1994) 13(4):601–9. doi: 10.1109/42.363108
159. Cherry SR, Sorenson JA, Phelps ME. Chapter 16—tomographic reconstruction in nuclear medicine. In: Physics in Nuclear Medicine (Fourth Edition). Philadelphia: W.B. Saunders (2012). p. 253–77.
160. Kemin H. Impact of different reconstruction algorithms and OSEM reconstruction parameters on quantitative results in SPECT/CT. J Nucl Med. (2018) 59(supplement 1):1800–1800.
161. Dewaraja YK, Wilderman SJ, Ljungberg M, Koral KF, Zasadny K, Kaminiski MS. Accurate dosimetry in 131I radionuclide therapy using patient-specific, 3-dimensional methods for SPECT reconstruction and absorbed dose calculation. J Nucl Med. (2005) 46(5):840–9.15872359
162. Gedik GK, Sari O. Influence of single photon emission computed tomography (SPECT) reconstruction algorithm on diagnostic accuracy of parathyroid scintigraphy: comparison of iterative reconstruction with filtered backprojection. Indian J Med Res. (2017) 145(4):479–87. doi: 10.4103/ijmr.IJMR_305_15
163. Cheng Z, Wen J, Huang G, Yan J. Applications of artificial intelligence in nuclear medicine image generation. Quant Imaging Med Surg. (2021) 11(6):2792–822. doi: 10.21037/qims-20-1078
164. Ravishankar S, Ye JC, Fessler JA. Image reconstruction: from sparsity to data-adaptive methods and machine learning. Proc IEEE Inst Electr Electron Eng. (2020) 108(1):86–109. doi: 10.1109/JPROC.2019.2936204
165. Patton JA, Turkington TG. SPECT/CT physical principles and attenuation correction. J Nucl Med Technol. (2008) 36(1):1–10. doi: 10.2967/jnmt.107.046839
166. Takahashi Y, Murase K, Higashino H, Mochizuki T, Motomura N. Attenuation correction of myocardial SPECT images with x-ray CT: effects of registration errors between x-ray CT and SPECT. Ann Nucl Med. (2002) 16(6):431–5. doi: 10.1007/BF02990083
167. Ritt P, Vija H, Hornegger J, Kuwert T. Absolute quantification in SPECT. Eur J Nucl Med Mol Imaging. (2011) 38(Suppl 1):S69–77. doi: 10.1007/s00259-011-1770-8
168. Prieto Canalejo MA, Palau San Pedro A, Geronazzo R, Minsky DM, Juárez-Orozco LE, Namías M. Synthetic attenuation correction maps for SPECT imaging using deep learning: a study on myocardial perfusion imaging. Diagnostics. (2023) 13(13):2214. doi: 10.3390/diagnostics13132214
169. Shanbhag AD, Miller RJH, Pieszko K, Lemley M, Kavanagh P, Feher A, et al. Deep learning-based attenuation correction improves diagnostic accuracy of cardiac SPECT. J Nucl Med. (2023) 64(3):472–8. doi: 10.2967/jnumed.122.264429
170. Shi L, Onofrey JA, Liu H, Liu YH, Liu C. Deep learning-based attenuation map generation for myocardial perfusion SPECT. Eur J Nucl Med Mol Imaging. (2020) 47(10):2383–95. doi: 10.1007/s00259-020-04746-6
171. Sakaguchi K, Kaida H, Yoshida S, Ishii K. Attenuation correction using deep learning for brain perfusion SPECT images. Ann Nucl Med. (2021) 35(5):589–99. doi: 10.1007/s12149-021-01600-z
172. Grimes J. Patient-specific Internal Dose Calculation Techniques for Clinical use in Targeted Radionuclide Therapy. Vancouver: University of British Columbia (2013).
173. Jaszczak RJ, Greer KL, Floyd CE, Harris CC, Coleman RE. Improved SPECT quantification using compensation for scattered photons. J Nucl Med. (1984) 25(8):893–900.6611390
174. Ogawa K, Harata Y, Ichihara T, Kubo A, Hashimoto S. A practical method for position-dependent compton-scatter correction in single photon emission CT. IEEE Trans Med Imaging. (1991) 10(3):408–12. doi: 10.1109/42.97591
175. Hutton BF, Buvat I, Beekman FJ. Review and current status of SPECT scatter correction. Phys Med Biol. (2011) 56(14):R85–112. doi: 10.1088/0031-9155/56/14/R01
176. Frey EC, Tsui BMW. New method for modeling the spatially-variant, object-dependent scatter response function in SPECT. IEEE Nuclear Science Symposium & Medical Imaging Conference, IEEE (1996).
177. Frey EC, Ju ZW, Tsui BMW. A fast projector-backprojector pair modeling the asymmetric, spatially varying scatter response function for scatter compensation in SPECT imaging. IEEE Trans Nucl Sci. (1993) 40(4):1192–7. doi: 10.1109/23.256735
178. Visvikis D, Lambin P, Beuschau Mauridsen K, Hustinx R, Lassmann M, Rischpler C, et al. Application of artificial intelligence in nuclear medicine and molecular imaging: a review of current status and future perspectives for clinical translation. Eur J Nucl Med Mol Imaging. (2022) 49(13):4452–63. doi: 10.1007/s00259-022-05891-w
179. Xiang H, Lim H, Fessler JA, Dewaraja YK. A deep neural network for fast and accurate scatter estimation in quantitative SPECT/CT under challenging scatter conditions. Eur J Nucl Med Mol Imaging. (2020) 47(13):2956–67. doi: 10.1007/s00259-020-04840-9
180. Yang J, Park D, Gullberg GT, Seo Y. Joint correction of attenuation and scatter in image space using deep convolutional neural networks for dedicated brain 18F-FDG PET. Phys Med Biol. (2019) 64(7):075019. doi: 10.1088/1361-6560/ab0606
181. McMillan AB, Bradshaw TJ. Artificial intelligence-based data corrections for attenuation and scatter in position emission tomography and single-photon emission computed tomography. PET Clin. (2021) 16(4):543–52. doi: 10.1016/j.cpet.2021.06.010
182. Zaidi H. Collimator-Detector response compensation in SPECT. In: Zaidi H, editor. Quantitative Analysis in Nuclear Medicine Imaging. New York: Springer US (2006). p. 141–2.
183. Dewaraja YK, Koral KF, Fessler JA. Regularized reconstruction in quantitative SPECT using CT side information from hybrid imaging. Phys Med Biol. (2010) 55(9):2523–39. doi: 10.1088/0031-9155/55/9/007
184. International Atomic Energy Agency. Partial volume correction. In: Buvat I, Frey EC, Green AJ, Ljungberg M, editors. IAEA Quantitative Nuclear Medicine Imaging: Concepts, Requirements and Methods, Human Health Report No.9. Vienna, Austria: IAEA (2014). p. 34.
185. Erlandsson K, Thomas B, Dickson J, Hutton BF. Partial volume correction in SPECT reconstruction with OSEM. Nucl Instrum Methods Phys Res Sect Accel Spectrometers Detect Assoc Equip. (2011) 648:S85–8. doi: 10.1016/j.nima.2010.12.106
186. Tran-Gia J, Salas-Ramirez M, Lassmann M. What you see is not what you get: on the accuracy of voxel-based dosimetry in molecular radiotherapy. J Nucl Med. (2020) 61(8):1178–86. doi: 10.2967/jnumed.119.231480
187. Bettinardi V, Castiglioni I, De Bernardi E, Gilardi MC. PET quantification: strategies for partial volume correction. Clin Transl Imaging. (2014) 2(3):199–218. doi: 10.1007/s40336-014-0066-y
188. Hoffman EJ, Huang SC, Phelps ME. Quantitation in positron emission computed tomography: 1. Effect of object size. J Comput Assist Tomogr. (1979) 3(3):299. doi: 10.1097/00004728-197906000-00001
189. Marquis H, Willowson K, Bailey D. Partial volume effect in SPECT & PET imaging and impact on radionuclide dosimetry estimates. Asia Ocean J Nucl Med Biol. (2023) 11(1):44–54. doi: 10.22038/AOJNMB.2022.63827.1448
190. Tran-Gia J, Lassmann M. Optimizing image quantification for Lu-177 SPECT/CT based on a 3D printed 2-compartment kidney phantom. J Nucl Med. (2017) 117:200170.
191. Grings A, Jobic C, Kuwert T, Ritt P. The magnitude of the partial volume effect in SPECT imaging of the kidneys: a phantom study. EJNMMI Phys. (2022) 9:1–16. doi: 10.1186/s40658-022-00446-2
192. Loening AM, Gambhir SS. AMIDE: a free software tool for multimodality medical image analysis. Mol Imaging. (2023) 1(2):131–7. Available online at: http://journals.sagepub.com/doi/abs/10.1162/15353500200303133
193. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. (2012) 9(7):676–82. doi: 10.1038/nmeth.2019
194. Sjögreen K, Ljungberg M, Strand SE. An activity quantification method based on registration of CT and whole-body scintillation camera images, with application to 131I. J Nucl Med. (2002) 43(7):972–82.
195. Adams MC, Turkington TG, Wilson JM, Wong TZ. A systematic review of the factors affecting accuracy of SUV measurements. AJR Am J Roentgenol. (2010) 195(2):310–20. doi: 10.2214/AJR.10.4923
196. van Sluis J, Noordzij W, de Vries EGE, Kok IC, de Groot DJA, Jalving M, et al. Manual versus artificial intelligence-based segmentations as a pre-processing step in whole-body PET dosimetry calculations. Mol Imaging Biol. (2023) 25(2):435–41. doi: 10.1007/s11307-022-01775-5
197. Loevinger R, Budinger TF, Watson EE. Mird Primer for Absorbed Dose Calculations. Revised, Subsequent Edition. New York, NY: Society of Nuclear Medicine (1991). p. 128.
198. Ljungberg M, Sjogreen Gleisner K. Three-dimensional image-based dosimetry in radionuclide therapy. IEEE Trans Radiation Plasma Med Sci. (2018) 2(6):1–1. doi: 10.1109/TRPMS.2018.2860563
199. Staanum PF, Frellsen AF, Olesen ML, Iversen P, Arveschoug AK. Practical kidney dosimetry in peptide receptor radionuclide therapy using [177Lu]Lu-DOTATOC and [177Lu]Lu-DOTATATE with focus on uncertainty estimates. EJNMMI Phys. (2021) 8:78. doi: 10.1186/s40658-021-00422-2
200. Siegel JA, Thomas SR, Stubbs JB, Stabin MG, Hays MT, Koral KF, et al. MIRD pamphlet no. 16: techniques for quantitative radiopharmaceutical biodistribution data acquisition and analysis for use in human radiation dose estimates. J Nucl Med. (1999) 40(2):37S–61S.10025848
201. Metz CE, Atkins FB, Beck RN. The geometric transfer function component for scintillation camera collimators with straight parallel holes. Phys Med Biol. (1980) 25(6):1059–70. doi: 10.1088/0031-9155/25/6/003
202. Jackson PA, Hofman MS, Hicks RJ, Scalzo M, Violet J. Radiation dosimetry in 177Lu-PSMA-617 therapy using a single posttreatment SPECT/CT scan: a novel methodology to generate time- and tissue-specific dose factors. J Nucl Med. (2020) 61(7):1030–6. doi: 10.2967/jnumed.119.233411
203. Glatting G, Kletting P, Reske SN, Hohl K, Ring C. Choosing the optimal fit function: comparison of the akaike information criterion and the F-test. Med Phys. (2007) 34(11):4285–92. doi: 10.1118/1.2794176
204. Kletting P, Schimmel S, Hänscheid H, Luster M, Fernández M, Nosske D, et al. The NUKDOS software for treatment planning in molecular radiotherapy. Z Med Phys. (2015) 25(3):264–74. doi: 10.1016/j.zemedi.2015.01.001
205. Kletting P, Kull T, Reske SN, Glatting G. Comparing time activity curves using the Akaike information criterion. Phys Med Biol. (2009) 54(21):N501–507. doi: 10.1088/0031-9155/54/21/N01
206. Gleisner KS, Brolin G, Sundlöv A, Mjekiqi E, Östlund K, Tennvall J, et al. Long-term retention of 177Lu/177mlu-DOTATATE in patients investigated by γ-spectrometry and γ-camera imaging. J Nucl Med. (2015) 56(7):976–84. doi: 10.2967/jnumed.115.155390
207. Gear J, Cox M, Gustafsson J, Sjogreen Gleisner K, Murray I, Glatting G, et al. EANM practical guidance on uncertainty analysis for molecular radiotherapy absorbed dose calculations. Eur J Nucl Med Mol Imaging. (2018) 45:2456–74. doi: 10.1007/s00259-018-4136-7
208. Peterson AB, Mirando DM, Dewaraja YK. Accuracy and uncertainty analysis of reduced time point imaging effect on time-integrated activity for 177Lu-DOTATATE PRRT in patients and clinically realistic simulations. EJNMMI Res. (2023) 13(1):57. doi: 10.1186/s13550-023-01007-z
209. Stabin MG, Siegel JA. Physical models and dose factors for use in internal dose assessment. Health Phys. (2003) 85(3):294–310. doi: 10.1097/00004032-200309000-00006
210. Capala J, Graves SA, Scott A, Sgouros G, James SS, Zanzonico P, et al. Dosimetry for radiopharmaceutical therapy: current practices and commercial resources. J Nucl Med. (2021) 62(Suppl 3):3S–11S. doi: 10.2967/jnumed.121.262749
211. Cristy M, Eckerman KF. Specific Absorbed Fractions of Energy at Various Ages from Internal Photon Sources: 7, Adult Male. Oak Ridge National Lab. (1987). Report No.: ORNL/TM–8381/V7. Available online at: http://inis.iaea.org/Search/search.aspx?orig_q=RN:19012902
212. Snyder W, Ford M, Warner G. MIRD Pamphlet No 5. Revised. Estimated of Specific Absorbed Fractions for Photon Sources Uniformly Distributed in Various Organs of a Heterogeneous Phantom. (1978).
213. Stabin M, Brill A, Segars W, Emmons M, Gesner J, Milam R. Realistic phantom series for OLINDA/EXM 2.0. J Nucl Med. (2006) 47(suppl 1):156P–156P.
214. Stabin M, Emmons MA, Segars WP, Fernald M, Brill AB. ICRP-89 based adult and pediatric phantom series. J Nucl Med. (2008) 49(supplement 1):14P–14P.
215. Menzel HG, Clement C, DeLuca P. ICRP Publication 110. Realistic reference phantoms: an ICRP/ICRU joint effort. A report of adult reference computational phantoms. Ann ICRP. (2009) 39(2):1–164. doi: 10.1016/j.icrp.2009.07.001
216. Segars W, Tsui B, Lalush D, Frey E, King M, Manocha D. Development and application of the new dynamic nurbs-based cardiac-torso (NCAT) phantom. J Nucl Med. (2001) 42:23.
217. Bardiès M, Gear JI. Scientific developments in imaging and dosimetry for molecular radiotherapy. Clin Oncol (R Coll Radiol). (2021) 33(2):117–24. doi: 10.1016/j.clon.2020.11.005
218. Akhavanallaf A, Shiri I, Arabi H, Zaidi H. Whole-body voxel-based internal dosimetry using deep learning. Eur J Nucl Med Mol Imaging. (2021) 48(3):670–82. doi: 10.1007/s00259-020-05013-4
219. Kramer R, Cassola VF, Khoury HJ, Vieira JW, Lima VdM, Brown KR. FASH And MASH: female and male adult human phantoms based on polygon mesh surfaces: II. Dosimetric calculations. Phys Med Biol. (2010) 55(1):163–89. doi: 10.1088/0031-9155/55/1/010
220. van den Hoven AF, Rosenbaum CENM, Elias SG, de Jong HWAM, Koopman M, Verkooijen HM, et al. Insights into the dose-response relationship of radioembolization with resin 90Y-microspheres: a prospective cohort study in patients with colorectal cancer liver metastases. J Nucl Med. (2016) 57(7):1014–9. doi: 10.2967/jnumed.115.166942
221. Huizing DMV, de Wit-van der Veen BJ, Verheij M, Stokkel MPM. Dosimetry methods and clinical applications in peptide receptor radionuclide therapy for neuroendocrine tumours: a literature review. EJNMMI Res. (2018) 8:89. doi: 10.1186/s13550-018-0443-z
222. Gustafsson J, Brolin G, Cox M, Ljungberg M, Johansson L, Gleisner KS. Uncertainty propagation for SPECT/CT-based renal dosimetry in 177Lu peptide receptor radionuclide therapy. Phys Med Biol. (2015) 60(21):8329. doi: 10.1088/0031-9155/60/21/8329
223. Hippeläinen E, Tenhunen M, Sohlberg A. Fast voxel-level dosimetry for (177)Lu labelled peptide treatments. Phys Med Biol. (2015) 60(17):6685–700. doi: 10.1088/0031-9155/60/17/6685
224. Svensson J, Berg G, Wängberg B, Larsson M, Forssell-Aronsson E, Bernhardt P. Renal function affects absorbed dose to the kidneys and haematological toxicity during 177Lu-DOTATATE treatment. Eur J Nucl Med Mol Imaging. (2015) 42(6):947–55. doi: 10.1007/s00259-015-3001-1
225. Dieudonné A, Hobbs RF, Lebtahi R, Maurel F, Baechler S, Wahl RL, et al. Study of the impact of tissue density heterogeneities on 3-dimensional abdominal dosimetry: comparison between dose kernel convolution and direct Monte Carlo methods. J Nucl Med. (2013) 54(2):236–43. doi: 10.2967/jnumed.112.105825
226. Dieudonné A, Hobbs RF, Bolch WE, Sgouros G, Gardin I. Fine-resolution voxel S values for constructing absorbed dose distributions at variable voxel size. J Nucl Med. (2010) 51(10):1600–7. doi: 10.2967/jnumed.110.077149
227. Agostinelli S, Allison J, Amako K, Apostolakis J, Araujo H, Arce P, et al. Geant4—a simulation toolkit. Nuclear Instruments and Methods in Physics Research Section A. (2003) 506(3):250–303. doi: 10.1016/S0168-9002(03)01368-8
228. Kawrakow I, Mainegra-Hing E, Tessier F. The EGSnrc Code System: Monte Carlo Simulation of Electron and Photon Transport. Semantic Scholar. (2016). Available online at: https://www.semanticscholar.org/paper/The-EGSnrc-Code-System%3A-Monte-Carlo-Simulation-of-Kawrakow-Mainegra%E2%80%90Hing/dd0c3a1d3c33cf0679999e835465b450d284da27 (Accessed December 03, 2023).
229. Larsson E, Ljungberg M, Strand SE, Jönsson BA. Monte carlo calculations of absorbed doses in tumours using a modified MOBY mouse phantom for pre-clinical dosimetry studies. Acta Oncol. (2011) 50(6):973–80. doi: 10.3109/0284186X.2011.582517
230. Waters L, Mckinney G, Durkee J, Fensin M, Hendricks J, James M, et al. The MCNPX monte carlo radiation transport code. AIP Conf Proc. (2007) 896:81–90. doi: 10.1063/1.2720459
231. Lee MS, Hwang D, Kim JH, Lee JS. Deep-dose: a voxel dose estimation method using deep convolutional neural network for personalized internal dosimetry. Sci Rep. (2019) 9(1):10308. doi: 10.1038/s41598-019-46620-y
232. Giap HB, Macey DJ, Podoloff DA. Development of a SPECT-based three-dimensional treatment planning system for radioimmunotherapy. J Nucl Med. (1995) 36(10):1885–94.7562060
233. Götz T, Schmidkonz C, Lang EW, Maier A, Kuwert T, Ritt P. A comparison of methods for adapting 177Lu dose-voxel-kernels to tissue inhomogeneities. Phys Med Biol. (2019) 64(24):245011. doi: 10.1088/1361-6560/ab5b81
234. Berger MJ. Distribution of absorbed dose around point sources of electrons and beta particles in water and other media. J Nucl Med. (1971) 2(Suppl 5):5–23.
235. Pacilio M, Amato E, Lanconelli N, Basile C, Torres LA, Botta F, et al. Differences in 3D dose distributions due to calculation method of voxel S-values and the influence of image blurring in SPECT. Phys Med Biol. (2015) 60(5):1945–64. doi: 10.1088/0031-9155/60/5/1945
236. Sanchez-Garcia M, Gardin I, Lebtahi R, Dieudonné A. Implementation and validation of collapsed cone superposition for radiopharmaceutical dosimetry of photon emitters. Phys Med Biol. (2015) 60(20):7861–76. doi: 10.1088/0031-9155/60/20/7861
237. Götz T, Lang E, Schmidkonz C, Kuwert T, Ludwig B. Dose voxel kernel prediction with neural networks for radiation dose estimation. Zeitschrift für Medizinische Physik. (2020) 31:1–14. doi: 10.1016/j.zemedi.2020.09.005
238. Akhavanallaf A, Shiri I, Arabi H, Zaidi H. Deep learning-assisted whole-body voxel-based internal dosimetry. 2020 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC) (2020), Boston, MA, USA: IEEE. p. 1–3
239. Bardies M, Ljungberg M. Monte Carlo Codes in Radionuclide Therapy. IOP Publishing. Lund University. (2011). Available online at: https://portal.research.lu.se/en/publications/monte-carlocodes-in-radionuclide-therapy (Accessed December 8, 2023).
240. Lanconelli N, Pacilio M, Lo Meo S, Botta F, Di Dia A, Aroche AT, et al. A free database of radionuclide voxel S values for the dosimetry of nonuniform activity distributions. Phys Med Biol. (2012) 57(2):517–33. doi: 10.1088/0031-9155/57/2/517
241. Bolch WE, Bouchet LG, Robertson JS, Wessels BW, Siegel JA, Howell RW, et al. MIRD pamphlet No. 17: the dosimetry of nonuniform activity distributions—radionuclide S values at the voxel level. J Nucl Med. (1999) 40(1):11S–36S.9935083
242. Yorke ED, Williams LE, Demidecki AJ, Heidorn DB, Roberson PL, Wessels BW. Multicellular dosimetry for beta-emitting radionuclides: autoradiography, thermoluminescent dosimetry and three-dimensional dose calculations. Med Phys. (1993) 20(2 Pt 2):543–50. doi: 10.1118/1.597050
243. Xiao Y, Roncali E, Hobbs R, James SS, Bednarz B, Benedict S, et al. Toward individualized voxel-level dosimetry for radiopharmaceutical therapy. Int J Radiat Oncol Biol Phys. (2021) 109(4):902–4. doi: 10.1016/j.ijrobp.2020.08.026
244. Helisch A, Förster GJ, Reber H, Buchholz HG, Arnold R, Göke B, et al. Pre-therapeutic dosimetry and biodistribution of86Y-DOTA-Phe1-Tyr3-octreotide versus 111In-pentetreotide in patients with advanced neuroendocrine tumours. Eur J Nucl Med Mol Imaging. (2004) 31(10):1386–92. doi: 10.1007/s00259-004-1561-6
245. Stabin MG, Wendt RE, Flux GD. RADAR guide: standard methods for calculating radiation doses for radiopharmaceuticals, part 2—data analysis and dosimetry. J Nucl Med. (2022) 63(3):485–92. doi: 10.2967/jnumed.121.262034
246. RADAR Home. (cited Nov 21, 2023). Available online at: https://www.doseinfo-radar.com/
247. Stabin MG, Siegel JA. RADAR dose estimate report: a compendium of radiopharmaceutical dose estimates based on OLINDA/EXM version 2.0. J Nucl Med. (2018) 59(1):154–60. doi: 10.2967/jnumed.117.196261
248. Huizing DMV, Peters SMB, Versleijen MWJ, Martens E, Verheij M, Sinaasappel M, et al. A head-to-head comparison between two commercial software packages for hybrid dosimetry after peptide receptor radionuclide therapy. EJNMMI Phys. (2020) 7(1):36. doi: 10.1186/s40658-020-00308-9
249. Kupitz D, Wetz C, Wissel H, Wedel F, Apostolova I, Wallbaum T, et al. Software-assisted dosimetry in peptide receptor radionuclide therapy with 177Lutetium-DOTATATE for various imaging scenarios. PLoS One. (2017) 12(11):e0187570. doi: 10.1371/journal.pone.0187570
250. Maughan NM, Garcia-Ramirez J, Arpidone M, Swallen A, Laforest R, Goddu SM, et al. Validation of post-treatment PET-based dosimetry software for hepatic radioembolization of yttrium-90 microspheres. Med Phys. (2019) 46(5):2394–402. doi: 10.1002/mp.13444
251. Hippeläinen ET, Tenhunen MJ, Mäenpää HO, Heikkonen JJ, Sohlberg AO. Dosimetry software Hermes internal radiation dosimetry: from quantitative image reconstruction to voxel-level absorbed dose distribution. Nucl Med Commun. (2017) 38(5):357–65. doi: 10.1097/MNM.0000000000000662
252. Denis-Bacelar A, Chittenden S, Divoli A, Gear J, Flux G. Qdose, a 3D treatment planning system for molecular radiotherapy. EANM conference. Conference Record; Milan, Italy, 2012. p. S351 (2). https://www.quantitativedose.com/
Keywords: 177Lutetium, SPECT, PET, theranostics, radiopharmaceutical therapy, patient-specific dosimetry, activity quantification, absorbed dose
Citation: Ramonaheng K, Qebetu M, Ndlovu H, Swanepoel C, Smith L, Mdanda S, Mdlophane A and Sathekge M (2024) Activity quantification and dosimetry in radiopharmaceutical therapy with reference to 177Lutetium. Front. Nucl. Med. 4:1355912. doi: 10.3389/fnume.2024.1355912
Received: 14 December 2023; Accepted: 12 March 2024;
Published: 28 March 2024.
Edited by:
Hesham M. H. Zakaly, Ural Federal University, RussiaReviewed by:
Jim O' Doherty, Siemens Healthineers, United StatesWeibo Li, Helmholtz Association of German Research Centres (HZ), Germany
© 2024 Ramonaheng, Qebetu, Ndlovu, Swanepoel, Smith, Mdanda, Mdlophane and Sathekge. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Keamogetswe Ramonaheng, kamo.ramonaheng@sanumeri.co.za