Healthy and diseased coronary bifurcation geometries influence near-wall and intravascular flow: A computational exploration of the hemodynamic risk

Residence time in complex left main bifurcation disease after stenting

BACKGROUND

Data regarding the mean resident time (RT) after left main (LM) bifurcation stenting are scant. In the present study we performed a patient-specific computational fluid dynamic (CFD) analysis to investigate the different post-stenting mean RT values in LM patients treated with single-or double stenting techniques.

METHODS

Patients were identified after reviewing the local Optical Coherence Tomography (OCT) scans database. Overall, 27 patients (mean age 65.5 ± 12.4, 21 males) [10 patients treated with provisional cross-over stenting, 7 with the double kissing crush (DK crush) and 10 with the nano-inverted T (NIT) technique, respectively] with isolated and significant LM bifurcation disease were analyzed.

RESULTS

After LM bifurcation stenting, the NIT showed a higher averages WSS values at all bifurcation sites compared to DK crush and provisional cross-over stenting. Moreover, the mean RT resulted lower after NIT compared to provisional or DK crush. During the diastolic phase, the average RT of the entire LM bifurcation was 0.46 s, 0.38 s and 0.33 s after using the provisional stenting, DK crush and NIT, respectively. Moreover, the average RT in the LM bifurcation decreased by 17.1 % using the DK crush and by 28.2 % using the NIT compared to the Provisional.

CONCLUSION

The present OCT-derived CFD analysis revealed that, in patients with complex bifurcation LM disease, the provisional approach resulted in lower WSS values, while double stenting techniques, especially the NIT technique, resulted in a marked reduction of average RT compared to the provisional approach.

CONDENSED ABSTRACT

In the present study we performed a patient-specific Optical coherence tomography (OCT)-based computational fluid dynamic (CFD) analysis to investigate the different post-stenting mean RT values in 27 patients treated with provisional cross-over stenting, DK crush and Nano-inverted-T (NIT) stenting. The NIT showed a higher averages WSS values at all bifurcation sites compared to DK crush and Provisional. The mean RT resulted lower in NIT compared to Provisional or DK crush. During the entire diastolic phase, the average RT of the entire LM bifurcation was 0.46 s, 0.38 s and 0.33 s after using the provisional stenting, DK crush and NIT, respectively. Moreover, the average RT in the entire LM bifurcation decreased by 17.1 % using the DK crush and by 28.2 % using the NIT compared to the Provisional.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment.Methods. According to Taguchi's method of design of experiment (DoE), the left main stem (LMS) length (lLMS), curvature (kLMS), diameter (dLMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (αLAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed.Results. The proximal LAD (LADproximal) was identified to be the most atheroprone region of the left coronary artery due to higherh2,large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa-1). In both patient-specific cases, based onh2and TAWSS,dlmsis the dominant geometric parameter while based on OSI and RRT,αLAD-LCxis the dominant one influencing hemodynamic condition in proximal LAD (p< 0.05). Based on RRT, the rank of the geometrical factors is:αLAD-LCx>dLMS>lLMS>kLMS, indicating thatαLAD-LCxis the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis.Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

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.

Left main coronary artery morphological phenotypes and its hemodynamic properties

BACKGROUND

Atherosclerosis may be linked to morphological defects that lead to variances in coronary artery hemodynamics. Few objective strategies exit at present for generalizing morphological phenotypes of coronary arteries in terms of hemodynamics. We used unsupervised clustering (UC) to classify the morphology of the left main coronary artery (LM) and looked at how hemodynamic distribution differed between phenotypes.

METHODS

In this study, 76 LMs were obtained from 76 patients. After LMs were reconstructed with coronary computed tomography angiography, centerlines were used to extract the geometric characteristics. Unsupervised clustering was carried out using these characteristics to identify distinct morphological phenotypes of LMs. The time-averaged wall shear stress (TAWSS) for each phenotype was investigated by means of computational fluid dynamics (CFD) analysis of the left coronary artery.

RESULTS

We identified four clusters (i.e., four phenotypes): Cluster 1 had a shorter stem and thinner branches (n = 26); Cluster 2 had a larger bifurcation angle (n = 10); Cluster 3 had an ostium at an angulation to the coronary sinus and a more curved stem, and thick branches (n = 10); and Cluster 4 had an ostium at an angulation to the coronary sinus and a flatter stem (n = 14). TAWSS features varied widely across phenotypes. Nodes with low TAWSS (L-TAWSS) were typically found around the branching points of the left anterior descending artery (LAD), particularly in Cluster 2.

