The influence of boundary conditions on wall shear stress distribution in patients specific coronary trees

Evaluation of enhanced external counterpulsation for diabetic foot based on a patient-specific 0D-1D cardiovascular system model

BACKGROUND AND OBJECTIVE

Diabetic foot (DF) complications often lead to severe vascular issues. This study investigated the effectiveness of enhanced external counterpulsation (EECP) and its derived innovative compression strategies in addressing poor perfusion in DF. Although developing non-invasive and efficient treatment methods for DF is critical, the hemodynamic alterations during EECP remain underexplored despite promising outcomes in microcirculation. This research sought to address this gap by developing a patient-specific 0D-1D model based on clinical ultrasound data to identify potentially superior compression strategies that could substantially enhance blood flow in patients with DF complications.

METHODS

Data were gathered from 10 patients with DF utilizing ultrasound for blood flow rate and computed tomography angiography (CTA) to identify lower limb conditions. Clinical measurements during standard EECP, with varying cuff pressures, facilitated the creation of a patient-specific 0D-1D model through a two-step parameter estimation process. The accuracy of this model was verified via comparison with the clinical measurements. Four compression strategies were proposed and rigorously evaluated using this model: EECP-Simp-I (removing hip cuffs), EECP-Simp-II (further removing the cuffs around the lower leg), EECP-Impr-I (removing all cuffs around the affected side), and EECP-Impr-II (building a loop circulation from the healthy side to the affected side).

RESULTS

The predicted results under the rest and standard EECP states were generally closely aligned with clinical measurements. The patient-specific 0D-1D model demonstrated that EECP-Simp-I and EECP-Impr-I contributed similar enhancement to perfusion in the dorsal artery (DA) and were comparable to standard EECP, while EECP-Simp-II had the least effect and EECP-Impr-II displayed the most significant enhancement. Pressure at the aortic root (AO) remained consistent across strategies.

CONCLUSIONS

EECP-Simp-I is recommended for patients with DF, emphasizing device simplification. However, EECP-Simp-II is discouraged as it significantly diminished blood perfusion in this study, except in cases of limb fragility. EECP-Impr-II showed superior enhancement of blood perfusion in DA to all other strategies but required a more complex EECP device. Despite increased AO pressure in all the proposed compression strategies, safety could be guaranteed as the pressue remained within a safe range.

Influence of the pressure wire on the fractional flow reserve calculation: CFD analysis of an ideal vessel and clinical patients with stenosis

BACKGROUND AND OBJECTIVE

Fractional Flow Reserve (FFR) is generally considered the gold standard in hemodynamics to assess the impact of a stenosis on the blood flow. The standard procedure to measure involves the displacement of a pressure guide along the circulatory system until it is placed next to the lesion to be analyzed. The main objective of the present study is to analyze the influence of the pressure guide on the invasive FFR measurements and its implications in clinical practice.

METHODS

We studied the influence of pressure wires on the measurement of Fractional Flow Reserve (FFR) through a combination of Computational Fluid Dynamics (CFD) simulations using 45 clinical patient data with 58 lesions and ideal geometries. The analysis is conducted considering patients that were subjected to a computer tomography and also have direct measurements using a pressure guide. Influence of the stenosis severity, degree of occlusion and blood viscosity has also been studied.

RESULTS

The influence of pressure wires specifically affects severe stenosis with a lumen diameter reduction of 50 % or greater. This type of stenosis leads to reduced hyperemic flow and increased coronary pressure drop. Thus, we identified that the placement of wires during FFR measurements results in partial obstruction of the coronary artery lumen, leading to increased pressure drop and subsequent reduction in blood flow. The severity of low FFR values associated with severe stenosis may be prone to overestimation when compared to stenosis without severe narrowing. These results have practical implications, particularly in the interpretation of lesions falling within the "gray zone" (0,75-0,80).

CONCLUSIONS

The pressure wire's presence significantly alters the flow on severe lesions, which has an impact on the FFR calculation. In contrast, the impact of the pressure wire appears to be reduced when the FFR is larger than 0.8. The findings provide critical information for physicians, emphasizing the need for cautious interpretation of FFR values, particularly in severe stenosis. It also offers insights into improving the correlation between FFRct models and invasive measurements by incorporating the influence of pressure wires.

Systematic review and meta-analysis of Murray's law in the coronary arterial circulation

Murray's law has been viewed as a fundamental law of physiology. Relating blood flow ([Formula: see text]) to vessel diameter (D) ([Formula: see text]·∝·D3), it dictates minimum lumen area (MLA) targets for coronary bifurcation percutaneous coronary intervention (PCI). The cubic exponent (3.0), however, has long been disputed, with alternative theoretical derivations, arguing this should be closer to 2.33 (7/3). The aim of this meta-analysis was to quantify the optimum flow-diameter exponent in human and mammalian coronary arteries. We conducted a systematic review and meta-analysis of all articles quantifying an optimum flow-diameter exponent for mammalian coronary arteries within the Cochrane library, PubMed Medline, Scopus, and Embase databases on 20 March 2023. A random-effects meta-analysis was used to determine a pooled flow-diameter exponent. Risk of bias was assessed with the National Institutes of Health (NIH) quality assessment tool, funnel plots, and Egger regression. From a total of 4,772 articles, 18 were suitable for meta-analysis. Studies included data from 1,070 unique coronary trees, taken from 372 humans and 112 animals. The pooled flow diameter exponent across both epicardial and transmural arteries was 2.39 (95% confidence interval: 2.24-2.54; I2 = 99%). The pooled exponent of 2.39 showed very close agreement with the theoretical exponent of 2.33 (7/3) reported by Kassab and colleagues. This exponent may provide a more accurate description of coronary morphometric scaling in human and mammalian coronary arteries, as compared with Murray's original law. This has important implications for the assessment, diagnosis, and interventional treatment of coronary artery disease.

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Experimental validation and clinical feasibility of 3D reconstruction of coronary artery bifurcation stents using intravascular ultrasound

The structural morphology of coronary stents and the local hemodynamic environment following stent deployment in coronary arteries are crucial determinants of procedural success and subsequent clinical outcomes. High-resolution intracoronary imaging has the potential to facilitate geometrically accurate three-dimensional (3D) reconstruction of coronary stents. This work presents an innovative algorithm for the 3D reconstruction of coronary artery stents, leveraging intravascular ultrasound (IVUS) and angiography. The accuracy and reproducibility of our method were tested in stented patient-specific silicone models, with micro-computed tomography serving as a reference standard. We also evaluated the clinical feasibility and ability to perform computational fluid dynamics (CFD) studies in a clinically stented coronary bifurcation. Our experimental and clinical studies demonstrated that our proposed algorithm could reproduce the complex 3D stent configuration with a high degree of precision and reproducibility. Moreover, the algorithm was proved clinically feasible in cases with stents deployed in a diseased coronary artery bifurcation, enabling CFD studies to assess the hemodynamic environment. In combination with patient-specific CFD studies, our method can be applied to stenting optimization, training in stenting techniques, and advancements in stent research and development.

Drug Delivery Angle for Various Atherosclerosis and Aneurysm Percentages of the Carotid Artery

Stroke is the second cause of mortality among adult males and the first cause of death in adult females all around the world. It is also recognized as one of the most important causes of morbidity and dementia in adults. Stenosis or rupture of the only channels of the blood supply from the heart to the brain (carotid arteries) is among the main causes of stroke. In this regard, treatment of the lesions of carotid arteries, including atherosclerosis and aneurysm, could be a huge step in preventing stroke and improving brain performance. Targeted drug delivery by drug-carrying nanoparticles is the latest method for optimal delivery of drug to the damaged parts of the artery. In this study, a wide range of carotid artery lesions, including different percentages of atherosclerosis and aneurysm, were considered. After analyzing the dynamics of the fluid flow in different damaged regions and selecting the magnetic framework with proper ligand (Fe3O4@MOF) as the drug carrier, the size of the particles and their number per cycle were analyzed. Based on the results, the particle size of 100 nm and the use of 300 particles per injection at each cardiac cycle can result in maximum drug delivery to the target site. Then, the effect of the hospital bed angle on drug delivery was investigated. The results showed a unique optimal drug delivery angle for each extent of atherosclerosis or aneurysm. For example, in a 50% aneurysm, drug delivery at an angle of 30° is about 387% higher than that at an angle of 15°. Finally, simulation of real geometry indicated the effectiveness of simple geometry instead of real geometry for the simulation of carotid arteries, which can remarkably decrease the computational time and costs.

Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures

BACKGROUND

Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis.

METHODS

Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the timeaveraged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, which can be assessed by the topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage atheroma volume change at 1-year follow-up (Δplaque atheroma volume). Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography.

RESULTS

Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich content is exposed to high TSVI or low TAWSS.

CONCLUSIONS

Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

Blood Flow Energy Identifies Coronary Lesions Culprit of Future Myocardial Infarction

The present study establishes a link between blood flow energy transformations in coronary atherosclerotic lesions and clinical outcomes. The predictive capacity for future myocardial infarction (MI) was compared with that of established quantitative coronary angiography (QCA)-derived predictors. Angiography-based computational fluid dynamics (CFD) simulations were performed on 80 human coronary lesions culprit of MI within 5 years and 108 non-culprit lesions for future MI. Blood flow energy transformations were assessed in the converging flow segment of the lesion as ratios of kinetic and rotational energy values (KER and RER, respectively) at the QCA-identified minimum lumen area and proximal lesion sections. The anatomical and functional lesion severity were evaluated with QCA to derive percentage area stenosis (%AS), vessel fractional flow reserve (vFFR), and translesional vFFR (ΔvFFR). Wall shear stress profiles were investigated in terms of topological shear variation index (TSVI). KER and RER predicted MI at 5 years (AUC = 0.73, 95% CI 0.65-0.80, and AUC = 0.76, 95% CI 0.70-0.83, respectively; p < 0.0001 for both). The predictive capacity for future MI of KER and RER was significantly stronger than vFFR (p = 0.0391 and p = 0.0045, respectively). RER predictive capacity was significantly stronger than %AS and ΔvFFR (p = 0.0041 and p = 0.0059, respectively). The predictive capacity for future MI of KER and RER did not differ significantly from TSVI. Blood flow kinetic and rotational energy transformations were significant predictors for MI at 5 years (p < 0.0001). The findings of this study support the hypothesis of a biomechanical contribution to the process of plaque destabilization/rupture leading to MI.

The influence of flow distribution strategy for the quantification of pressure- and wall shear stress-derived parameters in the coronary artery: A CTA-based computational fluid dynamics analysis

For image-based computational fluid dynamics (CFD) analysis to characterize the local coronary hemodynamic environment, the accuracy depends on the flow rate which is in turn associated with outlet branches' morphology. A good flow distribution strategy is important to mitigate the effect when certain branches cannot be considered. In this study, stenotic coronary arteries from 13 patients were used to analyze the effect of missing branches and different flow distribution strategies. Pressure- and wall shear stress (WSS)-derived parameters around the stenotic region (ROI) were compared, including fractional flow reserve (CT-FFR), instantaneous wave-free ratio (CT-iFR), resting distal to aortic coronary pressure (CT-Pd/Pa), time-averaged WSS, oscillatory shear index (OSI) and relative residence time (RRT). Three flow distribution strategies were the Huo-Kassab model at distal outlets (Type I), flow distribution based on outlet resistances (Type II), and a developed algorithm distributing flow at each bifurcation until the final outlets (Type III). Results showed that Type III strategy for models with truncated branch(es) had a good agreement in both pressure- and WSS-related results (interquatile range less than 0.12% and 4.02%, respectively) with the baseline model around the ROI. The relative difference of pressure- and WSS-related results were correlated with the flow differences in the ROI to the baseline mode. Type III strategy had the best performance in maintaining the flow in intermediate branches. It is recommended for CFD analysis. Removal of branches distal to a stenosis can be undertaken with an improved performance and maintained accuracy, while those proximal to the ROI should be kept.

Modelling blood flow in coronary arteries: Newtonian or shear-thinning non-Newtonian rheology?

BACKGROUND

The combination of medical imaging and computational hemodynamics is a promising technology to diagnose/prognose coronary artery disease (CAD). However, the clinical translation of in silico hemodynamic models is still hampered by assumptions/idealizations that must be introduced in model-based strategies and that necessarily imply uncertainty. This study aims to provide a definite answer to the open question of how to properly model blood rheological properties in computational fluid dynamics (CFD) simulations of coronary hemodynamics.

METHODS

The geometry of the right coronary artery (RCA) of 144 hemodynamically stable patients with different stenosis degree were reconstructed from angiography. On them, unsteady-state CFD simulations were carried out. On each reconstructed RCA two different simulation strategies were applied to account for blood rheological properties, implementing (i) a Newtonian (N) and (ii) a shear-thinning non-Newtonian (non-N) rheological model. Their impact was evaluated in terms of wall shear stress (WSS magnitude, multidirectionality, topological skeleton) and helical flow (strength, topology) profiles. Additionally, luminal surface areas (SAs) exposed to shear disturbances were identified and the co-localization of paired N and non-N SAs was quantified in terms of similarity index (SI).

RESULTS

The comparison between paired N vs. shear-thinning non-N simulations revealed remarkably similar profiles of WSS-based and helicity-based quantities, independent of the adopted blood rheology model and of the degree of stenosis of the vessel. Statistically, for each paired N and non-N hemodynamic quantity emerged negligible bias from Bland-Altman plots, and strong positive linear correlation (r > 0.94 for almost all the WSS-based quantities, r > 0.99 for helicity-based quantities). Moreover, a remarkable co-localization of N vs. non-N luminal SAs exposed to disturbed shear clearly emerged (SI distribution 0.95 [0.93, 0.97]). Helical flow topology resulted to be unaffected by blood rheological properties.

CONCLUSIONS

This study, performed on 288 angio-based CFD simulations on 144 RCA models presenting with different degrees of stenosis, suggests that the assumptions on blood rheology have negligible impact both on WSS and helical flow profiles associated with CAD, thus definitively answering to the question "is Newtonian assumption for blood rheology adequate in coronary hemodynamics simulations?".

Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees

Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary models combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. Personalized flow distributions and model parameters are informed by clinical CT myocardial perfusion imaging and cardiac function using surrogate-based optimization. We reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.

Agreement of wall shear stress distribution between two core laboratories using three-dimensional quantitative coronary angiography

Wall shear stress (WSS) estimated in models reconstructed from intravascular imaging and 3-dimensional-quantitative coronary angiography (3D-QCA) data provides important prognostic information and enables identification of high-risk lesions. However, these analyses are time-consuming and require expertise, limiting WSS adoption in clinical practice. Recently, a novel software has been developed for real-time computation of time-averaged WSS (TAWSS) and multidirectional WSS distribution. This study aims to examine its inter-corelab reproducibility. Sixty lesions (20 coronary bifurcations) with a borderline negative fractional flow reserve were processed using the CAAS Workstation WSS prototype to estimate WSS and multi-directional WSS values. Analysis was performed by two corelabs and their estimations for the WSS in 3 mm segments across each reconstructed vessel was extracted and compared. In total 700 segments (256 located in bifurcated vessels) were included in the analysis. A high intra-class correlation was noted for all the 3D-QCA and TAWSS metrics between the estimations of the two corelabs irrespective of the presence (range: 0.90-0.92) or absence (range: 0.89-0.90) of a coronary bifurcation, while the ICC was good-moderate for the multidirectional WSS (range: 0.72-0.86). Lesion level analysis demonstrated a high agreement of the two corelabls for detecting lesions exposed to an unfavourable haemodynamic environment (WSS > 8.24 Pa, κ = 0.77) that had a high-risk morphology (area stenosis > 61.3%, κ = 0.71) and were prone to progress and cause events. The CAAS Workstation WSS enables reproducible 3D-QCA reconstruction and computation of WSS metrics. Further research is needed to explore its value in detecting high-risk lesions.

