Nitric oxide transport in an axisymmetric stenosis

Abnormal flow pattern of low wall shear stress and high oscillatory shear index in spontaneous vertebral artery dissection with vertebral artery hypoplasia

Introduction

Spontaneous vertebral artery dissection (sVAD) might tend to develop in vertebral artery hypoplasia (VAH) with hemodynamic dysfunction and it is crucial to assess hemodynamics in sVAD with VAH to investigate this hypothesis. This retrospective study aimed to quantify hemodynamic parameters in patients with sVAD with VAH.

Methods

Patients who had suffered ischemic stroke due to an sVAD of VAH were enrolled in this retrospective study. The geometries of 14 patients (28 vessels) were reconstructed using Mimics and Geomagic Studio software from CT angiography (CTA). ANSYS ICEM and ANSYS FLUENT were utilized for mesh generation, set boundary conditions, solve governing equations, and perform numerical simulations. Slices were obtained at the upstream area, dissection or midstream area and downstream area of each VA. The blood flow patterns were visualized through instantaneous streamline and pressure at peak systole and late diastole. The hemodynamic parameters included pressure, velocity, time-averaged blood flow, time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), endothelial cell action potential (ECAP), relative residence time (RRT) and time-averaged nitric oxide production rate (TAR_NO).

Results

Significant focal increased velocity was present in the dissection area of steno-occlusive sVAD with VAH compared to other nondissected areas (0.910 m/s vs. 0.449 vs. 0.566, p < 0.001), while focal slow flow velocity was observed in the dissection area of aneurysmal dilatative sVAD with VAH according to velocity streamlines. Steno-occlusive sVAD with VAH arteries had a lower time-averaged blood flow (0.499 cm³/s vs. 2.268, p < 0.001), lower TAWSS (1.115 Pa vs. 2.437, p = 0.001), higher OSI (0.248 vs. 0.173, p = 0.006), higher ECAP (0.328 Pa^-1 vs. 0.094, p = 0.002), higher RRT (3.519 Pa^-1 vs. 1.044, p = 0.001) and deceased TAR_NO (104.014 nM/s vs. 158.195, p < 0.001) than the contralateral VAs.

Conclusion

Steno-occlusive sVAD with VAH patients had abnormal blood flow patterns of focal increased velocity, low time-averaged blood flow, low TAWSS, high OSI, high ECAP, high RRT and decreased TAR_NO. These results provide a good basis for further investigation of sVAD hemodynamics and support the applicability of the CFD method in testing the hemodynamic hypothesis of sVAD. More detailed hemodynamic conditions with different stages of sVAD are warranted in the future.

Risk evaluation of adverse aortic events in patients with non-circular aortic annulus after transcatheter aortic valve implantation: a numerical study

Transcatheter aortic valve implantation (TAVI) is a micro-invasive surgery used to treat patients with aortic stenosis (AS) efficiently. However, the uneven valve expansion can cause a non-circular annulus, which is one of the main factors leading to complications after TAVI. As a preliminary work, the main purpose of this study was to evaluate the risk of adverse aortic events in patients with a non-circular aortic annulus after TAVI. This study numerically investigated the distribution of four wall shear stress (WSS)-based indicators and three helicity-based indicators in eight patient-specific aortas with different annulus including circular, type I elliptical and type II elliptical shapes. Both elliptical annulus features can significantly enhance the intensity of the helicity (h2) in the ascending aorta (p < 0.001). However, for the type I elliptical annulus, the spiral flow structure was changed into low-velocity and disturbed flow pattern close to the inner side of the aortic arch. For the type II elliptical annulus, the spiral flow remained but became skewed in distribution. The elliptical annulus feature could increase the general level WSS-based indicators, especially in the ascending aorta. However, due to the disturbance of spiral flow or second helical flow in ascending aortas, areas with low TAWSS accompanied by high oscillatory shear index (OSI) and cross flow index (CFI) were observed in all the ascending aortas with non-circular annulus. The elliptical annulus feature can change the hemodynamic environment in the aortic arch, especially in the ascending aorta. Although both elliptical annulus features enhanced the strength of helicity, the uniform distribution of the helical flow was disturbed, especially in the ascending aorta, indicating the potential risk of adverse aortic events may increase. Therefore, for the patients without paravalvular leak but elliptical annulus shape after TAVI treatment, surgeons may be needed to consider further dilatation to make the non-circular annulus become circular.

