Proof-of-concept of DosiTest: A virtual multicentric clinical trial for assessing uncertainties in molecular radiotherapy dosimetry.

Clinical dosimetry in molecular radiotherapy (MRT) is a multi-step procedure, prone to uncertainties at every stage of the dosimetric workflow. These are difficult to assess, especially as some are complex or even impossible to measure experimentally. The DosiTest project was initiated to assess the variability associated with clinical dosimetry, by setting up a 'virtual' multicentric clinical dosimetry trial based on Monte Carlo (MC) modelling. A reference patient model with a realistic geometry and activity input for a specific tracer is considered. Reference absorbed dose rate distribution maps are generated at various time-points from MC modelling, combining precise information on density and activity distributions (voxel wise). Then, centre-specific calibration and patient SPECT/CT datasets are modelled, on which the clinical centres can perform clinical (i.e. image-based) dosimetry. The results of this dosimetric analysis can be benchmarked against the reference dosimetry to assess the variability induced by implementing different clinical dosimetry approaches. The feasibility of DosiTest is presented here for a clinical situation of therapeutic administration of 177Lu-DOTATATE (Lutathera®) peptide receptor radionuclide therapy (PRRT). From a real patient dataset composed of 5 SPECT/CT images and associated calibrations, we generated the reference absorbed dose rate images with GATE. Then, simulated SPECT/CT image generation based on GATE was performed, both for a calibration phantom and virtual patient images. Based on this simulated dataset, image-based dosimetry could be performed, and compared with reference dosimetry. The good agreement, between real and simulated images, and between reference and image-based dosimetry established the proof of concept of DosiTest.

Clinical implementation of PLANET® Dose for dosimetric assessment after [177Lu]Lu-DOTA-TATE: comparison with Dosimetry Toolkit® and OLINDA/EXM® V1.0.

BACKGROUND: The aim of this study was to compare a commercial dosimetry workstation (PLANET® Dose) and the dosimetry approach (GE Dosimetry Toolkit® and OLINDA/EXM® V1.0) currently used in our department for quantification of the absorbed dose (AD) to organs at risk after peptide receptor radionuclide therapy with [177Lu]Lu-DOTA-TATE. METHODS: An evaluation on phantom was performed to determine the SPECT calibration factor variations over time and to compare the Time Integrated Activity Coefficients (TIACs) obtained with the two approaches. Then, dosimetry was carried out with the two tools in 21 patients with neuroendocrine tumours after the first and second injection of 7.2 ± 0.2 GBq of [177Lu]Lu-DOTA-TATE (40 dosimetry analyses with each software). SPECT/CT images were acquired at 4 h, 24 h, 72 h and 192 h post-injection and were reconstructed using the Xeleris software (General Electric). The liver, spleen and kidneys masses and TIACs were determined using Dosimetry Toolkit® (DTK) and PLANET® Dose. The ADs were calculated using OLINDA/EXM® V1.0 and the Local Deposition Method (LDM) or Dose voxel-Kernel convolution (DK) on PLANET® Dose. RESULTS: With the phantom, the 3D calibration factors showed a slight variation (0.8% and 3.3%) over time, and TIACs of 225.19 h and 217.52 h were obtained with DTK and PLANET® Dose, respectively. In patients, the root mean square deviation value was 8.9% for the organ masses, 8.1% for the TIACs, and 9.1% and 7.8% for the ADs calculated with LDM and DK, respectively. The Lin's concordance correlation coefficient was 0.99 and the Bland-Altman plot analysis estimated that the AD value difference between methods ranged from - 0.75 to 0.49 Gy, from - 0.20 to 0.64 Gy, and from - 0.43 to 1.03 Gy for 95% of the 40 liver, kidneys and spleen dosimetry analyses. The dosimetry method had a minor influence on AD differences compared with the image registration and organ segmentation steps. CONCLUSIONS: The ADs to organs at risk obtained with the new workstation PLANET® Dose are concordant with those calculated with the currently used software and in agreement with the literature. These results validate the use of PLANET® Dose in clinical routine for patient dosimetry after targeted radiotherapy with [177Lu]Lu-DOTA-TATE.

Scientific Developments in Imaging and Dosimetry for Molecular Radiotherapy.