Wall shear stress-related plaque growth of lipid-rich plaques in human coronary arteries: an near-infrared spectroscopy and optical coherence tomography study

AIMS

Low wall shear stress (WSS) is acknowledged to play a role in plaque development through its influence on local endothelial function. Also, lipid-rich plaques (LRPs) are associated with endothelial dysfunction. However, little is known about the interplay between WSS and the presence of lipids with respect to plaque progression. Therefore, we aimed to study the differences in WSS-related plaque progression between LRPs, non-LRPs, or plaque-free regions in human coronary arteries.

METHODS AND RESULTS

In the present single-centre, prospective study, 40 patients who presented with an acute coronary syndrome successfully underwent near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) and optical coherence tomography (OCT) of at least one non-culprit vessel at baseline and completed a 1-year follow-up. WSS was computed applying computational fluid dynamics to a three-dimensional reconstruction of the coronary artery based on the fusion of the IVUS-segmented lumen with a CT-derived centreline, using invasive flow measurements as boundary conditions. For data analysis, each artery was divided into 1.5 mm/45° sectors. Plaque growth based on IVUS-derived percentage atheroma volume change was compared between LRPs, non-LRPs, and plaque-free wall segments, as assessed by both OCT and NIRS. Both NIRS- and OCT-detected lipid-rich sectors showed a significantly higher plaque progression than non-LRPs or plaque-free regions. Exposure to low WSS was associated with a higher plaque progression than exposure to mid or high WSS, even in the regions classified as a plaque-free wall. Furthermore, low WSS and the presence of lipids had a synergistic effect on plaque growth, resulting in the highest plaque progression in lipid-rich regions exposed to low shear stress.

CONCLUSION

This study demonstrates that NIRS- and OCT-detected lipid-rich regions exposed to low WSS are subject to enhanced plaque growth over a 1-year follow-up. The presence of lipids and low WSS proves to have a synergistic effect on plaque growth.

Simulation of LDL permeation into multilayer wall of a coronary bifurcation using WSS-dependent model: effects of hemorheology

Atherosclerosis, due to the permeation of low-density lipoprotein (LDL) particles into the arterial wall, is one of the most common and deadly diseases in today's world. Due to its importance, numerous studies have been conducted on the factors affecting this disease. In this study, using numerical simulation, the effects of Wall Shear Stress (WSS), non-Newtonian behavior of blood, different values of hematocrit and blood pressure on LDL permeation into the arterial wall layers are investigated in a 4-layer wall model of a coronary bifurcation. To obtain the velocity and concentration fields in the fluid domain, the Navier-Stokes, Brinkman, and mass transfer equations are numerically solved in the lumen and wall layers. Results show that it is important to consider the effects of WSS on transport properties of endothelium layer in bifurcations and this leads to completely different concentration profiles compared to the constant properties model. Our computations show that a giant accumulation of LDL in the intima layer of the outer wall of the left anterior descending artery, especially in low WSS regions, may lead to atherosclerosis. It is also, necessary to consider the non-Newtonian behavior of blood in bifurcations due to its direct effect on WSS. A pressure-induced increase in the half-width of leaky junctions may be responsible for the higher risk of atherosclerosis in hypertension. In addition, it is shown that the dominant mechanism in LDL permeation into the wall is convection, and also, hypertension increases the effect of mass transfer by convection mechanism more than the diffusion mechanism. Furthermore, our results are consistent with various clinical studies.

Transient wall shear stress estimation in coronary bifurcations using convolutional neural networks

BACKGROUND AND OBJECTIVE

Haemodynamic metrics, such as blood flow induced shear stresses at the inner vessel lumen, are associated with the development and progression of coronary artery disease. Understanding these metrics may therefore improve the assessment of an individual's coronary disease risk. However, the calculation of such luminal Wall Shear Stress (WSS) using traditional Computational Fluid Dynamics (CFD) methods is relatively slow and computationally expensive. As a result, CFD based haemodynamic computation is not suitable for integrated and large-scale use in clinical settings.

METHODS

In this work, deep learning techniques are proposed as an alternative method to CFD, whereby luminal WSS magnitude can be predicted in coronary bifurcations throughout the cardiac cycle based on the steady state solution (which takes <120 seconds to calculate including preprocessing), vessel geometry and additional global features. The deep learning model is trained on a dataset of 101 patient-specific and 2626 synthetic left main bifurcation models with 26 separate patient-specific cases used as the test set.

RESULTS

The model showed high fidelity predictions with <5% (normalised against mean WSS magnitude) deviation to CFD derived values as the gold-standard method, while being orders of magnitude faster with on average <2 minutes versus 3 hours computation for transient CFD.

CONCLUSIONS

This method therefore offers a new approach to substantially reduce the computational cost involved in, for example, large-scale population studies of coronary haemodynamic metrics, and may therefore open the pathway for future clinical integration.

Challenges in Modeling Hemodynamics in Cerebral Aneurysms Related to Arteriovenous Malformations

PURPOSE

The significantly higher incidence of aneurysms in patients with arteriovenous malformations (AVMs) suggests a strong hemodynamic relationship between these lesions. The presence of an AVM alters hemodynamics in proximal vessels by drastically changing the distal resistance, thus affecting intra-aneurysmal flow. This study discusses the challenges associated with patient-specific modeling of aneurysms in the presence of AVMs.

METHODS

We explore how the presence of a generic distal AVM affects upstream aneurysms by examining the relationship between distal resistance and aneurysmal wall shear stress using physiologically realistic estimates for the influence of the AVM on hemodynamics. Using image-based computational models of aneurysms and surrounding vasculature, aneurysmal wall-shear stress is calculated for a range of distal resistances corresponding to the presence of AVMs of various sizes and compared with a control case representing the absence of an AVM.

RESULTS

In the patient cases considered, the alteration in aneurysmal wall shear stress due to the presence of an AVM is considerable, as much as 19 times the base case wall shear stress. Furthermore, the relationship between aneurysmal wall shear stress and distal resistance is shown to be highly geometry-dependent and nonlinear. In most cases, the range of physiologically realistic possibilities for AVM-related distal resistance are so large that patient-specific flow measurements are necessary for meaningful predictions of wall shear stress.

CONCLUSIONS

The presented work offers insight on the impact of distal AVMs on aneurysmal wall shear stress using physiologically realistic computational models. Patient-specific modeling of hemodynamics in aneurysms and associated AVMs has great potential for understanding lesion pathogenesis, surgical planning, and assessing the effect of treatment of one lesion relative to another. However, we show that modeling approaches cannot usually meaningfully quantify the impact of AVMs if based solely on imaging data from CT and X-ray angiography, currently used in clinical practice. Based on recent studies, it appears that 4D flow MRI is one promising approach to obtaining meaningful patient-specific flow boundary conditions that improve modeling fidelity.

Helical Flow in Healthy and Diseased Patient-specific Coronary Bifurcations

Helical flow (HF) exists in healthy and diseased coronary bifurcations and was found to have a protective atherosclerotic vascular effect in other vessels. However, the role of HF in patient-specific human coronary arteries still needs further study, and is therefore the objective of this study in both healthy and diseased bifurcations. Computational studies were conducted on 16 patient-specific coronary bifurcations, including eight healthy and eight identical cases with idealized narrowing to represent disease. In general, higher HF intensity may have a favorable effect as it corelated to the reduction of the percentage vessel area exposed to adverse time averaged wall shear stress (TAWSS%) in both healthy and diseased models. The HF intensity and distribution of each model varies due to the complex shape of patient-specific models. The presence of disease appears to have an important impact on the downstream HF patterns and the TAWSS distributions. Clinical Relevance- By understanding the relationship between HF and hemodynamics, HF may be used as a predictor for the formation and progression of atherosclerotic plaque in coronary arteries instead of near-wall WSS measures, which can be determined with higher accuracy in vivo.

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.