CONCLUSION

Our findings demonstrated that UC is a powerful technique for morphologically classifying LMs. Different LM phenotypes exhibited distinct hemodynamic characteristics in certain regions. This morphological clustering method could aid in identifying people at high risk for developing coronary atherosclerosis, hence facilitating early intervention.

Visualization of multiphase pulsatile blood over single phase blood flow in a patient specific stenosed left coronary artery using image processing technique

BACKGROUND

Coronary arteries disease has been reported as one of the principal roots of deaths worldwide.

OBJECTIVE

The aim of this study is to analyze the multiphase pulsatile blood flow in the left coronary artery tree with stenosis.

METHODS

The 3D left coronary artery model was reconstructed using 2D computerized tomography (CT) scan images. The Red Blood Cell (RBC) and varying hemodynamic parameters for single and multiphase blood flow conditions were analyzed.

RESULTS

Results asserted that the multiphase blood flow modeling has a maximum velocity of 1.017 m/s and1.339 m/s at the stenosed region during the systolic and diastolic phases respectively. The increase in Wall Shear Stress (WSS) observed at the stenosed region during the diastole phase as compared during the systolic phase. It was also observed that the highest Oscillatory Shear Index (OSI) regions are found in the downstream area of stenosis and across the bifurcations. The increase in RBCs velocity from 0.45 m/s to 0.6 m/s across the stenosis was also noticed.

CONCLUSION

The computational multiphase blood flow analysis improves the understanding and accuracy of the complex flow conditions of blood elements (RBC and Plasma) and provides the progression of the disease development in the coronary arteries. This study helps to enhance the diagnosis of the blocked (stenosed) arteries more precisely compared to the single-phase blood flow modeling.

Patient-Specific Image-Based Computational Fluid Dynamics Analysis of Abdominal Aorta and Branches

The complicated abdominal aorta and its branches are a portion of the circulatory system prone to developing atherosclerotic plaque and aneurysms. These disorders are closely connected to the changing blood flow environment that the area’s complicated architecture produces (between celiac artery and iliac artery bifurcation); this phenomenon is widespread at arterial bifurcations. Based on computed tomography angiography (CTA) scans, this current work offers a numerical analysis of a patient-specific reconstruction of the abdominal aorta and its branches to identify and emphasize the most likely areas to develop atherosclerosis. The simulations were run following the heart cycle and under physiological settings. The wall shear stress (WSS), velocity field, and streamlines were examined. According to the findings, complex flow is primarily present at the location of arterial bifurcations, where abnormal flow patterns create recirculation zones with low and fluctuating WSS (<0.5 Pa), which are known to affect endothelial homeostasis and cause adverse vessel remodeling. The study provides a patient-specific hemodynamic analysis model, which couples in vivo CT imaging with in silico simulation under physiological circumstances. The study offers quantitative data on the range fluctuations of important hemodynamic parameters, such as WSS and recirculation region expansion, which are directly linked to the onset and progression of atherosclerosis. The findings could also help drug targeting at this vascular level by understanding blood flow patterns in the abdominal aorta and its branches.

Relationship between Coronary Arterial Geometry and the Presence and Extend of Atherosclerotic Plaque Burden: A Review Discussing Methodology and Findings in the Era of Cardiac Computed Tomography Angiography

Coronary artery disease (CAD) represents a modern pandemic associated with significant morbidity and mortality. The multi-faceted pathogenesis of this entity has long been investigated, highlighting the contribution of systemic factors such as hyperlipidemia and hypertension. Nevertheless, recent research has drawn attention to the importance of geometrical features of coronary vasculature on the complexity and vulnerability of coronary atherosclerosis. Various parameters have been investigated so far, including vessel-length, coronary artery volume index, cross-sectional area, curvature, and tortuosity, using primarily invasive coronary angiography (ICA) and recently non-invasive cardiac computed tomography angiography (CCTA). It is clear that there is correlation between geometrical parameters and both the haemodynamic alterations augmenting the atherosclerosis-prone environment and the extent of plaque burden. The purpose of this review is to discuss the currently available literature regarding this issue and propose a potential non-invasive imaging biomarker, the geometric risk score, which could be of importance to allow the early detection of individuals at increased risk of developing CAD.

