A proof-of-concept pipeline to guide evaluation of tumor tissue perfusion by dynamic contrast-agent imaging: Direct simulation and inverse tracer-kinetic procedures.

Dynamic contrast-enhanced (DCE) perfusion imaging has shown great potential to non-invasively assess cancer development and its treatment by their characteristic tissue signatures. Different tracer kinetics models are being applied to estimate tissue and tumor perfusion parameters from DCE perfusion imaging. The goal of this work is to provide an in silico model-based pipeline to evaluate how these DCE imaging parameters may relate to the true tissue parameters. As histology data provides detailed microstructural but not functional parameters, this work can also help to better interpret such data. To this aim in silico vasculatures are constructed and the spread of contrast agent in the tissue is simulated. As a proof of principle we show the evaluation procedure of two tracer kinetic models from in silico contrast-agent perfusion data after a bolus injection. Representative microvascular arterial and venous trees are constructed in silico. Blood flow is computed in the different vessels. Contrast-agent input in the feeding artery, intra-vascular transport, intra-extravascular exchange and diffusion within the interstitial space are modeled. From this spatiotemporal model, intensity maps are computed leading to in silico dynamic perfusion images. Various tumor vascularizations (architecture and function) are studied and show spatiotemporal contrast imaging dynamics characteristic of in vivo tumor morphotypes. The Brix II also called 2CXM, and extended Tofts tracer-kinetics models common in DCE imaging are then applied to recover perfusion parameters that are compared with the ground truth parameters of the in silico spatiotemporal models. The results show that tumor features can be well identified for a certain permeability range. The simulation results in this work indicate that taking into account space explicitly to estimate perfusion parameters may lead to significant improvements in the perfusion interpretation of the current tracer-kinetics models.

Assessing Early Cardiac Outflow Tract Adaptive Responses Through Combined Experimental-Computational Manipulations.

Mechanical forces are essential for proper growth and remodeling of the primitive pharyngeal arch arteries (PAAs) into the great vessels of the heart. Despite general acknowledgement of a hemodynamic-malformation link, the direct correlation between hemodynamics and PAA morphogenesis remains poorly understood. The elusiveness is largely due to difficulty in performing isolated hemodynamic perturbations and quantifying changes in-vivo. Previous in-vivo arch artery occlusion/ablation experiments either did not isolate the effects of hemodynamics, did not analyze the results in a 3D context or did not consider the effects of varying degrees of occlusion. Here, we overcome these limitations by combining minimally invasive occlusion experiments in the avian embryo with 3D anatomical models of development and in-silico testing of experimental phenomenon. We detail morphological and hemodynamic changes 24 hours post vessel occlusion. 3D anatomical models showed that occlusion geometries had more circular cross-sectional areas and more elongated arches than their control counterparts. Computational fluid dynamics revealed a marked change in wall shear stress-morphology trends. Instantaneous (in-silico) occlusion models provided mechanistic insights into the dynamic vessel adaptation process, predicting pressure-area trends for a number of experimental occlusion arches. We follow the propagation of small defects in a single embryo Hamburger Hamilton (HH) Stage 18 embryo to a more serious defect in an HH29 embryo. Results demonstrate that hemodynamic perturbation of the presumptive aortic arch, through varying degrees of vessel occlusion, overrides natural growth mechanisms and prevents it from becoming the dominant arch of the aorta.

A whole lung in silico model to estimate age dependent particle dosimetry.

Anatomical and physiological changes alter airflow characteristics and aerosol distribution in the developing lung. Correlation between age and aerosol dosimetry is needed, specifically because youth are more susceptible to medication side effects. In this study, we estimate aerosol dosages (particle diameters of 1, 3, and 5 [Formula: see text]m) in a 3 month-old infant, a 6 year-old child, and a 36 year-old adult by performing whole lung subject-specific particle simulations throughout respiration. For 3 [Formula: see text]m diameter particles we estimate total deposition as 88, 73, and [Formula: see text] and the conducting versus respiratory deposition ratios as 4.0, 0.5, and 0.4 for the infant, child, and adult, respectively. Due to their lower tidal volumes and functional residual capacities the deposited mass is smaller while the tissue concentrations are larger in the infant and child subjects, compared to the adult. Furthermore, we find that dose cannot be predicted by simply scaling by tidal volumes. These results highlight the need for additional clinical and computational studies that investigate the efficiency of treatment, while optimizing dosage levels in order to alleviate side effects, in youth.

