Effect of geometrical assumptions on numerical modeling of coronary blood flow under normal and disease conditions

Comparison of stable carotid plaques in patients with mild-to-moderate carotid stenosis with vulnerable plaques in patients with significant carotid stenosis

To compare the characteristics of stable and vulnerable carotid plaques, and investigate the diagnostic performance of wall shear stress (WSS) based on magnetic resonance plaque imaging in carotid plaques. Retrospectively analyzed and divided 64 atherosclerotic plaques into stable carotid plaque groups with mild-to-moderate stenosis and vulnerable carotid plaque groups with significant stenosis. Computational fluid dynamics simulations were performed to calculate WSS parameters by using three-dimensional wall geometry based on high-resolution magnetic resonance plaque imaging of carotid bifurcation and patient specific boundary conditions obtained through color Doppler ultrasound. WSS parameters including upstream (WSSup), downstream (WSSdown), and core (WSScore) of plaque. The WSS parameters values were compared between the stable and vulnerable carotid plaque groups. Receiver operating characteristic curves and area under the curve (ROC-AUC) and Python were used to evaluate discriminative efficacy of WSS. WSSdown exhibited significant decrease in the vulnerable carotid plaque group (2.88 ± 0.41 Pa) compared to the stable carotid plaque group (4.47 ± 0.84 Pa) (P = .003). The difference of WSSup (3.28 ± 0.85 Pa vs 4.02 ± 0.74 Pa) and WSScore (1.12 ± 0.18 Pa vs 1.38 ± 0.38 Pa) between the two groups were also pronounced (P = .02, 0.01, respectively). The ROC-AUC values for WSSup, WSSdown, WSScore were 0.75 (95% CI, 0.58-0.93), 0.96 (95% CI, 0.79-1.14), 0.69 (95% CI, 0.56-0.83) respectively. When the value of WSSdown was 3.5 Pa, the sensitivity was 93.7% (95% CI, 76.1-111), specificity and accuracy was 87.5% (95% CI, 70.0-105), 88.4% (95% CI, 70.6-105) respectively. Notably, among these parameters, WSSdown demonstrated the highest discriminative efficiency with a F1 Score of 0.90, Diagnostic Odds Ratio of 105.0 and Matthews Correlation Coefficient of 0.81. Vulnerable carotid plaques with significant stenosis have lower WSS compared to stable plaques with mild-to-moderate stenosis, and downstream WSS showing the highest diagnostic efficacy.

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.

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.

On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis

BACKGROUND AND OBJECTIVE

In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied.

METHODS

We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology.

RESULTS

22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately.

CONCLUSIONS

These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.

Hemodynamic and fluid flow analysis of a cerebral aneurysm: a CFD simulation

In this study, we investigate the hemodynamics parameters and their impact on the aneurysm rupture. The simulations are performed on an ideal (benchmark) and realistic model for the intracranial aneurysm that appears at the anterior communicating artery. The realistic geometry was reconstructed from patient-specific cerebral arteries. The computational fluid dynamics simulations are utilized to investigate the hemodynamic parameters such as flow recirculation, wall shear stress, and wall pressure. The boundary conditions are measured from the patient using ultrasonography. The solution of the governing equations is obtained by using the ANSYS-FLUENT 19.2 package. The CFD results indicate that the flow recirculation appears in the aneurysms zone. The effect of the flow recirculation on the bulge hemodynamics wall parameters is discussed to identify the rupture zone.

Segmenting 3D geometry of left coronary artery from coronary CT angiography using deep learning for hemodynamic evaluation

While coronary CT angiography (CCTA) is crucial for detecting several coronary artery diseases, it fails to provide essential hemodynamic parameters for early detection and treatment. These parameters can be easily obtained by performing computational fluid dynamic (CFD) analysis on the 3D artery geometry generated by CCTA image segmentation. As the coronary artery is small in size, manually segmenting the left coronary artery from CCTA scans is a laborious, time-intensive, error-prone, and complicated task which also requires a high level of expertise. Academics recently proposed various automated segmentation techniques for combatting these issues. To further aid in this process, we present CoronarySegNet, a deep learning-based framework, for autonomous and accurate segmentation as well as generation of 3D geometry of the left coronary artery. The design is based on the original U-net topology and includes channel-aware attention blocks as well as deep residual blocks with spatial dropout that contribute to feature map independence by eliminating 2D feature maps rather than individual components. We trained, tested, and statistically evaluated our model using CCTA images acquired from various medical centers across Bangladesh and the Rotterdam Coronary Artery Algorithm Evaluation challenge dataset to improve generality. In empirical assessment, CoronarySegNet outperforms several other cutting-edge segmentation algorithms, attaining dice similarity coefficient of 0.78 on an average while being highly significant (p < 0.05). Additionally, both the 3D geometries generated by machine learning and semi-automatic method were statistically similar. Moreover, hemodynamic evaluation performed on these 3D geometries showed comparable results.