Predicting hemodynamic indices in coronary artery aneurysms using response surface method: An application in Kawasaki disease

BACKGROUND AND OBJECTIVES

Coronary artery aneurysms (CAA), such as those in Kawasaki Disease (KD), induce hemodynamic alterations associated with thrombosis and atherosclerosis. Current clinical routines assess the risk level of the CAA cases based on the Z-Score, which considers the body surface area (BSA) and the CAA's diameter. A full geometric characterization and impact on hemodynamic metrics and their correlation with thrombotic risks have not been systematically investigated. The goal of this study was to investigate the effect of CAA shape indices on local hemodynamics using the response surface method (RSM) through considering KD applications.

METHODS

Transient computational fluid dynamics (CFD) simulations have been performed on idealized CAA geometries defined by geometrical ratios combining neck diameter, CAA diameter and CAA length. The results were used to develop full quadratic regression models of the indices using the response surface method (RSM). Validation using patient-specific KD models was performed.

RESULTS

The results indicated that the aneurysm diameter is the main determining factor in the thrombotic risk of CAA patients, which is consistent with clinical guidelines. Furthermore, it was observed that in most CAA cases having the same diameter, the one with the shorter length experiences higher RRT values, indicating flow stagnation and circulation.

CONCLUSIONS

The developed regression models can be used to ultimately assess the thrombotic risk of CAA cases from the hemodynamic perspective. The applicability of these models was tested on 2 KD patient specific models, with close values achieved between the models and the patient-specific results.

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.

A new and automated risk prediction of coronary artery disease using clinical endpoints and medical imaging-derived patient-specific insights: protocol for the retrospective GeoCAD cohort study

INTRODUCTION

Coronary artery disease (CAD) is the leading cause of death worldwide. More than a quarter of cardiovascular events are unexplained by current absolute cardiovascular disease risk calculators, and individuals without clinical risk factors have been shown to have worse outcomes. The 'anatomy of risk' hypothesis recognises that adverse anatomical features of coronary arteries enhance atherogenic haemodynamics, which in turn mediate the localisation and progression of plaques. We propose a new risk prediction method predicated on CT coronary angiography (CTCA) data and state-of-the-art machine learning methods based on a better understanding of anatomical risk for CAD. This may open new pathways in the early implementation of personalised preventive therapies in susceptible individuals as a potential key in addressing the growing burden of CAD.

METHODS AND ANALYSIS

GeoCAD is a retrospective cohort study in 1000 adult patients who have undergone CTCA for investigation of suspected CAD. It is a proof-of-concept study to test the hypothesis that advanced image-derived patient-specific data can accurately predict long-term cardiovascular events. The objectives are to (1) profile CTCA images with respect to variations in anatomical shape and associated haemodynamic risk expressing, at least in part, an individual's CAD risk, (2) develop a machine-learning algorithm for the rapid assessment of anatomical risk directly from unprocessed CTCA images and (3) to build a novel CAD risk model combining traditional risk factors with these novel anatomical biomarkers to provide a higher accuracy CAD risk prediction tool.

ETHICS AND DISSEMINATION

The study protocol has been approved by the St Vincent's Hospital Human Research Ethics Committee, Sydney-2020/ETH02127 and the NSW Population and Health Service Research Ethics Committee-2021/ETH00990. The project outcomes will be published in peer-reviewed and biomedical journals, scientific conferences and as a higher degree research thesis.

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.

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.

The effects of plaque morphological characteristics on the post-stenotic flow in left main coronary artery bifurcation