Molecular radiotherapy is a rapidly developing field with new vector and isotope combinations continually added to market. As with any radiotherapy treatment, it is vital that the absorbed dose and toxicity profile are adequately characterised. Methodologies for absorbed dose calculations for radiopharmaceuticals were generally developed to characterise stochastic effects and not suited to calculations on a patient-specific basis. There has been substantial scientific and technological development within the field of molecular radiotherapy dosimetry to answer this challenge. The development of imaging systems and advanced processing techniques enable the acquisition of accurate measurements of radioactivity within the body. Activity assessment combined with dosimetric models and radiation transport algorithms make individualised absorbed dose calculations not only feasible, but commonplace in a variety of commercially available software packages. The development of dosimetric parameters beyond the absorbed dose has also allowed the possibility to characterise the effect of irradiation by including biological parameters that account for radiation absorbed dose rates, gradients and spatial and temporal energy distribution heterogeneities. Molecular radiotherapy is in an exciting time of its development and the application of dosimetry in this field can only have a positive influence on its continued progression.

An innovative in vitro device providing continuous low doses of γ-rays mimicking exposure to the space environment: A dosimetric study.

Astronauts are exposed to microgravity and chronic irradiation but experimental conditions combining these two factors are difficult to reproduce on earth. We have created an experimental device able to combine chronic irradiation and altered gravity that may be used for cell cultures or plant models in a ground based facility. Irradiation was provided by thorium nitrate powder, conditioned so as to constitute a sealed source that could be placed in an incubator. Cell plates or plant seedlings could be placed in direct contact with the source or at various distances above it. Moreover, a random positioning machine (RPM) could be positioned on the source to simulate microgravity. The activity of the source was established using the Bateman formula. The spectrum of the source, calculated according to the natural decrease of radioactivity and the gamma spectrometry, showed very good adequacy. The experimental fluence was close to the theoretical fluence evaluation, attesting its uniform distribution. A Monte Carlo model of the irradiation device was processed by GATE code. Dosimetry was performed with radiophotoluminescent dosimeters exposed for one month at different locations (x and y axes) in various cell culture conditions. Using the RPM placed on the source, we reached a mean absorbed dose of gamma rays of (0.33 ± 0.17) mSv per day. In conclusion, we have elaborated an innovative device allowing chronic radiation exposure to be combined with altered gravity. Given the limited access to the International Space Station, this device could be useful to researchers interested in the field of space biology.

Optimization of GATE simulations for whole-body planar scintigraphic acquisitions using the XCAT male phantom with 177Lu-DOTATATE biokinetics in a Siemens Symbia T2.

Simulations of planar whole body acquisitions in therapeutic procedures are often extensively time-consuming and therefore rarely used. However, optimising tools and variance reduction techniques can be employed to overcome this problem. In this paper, a variety of features available in GATE are explored and their capabilities to reduce simulation time are evaluated. For this purpose, the male XCAT phantom was used as a virtual patient with 177Lu-DOTATATE pharmacokinetic for whole body planar acquisition simulations in a Siemens Symbia T2 model. Activity distribution was divided into 8 compartments that were simulated separately. GATE optimization techniques included reducing the amount of time spent in both voxel and detector tracking. Some acceleration techniques led to a decrease of CPU-time by a factor of 167, while image statistics were kept constant. In that context, the simulation of therapeutic procedure imaging would still require 46days on a single CPU, but this could be reduced to hours on a dedicated cluster.

Monte Carlo dose calculation in presence of low-density media: Application to lung SBRT treated during DIBH.

PURPOSE: Commercial algorithms used in Radiotherapy include approximations that are generally acceptable. However their limits can be seen when confronted with small fields and low-density media. These conditions exist during the treatment of lung cancers with Stereotactic Body Radiation Therapy (SBRT) achieved with the "Deep Inspiration Breath Hold" (DIBH) technique. A Monte Carlo (MC) model of a linear accelerator was used to assess the performance of two algorithms (Varian Acuros and AAA) in these conditions. This model is validated using phantoms with different densities. Lastly, results for SBRT cases are compared to both Acuros and AAA. METHODS: A Varian TrueBeam linac was modeled using GATE/Geant4 and validated by comparing dose distributions for simple fields to measurements in water and in heterogeneous phantoms composed of PMMA and two types of cork (corresponding to lung densities during free-breathing and DIBH). Experimental measurements are also compared to AAA and Acuros. Finally, results of Acuros/AAA are compared to MC for a clinical case (SBRT during DIBH). RESULTS: Based on 1D gamma index comparisons with measurements in water, the TrueBeam model was validated (>97% of points passed this test). In heterogeneous phantoms, and in particular for small field sizes, very low density (0.12g.cm-3) and at the edge of the field, MC model was still in good agreement with measurements whilst AAA and Acuros showed discrepancies. With the patient CT, similar differences between MC and AAA/Acuros were observed for static fields but disappeared using an SBRT arc field. CONCLUSIONS: Our MC model is validated and limits of commercial algorithms are shown in very low densities.