Blood flow CFD simulation on a cerebral artery of a stroke patient

Abstract

The purpose of this paper is to conduct a numerical simulation of the stroke patient's cerebral arteries and investigate the flow parameters due to the presence of stenosis. The computational fluid dynamics (CFD) simulations are based on simplified and realistic cerebral artery models. The seven simplified models (benchmarks) include straight cylindrical vessels with idealized stenosis with variable d/D (0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1). The realistic model of the cerebral artery is based on magnetic resonance imaging (MRI) for patient-specific cerebral arteries. The simulation for the realistic model of the cerebral artery is performed at boundary conditions measured by ultrasonography of the input and the output flow profiles (velocity and pressure). The obtained CFD results of the benchmarks are validated with actual data from the literature. Furthermore, a previous vascular contraction is assumed to be exist and the effect of this contraction area ratio on the blood flow regime is discussed and highlighted. Furthermore, CFD results show that a certain vascular contraction area critically affects the blood flow which shows increasing the wall shear stress WSS at the stenosis site. An increase in the blood velocity and vortex appears after the contraction zone, this lead to vessel occlusion and strokes.

Article highlights

Numerical simulation of pulsatile blood flow characteristics in a multi stenosed coronary artery

BACKGROUND

Coronary artery disease is reported as one of the most common sources of death all over the world. The presence of stenosis (plaque) in the coronary arteries results in the restriction of blood supply, which leads to myocardial infarction.

OBJECTIVE

The aim of this study was to investigate the effect of multi stenosis on hemodynamics parameters in idealized coronary artery models with varying degrees of stenosis and interspace distance between the stenosis.

METHODS

A finite volume-based software package (Ansys CFX version 17.2) was employed to model the blood flow. The hemodynamic stenosis parameters of blood, such as the pressure, velocity, and wall shear stress were obtained.

RESULTS

The computed results showed that the pressure drop is maximum across the 90% area stenosis (AS). The pressure drop is increased as the distance between the proximal and distal stenosis is decreased across the proximal stenosis for the model P70_D70 during the systolic period of the cardiac cycle. A recirculation zone is formed behind the stenosis and is restricted by the occurrence of distal stenosis as the interspacing distance decreases, which could lead to further progression of stenosis in the flow-disturbed area. The wall shear stress was found to increase as the distance between the proximal and distal stenosis is increased across the distal stenosis. The maximum wall shear stress was found at 90% AS.

CONCLUSIONS

In the clinical diagnosis, an overestimation of distal stenosis severity could be possible. Furthermore, the low wall shear stress zone in between the proximal and distal stenosis may help atherosclerotic growth or merge adjacent stenosis.

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.

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.

Platelet adhesion potential estimation in a normal and diseased coronary artery model: effects of shear stress magnitude versus shear stress history

Platelets play a salient role in the pathogenesis of coronary diseases; primarily through their adhesion to other platelets, endothelial cells and plasma proteins. It is necessary for platelets to activate in order for them to adhere to these different substrates. One of the key regulatory mechanical factors in platelet activation is shear stress, which has been shown to alter multiple platelet functions through the activation of mechanoreceptors. Our goal was to investigate how different numerical shear stress tracking techniques affect platelet adhesion estimates within physiologically relevant computational models. Previously, we developed a physiological coronary artery computational fluid dynamics model. Shear stress waveforms, obtained from these models, were used to monitor in vitro platelet and endothelial cell adhesion marker expression. In this work, the adhesion marker expression data was regressed to obtain numerical functions for receptor expression predictions. These functions were input into a customized adhesion model utilizing different shear stress tracking techniques. For the normal vascular conditions and minimal pathological disease models, shear stress tracking did not significantly affect the adhesion estimates. However, for the severe pathological model, the two shear stress tracking methods had vastly different estimates. Therefore, shear stress tracking methods must be chosen accurately to predict platelet adhesion potentials for accurate modeling techniques.

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.

