Does low and oscillatory wall shear stress correlate spatially with early atherosclerosis? A systematic review

Computational haemodynamics for pulmonary valve replacement by means of a reduced fluid-structure interaction model

Pulmonary valve replacement (PVR) consists of substituting a patient's original valve with a prosthetic one, primarily addressing pulmonary valve insufficiency, which is crucially relevant in Tetralogy of Fallot repairment. While extensive clinical and computational literature on aortic and mitral valve replacements is available, PVR's post-procedural haemodynamics in the pulmonary artery and the impact of prosthetic valve dynamics remain significantly understudied. Addressing this gap, we introduce a reduced Fluid-Structure Interaction (rFSI) model, applied for the first time to the pulmonary valve. This model couples a three-dimensional computational representation of pulmonary artery haemodynamics with a one-degree-of-freedom model to account for valve structural mechanics. Through this approach, we analyse patient-specific haemodynamics pre and post PVR. Patient-specific geometries, reconstructed from CT scans, are virtually equipped with a template valve geometry. Boundary conditions for the model are established using a lumped-parameter model, fine-tuned based on clinical patient data. Our model accurately reproduces patient-specific haemodynamic changes across different scenarios: pre-PVR, six months post-PVR, and a follow-up condition after a decade. It effectively demonstrates the impact of valve implantation on sustaining the diastolic pressure gradient across the valve. The numerical results indicate that our valve model is able to reproduce overall physiological and/or pathological conditions, as preliminary assessed on two different patients. This promising approach provides insights into post-PVR haemodynamics and prosthetic valve effects, shedding light on potential implications for patient-specific outcomes.

Numerical analysis of blood flow in the abdominal aorta under simulated weightlessness and earth conditions

Blood flow through the abdominal aorta and iliac arteries is a crucial area of research in hemodynamics and cardiovascular diseases. To get in to the problem, this study presents detailed analyses of blood flow through the abdominal aorta, together with left and right iliac arteries, under Earth gravity and weightless conditions, both at the rest stage, and during physical activity. The analysis were conducted using ANSYS Fluent software. The results indicate, that there is significantly less variation in blood flow velocity under weightless conditions, compared to measurement taken under Earth Gravity conditions. Study presents, that the maximum and minimum blood flow velocities decrease and increase, respectively, under weightless conditions. Our model for the left iliac artery revealed higher blood flow velocities during the peak of the systolic phase (systole) and lower velocities during the early diastolic phase (diastole). Furthermore, we analyzed the shear stress of the vessel wall and the mean shear stress over time. Additionally, the distribution of oscillatory shear rate, commonly used in hemodynamic analyses, was examined to assess the effects of blood flow on the blood vessels. Countermeasures to mitigate the negative effects of weightlessness on astronauts health are discussed, including exercises performed on the equipment aboard the space station. These exercises aim to maintain optimal blood flow, prevent the formation of atherosclerotic plaques, and reduce the risk of cardiovascular complications.

Exploring the relationship between hemodynamics and the immune microenvironment in carotid atherosclerosis: Insights from CFD and CyTOF technologies

Carotid atherosclerosis is a major cause of stroke. Hemodynamic forces, such as shear stress and oscillatory shear, play an important role in the initiation and progression of atherosclerosis. The alteration of the immune microenvironment is the fundamental pathological mechanism by which diverse external environmental factors impact the formation and progression of plaques. However, Current research on the relationship between hemodynamics and immunity in atherosclerosis still lack of comprehensive understanding. In this study, we combined computational fluid dynamics (CFD) and Mass cytometry (CyTOF) technologies to explore the changes in the immune microenvironment within plaques under different hemodynamic conditions. Our results indicated that neutrophils were enriched in adverse flow environments. M2-like CD163+CD86+ macrophages were predominantly enriched in high WSS and low OSI environments, while CD163-CD14+ macrophages were enriched in low WSS and high OSI environments. Functional analysis further revealed T cell pro-inflammatory activation and dysregulation in modulation, along with an imbalance in M1-like/M2-like macrophages, suggesting their potential involvement in the progression of atherosclerotic lesions mediated by adverse flow patterns. Our study elucidated the potential mechanisms by which hemodynamics regulated the immune microenvironment within plaques, providing intervention targets for future precision therapies.

A dual-targeting therapeutic nanobubble for imaging-guided atherosclerosis treatment

Atherosclerosis is a cardiovascular disease that seriously endangers human health. Low shear stress (LSS) is recognized as a vital factor in causing chronic inflammatory and further inducing the occurrence and development of atherosclerosis. Targeting imaging and treatment are of substantial significance for the diagnosis and therapy of atherosclerosis. On this ground, a kind of ultrasound (US) imaging-guided therapeutic polymer nanobubbles (NBs) with dual targeting of magnetism and antibody was rationally designed and constructed for the efficiently treating LSS-mediated atherosclerosis. Under the combined targeting effect of an external magnetic field and antibodies, the drug-loaded therapeutic NBs can be effectively accumulated in the inflammatory area caused by LSS. Upon US irradiation, the NBs can be selectively disrupted, leading to the rapid release of the loaded drugs at the targeted site. Notably, the US irradiation generates a cavitation effect that induces repairable micro gaps in nearby cells, thereby enhancing the uptake of released drugs and further improving the therapeutic effect. The prominent US imaging, efficient anti-inflammatory effect and treatment outcome of LSS-mediated atherosclerosis had been verified in vivo on a surgically constructed LSS-atherosclerosis animal model. This work showcased the potential of the designed NBs with multifunctionality for in vivo imaging, dual-targeting, and drug delivery in the treatment of atherosclerosis.

Low shear stress-induced blockage of autophagic flux impairs endothelial barrier and facilitates atherosclerosis in mice

Atherosclerosis preferentially occurs in areas with low shear stress (LSS) and oscillatory flow. LSS has been demonstrated to correlate with the development of atherosclerosis. The sphingosine 1-phosphate receptor 1 (S1PR1), involving intravascular blood flow sensing, regulates vascular development and vascular barrier function. However, whether LSS affects atherosclerosis via regulating S1PR1 remains incompletely clear. In this study, immunostaining results of F-actin, β-catenin, and VE-cadherin indicated that LSS impaired endothelial barrier function in human umbilical vein endothelial cells (HUVECs). Western blot analysis showed that LSS resulted in blockage of autophagic flux in HUVECs. In addition, autophagy agonist Rapamycin (Rapa) antagonized LSS-induced endothelial barrier dysfunction, whereas autophagic flux inhibitor Bafilomycin A1 (BafA1) exacerbated it, indicating that LSS promoted endothelial barrier dysfunction by triggering autophagic flux blockage. Notably, gene expression analysis revealed that LSS downregulated S1PR1 expression, which was antagonized by Rapa. Selective S1PR1 antagonist W146 impaired endothelial barrier function of HUVECs under high shear stress (HSS) conditions. Moreover, our data showed that expression of GAPARAPL2, a member of autophagy-related gene 8 (Atg8) proteins, was decreased in HUVECs under LSS conditions. Autophagic flux blockage induced by GAPARAPL2 knockdown inhibited S1PR1, aggravated endothelial barrier dysfunction of HUVECs in vitro, and promoted aortic atherosclerosis in ApoE-/- mice in vivo. Our study demonstrates that autophagic flux blockage induced by LSS downregulates S1PR1 expression and impairs endothelial barrier function. GABARAPL2 inhibition is involved in LSS-induced autophagic flux blockage, which impairs endothelial barrier function via downregulation of S1PR1.

Structural and Mechanical Properties of Human Superficial Femoral and Popliteal Arteries

The femoropopliteal artery (FPA) is the main artery in the lower limb. It supplies blood to the leg muscles and undergoes complex deformations during limb flexion. Atherosclerotic disease of the FPA (peripheral arterial disease, PAD) is a major public health burden, and despite advances in surgical and interventional therapies, the clinical outcomes of PAD repairs continue to be suboptimal, particularly in challenging calcified lesions and biomechanically active locations. A better understanding of human FPA mechanical and structural characteristics in relation to age, risk factors, and the severity of vascular disease can help develop more effective and longer-lasting treatments through computational modeling and device optimization. This review aims to summarize recent research on the main biomechanical and structural properties of human superficial femoral and popliteal arteries that comprise the FPA and describe their anatomy, composition, and mechanical behavior under different conditions.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstructs flow and increases pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the small pulmonary arteries due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we applied a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) approach to investigate the impact of CTEPH on pulmonary artery stiffening, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) in extralobar (main, right, and left) pulmonary arteries and intralobar (distal to the extralobar) arteries. Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar and intralobar arteries. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI, quantified by the novel parameter φ , which is related to thrombogenicity. Our analysis reveals a strong positive correlation between increases in mean pulmonary artery pressure (mPAP) and φ from baseline to CTEPH in individual subjects, which supports the suggestion that increased φ drives disease severity. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

