Doppler derived quantitative flow estimate in coronary artery bypass graft: a computational multiscale model for the evaluation of the current clinical procedure

Click method preserves but EDC method compromises the therapeutic activities of the peptide-activated hydrogels for critical ischemic vessel regeneration

Peptide-functionalized hydrogel is one of commonly used biomaterials to introduce hydrogel-induced vessel regeneration. Despite many reports about the discoveries of high-active peptides (or ligands) for regeneration, the study on the conjugating methods for the hydrogel functionalization with peptides is limited. Here, we compared the vasculogenic efficacy of the peptide-functionalized hydrogels prepared by two commonly used conjugating methods, 1-ethyl-3-(3-dimethylamino propyl) carbodiimide (EDC) and Click methods, through cell models, organ-on-chips models, animal models, and RNA sequencing analysis. Two vascular-related cell types, the human umbilical vein endothelial cells (HUVECs) and the adipose-derived stem cells (ADSCs), have been cultured on the hydrogel surfaces prepared by EDC/Click methods. It showed that the hydrogels prepared by Click method supported the higher vasculogenic activities while the ones made by EDC method compromised the peptide activities on hydrogels. The vasculogenesis assays further revealed that hydrogels prepared by Click method promoted a better vascular network formation. In a critical ischemic hindlimb model, only the peptide-functionalized hydrogels prepared by Click method successfully salvaged the ischemic limb, significantly improved blood perfusion, and enhanced the functional recoveries (through gait analysis and animal behavior studies). RNA sequencing studies revealed that the hydrogels prepared by Click method significantly promoted the PI3K-AKT pathway activation compared to the hydrogels prepared by EDC method. All the results suggested that EDC method compromised the functions of the peptides, while Click method preserved the vascular regenerating capacities of the peptides on the hydrogels, illustrating the importance of the conjugating method during the preparation of the peptide-functionalized hydrogels.

Modelling coronary flows: impact of differently measured inflow boundary conditions on vessel-specific computational hemodynamic profiles

BACKGROUND AND OBJECTIVES

The translation of hemodynamic quantities based on wall shear stress (WSS) or intravascular helical flow into clinical biomarkers of coronary atherosclerotic disease is still hampered by the assumptions/idealizations required by the computational fluid dynamics (CFD) simulations of the coronary hemodynamics. In the resulting budget of uncertainty, inflow boundary conditions (BCs) play a primary role. Accordingly, in this study we investigated the impact of the approach adopted for in vivo coronary artery blood flow rate assessment on personalized CFD simulations where blood flow rate is used as inflow BC.

METHODS

CFD simulations were carried out on coronary angiograms by applying personalized inflow BCs derived from four different techniques assessing in vivo surrogates of flow rate: continuous thermodilution, intravascular Doppler, frame count-based 3D contrast velocity, and diameter-based scaling law. The impact of inflow BCs on coronary hemodynamics was evaluated in terms of WSS- and helicity-based quantities.

RESULTS

As main findings, we report that: (i) coronary flow rate values may differ based on the applied flow derivation technique, as continuous thermodilution provided higher flow rate values than intravascular Doppler and diameter-based scaling law (p = 0.0014 and p = 0.0023, respectively); (ii) such intrasubject differences in flow rate values lead to different surface-averaged values of WSS magnitude and helical blood flow intensity (p<0.0020); (iii) luminal surface areas exposed to low WSS and helical flow topological features showed robustness to the flow rate values.

CONCLUSIONS

Although the absence of a clinically applicable gold standard approach prevents a general recommendation for one coronary blood flow rate derivation technique, our findings indicate that the inflow BC may impact computational hemodynamic results, suggesting that a standardization would be desirable to provide comparable results among personalized CFD simulations of the coronary hemodynamics.

Computational Fluid Dynamics Modeling of the Human Pulmonary Arteries with Experimental Validation