Influence of bifurcation angle in left coronary artery with stenosis: A CFD analysis

BACKGROUND: The left coronary artery commonly known as LCA gets divided into two branches, such as the left circumflex (LCX) and left anterior descending (LAD) at a particular angle. This angle is varies from person to person. The present computational study contributes remarkable expertise about the influence of this angle variation on the hemodynamic parameters in the presence of 80% area stenosis at the LAD branch. OBJECTIVE: This study aimed to compare the effect of the bifurcation angle on hemodynamic parameters in the left coronary artery with 80% stenosis. METHOD: Computational models of left coronary bifurcation angles of 30°, 60°, 90°, 120° were developed to understand the flow behavior of left coronary artery branches. The 80% area stenosis (AS) is considered at the LAD branch immediate to bifurcation. RESULTS: Measurements of pressure, velocity and wall shear stress were carried out corresponding to various bifurcation angles. It was found that the drop-in pressure increases as the angle increases from narrow to wider. A slight elevation in the velocity at the stenosis was observed. In addition, the obtained results further reveal a recirculation region immediately after the plaque, which leads to more deposition of plaque in the flow obstructed area. It is known that the shear stress at the arterial wall across the stenosis increases as the angle of bifurcation increases from narrow to wider. CONCLUSIONS: The bifurcation of the left coronary artery and size of the stenosis have a notable impact on the pressure and wall shear stress. These two factors should be given due consideration by cardiologists to assess the complexity of stenosis in the LCA branches.

Numerical computation of blood hemodynamic through constricted human left coronary artery: Pulsatile simulations

BACKGROUND AND OBJECTIVE

The accumulation of plaque in the coronary artery of the human heart restricts the path of blood flow in that region and leads to Coronary Artery Disease. This study's goal is to present the pulsatile blood flow conduct through four different levels of constrictions, i.e., healthy, 25%, 50%, and 75% in human left coronary arteries.

METHODS

Using CT scan data of a healthy person, the two-dimensional coronary model is constructed. A non-Newtonian Carreau model is used to study the maximum flow velocity, streamline effect, and maximum Wall Shear Stress at the respective constricted areas over the entire cardiac cycle. Finite Volume Method is executed for solving the governing equations. The fluctuating Wall Shear Stress (WSS) at different levels was assessed using Computational Fluid Dynamics (CFD).

RESULTS

The comparative study of the diseased arteries showcases that at the systolic phase, the 75% blocked artery attains the maximum velocity of 0.14 m/s and 0.53 m/s at t=0.005 s and t=0.115 s, respectively. While the maximum velocity takes a significant drop at t=0.23 s and t=0.345 s, this marks the diastolic phase. The streamline contour showcased the blood flow conduct at different phases of the cardiac cycle. At the peak systolic phase, a dense flow separation was observed near the blocked regions. It highlights the disturbed flow in that particular region. The most severely diseased artery acquires the maximum WSS of 18.81 Pa at the peak systolic phase, i.e., at t=0.115 s.

CONCLUSIONS

The computational study of the hemodynamic parameters can aid in the early anticipation of the degree of the severity of the diseased arteries. This study, in a way, could benefit doctors/surgeons to plan an early treatment/surgery on the grounds of the severity of the disease. Thus, a before time prognosis could restrain the number of deaths caused due to Coronary Artery Disease.

A comprehensive fluid-structure interaction model of the left coronary artery

A fluid structure interaction model of a left anterior descending (LAD) coronary artery was developed, incorporating transient blood flow, cyclic bending motion of the artery, and myocardial contraction. The 3D geometry was constructed based on a patient's computed tomography angiography data. To simulate disease conditions, a plaque was placed within the LAD to create a 70% stenosis. The bending motion of the blood vessel was prescribed based on the LAD spatial information. The pressure induced by myocardial contraction was applied to the outside of the blood vessel wall. The fluid domain was solved using the Navier-Stokes equations. The arterial wall was defined as a nonlinear elastic, anisotropic, and incompressible material, and the mechanical behavior was described using the modified hyper-elastic Mooney-Rivlin model. The fluid (blood) and solid (vascular wall) domains were fully coupled. The simulation results demonstrated that besides vessel bending/stretching motion, myocardial contraction had a significant effect on local hemodynamics and vascular all stress/strain distribution. It not only transiently increased blood flow velocity and fluid wall shear stress, but also changed shear stress patterns. The presence of the plaque significantly reduced vascular wall tensile strain. Compared to the coronary artery models developed previously, the current model had improved physiological relevance.

