Towards a patient-specific hepatic arterial modeling for microspheres distribution optimization in SIRT protocol

Exploring software navigation tools for liver tumour angiography: a scoping review

INTRODUCTION

Liver cancer presents a growing global health concern, necessitating advanced approaches for intervention. This review investigates the use and effectiveness of software navigation in interventional radiology for liver tumour procedures.

METHODS

In accordance with Preferred Reporting Items for Systematic reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR) guidelines, a scoping review was conducted of the literature published between 2013 and 2023 sourcing articles through MEDLINE, Scopus, CINAHL and Embase. Eligible studies focused on liver cancer, utilised cone-beam computed tomography (CBCT), and employed software for intervention. Twenty-one articles were deemed eligible for data extraction and analysis.

RESULTS

Categorised by type, software applications yielded diverse benefits. Feeder detection software significantly enhanced vessel identification, reducing non-target embolisation by up to 43%. Motion correction software demonstrated a 20% enhancement in image quality, effectively mitigating breathing-induced motion artefacts. Liver perfusion software facilitated efficient tumour targeting while simultaneously reducing the occurrence of side effects. Needle guide software enabled precise radiofrequency ablation needle placement. Additionally, these software applications provided detailed anatomical simulations. Overall, software integration resulted in shorter procedures, reduced radiation exposure and decreased contrast media usage.

CONCLUSION

This scoping review highlights the innovative yet relatively underexplored role of software navigation for liver tumour procedures. The integration of software applications not only enhances procedural efficiency but also bolsters operator confidence, and contributes to improved patient outcomes. Despite the current lack of uniformity and standardisation, these software-driven advancements hold significant promise for transforming liver tumour interventions. To realise these benefits, further research is needed to explore the clinical impact and optimal utilisation of software navigation tools in interventional radiology.

Spatiotemporal Analysis of Particle Spread to Assess the Hybrid Particle-Flow CFD Model of Radioembolization of HCC Tumors

OBJECTIVE

Computational fluid dynamics (CFD) models can potentially aid in pre-operative planning of transarterial radioactive microparticle injections to treat hepatocellular carcinoma, but these models are computationally very costly. Previously, we introduced the hybrid particle-flow model as a surrogate, less costly modelling approach for the full particle distribution in truncated hepatic arterial trees. We hypothesized that higher cross-sectional particle spread could increase the match between flow and particle distribution. Here, we investigate whether truncation is still reliable for selective injection scenarios, and if spread is an important factor to consider for reliable truncation.

METHODS

Moderate and severe up- and downstream truncation for selective injection served as input for the hybrid model to compare downstream particle distributions with non-truncated models. In each simulation, particle cross-sectional spread was quantified for 5-6 planes.

RESULTS

Severe truncation gave maximum differences in particle distribution of ∼4-11% and ∼8-9% for down- and upstream truncation, respectively. For moderate truncation, these differences were only ∼1-1.5% and ∼0.5-2%. Considering all particles, spread increased downstream of the tip to 80-90%. However, spread was found to be much lower at specific timepoints, indicating high time-dependency.

CONCLUSION

Combining domain truncation with hybrid particle-flow modelling is an effective method to reduce computational complexity, but moderate truncation is more reliable than severe truncation. Time-dependent spread measures show where differences might arise between flow and particle modelling.

SIGNIFICANCE

The hybrid particle-flow model cuts down computational time significantly by reducing the physical domain, paving the way towards future clinical applications.

Development of a scatter correction technique for planar 99mTc-MAA imaging to improve accuracy in lung shunt fraction estimation

PURPOSE

Prior to ⁹⁰Y selective internal radiation therapy (SIRT) treatment, ^99mTc-MAA scintigraphy imaging is used in the estimation of the lung shunt fraction (LSF). Planar imaging is recommended for determining a LSF ratio. However, the estimate may be affected by scatter contributions, attenuation and respiratory motion. The objective of this study was to correct for the effects of scatter in the LSF, towards the determination of a more accurate estimation method of LSF derived from planar scintigraphy imaging, which is recommended by international guidelines.

METHODS

The open access SIMIND Monte Carlo modelling software was used to estimate an optimum scatter window (SW) for scatter correction. The uncertainties associated with scatter and scatter contributions from the liver on the LSF were evaluated using an anthropomorphic thorax phantom and a virtual Vox-Man phantom. A brief retrospective examination of patient scans and tumour location investigated the impact that the inclusion of the simulated scatter corrections had on the LSF estimation.

RESULTS

The percentage overestimation of the manufacturer recommended method of LSF estimation was 192%. SW corrections improved the uncertainty to within 19% for the range of known LSFs. Similar findings were observed for our patient and tumour location studies.