Pulmonary hypertension (PH) is a chronic progressive disease characterized by elevated pulmonary arterial pressure, caused by an increase in pulmonary arterial impedance. Computational fluid dynamics (CFD) can be used to identify metrics representative of the stage of PH disease. However, experimental validation of CFD models is often not pursued due to the geometric complexity of the model or uncertainties in the reproduction of the required flow conditions. The goal of this work is to validate experimentally a CFD model of a pulmonary artery phantom using a particle image velocimetry (PIV) technique. Rapid prototyping was used for the construction of the patient-specific pulmonary geometry, derived from chest computed tomography angiography images. CFD simulations were performed with the pulmonary model with a Reynolds number matching those of the experiments. Flow rates, the velocity field, and shear stress distributions obtained with the CFD simulations were compared to their counterparts from the PIV flow visualization experiments. Computationally predicted flow rates were within 1% of the experimental measurements for three of the four branches of the CFD model. The mean velocities in four transversal planes of study were within 5.9 to 13.1% of the experimental mean velocities. Shear stresses were qualitatively similar between the two methods with some discrepancies in the regions of high velocity gradients. The fluid flow differences between the CFD model and the PIV phantom are attributed to experimental inaccuracies and the relative compliance of the phantom. This comparative analysis yielded valuable information on the accuracy of CFD predicted hemodynamics in pulmonary circulation models.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

A simple framework to generate 3D patient-specific model of coronary artery bifurcation from single-plane angiographic images

Although computer-based simulations, such as structural finite element analysis, have proven their usefulness to support procedural planning of coronary stenting, the link between the clinical practice and these engineering techniques is still limited to research test-cases. A key point to further promote such an interaction is to generate in a fast and effective manner the computational grids from the medical images. Hence, the present study proposes a simple framework to generate 3D meshes of coronary bifurcations from a pair of planar angiographic images obtained by X-ray angiography, which is the gold standard technique for the diagnosis of coronary artery stenosis.

Considerations for numerical modeling of the pulmonary circulation--a review with a focus on pulmonary hypertension

Both in academic research and in clinical settings, virtual simulation of the cardiovascular system can be used to rapidly assess complex multivariable interactions between blood vessels, blood flow, and the heart. Moreover, metrics that can only be predicted with computational simulations (e.g., mechanical wall stress, oscillatory shear index, etc.) can be used to assess disease progression, for presurgical planning, and for interventional outcomes. Because the pulmonary vasculature is susceptible to a wide range of pathologies that directly impact and are affected by the hemodynamics (e.g., pulmonary hypertension), the ability to develop numerical models of pulmonary blood flow can be invaluable to the clinical scientist. Pulmonary hypertension is a devastating disease that can directly benefit from computational hemodynamics when used for diagnosis and basic research. In the present work, we provide a clinical overview of pulmonary hypertension with a focus on the hemodynamics, current treatments, and their limitations. Even with a rich history in computational modeling of the human circulation, hemodynamics in the pulmonary vasculature remains largely unexplored. Thus, we review the tasks involved in developing a computational model of pulmonary blood flow, namely vasculature reconstruction, meshing, and boundary conditions. We also address how inconsistencies between models can result in drastically different flow solutions and suggest avenues for future research opportunities. In its current state, the interpretation of this modeling technology can be subjective in a research environment and impractical for clinical practice. Therefore, considerations must be taken into account to make modeling reliable and reproducible in a laboratory setting and amenable to the vascular clinic. Finally, we discuss relevant existing models and how they have been used to gain insight into cardiopulmonary physiology and pathology.

Computational fluid dynamics simulations of contrast agent bolus dispersion in a coronary bifurcation: impact on MRI-based quantification of myocardial perfusion

Contrast-enhanced first-pass magnetic resonance imaging (MRI) in combination with a tracer kinetic model, for example, MMID4, can be used to determine myocardial blood flow (MBF) and myocardial perfusion reserve (MPR). Typically, the arterial input function (AIF) required for this methodology is estimated from the left ventricle (LV). Dispersion of the contrast agent bolus might occur between the LV and the myocardial tissue. Negligence of bolus dispersion could cause an error in MBF determination. The aim of this study was to investigate the influence of bolus dispersion in a simplified coronary bifurcation geometry including one healthy and one stenotic branch on the quantification of MBF and MPR. Computational fluid dynamics (CFD) simulations were combined with MMID4. Different inlet boundary conditions describing pulsatile and constant flows for rest and hyperemia and differing outflow conditions have been investigated. In the bifurcation region, the increase of the dispersion was smaller than inside the straight vessels. A systematic underestimation of MBF values up to -16.1% for pulsatile flow and an overestimation of MPR up to 7.5% were found. It was shown that, under the conditions considered in this study, bolus dispersion can significantly influence the results of quantitative myocardial MR-perfusion measurements.

