The impact of helical flow on coronary atherosclerotic plaque development

Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries

Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.

Aortic geometry and long-term outcome in patients with a repaired coarctation

OBJECTIVE

This study aims to compare aortic morphology between repaired coarctation patients and controls, and to identify aortic morphological risk factors for hypertension and cardiovascular events (CVEs) in coarctation patients.

METHODS

Repaired coarctation patients with computed tomography angiography (CTA) or magnetic resonance angiography (MRA) were included, followed-up and compared with sex-matched and age-matched controls. Three-dimensional aortic shape was reconstructed using patients' CTA or MRA, or four-dimensional flow cardiovascular magnetic resonance in controls, and advanced geometrical characteristics were calculated and visualised using statistical shape modelling. In patients, we examined the association of geometrical characteristics with (1) baseline hypertension, using multivariable logistic regression; and (2) cardiovascular events (CVE, composite of aortic complications, coronary artery disease, ventricular arrhythmias, heart failure hospitalisation, stroke, transient ischaemic attacks and cardiovascular death), using multivariable Cox regression. The least absolute shrinkage and selection operator (LASSO) method selected the most informative multivariable model.

RESULTS

Sixty-five repaired coarctation patients (23 years (IQR 19-38)) were included, of which 44 (68%) patients were hypertensive at baseline. After a median follow-up of 8.7 years (IQR 4.8-15.4), 27 CVEs occurred in 20 patients. Aortic arch dimensions were smaller in patients compared with controls (diameter p<0.001, wall surface area p=0.026, volume p=0.007). Patients had more aortic arch torsion (p<0.001) and a higher curvature (p<0.001). No geometrical characteristics were associated with hypertension. LASSO selected left ventricular mass, male sex, tortuosity and age for the multivariable model. Left ventricular mass (p=0.014) was independently associated with CVE, and aortic tortuosity showed a trend towards significance (p=0.070).

CONCLUSION

Repaired coarctation patients have a smaller aortic arch and a more tortuous course of the aorta compared with controls. Besides left ventricular mass index, geometrical features might be of importance in long-term risk assessment in coarctation patients.

Patient-Specific Numerical Simulations of Coronary Artery Hemodynamics and Biomechanics: A Pathway to Clinical Use

PURPOSE

Numerical models that simulate the behaviors of the coronary arteries have been greatly improved by the addition of fluid-structure interaction (FSI) methods. Although computationally demanding, FSI models account for the movement of the arterial wall and more adequately describe the biomechanical conditions at and within the arterial wall. This offers greater physiological relevance over Computational Fluid Dynamics (CFD) models, which assume the walls do not move or deform. Numerical simulations of patient-specific cases have been greatly bolstered by the use of imaging modalities such as Computed Tomography Angiography (CTA), Magnetic Resonance Imaging (MRI), Optical Coherence Tomography (OCT), and Intravascular Ultrasound (IVUS) to reconstruct accurate 2D and 3D representations of artery geometries. The goal of this study was to conduct a comprehensive review on CFD and FSI models on coronary arteries, and evaluate their translational potential.

METHODS

This paper reviewed recent work on patient-specific numerical simulations of coronary arteries that describe the biomechanical conditions associated with atherosclerosis using CFD and FSI models. Imaging modality for geometry collection and clinical applications were also discussed.

RESULTS

Numerical models using CFD and FSI approaches are commonly used to study biomechanics within the vasculature. At high temporal and spatial resolution (compared to most cardiac imaging modalities), these numerical models can generate large amount of biomechanics data.

CONCLUSIONS

Physiologically relevant FSI models can more accurately describe atherosclerosis pathogenesis, and help to translate biomechanical assessment to clinical evaluation.

Impact of Cross-Limb Stent-Graft Configuration on Hemodynamics in Abdominal Aortic Aneurysm Interventional Therapy

PURPOSE

The cross-limb (CL) technique is a commonly used endovascular treatment for addressing unfavorable anatomical features in patients with abdominal aortic aneurysm (AAA). The configuration of CL stent-graft plays a critical role in determining the postoperative hemodynamic properties and physiological behaviors, which ultimately impact the efficacy and safety of endovascular AAA treatment. This study aims to investigate the relationship between hemodynamics and CL stent-graft configuration from a hemodynamic perspective.