Local post-stenotic hemodynamics has critical influence in the atherosclerotic plaque progression occurring in susceptible arterial sites, in particular the left main coronary artery (LMCA) bifurcation. Understanding the effects of plaque morphological characteristics: stenosis severity (SS), eccentricity index (EI) and lesion length (LL) on the post-stenotic flow behavior can significantly improve treatment planning. In order to investigate these effects, we have employed computational fluid dynamics (CFD) simulations in twenty computer-generated and five patient-specific LMCA models and the hemodynamic parameters: velocity, pressure (P), wall pressure gradient (WPG), wall shear stress (WSS), time averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) and helicity intensity (h₂) were analyzed. Our results revealed that the effect of stenosis eccentricity varied significantly for different values of stenosis severity and lesion length. Regions with low WSS, low TAWSS and high RRT were more prominent in models having higher stenosis severity. For smaller lesion length, at low and moderate stenosis severity, surface area with low TAWSS and high RRT decreased with increasing eccentricity index, whereas for high stenosis severity models, low TAWSS region and average RRT values increased with eccentricity. However, for models with longer lesion length, regions with high OSI and RRT overall increased gradually with eccentricity. The helicity intensity (h₂) of all models remained very low except at the most eccentric model with longer lesion length. The presence of very high helical flow in this model suggests the possibility of atheroprotective flow. It can be concluded that all plaque morphological characteristics covered under this investigation play an important role in plaque progression.

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.

Effect of transport parameters on atherosclerotic lesion growth: A parameter sensitivity analysis

BACKGROUND

Atherosclerosis is a degenerative disease of the arterial wall. It results in the formation of progressively growing plaque lesions that can harden and narrow their host arteries. Current computational models of the inflammatory process that govern atherosclerosis growth are reliant on a number of parameters that can freely vary and whose precise values are not well known.

METHODS

To identify the significance of variation in such parameters, a parametric sensitivity study had been conducted on the blood density, blood viscosity, plasma viscosity and bulk flow low density lipoprotein (LDL) concentration. Using computational modeling, the significance of variation in these parameters was assessed on the transport of LDL. The simulation was performed via the 2^k factorial experimental design, which was conducted to identify the significance of the select parameters on the intima LDL concentration and endothelial LDL coverage area.

RESULTS

Results identified the blood viscosity and bulk flow LDL concentration are the dominant parameters for the atherosclerotic lesion growth. The coverage of LDL on the arterial wall surface was strongly dependent on the blood viscosity. The significance of these findings was discussed.

CONCLUSION

This statistical study identifies two dominating blood factors, LDL concentration and blood viscosity, and how they influence atherosclerosis which will serves as a guideline for further investigation on the atherosclerosis topic.

Risk prediction of in-stent restenosis among patients with coronary drug-eluting stents: current clinical approaches and challenges

Introduction: In-stent restenosis (ISR) has been one of the biggest limitations to the success of percutaneous coronary intervention for the treatment of coronary artery disease (CAD). The introduction of drug-eluting stent (DES) was a revolution in the treatment of CAD because these devices drastically reduced ISR to very low levels (<5%). Subsequently, newer generation DES treatments have overcome the drawbacks of first-generation DES, i.e. delayed endothelialization, and late stent thrombosis. However, the issue of late ISR, including neoatherosclerosis after DES implantation especially in high-risk patients and complex lesions, still exists as a challenge to be overcome.Areas covered: We discuss the mechanisms of ISR development including neoatherosclerosis, past and current clinical status of ISR, and methods to predict and overcome this issue from pathological and clinical points of view.Expert opinion: The initial drawbacks of first-generation DES, such as delayed endothelial healing and subsequent risk of late stent thrombosis, have been improved upon by the current generation DES. To achieve better long-term clinical outcomes, further titration of drug-release and polymer degradation profile, strut thickness as well as material innovation are needed.

High precision invasive FFR, low-cost invasive iFR, or non-invasive CFR?: optimum assessment of coronary artery stenosis based on the patient-specific computational models

The objective of this paper is to apply computational fluid dynamic (CFD) as a complementary tool for clinical tests to not only predict the present and future status of left coronary artery stenosis but also to evaluate some clinical hypotheses. In order to assess the present status of the coronary artery stenosis severity, and thereby selecting the most appropriate type of treatment for each patient, fractional flow reserve (FFR), instantaneous wave free-ratio (iFR), and coronary flow reserve (CFR) are calculated. To examine FFR, iFR, and CFR results, the effect of geometric features of stenoses, including diameter reduction (%), lesion length (LL), and minimum lumen diameter (MLD), is studied on them. It is observed that FFR is a more conservative index than iFR and CFR to assess the severity of coronary stenosis. In addition, it is seen that FFR, iFR, and CFR decrease by increasing LL and decreasing MLD. Therefore, the morphological indices, LL/MLD and LL/MLD̂4, with the calculated conservative cut-off values equal to 5.5 and 3.6, are considered. Next, some controversial clinical hypotheses about the assessment of the severity of coronary stenosis are evaluated numerically. These include the examination of FFR, iFR, and CFR accuracies, investigating the effect of coronary hyperemia on iFR, as well as the reliability of the hybrid iFR-FFR decision-making strategy. The presented numerical model can also be used as a predictive tool to identify the atherosuseptible sites of arteries by calculating the time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT).

