Clinical Positron Emission Tomography Imaging Research Articles

FAST (fast analytical simulator of tracer)-PET: an accurate and efficient PET analytical simulation tool

Objective. Simulation of positron emission tomography (PET) images is an essential tool in the development and validation of quantitative imaging workflows and advanced image processing pipelines. Existing Monte Carlo or analytical PET simulators often compromise on either efficiency or accuracy. We aim to develop and validate fast analytical simulator of tracer (FAST)-PET, a novel analytical framework, to simulate PET images accurately and efficiently. Approach. FAST-PET simulates PET images by performing precise forward projection, scatter, and random estimation that match the scanner geometry and statistics. Although the same process should be applicable to other scanner models, we focus on the Siemens Biograph Vision-600 in this work. Calibration and validation of FAST-PET were performed through comparison with an experimental scan of a National Electrical Manufacturers Association (NEMA) Image Quality (IQ) phantom. Further validation was conducted between FAST-PET and Geant4 Application for Tomographic Emission (GATE) quantitatively in clinical image simulations in terms of intensity-based and texture-based features and task-based tumor segmentation. Main results. According to the NEMA IQ phantom simulation, FAST-PET’s simulated images exhibited partial volume effects and noise levels comparable to experimental images, with a relative bias of the recovery coefficient RC within 10% for all spheres and a coefficient of variation for the background region within 6% across various acquisition times. FAST-PET generated clinical PET images exhibit high quantitative accuracy and texture comparable to GATE (correlation coefficients of all features over 0.95) but with ∼100-fold lower computation time. The tumor segmentation masks comparison between both methods exhibited significant overlap and shape similarity with high concordance CCC > 0.97 across measures. Significance. FAST-PET generated PET images with high quantitative accuracy comparable to GATE, making it ideal for applications requiring extensive PET image simulations such as virtual imaging trials, and the development and validation of image processing pipelines.

Deep denoiser prior driven relaxed iterated Tikhonov method for low-count PET image restoration

Objective. Low-count positron emission tomography (PET) imaging is an efficient way to promote more widespread use of PET because of its short scan time and low injected activity. However, this often leads to low-quality PET images with clinical image reconstruction, due to high noise and blurring effects. Existing PET image restoration (IR) methods hinder their own restoration performance due to the semi-convergence property and the lack of suitable denoiser prior. Approach. To overcome these limitations, we propose a novel deep plug-and-play IR method called Deep denoiser Prior driven Relaxed Iterated Tikhonov method (DP-RI-Tikhonov). Specifically, we train a deep convolutional neural network denoiser to generate a flexible deep denoiser prior to handle high noise. Then, we plug the deep denoiser prior as a modular part into a novel iterative optimization algorithm to handle blurring effects and propose an adaptive parameter selection strategy for the iterative optimization algorithm. Main results. Simulation results show that the deep denoiser prior plays the role of reducing noise intensity, while the novel iterative optimization algorithm and adaptive parameter selection strategy can effectively eliminate the semi-convergence property. They enable DP-RI-Tikhonov to achieve an average quantitative result (normalized root mean square error, structural similarity) of (0.1364, 0.9574) at the stopping iteration, outperforming a conventional PET IR method with an average quantitative result of (0.1533, 0.9523) and a state-of-the-art deep plug-and-play IR method with an average quantitative result of (0.1404, 0.9554). Moreover, the advantage of DP-RI-Tikhonov becomes more obvious at the last iteration. Experiments on six clinical whole-body PET images further indicate that DP-RI-Tikhonov successfully reduces noise intensity and recovers fine details, recovering sharper and more uniform images than the comparison methods. Significance. DP-RI-Tikhonov’s ability to reduce noise intensity and effectively eliminate the semi-convergence property overcomes the limitations of existing methods. This advancement may have substantial implications for other medical IR.

DOCK-PET: database of CNS kinetic parameters in the healthy human brain for existing PET tracers.