METHODS

Five distinct geometric models of cross-limb (CL) stent-graft configurations were constructed by optimizing the real clinical computed tomography angiography (CTA) data. These models varied in main body lengths and cross angles and were used to perform numerical simulations to analyze various hemodynamic parameters. Flow pattern, distribution of wall shear stress (WSS)-related parameters, localized normalized helicity (LNH), pressure drop, and the displacement force of all models were examined in this paper.

RESULTS

In patient-specific cases, helical flow and WSS increase with the main body. However, it also generated secondary flow in localized areas, leading to increased oscillation in the WSS direction. Notably, increasing the stent graft's main body length or decreasing the cross angle reduced the displacement force exerted on the stent-graft. Reducing the cross angle did not significantly alter the hemodynamic characteristics.

CONCLUSION

In the clinical practice of CL deployment, it is crucial to carefully consider the stent-graft configuration and the patient specific to achieve optimal postoperative outcomes. This study provides valuable insights for guiding stent selection and treatment planning in patients with abdominal aortic aneurysm undergoing CL techniques, from a hemodynamic perspective.

Analyzing the effects of helical flow in blood vessels using acoustofluidic-based dynamic flow generator

The effects of helical flow in a blood vessel are investigated in a dynamic flow generator using surface acoustic wave (SAW) in the microfluidic device. The SAW, generated by an interdigital transducer (IDT), induces acoustic streaming, resulting in a stable and consistent helical flow pattern in microscale channels. This approach allows rapid development of helical flow within the channel without directly contacting the medium. The precise design of the window enables the creation of distinct unidirectional vortices, which can be controlled by adjusting the amplitude of the SAW. Within this device, optimal operational parameters of the dynamic flow generator to preserve the integrity of endothelial cells are found, and in such settings, the actin filaments within the cells are aligned to the desired state. Our findings reveal that intracellular Ca2+ concentrations vary in response to flow conditions. Specifically, comparable maximum intensity and graphical patterns were observed between low-flow rate helical flow and high-flow rate Hagen-Poiseuille flow. These suggest that the cells respond to the helical flow through mechanosensitive ion channels. Finally, adherence of monocytes is effectively reduced under helical flow conditions in an inflammatory environment, highlighting the atheroprotective role of helical flow. STATEMENT OF SIGNIFICANCE: Helical flow in blood vessels is well known to prevent atherosclerosis. However, despite efforts to replicate helical flow in microscale channels, there is still a lack of in vitro models which can generate helical flow for analyzing its effects on the vascular system. In this study, we developed a method for generating steady and constant helical flow in microfluidic channel using acoustofluidic techniques. By utilizing this dynamic flow generator, we were able to observe the atheroprotective aspects of helical flow in vitro, including the enhancement of calcium ion flux and reduction of monocyte adhesion. This study paves the way for an in vitro model of dynamic cell culture and offers advanced investigation into helical flow in our circulatory system.

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.

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?".

Computational fluid dynamics as supporting technology for coronary artery disease diagnosis and treatment: an international survey

Background

Computational fluid dynamics (CFD) is emerging as an effective technology able to improve procedural outcomes and enhance clinical decision-making in patients with coronary artery disease (CAD). The present study aims to assess the state of knowledge, use and clinical acceptability of CFD in the diagnosis and treatment of CAD.

Methods

We realized a 20-questions international, anonymous, cross-sectional survey to cardiologists to test their knowledge and confidence on CFD as a technology applied to patients suffering from CAD. Responses were recorded between May 18, 2022, and June 12, 2022.

Results

A total of 466 interventional cardiologists (mean age 48.4 ± 8.3 years, males 362), from 42 different countries completed the survey, for a response rate of 45.9%. Of these, 66.6% declared to be familiar with the term CFD, especially for optimization of existing interventional techniques (16.1%) and assessment of hemodynamic quantities related with CAD (13.7%). About 30% of respondents correctly answered to the questions exploring their knowledge on the pathophysiological role of some CFD-derived quantities such as wall shear stress and helical flow in coronary arteries. Among respondents, 85.9% would consider patient-specific CFD-based analysis in daily interventional practice while 94.2% declared to be interested in receiving a brief foundation course on the basic CFD principles. Finally, 87.7% of respondents declared to be interested in a cath-lab software able to conduct affordable CFD-based analyses at the point-of-care.