Effect of PCI on ophthalmic artery hemodynamics in patients with acute coronary syndrome

Purpose

We aimed to explore the effects of percutaneous coronary intervention (PCI) on the ophthalmic artery (OA) hemodynamics in patients with acute coronary syndrome (ACS).

Methods

A total of 73 participants (Group0: healthy controls, Group1: Patients with ACS underwent PCI < 3 months, Group2: Patients with ACS underwent PCI ≥ 3 months) were enrolled. Computed tomographic angiography images were used to construct three-dimensional models of participants' OAs. Numerical simulations based on computational fluid dynamics were used to acquire hemodynamic parameters.

Results

The angle between the OA and internal carotid artery in Group2 was significantly larger compared with Group0 and Group1 (P = 0.003 and P = 0.044). Hemodynamic simulation showed a significantly slower OA blood velocity in Group1 than in the control (P < 0.001) and Group2 (P = 0.033). Lower wall shear stress was found in Group1 than that in control (P = 0.040). Patients after PCI had a higher wall pressure than healthy controls (P = 0.012 and P = 0.004). Mass flow ratios were decreased in Group1 and Group2 (P = 0.021 and P = 0.002). The hemodynamic parameters of OA were correlated with several clinical indicators.

Conclusions

The OA blood flow velocity of patients with ACS after PCI initially slowed down, which increased the risk of plaque formation, and then showed an increasing trend. There was a correlation between OA hemodynamic parameters and clinical indexes related to cardiac stress. Ischemia-reperfusion injury and changes in blood flow status after PCI may affect OA morphology and hemodynamics, leading to ocular lesions.

Trial registration

ChiCTR2100050428.

Longitudinal wall shear stress evaluation using centerline projection approach in the numerical simulations of the patient-based carotid artery

In this numerical study, areas of the carotid bifurcation and of a distal stenosis in the internal carotid artery are closely observed to evaluate the patient's current risks of ischemic stroke. An indicator for the vessel wall defects is the stress exerted by blood on the vessel tissue, typically expressed by the amplitude of the wall shear stress vector (WSS) and its oscillatory shear index. To detect negative shear stresses corresponding with reversal flow, we perform orientation-based shear evaluation. We investigate the longitudinal component of the wall shear vector, where tangential vectors aligned longitudinally with the vessel are necessary. However, resulting from imaging segmentation resolution of patients' computed tomography angiography scans and stenotic regions, the geometry model's mesh is non-smooth on its surface areas and the automatically generated tangential vector field is discontinuous and multi-directional, making an interpretation of our orientation-based risk indicators unreliable. We improve the evaluation of longitudinal shear stress by applying the projection of the vessel's centerline to the surface to construct smooth tangential field aligned longitudinally with the vessel. We validate our approach for the longitudinal WSS component and the corresponding oscillatory index by comparing them to results obtained using automatically generated tangents in both rigid and elastic vessel modeling and to amplitude-based indicators. We present the major benefit of our longitudinal WSS evaluation based on its directionality for the cardiovascular risk assessment, which is the detection of negative WSS indicating persistent reversal or transverse flow. This is impossible in the case of the amplitude-based WSS.

Sphingomyelinase-responsive nanomicelles for targeting atherosclerosis

Atherosclerosis, a leading cause of cardiovascular diseases requires approaches to enhance disease monitoring and treatment. Nanoparticles offer promising potential in this area by being customisable to target components or molecular processes within plaques, while carrying diagnostic and therapeutic agents. However, the number of biomarkers available to target this disease is limited. This study investigated the use of sphingomyelin-based nanomicelles triggered by sphingomyelinase (SMase) in atherosclerotic plaques. Accumulation of iron oxide-based nanomicelles in the plaque was demonstrated by fluorescence, MR imaging and electron microscopy. These findings demonstrate the possibility of utilising SMase as a mechanism to retain nanoprobes within plaques, thus opening up possibilities for future therapeutic interventions.

In silico hemodynamical simulations show additional benefits of artery wall softening induced by antihypertensive drugs

BACKGROUND AND OBJECTIVES

Age-related arterial stiffening increases peripheral resistance and decreases arterial distensibility, thus contributing to hypertension, an important risk factor of atherosclerosis. It causes abnormal blood flow, endothelial dysfunction, higher pulse wave velocity, and consequently elevated pressure wave amplitude.

METHODS

This paper presents the influence of these changes via multiscale 3D-0D transient computational fluid dynamics simulations of blood flow in five patient-specific geometries of human carotid bifurcation using archetypal flow waveforms for young and old subjects.

RESULTS

The proposed model shows a significant decrease in the time-averaged wall shear stress (TAWSS) for the old archetypal flow waveform. This is in good agreement with clinical data on a straight segment of common carotid arteries available for young and old subjects. Moreover, our study showed that the decrease of area-averaged TAWSS related to the old flow waveform is much more pronounced (2.5 ÷ 4.5 times higher) at risk areas (areas showing TAWSS below its threshold value of 0.48 Pa) than in straight segments commonly considered in clinical studies.

CONCLUSIONS

Since arterial stiffness can be lowered through long-term usage of any of the five basic groups of antihypertensives, possible benefits of such medical therapy could be not only lowering blood pressure and peripheral resistance but also in increasing the TAWSS and thus attenuating an important mechanism of the atherosclerotic process.

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.

Near-wall hemodynamic parameters of finger arteries altered by hand-transmitted vibration

BACKGROUND

Sustained exposure to high-level hand-transmitted vibrations may result in angioneurotic disorders, which partly originate from vibration-altered hemodynamics in the finger arteries when repeating these disturbances throughout working life. Hence, the aim of this study is to assess the most relevant hemodynamic descriptors in the digital arteries, determine the relationship between the latter and vibration features, and gain better understanding of the physiological mechanisms involved.

METHODS

An experimental setup, mainly comprised of an ultra-high frequency ultrasound scanner and a vibration shaker, was used to image the digital proper volar arteries of the forefinger. Raw ultrasound data were post-processed by custom-made numerical routines to supply a pulsatile fluid mechanics model for computing the hemodynamic descriptors. Twenty-four healthy volunteers participated in the measurement campaign. Classical statistical methods were then applied to the dataset and also the wavelet transform for calculating the signal power in the frequency bands matching cardiac, respiratory, myogenic and neurogenic activities.

RESULTS

The artery diameter, the wall shear stress - WSS - and the WSS temporal gradient - WSSTG - were found to be the most relevant descriptors. Vibration-induced WSS was divided by three compared to its basal value whatever the vibration frequency and it was proportional to log2 of the acceleration level. Marked increases in WSSTG when stopping vibration might also lead to adverse health effects. Vibration caused a drop in WSS power for the frequency band associated with the neurogenic activity of the sympathetic nervous system.

CONCLUSION

This study may pave the way for a new framework to prevent vibration-induced vascular risk.

Integrated systems biology approach identifies gene targets for endothelial dysfunction

Endothelial dysfunction (ED) is critical in the development and progression of cardiovascular (CV) disorders, yet effective therapeutic targets for ED remain elusive due to limited understanding of its underlying molecular mechanisms. To address this gap, we employed a systems biology approach to identify potential targets for ED. Our study combined multi omics data integration, with siRNA screening, high content imaging and network analysis to prioritise key ED genes and identify a pro- and anti-ED network. We found 26 genes that, upon silencing, exacerbated the ED phenotypes tested, and network propagation identified a pro-ED network enriched in functions associated with inflammatory responses. Conversely, 31 genes ameliorated ED phenotypes, pointing to potential ED targets, and the respective anti-ED network was enriched in hypoxia, angiogenesis and cancer-related processes. An independent screen with 17 drugs found general agreement with the trends from our siRNA screen and further highlighted DUSP1, IL6 and CCL2 as potential candidates for targeting ED. Overall, our results demonstrate the potential of integrated system biology approaches in discovering disease-specific candidate drug targets for endothelial dysfunction.