CONCLUSION

The incorporated scatter corrections can significantly improve the accuracy of the LSF estimation, thereby providing a robust gamma camera, patient and tumour depth specific correction which is easily implementable. This is supported by Monte Carlo, phantom and preliminary patient studies.

The Impact of Injection Distance to Bifurcations on Yttrium-90 Distribution in Liver Cancer Radioembolization

PURPOSE

To model the effect of the injection location on the distribution of yttrium-90 (⁹⁰Y) microspheres in the liver during radioembolization using computational simulation and to determine the potential effects of radial movements of the catheter tip.

MATERIALS AND METHODS

Numerical studies were conducted using images from a representative patient with hepatocellular carcinoma. The right hepatic artery (RHA) was segmented from contrast-enhanced cone-beam computed tomography scans. The blood flow was investigated in the trunk of the RHA using numerical simulations for 6 injection position scenarios at 2 sites located at a distance of approximately 5 and 20 mm upstream of the first bifurcation (RHA diameters of approximately 4.6 mm). The ⁹⁰Y delivery to downstream vessels was calculated from the simulated hepatic artery hemodynamics.

RESULTS

Varying the injection location along the RHA and across the vessel cross-section resulted in different simulated microsphere distributions in the downstream vascular bed. When the catheter tip was 5 mm upstream of the bifurcation, ⁹⁰Y distribution in the downstream branches varied by as much as 53% with a 1.5-mm radial movement of the tip. However, the catheter radial movement had a weaker effect on the microsphere distribution when the injection plane was farther from the first bifurcation (20 mm), with a maximum delivery variation of 9% to a downstream branch.

CONCLUSIONS

An injection location far from bifurcations is recommended to minimize the effect of radial movements of the catheter tip on the microsphere distribution.

Computational Fluid Dynamics Modeling of Liver Radioembolization: A Review

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.

A new software tool for planning interventional procedures in liver cancer

INTRODUCTION

Intra-arterial therapy is an effective way of performing chemotherapy or radiation therapy in patients with primary liver cancer (i.e. hepatocellular carcinoma). Although this minimally invasive approach is now an established treatment option, support tools for pre-operative planning and intra-operative assistance might be helpful.

MATERIAL AND METHODS

We developed an approach for semi-automatic segmentation of computed tomography angiography images of the main arterial branches (required for access path to the treatment site), automatic segmentation of the liver, arterial and venous tree, and interactive segmentation of the tumors (required for procedure-specific planning). This approach was then integrated into a liver-specific workflow within EndoSize^® solution, a planning software for endovascular procedures. The main branches extraction approach was qualitatively evaluated inside the software, while the automatic segmentation methods were quantitatively assessed.

RESULTS

Main branches extraction provides a success rate of 85% (i.e. all arteries correctly extracted) in a dataset of 172 patients. On public databases, a mean DICE of 0.91, 0.47 and 0.92 was obtained for liver, venous and arterial trees segmentation, respectively.

CONCLUSIONS

This pipeline is suitable for directly accessing the treatment site, giving anatomic measurements, and visualizing the hepatic trees, liver, and surrounding arteries during the pre-operative planning.

ABBREVIATIONS

HCC: hepatocellular carcinoma; TACE: transarterial chemoembolization; SIRT: selective internal radiation therapy; CT: computed tomography; CTA: computed tomography angiography; AMS: superior mesenteric artery; LGA: left gastric artery; RHA: right hepatic artery; LHA: left hepatic artery; rbHA: right branch of the hepatic artery; lbHA: left branch of the hepatic artery; GDA: gastroduodenal artery; VOI: volume of interest; SD: standard deviation; MICCAI: medical image computing and computer assisted interventions; MR: magnetic resonance.

A proof-of-concept study of the in-vivo validation of a computational fluid dynamics model of personalized radioembolization

Radioembolization (RE) with yttrium-90 (90Y) microspheres, a transcatheter intraarterial therapy for patients with liver cancer, can be modeled computationally. The purpose of this work was to correlate the results obtained with this methodology using in vivo data, so that this computational tool could be used for the optimization of the RE procedure. The hepatic artery three-dimensional (3D) hemodynamics and microsphere distribution during RE were modeled for six 90Y-loaded microsphere infusions in three patients with hepatocellular carcinoma using a commercially available computational fluid dynamics (CFD) software package. The model was built based on in vivo data acquired during the pretreatment stage. The results of the simulations were compared with the in vivo distribution assessed by 90Y PET/CT. Specifically, the microsphere distribution predicted was compared with the actual 90Y activity per liver segment with a commercially available 3D-voxel dosimetry software (PLANET Dose, DOSIsoft). The average difference between the CFD-based and the PET/CT-based activity distribution was 2.36 percentage points for Patient 1, 3.51 percentage points for Patient 2 and 2.02 percentage points for Patient 3. These results suggest that CFD simulations may help to predict 90Y-microsphere distribution after RE and could be used to optimize the RE procedure on a patient-specific basis.