Haemodynamic mechanism of formation and distribution of coronary atherosclerosis: A lesion-specific model

Coronary arterial disease, as the most devastated cardiovascular disease, is caused by the atherosclerosis in the coronary arteries, which blocks the blood flow to the heart, resulting in the deficient supply of oxygen and nutrition to the heart, and eventually leading to heart failure. To date, haemodynamic mechanisms for atherosclerosis development are not fully understood although it is believed that the haemodynamic disturbance at the region of the arterial bifurcation, particular, bifurcation angle, plays an important role in the atherosclerosis development. In this study, two types of computational fluid dynamics models, lesion-specific and idealized models, combined with the computer tomography imaging techniques, are used to explore the mechanism of formation and distribution of the atherosclerosis around the bifurcation of left coronary artery and its association with the bifurcation angle. The lesion-specific model is used to characterize the effect of personalized features on the haemodynamic performance, while the idealized model is focusing on the effect of single factor, bifurcation angle, on the haemodynamic performance. The simulated results from both types of the models, combined with the clinical observation, revealed that the three key areas around the bifurcations are prone to formation of the atherosclerosis. Unlike the idealized models, lesion-specific modelling results did not show the significant correlation between the wall shear stress and bifurcation angle, although the mean value of the wall shear stress in smaller bifurcation angles (less than 90°) is higher than that with larger bifurcation angles (greater than 90°). In conclusion, lesion-specific computational fluid dynamics modelling is an efficient and convenient way to predict the haemodynamic performance around the bifurcation region, allowing the comprehensive information for the clinicians to predict the atherosclerosis development. The idealized models, which only focus on single parameter, may not provide the sufficient and reliable information for the clinical application. A novel multi-parameters modelling technique, therefore, is suggested to be developed in future, allowing the effects of many parameters on the haemodynamic performance to be evaluated.

Assessment of surface roughness and blood rheology on local coronary haemodynamics: a multi-scale computational fluid dynamics study

The surface roughness of the coronary artery is associated with the onset of atherosclerosis. The study applies, for the first time, the micro-scale variation of the artery surface to a 3D coronary model, investigating the impact on haemodynamic parameters which are indicators for atherosclerosis. The surface roughness of porcine coronary arteries have been detailed based on optical microscopy and implemented into a cylindrical section of coronary artery. Several approaches to rheology are compared to determine the benefits/limitations of both single and multiphase models for multi-scale geometry. Haemodynamic parameters averaged over the rough/smooth sections are similar; however, the rough surface experiences a much wider range, with maximum wall shear stress greater than 6 Pa compared to the approximately 3 Pa on the smooth segment. This suggests the smooth-walled assumption may neglect important near-wall haemodynamics. While rheological models lack sufficient definition to truly encompass the micro-scale effects occurring over the rough surface, single-phase models (Newtonian and non-Newtonian) provide numerically stable and comparable results to other coronary simulations. Multiphase models allow for phase interactions between plasma and red blood cells which is more suited to such multi-scale models. These models require additional physical laws to govern advection/aggregation of particulates in the near-wall region.

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.

Exploring wall shear stress spatiotemporal heterogeneity in coronary arteries combining correlation-based analysis and complex networks with computational hemodynamics

Atherosclerosis at the early stage in coronary arteries has been associated with low cycle-average wall shear stress magnitude. However, parallel to the identification of an established active role for low wall shear stress in the onset/progression of the atherosclerotic disease, a weak association between lesions localization and low/oscillatory wall shear stress has been observed. In the attempt to fully identify the wall shear stress phenotype triggering early atherosclerosis in coronary arteries, this exploratory study aims at enriching the characterization of wall shear stress emerging features combining correlation-based analysis and complex networks theory with computational hemodynamics. The final goal is the characterization of the spatiotemporal and topological heterogeneity of wall shear stress waveforms along the cardiac cycle. In detail, here time-histories of wall shear stress magnitude and wall shear stress projection along the main flow direction and orthogonal to it (a measure of wall shear stress multidirectionality) are analyzed in a representative dataset of 10 left anterior descending pig coronary artery computational hemodynamics models. Among the main findings, we report that the proposed analysis quantitatively demonstrates that the model-specific inlet flow-rate shapes wall shear stress time-histories. Moreover, it emerges that a combined effect of low wall shear stress magnitude and of the shape of the wall shear stress–based descriptors time-histories could trigger atherosclerosis at its earliest stage. The findings of this work suggest for new experiments to provide a clearer determination of the wall shear stress phenotype which is at the basis of the so-called arterial hemodynamic risk hypothesis in coronary arteries.