Implementation of new physics models for low energy electrons in liquid water in Geant4-DNA.

A new alternative set of elastic and inelastic cross sections has been added to the very low energy extension of the Geant4 Monte Carlo simulation toolkit, Geant4-DNA, for the simulation of electron interactions in liquid water. These cross sections have been obtained from the CPA100 Monte Carlo track structure code, which has been a reference in the microdosimetry community for many years. They are compared to the default Geant4-DNA cross sections and show better agreement with published data. In order to verify the correct implementation of the CPA100 cross section models in Geant4-DNA, simulations of the number of interactions and ranges were performed using Geant4-DNA with this new set of models, and the results were compared with corresponding results from the original CPA100 code. Good agreement is observed between the implementations, with relative differences lower than 1% regardless of the incident electron energy. Useful quantities related to the deposited energy at the scale of the cell or the organ of interest for internal dosimetry, like dose point kernels, are also calculated using these new physics models. They are compared with results obtained using the well-known Penelope Monte Carlo code.

Realistic multi-cellular dosimetry for 177Lu-labelled antibodies: model and application.

Current preclinical dosimetric models often fail to take account of the complex nature of absorbed dose distribution typical of in vitro clonogenic experiments in targeted radionuclide therapy. For this reason, clonogenic survival is often expressed as a function of added activity rather than the absorbed dose delivered to cells/cell nuclei. We designed a multi-cellular dosimetry model that takes into account the realistic distributions of cells in the Petri dish, for the establishment of survival curves as a function of the absorbed dose. General-purpose software tools were used for the generation of realistic, randomised 3D cell culture geometries based on experimentally determined parameters (cell size, cell density, cluster density, average cluster size, cell cumulated activity). A mixture of Monte Carlo and analytical approaches was implemented in order to achieve as accurate as possible results while reducing calculation time. The model was here applied to clonogenic survival experiments carried out to compare the efficacy of Betalutin®, a novel 177Lu-labelled antibody radionuclide conjugate for the treatment of non-Hodgkin lymphoma, to that of 177Lu-labelled CD20-specific (rituximab) and non-specific antibodies (Erbitux) on lymphocyte B cells. The 3D cellular model developed allowed a better understanding of the radiative and non-radiative processes associated with cellular death. Our approach is generic and can also be applied to other radiopharmaceuticals and cell distributions.

Model-based versus specific dosimetry in diagnostic context: comparison of three dosimetric approaches.

PURPOSE: The dosimetric assessment of novel radiotracers represents a legal requirement in most countries. While the techniques for the computation of internal absorbed dose in a therapeutic context have made huge progresses in recent years, in a diagnostic scenario the absorbed dose is usually extracted from model-based lookup tables, most often derived from International Commission on Radiological Protection (ICRP) or Medical Internal Radiation Dose (MIRD) Committee models. The level of approximation introduced by these models may impact the resulting dosimetry. The aim of this work is to establish whether a more refined approach to dosimetry can be implemented in nuclear medicine diagnostics, by analyzing a specific case. METHODS: The authors calculated absorbed doses to various organs in six healthy volunteers administered with flutemetamol ((18)F) injection. Each patient underwent from 8 to 10 whole body 3D PET/CT scans. This dataset was analyzed using a Monte Carlo (MC) application developed in-house using the toolkit gate that is capable to take into account patient-specific anatomy and radiotracer distribution at the voxel level. They compared the absorbed doses obtained with GATE to those calculated with two commercially available software: OLINDA/EXM and STRATOS implementing a dose voxel kernel convolution approach. RESULTS: Absorbed doses calculated with gate were higher than those calculated with OLINDA. The average ratio between gate absorbed doses and OLINDA's was 1.38 ± 0.34 σ (from 0.93 to 2.23). The discrepancy was particularly high for the thyroid, with an average GATE/OLINDA ratio of 1.97 ± 0.83 σ for the six patients. Differences between STRATOS and GATE were found to be higher. The average ratio between GATE and STRATOS absorbed doses was 2.51 ± 1.21 σ (from 1.09 to 6.06). CONCLUSIONS: This study demonstrates how the choice of the absorbed dose calculation algorithm may introduce a bias when gamma radiations are of importance, as is the case in nuclear medicine diagnostics.

Multi-scale hybrid models for radiopharmaceutical dosimetry with Geant4.