Simulation of a detoxifying organ function: Focus on hemodynamics modeling and convection-reaction numerical simulation in microcirculatory networks.

When modeling a detoxifying organ function, an important component is the impact of flow on the metabolism of a compound of interest carried by the blood. We here study the effects of red blood cells (such as the Fahraeus-Lindqvist effect and plasma skimming) on blood flow in typical microcirculatory components such as tubes, bifurcations and entire networks, with particular emphasis on the liver as important representative of detoxifying organs. In one of the plasma skimming models, under certain conditions, oscillations between states are found and analyzed in a methodical study to identify their causes and influencing parameters. The flow solution obtained is then used to define the velocity at which a compound would be transported. A convection-reaction equation is studied to simulate the transport of a compound in blood and its uptake by the surrounding cells. Different types of signal sharpness have to be handled depending on the application to address different temporal compound concentration profiles. To permit executing the studied models numerically stable and accurate, we here extend existing transport schemes to handle converging bifurcations, and more generally multi-furcations. We study the accuracy of different numerical schemes as well as the effect of reactions and of the network itself on the bolus shape. Even though this study is guided by applications in liver micro-architecture, the proposed methodology is general and can readily be applied to other capillary network geometries, hence to other organs or to bioengineered network designs.

Predicting the risk of post-hepatectomy portal hypertension using a digital twin: A clinical proof of concept.

BACKGROUND & AIMS: Despite improvements in medical and surgical techniques, post-hepatectomy liver failure (PHLF) remains the leading cause of postoperative death. High postoperative portal vein pressure (PPV) and portocaval gradient (PCG), which cannot be predicted by current tools, are the most important determinants of PHLF. Therefore, we aimed to evaluate a digital twin to predict the risk of postoperative portal hypertension (PHT). METHODS: We prospectively included 47 patients undergoing major hepatectomy. A mathematical (0D) model of the entire blood circulation was assessed and automatically calibrated from patient characteristics. Hepatic flows were obtained from preoperative flow MRI (n = 9), intraoperative flowmetry (n = 16), or estimated from cardiac output (n = 47). Resection was then simulated in these 3 groups and the computed PPV and PCG were compared to intraoperative data. RESULTS: Simulated post-hepatectomy pressures did not differ between the 3 groups, comparing well with collected data (no significant differences). In the entire cohort, the correlation between measured and simulated PPV values was good (r = 0.66, no adjustment to intraoperative events) or excellent (r = 0.75) after adjustment, as well as for PCG (respectively r = 0.59 and r = 0.80). The difference between simulated and measured post-hepatectomy PCG was ≤3 mmHg in 96% of cases. Four patients suffered from lethal PHLF for whom the model satisfactorily predicted their postoperative pressures. CONCLUSIONS: We demonstrated that a 0D model could correctly anticipate postoperative PHT, even using estimated hepatic flow rates as input data. If this major conceptual step is confirmed, this algorithm could change our practice toward more tailor-made procedures, while ensuring satisfactory outcomes. LAY SUMMARY: Post-hepatectomy portal hypertension is a major cause of liver failure and death, but no tool is available to accurately anticipate this potentially lethal complication for a given patient. Herein, we propose using a mathematical model to predict the portocaval gradient at the end of liver resection. We tested this model on a cohort of 47 patients undergoing major hepatectomy and demonstrated that it could modify current surgical decision-making algorithms.

Verification of the coupled-momentum method with Womersley's Deformable Wall analytical solution.