Refining Our Understanding of the Flow Through Coronary Artery Branches; Revisiting Murray's Law in Human Epicardial Coronary Arteries

Background: Quantification of coronary blood flow is used to evaluate coronary artery disease, but our understanding of flow through branched systems is poor. Murray’s law defines coronary morphometric scaling, the relationship between flow (Q) and vessel diameter (D) and is the basis for minimum lumen area targets when intervening on bifurcation lesions. Murray’s original law (Q α D^P) dictates that the exponent (P) is 3.0, whilst constant blood velocity throughout the system would suggest an exponent of 2.0. In human coronary arteries, the value of Murray’s exponent remains unknown.

Aim: To establish the exponent in Murray’s power law relationship that best reproduces coronary blood flows (Q) and microvascular resistances (Rmicro) in a bifurcating coronary tree.

Methods and Results: We screened 48 cases, and were able to evaluate inlet Q and Rmicro in 27 branched coronary arteries, taken from 20 patients, using a novel computational fluid dynamics (CFD) model which reconstructs 3D coronary anatomy from angiography and uses pressure-wire measurements to compute Q and Rmicro distribution in the main- and side-branches. Outputs were validated against invasive measurements using a Rayflow™ catheter. A Murray’s power law exponent of 2.15 produced the strongest correlation and closest agreement with inlet Q (zero bias, r = 0.47, p = 0.006) and an exponent of 2.38 produced the strongest correlation and closest agreement with Rmicro (zero bias, r = 0.66, p = 0.0001).

Conclusions: The optimal power law exponents for Q and Rmicro were not 3.0, as dictated by Murray’s Law, but 2.15 and 2.38 respectively. These data will be useful in assessing patient-specific coronary physiology and tailoring revascularisation decisions.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.

Automating fractional flow reserve (FFR) calculation from CT scans: A rapid workflow using unsupervised learning and computational fluid dynamics

Fractional flow reserve (FFR) provides the functional relevance of coronary atheroma. The FFR-guided strategy has been shown to reduce unnecessary stenting, improve overall health outcome, and to be cost-saving. The non-invasive, coronary computerised tomography (CT) angiography-derived FFR (cFFR) is an emerging method in reducing invasive catheter based measurements. This computational fluid dynamics-based method is laborious as it requires expertise in multidisciplinary analysis of combining image analysis and computational mechanics. In this work, we present a rapid method, powered by unsupervised learning, to automatically calculate cFFR from CT scans without manual intervention.

Accuracy of vascular tortuosity measures using computational modelling

Severe coronary tortuosity has previously been linked to low shear stresses at the luminal surface, yet this relationship is not fully understood. Several previous studies considered different tortuosity metrics when exploring its impact of on the wall shear stress (WSS), which has likely contributed to the ambiguous findings in the literature. Here, we aim to analyze different tortuosity metrics to determine a benchmark for the highest correlating metric with low time-averaged WSS (TAWSS). Using Computed Tomography Coronary Angiogram (CTCA) data from 127 patients without coronary artery disease, we applied all previously used tortuosity metrics to the left main coronary artery bifurcation, and to its left anterior descending and left circumflex branches, before modelling their TAWSS using computational fluid dynamics (CFD). The tortuosity measures included tortuosity index, average absolute-curvature, root-mean-squared (RMS) curvature, and average squared-derivative-curvature. Each tortuosity measure was then correlated with the percentage of vessel area that showed a < 0.4 Pa TAWSS, a threshold associated with altered endothelial cell cytoarchitecture and potentially higher disease risk. Our results showed a stronger correlation between curvature-based versus non-curvature-based tortuosity measures and low TAWSS, with the average-absolute-curvature showing the highest coefficient of determination across all left main branches (p < 0.001), followed by the average-squared-derivative-curvature (p = 0.001), and RMS-curvature (p = 0.002). The tortuosity index, the most widely used measure in literature, showed no significant correlation to low TAWSS (p = 0.86). We thus recommend the use of average-absolute-curvature as a tortuosity measure for future studies.

Coronary Artery Stenting Affects Wall Shear Stress Topological Skeleton

Despite the important advancements in the stent technology for the treatment of diseased coronary arteries, major complications still affect the post-operative long-term outcome. The stent-induced flow disturbances, and especially the altered wall shear stress (WSS) profile at the strut level, play an important role in the pathophysiological mechanisms leading to stent thrombosis (ST) and in-stent restenosis (ISR). In this context, the analysis of the WSS topological skeleton is gaining more and more interest by extending the current understanding of the association between local hemodynamics and vascular diseases. The present study aims to analyze the impact that a deployed coronary stent has on the WSS topological skeleton. Computational fluid dynamics simulations were performed in three stented human coronary artery geometries reconstructed from clinical images. The selected cases presented stents with different designs (i.e., two contemporary drug eluting stents and one bioresorbable scaffold) and included regions with stent malapposition or overlapping. A recently proposed Eulerian-based approach was applied to analyze the WSS topological skeleton features. The results highlighted that the presence of single or multiple stents within a coronary artery markedly impacts the WSS topological skeleton. In particular, repetitive patterns of WSS divergence were observed at the luminal surface, highlighting a WSS contraction action proximal to the struts and a WSS expansion action distal to the struts. This WSS action pattern was independent from the stent design. In conclusions, these findings could contribute to a deeper understanding of the hemodynamic-driven processes underlying ST and ISR.

Secondary flow in bifurcations - Important effects of curvature, bifurcation angle and stents

Coronary bifurcations have complex flow patterns including secondary flow zones and helical flow, which directly affect pathophysiological mechanisms such as the development of atherosclerosis. The objective of this study was to generate insights into the effects of curvature, bifurcation angle and the presence of stents on flow patterns and resulting haemodynamics in coronary left main bifurcations. The blood flow and associated metrics were modelled in both idealised and patient-specific bifurcations with varying curvature and bifurcation angles with and without stents, resulting in a total of 128 geometries considered. The results showed that larger curvature of bifurcating vessels has a significant influence on secondary flow, especially with distance to the bifurcation region, causing a skew, spin and asymmetry of Dean vortices, an increase in helical flow intensity with symmetry loss, and a decrease in adversely low time-average wall shear stress (TAWSS). Generally, asymmetric flow patterns coincided with adversely low TAWSS regions. In identical stented geometries, the presence of the stents induced local recirculation immediately adjacent to the stent struts, thus generating adversely low TAWSS in these areas, with some effect on the overall secondary flow. Overall, the effect of stents outweighed the effect of curvature and BA. This new knowledge contributes to a better understanding of the joint effects of curvature, bifurcation angle, and stents on flow patterns and haemodynamics in coronary bifurcations.

Risk of myocardial infarction based on endothelial shear stress analysis using coronary angiography

BACKGROUND AND AIMS

Wall shear stress (WSS) has been associated with atherogenesis and plaque progression. The present study assessed the value of WSS analysis derived from conventional coronary angiography to detect lesions culprit for future myocardial infarction (MI).

METHODS AND RESULTS

Three-dimensional quantitative coronary angiography (3DQCA), was used to calculate WSS and pressure drop in 80 patients. WSS descriptors were compared between 80 lesions culprit of future MI and 108 non-culprit lesions (controls). Endothelium-blood flow interaction was assessed by computational fluid dynamics (10.8 ± 1.41 min per vessel). Median time between baseline angiography and MI was 25.9 (21.9-29.8) months. Mean patient age was 70.3 ± 12.7. Clinical presentation was STEMI in 35% and NSTEMI in 65%. Culprit lesions showed higher percent area stenosis (%AS), translesional vFFR difference (ΔvFFR), time-averaged WSS (TAWSS) and topological shear variation index (TSVI) compared to non-culprit lesions (p < 0.05 for all). TSVI was superior to TAWSS in predicting MI (AUC-TSVI = 0.77, 95%CI 0.71-0.84 vs. AUC-TAWSS = 0.61, 95%CI 0.53-0.69, p < 0.001). The addition of TSVI increased predictive and reclassification abilities compared to a model based on %AS and ΔvFFR (NRI = 1.04, p < 0.001, IDI = 0.22, p < 0.001).

