A methodology for numerically analysing the hepatic artery haemodynamics during B-TACE: a proof of concept

Two-dimensional perfusion angiography permits direct visualization of redistribution of flow in hepatocellular carcinoma during b-TACE

OBJECTIVES

To demonstrate in vivo redistribution of the blood flow towards HCC's lesions by utilizing two-dimensional perfusion angiography in b-TACE procedures.

MATERIAL AND METHODS

In total, 30 patients with 35 HCC nodules treated in the period between January 2019 and November 2021. For each patient, a post-processing software leading to a two-dimensional perfusion angiography was applied on each angiography performed via balloon microcatheter, before and after inflation. On the colour map obtained, reflecting the evolution of contrast intensity change over time, five regions of interests (ROIs) were assessed: one on the tumour (ROI-t), two in the immediate peritumoural healthy liver parenchyma (ROI-ihl) and two in the peripheral healthy liver parenchyma (ROI-phl). The results have been interpreted with a novel in silico model that simulates the hemodynamics of the hepatic arterial system.

RESULTS

Among the ROIs drawn inside the same segment of target lesion, the time-to-peak of the ROI-t and of the ROI-ihl have a significantly higher mean value when the balloon was inflated compared with the ROIs obtained with deflated balloon (10.33 ± 3.66 s vs 8.87 ± 2.60 s (p = 0.015) for ROI-t; 10.50 ± 3.65 s vs 9.23 ± 2.70 s (p = 0.047) for ROI-ihl). The in silico model prediction time-to-peak delays when balloon was inflated, match with those observed in vivo. The numerical flow analysis shows how time-to-peak delays are caused by the obstruction of the balloon-occluded artery and the opening of intra-hepatic collateral.

CONCLUSION

The measurements identify predictively the flow redistribution in the hepatic arteries during b-TACE, supporting a proper positioning of the balloon microcatheter.

In silico modeling for the hepatic circulation and transport: From the liver organ to lobules

The function of the liver depends critically on its blood supply. Numerous in silico models have been developed to study various aspects of the hepatic circulation, including not only the macro-hemodynamics at the organ level, but also the microcirculation at the lobular level. In addition, computational models of blood flow and bile flow have been used to study the transport, metabolism, and clearance of drugs in pharmacokinetic studies. These in silico models aim to provide insights into the liver organ function under both healthy and diseased states, and to assist quantitative analysis for surgical planning and postsurgery treatment. The purpose of this review is to provide an update on state-of-the-art in silico models of the hepatic circulation and transport processes. We introduce the numerical methods and the physiological background of these models. We also discuss multiscale frameworks that have been proposed for the liver, and their linkage with the large context of systems biology, systems pharmacology, and the Physiome project. This article is categorized under: Metabolic Diseases > Computational Models Metabolic Diseases > Biomedical Engineering Cardiovascular Diseases > Computational Models.

In Vitro Model for Simulating Drug Delivery during Balloon-Occluded Transarterial Chemoembolization

Background: Balloon-occluded transarterial chemoembolization (B-TACE) has emerged as a safe and effective procedure for patients with liver cancer, which is one of the deadliest types of cancer worldwide. B-TACE consist of the transcatheter intraarterial infusion of chemotherapeutic agents, followed by embolizing particles, and it is performed with a microballoon catheter that temporarily occludes a hepatic artery. B-TACE relies on the blood flow redistribution promoted by the balloon-occlusion. However, flow redistribution phenomenon is not yet well understood. Methods: This study aims to present a simple in vitro model (IVM) where B-TACE can be simulated. Results: By visually analyzing the results of various clinically-realistic experiments, the IVM allows for the understanding of balloon-occlusion-related hemodynamic changes and the importance of the occlusion site. Conclusion: The IVM can be used as an educational tool to help clinicians better understand B-TACE treatments. This IVM could also serve as a base for a more sophisticated IVM to be used as a research tool.

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.

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.