In this paper, we perform a verification study of the Coupled-Momentum Method (CMM), a 3D fluid-structure interaction (FSI) model which uses a thin linear elastic membrane and linear kinematics to describe the mechanical behavior of the vessel wall. The verification of this model is done using Womersley's deformable wall analytical solution for pulsatile flow in a semi-infinite cylindrical vessel. This solution is, under certain premises, the analytical solution of the CMM and can thus be used for model verification. For the numerical solution, we employ an impedance boundary condition to define a reflection-free outflow boundary condition and thus mimic the physics of the analytical solution, which is defined on a semi-infinite domain. We first provide a rigorous derivation of Womersley's deformable wall theory via scale analysis. We then illustrate different characteristics of the analytical solution such as space-time wave periodicity and attenuation. Finally, we present the verification tests comparing the CMM with Womersley's theory.

Cohort-based multiscale analysis of hemodynamic-driven growth and remodeling of the embryonic pharyngeal arch arteries.

Growth and remodeling of the primitive pharyngeal arch artery (PAA) network into the extracardiac great vessels is poorly understood but a major source of clinically serious malformations. Undisrupted blood flow is required for normal PAA development, yet specific relationships between hemodynamics and remodeling remain largely unknown. Meeting this challenge is hindered by the common reductionist analysis of morphology to single idealized models, where in fact structural morphology varies substantially. Quantitative technical tools that allow tracking of morphological and hemodynamic changes in a population-based setting are essential to advancing our understanding of morphogenesis. Here, we have developed a methodological pipeline from high-resolution nano-computed tomography imaging and live-imaging flow measurements to multiscale pulsatile computational models. We combine experimental-based computational models of multiple PAAs to quantify hemodynamic forces in the rapidly morphing Hamburger Hamilton (HH) stage HH18, HH24 and HH26 embryos. We identify local morphological variation along the PAAs and their association with specific hemodynamic changes. Population-level mechano-morphogenic variability analysis is a powerful strategy for identifying stage-specific regions of well and poorly tolerated morphological and/or hemodynamic variation that may protect or initiate cardiovascular malformations.

A Lumped Parameter Model to Study Atrioventricular Valve Regurgitation in Stage 1 and Changes Across Stage 2 Surgery in Single Ventricle Patients.

GOAL: This manuscript evaluates atrioventric-ular valve regurgitation (AVVR) in babies born with an already very challenging heart condition, i.e., with single ventricle physiology. Although the second surgery that single ventricle patients undergo is thought to decrease AVVR, there is much controversy in the clinical literature about AVVR treatment. METHODS: The effect of AVVR on Stage 1 haemodynamics and resulting acute changes from conversion to Stage 2 circulation in single ventricle patients are analyzed through lumped parameter models. Several degrees of AVVR severity are analyzed, for two types of valve regurgitation: incomplete leaflet closure and valve prolapse. RESULTS: The models show that increasing AVVR in Stage 1 induces the following effects: first, higher stroke volume and associated decrease in ventricular end-systolic volume; second, increase in atrial volumes with V-loop enlargement in pressure-volume curves; third, pulmonary venous hypertension. The Stage 2 surgery results in volume unloading of the ventricle, thereby, driving a decrease in AVVR. However, this effect is offset by an increase in ventricular pressures resulting in a net increase in regurgitation fraction (RF) of approximately 0.1 (for example, in severe AVVR, the preoperative RF increases from 60% to 70% postoperatively). Moreover, despite some improvements to sarcomere function early after Stage 2 surgery, it may deteriorate in cases of severe AVVR. CONCLUSION: In patients with moderate to severe AVVR, restoration of atrioventricular valve competence prior to, or at the time of, Stage 2 surgery would likely lead to improved haemodynamics and clinical outcome as the models suggest that uncorrected AVVR can worsen across Stage 2 surgery. This was found to be independent of the AVVR degree and mechanisms.

Closed-Loop Lumped Parameter Modeling of Hemodynamics During Cirrhogenesis in Rats.