CONCLUSIONS

A 3DQCA-based WSS analysis was feasible and can identify lesions culprit for future MI. The combination of area stenoses, pressure gradients and WSS predicted the occurrence of MI. TSVI, a novel WSS descriptor, showed strong predictive capacity to detect lesions prone to cause MI.

An in silico trials platform for the evaluation of effect of the arterial anatomy configuration on stent implantation. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:4213-4217. [PMID: 34892153 DOI: 10.1109/embc46164.2021.9629950] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The introduction of Bioresorbable Vascular Scaffolds (BVS) has revolutionized the treatment of atherosclerosis. InSilc is an in silico clinical trial (ISCT) platform in a Cloud-based environment used for the design, development and evaluation of BVS. Advanced multi-disciplinary and multiscale models are integrated in the platform towards predicting the short/acute and medium/long term scaffold performance. In this study, InSilc platform is employed in a use case scenario and demonstrates how the whole in silico pipeline allows the interpretation of the effect of the arterial anatomy configuration on stent implantation.

An optimal control approach to determine resistance-type boundary conditions from in-vivo data for cardiovascular simulations

The choice of appropriate boundary conditions is a fundamental step in computational fluid dynamics (CFD) simulations of the cardiovascular system. Boundary conditions, in fact, highly affect the computed pressure and flow rates, and consequently haemodynamic indicators such as wall shear stress (WSS), which are of clinical interest. Devising automated procedures for the selection of boundary conditions is vital to achieve repeatable simulations. However, the most common techniques do not automatically assimilate patient-specific data, relying instead on expensive and time-consuming manual tuning procedures. In this work, we propose a technique for the automated estimation of outlet boundary conditions based on optimal control. The values of resistive boundary conditions are set as control variables and optimized to match available patient-specific data. Experimental results on four aortic arches demonstrate that the proposed framework can assimilate 4D-Flow MRI data more accurately than two other common techniques based on Murray's law and Ohm's law.

Comparison of Swine and Human Computational Hemodynamics Models for the Study of Coronary Atherosclerosis

Coronary atherosclerosis is a leading cause of illness and death in Western World and its mechanisms are still non completely understood. Several animal models have been used to 1) study coronary atherosclerosis natural history and 2) propose predictive tools for this disease, that is asymptomatic for a long time, aiming for a direct translation of their findings to human coronary arteries. Among them, swine models are largely used due to the observed anatomical and pathophysiological similarities to humans. However, a direct comparison between swine and human models in terms of coronary hemodynamics, known to influence atherosclerotic onset/development, is still lacking. In this context, we performed a detailed comparative analysis between swine- and human-specific computational hemodynamic models of coronary arteries. The analysis involved several near-wall and intravascular flow descriptors, previously emerged as markers of coronary atherosclerosis initiation/progression, as well as anatomical features. To do that, non-culprit coronary arteries (18 right–RCA, 18 left anterior descending–LAD, 13 left circumflex–LCX coronary artery) from patients presenting with acute coronary syndrome were imaged by intravascular ultrasound and coronary computed tomography angiography. Similarly, the three main coronary arteries of ten adult mini-pigs were also imaged (10 RCA, 10 LAD, 10 LCX). The geometries of the imaged coronary arteries were reconstructed (49 human, 30 swine), and computational fluid dynamic simulations were performed by imposing individualized boundary conditions. Overall, no relevant differences in 1) wall shear stress-based quantities, 2) intravascular hemodynamics (in terms of helical flow features), and 3) anatomical features emerged between human- and swine-specific models. The findings of this study strongly support the use of swine-specific computational models to study and characterize the hemodynamic features linked to coronary atherosclerosis, sustaining the reliability of their translation to human vascular disease.

Patient-specific computational simulation of coronary artery bifurcation stenting

Patient-specific and lesion-specific computational simulation of bifurcation stenting is an attractive approach to achieve individualized pre-procedural planning that could improve outcomes. The objectives of this work were to describe and validate a novel platform for fully computational patient-specific coronary bifurcation stenting. Our computational stent simulation platform was trained using n = 4 patient-specific bench bifurcation models (n = 17 simulations), and n = 5 clinical bifurcation cases (training group, n = 23 simulations). The platform was blindly tested in n = 5 clinical bifurcation cases (testing group, n = 29 simulations). A variety of stent platforms and stent techniques with 1- or 2-stents was used. Post-stenting imaging with micro-computed tomography (μCT) for bench group and optical coherence tomography (OCT) for clinical groups were used as reference for the training and testing of computational coronary bifurcation stenting. There was a very high agreement for mean lumen diameter (MLD) between stent simulations and post-stenting μCT in bench cases yielding an overall bias of 0.03 (− 0.28 to 0.34) mm. Similarly, there was a high agreement for MLD between stent simulation and OCT in clinical training group [bias 0.08 (− 0.24 to 0.41) mm], and clinical testing group [bias 0.08 (− 0.29 to 0.46) mm]. Quantitatively and qualitatively stent size and shape in computational stenting was in high agreement with clinical cases, yielding an overall bias of < 0.15 mm. Patient-specific computational stenting of coronary bifurcations is a feasible and accurate approach. Future clinical studies are warranted to investigate the ability of computational stenting simulations to guide decision-making in the cardiac catheterization laboratory and improve clinical outcomes.

Three dimensional reconstruction of coronary artery stents from optical coherence tomography: experimental validation and clinical feasibility

The structural morphology of coronary stents (e.g. stent expansion, lumen scaffolding, strut apposition, tissue protrusion, side branch jailing, strut fracture), and the local hemodynamic environment after stent deployment are key determinants of procedural success and subsequent clinical outcomes. High-resolution intracoronary imaging has the potential to enable the geometrically accurate three-dimensional (3D) reconstruction of coronary stents. The aim of this work was to present a novel algorithm for 3D stent reconstruction of coronary artery stents based on optical coherence tomography (OCT) and angiography, and test experimentally its accuracy, reproducibility, clinical feasibility, and ability to perform computational fluid dynamics (CFD) studies. Our method has the following steps: 3D lumen reconstruction based on OCT and angiography, stent strut segmentation in OCT images, packaging, rotation and straightening of the segmented struts, planar unrolling of the segmented struts, planar stent wireframe reconstruction, rolling back of the planar stent wireframe to the 3D reconstructed lumen, and final stent volume reconstruction. We tested the accuracy and reproducibility of our method in stented patient-specific silicone models using micro-computed tomography (μCT) and stereoscopy as references. The clinical feasibility and CFD studies were performed in clinically stented coronary bifurcations. The experimental and clinical studies showed that our algorithm (1) can reproduce the complex spatial stent configuration with high precision and reproducibility, (2) is feasible in 3D reconstructing stents deployed in bifurcations, and (3) enables CFD studies to assess the local hemodynamic environment within the stent. Notably, the high accuracy of our algorithm was consistent across different stent designs and diameters. Our method coupled with patient-specific CFD studies can lay the ground for optimization of stenting procedures, patient-specific computational stenting simulations, and research and development of new stent scaffolds and stenting techniques.

Identification of the haemodynamic environment permissive for plaque erosion