Towards Additive Manufacture of Functional, Spline-Based Morphometric Models of Healthy and Diseased Coronary Arteries: In Vitro Proof-of-Concept Using a Porcine Template

The aim of this study is to assess the additive manufacture of morphometric models of healthy and diseased coronary arteries. Using a dissected porcine coronary artery, a model was developed with the use of computer aided engineering, with splines used to design arteries in health and disease. The model was altered to demonstrate four cases of stenosis displaying varying severity, based on published morphometric data available. Both an Objet Eden 250 printer and a Solidscape 3Z Pro printer were used in this analysis. A wax printed model was set into a flexible thermoplastic and was valuable for experimental testing with helical flow patterns observed in healthy models, dominating the distal LAD (left anterior descending) and left circumflex arteries. Recirculation zones were detected in all models, but were visibly larger in the stenosed cases. Resin models provide useful analytical tools for understanding the spatial relationships of blood vessels, and could be applied to preoperative planning techniques, but were not suitable for physical testing. In conclusion, it is feasible to develop blood vessel models enabling experimental work; further, through additive manufacture of bio-compatible materials, there is the possibility of manufacturing customized replacement arteries.

Fiber-reinforced computational model of the aortic root incorporating thoracic aorta and coronary structures

Cardiovascular diseases are still the leading causes of death in the developed world. The decline in the mortality associated with circulatory system diseases is accredited to development of new diagnostic and prognostic tools. It is well known that there is an inter relationship between the aortic valve impairment and pathologies of the aorta and coronary vessels. However, due to the limitations of the current tools, the possible link is not fully elucidated. Following our previous model of the aortic root including the coronaries, in this study, we have further developed the global aspect of the model by incorporating the anatomical structure of the thoracic aorta. This model is different from all the previous studies in the sense that inclusion of the coronary structures and thoracic aorta into the natural aortic valve introduces the notion of globality into the model enabling us to explore the possible link between the regional pathologies. The developed model was first validated using the available data in the literature under physiological conditions. Then, to provide a support for the possible association between the localized cardiovascular pathologies and global variations in hemodynamic conditions, we simulated the model for two pathological conditions including moderate and severe aortic valve stenoses. The findings revealed that malformations of the aortic valve are associated with development of low wall shear stress regions and helical blood flow in thoracic aorta that are considered major contributors to aortic pathologies.

3D physiological model of the aortic valve incorporating small coronary arteries

The diseases of the coronary arteries and the aortic root are still the leading causes of mortality and morbidity worldwide. In this study, a 3D global fluid-structure interaction of the aortic root with inclusion of anatomically inspired small coronary arteries using the finite element method is presented. This innovative model allows to study the impact and interaction of root biomechanics on coronary hemodynamics and brings a new understanding to small coronary vessels hemodynamics. For the first time, the velocity profiles and shear stresses are reported in distal coronary arteries as a result of the aortic flow conditions in a global fluid-structure interaction model.

The impact of the aortic valve impairment on the distant coronary arteries hemodynamics: a fluid-structure interaction study

Atherosclerosis is still the leading cause of death in the developed world. Although its initiation and progression is a complex multifactorial process, it is well known that blood flow-induced wall shear stress (WSS) is an important factor involved in early atherosclerotic plaque initiation. In recent clinical studies, it was established that the regional pathologies of the aortic valve can be involved in the formation of atherosclerotic plaques. However, the impact of hemodynamic effects is not yet fully elucidated for disease initiation and progression. In this study, our developed 3D global fluid-structure interaction model of the aortic root incorporating coronary arteries is used to investigate the possible interaction between coronary arteries and aortic valve pathologies. The coronary hemodynamics was examined and quantified for different degrees of aortic stenosis varying from nonexistent to severe. For the simulated healthy model, the calculated WSS varied between 0.41 and 1.34 Pa which is in the atheroprotective range. However, for moderate and severe aortic stenoses, wide regions of the coronary structures, especially the proximal sections around the first bifurcation, were exposed to lower values of WSS and therefore they were prone to atherosclerosis even in the case of healthy coronary arteries.