Comparison of bioresorbable vs durable polymer drug-eluting stents in unprotected left main (from the RAIN-CARDIOGROUP VII Study)

BACKGROUND

There are limited data regarding the impact of bioresorbable polymer drug eluting stent (BP-DES) compared to durable polymer drug eluting stent (DP-DES) in patients treated with percutaneous coronary intervention using ultrathin stents in left main or bifurcations.

METHODS

In the RAIN registry (ClinicalTrials NCT03544294, june 2018 retrospectively registered) patients with a ULM or bifurcation stenosis treated with PCI using ultrathin stents (struts thinner than 81 μm) were enrolled. The primary endpoint was the rate of target lesion revascularization (TLR); major adverse cardiovascular events (MACE, a composite of all-cause death, myocardial infarction, TLR and stent thrombosis) and its components, along with target vessel revascularization (TVR) were the secondary ones. A propensity score with matching analysis to compare patients treated with BP-DES versus DP-DES was also assessed.

RESULTS

From 3001 enrolled patients, after propensity score analysis 1400 patients (700 for each group) were selected. Among them, 352 had ULM disease and 1048 had non-LM bifurcations. At 16 months (12-22), rates of TLR (3.7% vs 2.9%, p = 0.22) and MACE were similar (12.3% vs. 11.6%, p = 0.74) as well as for the other endpoints. Sensitivity analysis of outcomes after a two-stents strategy, showed better outcome in term of MACE (20.4% vs 10%, p = 0.03) and TVR (12% vs 4.6%, p = 0.05) and a trend towards lower TLR in patients treated with BP-DES.

CONCLUSION

In patients with bifurcations or ULM treated with ultrathin stents BP-DES seems to perform similarly to DP-DES: the trends toward improved clinical outcomes in patients treated with the BP-DES might potentially be of value for speculating the stent choice in selected high-risk subgroups of patients at increased risk of ischemic events.

TRIAL REGISTRATION

ClinicalTrials.gov Identifier: NCT03544294. Retrospectively registered June 1, 2018.

A framework for automated and objective modification of tubular structures: Application to the internal carotid artery

Patient-specific medical image-based computational fluid dynamics has been widely used to reveal fundamental insight into mechanisms of cardiovascular disease, for instance, correlating morphology to adverse vascular remodeling. However, segmentation of medical images is laborious, error-prone, and a bottleneck in the development of large databases that are needed to capture the natural variability in morphology. Instead, idealized models, where morphological features are parameterized, have been used to investigate the correlation with flow features, but at the cost of limited understanding of the complexity of cardiovascular flows. To combine the advantages of both approaches, we developed a tool that preserves the patient-specificness inherent in medical images while allowing for parametric alteration of the morphology. In our open-source framework morphMan we convert the segmented surface to a Voronoi diagram, modify the diagram to change the morphological features of interest, and then convert back to a new surface. In this paper, we present algorithms for modifying bifurcation angles, location of branches, cross-sectional area, vessel curvature, shape of bends, and surface roughness. We show qualitative and quantitative validation of the algorithms, performing with an accuracy exceeding 97% in general, and proof-of-concept on combining the tool with computational fluid dynamics. By combining morphMan with appropriate clinical measurements, one could explore the morphological parameter space and resulting hemodynamic response using only a handful of segmented surfaces, effectively minimizing the main bottleneck in image-based computational fluid dynamics.

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.

Computational fluid dynamic-derived wall shear stress of non-significant left main bifurcation disease may predict acute vessel thrombosis at 3-year follow-up