Information about developed positron emission tomography (PET) tracers and obtained clinical PET images is publicly available in a database. However, findings regarding the kinetic parameters of PET tracers are yet to be summarized. Therefore, in this study, we created an open-access database of central nervous system (CNS) kinetic parameters in the healthy human brain for existing PET tracers (DOCK-PET). Our database includes information on the kinetic parameters and compounds of existing CNS-PET tracers. The kinetic parameter dataset comprises the analysis methods, VT, BPND, K parameters, relevant literature, and study details. The list of PET tracers and kinetic parameter information was compiled through keyword-based searches of PubMed and the Molecular Imaging and Contrast Agent Database (MICAD). The kinetic parameters obtained, including VT, BPND, and K parameters, were reorganized based on the defined brain anatomical regions. All data were rigorously double-checked before being summarized in Microsoft Excel and JavaScript Object Notation (JSON) formats. Of the 247 PET tracers identified through searches using the PubMed and MICAD websites, the kinetic parameters of 120 PET tracers were available. Among the 120 PET tracers, compound structures with chemical and physical properties were obtained from the PubChem website or the ChemDraw software. Furthermore, the affinity information of the 104 PET tracers was gathered from PubChem or extensive literature surveys of the 120 PET tracers. We developed a comprehensive open-access database, DOCK-PET, that includes both kinetic parameters of healthy humans and compound information for existing CNS-PET tracers.

Fully Automated Cassette-Based Synthesis of 2-Deoxy-2-[18F]Fluorocellobiose Using Trasis AllInOne Module.

Due to the continuous rise in global incidence and severity of invasive fungal infections (IFIs), particularly among immunocompromised and immunodeficient patients, there is an urgent demand for swift and accurate fungal pathogen diagnosis. Therefore, the need for fungal-specific positron emission tomography (PET) imaging agents that can detect the infection in the early stages is increasing. Cellobiose, a disaccharide, is readily metabolized by fungal pathogens such as Aspergillus species. Recently, our group reported fluorine-18 labeled cellobiose, 2-deoxy-2-[18F]fluorocellobiose ([18F]FCB), for specific imaging of Aspergillus infection. The positive imaging findings with very low background signal on delayed imaging make this ligand a promising fungal-specific imaging ligand. Inspired by this result, the decision was made to automate the radiolabeling procedure for better reproducibility and to facilitate clinical translation. A Trasis AllInOne (Trasis AIO) automated module was used for this purpose. The reagent vials contain commercially available 2-deoxy-2-[18F]fluoroglucose ([18F]FDG), glucose-1-phosphate, and enzyme (cellobiose phosphorylase). A Sep-Pak cartridge was used to purify the tracer. The overall radiochemical yield was 50%-70% (n = 6, decay corrected) in 75-min synthesis time with a radiochemical purity of > 98%. This is a highly reliable protocol to produce current good manufacturing practice (cGMP)-compliant [18F]FCB for clinical PET imaging.

Learning Without Real Data Annotations to Detect Hepatic Lesions in PET Images.

Deep neural networks have been recently applied to lesion identification in fluorodeoxyglucose (FDG) positron emission tomography (PET) images, but they typically rely on a large amount of well-annotated data for model training. This is extremely difficult to achieve for neuroendocrine tumors (NETs), because of low incidence of NETs and expensive lesion annotation in PET images. The objective of this study is to design a novel, adaptable deep learning method, which uses no real lesion annotations but instead low-cost, list mode-simulated data, for hepatic lesion detection in real-world clinical NET PET images. We first propose a region-guided generative adversarial network (RG-GAN) for lesion-preserved image-to-image translation. Then, we design a specific data augmentation module for our list-mode simulated data and incorporate this module into the RG-GAN to improve model training. Finally, we combine the RG-GAN, the data augmentation module and a lesion detection neural network into a unified framework for joint-task learning to adaptatively identify lesions in real-world PET data. The proposed method outperforms recent state-of-the-art lesion detection methods in real clinical 68Ga-DOTATATE PET images, and produces very competitive performance with the target model that is trained with real lesion annotations. With RG-GAN modeling and specific data augmentation, we can obtain good lesion detection performance without using any real data annotations. This study introduces an adaptable deep learning method for hepatic lesion identification in NETs, which can significantly reduce human effort for data annotation and improve model generalizability for lesion detection with PET imaging.

Sensitive Positron Emission Tomography Imaging of PD-L1 Expression in Human Breast and Lung Carcinoma Xenografts Using the Radiometalated Peptide Ga-68-TRAP-WL12.