Age and the Residual Risk of Cardiovascular Disease following Low Density Lipoprotein-Cholesterol Exposure

We believe that there is sufficient evidence from basic science, longitudinal cohort studies and randomised controlled trials which validates the low-density lipoprotein cholesterol (LDL-C) or lipid hypothesis. It is important that we can communicate details of the cardiovascular disease (CVD) risk reduction that the average patient could expect depending on the scale of LDL-C decrease following lipid lowering therapy. It is also essential that residual risk (ResR) of CVD be highlighted. To achieve this aim by using existing trial evidence, we developed mathematical models initially for relative risk reduction (RRR) and absolute risk (AR) reduction and then showed that despite optimising LDL-C levels, a considerable degree of ResR remains that is dependent on AR. Age is significantly associated with AR (odds ratio: 1.02, 95% confidence intervals: 1.01-1.04) as was previously demonstrated by analysing the Whickham study cohort using a logistic regression model (age remaining significant even when all the other significant risk factors such as sex, smoking, systolic blood pressure, diabetes and family history were included in the regression model). A discussion of a paper by Ference et al. provided detailed evidence of the relationship between age and AR, based on lifetime LDL-C exposure. Finally, we discussed non-traditional CVD risk factors that may contribute to ResR based on randomised controlled trials investigating drugs improving inflammation, thrombosis, metabolic and endothelial status.

Mechanical wall stress and wall shear stress are associated with atherosclerosis development in non-calcified coronary segments

BACKGROUND AND AIMS

Atherosclerotic plaque onset and progression are known to be affected by local biomechanical factors. While the role of wall shear stress (WSS) has been studied, the impact of another biomechanical factor, namely mechanical wall stress (MWS), remains poorly understood. In this study, we investigated the association of MWS, independently and combined with WSS, towards atherosclerosis in coronary arteries.

METHODS

Thirty-four human coronary arteries were analyzed using near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) and optical coherence tomography (OCT) at baseline and after 12 months. Baseline WSS and MWS were calculated using computational models, and wall thickness (ΔWT) and lipid-rich necrotic core size (ΔLRNC) change were measured in non-calcified coronary segments. The arteries were further divided into 1.5 mm/45° sectors and categorized as plaque-free or plaque sectors. For each category, associations between biomechanical factors (WSS & MWS) and changes in coronary wall (ΔWT & ΔLRNC) were studied using linear mixed models.

RESULTS

In plaque-free sectors, higher MWS (p < 0.001) was associated with greater vessel wall growth. Plaque sectors demonstrated wall thickness reduction over time, likely due to medical therapy, where higher levels of WSS and WMS, individually and combined, (p < 0.05) were associated with a greater reduction. Sectors with low MWS combined with high WSS demonstrated the highest LRNC increase (p < 0.01).

CONCLUSIONS

In this study, we investigated the association of the (largely-overlooked) biomechanical factor MWS with coronary atherosclerosis, individually and combined with WSS. Our results demonstrated that both MWS and WSS significantly correlate with atherosclerotic plaque initiation and development.

Microarray Expression Profile and Bioinformatic Analysis of Circular RNA in Human Arteriosclerosis Obliterans

Background

Arteriosclerosis obliterans (ASO) is the leading cause of nontraumatic lower-extremity amputations. Multiple researches have suggested that circular RNAs (circRNAs) played vital regulatory functions in cancer and cardiovascular disease. Nevertheless, the underlying effect and pathological mechanism of circRNAs in the formation and progression of ASO are still indistinct.

Methods and Results

This study used microarray analysis to investigate the expression portrait of circRNAs in normal lower extremity arteries and ASO arteries. Bioinformatics analysis was conducted using the KEGG database to study the enrichment of differentially expressed circRNAs (DE circRNAs) and predict their functions. The accuracy of microarray assay was verified by evaluating expression of the top 5 upregulated and 5 downregulated circRNAs (raw density of normal group ≥200) using RT-qPCR. A circRNA-miRNA-mRNA interaction network was further predicted using software. Compared to the normal lower extremity group, the ASO arteries with HE and EVG staining presented hyperplastic fibrous membrane and luminal stenosis. A total of 12,735 circRNAs were identified, including 1196 DE circRNAs with 276 upregulated and 920 downregulated in ASO group based on |log2(FC)| > 1 and padj < 0.05. Among selected 10 circRNAs, RT-qPCR confirmed that hsa_circ_0003266, hsa_circ_0118936 and hsa_circ_0067161 were upregulated while hsa_circ_0091934 and hsa_circ_0092022 were downregulated in ASO group (p < 0.05). GO analysis presented that the DE circRNAs were primarily enriched in protein binding, intracellular part and organelle organization. KEGG pathway analysis indicated that MAPK signaling pathway, human T-cell leukemia virus 1 infection, proteoglycans in cancer were associated with the DE circRNAs. The circRNA-miRNA-mRNA interactive network revealed that both mRNAs and miRNAs linked to circRNAs played an indispensable role in ASO.

Conclusion

This study described the expression portrait of circRNAs in human ASO arteries, and revealed the molecular background for further investigations of the circRNA regulatory mechanism in the formation and progression of ASO.

Mechanosignalling pathways that regulate endothelial barrier function

Blood vessels are lined by a single layer of endothelial cells that provide a barrier between circulating plasma and the underlying tissue. Permeability of endothelial cells is tightly regulated, and increased permeability is associated with a number of diseases including atherosclerosis. Endothelial cells are continuously exposed to mechanical forces exerted by flowing blood and are particularly sensitive to shear stress, which is a key determinant of endothelial function. Undisturbed flow promotes endothelial resilience and reduces permeability to macromolecules whereas disturbed flow promotes endothelial dysfunction and barrier disruption. This review will outline recent advances in our understanding of how disturbed and undisturbed flow regulate paracellular and transcellular permeability and will highlight potential cellular targets that could form the basis of therapies to limit the development of cardiovascular disease.

The Impact of Left Ventricular Assist Device Outflow Graft Positioning on Aortic Hemodynamics: Improving Flow Dynamics to Mitigate Aortic Insufficiency

Heart failure is a global health concern with significant implications for healthcare systems. Left ventricular assist devices (LVADs) provide mechanical support for patients with severe heart failure. However, the placement of the LVAD outflow graft within the aorta has substantial implications for hemodynamics and can lead to aortic insufficiency during long-term support. This study employs computational fluid dynamics (CFD) simulations to investigate the impact of different LVAD outflow graft locations on aortic hemodynamics. The introduction of valve morphology within the aorta geometry allows for a more detailed analysis of hemodynamics at the aortic root. The results demonstrate that the formation of vortex rings and subsequent vortices during the high-velocity jet flow from the graft interacted with the aortic wall. Time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI) indicate that modification of the outflow graft location changes mechanical states within the aortic wall and aortic valve. Among the studied geometric factors, both the height and inclination angle of the LVAD outflow graft are important in controlling retrograde flow to the aortic root, while the azimuthal angle primarily determines the rotational direction of blood flow in the aortic arch. Thus, precise positioning of the LVAD outflow graft emerges as a critical factor in optimizing patient outcomes by improving the hemodynamic environment.

Considering the Influence of Coronary Motion on Artery-Specific Biomechanics Using Fluid-Structure Interaction Simulation

The endothelium in the coronary arteries is subject to wall shear stress and vessel wall strain, which influences the biology of the arterial wall. This study presents vessel-specific fluid-structure interaction (FSI) models of three coronary arteries, using directly measured experimental geometries and boundary conditions. FSI models are used to provide a more physiologically complete representation of vessel biomechanics, and have been extended to include coronary bending to investigate its effect on shear and strain. FSI both without- and with-bending resulted in significant changes in all computed shear stress metrics compared to CFD (p = 0.0001). Inclusion of bending within the FSI model produced highly significant changes in Time Averaged Wall Shear Stress (TAWSS) + 9.8% LAD, + 8.8% LCx, - 2.0% RCA; Oscillatory Shear Index (OSI) + 208% LAD, 0% LCx, + 2600% RCA; and transverse wall Shear Stress (tSS) + 180% LAD, + 150% LCx and + 200% RCA (all p < 0.0001). Vessel wall strain was homogenous in all directions without-bending but became highly anisotropic under bending. Changes in median cyclic strain magnitude were seen for all three vessels in every direction. Changes shown in the magnitude and distribution of shear stress and wall strain suggest that bending should be considered on a vessel-specific basis in analyses of coronary artery biomechanics.