OBJECTIVE: Cirrhosis is the common end stage of any given chronic liver disease, developing after persistent destruction and regeneration of parenchymal liver cells. The associated architectural distortion increases the intrahepatic vascular resistance, leading to portal hypertension and systemic circulatory disorders. This study investigates the impact of the changing vascular resistances on the hepatic and global circulation hemodynamics during cirrhogenesis. METHODS: Cirrhogenesis was revisited using the thioacetamide rat model (N = 20). Rats were sacrificed at weeks 0, 6, 12, and 18. For each time-point, three-dimensional vascular geometries were created by combining hepatic vascular corrosion casting with µCT imaging. Morphological quantification of the trees branching topology provided the input for a lobe-specific lumped parameter model of the liver that was coupled to a closed-loop model of the entire circulation of the rat. Hemodynamics was simulated in physiological and pathological circumstances. RESULTS: The simulations showed the effect of the liver vascular resistances (driven by the hepatic venous resistance increase) on liver hemodynamics with portal hypertension observed after 12 weeks. The closed-loop model was further adapted to account for systemic circulatory compensation mechanisms and disorders frequently observed in cirrhosis and simulated their impact on the hepatic, systemic, and pulmonary hemodynamics. CONCLUSION: The simulations explain how vascular changes due to cirrhosis severely disrupt both hepatic and global hemodynamics. SIGNIFICANCE: This study is a priori the first to model the rat's entire blood circulation during cirrhogenesis. Since it is able to simulate cirrhosis main characteristics, the model may be translated to humans for the assessment of liver interventions.

Model and methods to assess hepatic function from indocyanine green fluorescence dynamical measurements of liver tissue.

The indocyanine green (ICG) clearance, presented as plasma disappearance rate is, presently, a reliable method to estimate the hepatic "function". However, this technique is not instantaneously available and thus cannot been used intra-operatively (during liver surgery). Near-infrared spectroscopy enables to assess hepatic ICG concentration over time in the liver tissue. This article proposes to extract more information from the liver intensity dynamics by interpreting it through a dedicated pharmacokinetics model. In order to account for the different exchanges between the liver tissues, the proposed model includes three compartments for the liver model (sinusoids, hepatocytes and bile canaliculi). The model output dependency to parameters is studied with sensitivity analysis and solving an inverse problem on synthetic data. The estimation of model parameters is then performed with in-vivo measurements in rabbits (El-Desoky et al. 1999). Parameters for different liver states are estimated, and their link with liver function is investigated. A non-linear (Michaelis-Menten type) excretion rate from the hepatocytes to the bile canaliculi was necessary to reproduce the measurements for different liver conditions. In case of bile duct ligation, the model suggests that this rate is reduced, and that the ICG is stored in the hepatocytes. Moreover, the level of ICG remains high in the blood following the ligation of the bile duct. The percentage of retention of indocyanine green in blood, which is a common test for hepatic function estimation, is also investigated with the model. The impact of bile duct ligation and reduced liver inflow on the percentage of ICG retention in blood is studied. The estimation of the pharmacokinetics model parameters may lead to an evaluation of different liver functions.

Computational simulation of postoperative pulmonary flow distribution in Alagille patients with peripheral pulmonary artery stenosis.

BACKGROUND: Up to 90% of individuals with Alagille syndrome have congenital heart diseases. Peripheral pulmonary artery stenosis (PPS), resulting in right ventricular hypertension and pulmonary flow disparity, is one of the most common abnormalities, yet the hemodynamic effects are ill-defined, and optimal patient management and treatment strategies are not well established. The purpose of this pilot study is to use recently refined computational simulation in the setting of multiple surgical strategies, to examine the influence of pulmonary artery reconstruction on hemodynamics in this population. MATERIALS AND METHODS: Based on computed tomography angiography and cardiac catheterization data, preoperative pulmonary artery models were constructed for 4 patients with Alagille syndrome with PPS (all male, age range: 0.6-2.9 years), and flow simulations with deformable walls were performed. Surgeon directed virtual surgery, mimicking the surgical procedure, was then performed to derive postoperative models. Postoperative simulation-derived hemodynamics and blood flow distribution were then compared with the clinical results. RESULTS: Simulations confirmed substantial resistance, resulting from preoperative severe ostial stenoses, and the use of newly developed adaptive outflow boundary conditions led to excellent agreement with in vivo measurements. Relief of PPS decreased pulmonary artery pressures and improved pulmonary flow distribution both in vivo and in silico with good correlation. CONCLUSIONS: Using adaptive outflow boundary conditions, computational simulations can estimate postoperative overall pulmonary flow distribution in patients with Alagille syndrome after pulmonary artery reconstruction. Obstruction relief along with pulmonary artery vasodilation determines postoperative pulmonary flow distribution and newer methods can incorporate these physiologic changes. Evolving blood flow simulations may be useful in surgical or transcatheter planning and in understanding the complex interplay among various obstructions in patients with peripheral pulmonary stenosis.