A NEW METHOD TO STUDY ATHEROSCLEROSIS AND ASSESS THE EFFECTIVENESS OF PERCUTANEOUS CORONARY INTERVENTION BY COMPUTATIONAL FLUID DYNAMICS SIMULATION

High-pitch spiral computed tomography coronary angiography (CTCA) is able to perform a whole-heart scan within one heartbeat, resulting in high-quality images with high spatial and temporal resolution. To investigate the performance of high-quality CTCA images, an anatomic stenosis evaluation by digital subtracted angiography (DSA) was compared to a functional stenosis evaluation by CTCA-derived fraction flow reserve (FFR). A total of 54 arterial segments with stenosis were collected from 23 patients, and three-dimensional (3D) geometrical models were reconstructed. The computational fluid dynamics (CFDs) analysis was used to calculate the pressure distributions and FFR values. The correlation between anatomic and functional evaluation factors was assessed with either the ratio of anatomic reduction or CTCA-derived FFR values at the corresponding anatomic locations. Pearson correlation analysis was performed, and a significant correlation was found relating to the diameter ([Formula: see text]) and the cross-sectional area ([Formula: see text]). A significant correlation was also found in the functional evaluation relating to the diameter ([Formula: see text]) and the cross-sectional area ([Formula: see text]). High-quality CT images greatly reduce the time needed for geometric reconstruction. Significant advances in the accuracy of the reconstruction have resulted in more accurate CFD analysis, which can help to improve clinical diagnoses. The results of this study show that the CFD method can be a feasible tool for the clinic diagnosis of stenosis and for determining whether a patient requires percutaneous coronary intervention (PCI).

Platelets modulate endothelial cell response to dynamic shear stress through PECAM-1

INTRODUCTION

Both vascular endothelial cells and platelets are sensitive to blood flow induced shear stress. We have recently reported that platelet-endothelial cell interaction could greatly affect platelet activation under flow. In the present study, we aimed to investigate how platelet-endothelial cell interaction affected endothelial cell inflammatory responses under flow.

MATERIALS AND METHODS

Human coronary artery endothelial cells were exposed to normal or low pulsatile shear stress with or without the presence of platelets. Following shear exposure, endothelial cell ICAM-1 expression was measured using ELISA, Western blot and PCR; cell surface PECAM-1 expression/phosphorylation was measured using ELISA. Platelet adhesion to endothelial cells was quantified using immunofluorescence microscopy. To determine the role of PECAM-1 in platelet-endothelial cell interaction, endothelial cell PECAM-1 expression was suppressed using siRNA.

RESULTS

Pathological low shear stress induced a significant increase in endothelial cell ICAM-1 expression, at both protein and mRNA levels. Platelet adhesion to endothelial cells increased significantly under low shear stress, co-localizing with PECAM-1 at endothelial cell junctions. The presence of platelets inhibited low shear stress-induced ICAM-1 upregulation. When endothelial cell PECAM-1 expression was suppressed, platelet adhesion to endothelial cells under low shear stress decreased significantly; endothelial cell ICAM-1 expression was not affected by shear stress, with or without platelets.

CONCLUSIONS

These results suggested that PECAM-1 could mediate platelet adhesion to endothelial cells under shear stress. Platelets binding to endothelial cells interfered with endothelial cell mechanotransduction through PECAM-1, affecting endothelial cell inflammatory responses towards pathological shear flow.

Patient-Specific Computational Models of Coronary Arteries Using Monoplane X-Ray Angiograms

Coronary artery disease (CAD) is the most common type of heart disease in western countries. Early detection and diagnosis of CAD is quintessential to preventing mortality and subsequent complications. We believe hemodynamic data derived from patient-specific computational models could facilitate more accurate prediction of the risk of atherosclerosis. We introduce a semiautomated method to build 3D patient-specific coronary vessel models from 2D monoplane angiogram images. The main contribution of the method is a robust segmentation approach using dynamic programming combined with iterative 3D reconstruction to build 3D mesh models of the coronary vessels. Results indicate the accuracy and robustness of the proposed pipeline. In conclusion, patient-specific modelling of coronary vessels is of vital importance for developing accurate computational flow models and studying the hemodynamic effects of the presence of plaques on the arterial walls, resulting in lumen stenoses, as well as variations in the angulations of the coronary arteries.