'Heart-like' cross-sectional shape can better improve the hemodynamics in spiral laminar flow graft for small-caliber bypass application: a numerical study

Small-caliber grafts remain disappointed in the long-term bypass surgeries of coronary and peripheral arterial diseases. In order to improve the hemodynamics in small-caliber artery bypass grafts (ABGs), an improved spiral laminar flow (improved-SLF) graft with a 'heart-like' cross-sectional shape was proposed and verified by computational fluid dynamics simulation in this study. The results show that such graft can indeed induce a spiral flow and enhance the WSS distribution on the graft section. Furthermore, the helically distributed ribbon of unfavorable WSS observed in the original SLF graft was eliminated in the improved-SLF graft due to its smoothed and gentle helical ridge. On the other hand, improved-SLF ABG improved the WSS distribution in the distal anastomosis as well, because it maintained the strength of spiral flow when entering the anastomosis region. Finally, the improved-SLF ABG slightly increased the pressure drop along the bypass due to its small change of the general graft structure. As a proof-of-concept study, it can be concluded that improved-SLF graft can not only evenly enhance the WSS distribution in the graft section, but also improve the hemodynamic environment in the distal anastomosis without significantly increasing the pressure drop along the bypass, indicating such new helical-type graft may be more suitable to be used in the small-caliber graft bypass surgeries.

Delivery of Nitric Oxide in the Cardiovascular System: Implications for Clinical Diagnosis and Therapy

Nitric oxide (NO) is a key molecule in cardiovascular homeostasis and its abnormal delivery is highly associated with the occurrence and development of cardiovascular disease (CVD). The assessment and manipulation of NO delivery is crucial to the diagnosis and therapy of CVD, such as endothelial dysfunction, atherosclerotic progression, pulmonary hypertension, and cardiovascular manifestations of coronavirus (COVID-19). However, due to the low concentration and fast reaction characteristics of NO in the cardiovascular system, clinical applications centered on NO delivery are challenging. In this tutorial review, we first summarized the methods to estimate the in vivo NO delivery process, based on computational modeling and flow-mediated dilation, to assess endothelial function and vulnerability of atherosclerotic plaque. Then, emerging bioimaging technologies that have the potential to experimentally measure arterial NO concentration were discussed, including Raman spectroscopy and electrochemical sensors. In addition to diagnostic methods, therapies aimed at controlling NO delivery to regulate CVD were reviewed, including the NO release platform to treat endothelial dysfunction and atherosclerosis and inhaled NO therapy to treat pulmonary hypertension and COVID-19. Two potential methods to improve the effectiveness of existing NO therapy were also discussed, including the combination of NO release platform and computational modeling, and stem cell therapy, which currently remains at the laboratory stage but has clinical potential for the treatment of CVD.

Flow shear stress controls the initiation of neovascularization via heparan sulfate proteoglycans within a biomimetic microfluidic model

Endothelial cells (ECs) in vivo are subjected to three forms of shear stress induced by luminal blood flow, transendothelial flow and interstitial flow simultaneously. It is controversial that shear stress, especially the component induced by luminal flow, was thought to inhibit the initialization of angiogenesis and trigger arteriogenesis. Here, we combined microfabrication techniques and delicate numerical simulations to reconstruct the initial physiological microenvironment of neovascularization in vitro, where ECs experience high luminal shear stress, physiological transendothelial flow and various vascular endothelial growth factor (VEGF) distributions simultaneously. With the biomimetic microfluidic model, cell alignment and endothelial sprouting assays were carried out. We found that luminal shear stress inhibits endothelial sprouting and tubule formation in a dose-dependent manner. Although a high concentration of VEGF increases EC sprouting, neither a positive nor a negative VEGF gradient additionally affects the degree of sprouting, and luminal shear stress significantly attenuates neovascularization even in the presence of VEGF. Heparinase was used to selectively degrade the heparan sulfate proteoglycan (HSPG) coating on ECs and messenger RNA profiles in ECs were analyzed. It turned out that HSPGs could act as a mechanosensor to sense the change of fluid shear stress, modulate multiple EC gene expressions, and hence affect neovascularization. In summary, distraction from the stabilized state, such as decreased luminal shear stress, increased VEGF and the destructed mechanotransduction of HSPGs would induce the initiation of neovascularization. Our study highlights the key role of the magnitude and forms of shear stress in neovascularization.

