Impact of fluid-structure interaction on direct tumor-targeting in a representative hepatic artery system. Ann Biomed Eng. 2014;42:461-474. [PMID: 24048712 DOI: 10.1007/s10439-013-0910-7]

Improved magnetic drug targeting with maximized magnetic forces and limited particle spreading

BACKGROUND

In magnetic drug targeting (MDT), micro- or nanoparticles are injected into the human body to locally deliver therapeutics. These magnetic particles can be guided from a distance by external magnetic fields and gradients from electromagnets.

PURPOSE

During the particles' movement through the vascular network, they are affected by magnetic forces, fluid (drag) forces, particle interactions, diffusion, etc. Adequate targeting is hindered when drag forces overcome the magnetic forces and particles present in vessels are carried away from the targeted region. Moreover, the magnetic force directions and diffusion mechanisms can cause particles to scatter, while they should remain together for an effective targeting performance. In this work, these adverse effects are tackled using optimization methods.

METHODS

We formulate an optimization problem with respect to the currents in surrounding electromagnets that aims to maximize the magnetic force on a particle along a predefined direction. A boundary on the magnetic force divergence is introduced as a constraint to limit particle spreading. We also consider particles to be moved from an initial to a target location in a finite-time interval. To this end dynamic optimization is applied.

RESULTS

Simulations for particles in a bifurcated vessel show an increase of particle speed by 20% and a successful movement towards the targeted regions without spreading. For the dynamic optimization, simulation results demonstrate that particle collections are accurately guided with 10 times less scattering and 10 times more particles at the target than without the divergence constraint.

CONCLUSIONS

The proposed methods significantly improve the steering and capturing of particles in a region of interest. They are applicable to any magnetic drug targeting configuration with electromagnets.

A Hybrid Particle-Flow CFD Modeling Approach in Truncated Hepatic Arterial Trees for Liver Radioembolization: A Patient-specific Case Study

Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer. At its intermediate, unresectable stage, HCC is typically treated by local injection of embolizing microspheres in the hepatic arteries to selectively damage tumor tissue. Interestingly, computational fluid dynamics (CFD) has been applied increasingly to elucidate the impact of clinically variable parameters, such as injection location, on the downstream particle distribution. This study aims to reduce the computational cost of such CFD approaches by introducing a novel truncation algorithm to simplify hepatic arterial trees, and a hybrid particle-flow modeling approach which only models particles in the first few bifurcations. A patient-specific hepatic arterial geometry was pruned at three different levels, resulting in three trees: Geometry 1 (48 outlets), Geometry 2 (38 outlets), and Geometry 3 (17 outlets). In each geometry, 1 planar injection and 3 catheter injections (each with different tip locations) were performed. For the truncated geometries, it was assumed that, downstream of the truncated outlets, particles distributed themselves proportional to the blood flow. This allowed to compare the particle distribution in all 48 "outlets" for each geometry. For the planar injections, the median difference in outlet-specific particle distribution between Geometry 1 and 3 was 0.21%; while the median difference between outlet-specific flow and particle distribution in Geometry 1 was 0.40%. Comparing catheter injections, the maximum median difference in particle distribution between Geometry 1 and 3 was 0.24%, while the maximum median difference between particle and flow distribution was 0.62%. The results suggest that the hepatic arterial tree might be reliably truncated to estimate the particle distribution in the full-complexity tree. In the resulting hybrid particle-flow model, explicit particle modeling was only deemed necessary in the first few bifurcations of the arterial tree. Interestingly, using flow distribution as a surrogate for particle distribution in the entire tree was considerably less accurate than using the hybrid model, although the difference was much higher for catheter injections than for planar injections. Future work should focus on replicating and experimentally validating these results in more patient-specific geometries.

Particle Distribution in Embolotherapy, How Do They Get There? A Critical Review of the Factors Affecting Arterial Distribution of Embolic Particles

Embolization has tremendously evolved in recent years and has expanded to treatment of a variety of pathologic processes. There has been emerging evidence that the level of arterial occlusion and the distribution of embolic particles may play an important role in the clinical outcome. This is a comprehensive literature review to identify variables that play important role in determination of level of occlusion of blood vessels and distribution of embolic particles. The literature searches between 1996 to 2020 through PubMed and Ovid-MEDLINE yielded over 1030 articles of which 30 studies providing details on the level of occlusion are reviewed here. We divided the playing factors into characteristics of the particles, solution/injection and vascular bed. Accordingly, particle size, type and aggregation, compressibility/deformability, and biodegradability are categorized as the factors involving particles' behavioral nature. Infusion rate and concentration/dilution of the medium are related to the carrying solution. Hemodynamics and the arterial resistance are characteristics of the vascular bed that also play an important role in the distribution of embolic particles. Understanding and predicting the level of embolization is a complex multi-factor problem that requires more evidence, warranting further randomized controlled trials, and powered human and animal studies.

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.