The accuracy of radiopharmaceutical absorbed dose distributions computed through Monte Carlo (MC) simulations is mostly limited by the low spatial resolution of 3D imaging techniques used to define the simulation geometry. This issue also persists with the implementation of realistic hybrid models built using polygonal mesh and/or NURBS as they require to be simulated in their voxel form in order to reduce computation times. The existing trade-off between voxel size and simulation speed leads on one side, in an overestimation of the size of small radiosensitive structures such as the skin or hollow organs walls and, on the other, to unnecessarily detailed voxelization of large, homogeneous structures.We developed a set of computational tools based on VTK and Geant4 in order to build multi-resolution organ models. Our aim is to use different voxel sizes to represent anatomical regions of different clinical relevance: the MC implementation of these models is expected to improve spatial resolution in specific anatomical structures without significantly affecting simulation speed. Here we present the tools developed through a proof of principle example. Our approach is validated against the standard Geant4 technique for the simulation of voxel geometries.

Dose point kernels in liquid water: an intra-comparison between GEANT4-DNA and a variety of Monte Carlo codes.

Modeling the radio-induced effects in biological medium still requires accurate physics models to describe the interactions induced by all the charged particles present in the irradiated medium in detail. These interactions include inelastic as well as elastic processes. To check the accuracy of the very low energy models recently implemented into the GEANT4 toolkit for modeling the electron slowing-down in liquid water, the simulation of electron dose point kernels remains the preferential test. In this context, we here report normalized radial dose profiles, for mono-energetic point sources, computed in liquid water by using the very low energy "GEANT4-DNA" physics processes available in the GEANT4 toolkit. In the present study, we report an extensive intra-comparison of profiles obtained by a large selection of existing and well-documented Monte-Carlo codes, namely, EGSnrc, PENELOPE, CPA100, FLUKA and MCNPX.

Hybrid microPET imaging for dosimetric applications in mice: improvement of activity quantification in dynamic microPET imaging for accelerated dosimetry applied to 6-[18 F]fluoro-L-DOPA and 2-[18 F]fluoro-L-tyrosine.

PURPOSE: Dynamic microPET imaging has advantages over traditional organ harvesting, but is prone to quantification errors in small volumes. Hybrid imaging, where microPET activities are cross-calibrated using post scan harvested organs, can improve quantification. Organ harvesting, dynamic imaging and hybrid imaging were applied to determine the human and mouse radiation dosimetry of 6-[18 F]fluoro-L-DOPA and 2-[18 F]fluoro-L-tyrosine and compared. PROCEDURES: Two-hour dynamic microPET imaging was performed with both tracers in four separate mice for 18 F-FDOPA and three mice for 18 F-FTYR. Organ harvesting was performed at 2, 5, 10, 30, 60 and 120 min post tracer injection with n = 5 at each time point for 18 F-FDOPA and n = 3 at each time point for 18 F-FTYR. Human radiation dosimetry projected from animal data was calculated for the three different approaches for each tracer using OLINDA/EXM. S-factors for the MOBY phantom were used to calculate the animal dosimetry. RESULTS: Correlations between dose estimates based on organ harvesting and imaging was improved from r = 0.997 to r = 0.999 for 18 F-FDOPA and from r = 0.985 to r = 0.996 (p < 0.0001 for all) for 18 F-FTYR by using hybrid imaging. CONCLUSION: Hybrid imaging yields comparable results to traditional organ harvesting while partially overcoming the limitations of pure imaging. It is an advantageous technique in terms of number of animals needed and labour involved.

The assessment and management of risks associated with exposures to short-range Auger- and beta-emitting radionuclides. State of the art and proposals for lines of research.