Subject-specific one-dimensional fluid dynamics model of chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) develops due to the accumulation of blood clots in the lung vasculature that obstruct flow and increase pressure. The mechanobiological factors that drive progression of CTEPH are not understood, in part because mechanical and hemodynamic changes in the pulmonary vasculature due to CTEPH are not easily measurable. Using previously published hemodynamic measurements and imaging from a large animal model of CTEPH, we developed a subject-specific one-dimensional (1D) computational fluid dynamic (CFD) models to investigate the impact of CTEPH on pulmonary artery stiffening, time averaged wall shear stress (TAWSS), and oscillatory shear index (OSI). Our results demonstrate that CTEPH increases pulmonary artery wall stiffness and decreases TAWSS in extralobar (main, right and left pulmonary arteries) and intralobar vessels. Moreover, CTEPH increases the percentage of the intralobar arterial network with both low TAWSS and high OSI. This subject-specific experimental-computational framework shows potential as a predictor of the impact of CTEPH on pulmonary arterial hemodynamics and pulmonary vascular mechanics. By leveraging advanced modeling techniques and calibrated model parameters, we predict spatial distributions of flow and pressure, from which we can compute potential physiomarkers of disease progression, including the combination of low mean wall shear stress with high oscillation. Ultimately, this approach can lead to more spatially targeted interventions that address the needs of individual CTEPH patients.

Hemodynamic parameters impact the stability of distal stent graft-induced new entry

Stent graft-induced new entry tear (SINE) is a serious complication in aortic dissection patients caused by the stent-graft itself after thoracic endovascular aortic repair (TEVAR). The stability of SINE is a key indicator for the need and timing of reinterventions. This study aimed to understand the role of hemodynamics in SINE stability by means of computational fluid dynamics (CFD) analysis based on patient-specific anatomical information. Four patients treated with TEVAR who developed a distal SINE (dSINE) were included; two patients had a stable dSINE and two patients experienced expansion of the dSINE upon follow-up examinations. CFD simulations were performed on geometries reconstructed from computed tomography scans acquired upon early detection of dSINE in these patients. Computational results showed that stable dSINEs presented larger regions with low time-averaged wall shear stress (TAWSS) and high relative residence time (RRT), and partial thrombosis was observed at subsequent follow-ups. Furthermore, significant systolic antegrade flow was observed in the unstable dSINE which also had a larger retrograde flow fraction (RFF) on the SINE plane. In conclusion, this pilot study suggested that high RRT and low TAWSS may indicate stable dSINE by promoting thrombosis, whereas larger RFF and antegrade flows inside dSINE might be associated with its expansion.

The relationship between geometry and hemodynamics of the stenotic carotid artery based on computational fluid dynamics

OBJECTIVE

The purpose of this work was to investigate the relationship between the geometric factors and the hemodynamics of the stenotic carotid artery.

METHODS

We retrospectively reviewed data of patients with carotid stenosis (40%-95%). The Navier-Stokes equations were solved using ANSYS CFX 18.0. Correlation analysis was based on Spearman's test. Geometric variables (p < 0.1 in the univariate analysis) were entered into the logistical regression. A receiver-operating characteristics analysis was used to detect hemodynamically significant lesions.

RESULTS

81 patients (96 arteries) were evaluated. The logistic regression analysis revealed that the translesional pressure ratio was significantly correlated with the stenosis degree (OR = 1.147, p < 0.001) and the angle between internal carotid artery and external carotid artery (angle γ) (OR = 0.933, p = 0.01). The translesional wall shear stress ratio was significantly correlated with stenosis degree (OR = 1.094, p < 0.001), lesion length (OR = 0.873, p = 0.01), lumen area of internal carotid artery (OR = 0.867, p = 0.002), and lumen area of common carotid artery (OR = 1.058, p = 0.01). For predicting low translesional pressure ratio, the AUC was 0.71 (p < 0.001) for angle γ, and was 0.87 (p < 0.001) for stenosis degree. For predicting high translesional wall shear stress ratio, the AUC was 0.62 (p = 0.04) for lumen area of internal carotid artery, and was 0.77 (p < 0.001) for stenosis degree.

CONCLUSIONS

Apart from stenosis degree, other geometric characteristics of lesions may also have an influence on hemodynamics of the stenotic carotid artery.

Near-wall hemodynamic changes in subclavian artery perfusion induced by retrograde inner branched thoracic endograft implantation

Objective

Left subclavian artery (LSA)-branched endografts with retrograde inner branch configuration (thoracic branch endoprosthesis [TBE]) offer a complete endovascular solution when LSA preservation is required during zone 2 thoracic endovascular aortic repair. However, the hemodynamic consequences of the TBE have not been well-investigated. We compared near-wall hemodynamic parameters before and after the TBE implantation using computational fluid dynamic simulations.

Methods

Eleven patients who had undergone TBE implantation were included. Three-dimensional aortic arch geometries were constructed from the pre- and post-TBE implantation computed tomography images. The resulting 22 three-dimensional aortic arch geometries were then discretized into finite element meshes for computational fluid dynamic simulations. Inflow boundary conditions were prescribed using normal physiological pulsatile circulation. Outlet boundary conditions consisted of Windkessel models with previously published values. Blood flow, modeled as Newtonian fluid, simulations were performed with rigid wall assumptions using SimVascular's incompressible Navier-Stokes solver. We compared well-established hemodynamic descriptors: pressure, flow rate, time-averaged wall shear stress (TAWSS), the oscillatory shear index (OSI), and percent area with an OSI of >0.2. Data were presented on the stented portion of the LSA.

Results

TBE implantation was associated with a small decrease in peak LSA pressure (153 mm Hg; interquartile range [IQR], 151-154 mm Hg vs 159 mm Hg; IQR, 158-160 mm Hg; P = .005). No difference was observed in peak LSA flow rates before and after implantation: 40.4 cm³/ (IQR, 39.5-41.6 cm³/s) vs 41.3 cm³/s (IQR, 37.2-44.8 cm³/s; P = .59). There was a significant postimplantation increase in TAWSS (15.2 dynes/cm² [IQR, 12.2-17.7 dynes/cm²] vs 6.2 dynes/cm² [IQR, 5.7-10.3 dynes/cm²]; P = .003), leading to decreases in both the OSI (0.088 [IQR, 0.063 to -0.099] vs 0.1 [IQR, 0.096-0.16]; P = .03) and percentage of area with an OSI of >0.2 (10.4 [IQR, 5.8-15.8] vs 15.7 [IQR, 10.7-31.9]; P = .13). Neither LSA side branch angulation (median, 81°, IQR, 77°-109°) nor moderate compression (16%-58%) seemed to have an impact on the pressure, flow rate, TAWSS, or percentage of area with an OSI of >0.2 in the stented LSA.

Conclusions

The implantation of TBE produces modest hemodynamic disturbances that are unlikely to result in clinically relevant changes.

Aortic haemodynamics and wall stress analysis following arch aneurysm repair using a single-branched endograft

Introduction

Thoracic endovascular aortic repair (TEVAR) of the arch is challenging given its complex geometry and the involvement of supra-aortic arteries. Different branched endografts have been designed for use in this region, but their haemodynamic performance and the risk for post-intervention complications are not yet clear. This study aims to examine aortic haemodynamics and biomechanical conditions following TVAR treatment of an aortic arch aneurysm with a two-component single-branched endograft.

Methods

Computational fluid dynamics and finite element analysis were applied to a patient-specific case at different stages: pre-intervention, post-intervention and follow-up. Physiologically accurate boundary conditions were used based on available clinical information.

Results

Computational results obtained from the post-intervention model confirmed technical success of the procedure in restoring normal flow to the arch. Simulations of the follow-up model, where boundary conditions were modified to reflect change in supra-aortic vessel perfusion observed on the follow-up scan, predicted normal flow patterns but high levels of wall stress (up to 1.3M MPa) and increased displacement forces in regions at risk of compromising device stability. This might have contributed to the suspected endoleaks or device migration identified at the final follow up.