Noninvasive imaging of the immune checkpoint protein programmed death ligand 1 (PD-L1; synonyms: CD274, B7-H1) holds great promise to improve patient selection and, thus, response rates for immune checkpoint therapy (ICT) with monoclonal antibodies targeting the PD1/PD-L1 axis. The PD-L1 specific peptide WL12 (cyclo(AcY-(NMe)A-N-P-H-L-Hyp-W-S-W(Me)-(NMe)Nle-(NMe)Nle-O-C)-G-NH2) was functionalized with the Gallium-68 chelator TRAP by means of click chemistry (CuAAC). The resulting conjugate TRAP-WL12 was labeled with Gallium-68 within 16 min, with approximately 90% radiochemical yield and 99% radiochemical purity, affording Ga-68-TRAP-WL12 with molar activities typically exceeding 100 MBq/nmol. This radiotracer was characterized by positron emission tomography (PET) imaging and ex vivo biodistribution in murine xenografts of nontransfected PD-L1 expressing tumor cell lines, MDA-MB-231 (human breast carcinoma), and H2009 (human lung adenocarcinoma). It showed a favorable biodistribution profile with rapid renal clearance and low background (tumor-to-blood ratio = 26.6, 3 h p.i.). Conjugation of the Ga-68-TRAP moiety to WL12 successfully mitigated the nonspecific uptake of this peptide in organs, particularly the liver. This was demonstrated by comparing Ga-68-TRAP-WL12 with the archetypical Ga-68-DOTA-WL12, for which tumor-to-liver ratios of 1.4 and 0.5, respectively, were found. Although immunohistochemistry (IHC) revealed a low PD-L1 expression in MDA-MB-213 and H2009 xenografts that corresponds well to the clinical situation, PET showed high tumor uptakes (6.6 and 7.3% injected activity per gram of tissue (iA/g), respectively) for Ga-68-TRAP-WL12. Thus, this tracer has the potential for routine clinical PD-L1 PET imaging because it detects even very low PD-L1 expression densities with high sensitivity and may open an avenue to replace PD-L1 IHC of biopsies as the standard means to select potential responders for ICT.

Preclinical Evaluation and a Pilot Clinical Positron Emission Tomography Imaging Study of [68Ga]Ga-FAPI-FUSCC-II.

Fibroblast activation protein (FAP), a type II integral membrane serine protease, is a promising target for tumor diagnosis and therapy. OncoFAP has been recently discovered for PET imaging procedures for various solid malignancies. In this study, we presented the development of manual radiolabeling procedures for the preparation of OncoFAP-based radiopharmaceuticals for cancer imaging. A novel series of [68Ga/177Lu]Ga/Lu-FAPI-FUSCC-I/II were produced with high radiochemical yields. [68Ga]Ga-FAPI-FUSCC-I/II and [177Lu]Lu-FAPI-FUSCC-I/II were stable in phosphate-buffered saline, fetal bovine serum, and human serum for at least 3 h. In vitro cellular uptake and blocking experiments implied that they had specificity to FAP. Additionally, the low nanomolar IC50 values of FAPI-FUSCC-II indicated that it had a high target affinity to FAP. The in vivo biodistribution and blocking study in mice bearing HT-1080-FAP tumors showed that both exhibited specific tumor uptake. [68Ga]Ga-FAPI-FUSCC-II showed a higher tumor uptake and a higher tumor/nontarget ratio than [68Ga]Ga-FAPI-FUSCC-I and [68Ga]Ga-FAPI-04. The results of ex vivo biodistribution were in accordance with the biodistribution results. Clinical [68Ga]Ga-FAPI-FUSCC-II-PET/CT imaging further demonstrated its favorable biodistribution and kinetics with elevated and reliable uptake by primary tumors (maximum standardized uptake value (SUVmax), 12.17 ± 6.67) and distant metastases (SUVmax, 9.24 ± 4.28). In summary, [68Ga]Ga-FAPI-FUSCC-II displayed increased tumor uptake and retention compared to [68Ga]Ga-FAPI-04, giving it potential as a promising tracer for the diagnostic imaging of malignant tumors with positive FAP expression.

A multivariate MRI model of Alzheimer’s disease risk is associated with clinical diagnosis, PET imaging, and plasma biomarkers in a mixed dementia sample