Synthetic dataset generation for the analysis and the evaluation of image-based hemodynamics of the human aorta

Here, we consider the issue of generating a suitable controlled environment for the evaluation of phase contrast (PC) MRI measurements. The computational framework, tailored to build synthetic datasets, is based on a two-step approach, i.e., define and implement (1) an accurate CFD model and (2) an image generator able to mime the overall outcomes of a PC MRI acquisition starting from datasets retrieved by the computational model. About 20 different datasets were built by changing relevant image parameters (pixel size, slice thickness, time frames per cardiac cycle). Focusing our attention on the thoracic aorta, synthetic images were processed in order to: (1) verify to which extent the fluid dynamics into the aortic arch is influenced by the image parameters; (2) establish the effect of spatial and temporal interpolation. Our study demonstrates that the integral scale of the aortic bulk flow could be described satisfactorily even when using images which are nowadays acquirable with MRI scanners. However, attention must be paid to near-wall velocities that can be affected by large inaccuracy. In detail, in bulk flow regions error values are well bounded (below 5% for most of the analyzed resolutions), while errors greater than 100% are systematically present at the vessel's wall. Moreover, also the data interpolation process can be responsible for large inaccuracies in new data generation, due to the inherent complexity of the flow field in some connected regions.

Full-hexahedral structured meshing for image-based computational vascular modeling

Image-based computational modeling offers a virtual access to spatially and temporally high resolution flow and structural mechanical data in vivo. Due to inter-subject morphological variability, mesh generation represents a critical step in modeling the patient-specific geometry and is usually performed using unstructured tetrahedral meshing algorithms. Although hexahedral structured meshes are known to provide higher accuracy and reduce the computational costs both for Finite Element Analysis and Computational Fluid Dynamics, their application in computational cardiovascular studies is challenging due to the complex 3D and branching topology of vascular territories. In this study, we propose a robust procedure for structured mesh generation, tailoring the mesh structure to the subject-specific vessel topology. The proposed methodology is based on centerline-based synthetic descriptors (i.e. centerlines, radii and centerlines' normals) which are used to solve the meshing problem following a bottom-up approach. First, topologically equivalent block-structures are placed inside and outside the lumen domain. Then, a projection operation is performed, returning a parametric volume mesh which fits the original triangulated model with sub-micrometric accuracy. Additionally, a three-layered arterial wall (resembling the intima, media and adventitia) is artificially generated, with the possibility of setting variable thickness (e.g. proximal-to-distal tapering) and material anisotropy (e.g. position-dependent collagen-fibers' orientation). This new meshing procedure, implemented using open-source software packages only, is demonstrated on two challenging human cases, being an aortic arch and an abdominal aortic aneurysm. High-quality meshes are generated in both cases, according to shape-quality metrics. By increasing the computation accuracy, the developed meshing tool has the potential to further add "confidence" to the use of computational methods in vascular applications.

Theoretical models for coronary vascular biomechanics: progress & challenges

A key aim of the cardiac Physiome Project is to develop theoretical models to simulate the functional behaviour of the heart under physiological and pathophysiological conditions. Heart function is critically dependent on the delivery of an adequate blood supply to the myocardium via the coronary vasculature. Key to this critical function of the coronary vasculature is system dynamics that emerge via the interactions of the numerous constituent components at a range of spatial and temporal scales. Here, we focus on several components for which theoretical approaches can be applied, including vascular structure and mechanics, blood flow and mass transport, flow regulation, angiogenesis and vascular remodelling, and vascular cellular mechanics. For each component, we summarise the current state of the art in model development, and discuss areas requiring further research. We highlight the major challenges associated with integrating the component models to develop a computational tool that can ultimately be used to simulate the responses of the coronary vascular system to changing demands and to diseases and therapies.