Discussion

Our study demonstrated that detailed haemodynamic and biomechanical analysis can help identify possible causes for post-TEVAR complications in a patient-specific setting. Further refinement and validation of the computational workflow will allow personalised assessment to aid in surgical planning and clinical decision making.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

Shear stress metrics associated with pro-atherogenic high-risk anatomical features in a carotid artery bifurcation model

BACKGROUND

Diseases associated with atherosclerotic plaques in the carotid artery are a major cause of deaths in the United States. Blood-flow-induced shear-stresses are known to trigger plaque formation. Prior literature suggests that the internal carotid artery sinus is prone to atherosclerosis, but there is limited understanding of why only certain patients are predisposed towards plaque formation.

METHODS

We computationally investigate the effect of vessel geometry on wall-shear-stress distribution by comparing flowfields and wall-shear-stress-metrics between a low-risk and a novel predisposed high-risk carotid artery bifurcation anatomy. Both models were developed based on clinical risk estimations and patient-averaged anatomical features. The high-risk geometry has a larger internal carotid artery branching angle and a lower internal-to-carotid-artery-diameter-ratio. A patient-averaged physiological carotid artery inflow waveform is used.

FINDINGS

The high-risk geometry experiences stronger flow separation in the sinus. Furthermore, it experiences a more equal flow split at the bifurcation, thereby reducing internal carotid artery flowrate and increasing atherosclerosis-prone low-velocity areas. Lowest time-averaged-wall-shear-stresses are present at the sinus outer wall, where plaques are often found, for both geometries. The high-risk geometry has significantly high, unfavorable oscillatory-shear-index values not found in the low-risk geometry. High oscillatory-shear-index areas are located at the vessels outside walls distal to the bifurcation and on the sinus wall.

INTERPRETATION

These results highlight the effectiveness of oscillatory-shear-index, to augment classical time-averaged-wall-shear-stress, in evaluating pro-atherogenic geometry features. Furthermore, the flow split at the bifurcation is a promising clinical indicator for atherosclerosis risk as it can be directly accessed using clinical imaging, whereas shear-stress-metrics cannot.

Modification of the swirling well cell culture model to alter shear stress metrics

Effects of hemodynamic shear stress on endothelial cells have been extensively investigated using the "swirling well" method, in which cells are cultured in dishes or multiwell plates placed on an orbital shaker. A wave rotates around the well, producing complex patterns of shear. The method allows chronic exposure to flow with high throughput at low cost but has two disadvantages: a number of shear stress characteristics change in a broadly similar way from the center to the edge of the well, and cells at one location in the well may release mediators into the medium that affect the behavior of cells at other locations, exposed to different shears. These properties make it challenging to correlate cell properties with shear. The present study investigated simple alterations to ameliorate these issues. Flows were obtained by numerical simulation. Increasing the volume of fluid in the well-altered dimensional but not dimensionless shear metrics. Adding a central cylinder to the base of the well-forced fluid to flow in a square toroidal channel and reduced multidirectionality. Conversely, suspending a cylinder above the base of the well made the flow highly multidirectional. Increasing viscosity in the latter model increased the magnitude of dimensional but not dimensionless metrics. Finally, tilting the well changed the patterns of different wall shear stress metrics in different ways. Collectively, these methods allow similar flows over most of the cells cultured and/or allow the separation of different shear metrics. A combination of the methods overcomes the limitations of the baseline model.

In silico experiments of intimal hyperplasia development: disendothelization in an axisymmetric idealized artery

We use in silico experiments to study the role of the hemodynamics and of the type of disendothelization on the physiopathology of intimal hyperplasia. We apply a multiscale bio-chemo-mechanical model of intimal hyperplasia on an idealized axisymmetric artery that suffers two kinds of disendothelizations. The model predicts the spatio-temporal evolution of the lesions development, initially localized at the site of damages, and after few days displaced downstream of the damaged zones, these two stages being observed whatever the kind of damage. Considering macroscopic quantities, the model sensitivity to pathology-protective and pathology-promoting zones is qualitatively consistent with experimental findings. The simulated pathological evolutions demonstrate the central role of two parameters: (a) the initial damage shape on the morphology of the incipient stenosis, and (b) the local wall shear stresses on the overall spatio-temporal dynamics of the lesion.

Pulmonary artery blood flow dynamics in chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension is caused by incomplete resolution and organization of thrombi. Blood flow dynamics are involved in thrombus formation; however, only a few studies have reported on pulmonary artery blood flow dynamics in patients with chronic thromboembolic pulmonary hypertension. Furthermore, the effects of treatment interventions on pulmonary artery blood flow dynamics are not fully understood. The aim of the study was to evaluate pulmonary artery blood flow dynamics in patients with chronic thromboembolic pulmonary hypertension before and after pulmonary endarterectomy and balloon pulmonary angioplasty, using computational fluid dynamics. We analyzed patient-specific pulmonary artery models of 10 patients with chronic thromboembolic pulmonary hypertension and three controls using computational fluid dynamics. In patients with chronic thromboembolic pulmonary hypertension, flow velocity and wall shear stress in the pulmonary arteries were significantly decreased, and the oscillatory shear index and blood stagnation volume were significantly increased than in controls. Pulmonary endarterectomy induced redistribution of pulmonary blood flow and improved blood flow dynamics in the pulmonary artery. Balloon pulmonary angioplasty improved pulmonary blood flow disturbance, decreased blood flow stagnation, and increased wall shear stress, leading to vasodilatation of the distal portion of the pulmonary artery following balloon pulmonary angioplasty treatment.

Multi-Objective Optimisation of a Novel Bypass Graft with a Spiral Ridge

The low long-term patency of bypass grafts is a major concern for cardiovascular treatments. Unfavourable haemodynamic conditions in the proximity of distal anastomosis are closely related to thrombus creation and lumen lesions. Modern graft designs address this unfavourable haemodynamic environment with the introduction of a helical component in the flow field, either by means of out-of-plane helicity graft geometry or a spiral ridge. While the latter has been found to lack in performance when compared to the out-of-plane helicity designs, recent findings support the idea that the existing spiral ridge grafts can be further improved in performance through optimising relevant design parameters. In the current study, robust multi-objective optimisation techniques are implemented, covering a wide range of possible designs coupled with proven and well validated computational fluid dynamics (CFD) algorithms. It is shown that the final set of suggested design parameters could significantly improve haemodynamic performance and therefore could be used to enhance the design of spiral ridge bypass grafts.

Pulse wave and vector flow Imaging for atherosclerotic disease progression in hypercholesterolemic swine

Non-invasive monitoring of atherosclerosis remains challenging. Pulse Wave Imaging (PWI) is a non-invasive technique to measure the local stiffness at diastolic and end-systolic pressures and quantify the hemodynamics. The objective of this study is twofold, namely (1) to investigate the capability of (adaptive) PWI to assess progressive change in local stiffness and homogeneity of the carotid in a high-cholesterol swine model and (2) to assess the ability of PWI to monitor the change in hemodynamics and a corresponding change in stiffness. Nine (n=9) hypercholesterolemic swine were included in this study and followed for up to 9 months. A ligation in the left carotid was used to cause a hemodynamic disturbance. The carotids with detectable hemodynamic disturbance showed a reduction in wall shear stress immediately after ligation (2.12 ± 0.49 to 0.98 ± 0.47 Pa for 40-90% ligation (Group B) and 1.82 ± 0.25 to 0.49 ± 0.46 Pa for >90% ligation (Group C)). Histology revealed subsequent lesion formation after 8-9 months, and the type of lesion formation was dependent on the type of the induced ligation, with more complex plaques observed in the carotids with a more significant ligation (C: >90%). The compliance progression appears differed for groups B and C, with an increase in compliance to 2.09 ± 2.90×10-10 m2 Pa-1 for group C whereas the compliance of group B remained low at 8 months (0.95 ± 0.94×10-10 m2 Pa-1). In summary, PWI appeared capable of monitoring a change in wall shear stress and separating two distinct progression pathways resulting in distinct compliances.

Wall shear stress and its role in atherosclerosis