AbstractBackgroundWe sought to validate two structural MRI‐based brain health models: Spatial Pattern of Atrophy for REcognition of Alzheimer’s disease (SPARE‐AD) and SPARE‐Brain Age Gap (SPARE‐BAG), which estimates the discrepancy between biological brain age and chronological age. These models' generalization to non‐amnestic AD (naAD) syndromes and associations with positron emission tomography (PET) and plasma biomarkers remain underinvestigated. We hypothesized that SPARE‐AD scores in naAD would be lower (less abnormal) than in amnestic AD (aAD) but elevated (more abnormal) relative to controls; and that SPARE‐AD but not SPARE‐BAG scores would be associated with PET and plasma biomarkers of AD.MethodSPARE‐AD and SPARE‐BAG scores were estimated for each of 1607 MRI scans from 277 cognitively normal adults (62.8% female; mean age = 69.5 years) and 902 clinical participants (54.2% female; mean age = 73.1 years), including 585 with amnestic mild cognitive impairment (aMCI) or aAD; 56 with naAD syndromes; 28 with vascular disease; 28 with Lewy body disorders (LBD); 18 with suspected frontotemporal lobar degeneration (FTLD); 56 with non‐amnestic MCI of uncertain etiology; and 131 with other clinical diagnoses. We used mixed effects models adjusting for age, sex, and global cognition to test group differences in SPARE‐AD and SPARE‐BAG. Additionally, we assessed associations with tau (n = 122) and amyloid‐beta (n = 199) PET imaging and plasma‐based phosphorylated tau‐181 (p‐tau181) and glial fibrillary acidic protein (GFAP) (n = 358).ResultSPARE‐AD scores were elevated vs. controls in aMCI/aAD, naAD, vascular disease, FTLD, and naMCI (all p<0.003) but not in LBD or other diagnoses (p>0.1). NaAD patients had lower SPARE‐AD (p<0.02) but higher SPARE‐BAG scores (p<0.02) than aMCI/aAD patients. SPARE‐AD scores were positively associated with mean cortical tau PET uptake, amyloid PET positivity, and plasma p‐tau181 and GFAP (all p<0.001). In contrast, SPARE‐BAG scores were not associated with PET or plasma biomarkers (all p>0.05).ConclusionSPARE‐AD values for naAD syndromes were intermediate between those of controls and aMCI/AD patients, likely reflecting divergent, partially overlapping cortical atrophy relative to aAD. PET and plasma biomarkers of AD neuropathologic change were associated with SPARE‐AD but not with SPARE‐BAG scores, suggesting SPARE‐BAG reflects age‐related atrophy that is largely independent of AD.

Astrocyte Signature in Alzheimer's Disease Continuum through a Multi-PET Tracer Imaging Perspective.

Reactive astrogliosis is an early event in the continuum of Alzheimer's disease (AD). Current advances in positron emission tomography (PET) imaging provide ways of assessing reactive astrogliosis in the living brain. In this review, we revisit clinical PET imaging and in vitro findings using the multi-tracer approach, and point out that reactive astrogliosis precedes the deposition of Aβ plaques, tau pathology, and neurodegeneration in AD. Furthermore, considering the current view of reactive astrogliosis heterogeneity-more than one subtype of astrocyte involved-in AD, we discuss how astrocytic body fluid biomarkers might fit into trajectories different from that of astrocytic PET imaging. Future research focusing on the development of innovative astrocytic PET radiotracers and fluid biomarkers may provide further insights into the heterogeneity of reactive astrogliosis and improve the detection of AD in its early stages.

Thermodynamic and Kinetic Stabilities of Al(III) Complexes with N2O3 Pentadentate Ligands.

Al(III) complexes have been recently investigated for their potential use in imaging with positron emission tomography (PET) by formation of ternary complexes with the radioisotope fluorine-18 (18F). Although the derivatives of 1,4,7-triazacyclononane-1,4,7-triacetic acid (NOTA) are the most applied chelators for [Al18F]2+ labelling and (pre)clinical PET imaging, non-macrocyclic, semi-rigid pentadentate chelators having two N- and three O-donor atoms such as RESCA1 and AMPDA-HB have been proposed with the aim to allow room temperature labelling of temperature-sensitive biomolecules. The paucity of stability data on Al(III) complexes used for PET imaging instigated a complete thermodynamic and kinetic solution study on Al(III) complexes with aminomethylpiperidine (AMP) derivatives AMPTA and AMPDA-HB and the comparison with a RESCA1-like chelator CD3A-Bn (trans-1,2-diaminocyclohexane-N-benzyl-N,N',N'-triacetic acid). The stability constant of [Al(AMPDA-HB)] is about four orders of magnitude higher than that of [Al(AMPTA)] and [Al(CD3A-Bn)], highlighting the greater affinity of phenolates with respect to acetate O-donors. On the other hand, the kinetic inertness of the complexes, determined by following the Cu2+-mediated transmetallation reactions in the 7.5-10.5 pH range, resulted in a spontaneous and hydroxide-assisted dissociation slightly faster for [Al(AMPTA)] than for the other two complexes (t1/2 = 4.5 h for [Al(AMPTA)], 12.4 h for [Al(AMPDA-HB)], and 24.1 h for [Al(CD3A-Bn)] at pH 7.4 and 25 °C). Finally, the [AlF]2+ ternary complexes were prepared and their stability in reconstituted human serum was determined by 19F NMR experiments.

Experimental study of medical isotopes 62,64Cu and 68Ga production using intense picosecond laser pulse