CFD Simulations of Radioembolization: A Proof-of-Concept Study on the Impact of the Hepatic Artery Tree Truncation

Radioembolization (RE) is a treatment for patients with liver cancer, one of the leading cause of cancer-related deaths worldwide. RE consists of the transcatheter intraarterial infusion of radioactive microspheres, which are injected at the hepatic artery level and are transported in the bloodstream, aiming to target tumors and spare healthy liver parenchyma. In paving the way towards a computer platform that allows for a treatment planning based on computational fluid dynamics (CFD) simulations, the current simulation (model preprocess, model solving, model postprocess) times (of the order of days) make the CFD-based assessment non-viable. One of the approaches to reduce the simulation time includes the reduction in size of the simulated truncated hepatic artery. In this study, we analyze for three patient-specific hepatic arteries the impact of reducing the geometry of the hepatic artery on the simulation time. Results show that geometries can be efficiently shortened without impacting greatly on the microsphere distribution.

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.

In Vitro Study of Particle Transport in Successively Bifurcating Vessels

To reach a predictive understanding of how particles travel through bifurcating vessels is of paramount importance in many biomedical settings, including embolization, thromboembolism, and drug delivery. Here we utilize an in vitro model in which solid particles are injected through a rigid vessel that symmetrically bifurcates in successive branching generations. The geometric proportion and fluid dynamics parameters are relevant to the liver embolization. The volumetric flow field is reconstructed via phase-contrast magnetic resonance imaging, from which the particle trajectories are calculated for a range of size and density using the particle equation of motion. The method is validated by directly tracking the injected particles via optical imaging. The results indicate that, opposite to the common assumption, the particles distribution is fundamentally different from the volumetric flow partition. In fact, the amount of delivered particles vary substantially between adjacent branches even when the flow is uniformly distributed. This is not due to the inertia of the particles, nor to gravity. The particle distribution is rather rooted in their different pathways, which in turn are linked to their release origin along the main vessel cross-section. Therefore, the tree geometry and the associated flow streamlines are the prime determinant of the particle fate, while local changes of volumetric flow rate to selected branches do not generally produce proportional changes of particle delivery.

The role of angled-tip microcatheter and microsphere injection velocity in liver radioembolization: A computational particle-hemodynamics study

Liver radioembolization is a promising treatment option for combating liver tumors. It is performed by placing a microcatheter in the hepatic artery and administering radiation-emitting microspheres through the arterial bloodstream so that they get lodged in the tumoral bed. In avoiding nontarget radiation, the standard practice is to conduct a pretreatment, in which the microcatheter location and injection velocity are decided. However, between pretreatment and actual treatment, some of the parameters that influence the particle distribution in the liver can vary, resulting in radiation-induced complications. The present study aims to analyze the influence of a commercially available microcatheter with an angled tip and particle injection velocity in terms of segment-to-segment particle distribution. Specifically, 4 tip orientations and 2 injection velocities are combined to yield a set of 8 numerical simulations of the particle-hemodynamics in a patient-specific truncated hepatic artery. For each simulation, 4 cardiac pulses are simulated. Particles are injected during the first cycle, and the remaining pulses enable the majority of the injected particles to exit the computational domain. Results indicate that, in terms of injection velocity, particles are more spread out in the cross-sectional lumen areas as the injection velocity increases. The tip's orientation also plays a role because it influences the near-tip hemodynamics, therefore altering the particle travel through the hepatic artery. However, results suggest that particle distribution tries to match the blood flow split, therefore particle injection velocity and microcatheter tip orientation playing a minor role in segment-to-segment particle distribution.

Computational particle-haemodynamics analysis of liver radioembolization pretreatment as an actual treatment surrogate

Liver radioembolization (RE) is a treatment option for patients with unresectable and chemorefractory primary and metastatic liver tumours. RE consists of intra-arterially administering via catheter radioactive microspheres that locally attack the tumours, sparing healthy tissue. Prior to RE, the standard practice is to conduct a treatment-mimicking pretreatment assessment via the infusion of ^99m Tc-labelled macroaggregated albumin microparticles. The usefulness of this pretreatment has been debated in the literature, and thus, the aim of the present study is to shed light on this issue by numerically simulating the liver RE pretreatment and actual treatment particle-haemodynamics in a patient-specific hepatic artery under two different literature-based cancer scenarios and two different placements of a realistic end-hole microcatheter in the proper hepatic artery. The parameters that are analysed are the following: microagent quantity and size (accounting for RE pretreatment and treatment), catheter-tip position (near the proper hepatic artery bifurcation and away from it), and cancer burden (10% and 30% liver involvement). The conclusion that can be reached from the simulations is that when it comes to mimicking RE in terms of delivering particles to tumour-bearing segments, the catheter-tip position is much more important (because of the importance of local haemodynamic pattern alteration) than the infused microagents (i.e. quantity and size). Cancer burden is another important feature because the increase in blood flow rate to tumour-bearing segments increases the power to drag particles. These numerical simulation-based conclusions are in agreement with clinically observed events reported in the literature. Copyright © 2016 John Wiley & Sons, Ltd.