Conclusions

Interventional cardiologists reported to be profoundly interested in adopting CFD simulations as a technology supporting decision making in the treatment of CAD in daily practice.

Vascular Cast to Program Antistenotic Hemodynamics and Remodeling of Vein Graft

The structural stability of medical devices is established by managing stress distribution in response to organ movement. Veins abruptly dilate upon arterial grafting due to the mismatched tissue property, resulting in flow disturbances and consequently stenosis. Vascular cast is designed to wrap the vein-artery grafts, thereby adjusting the diameter and property mismatches by relying on the elastic fixity. Here, a small bridge connection in the cast structure serves as an essential element to prevent stress concentrations due to the improved elastic fixity. Consequently, the vein dilation is efficiently suppressed, healthy (laminar and helical) flow is induced effectively, and the heathy functions of vein grafting are promoted, as indicated by the flow directional alignment of endothelial cells with arterialization, muscle expansion, and improved contractility. Finally, collaborative effects of the bridge drastically suppress stenosis with patency improvement. As a key technical point, the advantages of the bridge addition are validated via the computational modeling of fluid-structure interaction, followed by a customized ex vivo set-up and analyses. The calculated effects are verified using a series of cell, rat, and canine models towards translation. The bridge acted like "Little Dutch boy" who saved the big mass using one finger by supporting the cast function.

A review on non-Newtonian fluid models for multi-layered blood rheology in constricted arteries

Haemodynamics is a branch of fluid mechanics which investigates the features of blood when it flows not only via blood vessels of smaller/larger diameter, but also under normal as well as abnormal flow states, such as in the presence of stenosis, aneurysm, and thrombosis. This review aims to discuss the rheological properties of blood, geometry of constrictions, dilations and the emergence of single-layered fluid to four-layered fluid models. To discuss further the influence of the aforesaid parameters on the physiologically important flow quantities, the mathematical formulation and solution methodology of the two-layered and four layered arterial blood flow problems studied by the authors (Afiqah and Sankar in ARPN J Eng Appl Sci 15:1129--1143, 2020, Comput Methods Programs Biomed 199:105907, 2021. 10.1016/j.cmpb.2020.105907) are recalled. It should be pointed out that the increasing resistive impedance to flow in three distinct states encompassing healthy, anaemic, and diabetic demonstrates that the greater the restriction in the artery, very few blood is carried to the pathetic organs, leading to subjects' death. It is also discovered that the pulsatile nature of blood movement produces a dynamic environment that poses a slew of intriguing and unstable fluid mechanical state. It is hoped that the intriguing results gathered from this literature survey and review conducted may help the medical practitioners to forecast blood behaviour mobility in stenotic arteries. Furthermore, the physiological information gathered from the available clinical data from the literature on patients diagnosed with diabetes and anaemia may be beneficial to doctors in deciding the therapeutic procedure for treating some particular cardiovascular disease.

Computational Study of Hemodynamic Field of an Occluded Artery Model with Anastomosis

In this research work, the hemodynamic field of an occluded artery with anastomosis by means of computational simulation has been studied. The main objective of the current study is the investigation of 3D flow field phenomena in the by-pass region and the effect of the bypass graft to stenosis volume flow ratio on their formation. The anastomosis type was end-to-side with a 45° angle, while stenosis imposed a 75% area blockage of the aorta vessel and the total volume flow was 220 lt/h. The computational study of the flow field was utilized via a laminar flow model and three turbulence models (k-ε RNG, standard k-ω, and k-ω SST). Numerical results were compared qualitatively with experimental visualizations carried out under four different flow conditions, varying according to the flow ratio between the stenosis and the anastomotic graft. Comparison between computational results and experimental visualization findings exhibited a good agreement. Results showed that SST k-ω turbulence models reproduce better visually obtained flow patterns. Furthermore, cross-sectional velocity distributions demonstrated two distinct flow patterns down the bypass graft, depending on the flow ratio. Low values of flow ratio are characterized by fluid rolling up, whereas for high values fluid volume twisting was observed. Finally, areas with low wall shear stresses were mapped, as these are more prone to postoperative degradation of the bypass graft due to the development of subendothelial hyperplasia.

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.