Radioisotopes are indispensable agents in medical diagnosis and treatment, among which copper-62, 64 (62,64Cu) and gallium-68 (68Ga) are medical isotopes widely used in positron emission tomography (PET) imaging. Experiments that generate these radioisotopes via laser-induced photonuclear reactions were performed on the XingGuangIII laser facility of the Laser Fusion Research Center in Mianyang. High-charge ( ∼ 40 nC) MeV electron (e− ) beams were generated with 100 terawatts, picosecond laser pulses. The e− beams were then impinged on a metal stack composed of Ta foil and activation plates (natural Cu and Ga2O3), producing high-energy bremsstrahlung x-rays and medical isotopes 62,64Cu and 68Ga, respectively. The characteristic emissions of the produced 62,64Cu and 68Ga were detected off-line , and the production yields of 62,64Cu and 68Ga were obtained to be of the order of 106 per laser shot. For electrons with energy higher than 8 MeV, the dependence of isotope production efficiency (per e− ) on electron temperature () is investigated through Geant4 simulations. It is found that the production efficiency increases with the . At ∼ 10 MeV, the values are 10−4 for 62Cu and 2 × 10 − 5 for 68Ga, respectively. The prospect of producing the medical isotopes 62,64Cu and 68Ga are further evaluated using a table-top femtosecond laser system of high repetition. With a repetition rate of 100 Hz, their activity is expected to reach 0.2 GBq for 62Cu, 0.1 GBq for 64Cu and 0.05 GBq for 68Ga, respectively. Such activity would meet the required dose for clinical PET imaging, indicating the great potential to produce medical isotopes with an all-optical, high-repetition laser system.

Publicly available framework for simulating and experimentally validating clinical PET systems.