Atherosclerosis (AS) is the major form of cardiovascular disease and the leading cause of morbidity and mortality in countries around the world. Atherosclerosis combines the interactions of systemic risk factors, haemodynamic factors, and biological factors, in which biomechanical and biochemical cues strongly regulate the process of atherosclerosis. The development of atherosclerosis is directly related to hemodynamic disorders and is the most important parameter in the biomechanics of atherosclerosis. The complex blood flow in arteries forms rich WSS vectorial features, including the newly proposed WSS topological skeleton to identify and classify the WSS fixed points and manifolds in complex vascular geometries. The onset of plaque usually occurs in the low WSS area, and the plaque development alters the local WSS topography. low WSS promotes atherosclerosis, while high WSS prevents atherosclerosis. Upon further progression of plaques, high WSS is associated with the formation of vulnerable plaque phenotype. Different types of shear stress can lead to focal differences in plaque composition and to spatial variations in the susceptibility to plaque rupture, atherosclerosis progression and thrombus formation. WSS can potentially gain insight into the initial lesions of AS and the vulnerable phenotype that gradually develops over time. The characteristics of WSS are studied through computational fluid dynamics (CFD) modeling. With the continuous improvement of computer performance-cost ratio, WSS as one of the effective parameters for early diagnosis of atherosclerosis has become a reality and will be worth actively promoting in clinical practice. The research on the pathogenesis of atherosclerosis based on WSS is gradually an academic consensus. This article will comprehensively review the systemic risk factors, hemodynamics and biological factors involved in the formation of atherosclerosis, and combine the application of CFD in hemodynamics, focusing on the mechanism of WSS and the complex interactions between WSS and plaque biological factors. It is expected to lay a foundation for revealing the pathophysiological mechanisms related to abnormal WSS in the progression and transformation of human atherosclerotic plaques.

EVA1A (Eva-1 Homolog A) Promotes Endothelial Apoptosis and Inflammatory Activation Under Disturbed Flow Via Regulation of Autophagy

BACKGROUND

Hemodynamic wall shear stress (WSS) exerted on the endothelium by flowing blood determines the spatial distribution of atherosclerotic lesions. Disturbed flow (DF) with a low WSS magnitude and reversing direction promotes atherosclerosis by regulating endothelial cell (EC) viability and function, whereas un-DF which is unidirectional and of high WSS magnitude is atheroprotective. Here, we study the role of EVA1A (eva-1 homolog A), a lysosome and endoplasmic reticulum-associated protein linked to autophagy and apoptosis, in WSS-regulated EC dysfunction.

METHODS

The effect of WSS on EVA1A expression was studied using porcine and mouse aortas and cultured human ECs exposed to flow. EVA1A was silenced in vitro in human ECs and in vivo in zebrafish using siRNA (small interfering RNA) and morpholinos, respectively.

RESULTS

EVA1A was induced by proatherogenic DF at both mRNA and protein levels. EVA1A silencing resulted in decreased EC apoptosis, permeability, and expression of inflammatory markers under DF. Assessment of autophagic flux using the autolysosome inhibitor, bafilomycin coupled to the autophagy markers LC3-II (microtubule-associated protein 1 light chain 3-II) and p62, revealed that EVA1A knockdown promotes autophagy when ECs are exposed to DF, but not un-DF . Blocking autophagic flux led to increased EC apoptosis in EVA1A-knockdown cells exposed to DF, suggesting that autophagy mediates the effects of DF on EC dysfunction. Mechanistically, EVA1A expression was regulated by flow direction via TWIST1 (twist basic helix-loop-helix transcription factor 1). In vivo, knockdown of EVA1A orthologue in zebrafish resulted in reduced EC apoptosis, confirming the proapoptotic role of EVA1A in the endothelium.

CONCLUSIONS

We identified EVA1A as a novel flow-sensitive gene that mediates the effects of proatherogenic DF on EC dysfunction by regulating autophagy.

Time-resolved simulation of blood flow through left anterior descending coronary artery: effect of varying extent of stenosis on hemodynamics

BACKGROUND AND OBJECTIVES

Real-time blood flow variation is crucial for understanding the dynamic development of coronary atherosclerosis. The main objective of this study is to investigate the effect of varying extent of stenosis on the hemodynamic features in left anterior descending coronary artery.

METHODS

Various Computational fluid dynamics (CFD) models were constructed with patient-specific CT image data, using actual fractional flow reserve (FFR) as boundary conditions to provide a real-time quantitative description of hemodynamic properties. The hemodynamic parameters, such as the local and instantaneous wall shear stress (WSS), oscillating shear index (OSI) and relative residence time (RRT), blood flow velocity and pressure drop during various phases of cardiac cycle were provided in detail.

RESULTS

There was no evident variation in hemodynamic parameters in the cases of less than 50% stenosis while there were abrupt and dramatic changes in hemodynamics when the stenosis aggravated from 60 to 70%. Furthermore, when the stenosis was beyond 70%, there existed substantial pressure difference, WSS, and blood flow velocity in the center of the stenosis. Although OSI and RRT increased along with the aggravation of stenosis, they appeared with obvious abnormalities across all cases, even in mild stenosis.

CONCLUSION

The simulation could present a dynamic and comprehensive profile of how hemodynamic parameters vary in accordance with divergent severities of stenosis, which could serve as an effective reference for the clinicians to have a deeper insight into the pathological mechanism of coronary atherosclerosis and stenosis.

Differences in Carotid Artery Geometry and Flow Caused by Body Postural Changes and Physical Exercise

Different body postures and physical exercises may lead to changes in arterial geometry and hemodynamics, which may be associated with the distribution of atherosclerosis lesions. This study was aimed at investigating potential geometric and hemodynamic changes of the carotid bifurcation in different body postures and after high-intensity interval training (HIIT) workouts. Three-dimensional vascular ultrasound (3DVUS) and Doppler ultrasound images were acquired for 21 healthy participants (aged 29 ± 6 y, 14 men and 7 women) in different body postures (sitting and three sleeping postures [supine, left lateral and right lateral]) and after physical exercises. The common carotid artery (CCA) and internal carotid artery (ICA) diameters of the left carotid artery were found to increase significantly from supine to left lateral (both p <0.05). CCA diameters (p < 0.05) and ICA/CCA diameter ratio (p < 0.01) of the left carotid artery changed significantly from supine to sitting. Significant differences in CCA peak systolic velocity (CCA PSV, p < 0.001), CCA end-diastolic velocity (CCA EDV, p < 0.001), CCA pulsatility index (CCA PI, p < 0.001) and maximum velocity-based wall shear stress at the CCA (WSS_(max) at the CCA, p < 0.001) were identified in different postures. After physical exercises, significant increases were observed in the CCA diameter (p < 0.001), CCA PSV (p < 0.001), ICA PSV (p < 0.05), WSS_(max) at the CCA (p < 0.001) and WSS_(max) at the ICA (p < 0.05), as were significantly lower values of the CCA EDV (p < 0.01) and ICA/CCA PSV ratio (p < 0.05). Side-to-side differences were also detected in different postural change scenarios and after physical exercise; more significant differences were found to occur only in the left-sided carotid artery. Significant differences were identified under postural change and after physical exercise among healthy adults, suggesting that daily activity has an effect on the carotid bifurcation. These changes may be associated with formation and development of carotid atherosclerosis. Moreover, these side differences might be severe for patients and worth further attention in clinical practice.

A spatiotemporal analysis of the left coronary artery biomechanics using fluid-structure interaction models

Biomechanics plays a critical role in coronary artery disease development. FSI simulation is commonly used to understand the hemodynamics and mechanical environment associated with atherosclerosis pathology. To provide a comprehensive characterization of patient-specific coronary biomechanics, an analysis of FSI simulation in the spatial and temporal domains was performed. In the current study, a three-dimensional FSI model of the LAD coronary artery was built based on a patient-specific geometry using COMSOL Multiphysics. The effect of myocardial bridging was simulated. Wall shear stress and its derivatives including time-averaged wall shear stress, wall shear stress gradient, and OSI were calculated across the cardiac cycle in multiple locations. Arterial wall strain (radial, circumferential, and longitudinal) and von Mises stress were calculated. To assess perfusion, vFFR was calculated. The results demonstrated the FSI model could identify regional and transient differences in biomechanical parameters within the coronary artery. The addition of myocardial bridging caused a notable change in von Mises stress and an increase in arterial strain during systole. The analysis performed in this manner takes greater advantage of the information provided in the space and time domains and can potentially assist clinical evaluation.