The assessment and management of risks associated with exposures to ionising radiation are defined by the general radiological protection system, proposed by the International Commission on Radiological Protection (ICRP). This system is regarded by a large majority of users as a robust system although there are a number of dissenting voices, claiming that it is not suitable for estimating the risks resulting from internal exposures. One of the specific issues of internal exposure involves short-range radiations such as Auger and beta particles. Auger- and beta-emitting radionuclides can be distributed preferentially in certain tissue structures and even in certain cellular organelles, according to their chemical nature and the vector with which they are associated. Given the limited range of the low-energy electrons in biological matter, this heterogeneous distribution can generate highly localised energy depositions and exacerbate radiotoxic responses at cellular level. These particularities in energy distribution and cellular responses are not taken into account by the conventional methods for the assessment of risk.Alternative systems have been proposed, based on dosimetry conducted at the cellular or even molecular level, whose purpose is to determine the energy deposition occurring within the DNA molecule. However, calculation of absorbed doses at the molecular level is not sufficient to ensure a better assessment of the risks incurred. Favouring such a microdosimetric approach for the risk assessments would require a comprehensive knowledge of the biological targets of radiation, the dose-response relationships at the various levels of organisation, and the mechanisms leading from cellular energy deposition to the appearance of a health detriment. The required knowledge is not fully available today and it is not yet possible to link an intracellular energy deposition to a probability of occurrence of health effects or to use methods based on cellular dosimetry directly.The imperfections of the alternative approaches proposed so far should not discourage efforts. Protection against exposure to Auger and low-energy beta emitters would benefit from mechanistic studies, dedicated to the study of energy depositions of the radionuclides in various cellular structures, but also from radiotoxicological studies to define the relative biological effectiveness of the various Auger emitters used in medicine and of certain low-energy beta emitters, whose behaviour may depend greatly on their chemical form during intake. The scientific expertise, as well as the human and physical resources needed to conduct these studies, is available. They could be now mobilised into international low-dose research programmes, in order to ultimately improve the protection of people exposed to these specific radionuclides.

Comparisons of dosimetric approaches for fractionated radioimmunotherapy of non-Hodgkin lymphoma.

AIM: The aim of this study was to compare different dosimetric approaches on therapy naïve patients enrolled in a multicentre fractionated radioimmunotherapy trial, to determine which methodological approach correlates with bone marrow toxicity. METHODS: Twenty-height non-Hodgkin lymphoma patients were treated with one or two fractions of 90Y-Ibritumomab-Tiuxetan (11.1 MBq/kg) 8 to 12 weeks apart in four different institutions. Quantitative imaging with 111In-Ibritumomab-Tiuxetan (185 MBq) was performed at 0, 1, 4 and 7 days after infusion, starting two weeks before the therapeutic administration. A whole-body (WB) CT scan was also acquired prior to the 111In-Ibritumomab injection, for attenuation correction purposes and was segmented to derive patient-specific organ masses. All dosimetry processing was centralized in a single institution. The first method (M_2D) was based on geometric mean WB scans, corrected for attenuation, scatter and organs superposition. The second method (M_2.5D) was based on the computed assisted matrix inversion approach and used segmented CT scans. The third method (M_3D) used iterative reconstruction of tomographic scans, corrected for attenuation, scatter and collimator response. Absorbed doses were estimated for lungs, liver, kidneys and spleen using MIRD S values adjusted for organ masses. Bone marrow (BM) absorbed doses were evaluated according to imaging methods (3) and compared to blood-based approaches. RESULTS: For some patients, organ masses such as liver or spleen significantly differed from male/female reference masses, whereas lungs and kidneys masses were relatively constant. Except for lungs, absorbed doses estimated by M_2D were higher than those from M_2.5D and these, in turn, were higher that those calculated from M_3D (Wilcoxon P<8.6e-4). Median organ absorbed dose estimates were equivalent for both fractions except for the spleen. In fact, spleen absorbed doses for the second fraction were lower than those for the first fraction, regardless of the approach. Possible explanations are that patient spleen masses were kept constant for analysis of both fractions and/or that spleen uptake was lowered after the first fraction. Estimation of BM absorbed doses from blood sampling was unable to predict platelet toxicity, but image-based methods performed better. Additionally, for most organs, the absorbed dose delivered by the first fraction could predict that delivered by the second fraction. CONCLUSION: These results confirm that different acquisition/processing protocols will lead to statistically different absorbed doses. Additionally, image-based dosimetric approaches are needed in order to correlate absorbed dose to bone marrow toxicity.

EANM Dosimetry Committee guidance document: good practice of clinical dosimetry reporting.

Many recent publications in nuclear medicine contain data on dosimetric findings for existing and new diagnostic and therapeutic agents. In many of these articles, however, a description of the methodology applied for dosimetry is lacking or important details are omitted. The intention of the EANM Dosimetry Committee is to guide the reader through a series of suggestions for reporting dosimetric approaches. The authors are aware of the large amount of data required to report the way a given clinical dosimetry procedure was implemented. Another aim of this guidance document is to provide comprehensive information for preparing and submitting publications and reports containing data on internal dosimetry. This guidance document also contains a checklist which could be useful for reviewers of manuscripts submitted to scientific journals or for grant applications. In addition, this document could be used to decide which data are useful for a documentation of dosimetry results in individual patient records. This may be of importance when the approval of a new radiopharmaceutical by official bodies such as EMA or FDA is envisaged.