Monte Carlo (MC) simulations are a powerful tool to model medical imaging systems. However, before simulations can be considered the ground truth, they have to be validated with experiments. To provide a pipeline that models a clinical positron emission tomography (PET)/CT system using MC simulations after extensively validating the results against experimental measurements. A clinical four-ring PET imaging system was modeled using Geant4 application for tomographic emission (v. 9.0). To validate the simulations, PET images were acquired of a cylindrical phantom, point source, and image quality phantom with the modeled system and the simulations of the experimental procedures. For the purpose of validating the quantification capabilities and image quality provided by the simulation pipeline, the simulations were compared against the measurements in terms of their count rates and sensitivity as well as their image uniformity, resolution, recovery coefficients (RCs), coefficients of variation, contrast, and background variability. When compared to the measured data, the number of true detections in the MC simulations was within 5%. The scatter fraction was found to be 30.0%±2.2% and 28.8%±1.7% in the measured and simulated scans, respectively. Analyzing the measured and simulated sinograms, the sensitivities were found to be 8.2 and 7.8cps/kBq, respectively. The fraction of random coincidences were 19% in the measured data and 25% in the simulation. When calculating the image uniformity within the axial slices, the measured image exhibited a uniformity of 0.015±0.005, whereas the simulated image had a uniformity of 0.029±0.011. In the axial direction, the uniformity was measured to be 0.024±0.006 and 0.040±0.015 for the measured and simulated data, respectively. Comparing the image resolution, an average percentage difference of 2.9% was found between the measurements and simulations. The RCs calculated in both the measured and simulated images were found to be within the EARL ranges, except for that of the simulation of the smallest sphere. The coefficients of variation for the measured and simulated images were found to be 12% and 13%, respectively. Lastly, the background variability was consistent between the measurements and simulations, whereas the average percentage difference in the sphere contrasts was found to be 8.8%. The clinical PET/CT system was modeled and validated to provide a simulation pipeline for the community. The pipeline and the validation procedures have been made available (https://github.com/teaghan/PET_MonteCarlo).

Evaluation of the performance of a high-resolution time-of-flight PET system dedicated to the head and breast according to NEMA NU 2-2012 standard

BackgroundThis study evaluated the physical performance of a positron emission tomography (PET) system dedicated to the head and breast according to the National Electrical Manufacturers Association (NEMA) NU2-2012 standard.MethodsThe spatial resolution, sensitivity, scatter fraction, count rate characteristics, corrections for count losses and randoms, and image quality of the system were determined. All measurements were performed according to the NEMA NU2-2012 acquisition protocols, but image quality was assessed using a brain-sized phantom. Furthermore, scans of the three-dimensional (3D) Hoffmann brain phantom and mini-Derenzo phantom were acquired to allow visual evaluation of the imaging performance for small structures.ResultsThe tangential, radial, and axial full width at half maximum (FWHM) at a 10-mm offset in half the axial field of view were measured as 2.3, 2.5, and 2.9 mm, respectively. The average system sensitivity at the center of the field of view and at a 10-cm radial offset was 7.18 and 8.65 cps/kBq, respectively. The peak noise-equivalent counting rate was 35.2 kcps at 4.8 kBq/ml. The corresponding scatter fraction at the peak noise-equivalent counting rate was 46.8%. The peak true rate and scatter fraction at 8.6 kBq/ml were 127.8 kcps and 54.3%, respectively. The percent contrast value for a 10-mm sphere was approximately 50%. On the 3D Hoffman brain phantom image, the structures of the thin layers composing the phantom were visualized on the sagittal and coronal images. On the mini-Derenzo phantom, each of the 1.6-mm rods was clearly visualized.ConclusionTaken together, these results indicate that the head- and breast-dedicated PET system has high resolution and is well suited for clinical PET imaging.

Hypometabolic and hypermetabolic brain regions in patients with ALS-FTD show distinct patterns of grey and white matter degeneration: A pilot multimodal neuroimaging study

BackgroundUp to 50% of amyotrophic lateral sclerosis (ALS) patients develop some degree of cognitive dysfunction and a small proportion of these develop frontotemporal dementia (FTD). Non-invasive techniques of magnetic resonance imaging (MRI) and [18F]-fluoro-2-deoxy-d-glucose (18F-FDG) positron emission tomography (PET) have demonstrated structural and metabolic abnormalities, respectively, in the brains of such patients with ALS-FTD. Although initial 18F-FDG PET studies in ALS patients showed only hypometabolism of motor and extramotor brain regions, subsequent studies have demonstrated hypermetabolic changes as well. Such contrasting findings prompted us to hypothesize that hypo- and hypermetabolic brain regions in ALS-FTD patients are associated with divergent degeneration of structural grey matter (GM) and white matter (WM). MethodsCerebral glucose metabolic rate (CMRglc), cortical thickness (CT), fractal dimension (FD), and graph theory WM network analyses were performed on clinical MRI and 18F-FDG PET images from 8 ALS-FTD patients and 14 neurologic controls to explore the relationship between GM-WM degeneration and hypo- and hypermetabolic brain regions. ResultsCMRglc revealed significant hypometabolism in frontal and precentral gyrus brain regions, with hypermetabolism in temporal, occipital and cerebellar regions. Cortical thinning was noted in both hypo- and hypermetabolic brain areas. Unlike CT, FD did not reveal widespread GM degeneration in hypo- and hypermetabolic brain regions of ALS-FTD patients. Graph theory analysis showed severe WM degeneration in hypometabolic but not hypermetabolic areas, especially in the right hemisphere. ConclusionOur multimodal MRI-PET study provides insights into potentially differential pathophysiological mechanisms between hypo- and hypermetabolic brain regions of ALS-FTD patients.

Fast beta-emitter Monte Carlo simulations and full patient dose calculations of targeted radionuclide therapy: introducing egs_mird.

Targeted radionuclide therapy (TRT) is a fast-growing field garnering much interest, with several clinical trials currently underway, that has a steady increase in development of treatment techniques. Unfortunately, within the field and many clinical trials, the dosimetry calculation techniques used remain relatively simple, often using a mix of S-value calculations and kernelconvolutions. The common TRT calculation techniques, although very quick, can often ignore important aspects of patient anatomy and radionuclide distribution, as well as the interplay there-in. This paper introduces egs_mird, a new Monte Carlo (MC) application built in EGSnrc which allows users to model full patient tissue and density (using clinical CT images) and radionuclide distribution (using clinical PET images) for fast and detailed dose TRT calculation. The novel application egs_mird is introduced along with a general outline of the structure of egs_mird simulations. The general structure of the code, and the track-length (TL) estimator scoring implementation for variance reduction, is described. A new egs++ source class egs_internal_source, created to allow detailed patient-wide source distribution, and a modified version of egs_radionuclide_source, changed to be able to work with egs_internal_source, are also described. The new code is compared to other MC calculations of S-values kernels of 131 I, 90 Y, and 177 Lu in the literature, along with further self-validation using a histogram variant of the electron Fano test. Several full patient prostate 177 Lu TRT prostate cancer treatment simulations are performed using a single set of patient DICOM CT and [18 F]-DCFPyL PETdata. Good agreement is found between S-value kernels calculated using egs_mird with egs_internal_source and those found in the literature. Calculating 1000 doses (individual voxel uncertainties of ∼0.05%) in a voxel grid Fano test for monoenergetic 500keV electrons and 177 Lu electrons results in 94% and 99% of the doses being within 0.1% of the expected dose, respectively. For a hypothetical 177 Lu treatment, patient prostate, rectum, bone marrow, and bladder dose volume histograms (DVHs) results did not vary significantly when using the TL estimator and not modeling electron transport, modeling bone marrow explicitly (rather than using generic tissue compositions), and reducing activity to voxels containing partial or full calcifications to half or none, respectively. Dose profiles through different regions demonstrate there are some differences with model choices not seen in the DVH. Simulations using the TL estimator can be completed in under 15 min (∼30 min when using standard interaction scoring). This work shows egs_mird to be a reliable MC code for computing TRT doses as realistically as the patient Computed Tomography (CT) and Positron Emission Tomography (PET) data allow. Furthermore, the code can compute doses to sub-1% uncertainty within 15 min, with little to no optimization. Thus, this work supports the use of egs_mird for dose calculations inTRT.

Combined multi-scale mesh and full-matrix inversion for enhancing time-domain breast diffuse optical tomography.

Time-domain diffuse optical tomography can efficiently reconstruct both absorption and reduced scattering coefficients but is heavily limited by the ill-posedness in its inverse problem and low spatial resolution. To deal with these adversities, the truncated singular value decomposition (TSVD)-based whole-weighting-matrix inversion scheme can be a particularly suitable implementation. Unfortunately, TSVD is subject to a storage challenge for three-dimensional imaging of a bulk region, such as breast. In this paper, a multi-scale mesh strategy based on computed tomography (CT) anatomical geometry is adopted to solve the storage challenge, where a fine mesh is used in forward calculation to ensure accuracy, and a coarse mesh in the inversion process to enable TSVD-based inversion of the whole-weighting matrix. We validate the proposed strategy using simulated data for a single lesion model from clinical positron emission tomography images of a breast cancer patient, and further, for a complex model that is constructed by setting dual lesions at different separations in the CT breast geometry.

Cycle-consistent learning-based hybrid iterative reconstruction for whole-body PET imaging

Objective. To develop a cycle-consistent learning-based hybrid iterative reconstruction (IR) method that takes only slightly longer than analytic reconstruction, while pursuing the image resolution and tumor quantification achievable by IR for whole-body PET imaging. Approach. We backproject the raw positron emission tomography (PET) data to generate a blurred activity distribution. From the backprojection to the IR label, a reconstruction mapping that approximates the deblurring filters for the point spread function and the physical effects of the PET system is unrolled to a neural network with stacked convolutional layers. By minimizing the cycle-consistent loss, we train the reconstruction and inverse mappings simultaneously. Main results. In phantom study, the proposed method results in an absolute relative error (RE) of the mean activity of 4.0% ± 0.7% in the largest hot sphere, similar to the RE of the full-count IR and significantly smaller than that obtained by CycleGAN postprocessing. Achieving a noise reduction of 48.1% ± 0.5% relative to the low-count IR, the proposed method demonstrates advantages over the low-count IR and CycleGAN in terms of resolution maintenance, contrast recovery, and noise reduction. In patient study, the proposed method obtains a noise reduction of 44.6% ± 8.0% for the lung and the liver, while maintaining the regional mean activity in both simulated lesions and real tumors. The run time of the proposed method is only half that of the conventional IR. Significance. The proposed cycle-consistent learning from the backprojection rather than the raw PET data or an IR result enables improved reconstruction accuracy, reduced memory requirements, and fast implementation speeds for clinical whole-body PET imaging.

Synthesis of 68Ga-radiopharmaceuticals using both generator-derived and cyclotron-produced 68Ga as exemplified by [68Ga]Ga-PSMA-11 for prostate cancer PET imaging.

[68Ga]Ga-PSMA-11, a urea-based peptidomimetic, is a diagnostic radiopharmaceutical for positron emission tomography (PET) imaging that targets the prostate-specific membrane antigen (PSMA). The recent Food and Drug Administration approval of [68Ga]Ga-PSMA-11 for PET imaging of patients with prostate cancer, expected follow-up approval of companion radiotherapeutics (e.g., [177Lu]Lu-PSMA-617, [225Ac]Ac-PSMA-617) and large prostate cancer patient volumes requiring access are poised to create an unprecedented demand for [68Ga]Ga-PSMA-11 in nuclear medicine clinics around the world. Meeting this global demand is going to require a variety of synthesis methods compatible with 68Ga eluted from a generator or produced on a cyclotron. To address this urgent need in the PET radiochemistry community, herein we report detailed protocols for the synthesis of [68Ga]Ga-PSMA-11, (also known as HBED-CC, Glu-urea-Lys(Ahx)-HBED-CC and PSMA-HBED-CC) using both generator-eluted and cyclotron-produced 68Ga and contrast the pros and cons of each method. The radiosyntheses are automated and have been validated for human use at two sites (University of Michigan (UM), United States; Royal Prince Alfred Hospital (RPA), Australia) and used to produce [68Ga]Ga-PSMA-11 for patient use in good activity yields (single generator, 0.52 GBq (14 mCi); dual generators, 1.04-1.57 GBq (28-42 mCi); cyclotron method (single target), 1.47-1.89 GBq (40-51 mCi); cyclotron method (dual target), 3.63 GBq (98 mCi)) and high radiochemical purity (99%) (UM, n = 645; RPA, n > 600). Both methods are appropriate for clinical production but, in the long term, the method employing cyclotron-produced 68Ga is the most promising for meeting high patient volumes. Quality control testing (visual inspection, pH, radiochemical purity and identity, radionuclidic purity and identity, sterile filter integrity, bacterial endotoxin content, sterility, stability) confirmed doses are suitable for clinical use, and there is no difference in clinical prostate cancer PET imaging using [68Ga]Ga-PSMA-11 prepared using the two production methods.

Clinical evaluation of [18F]pitavastatin for quantitative analysis of hepatobiliary transporter activity

It is widely accepted that uptake and efflux transporters on clearance organs play crucial roles in drug disposition. Although in vitro transporter assay system can identify the intrinsic properties of the target transporters, it is not so easy to precisely predict in vivo pharmacokinetic parameters from in vitro data. Positron emission tomography (PET) imaging is a useful tool to directly assess the activity of drug transporters in humans. We recently developed a practical synthetic method for fluorine-18-labeled pitavastatin ([18F]PTV) as a PET probe for quantitative evaluation of hepatobiliary transport. In the present study, we conducted clinical PET imaging with [18F]PTV and compared the pharmacokinetic properties of the probe for healthy subjects with or without rifampicin pretreatment. Rifampicin pretreatment significantly suppressed the hepatic maximum concentration and biliary excretion of the probe to 52% and 34% of the control values, respectively. Rifampicin treatment markedly decreased hepatic uptake clearance (21% of the control), and moderately canalicular efflux clearance with regard to hepatic concentration (52% of the control). These results demonstrate that [18F]PTV is a useful probe for clinical investigation of the activities of hepatobiliary uptake/efflux transporters in humans.

Optimizing Immuno-PET Imaging of Tumor PD-L1 Expression: Pharmacokinetic, Biodistribution, and Dosimetric Comparisons of 89Zr-Labeled Anti-PD-L1 Antibody Formats.

PET imaging of programmed cell death ligand 1 (PD-L1) may help to noninvasively predict and monitor responses to anti-programmed cell death 1/anti-PD-L1 immunotherapies. In this study, we compared the imaging characteristics of 3 radioligands derived from the anti-PD-L1 IgG1 complement 4 (C4). In addition to the IgG C4, we produced a fragment antigen-binding (Fab) C4, as well as a double-mutant IgG C4 (H310A/H435Q) with minimal affinity for the murine neonatal Fc receptor. Methods: The pharmacokinetics, biodistribution, and dosimetry of the 3 89Zr-labeled C4 ligands were compared by longitudinal PET/CT imaging in nude mice bearing subcutaneous human non-small cell lung cancer xenografts with positive (H1975 model) or negative (A549 model) endogenous PD-L1 expression. Results: The C4 radioligands substantially accumulated in PD-L1-positive tumors but not in PD-L1-negative tumors or in blocked PD-L1-positive tumors, confirming their PD-L1-specific tumor targeting. 89Zr-Fab C4 and 89Zr-IgG C4 (H310A/H435Q) were rapidly eliminated compared with 89Zr-IgG C4. Consequently, maximal tumor-to-muscle ratios were obtained earlier, at 4 h after injection for 89Zr-Fab C4 (ratio, ∼6) and 24 h after injection for 89Zr-IgG C4 (H310A/H435Q) (ratio, ∼9), versus 48 h after injection for 89Zr-IgG C4 (ratio, ∼8). Background activity in nontumor tissues was low, except for high kidney retention of 89Zr-Fab C4 and persistent liver accumulation of 89Zr-IgG C4 (H310A/H435Q) compared with 89Zr-IgG C4. Dosimetry estimates suggested that the C4 radioligands would yield organ-absorbed doses tolerable for repeated clinical PET imaging studies. Conclusion: This study highlights the potential of designing radioligands with shorter pharmacokinetics for PD-L1 immuno-PET imaging in a preclinical model and encourages further clinical translation of such radioligands.