Distinct Temporal Pattern of the Prediction of Lumen Remodeling of Lower Extremity Vein Bypass Grafts by Initial Local Hemodynamics

We predicted human lower extremity vein bypass graft remodeling by hemodynamics. Computed tomography and duplex ultrasound scans of 55 patients were performed at 1 week and 1, 6, and 12 months post-implantation to obtain wall shear stress (WSS) and oscillatory shear index (OSI) at 1-mm intervals via computational fluid dynamics simulations. Graft remodeling was quantified by computed tomography-measured lumen diameter changes in the early (1 week-1 month), intermediate (1-6 months), and late (6-12 months) periods. Linear mixed-effect models were constructed to examine the overall relationship between remodeling and initial hemodynamics using the average data of all cross sections within the same graft. A significant association of graft remodeling with WSS (p < 0.001) and time (p = 0.001) was found; however, the effect size decreased with time (every 2.7 dyne/cm² increase of WSS was associated with a 0.39, 0.35, 0.002 mm diameter increase in the three periods, respectively). The association of remodeling with OSI was significant only in the intermediate period (every 0.1 increase of OSI was associated with a 0.25 mm lumen diameter decrease, p = 0.004). Therefore, the association of graft lumen remodeling with local hemodynamics has a distinct temporal pattern; WSS and OSI are predictive of remodeling only in certain postoperative periods.

Potential relationship between high wall shear stress and plaque rupture causing acute coronary syndrome

The relationship between high wall shear stress (WSS) and plaque rupture (PR) in longitudinal and circumferential locations remains uncertain. Overall, 100 acute coronary syndrome patients whose culprit lesions had PR, documented by optical coherence tomography (OCT), were enrolled. Lesion-specific three-dimensional coronary artery models were created using OCT data. WSS was computed with computational fluid dynamics analysis. PR was classified into upstream-PR, minimum lumen area-PR, and downstream-PR according to the PR's longitudinal location, and into central-PR and lateral-PR according to the disrupted fibrous cap circumferential location. In the longitudinal 3-mm segmental analysis, multivariate analysis demonstrated that higher WSS in the upstream segment was independently associated with upstream-PR, and thinner fibrous cap was independently associated with downstream-PR. In the PR cross-sections, the PR region had a significantly higher average WSS than non-PR region. In the cross-sectional analysis, the in-lesion peak WSS was frequently observed in the lateral (66.7%) and central regions (70%) in lateral-PR and central-PR, respectively. Multivariate analysis demonstrated that the presence of in-lesion peak WSS at the lateral region, thinner broken fibrous cap, and larger lumen area were independently associated with lateral-PR, while the presence of in-lesion peak WSS at the central region and thicker broken fibrous cap were independently associated with central-PR. In conclusion, OCT-based WSS simulation revealed that high WSS might be related to the longitudinal and circumferential locations of PR.

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

AIMS

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

METHODS AND RESULTS

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

CONCLUSION

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

The Need to Shift from Morphological to Structural Assessment for Carotid Plaque Vulnerability

Degree of luminal stenosis is generally considered to be an important indicator for judging the risk of atherosclerosis burden. However, patients with the same or similar degree of stenosis may have significant differences in plaque morphology and biomechanical factors. This study investigated three patients with carotid atherosclerosis within a similar range of stenosis. Using our developed fluid-structure interaction (FSI) modelling method, this study analyzed and compared the morphological and biomechanical parameters of the three patients. Although their degrees of carotid stenosis were similar, the plaque components showed a significant difference. The distribution range of time-averaged wall shear stress (TAWSS) of patient 2 was wider than that of patient 1 and patient 3. Patient 2 also had a much smaller plaque stress compared to the other two patients. There were significant differences in TAWSS and plaque stresses among three patients. This study suggests that plaque vulnerability is not determined by a single morphological factor, but rather by the combined structure. It is necessary to transform the morphological assessment into a structural assessment of the risk of plaque rupture.

Transcriptomic analysis of a 3D blood-brain barrier model exposed to disturbed fluid flow

Cerebral aneurysms are more likely to form at bifurcations in the vasculature, where disturbed fluid is prevalent due to flow separation at sufficiently high Reynolds numbers. While previous studies have demonstrated that altered shear stress exerted by disturbed flow disrupts endothelial tight junctions, less is known about how these flow regimes alter gene expression in endothelial cells lining the blood-brain barrier. Specifically, the effect of disturbed flow on expression of genes associated with cell-cell and cell-matrix interaction, which likely mediate aneurysm formation, remains unclear. RNA sequencing of immortalized cerebral endothelial cells isolated from the lumen of a 3D blood-brain barrier model reveals distinct transcriptional changes in vessels exposed to fully developed and disturbed flow profiles applied by both steady and physiological waveforms. Differential gene expression, validated by qRT-PCR and western blotting, reveals that lumican, a small leucine-rich proteoglycan, is the most significantly downregulated gene in endothelial cells exposed to steady, disturbed flow. Knocking down lumican expression reduces barrier function in the presence of steady, fully developed flow. Moreover, adding purified lumican into the hydrogel of the 3D blood-brain barrier model recovers barrier function in the region exposed to fully developed flow. Overall, these findings emphasize the importance of flow regimes exhibiting spatial and temporal heterogeneous shear stress profiles on cell-matrix interaction in endothelial cells lining the blood-brain barrier, while also identifying lumican as a contributor to the formation and maintenance of an intact barrier.

Differential hemodynamics between arteriovenous fistulas with or without intervention before successful use

A significant number of arteriovenous fistulas (AVFs) fail to maturate for dialysis. Although interventions promote maturation, functional primary patency loss is higher for AVFs with interventions (assisted maturation) than AVFs without interventions (un-assisted maturation). Although blood flow-associated hemodynamics have long been proposed to affect AVF remodeling, the optimal hemodynamic parameters for un-assisted maturation are unclear. Additionally, AVF maturation progress is generally not investigated until 6 weeks after AVF creation, and the examination is focused on the AVF's venous limb. In this exploratory study, patients (n = 6) underwent magnetic resonance imaging (MRI) at 1 day, 6 weeks, and 6 months after AVF creation surgery. Before successful use for hemodialysis, three AVFs required intervention and three did not. MRI of the AVFs were used to calculate lumen cross-sectional area (CSA) and perform computational fluid dynamics (CFD) to analyze hemodynamics, including velocity, wall shear stress (WSS), and vorticity. For the venous limb, the no-intervention group and intervention group had similar pre-surgery vein diameter and 1-day post-surgery venous CSA. However, the no-intervention group had statistically larger 1-day venous velocity (0.97 ± 0.67 m/s; mean ± SD), WSS (333 ± 336 dyne/cm2) and vorticity (1709 ± 1290 1/s) than the intervention group (velocity = 0.23 ± 0.10 m/s; WSS = 49 ± 40 dyne/cm2; vorticity = 493.1 ± 227 1/s) (P < 0.05). At 6 months, the no-intervention group had statistically larger venous CSA (43.5 ± 27.4 mm2) than the intervention group (15.1 ± 6.2 mm2) (P < 0.05). Regarding the arterial limb, no-intervention AVF arteries also had statistically larger 1-day velocity (1.17 ± 1.0 m/s), WSS (340 ± 423 dyne/cm2), vorticity (1787 ± 1694 1/s), and 6-month CSA (22.6 ± 22.7 mm2) than the intervention group (velocity = 0.64 ± 0.36 m/s; WSS = 104 ± 116 dyne/cm2, P < 0.05; vorticity = 867 ± 4551/s; CSA = 10.7 ± 6.0 mm2, P < 0.05). Larger venous velocity, WSS, and vorticity immediately after AVF creation surgery may be important for later lumen enlargement and AVF maturation, with the potential to be used as a tool to help diagnose poor AVF maturation earlier. However, future studies using a larger cohort are needed to validate this finding and determine cut off values, if any.

Atherogenic potential of microgravity hemodynamics in the carotid bifurcation: a numerical investigation