Spatiotemporal changes of local hemodynamics and plaque components during atherosclerotic progression in rabbit

BACKGROUND AND OBJECTIVE

Recent evidence demonstrates that the atherogenic process is discontinuous. Our goal is to study changes of plaque components and local hemodynamics during atherosclerotic progression.

METHODS

The histological and immunohistochemical staining of high-fat diet rabbit aorta were evaluated at 0, 8, 10 and 12 weeks, respectively. In addition, the blood flow and LDL transport were simulated at the above four time points.

RESULTS

The plaque thickness at different characteristic regions increased at different rates. The collagen continued to increase, while the elastin, fibronectin, macrophages and smooth muscle cells increased first and then decreased. The relative surface LDL concentration decreased at 8 weeks, and then it increased first and decreased slightly. Meanwhile, the hemodynamic environment became better firstly at 8 weeks, then got slightly worse and lastly improved again.

CONCLUSIONS

The local hemodynamics and plaque components vary nonlinearly during atherosclerotic progression in rabbit aorta.

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 definition of low wall shear stress and its effect on plaque progression estimation in human coronary arteries

Wall shear stress (WSS), the frictional force of the blood on the vessel wall, plays a crucial role in atherosclerotic plaque development. Low WSS has been associated with plaque growth, however previous research used different approaches to define low WSS to investigate its effect on plaque progression. In this study, we used four methodologies to allocate low, mid and high WSS in one dataset of human coronary arteries and investigated the predictive power of low WSS for plaque progression. Coronary reconstructions were based on multimodality imaging, using intravascular ultrasound and CT-imaging. Vessel-specific flow was measured using Doppler wire and computational fluid dynamics was performed to calculate WSS. The absolute WSS range varied greatly between the coronary arteries. On the population level, the established pattern of most plaque progression at low WSS was apparent in all methodologies defining the WSS categories. However, for the individual patient, when using measured flow to determine WSS, the absolute WSS values range so widely, that the use of absolute thresholds to determine low WSS was not appropriate to identify regions at high risk for plaque progression.

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.

In Vivo Based Fluid-Structure Interaction Biomechanics of the Left Anterior Descending Coronary Artery

A fluid-structure interaction-based biomechanical model of the entire left anterior descending coronary artery is developed from in vivo imaging via the finite element method in this paper. Included in this investigation is ventricle contraction, three-dimensional motion, all angiographically visible side branches, hyper/viscoelastic artery layers, non-Newtonian and pulsatile blood flow, and the out-of-phase nature of blood velocity and pressure. The fluid-structure interaction model is based on in vivo angiography of an elite athlete's entire left anterior descending coronary artery where the influence of including all alternating side branches and the dynamical contraction of the ventricle is investigated for the first time. Results show the omission of side branches result in a 350% increase in peak wall shear stress and a 54% decrease in von Mises stress. Peak von Mises stress is underestimated by up to 80% when excluding ventricle contraction and further alterations in oscillatory shear indices are seen, which provide an indication of flow reversal and has been linked to atherosclerosis localization. Animations of key results are also provided within a video abstract. We anticipate that this model and results can be used as a basis for our understanding of the interaction between coronary and myocardium biomechanics. It is hoped that further investigations could include the passive and active components of the myocardium to further replicate in vivo mechanics and lead to an understanding of the influence of cardiac abnormalities, such as arrythmia, on coronary biomechanical responses.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Mechanosensitive Piezo1 Channel Evoked-Mechanical Signals in Atherosclerosis