Partial hepatectomy hemodynamics changes: Experimental data explained by closed-loop lumped modeling

The liver function may be degraded after partial liver ablation surgery. Adverse liver hemodynamics have been shown to be associated to liver failure. The link between these hemodynamics changes and ablation size is however poorly understood. This article proposes to explain with a closed-loop lumped model the hemodynamics changes observed during twelve surgeries in pigs. The portal venous tree is modeled with a pressure-dependent variable resistor. The variables measured, before liver ablation, are used to tune the model parameters. Then, the liver partial ablation is simulated with the model and the simulated pressures and flows are compared with post-operative measurements. Fluid infusion and blood losses occur during the surgery. The closed-loop model presented accounts for these blood volume changes. Moreover, the impact of blood volume changes and the liver lobe mass estimations on the simulated variables is studied. The typical increase of portal pressure, increase of liver pressure loss, slight decrease of portal flow and major decrease in arterial flow are quantitatively captured by the model for a 75% hepatectomy. It appears that the 75% decrease in hepatic arterial flow can be explained by the resistance increase induced by the surgery, and that no hepatic arterial buffer response (HABR) mechanism is needed to account for this change. The different post-operative states, observed in experiments, are reproduced with the proposed model. Thus, an explanation for inter-subjects post-operative variability is proposed. The presented framework can easily be adapted to other species circulations and to different pathologies for clinical hepatic applications.

Liver cancer arterial perfusion modelling and CFD boundary conditions methodology: a case study of the haemodynamics of a patient-specific hepatic artery in literature-based healthy and tumour-bearing liver scenarios

Some of the latest treatments for unresectable liver malignancies (primary or metastatic tumours), which include bland embolisation, chemoembolisation, and radioembolisation, among others, take advantage of the increased arterial blood supply to the tumours to locally attack them. A better understanding of the factors that influence this transport may help improve the therapeutic procedures by taking advantage of flow patterns or by designing catheters and infusion systems that result in the injected beads having increased access to the tumour vasculature. Computational analyses may help understand the haemodynamic patterns and embolic-microsphere transport through the hepatic arteries. In addition, physiological inflow and outflow boundary conditions are essential in order to reliably represent the blood flow through arteries. This study presents a liver cancer arterial perfusion model based on a literature review and derives boundary conditions for tumour-bearing liver-feeding hepatic arteries based on the arterial perfusion characteristics of normal and tumorous liver segment tissue masses and the hepatic artery branching configuration. Literature-based healthy and tumour-bearing realistic scenarios are created and haemodynamically analysed for the same patient-specific hepatic artery. As a result, this study provides boundary conditions for computational fluid dynamics simulations that will allow researchers to numerically study, for example, various intravascular devices used for liver disease intra-arterial treatments with different cancer scenarios. Copyright © 2016 John Wiley & Sons, Ltd.

Drug-targeting methodologies with applications: A review

Targeted drug delivery to solid tumors is a very active research area, focusing mainly on improved drug formulation and associated best delivery methods/devices. Drug-targeting has the potential to greatly improve drug-delivery efficacy, reduce side effects, and lower the treatment costs. However, the vast majority of drug-targeting studies assume that the drug-particles are already at the target site or at least in its direct vicinity. In this review, drug-delivery methodologies, drug types and drug-delivery devices are discussed with examples in two major application areas: (1) inhaled drug-aerosol delivery into human lung-airways; and (2) intravascular drug-delivery for solid tumor targeting. The major problem addressed is how to deliver efficiently the drug-particles from the entry/infusion point to the target site. So far, most experimental results are based on animal studies. Concerning pulmonary drug delivery, the focus is on the pros and cons of three inhaler types, i.e., pressurized metered dose inhaler, dry powder inhaler and nebulizer, in addition to drug-aerosol formulations. Computational fluid-particle dynamics techniques and the underlying methodology for a smart inhaler system are discussed as well. Concerning intravascular drug-delivery for solid tumor targeting, passive and active targeting are reviewed as well as direct drug-targeting, using optimal delivery of radioactive microspheres to liver tumors as an example. The review concludes with suggestions for future work, considereing both pulmonary drug targeting and direct drug delivery to solid tumors in the vascular system.