Wall shear stress (WSS) plays a pivotal role on plaque progression in coronary artery disease. We assess the prognostic role of baseline mean WSS in developing a bifurcation-located myocardial infarction (B-MI) over the following 3 years in angiographically non-significant LM bifurcation disease. For this purpose, we retrospectively reviewed the procedural and medical records of consecutive patients evaluated in our center from 1st January 2014 to 1st January 2019 who had a non-significant LM bifurcation disease as evaluated at coronary computed tomography angiography (CCTA) and confirmed by coronary angiography. Each bifurcation model was reconstructed on the patient-specific geometries derived from the CCTA. The population was divided into two groups: patients with (n = 12) and without B-MI (n = 20) over the following 3 years. Both the mean WSSprox of each branch and the WSSentire_lesion of each vessel, adjusted for the respective mean lesions lengths and 3-dimensional percentage of stenosis (DS%), resulted in independent predictors of 3-year B-MI. Multivariate Cox-regression analysis confirmed that a baseline mean WSSentire_model ≥ 5.05 Pa (HR 1.98, 95% CI 1.83-2.10, p = 0.001) was a predictor of 3-year B-MI independently from the entire mean lesions lengths (HR 1.56. 95% CI 1.43.1.68, p = 0.002) and DS% (HR 1.26, 95% CI 1.18-1.37, p = 0.03). In conclusion, in patients with angiographically non-significant LM bifurcation disease, both the mean WSSprox of each branch and WSSentire_lesion of each stenotic vessel predicted the occurrence of B-MI over the following 3 years. Moreover, the WSSentire_bifurcation ≥ 5.05 Pa seems to be a predictor of 3-year B-MI independently from the DS% and lesions lengths.

Histopathologic and physiologic effect of bifurcation stenting: current status and future prospects

Introduction: Coronary bifurcation lesions are involved in up to 20% of all percutaneous coronary interventions (PCI). However, bifurcation lesion intervention is associated with a high complication rate, and optimal treatment of coronary bifurcation is an ongoing debate.Areas covered: Both different stenting techniques and a variety of devices have been suggested for bifurcation treatment, including the use of conventional coronary stents, bioresorbable vascular scaffolds (BVS), drug-eluting balloons (DEB), and stents dedicated to bifurcations. This review will summarize different therapeutic approaches with their advantages and shortcomings, with special emphasis on histopathologic and physiologic effects of each treatment strategy.Expert opinion: Histopathology and clinical data have shown that a more simple treatment strategy is beneficial in bifurcation lesions, achieving superior results. Bifurcation interventions through balloon angioplasty or placement of stents can importantly alter the bifurcation's geometry and accordingly modify local flow conditions. Computational fluid dynamics (CFD) studies have shown that the outcome of bifurcation interventions is governed by local hemodynamic shear conditions. Minimizing detrimental flow conditions as much as possible should be the ultimate strategy to achieve long-term success of bifurcation interventions.

Local fluid dynamics in patients with bifurcated coronary lesions undergoing percutaneous coronary interventions

Although the coronary arteries are uniformly exposed to systemic cardiovascular risk factors, atherosclerosis development has a non-random distribution, which follows the local mechanical stresses including flow-related hemodynamic forces. Among these, wall shear stress plays an essential role and it represents the major flow-related factor affecting the distribution of atherosclerosis in coronary bifurcations. Furthermore, an emerging body of evidence suggests that hemodynamic factors such as low and oscillating wall shear stress may facilitate the development of in-stent restenosis and stent thrombosis after successful drug-eluting stent implantation. Drug-eluting stent implantation represents the gold standard for bifurcation interventions. In this specific setting of interventions on bifurcated lesions, the impact of fluid dynamics is expected to play a major role and constitutes substantial opportunity for future technical improvement. In the present review, available data is summarized regarding the role of local fluid dynamics in the clinical outcome of patients with bifurcated lesions.