Tumor Cell Load and Heterogeneity Estimation From Diffusion-Weighted MRI Calibrated With Histological Data: an Example From Lung Cancer.

Diffusion-weighted magnetic resonance imaging (DWI) is a key non-invasive imaging technique for cancer diagnosis and tumor treatment assessment, reflecting Brownian movement of water molecules in tissues. Since densely packed cells restrict molecule mobility, tumor tissues produce usually higher signal (a.k.a. less attenuated signal) on isotropic maps compared with normal tissues. However, no general quantitative relation between DWI data and the cell density has been established. In order to link low-resolution clinical cross-sectional data with high-resolution histological information, we developed an image processing and analysis chain, which was used to study the correlation between the diffusion coefficient (D value) estimated from DWI and tumor cellularity from serial histological slides of a resected non-small cell lung cancer tumor. Color deconvolution followed by cell nuclei segmentation was performed on digitized histological images to determine local and cell-type specific 2d (two-dimensional) densities. From these, the 3d cell density was inferred by a model-based sampling technique, which is necessary for the calculation of local and global 3d tumor cell count. Next, DWI sequence information was overlaid with high-resolution CT data and the resected histology using prominent anatomical hallmarks for co-registration of histology tissue blocks and non-invasive imaging modalities' data. The integration of cell numbers information and DWI data derived from different tumor areas revealed a clear negative correlation between cell density and D value. Importantly, spatial tumor cell density can be calculated based on DWI data. In summary, our results demonstrate that tumor cell count and heterogeneity can be predicted from DWI data, which may open new opportunities for personalized diagnosis and therapy optimization.

Modulating Portal Hemodynamics With Vascular Ring Allows Efficient Regeneration After Partial Hepatectomy in a Porcine Model.

OBJECTIVE: To investigate safety and efficacy of temporary portal hemodynamics modulation with a novel percutaneously adjustable vascular ring (MID-AVR) onto a porcine model of 75% hepatectomy. BACKGROUND: Postoperative liver failure is a leading cause of mortality after major hepatectomy. Portal flow modulation is an increasingly accepted concept to prevent postoperative liver failure. Nonetheless, the current strategies have shortcomings. METHODS: Resection was performed under hemodynamic monitoring in 17 large, white pigs allocated into 2 groups. Eight pigs had ring around the portal vein for 3 days with the aim of reducing changes in hemodynamics due to hepatectomy. Analysis of hemodynamics, laboratory, and histopathological parameters was performed. RESULTS: Percutaneous inflation, deflation, and removal of the MID-AVR were safe. Two (25%) pigs in the MID-AVR group and 4 (45%) controls died before day 3 (P = NS). A moderate increase of portal flow rate per liver mass after resection was associated with better survival (P = 0.017). The portocaval pressure gradient was lower after hepatectomy in the MID-AVR group (P = 0.001). Postoperative serum bilirubin levels were lower in the MID-AVR group (P = 0.007 at day 5). In the MID-AVR group, the Ki67 index was significantly higher on day 3 (P = 0.043) and the architectural derangement was lower (P < 0.05). Morphometric quantification of the bile canaliculi revealed a significantly lower number of intersection branches (P < 0.05) and intersection nodes (P < 0.001) on day 7 compared with the preoperative specimen, in the control group. These differences were not found in the ring group. CONCLUSIONS: MID-AVR is safe for portal hemodynamics modulation. It might improve liver regeneration by protecting liver microarchitecture.