Recently, more and more works have focused and used extensive resources on atherosclerosis research, which is one of the major causes of death globally. Alongside traditional risk factors, such as hyperlipidemia, smoking, hypertension, obesity, and diabetes, mechanical forces, including shear stress, pressure and stretches exerted on endothelial cells by flow, is proved to be crucial in atherosclerosis development. Studies have recognized the mechanosensitive Piezo1 channel as a special sensor and transducer of various mechanical forces into biochemical signals, and recent studies report its role in atherosclerosis through different mechanical forces in pressure, stretching and turbulent shear stress. Based on our expertise in this field and considering the recent advancement of atherosclerosis research, we will be focusing on the function of Piezo1 and its involvement in various cellular mechanisms and consequent involvement in the development of atherosclerosis in this review. Also, we will discuss various functions of Piezo1 involvement in atherosclerosis and come up with new mechanistic insight for future research. Based on the recent findings, we suggest Piezo1 as a valid candidate for novel therapeutic innovations, in which deep exploration and translating its findings into the clinic will be a new therapeutic strategy for cardiovascular diseases, particularly atherosclerosis.

Contribution of atherosclerotic plaque location and severity to the near-wall hemodynamics of the carotid bifurcation: an experimental study and FSI modeling

Atherosclerosis is initiated by endothelial injury that is related to abnormal values of hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI) and stress phase angle (SPA), which are more common in arterial bifurcations due to the complex structure. An experimental model of human carotid bifurcation with accurate geometrical and mechanical features was set up, and using realistic pulsatile flow rates, the inlet and outlet pressure pulses were measured for normal and stenosed models with 40% and 80% severities at common carotid (CCA), internal carotid (ICA) and external carotid (ECA) arteries. Based on the obtained experimental data, fluid-structure models were developed to obtain WSS, OSI, and SPA and evaluate pathological consequences at different locations. Mild severity had minor impact, however, inducing severe 80% stenosis in each branch led to considerable localized changes of hemodynamic parameters both in the stenosis site and other locations. This included sharp increases in WSS values accompanied by very low values close to zero before and after the peaks. Severe stenosis not only caused significant changes in the local artery, but also in other branches. OSI and SPA were less sensitive to stenosis, although high peaks were observed on bifurcation site for the stenosis at ECA. The interconnection of arteries at carotid bifurcation results in altered pressure/flow patterns in all branches when a stenosis is applied in any site. Such effect confirms pathological findings that atherosclerotic plaques are observed simultaneously in different carotid branches, although with different degrees of plaque growth and severity.

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.

Quantitative assessment of coronary plaque volume change related to triglyceride glucose index: The Progression of AtheRosclerotic PlAque DetermIned by Computed TomoGraphic Angiography IMaging (PARADIGM) registry

Background

The association between triglyceride glucose (TyG) index and coronary atherosclerotic change remains unclear. We aimed to evaluate the association between TyG index and coronary plaque progression (PP) using serial coronary computed tomography angiography (CCTA).

Methods

A total of 1143 subjects (aged 60.7 ± 9.3 years, 54.6% male) who underwent serial CCTA with available data on TyG index and diabetic status were analyzed from The Progression of AtheRosclerotic PlAque DetermIned by Computed TomoGraphic Angiography IMaging (PARADIGM) registry. PP was defined as plaque volume (PV) (mm³) at follow-up minus PV at index > 0. Annual change of PV (mm³/year) was defined as PV change divided by inter-scan period. Rapid PP was defined as the progression of percent atheroma volume (PV divided by vessel volume multiplied by 100) ≥ 1.0%/year.

Results

The median inter-scan period was 3.2 (range 2.6–4.4) years. All participants were stratified into three groups based on TyG index tertiles. The overall incidence of PP was 77.3%. Baseline total PV (group I [lowest]: 30.8 (0.0–117.7), group II: 47.2 (6.2–160.4), and group III [highest]: 57.5 (8.4–154.3); P < 0.001) and the annual change of total PV (group I: 5.7 (0.0–20.2), group II: 7.6 (0.5–23.5), and group III: 9.4 (1.4–27.7); P = 0.010) were different among all groups. The risk of PP (odds ratio [OR] 1.648; 95% confidence interval [CI] 1.167–2.327; P = 0.005) and rapid PP (OR 1.777; 95% CI 1.288–2.451; P < 0.001) was increased in group III compared to that in group I. TyG index had a positive and significant association with an increased risk of PP and rapid PP after adjusting for confounding factors.

Conclusion

TyG index is an independent predictive marker for the progression of coronary atherosclerosis.

Clinical registration ClinicalTrials.gov NCT02803411

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.