Tobacco and e-cigarette products initiate Kupffer cell inflammatory responses

Kupffer cells are liver resident macrophages that are responsible for screening and clearing blood of pathogens and foreign particles. It has recently been shown that Kupffer cells interact with platelets, through an adhesion based mechanism, to aid in pathogen clearance and then these platelets re-enter the general systemic circulation. Thus, a mechanism has been identified that relates liver inflammation to possible changes in the systemic circulation. However, the role that Kupffer cells play in cardiovascular disease initiation/progression has not been elucidated. Thus, our objective was to determine whether or not Kupffer cells are responsive to a classical cardiovascular risk factor and if these changes can be transmitted into the general systemic circulation. If Kupffer cells initiate inflammatory responses after exposure to classical cardiovascular risk factors, then this provides a potential alternative/synergistic pathway for cardiovascular disease initiation. We aimed to elucidate the prevalence of this potential pathway. We hypothesized that Kupffer cells would initiate a robust inflammatory response after exposure to tobacco cigarette or e-cigarette products and that the inflammatory response would have the potential to antagonize other salient cells for cardiovascular disease progression. To test this, Kupffer cells were incubated with tobacco smoke extracts, e-cigarette vapor extracts or pure nicotine. Complement deposition onto Kupffer cells, Kupffer cell complement receptor expression, oxidative stress production, cytokine release and viability and density were assessed after the exposure. We observed a robust inflammatory response, oxidative stress production and cytokine release after Kupffer cells were exposed to tobacco or e-cigarette extracts. We also observed a marginal decrease in cell viability coupled with a significant decrease in cell density. In general, this was not a function of the extract formulation (e.g. tobacco vs. e-cigarette products or the formulation of the cigarette product). These results indicate that Kupffer cells are responsive to classical cardiovascular risk factors and that an inflammatory response is initiated that may pass into the general systemic circulation.

Numerical investigation of blood flow in three-dimensional porcine left anterior descending artery with various stenoses

Coronary heart disease causes obstruction of coronary blood flow and is the leading cause of death worldwide. The effect of focal stenosis on downstream flow pattern in the coronary arterial tree is not well understood. Here, the blood flows in normal and diseased porcine left anterior descending (LAD) arterial tree were modeled and compared to determine the effects of stenosis on the blood flow distribution and hemodynamic parameters. The anatomical model of LAD was extracted from a porcine heart by computed tomography (CT), which was comprised of a main trunk and nine side branches. Stenoses with various severities were imposed into the main trunk between the first and second side branches, and the boundary condition at each outlet accounted for the effect of stenosis on the flow rate in the downstream vasculature. It was found that only significant stenosis (≥75% area reduction) considerably altered pressure drop and total flow rate distribution in branches and at each bifurcation. The effect of significant stenosis on bifurcations, however, diminished at downstream locations. As demonstrated by distributions of oscillatory shear index and relative residence time, non-significant stenosis (<75% area reduction) has the potential to induce atherosclerosis near the ostium of downstream side branch, while significant stenosis can promote atherosclerosis in its wake.

Coronary CT angiography: Beyond morphological stenosis analysis

Rapid technological developments in computed tomography (CT) imaging technique have made coronary CT angiography an attractive imaging tool in the detection of coronary artery disease. Despite visualization of excellent anatomical details of the coronary lumen changes, coronary CT angiography does not provide hemodynamic changes caused by presence of plaques. Computational fluid dynamics (CFD) is a widely used method in the mechanical engineering field to solve complex problems through analysing fluid flow, heat transfer and associated phenomena by using computer simulations. In recent years, CFD is increasingly used in biomedical research due to high performance hardware and software. CFD techniques have been used to study cardiovascular hemodynamics through simulation tools to assist in predicting the behaviour of circulatory blood flow inside the human body. Blood flow plays a key role in the localization and progression of coronary artery disease. CFD simulation based on 3D luminal reconstructions can be used to analyse the local flow fields and flow profiling due to changes of vascular geometry, thus, identifying risk factors for development of coronary artery disease. The purpose of this article is to provide an overview of the coronary CT-derived CFD applications in coronary artery disease.

Effects of stenosis on the porcine left anterior descending arterial tree