The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction

Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.

Non-Newtonian Effects on Patient-Specific Modeling of Fontan Hemodynamics

The Fontan procedure is a common palliative surgery for congenital single ventricle patients. In silico and in vitro patient-specific modeling approaches are widely utilized to investigate potential improvements of Fontan hemodynamics that are related to long-term complications. However, there is a lack of consensus regarding the use of non-Newtonian rheology, warranting a systematic investigation. This study conducted in silico patient-specific modeling for twelve Fontan patients, using a Newtonian and a non-Newtonian model for each patient. Differences were quantified by examining clinically relevant metrics: indexed power loss (iPL), indexed viscous dissipation rate (iVDR), hepatic flow distribution (HFD), and regions of low wall shear stress (A_WSS). Four sets of "non-Newtonian importance factors" were calculated to explore their effectiveness in identifying the non-Newtonian effect. No statistical differences were observed in iPL, iVDR, and HFD between the two models at the population-level, but large inter-patient variations exist. Significant differences were detected regarding A_WSS, and its correlations with non-Newtonian importance factors were discussed. Additionally, simulations using the non-Newtonian model were computationally faster than those using the Newtonian model. These findings distinguish good importance factors for identifying non-Newtonian rheology and encourage the use of a non-Newtonian model to assess Fontan hemodynamics.

EFFECT OF STENOSIS SEVERITY ON SHEAR-INDUCED DIFFUSION OF RED BLOOD CELLS IN CORONARY ARTERIES

In large blood vessels, migration of red blood cells (RBCs) affects the concentration of platelets and the transport of oxygen to the arterial endothelial cells. In this work, we investigate the locations where hydrodynamic diffusion of RBCs occurs and the effects of stenosis severity on shear-induced diffusion (SID) of RBCs, concentration distribution and wall shear stress (WSS). For the first time, multiphase mixture theory approach with Phillips shear-induced diffusive flux model coupled with Quemada non-Newtonian viscosity model has been applied to numerically simulate the RBCs macroscopic behavior in four different degrees of stenosis (DOS) geometries, viz., 30%, 50%, 70% and 85%. Considering SID of RBCs, the calculated average WSS increased by 77.70% which emphasises the importance of SID in predicting hemodynamic parameters. At the stenosis throat, it was observed that 85% DOS model had the lowest concentration of RBCs near the wall and highest concentration at the center. For the stenosis models with 70% and 85% DOS, the RBC lumen wall concentration at the distal section of stenosis becomes inhomogeneous with the maximum fluctuation of 1.568%. Finally, the wall regions with low WSS and low RBC concentrations correlate well with the atherosclerosis sites observed clinically.

Numerical simulation of haemodynamics and low-density lipoprotein transport in the rabbit aorta and their correlation with atherosclerotic plaque thickness

Two mechanisms of shear stress and mass transport have been recognized to play an important role in the development of localized atherosclerosis. However, their relationship and roles in atherogenesis are still obscure. It is necessary to investigate quantitatively the correlation among low-density lipoproteins (LDL) transport, haemodynamic parameters and plaque thickness. We simulated blood flow and LDL transport in rabbit aorta using computational fluid dynamics and evaluated plaque thickness in the aorta of a high-fat-diet rabbit. The numerical results show that regions with high luminal LDL concentration tend to have severely negative haemodynamic environments (HEs). However, for regions with moderately and slightly high luminal LDL concentration, the relationship between LDL concentration and the above haemodynamic indicators is not clear cut. Point-by-point correlation with experimental results indicates that severe atherosclerotic plaque corresponds to high LDL concentration and seriously negative HEs, less severe atherosclerotic plaque is related to either moderately high LDL concentration or moderately negative HEs, and there is almost no atherosclerotic plaque in regions with both low LDL concentration and positive HEs. In conclusion, LDL distribution is closely linked to blood flow transport, and the synergetic effects of luminal surface LDL concentration and wall shear stress-based haemodynamic indicators may determine plaque thickness.

Synergy between shear-induced migration and secondary flows on red blood cells transport in arteries: considerations on oxygen transport