Multiscale Computational Fluid Dynamics Modeling for Personalized Liver Cancer Radioembolization Dosimetry

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.

Transarterial drug delivery for liver cancer: numerical simulations and experimental validation of particle distribution in patient-specific livers

Background: Transarterial therapies are routinely used for the locoregional treatment of unresectable hepatocellular carcinoma (HCC). However, the impact of clinical parameters (i.e. injection location, particle size, particle density etc.) and patient-specific conditions (i.e. hepatic geometry, cancer burden) on the intrahepatic particle distribution (PD) after transarterial injection of embolizing microparticles is still unclear. Computational fluid dynamics (CFD) may help to better understand this impact.Methods: Using CFD, both the blood flow and microparticle mass transport were modeled throughout the 3D-reconstructed arterial vasculature of a patient-specific healthy and cirrhotic liver. An experimental feasibility study was performed to simulate the PD in a 3D-printed phantom of the cirrhotic arterial network.Results: Axial and in-plane injection locations were shown to be effective parameters to steer particles toward tumor tissue in both geometries. Increasing particle size or density made it more difficult for particles to exit the domain. As cancer burden increased, the catheter tip location mattered less. The in vitro study and numerical results confirmed that PD largely mimics flow distribution, but that significant differences are still possible.Conclusions: Our findings highlight that optimal parameter choice can lead to selective targeting of tumor tissue, but that targeting potential highly depends on patient-specific conditions.

Roudsari B, Vu CT, Roncali E. Computational Modeling of the Liver Arterial Blood Flow for Microsphere Therapy: Effect of Boundary Conditions

Transarterial embolization is a minimally invasive treatment for advanced liver cancer using microspheres loaded with a chemotherapeutic drug or radioactive yttrium-90 (90Y) that are injected into the hepatic arterial tree through a catheter. For personalized treatment, the microsphere distribution in the liver should be optimized through the injection volume and location. Computational fluid dynamics (CFD) simulations of the blood flow in the hepatic artery can help estimate this distribution if carefully parameterized. An important aspect is the choice of the boundary conditions imposed at the inlet and outlets of the computational domain. In this study, the effect of boundary conditions on the hepatic arterial tree hemodynamics was investigated. The outlet boundary conditions were modeled with three-element Windkessel circuits, representative of the downstream vasculature resistance. Results demonstrated that the downstream vasculature resistance affected the hepatic artery hemodynamics such as the velocity field, the pressure field and the blood flow streamline trajectories. Moreover, the number of microspheres received by the tumor significantly changed (more than 10% of the total injected microspheres) with downstream resistance variations. These findings suggest that patient-specific boundary conditions should be used in order to achieve a more accurate drug distribution estimation with CFD in transarterial embolization treatment planning.

Personalized Dosimetry for Liver Cancer Y-90 Radioembolization Using Computational Fluid Dynamics and Monte Carlo Simulation

Yttrium-90 (Y-90) transarterial radioembolization uses radioactive microspheres injected into the hepatic artery to irradiate liver tumors internally. One of the major challenges is the lack of reliable dosimetry methods for dose prediction and dose verification. We present a patient-specific dosimetry approach for personalized treatment planning based on computational fluid dynamics (CFD) simulations of the microsphere transport combined with Y-90 physics modeling called CFDose. The ultimate goal is the development of a software to optimize the amount of activity and injection point for optimal tumor targeting. We present the proof-of-concept of a CFD dosimetry tool based on a patient's angiogram performed in standard-of-care planning. The hepatic arterial tree of the patient was segmented from the cone-beam CT (CBCT) to predict the microsphere transport using multiscale CFD modeling. To calculate the dose distribution, the predicted microsphere distribution was convolved with a Y-90 dose point kernel. Vessels as small as 0.45 mm were segmented, the microsphere distribution between the liver segments using flow analysis was predicted, the volumetric microsphere and resulting dose distribution in the liver volume were computed. The patient was imaged with positron emission tomography (PET) 2 h after radioembolization to evaluate the Y-90 distribution. The dose distribution was found to be consistent with the Y-90 PET images. These results demonstrate the feasibility of developing a complete framework for personalized Y-90 microsphere simulation and dosimetry using patient-specific input parameters.