In this study, hemodynamics in the normal and diseased left anterior descending (LAD) arterial trees were compared to estimate the effects of stenosis on the flow rate distribution, detailed flow field and wall shear stress (WSS) distributions on both transient and steady flow conditions. The anatomical model was reconstructed from computed tomography (CT) images of a normal porcine heart. Stenosis was imposed at the trunk of LAD tree with different area stenosis ratios (25%, 50%, 75% and 90%) to study the effect of stenosis. It has been found that the steady-state simulation is competent to demonstrate the prominent flow features of transit flow. Only significant stenosis altered the flow rate considerably, and its effects diminished at the downstream bifurcations. Based on the WSS distribution, non-significant stenosis potentially led to the development of atherosclerosis near the ostium of the downstream side branch and significant stenosis might promote atherosclerosis in its wake.

Effects of Cyclic Motion on Coronary Blood Flow

The goal of this study was to establish a computational fluid dynamics model to investigate the effect of cyclic motion (i.e., bending and stretching) on coronary blood flow. The three-dimensional (3D) geometry of a 50-mm section of the left anterior descending artery (normal or with a 60% stenosis) was constructed based on anatomical studies. To describe the bending motion of the blood vessel wall, arbitrary Lagrangian–Eularian methods were used. To simulate artery bending and blood pressure change induced stretching, the arterial wall was modeled as an anisotropic nonlinear elastic solid using the five-parameter Mooney–Rivlin hyperelastic model. Employing a laminar model, the flow field was solved using the continuity equations and Navier–Stokes equations. Blood was modeled as an incompressible Newtonian fluid. A fluid–structure interaction approach was used to couple the fluid domain and the solid domain iteratively, allowing force and total mesh displacement to be transferred between the two domains. The results demonstrated that even though the bending motion of the coronary artery could significantly affect blood cell trajectory, it had little effect on flow parameters, i.e., blood flow velocity, blood shear stress, and wall shear stress. The shape of the stenosis (asymmetric or symmetric) hardly affected flow parameters either. However, wall normal stresses (axial, circumferential, and radial stress) can be greatly affected by the blood vessel wall motion. The axial wall stress was significantly higher than the circumferential and radial stresses, as well as wall shear stress. Therefore, investigation on effects of wall stress on blood vessel wall cellular functions may help us better understand the mechanism of mechanical stress induced cardiovascular disease.

Towards patient-specific modeling of coronary hemodynamics in healthy and diseased state

A model describing the primary relations between the cardiac muscle and coronary circulation might be useful for interpreting coronary hemodynamics in case multiple types of coronary circulatory disease are present. The main contribution of the present study is the coupling of a microstructure-based heart contraction model with a 1D wave propagation model. The 1D representation of the vessels enables patient-specific modeling of the arteries and/or can serve as boundary conditions for detailed 3D models, while the heart model enables the simulation of cardiac disease, with physiology-based parameter changes. Here, the different components of the model are explained and the ability of the model to describe coronary hemodynamics in health and disease is evaluated. Two disease types are modeled: coronary epicardial stenoses and left ventricular hypertrophy with an aortic valve stenosis. In all simulations (healthy and diseased), the dynamics of pressure and flow qualitatively agreed with observations described in literature. We conclude that the model adequately can predict coronary hemodynamics in both normal and diseased state based on patient-specific clinical data.

Hemodynamic impacts of various types of stenosis in the left coronary artery bifurcation: a patient-specific analysis

This study investigates the hemodynamic changes to various types of coronary stenosis in the left coronary artery bifurcation, based on a patient-specific analysis. Twenty two patients with left coronary artery disease were included in this study. All stenoses involving the left coronary artery bifurcation were classified into four types, according to their locations: A) left circumflex (LCx) and left anterior descending (LAD), B) LCx only, C) left main stem only, and D) LAD only. Computational fluid dynamics (CFD) was performed to analyze the flow and wall shear stress (WSS) changes in all reconstructed left coronary geometries. Our results showed that the flow velocity and WSS were significantly increased at stenotic locations. High WSS was found at >70% lumen stenosis, which ranged from 2.5 Pa to 3.5 Pa. This study demonstrates that in patients with more than 50% stenosis in the left coronary artery bifurcation, WSS plays an important role in providing information about the extent of coronary atherosclerosis in the left coronary artery branch.

Impact of plaques in the left coronary artery on wall shear stress and pressure gradient in coronary side branches