Womersley number-based estimates of blood flow rate in Doppler analysis: in vivo validation by means of phase-contrast MRI

A common clinical practice during single-point Doppler analysis is to measure the centerline maximum velocity and to recover the time-averaged flow rate by exploiting an assumption on the shape of velocity profile (a priori formula), either a parabolic or a flat one. In a previous study, we proposed a new formula valid for the peak instant linking the maximum velocity and the flow rate by including a well-established dimensionless fluid-dynamics parameter (the Womersley number), in order to account for the hemodynamics conditions (Womersley number-based formula). Several in silico tests confirmed the reliability of the new formula. Nevertheless, an in vivo confirmation is missing limiting the clinical applicability of the formula. An experimental in vivo protocol using cine phase-contrast MRI (2-D PCMRI) technique has been designed and applied to ten healthy young volunteers in three different arterial districts: the abdominal aorta, the common carotid artery, and the brachial artery. Each PCMRI dataset has been used twice: 1) to compute the value of the blood flow rate used as a gold standard and 2) to estimate the flow rate by measuring directly the maximum velocity and the diameter (i.e., emulating the intravascular Doppler data acquisition) and by applying to these data the a priori and the Womersley number-based formulae. All the in vivo results have confirmed that the Womersley number-based formula provides better estimates of the flow rate at the peak instant with respect to the a priori formula. More precisely, mean performances of the Womersley number-based formula are about three times better than the a priori results in the abdominal aorta, five times better in the common carotid artery, and two times better in the brachial artery.

Outflow Conditions for Image-Based Hemodynamic Models of the Carotid Bifurcation: Implications for Indicators of Abnormal Flow

Computational fluid dynamics (CFD) models have become very effective tools for predicting the flow field within the carotid bifurcation, and for understanding the relationship between local hemodynamics, and the initiation and progression of vascular wall pathologies. As prescribing proper boundary conditions can affect the solutions of the equations governing blood flow, in this study, we investigated the influence to assumptions regarding the outflow boundary conditions in an image-based CFD model of human carotid bifurcation. Four simulations were conducted with identical geometry, inlet flow rate, and fluid parameters. In the first case, a physiological time-varying flow rate partition at branches along the cardiac cycle was obtained by coupling the 3D model of the carotid bifurcation at outlets with a lumped-parameter model of the downstream vascular network. Results from the coupled model were compared with those obtained by imposing three fixed flow rate divisions (50/50, 60/40, and 70/30) between the two branches of the isolated 3D model of the carotid bifurcation. Three hemodynamic wall parameters were considered as indicators of vascular wall dysfunction. Our findings underscore that the overall effect of the assumptions done in order to simulate blood flow within the carotid bifurcation is mainly in the hot-spot modulation of the hemodynamic descriptors of atherosusceptible areas, rather than in their distribution. In particular, the more physiological, time-varying flow rate division deriving from the coupled simulation has the effect of damping wall shear stress (WSS) oscillations (differences among the coupled and the three fixed flow partition models are up to 37.3% for the oscillating shear index). In conclusion, we recommend to adopt more realistic constraints, for example, by coupling models at different scales, as in this study, when the objective is the outcome prediction of alternate therapeutic interventions for individual patients, or to test hypotheses related to the role of local fluid dynamics and other biomechanical factors in vascular diseases.

Womersley number-based estimation of flow rate with Doppler ultrasound: sensitivity analysis and first clinical application

In this paper we continue in investigating the approach we have proposed in a paper recently published, for a reliable estimate of (peak systolic) blood flow rate from velocity Doppler measurements. Basic features of this approach together with some in silico test cases were discussed in that work. Here, we provide more insights of this approach by performing a sensitivity analysis of the formulas relating blood flow rate to velocity. In particular we analyze how our estimates are affected by perturbation or errors in measurements in comparison with a standard method for catheter based estimates based on the assumption of a parabolic velocity profile. A first glance to in vivo clinical applications is given as well.