Shear-induced migration of red blood cells (RBCs) is a well-known phenomenon characterizing blood flow in the small vessels (micrometre to millimetre size) of the cardiovascular system. In large vessels, like the abdominal aorta and the carotid artery (millimetre to centimetre size), the extent of this migration and its interaction with secondary flows has not been fully elucidated. RBC migration exerts its influence primarily on platelet concentration, oxygen transport and oxygen availability at the luminal surface, which could influence vessel wall disease processes in and adjacent to the intima. Phillips' shear-induced particle migration model, coupled to the Quemada viscosity model, was employed to simulate the macroscopic behaviour of RBCs in four patient-specific geometries: a normal abdominal aorta, an abdominal aortic aneurysm (AAA), a normal carotid bifurcation and a stenotic carotid bifurcation. Simulations show a migration of RBCs from the near-wall region with a lowering of wall haematocrit (volume fraction of RBCs) on the posterior side of the normal aorta and on the lateral-external side of the iliac arteries. A marked migration is observed on the outer wall of the carotid sinus, along the common carotid artery and in the carotid stenosis. No significant migration is observed in the AAA. The spatial and temporal patterns of wall haematocrit are correlated with the near-wall shear layer and with the secondary flows induced by the vessel curvature. In particular, secondary flows accentuate the initial lowering in RBC near-wall concentration by convecting RBCs from the inner curvature side to the outer curvature side. The results reinforce data in literature showing a decrease in oxygen partial pressure on the inner curvature wall of the carotid sinus induced by the presence of secondary flows. The lowering of wall haematocrit is postulated to induce a decrease in oxygen availability at the luminal surface through a diminished concentration of oxyhaemoglobin, hence contributing, with the reported lowered oxygen partial pressure, to local hypoxia.

Nitric oxide transport in normal human thoracic aorta: effects of hemodynamics and nitric oxide scavengers

Despite the crucial role of nitric oxide (NO) in the homeostasis of the vasculature, little quantitative information exists concerning NO transport and distribution in medium and large-sized arteries where atherosclerosis and aneurysm occur and hemodynamics is complex. We hypothesized that local hemodynamics in arteries may govern NO transport and affect the distribution of NO in the arteries, hence playing an important role in the localization of vascular diseases. To substantiate this hypothesis, we presented a lumen/wall model of the human aorta based on its MRI images to simulate the production, transport and consumption of NO in the arterial lumen and within the aortic wall. The results demonstrated that the distribution of NO in the aorta was quite uneven with remarkably reduced NO bioavailability in regions of disturbed flow, and local hemodynamics could affect NO distribution mainly via flow dependent NO production rate of endothelium. In addition, erythrocytes in the blood could moderately modulate NO concentration in the aorta, especially at the endothelial surface. However, the reaction of NO within the wall could only slightly affect NO concentration on the luminal surface, but strongly reduce NO concentration within the aortic wall. A strong positive correlation was revealed between wall shear stress and NO concentration, which was affected by local hemodynamics and NO reaction rate. In conclusion, the distribution of NO in the aorta may be determined by local hemodynamics and modulated differently by NO scavengers in the lumen and within the wall.

Simulation of contrast agent transport in arteries with multilayer arterial wall: impact of arterial transmural transport on the bolus delay and dispersion

One assumption of DSC-MRI is that the injected contrast agent is kept totally intravascular and the arterial wall is impermeable to contrast agent. The assumption is unreal for such small contrast agent as Gd-DTPA can leak into the arterial wall. To investigate whether the unreal assumption is valid for the estimation of the delay and dispersion of the contrast agent bolus, we simulated flow and Gd-DTPA transport in a model with multilayer arterial wall and analyzed the bolus delay and dispersion qualified by mean vascular transit time (MVTT) and the variance of the vascular transport function. Factors that may affect Gd-DTPA transport hence the delay and dispersion were further investigated, such as integrity of endothelium and disturbed flow. The results revealed that arterial transmural transport would slightly affect MVTT and moderately increase the variance. In addition, although the integrity of endothelium can significantly affect the accumulation of contrast agent in the arterial wall, it had small effects on the bolus delay and dispersion. However, the disturbed flow would significantly increase both MVTT and the variance. In conclusion, arterial transmural transport may have a small effect on the bolus delay and dispersion when compared to the flow pattern in the artery.