Inverse problems in reduced order models of cardiovascular haemodynamics: aspects of data assimilation and heart rate variability.

Inverse problems in cardiovascular modelling have become increasingly important to assess each patient individually. These problems entail estimation of patient-specific model parameters from uncertain measurements acquired in the clinic. In recent years, the method of data assimilation, especially the unscented Kalman filter, has gained popularity to address computational efficiency and uncertainty consideration in such problems. This work highlights and presents solutions to several challenges of this method pertinent to models of cardiovascular haemodynamics. These include methods to (i) avoid ill-conditioning of the covariance matrix, (ii) handle a variety of measurement types, (iii) include a variety of prior knowledge in the method, and (iv) incorporate measurements acquired at different heart rates, a common situation in the clinic where the patient state differs according to the clinical situation. Results are presented for two patient-specific cases of congenital heart disease. To illustrate and validate data assimilation with measurements at different heart rates, the results are presented on a synthetic dataset and on a patient-specific case with heart valve regurgitation. It is shown that the new method significantly improves the agreement between model predictions and measurements. The developed methods can be readily applied to other pathophysiologies and extended to dynamical systems which exhibit different responses under different sets of known parameters or different sets of inputs (such as forcing/excitation frequencies).

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.

Aerosol transport throughout inspiration and expiration in the pulmonary airways.

Little is known about transport throughout the respiration cycle in the conducting airways. It is challenging to appropriately describe the time-dependent number of particles entering back into the model during exhalation. Modeling the entire lung is not feasible; therefore, multidomain methods must be used. Here, we present a new framework that is designed to simulate particles throughout the respiration cycle, incorporating realistic airway geometry and respiration. This framework is applied for a healthy rat lung exposed to  âˆ¼ 1µm diameter particles, chosen to facilitate parameterization and validation. The flow field is calculated in the conducting airways (3D domain) by solving the incompressible Navier-Stokes equations with experimentally derived boundary conditions. Particles are tracked throughout inspiration by solving a modified Maxey-Riley equation. Next, we pass the time-dependent particle concentrations exiting the 3D model to the 1D volume conservation and advection-diffusion models (1D domain). Once the 1D models are solved, we prescribe the time-dependent number of particles entering back into the 3D airways to again solve for 3D transport. The coupled simulations highlight that about twice as many particles deposit during inhalation compared to exhalation for the entire lung. In contrast to inhalation, where most particles deposit at the bifurcation zones, particles deposit relatively uniformly on the gravitationally dependent side of the 3D airways during exhalation. Strong agreement to previously collected regional experimental data is shown, as the 1D models account for lobe-dependent morphology. This framework may be applied to investigate dosimetry in other species and pathological lungs.

Adaptive outflow boundary conditions improve post-operative predictions after repair of peripheral pulmonary artery stenosis.

Peripheral pulmonary artery stenosis (PPS) is a congenital abnormality resulting in pulmonary blood flow disparity and right ventricular hypertension. Despite recent advance in catheter-based interventions, surgical reconstruction is still preferred to treat complex PPS. However optimal surgical strategies remain unclear. It would be of great benefit to be able to predict post-operative hemodynamics to assist with surgical planning toward optimizing outcomes. While image-based computational fluid dynamics has been used in cardiovascular surgical planning, most studies have focused on the impact of local geometric changes on hemodynamic performance. Previous experimental studies suggest morphological changes in the pulmonary arteries not only alter local hemodynamics but also lead to distal pulmonary adaptation. In this proof of concept study, a constant shear stress hypothesis and structured pulmonary trees are used to derive adaptive outflow boundary conditions for post-operative simulations. Patient-specific simulations showed the adaptive outflow boundary conditions by the constant shear stress model to provide better predictions of pulmonary flow distribution than the conventional strategy of maintaining outflow boundary conditions. On average, the relative difference, when compared to the gold standard clinical test, in blood flow distribution to the right lung is reduced from 20 to 4 %. This suggests adaptive outflow boundary conditions should be incorporated into post-operative modeling in patients with complex PPS.