Endothelial erosion of atherosclerotic plaques is the underlying cause of approximately 30% of acute coronary syndromes (ACS). As the vascular endothelium is profoundly affected by the haemodynamic environment to which it is exposed, we employed computational fluid dynamic (CFD) analysis of the luminal geometry from 17 patients with optical coherence tomography (OCT)-defined plaque erosion, to determine the flow environment permissive for plaque erosion. Our results demonstrate that 15 of the 17 cases analysed occurred on stenotic plaques with median 31% diameter stenosis (interquartile range 28–52%), where all but one of the adherent thrombi located proximal to, or within the region of maximum stenosis. Consequently, all flow metrics related to elevated flow were significantly increased (time averaged wall shear stress, maximum wall shear stress, time averaged wall shear stress gradient) with a reduction in relative residence time, compared to a non-diseased reference segment. We also identified two cases that did not exhibit an elevation of flow, but occurred in a region exposed to elevated oscillatory flow. Our study demonstrates that the majority of OCT-defined erosions occur where the endothelium is exposed to elevated flow, a haemodynamic environment known to evoke a distinctive phenotypic response in endothelial cells.

Effect of Material and Population on the Delivery of Nanoparticles to an Atherosclerotic Plaque: A Patient-specific In Silico Study

Coronary artery disease (CAD) is the prevalent reason of mortality all around the world. Targeting CAD, specifically atherosclerosis, with controlled delivery of micro and nanoparticles, as drug carriers, is a very proficient approach. In this work, a patient-specific and realistic model of an atherosclerotic plaque in the left anterior descending (LAD) artery was created by image-processing of CT-scan images and implementing a finite-element mesh. Next, a fluid-solid interaction simulation considering the physiological boundary conditions was conducted. By considering the simulated force fields and particle-particle interactions, the correlation between injected particles at each cardiac cycle and the surface density of adhered particles over the atherosclerotic plaque (SDP) were examined. For large particles (800 and 1000 nm) the amount of SDP on the plaque increased significantly when the number of the injected particles became higher. However, by increasing the number of the injected particles, for the larger particles (800 and 1000 nm) the increase in SDP was about 50% greater than that of the smaller ones (400 and 600 nm). Furthermore, for constant number of particles, depending on their size, different trends in SDP were observed. Subsequently, the distribution and adhesion of metal-based nanoparticles including SiO₂, Fe₃O₄, NiO₂, silver and gold with different properties were simulated. The injection of metal particles with medium density among the considered particles resulted in the highest SDP. Remarkably, the affinity, the geometrical features, and the biophysical factors involved in the adhesion outweighed the effect of difference in the density of particles on the SDP. Finally, the consideration of the lift force in the simulations significantly reduced the SDP and consistently decreased the particle residence time in the studied domain.

Blood Flow Modeling in Coronary Arteries: A Review

Atherosclerosis is one of the main causes of cardiovascular events, namely, myocardium infarction and cerebral stroke, responsible for a great number of deaths every year worldwide. This pathology is caused by the progressive accumulation of low-density lipoproteins, cholesterol, and other substances on the arterial wall, narrowing its lumen. To date, many hemodynamic studies have been conducted experimentally and/or numerically; however, this disease is not yet fully understood. For this reason, the research of this pathology is still ongoing, mainly, resorting to computational methods. These have been increasingly used in biomedical research of atherosclerosis because of their high-performance hardware and software. Taking into account the attempts that have been made in computational techniques to simulate realistic conditions of blood flow in both diseased and healthy arteries, the present review aims to give an overview of the most recent numerical studies focused on coronary arteries, by addressing the blood viscosity models, and applied physiological flow conditions. In general, regardless of the boundary conditions, numerical studies have been contributed to a better understanding of the development of this disease, its diagnosis, and its treatment.

Physiology and coronary artery disease: emerging insights from computed tomography imaging based computational modeling

Improvements in spatial and temporal resolution now permit robust high quality characterization of presence, morphology and composition of coronary atherosclerosis in computed tomography (CT). These characteristics include high risk features such as large plaque volume, low CT attenuation, napkin-ring sign, spotty calcification and positive remodeling. Because of the high image quality, principles of patient-specific computational fluid dynamics modeling of blood flow through the coronary arteries can now be applied to CT and allow the calculation of local lesion-specific hemodynamics such as endothelial shear stress, fractional flow reserve and axial plaque stress. This review examines recent advances in coronary CT image-based computational modeling and discusses the opportunity to identify lesions at risk for rupture much earlier than today through the combination of anatomic and hemodynamic information.

3D reconstruction of coronary artery bifurcations from coronary angiography and optical coherence tomography: feasibility, validation, and reproducibility

The three-dimensional (3D) representation of the bifurcation anatomy and disease burden is essential for better understanding of the anatomical complexity of bifurcation disease and planning of stenting strategies. We propose a novel methodology for 3D reconstruction of coronary artery bifurcations based on the integration of angiography, which provides the backbone of the bifurcation, with optical coherence tomography (OCT), which provides the vessel shape. Our methodology introduces several technical novelties to tackle the OCT frame misalignment, correct positioning of the OCT frames at the carina, lumen surface reconstruction, and merging of bifurcation lumens. The accuracy and reproducibility of the methodology were tested in n = 5 patient-specific silicone bifurcations compared to contrast-enhanced micro-computed tomography (µCT), which was used as reference. The feasibility and time-efficiency of the method were explored in n = 7 diseased patient bifurcations of varying anatomical complexity. The OCT-based reconstructed bifurcation models were found to have remarkably high agreement compared to the µCT reference models, yielding r² values between 0.91 and 0.98 for the normalized lumen areas, and mean differences of 0.005 for lumen shape and 0.004 degrees for bifurcation angles. Likewise, the reproducibility of our methodology was remarkably high. Our methodology successfully reconstructed all the patient bifurcations yielding favorable processing times (average lumen reconstruction time < 60 min). Overall, our method is an easily applicable, time-efficient, and user-friendly tool that allows accurate and reproducible 3D reconstruction of coronary bifurcations. Our technique can be used in the clinical setting to provide information about the bifurcation anatomy and plaque burden, thereby enabling planning, education, and decision making on bifurcation stenting.

Lipid-rich Plaques Detected by Near-infrared Spectroscopy Are More Frequently Exposed to High Shear Stress

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques.

Spherization of red blood cells and platelet margination in COPD patients

Red blood cells (RBCs) in pathological situations undergo biochemical and conformational changes, leading to alterations in rheology involved in cardiovascular events. The shape of RBCs in volunteers and stable and exacerbated chronic obstructive pulmonary disease (COPD) patients was analyzed. The effects of RBC spherization on platelet transport (displacement in the flow field caused by their interaction with RBCs) were studied in vitro and by numerical simulations. RBC spherization was observed in COPD patients compared with volunteers. In in vitro experiments at a shear rate of 100 s^-1 , treatment of RBCs with neuraminidase induced greater sphericity, which mainly affected platelet aggregates without changing aggregate size. At 400 s^-1 , neuraminidase treatment changes both the size of the aggregates and the number of platelet aggregates. Numerical simulations indicated that RBC spherization induces an increase of the platelet mean square displacement, which is traditionally linked to the platelet diffusion coefficient. RBCs of COPD patients are more spherical than healthy volunteers. Experimentally, RBC spherization induces increased platelet transport to the wall. Additional studies are needed to understand the link between the effect of RBCs on platelet transport and the increased cardiovascular events observed in COPD patients.

Endothelial shear stress and vascular remodeling in bioresorbable scaffold and metallic stent

BACKGROUND AND AIMS

The impact of endothelial shear stress (ESS) on vessel remodeling in vessels implanted with bioresorbable scaffold (BRS) as compared to metallic drug-eluting stent (DES) remains elusive. The aim of this study was to determine whether the relationship between ESS and remodeling patterns differs in BRS from those seen in metallic DES at 3-year follow-up.

METHODS

In the ABSORB II randomized trial, lesions were investigated by serial coronary angiography and intravascular ultrasound (IVUS). Three-dimensional reconstructions of coronary arteries post-procedure and at 3 years were performed. ESS was quantified using non-Newtonian steady flow simulation. IVUS cross-sections in device segment were matched using identical landmarks.

RESULTS