In this study, we investigate plaques located at the left coronary bifurcation. We focus on the effect that the resulting changes in wall shear stress (WSS) and wall pressure stress gradient (WPSG) have on atherosclerotic progress in coronary artery disease. Coronary plaques were simulated and placed at the left main stem and the left anterior descending to produce >50% narrowing of the coronary lumen. Computational fluid dynamics analysis was carried out, simulating realistic physiological conditions that show the in vivo cardiac haemodynamic. WSS and WPSG in the left coronary artery were calculated and compared in the left coronary models, with and without the presence of plaques during cardiac cycles. Our results showed that WSS decreased while WPSG was increased in coronary side branches due to the presence of plaques. There is a direct correlation between coronary plaques and subsequent WSS and WPSG variations based on the bifurcation plaques simulated in the realistic coronary models.

Computational fluid dynamics analysis of the effect of plaques in the left coronary artery

This study was to investigate the hemodynamic effect of simulated plaques in left coronary artery models, which were generated from a sample patient's data. Plaques were simulated and placed at the left main stem and the left anterior descending (LAD) to produce at least 60% coronary stenosis. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics, and comparison of wall shear stress (WSS) between Newtonian and non-Newtonian fluid models was performed. The pressure gradient (PSG) and flow velocities in the left coronary artery were measured and compared in the left coronary models with and without presence of plaques during cardiac cycle. Our results showed that the highest PSG was observed in stenotic regions caused by the plaques. Low flow velocity areas were found at postplaque locations in the left circumflex, LAD, and bifurcation. WSS at the stenotic locations was similar between the non-Newtonian and Newtonian models although some more details were observed with non-Newtonian model. There is a direct correlation between coronary plaques and subsequent hemodynamic changes, based on the simulation of plaques in the realistic coronary models.

Effects of physiologically relevant dynamic shear stress on platelet complement activation

Disturbed shear stress, commonly found in cardiovascular diseases, plays important roles in platelet activation and functions. It has been reported that when activated by elevated shear stress, platelets were able to support complement activation to completion. In this study, through a dynamic cone and plate shearing device, three physiologically relevant shear stresses were applied to platelets, mimicking the shear conditions when platelets pass through a normal left coronary artery (0.05-1 Pa), a 60% stenosis (elevated shear stress at 6.5 Pa for less than 0.1 s), and when platelets are trapped in a recirculation zone past a stenosis (<0.5 Pa). After shear exposure, platelet-surface complement activation (C1q, C4d, iC3b, and SC5b-9 depositions) was measured using a solid-phase ELISA approach and flow cytometry. Production of complement regulatory proteins - C1-inhibitor (C1-INH) and complement receptor 1 (CR1), was also measured. Results demonstrated that low-pulsatile shear stress (recirculation) was able to initiate platelet complement activation, by increasing C1q deposition significantly. Both pathological shear stresses triggered significant increases in C1 inhibitor generation and noticeable changes in CR1 production, effectively preventing complement activation from completion. These results suggested that for platelets, low-pulsatile shear stress may be more pro-atherogenic, compared to elevated shear stress, especially when the shear stress exposure time is short.

Computation of hemodynamics in the left coronary artery with variable angulations

The purpose of this study was to investigate the hemodynamic effect of variations in the angulations of the left coronary artery, based on simulated and realistic coronary artery models. Twelve models consisting of four realistic and eight simulated coronary artery geometries were generated with the inclusion of left main stem, left anterior descending and left circumflex branches. The simulated models included various coronary artery angulations, namely, 15°, 30°, 45°, 60°, 75°, 90°, 105° and 120°. The realistic coronary angulations were based on selected patient's data with angles ranging from narrow angles of 58° and 73° to wide angles of 110° and 120°. Computational fluid dynamics analysis was performed to simulate realistic physiological conditions that reflect the in vivo cardiac hemodynamics. The wall shear stress, wall shear stress gradient, velocity flow patterns and wall pressure were measured in simulated and realistic models during the cardiac cycle. Our results showed that a disturbed flow pattern was observed in models with wider angulations, and wall pressure was found to reduce when the flow changed from the left main stem to the bifurcated regions, based on simulated and realistic models. A low wall shear stress gradient was demonstrated at left bifurcations with wide angles. There is a direct correlation between coronary angulations and subsequent hemodynamic changes, based on realistic and simulated models. Further studies based on patients with different severities of coronary artery disease are required to verify our results.