Long-duration spaceflight poses multiple hazards to human health, including physiological changes associated with microgravity. The hemodynamic adaptations occurring upon entry into weightlessness have been associated with retrograde stagnant flow conditions and thromboembolic events in the venous vasculature but the impact of microgravity on cerebral arterial hemodynamics and function remains poorly understood. The objective of this study was to quantify the effects of microgravity on hemodynamics and wall shear stress (WSS) characteristics in 16 carotid bifurcation geometries reconstructed from ultrasonography images using computational fluid dynamics modeling. Microgravity resulted in a significant 21% increase in flow stasis index, a 22-23% decrease in WSS magnitude and a 16-26% increase in relative residence time in all bifurcation branches, while preserving WSS unidirectionality. In two anatomies, however, microgravity not only promoted flow stasis but also subjected the convex region of the external carotid arterial wall to a moderate increase in WSS bidirectionality, which contrasted with the population average trend. This study suggests that long-term exposure to microgravity has the potential to subject the vasculature to atheroprone hemodynamics and this effect is modulated by subject-specific anatomical features. The exploration of the biological impact of those microgravity-induced WSS aberrations is needed to better define the risk posed by long spaceflights on cardiovascular health.

An Oscillatory Shear Index-Based Model to Describe Progressive Carotid Artery Stenosis

Background and aims: This study describes and demonstrates the applicability of a novel in silico method for modeling progressive carotid artery stenosis using the oscillatory shear index (OSI) as the basis of stenosis. Methods: Three-dimensional reconstructions of 11 carotid arteries were generated using patient-derived magnetic resonance angiography and duplex ultrasound data. Computational fluid dynamic simulations were sequentially generated following computational stenosis assessment, and corresponding changes in OSI were observed and used as measure of morphological stabilization. Results: 6 carotid models showed progressive stenosis with statistically significant increases in regions of high OSI (OSI >.2, P < .05) with eventual carotid occlusion in 1 of the cases. Three models remained free or nearly free of increased OSI, whereas 1 model showed an overall decrease in high OSI regions (P < .05) and another trended in that direction but did not achieve statistical significance (P = .145). Conclusions: To our knowledge, this is the first computational model describing progressive stenosis in any peripheral artery including the carotid. Taken together, this study provides a novel framework for computational hemodynamic investigations on progressive atherosclerosis in the carotid artery.

A New Method for Quantifying Abdominal Aortic Wall Shear Stress Using Phase Contrast Magnetic Resonance Imaging and the Womersley Solution

Wall shear stress (WSS) is an important mediator of cardiovascular pathologies and there is a need for its reliable evaluation as a potential prognostic indicator. The purpose of this work was to develop a method that quantifies WSS from two-dimensional (2D) phase contrast magnetic resonance (PCMR) imaging derived flow waveforms, apply this method to PCMR data acquired in the abdominal aorta of healthy volunteers, and to compare PCMR-derived WSS values to values predicted from a computational fluid dynamics (CFD) simulation. The method uses PCMR-derived flow versus time waveforms constrained by the Womersley solution for pulsatile flow in a cylindrical tube. The method was evaluated for sensitivity to input parameters, intrastudy repeatability and was compared with results from a patient-specific CFD simulation. 2D-PCMR data were acquired in the aortas of healthy men (n = 12) and women (n = 15) and time-averaged WSS (TAWSS) was compared. Agreement was observed when comparing TAWSS between CFD and the PCMR flow-based method with a correlation coefficient of 0.88 (CFD: 15.0 ± 1.9 versus MRI: 13.5 ± 2.4 dyn/cm2) though comparison of WSS values between the PCMR-based method and CFD predictions indicate that the PCMR method underestimated instantaneous WSS by 3.7 ± 7.6 dyn/cm2. We found no significant difference in TAWSS magnitude between the sexes; 8.19 ± 2.25 versus 8.07 ± 1.71 dyn/cm2, p = 0.16 for men and women, respectively.

Probing the possibility of lesion formation/progression in vicinity of a primary atherosclerotic plaque: A fluid-solid interaction study and angiographic evidences

It is shown that certain locations in the arterial tree, such as coronary and cerebral arteries, are more prevalent to plaque formation. Endothelial activation and consequent plaque development are attributed to local hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI), relative residence time (RRT), and stress phase angle. After a certain level of plaque progression, these hemodynamic parameters are disturbed before and after the plaque. In the current study, it is hypothesized that the vicinity of a primary lesion is susceptible for further degeneration and second plaque formation. A fluid-solid interaction (FSI) model of the coronary artery with different levels of asymmetric constriction, is simulated and the trend of hemodynamic parameters were studied in both of the plaque side (PS) and the opposite wall (facing the plaque [PF]). Also, a novel factor is introduced that can identify the high-risk regions associated with WSS oscillations to negative values. Our results indicate that when more than half of the artery is constricted, the downstream of the plaque is highly exposed to endothelial pathogenesis the PS, such that negative WSS, and as well, critical values of OSI and RRT, that is, -1.2 Pa, 0.42 and 6.5 s, respectively arise in this region. PS endothelial cells in this region exposed to the highest risk of atherosclerosis based on the proposed index (3 out of 3). As well, three cases of angiographic images are provided that confirms existence of secondary lesion close to the primary one as predicted by our computational simulations.

High intensity exercise and passive hot water immersion cause similar post intervention changes in peripheral and cerebral shear

Passive hot water immersion (PHWI) provides a peripheral vasculature shear stimulus comparable to low intensity exercise within the active skeletal muscle, whereas moderate and high intensity exercise elicit substantially greater shear rates in the peripheral vasculature, likely conferring greater vascular benefits. Few studies have compared post intervention shear rates in the peripheral and cerebral vasculature following high intensity exercise and PHWI, especially considering that the post intervention recovery period represents a key window in which adaptation occurs. Therefore, we aimed to compare shear rates in the internal carotid artery (ICA), vertebral artery (VA) and common femoral artery (CFA) between high intensity exercise and PHWI for up to 80 minutes post intervention. Fifteen healthy (27 ± 4 years), moderately trained individuals underwent three-time matched interventions in a randomised order which included 30 minutes of whole-body immersion in a 42°C hot bath, 30 minutes of treadmill running and 5x4 minute high intensity intervals (HIIE). There were no differences in ICA (P= 0.4643) and VA (P=0.1940) shear rates between PHWI and exercise (both continuous and HIIE) post intervention. All three interventions elicited comparable increases in CFA shear rate post intervention (P=0.0671), however, CFA shear rate was slightly higher 40 minutes post threshold running (P=0.0464) and, slightly higher, although not statically for HIIE (P=0.0565) compared with PHWI. Our results suggest that time and core temperature matched high intensity exercise and PHWI elicit limited changes in cerebral shear and comparable increases in peripheral vasculature shear rates when measured for up to 80 minutes post intervention.

An automated software for real-time quantification of wall shear stress distribution in quantitative coronary angiography data

BACKGROUND

Wall shear stress (WSS) estimated in 3D-quantitative coronary angiography (QCA) models appears to provide useful prognostic information and identifies high-risk patients and lesions. However, conventional computational fluid dynamics (CFD) analysis is cumbersome limiting its application in the clinical arena. This report introduces a user-friendly software that allows real-time WSS computation and examines its reproducibility and accuracy in assessing WSS distribution against conventional CFD analysis.

METHODS

From a registry of 414 patients with borderline negative fractional flow reserve (0.81-0.85), 100 lesions were randomly selected. 3D-QCA and CFD analysis were performed using the conventional approach and the novel CAAS Workstation WSS software, and QCA as well as WSS estimations of the two approaches were compared. The reproducibility of the two methodologies was evaluated in a subgroup of 50 lesions.

RESULTS

A good agreement was noted between the conventional approach and the novel software for 3D-QCA metrics (ICC range: 0.73-0-93) and maximum WSS at the lesion site (ICC: 0.88). Both methodologies had a high reproducibility in assessing lesion severity (ICC range: 0.83-0.97 for the conventional approach; 0.84-0.96 for the CAAS Workstation WSS software) and WSS distribution (ICC: 0.85-0.89 and 0.83-0.87, respectively). Simulation time was significantly shorter using the CAAS Workstation WSS software compared to the conventional approach (4.13 ± 0.59 min vs 23.14 ± 2.56 min, p < 0.001).

CONCLUSION

CAAS Workstation WSS software is fast, reproducible, and accurate in assessing WSS distribution. Therefore, this software is expected to enable the broad use of WSS metrics in the clinical arena to identify high-risk lesions and vulnerable patients.