Paired ESS calculations post-procedure and at 3 years were feasible in 57 lesions in 56 patients. Post-procedure, median ESS at frame level was higher in BRS than in DES, with marginal statistical significance (0.97 ± 0.48 vs. 0.75 ± 0.39 Pa, p = 0.063). In the BRS arm, vessel area and lumen area showed larger increases in the highest tercile of median ESS post-procedure as compared to the lowest tercile. In contrast, in DES, no significant relationship between median ESS post-procedure and remodeling was observed. In multivariate analysis, smaller vessel area, larger lumen area, higher plaque burden post-procedure, and higher median ESS post-procedure were independently associated with expansive remodeling in matched frames. Only in BRS, younger age was an additional significant predictor of expansive remodeling.

CONCLUSIONS

In a subset of lesions with large plaque burden, shear stress could be associated with expansive remodeling and late lumen enlargement in BRS, while ESS had no impact on vessel dimension in metallic DES.

Usefulness of Proximal Coronary Wave Speed for Wave Intensity Analysis in Diseased Coronary Vessels

Background: Wave speed is needed to separate net wave intensity into forward and backward traveling components. However, wave speed in diseased coronary arteries cannot be assessed from hemodynamic measurements obtained distal to a stenosis. Wave speed inherently depends on arterial wall properties which should be similar proximal and distal to a stenosis. Our hypothesis is that proximal wave speed can be used to separate net wave intensity obtained distal to a stenosis.

Methods: We assessed coronary wave speed using the sum-of-squares single-point technique (SPc) based on simultaneous intracoronary pressure and flow velocity measurements in human coronary arteries. SPc at resting flow was determined in diseased coronary vessels of 12 patients both proximal and distal to the stenosis. In seven of these vessels, distal measurements were additionally obtained after revascularization by stent placement. SPc was also assessed at two axial locations in 14 reference vessels without a stenosis.

Results: (1) No difference in SPc was present between proximal and distal locations in the reference vessels. (2) In diseased vessels with a focal stenosis, SPc at the distal location was paradoxically larger than SPc proximal to the stenosis (28.4 ± 3.7 m/s vs. 18.3 ± 1.8 m/s, p < 0.02), despite the lower distending pressure downstream of the stenosis. The corresponding separated wave energy tended to be underestimated when derived from SPc at the distal compared with the proximal location. (3) After successful revascularization, SPc at the distal location no longer differed from SPc at the proximal location prior to revascularization (21.9 ± 2.0 m/s vs. 20.8 ± 1.9 m/s, p = 0.48). Accordingly, no significant difference in separated wave energy was observed for forward or backward waves.

Conclusion: In diseased coronary vessels, SPc assessed from distal hemodynamic signals is erroneously elevated. Our findings suggest that proximal wave speed can be used to separate wave intensity profiles obtained downstream of a stenosis. This approach may extend the application of wave intensity analysis to diseased coronary vessels.

The effects of clinically-derived parametric data uncertainty in patient-specific coronary simulations with deformable walls

Cardiovascular simulations are increasingly used for noninvasive diagnosis of cardiovascular disease, to guide treatment decisions, and in the design of medical devices. Quantitative assessment of the variability of simulation outputs due to input uncertainty is a key step toward further integration of cardiovascular simulations in the clinical workflow. In this study, we present uncertainty quantification in computational models of the coronary circulation to investigate the effect of uncertain parameters, including coronary pressure waveform, intramyocardial pressure, morphometry exponent, and the vascular wall Young's modulus. We employ a left coronary artery model with deformable vessel walls, simulated via an Arbitrary-Lagrangian-Eulerian framework for fluid-structure interaction, with a prescribed inlet pressure and open-loop lumped parameter network outlet boundary conditions. Stochastic modeling of the uncertain inputs is determined from intra-coronary catheterization data or gathered from the literature. Uncertainty propagation is performed using several approaches including Monte Carlo, Quasi Monte Carlo sampling, stochastic collocation, and multi-wavelet stochastic expansion. Variabilities in the quantities of interest, including branch pressure, flow, wall shear stress, and wall deformation are assessed. We find that uncertainty in inlet pressures and intramyocardial pressures significantly affect all resulting QoIs, while uncertainty in elastic modulus only affects the mechanical response of the vascular wall. Variability in the morphometry exponent used to distribute the total downstream vascular resistance to the single outlets, has little effect on coronary hemodynamics or wall mechanics. Finally, we compare convergence behaviors of statistics of QoIs using several uncertainty propagation methods on three model benchmark problems and the left coronary simulations. From the simulation results, we conclude that the multi-wavelet stochastic expansion shows superior accuracy and performance against Quasi Monte Carlo and stochastic collocation methods.

Does the inflow velocity profile influence physiologically relevant flow patterns in computational hemodynamic models of left anterior descending coronary artery?

Patient-specific computational fluid dynamics is a powerful tool for investigating the hemodynamic risk in coronary arteries. Proper setting of flow boundary conditions in computational hemodynamic models of coronary arteries is one of the sources of uncertainty weakening the findings of in silico experiments, in consequence of the challenging task of obtaining in vivo 3D flow measurements within the clinical framework. Accordingly, in this study we evaluated the influence of assumptions on inflow velocity profile shape on coronary artery hemodynamics. To do that, (1) ten left anterior descending coronary artery (LAD) geometries were reconstructed from clinical angiography, and (2) eleven velocity profiles with realistic 3D features such as eccentricity and differently shaped (single- and double-vortex) secondary flows were generated analytically and imposed as inflow boundary conditions. Wall shear stress and helicity-based descriptors obtained prescribing the commonly used parabolic velocity profile were compared with those obtained with the other velocity profiles. Our findings indicated that the imposition of idealized velocity profiles as inflow boundary condition is acceptable as long the results of the proximal vessel segment are not considered, in LAD coronary arteries. As a pragmatic rule of thumb, a conservative estimation of the length of influence of the shape of the inflow velocity profile on LAD local hemodynamics can be given by the theoretical entrance length for cylindrical conduits in laminar flow conditions.

Preclinical evaluation of a thin-strut bioresorbable scaffold (ArterioSorb): acute-phase invasive imaging assessment and hemodynamic implication

AIMS

The aim of this study was to assess the acute performance of the 95 µm ArterioSorb oriented poly L-lactic acid (PLLA) scaffold in comparison with the XIENCE metallic drug-eluting stent (DES) in porcine coronary arteries.

METHODS AND RESULTS

In 15 non-atherosclerotic Yucatan mini pigs, the ArterioSorb (3.0/14 mm) and XIENCE (3.0/15 mm) were implanted in 25 and 15 vessels, respectively. Acute performance was evaluated by using quantitative coronary angiography (QCA) and optical coherence tomography (OCT). Following three-dimensional reconstruction of the coronary arteries, endothelial shear stress (ESS) was quantified using non-Newtonian steady-flow simulation. Acute recoil measured by QCA was comparable in the two arms. Post-procedural flow and scaffold/stent area by OCT did not differ between the two devices. ESS post procedure was comparable between ArterioSorb and XIENCE (2.21±1.97 vs 2.25±1.71 Pa, p=0.314).

CONCLUSIONS

Acute recoil, luminal dimensions and ESS in the ArterioSorb oriented PLLA scaffold with thin struts of 95 µm were comparable to those in the XIENCE metallic DES.

A review study on blood in human coronary artery: Numerical approach

Computational fluid dynamics (CFD) study of blood flow in human coronary artery is one of the emerging fields of Biomed- ical engineering. In present review paper, Finite Volume Method with governing equations and boundary conditions are briefly discussed for different coronary models. Many researchers have come up with astonishing results related to the various factors (blood viscosity, rate of blood flow, shear stress on the arterial wall, Reynolds number, etc.) affecting the hemodynamic of blood in the right/left coronary artery. The aim of this paper is to present an overview of all those work done by the researchers to justify their work related to factors which hampers proper functioning of heart and lead to Coronary Artery Disease (CAD). Governing equations like Navier-stokes equations, continuity equations etc. are widely used and are solved using CFD solver to get a clearer view of coronary artery blockage. Different boundary conditions and blood properties published in the last ten years are summarized in the tabulated form. This table will help new researchers to work on this area.