The effects of stenting on hemorheological parameters: An in vitro investigation under various blood flow conditions

 Despite their wide clinical usage, stent functionality may be compromised by complications at the site of implantation, including early/late stent thrombosis and occlusion. Although several studies have described the effect of fluid-structure interaction on local haemodynamics, there is yet limited information on the effect of the stent presence on specific hemorheological parameters. The current work investigates the red blood cell (RBC) mechanical behavior and physiological changes as a result of flow through stented vessels. Blood samples from healthy volunteers were prepared as RBC suspensions in plasma and in phosphate buffer saline at 45% haematocrit. Self-expanding nitinol stents were inserted in clear perfluoroalkoxy alkane tubing which was connected to a syringe, and integrated in a syringe pump. The samples were tested at flow rates of 17.5, 35 and 70 ml/min, and control tests were performed in non-stented vessels. For each flow rate, the sample viscosity, RBC aggregation and deformability, and RBC lysis were estimated. The results indicate that the presence of a stent in a vessel has an influence on the hemorheological characteristics of blood. The viscosity of all samples increases slightly with the increase of the flow rate and exposure. RBC aggregation and elongation index (EI) decrease as the flow rate and exposure increases. RBC lysis for the extreme cases is evident. The results indicate that the stresses developed in the stent area for the extreme conditions could be sufficiently high to influence the integrity of the RBC membrane.

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.

Multiscale three-dimensional surface reconstruction and surface roughness of porcine left anterior descending coronary arteries

The aim of this study was to investigate the multiscale surface roughness characteristics of coronary arteries, to aid in the development of novel biomaterials and bioinspired medical devices. Porcine left anterior descending coronary arteries were dissected ex vivo, and specimens were chemically fixed and dehydrated for testing. Surface roughness was calculated from three-dimensional reconstructed surface images obtained by optical, scanning electron and atomic force microscopy, ranging in magnification from 10× to 5500×. Circumferential surface roughness decreased with magnification, and microscopy type was found to influence surface roughness values. Longitudinal surface roughness was not affected by magnification or microscopy types within the parameters of this study. This study found that coronary arteries exhibit multiscale characteristics. It also highlights the importance of ensuring consistent microscopy parameters to provide comparable surface roughness values.

Determination of hemodynamic risk for vascular disease in planar artery bifurcations

Understanding hemodynamics in blood circulation is crucial in order to unveil the mechanisms underlying the formation of stenosis and atherosclerosis. In fact, there are experimental evidences pointing out to the existence of some given vessel configurations that are more likely to develop the above mentioned pathologies. Along this manuscript, we performed an exhaustive investigation in a simplified model aiming to characterize by means of physical quantities those regions and configurations in vessel bifurcations that are more likely to develop such pathologies. The two-fold analysis is based, on the one hand, on numerical simulations (via CFD) and, on the other hand, on experiments realized in an ad-hoc designed polydimethylsiloxane (PDMS) channel with the appropriate parameters and appropriate fluid flows. The results obtained demonstrate that low velocity regions and low shear stress zones are located in the outer walls of bifurcations. In fact, we found that there is a critical range of bifurcation angles that is more likely to vascular disease than the others in correspondence with some experimental evidence. The effect of the inflow velocity on this critical range is also analyzed.

Effects of freezing, fixation and dehydration on surface roughness properties of porcine left anterior descending coronary arteries

BACKGROUND

To allow measurements of surface roughness to be made of coronary arteries using various imaging techniques, chemical processing, such as fixation and dehydration, is commonly used. Standard protocols suggest storing fresh biological tissue at -40°C. The aim of this study was to quantify the changes caused by freezing and chemical processing to the surface roughness measurements of coronary arteries, and to determine whether correction factors are needed for surface roughness measurements of coronary arteries following chemical processes typically used before imaging these arteries.

METHODS

Porcine left anterior descending coronary arteries were dissected ex vivo. Surface roughness was then calculated following three-dimensional reconstruction of surface images obtained using an optical microscope. Surface roughness was measured before and after a freeze cycle to assess changes during freezing, after chemical fixation, and again after dehydration, to determine changes during these steps of chemical processing.

RESULTS

No significant difference was caused due to the freeze cycle (p>0.05). There was no significant difference in the longitudinally measured surface roughness (Ra_L=0.99±0.39μm; p>0.05) of coronary arteries following fixation and dehydration either. However, the circumferentially measured surface roughness increased significantly following a combined method of processing (Ra_C=1.36±0.40, compared 1.98±0.27μm, respectively; p<0.05). A correction factor can compensate for the change Ra_Cβ=Ra_C1+0.46in Ra_C due to processing of tissue, Where Ra_Cβ, the corrected Ra_C, had a mean of 1.31±0.21μm.

CONCLUSIONS

Independently, freezing, fixation and dehydration do not alter the surface roughness of coronary arteries. Combined, however, fixation and dehydration significantly increase the circumferential, but not longitudinal, surface roughness of coronary arteries.