Personalized Dosimetry for Liver Cancer Y-90 Radioembolization Using Computational Fluid Dynamics and Monte Carlo Simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Roncali E
PROVIDER: S-EPMC7160004 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Yttrium-90 (90Y) radioembolization is a minimally invasive procedure increasingly used for advanced liver cancer treatment. In this method, radioactive microspheres are injected into the hepatic arterial bloodstream to target, irradiate, and kill cancer cells. Accurate and precise treatment planning can lead to more efficient and safer treatment by delivering a higher radiation dose to the tumor while minimizing the exposure of the surrounding liver parenchyma. Treatment planning primarily relies on the estimated radiation dose delivered to tissue. However, current methods used to estimate the dose are based on simplified assumptions that make the dosimetry results unreliable. In this work, we present a computational model to predict the radiation dose from the 90Y activity in different liver segments to provide a more realistic and personalized dosimetry. Computational fluid dynamics (CFD) simulations were performed in a 3D hepatic arterial tree model segmented from cone-beam CT angiographic data obtained from a patient with hepatocellular carcinoma (HCC). The microsphere trajectories were predicted from the velocity field. 90Y dose distribution was then calculated from the volumetric distribution of the microspheres. Two injection locations were considered for the microsphere administration, a lobar and a selective injection. Results showed that 22% and 82% of the microspheres were delivered to the tumor, after each injection, respectively, and the combination of both injections ultimately delivered 49% of the total administered 90Y microspheres to the tumor. Results also illustrated the nonhomogeneous distribution of microspheres between liver segments, indicating the importance of developing patient-specific dosimetry methods for effective radioembolization treatment.
| S-EPMC7580665 | biostudies-literature
Project description:BACKGROUND:Intrahepatic dosimetry is paramount to optimize radioembolization treatment accuracy using radioactive holmium-166 microspheres (166Ho). This requires a practical protocol that combines quantitative imaging of microsphere distribution with automated and robust delineation of the volumes of interest. To this end, we propose a dual isotope single photon emission computed tomography (SPECT) protocol based on 166Ho therapeutic microspheres and technetium-99?m (99mTc) stannous phytate, which accumulates in healthy liver tissue. This protocol may allow accurate and automatic estimation of tumor-absorbed dose and healthy liver-absorbed dose. The current study focuses on a Monte Carlo-based reconstruction framework that inherently corrects for scatter crosstalk between the 166Ho and 99mTc imaging. To demonstrate the feasibility of the method, it is evaluated with realistic phantom experiments and patient data. METHODS:The Utrecht Monte Carlo System (UMCS) was extended to include detailed modeling of crosstalk interactions between 99mTc and 166Ho. First, 99mTc images were reconstructed including energy window-based corrections for 166Ho downscatter. Next, 99mTc downscatter in the 81-keV 166Ho window was Monte Carlo simulated to allow quantitative reconstruction of the 166Ho images. The accuracy of the 99mTc-downscatter modeling was evaluated by comparing measurements with simulations. In addition, the ratio between 99mTc and 166Ho yielding the best 166Ho dose estimates was established and the quantitative accuracy was reported. RESULTS:Given the same level of activity, 99mTc contributes twice as many counts to the 81-keV window than 166Ho, and four times as many counts to the 140-keV window, applying a 166Ho/99mTc ratio of 5:1 yielded a high accuracy in both 166Ho and 99mTc reconstruction. Phantom experiments revealed that the accuracy of quantitative 166Ho activity recovery was reduced by 10% due to the presence of 99mTc. Twenty iterations (8 subsets) of the SPECT/CT reconstructions were considered feasible for clinical practice. Applicability of the proposed protocol was shown in a proof-of-concept case. CONCLUSION:A novel 166Ho/99mTc dual-isotope protocol for automatic dosimetry compensates accurately for downscatter and allows for the addition of 99mTc without compromising 166Ho SPECT image quality.
| S-EPMC7056760 | biostudies-literature
Project description:In recent years, yttrium-90 ((90)Y) microsphere radioembolization has been establishing itself as a safe and efficacious treatment for both primary and metastatic liver cancers. This extends to both first-line therapies as well as in the salvage setting. In addition, radioembolization appears efficacious for patients with portal vein thrombosis, which is currently a contraindication for surgery, transplantation and transarterial chemoembolization. This article reviews the efficacy and expanding use of (90)Y microsphere radioembolization with an added emphasis on recent advances in personalized dosimetry and interventional radiology techniques. Directions for future research into combination therapies with radioembolization and expansion into sites other than the liver are also explored.
| S-EPMC5258157 | biostudies-literature
Project description:Yttrium-90 radioembolization (RE) is a widely used transcatheter intraarterial therapy for patients with unresectable liver cancer. In the last decade, computer simulations of hepatic artery hemodynamics during RE have been performed with the aim of better understanding and improving the therapy. In this review, we introduce the concept of computational fluid dynamics (CFD) modeling with a clinical perspective and we review the CFD models used to study RE from the fluid mechanics point of view. Finally, we show what CFD simulations have taught us about the hemodynamics during RE, the current capabilities of CFD simulations of RE, and we suggest some future perspectives.
| S-EPMC8716346 | biostudies-literature
Project description:Mesh-type and voxel-based computational phantoms comprise the current state of the art for internal dose assessment via Monte Carlo simulations but excel in different aspects, with mesh-type phantoms offering advantages over their voxel counterparts in terms of their flexibility and realistic representation of detailed patient- or subject-specific anatomy. We have developed PARaDIM (pronounced "paradigm": Particle and Heavy Ion Transport Code System-Based Application for Radionuclide Dosimetry in Meshes), a freeware application for implementing tetrahedral mesh-type phantoms in absorbed dose calculations. It considers all medically relevant radionuclides, including α, β, γ, positron, and Auger/conversion electron emitters, and handles calculation of mean dose to individual regions, as well as 3-dimensional dose distributions for visualization and analysis in a variety of medical imaging software. This work describes the development of PARaDIM, documents the measures taken to test and validate its performance, and presents examples of its uses. Methods: Human, small-animal, and cell-level dose calculations were performed with PARaDIM and the results compared with those of widely accepted dosimetry programs and literature data. Several tetrahedral phantoms were developed or adapted using computer-aided modeling techniques for these comparisons. Results: For human dose calculations, agreement of PARaDIM with OLINDA 2.0 was good-within 10%-20% for most organs-despite geometric differences among the phantoms tested. Agreement with MIRDcell for cell-level S value calculations was within 5% in most cases. Conclusion: PARaDIM extends the use of Monte Carlo dose calculations to the broader community in nuclear medicine by providing a user-friendly graphical user interface for calculation setup and execution. PARaDIM leverages the enhanced anatomic realism provided by advanced computational reference phantoms or bespoke image-derived phantoms to enable improved assessments of radiation doses in a variety of radiopharmaceutical use cases, research, and preclinical development. PARaDIM can be downloaded freely at www.paradim-dose.org.
| S-EPMC6894378 | biostudies-literature
Project description:A model based on a specific phantom, called QuAArC, has been designed for the evaluation of planning and verification systems of complex radiotherapy treatments, such as volumetric modulated arc therapy (VMAT). This model uses the high accuracy provided by the Monte Carlo (MC) simulation of log files and allows the experimental feedback from the high spatial resolution of films hosted in QuAArC. This cylindrical phantom was specifically designed to host films rolled at different radial distances able to take into account the entrance fluence and the 3D dose distribution. Ionization chamber measurements are also included in the feedback process for absolute dose considerations. In this way, automated MC simulation of treatment log files is implemented to calculate the actual delivery geometries, while the monitor units are experimentally adjusted to reconstruct the dose-volume histogram (DVH) on the patient CT. Prostate and head and neck clinical cases, previously planned with Monaco and Pinnacle treatment planning systems and verified with two different commercial systems (Delta4 and COMPASS), were selected in order to test operational feasibility of the proposed model. The proper operation of the feedback procedure was proved through the achieved high agreement between reconstructed dose distributions and the film measurements (global gamma passing rates > 90% for the 2%/2 mm criteria). The necessary discretization level of the log file for dose calculation and the potential mismatching between calculated control points and detection grid in the verification process were discussed. Besides the effect of dose calculation accuracy of the analytic algorithm implemented in treatment planning systems for a dynamic technique, it was discussed the importance of the detection density level and its location in VMAT specific phantom to obtain a more reliable DVH in the patient CT. The proposed model also showed enough robustness and efficiency to be considered as a pre-treatment VMAT verification system.
| S-EPMC5117721 | biostudies-literature
Project description:This work is a comparative study of the dosimetry calculated by Plaque Simulator, a treatment planning system for eye plaque brachytherapy, to the dosimetry calculated using Monte Carlo simulation for an Eye Physics model EP917 eye plaque. Monte Carlo (MC) simulation using MCNPX 2.7 was used to calculate the central axis dose in water for an EP917 eye plaque fully loaded with 17 IsoAid Advantage (125)I seeds. In addition, the dosimetry parameters ?, gL(r), and F(r,?) were calculated for the IsoAid Advantage model IAI-125 (125)I seed and benchmarked against published data. Bebig Plaque Simulator (PS) v5.74 was used to calculate the central axis dose based on the AAPM Updated Task Group 43 (TG-43U1) dose formalism. The calculated central axis dose from MC and PS was then compared. When the MC dosimetry parameters for the IsoAid Advantage (125)I seed were compared with the consensus values, ? agreed with the consensus value to within 2.3%. However, much larger differences were found between MC calculated gL(r) and F(r,?) and the consensus values. The differences between MC-calculated dosimetry parameters are much smaller when compared with recently published data. The differences between the calculated central axis absolute dose from MC and PS ranged from 5% to 10% for distances between 1 and 12 mm from the outer scleral surface. When the dosimetry parameters for the (125)I seed from this study were used in PS, the calculated absolute central axis dose differences were reduced by 2.3% from depths of 4 to 12 mm from the outer scleral surface. We conclude that PS adequately models the central dose profile of this plaque using its defaults for the IsoAid model IAI-125 at distances of 1 to 7 mm from the outer scleral surface. However, improved dose accuracy can be obtained by using updated dosimetry parameters for the IsoAid model IAI-125 (125)I seed.
| S-EPMC5691011 | biostudies-literature
Project description:BackgroundThe purpose was to validate 90Y PET gradient-based tumor segmentation in phantoms and to evaluate the impact of the segmentation method on reported tumor absorbed dose (AD) and biological effective dose (BED) in 90Y microsphere radioembolization (RE) patients. A semi-automated gradient-based method was applied to phantoms and patient tumors on the 90Y PET with the initial bounding volume for gradient detection determined from a registered diagnostic CT or MR; this PET-based segmentation (PS) was compared with radiologist-defined morphologic segmentation (MS) on CT or MRI. AD and BED volume histogram metrics (D90, D70, mean) were calculated using both segmentations and concordance/correlations were investigated. Spatial concordance was assessed using Dice similarity coefficient (DSC) and mean distance to agreement (MDA). PS was repeated to assess intra-observer variability.ResultsIn phantoms, PS demonstrated high accuracy in lesion volumes (within 15%), AD metrics (within 11%), high spatial concordance relative to morphologic segmentation (DSC?>?0.86 and MDA?<?1.5 mm), and low intra-observer variability (DSC?>?0.99, MDA?<?0.2 mm, AD/BED metrics within 2%). For patients (58 lesions), spatial concordance between PS and MS was degraded compared to in-phantom (average DSC = 0.54, average MDA = 4.8 mm); the average mean tumor AD was 226?±?153 and 197?±?138 Gy, respectively for PS and MS. For patient AD metrics, the best Pearson correlation (r) and concordance correlation coefficient (ccc) between segmentation methods was found for mean AD (r?=?0.94, ccc?=?0.92), but worsened as the metric approached the minimum dose (for D90, r?=?0.77, ccc?=?0.69); BED metrics exhibited a similar trend. Patient PS showed low intra-observer variability (average DSC = 0.81, average MDA = 2.2 mm, average AD/BED metrics within 3.0%).Conclusions90Y PET gradient-based segmentation led to accurate/robust results in phantoms, and showed high concordance with MS for reporting mean tumor AD/BED in patients. However, tumor coverage metrics such as D90 exhibited worse concordance between segmentation methods, highlighting the need to standardize segmentation methods when reporting AD/BED metrics from post-therapy 90Y PET. Estimated differences in reported AD/BED metrics due to segmentation method will be useful for interpreting RE dosimetry results in the literature including tumor response data.
| S-EPMC6265358 | biostudies-literature
Project description:Conventional mammography can suffer from poor contrast between healthy and cancerous tissues due to the small difference in attenuation properties. Coherent scatter slot scan imaging is an imaging technique which provides additional information and is compatible with conventional mammography. A Monte Carlo simulation of coherent scatter slot scan imaging was performed to assess its performance and provide system optimization. Coherent scatter could be exploited using a system similar to conventional slot scan mammography system with antiscatter grids tilted at the characteristic angle of cancerous tissues. System optimization was performed across several parameters, including source voltage, tilt angle, grid distances, grid ratio, and shielding geometry. The simulated carcinomas were detectable for tumors as small as 5 mm in diameter, so coherent scatter analysis using a wide-slot setup could be promising as an enhancement for screening mammography. Employing coherent scatter information simultaneously with conventional mammography could yield a conventional high spatial resolution image with additional coherent scatter information.
| S-EPMC4999523 | biostudies-other
Project description:IntroductionThe use of high-dose-rate brachytherapy is currently a widespread practice worldwide. The most common isotope source is 192Ir, but 60Co is also becoming available for HDR. One of main advantages of 60Co compared to 192Ir is the economic and practical benefit because of its longer half-live, which is 5.27 years. Recently, Eckert & Ziegler BEBIG, Germany, introduced a new afterloading brachytherapy machine (MultiSource®); it has the option to use either the 60Co or 192Ir HDR source. The source for the Monte Carlo calculations is the new 60Co source (model Co0.A86), which is referred to as the new BEBIG 60Co HDR source and is a modified version of the 60Co source (model GK60M21), which is also from BEBIG.Objective and methodsThe purpose of this work is to obtain the dosimetry parameters in accordance with the AAPM TG-43U1 formalism with Monte Carlo calculations regarding the BEBIG 60Co high-dose-rate brachytherapy to investigate the required treatment-planning parameters. The geometric design and material details of the source was provided by the manufacturer and was used to define the Monte Carlo geometry. To validate the source geometry, a few dosimetry parameters had to be calculated according to the AAPM TG-43U1 formalism. The dosimetry studies included the calculation of the air kerma strength Sk, collision kerma in water along the transverse axis with an unbounded phantom, dose rate constant and radial dose function. The Monte Carlo code system that was used was EGSnrc with a new cavity code, which is a part of EGS++ that allows calculating the radial dose function around the source. The spectrum to simulate 60Co was composed of two photon energies, 1.17 and 1.33 MeV. Only the gamma part of the spectrum was used; the contribution of the electrons to the dose is negligible because of the full absorption by the stainless-steel wall around the metallic 60Co. The XCOM photon cross-section library was used in subsequent simulations, and the photoelectric effect, pair production, Rayleigh scattering and bound Compton scattering were included in the simulation. Variance reduction techniques were used to speed up the calculation and to considerably reduce the computer time. The cut-off energy was 10 keV for electrons and photons. To obtain the dose rate distributions of the source in an unbounded liquid water phantom, the source was immersed at the center of a cube phantom of 100 cm3. The liquid water density was 0.998 g/cm3, and photon histories of up to 1010 were used to obtain the results with a standard deviation of less than 0.5% (k = 1). The obtained dose rate constant for the BEBIG 60Co source was 1.108±0.001 cGyh-1U-1, which is consistent with the values in the literature. The radial dose functions were compared with the values of the consensus data set in the literature, and they are consistent with the published data for this energy range.
| S-EPMC4587958 | biostudies-literature