The improvement of the shear stress and oscillatory shear index of coronary arteries during Enhanced External Counterpulsation in patients with coronary heart disease

Background

Enhanced External Counterpulsation (EECP) can chronically relieve ischemic chest pain and improve the prognosis of coronary heart disease (CHD). Despite its role in mitigating heart complications, EECP and the mechanisms behind its therapeutic nature, such as its effects on blood flow hemodynamics, are still not fully understood. This study aims to elucidate the effect of EECP on significant hemodynamic parameters in the coronary arterial tree.

Methods

A finite volume method was used in conjunction with the inlet pressure wave (surrogated by the measured aortic pressure) before and during EECP and outlet flow resistance, assuming the blood as newtonian fluid. The time-average wall shear stress (TAWSS) and oscillatory shear index (OSI) were determined from the flow field.

Results

Regardless of the degree of vascular stenosis, hemodynamic conditions and flow patterns could be improved during EECP. In comparison with the original tree, the tree with a severe stenosis (75% area stenosis) demonstrated more significant improvement in hemodynamic conditions and flow patterns during EECP, with surface area ratio of TAWSS risk area (SAR-TAWSS) reduced from 12.3% to 6.7% (vs. SAR-TAWSS reduced from 7.2% to 5.6% in the original tree) and surface area ratio of OSI risk area (SAR-OSI) reduced from 6.8% to 2.5% (vs. SAR-OSI of both 0% before and during EECP in the original tree because of mild stenosis). Moreover, it was also shown that small ratio of diastolic pressure (D) and systolic pressure (S) (D/S) could only improve the hemodynamic condition mildly. The SAR-TAWSS reduction ratio significantly increased as D/S became larger.

Conclusions

A key finding of the study was that the improvement of hemodynamic conditions along the LMCA trees during EECP became more significant with the increase of D/S and the severity degree of stenoses at the bifurcation site. These findings have important implications on EECP as adjuvant therapy before or after percutaneous coronary intervention (PCI) in patients with diffuse atherosclerosis.

Topologic and Hemodynamic Characteristics of the Human Coronary Arterial Circulation

Background

Many processes contributing to the functional and structural regulation of the coronary circulation have been identified. A proper understanding of the complex interplay of these processes requires a quantitative systems approach that includes the complexity of the coronary network. The purpose of this study was to provide a detailed quantification of the branching characteristics and local hemodynamics of the human coronary circulation.

Methods

The coronary arteries of a human heart were filled post-mortem with fluorescent replica material. The frozen heart was alternately cut and block-face imaged using a high-resolution imaging cryomicrotome. From the resulting 3D reconstruction of the left coronary circulation, topological (node and loop characteristics), topographic (diameters and length of segments), and geometric (position) properties were analyzed, along with predictions of local hemodynamics (pressure and flow).

Results

The reconstructed left coronary tree consisted of 202,184 segments with diameters ranging from 30 μm to 4 mm. Most segments were between 100 μm and 1 mm long. The median segment length was similar for diameters ranging between 75 and 200 μm. 91% of the nodes were bifurcations. These bifurcations were more symmetric and less variable in smaller vessels. Most of the pressure drop occurred in vessels between 200 μm and 1 mm in diameter. Downstream conductance variability affected neither local pressure nor median local flow and added limited extra variation of local flow. The left coronary circulation perfused 358 cm³ of myocardium. Median perfused volume at a truncation level of 100 to 200 μm was 20 mm³ with a median perfusion of 5.6 ml/min/g and a high local heterogeneity.

Conclusion

This study provides the branching characteristics and hemodynamic analysis of the left coronary arterial circulation of a human heart. The resulting model can be deployed for further hemodynamic studies at the whole organ and local level.

Relationship Between Stent Diameter, Platelet Reactivity, and Thrombotic Events After Percutaneous Coronary Artery Revascularization

Small vessel diameter and residual platelet reactivity are independent predictors of thrombotic events after percutaneous coronary intervention (PCI). We sought to determine whether an interaction exists between residual platelet reactivity and stent diameter regarding the occurrence of stent thrombosis and other adverse events after PCI. We stratified patients in the prospective ADAPT-DES registry who underwent single-lesion PCI according to if they received a small diameter stent (SDS, defined as a stent with a diameter of 2.25 mm). Patients receiving an SDS were compared with patients receiving a stent ≥2.5 mm using Kaplan-Meier rates and multivariable Cox proportional hazards regression. We defined major adverse cardiac events (MACE) as the composite of cardiac death, myocardial infarction, and stent thrombosis (ST). Among 5,608 patients who underwent single-lesion PCI in ADAPT-DES, 222 (4.0%) patients received an SDS. Patients with an SDS were more likely than patients without an SDS to have 3-vessel disease but received, on average, fewer stents and were less likely to present with a thrombotic lesion. Receiving versus not receiving an SDS was associated with increased risk of ST (adjusted hazard ratio 4.35, 95% confidence interval 1.95 to 9.73, p <0.001) as well as MACE (adjusted hazard ratio 1.75, 95% confidence interval 1.11 to 2.75, p = 0.02). There was no statistical interaction between platelet reactivity and SDS regarding ST (p = 0.12) or MACE (p = 0.51). In conclusion, PCI with small drug-eluting stents is associated with a high risk of thrombotic events, including ST. Further studies should explore whether alternative treatment strategies are appropriate in small vessels.

Computational instantaneous wave-free ratio (IFR) for patient-specific coronary artery stenoses using 1D network models

In this work, we estimate the diagnostic threshold of the instantaneous wave-free ratio (iFR) through the use of a one-dimensional haemodynamic framework. To this end, we first compared the computed fractional flow reserve (cFFR) predicted from a 1D computational framework with invasive clinical measurements. The framework shows excellent promise and utilises minimal patient data from a cohort of 52 patients with a total of 66 stenoses. The diagnostic accuracy of the cFFR model was 75.76%, with a sensitivity of 71.43%, a specificity of 77.78%, a positive predictive value of 60%, and a negative predictive value of 85.37%. The validated model was then used to estimate the diagnostic threshold of iFR. The model determined a quadratic relationship between cFFR and the ciFR. The iFR diagnostic threshold was determined to be 0.8910 from a receiver operating characteristic curve that is in the range of 0.89 to 0.9 that is normally reported in clinical studies.

Transfer of Low-Density Lipoproteins in Coronary Artery Bifurcation Lesions with Stenosed Side Branch: Numerical Study

Evidence from clinical data suggests that the stenotic side branch (SB) is one of the key predictors for SB occlusion-based adverse events. In this study, we hypothesized that coronary bifurcations with stenotic SB might lead to severe concentration polarization of atherogenic lipids, such as the low-density lipoproteins (LDL), motivating the adverse events in the clinic. To confirm this hypothesis, this work numerically investigated the transport of LDL in different bifurcation lesions based on the Medina classification with various location and stenosis severities. The results showed that the coronary bifurcations with stenotic SB might be suffering more serious concentration polarization of LDL on the luminal surface of the SB due to higher level of LDL concentrations. Moreover, compared to the other bifurcation lesion types, the type (1,0,1) had the highest luminal surface LDL concentration along the SB and the highest degree of risk to enhance the process of atherosclerosis. In addition, this study also showed that the luminal surface LDL concentration increased with elevated stenosis severity. The type (1,0,1) with the severe stenosis (75% diameter reduction) had the highest concentration at the SB. In conclusion, these results suggested that both location of lesions and stenosis severities had great influence on the distribution of LDL on the luminal surface of the SB. Therefore, the estimation of disease severity and the interventional therapy should be carried out not only according to the stenosis severities in clinic. Moreover, compared to the other bifurcation lesion types, the type (1,0,1), rather than the type (1,1,1) as usually considered, had the highest luminal surface LDL concentration along the SB and the highest degree of risk to enhance the process of atherosclerosis.