Wall shear stresses remain elevated in mature arteriovenous fistulas: a case study

Maintaining vascular access (VA) patency continues to be the greatest challenge for dialysis patients. VA dysfunction, primarily due to venous neointimal hyperplasia development and stenotic lesion formation, is mainly attributed to complex hemodynamics within the arteriovenous fistula (AVF). The effect of VA creation and the subsequent geometrical remodeling on the hemodynamics and shear forces within a mature patient-specific AVF is investigated. A 3D reconstructed geometry of a healthy vein and a fully mature patient-specific AVF was developed from a series of 2D magnetic resonance image scans. A previously validated thresholding technique for region segmentation and lumen cross section contour creation was conducted in MIMICS 10.01, allowing for the creation of a 3D reconstructed geometry. The healthy vein and AVF computational models were built, subdivided, and meshed in GAMBIT 2.3. The computational fluid dynamic (CFD) code FLUENT 6.3.2 (Fluent Inc., Lebanon, NH) was employed as the finite volume solver to determine the hemodynamics and shear forces within the healthy vein and patient-specific AVF. Geometrical alterations were evaluated and a CFD analysis was conducted. Substantial geometrical remodeling was observed, following VA creation with an increase in cross-sectional area, out of plane curvature (maximum angle of curvature in AVF=30?deg), and angle of blood flow entry. The mean flow velocity entering the vein of the AVF is dramatically increased. These factors result in complex three-dimensional hemodynamics within VA junction (VAJ) and efferent vein of the AVF. Complex flow patterns were observed and the maximum and mean wall shear stress (WSS) magnitudes are significantly elevated. Flow reversal was found within the VAJ and efferent vein. Extensive geometrical remodeling during AVF maturation does not restore physiological hemodynamics to the VAJ and venous conduit of the AVF, and high WSS and WSS gradients, and flow reversal persist. It is theorized that the vessel remodelling and the continued non-physiological hemodynamics within the AVF compound to result in stenotic lesion development.     

Angiogenesis and vascular remodelling in normal and cancerous tissues

Vascular development and homeostasis are underpinned by two fundamental features: the generation of new vessels to meet the metabolic demands of under-perfused regions and the elimination of vessels that do not sustain flow. In this paper we develop the first multiscale model of vascular tissue growth that combines blood flow, angiogenesis, vascular remodelling and the subcellular and tissue scale dynamics of multiple cell populations. Simulations show that vessel pruning, due to low wall shear stress, is highly sensitive to the pressure drop across a vascular network, the degree of pruning increasing as the pressure drop increases. In the model, low tissue oxygen levels alter the internal dynamics of normal cells, causing them to release vascular endothelial growth factor (VEGF), which stimulates angiogenic sprouting. Consequently, the level of blood oxygenation regulates the extent of angiogenesis, with higher oxygenation leading to fewer vessels. Simulations show that network remodelling (and de novo network formation) is best achieved via an appropriate balance between pruning and angiogenesis. An important factor is the strength of endothelial tip cell chemotaxis in response to VEGF. When a cluster of tumour cells is introduced into normal tissue, as the tumour grows hypoxic regions form, producing high levels of VEGF that stimulate angiogenesis and cause the vascular density to exceed that for normal tissue. If the original vessel network is sufficiently sparse then the tumour may remain localised near its parent vessel until new vessels bridge the gap to an adjacent vessel. This can lead to metastable periods, during which the tumour burden is approximately constant, followed by periods of rapid growth.     

Radial artery wall shear stress evaluation in patients with arteriovenous fistula for hemodialysis access

It has been extensively documented that changes in blood flow induce vascular remodeling and this phenomenon seems to be correlated to the shear forces imposed on the vessel wall by motion of blood. Wall shear stress, the tractive force that acts on the endothelium, has been shown to influence endothelial cell function. To study changes in wall shear stress that develop on the vessel wall upon changes of blood flow, we set up a technique that allows estimation of shear stress in the radial artery of patients on chronic hemodialysis therapy. The technique is based on color-flow Doppler examination of the radial artery before and after surgical creation of radiocephalic fistula for hemodialysis. Calculation of time function wall shear stress and blood flow rate in the radial artery is performed on the basis of arterial diameter, center-line velocity waveform and blood viscosity, using a numerical method developed according to Womersley's theory for pulsatile flow in tubes. The results presented confirm that the model developed is suitable for calculation of the wall shear stress that develops in the radial artery of patients before and after surgical creation of an arteriovenous fistula for hemodialysis. This methodology was developed for characterization of wall shear stress in the radial artery but may be well applied to other vessels that can be examined by echo-Doppler technique.     

Branching exponent heterogeneity and wall shear stress distribution in vascular trees

A bifurcating arterial system with Poiseuille flow can function at minimum cost and with uniform wall shear stress if the branching exponent (z) = 3 [where z is defined by (D(1))(z) = (D(2))(z) + (D(3))(z); D(1) is the parent vessel diameter and D(2) and D(3) are the two daughter vessel diameters at a bifurcation]. Because wall shear stress is a physiologically transducible force, shear stress-dependent control over vessel diameter would appear to provide a means for preserving this optimal structure through maintenance of uniform shear stress. A mean z of 3 has been considered confirmation of such a control mechanism. The objective of the present study was to evaluate the consequences of a heterogeneous distribution of z values about the mean with regard to this uniform shear stress hypothesis. Simulations were carried out on model structures otherwise conforming to the criteria consistent with uniform shear stress when z = 3 but with varying distributions of z. The result was that when there was significant heterogeneity in z approaching that found in a real arterial tree, the coefficient of variation in shear stress was comparable to the coefficient of variation in z and nearly independent of the mean value of z. A systematic increase in mean shear stress with decreasing vessel diameter was one component of the variation in shear stress even when the mean z = 3. The conclusion is that the influence of shear stress in determining vessel diameters is not, per se, manifested in a mean value of z. In a vascular tree having a heterogeneous distribution in z values, a particular mean value of z (e.g., z = 3) apparently has little bearing on the uniform shear stress hypothesis.     

Flow and pressure distributions in vascular networks consisting of distensible vessels

We examine the influence of vessel distensibility on the fraction of the total network flow passing through each vessel of a model vascular network. An exact computational methodology is developed yielding an analytic proof. For a class of structurally heterogeneous asymmetric vascular networks, if all the individual vessels share a common distensibility relation when the total network flow is changed, this methodology proves that each vessel will continue to receive the same fraction of the total network flow. This constant flow partitioning occurs despite a redistribution of pressures, which may result in a decrease in the diameter of one and an increase in the diameter of the other of two vessels having a common diameter at a common pressure. This theoretical observation, taken along with published experimental observations on pulmonary vessel distensibilities, suggests that vessel diameter-independent distensibility in the pulmonary vasculature may be an evolutionary adaptation for preserving the spatial distribution of pulmonary blood flow in the face of large variations in cardiac output.     

Shear stress,vascular remodeling and neointimal formation

The role of shear stress in atherosclerosis has been well documented. However, its role in restenosis was underexposed. In this paper a novel in vivo measuring technique and several of its applications related to restenosis will be described. The technique consists of a combination of 3D reconstruction of blood vessels and computational fluid dynamics (CFD). The 3D imaging techniques use either of 3D intravascular ultrasound (IVUS) as a stand-alone technique or a fusion of biplane angiography and IVUS (ANGUS). CFD is applied in order to relate local shear stress distribution to the morphology of the vessel wall. In the applications of these techniques it will be demonstrated that shear stress plays a role in the prediction of neointimal formation in in-stent restenosis and in vascular remodeling after balloon angioplasty. Attempts to locally increase shear stress by a newly developed flow divider indicate that shear stress reduce in-stent neointimal formation by 50%.     

Mechanisms of blood flow-induced vascular enlargement

Chronic changes in wall shear stress lead to vascular remodeling, characterized by increased vascular wall diameter and thickness, to restore wall shear stress values to baseline. Release of nitric oxide from endothelial cells exposed to excessive shear is a fundamental step in the remodeling process, and potentially triggers a cascade of events, including growth factor induction and matrix metalloproteinase activation, that together contribute to restructuralization of the vessel wall. Understanding these processes could help explain how changes in blood vessel wall structure occur in the context of atherosclerosis or aortic aneurisms.     

A mathematical model of maternal vascular growth and remodeling and changes in maternal hemodynamics in uncomplicated pregnancy

The maternal vasculature undergoes tremendous growth and remodeling (G&R) that enables a?>?15-fold increase in blood flow through the uterine vasculature from conception to term. Hemodynamic metrics (e.g., uterine artery pulsatility index, UA-PI) are useful for the prognosis of pregnancy complications; however, improved characterization of the maternal hemodynamics is necessary to improve prognosis. The goal of this paper is to develop a mathematical framework to characterize maternal vascular G&R and hemodynamics in uncomplicated human pregnancies. A validated 1D model of the human vascular tree from the literature was adapted and inlet blood flow waveforms at the ascending aorta at 4 week increments from 0 to 40 weeks of gestation were prescribed. Peripheral resistances of each terminal vessel were adjusted to achieve target flow rates and mean arterial pressure at each gestational age. Vessel growth was governed by wall shear stress (and axial lengthening in uterine vessels), and changes in vessel distensibility were related to vessel growth. Uterine artery velocity waveforms generated from this model closely resembled ultrasound results from the literature. The literature UA-PI values changed significantly across gestation, increasing in the first month of gestation, then dramatically decreasing from 4 to 20 weeks. Our results captured well the time-course of vessel geometry, material properties, and UA-PI. This 1D fluid-G&R model captured the salient hemodynamic features across a broad range of clinical reports and across gestation for uncomplicated human pregnancy. While results capture available data well, this study highlights significant gaps in available data required to better understand vascular remodeling in pregnancy.

     

Circumferential vascular deformation after stent implantation alters wall shear stress evaluated with time-dependent 3D computational fluid dynamics models.

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that stent geometry may cause local alterations in wall shear stress (WSS) that have been associated with neointimal hyperplasia and subsequent restenosis. However, previous CFD studies have ignored histological evidence of vascular straightening between circumferential stent struts. We tested the hypothesis that consideration of stent-induced vascular deformation may more accurately predict alterations in indexes of WSS that may subsequently account for histological findings after stenting. We further tested the hypothesis that the severity of these alterations in WSS varies with the degree of vascular deformation after implantation. Steady-state and time-dependent simulations of three-dimensional CFD arteries based on canine coronary artery measurements of diameter and blood flow were conducted, and WSS and WSS gradients were calculated. Circumferential straightening introduced areas of high WSS between stent struts that were absent in stented vessels of circular cross section. The area of vessel exposed to low WSS was dependent on the degree of circumferential vascular deformation and axial location within the stent. Stents with four vs. eight struts increased the intrastrut area of low WSS in vessels, regardless of cross-sectional geometry. Elevated WSS gradients were also observed between struts in vessels with polygonal cross sections. The results obtained using three-dimensional CFD models suggest that changes in vascular geometry after stent implantation are important determinants of WSS distributions that may be associated with subsequent neointimal hyperplasia.     

A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations

We present a novel framework for investigating the role of vascular structure on arterial haemodynamics in large vessels, with a special focus on the human common carotid artery (CCA). The analysis is carried out by adopting a three-dimensional (3D) derived, fibre-reinforced, hyperelastic structural model, which is coupled with an axisymmetric, reduced order model describing blood flow. The vessel transmural pressure and lumen area are related via a Holzapfel–Ogden type of law, and the residual stresses along the thickness and length of the vessel are also accounted for. After a structural characterization of the adopted hyperelastic model, we investigate the link underlying the vascular wall response and blood-flow dynamics by comparing the proposed framework results against a popular tube law. The comparison shows that the behaviour of the model can be captured by the simpler linear surrogate only if a representative value of compliance is applied. Sobol’s multi-variable sensitivity analysis is then carried out in order to identify the extent to which the structural parameters have an impact on the CCA haemodynamics. In this case, the local pulse wave velocity (PWV) is used as index for representing the arterial transmission capacity of blood pressure waveforms. The sensitivity analysis suggests that some geometrical factors, such as the stress-free inner radius and opening angle, play a major role on the system’s haemodynamics. Subsequently, we quantified the differences in haemodynamic variables obtained from different virtual CCAs, tube laws and flow conditions. Although each artery presents a distinct vascular response, the differences obtained across different flow regimes are not significant. As expected, the linear tube law is unable to accurately capture all the haemodynamic features characterizing the current model. The findings from the sensitivity analysis are further confirmed by investigating the axial stretching effect on the CCA fluid dynamics. This factor does not seem to alter the pressure and flow waveforms. On the contrary, it is shown that, for an axially stretched vessel, the vascular wall exhibits an attenuation in absolute distension and an increase in circumferential stress, corroborating the findings of previous studies. This analysis shows that the new model offers a good balance between computational complexity and physics captured, making it an ideal framework for studies aiming to investigate the profound link between vascular mechanobiology and blood flow.

     

A model of a radially expanding and contracting lymphangion

The lymphatic system is an extensive vascular network featuring valves and contractile walls that pump interstitial fluid and plasma proteins back to the main circulation. Immune function also relies on the lymphatic system's ability to transport white blood cells. Failure to drain and pump this excess fluid results in edema characterized by fluid retention and swelling of limbs. It is, therefore, important to understand the mechanisms of fluid transport and pumping of lymphatic vessels. Unfortunately, there are very few studies in this area, most of which assume Poiseuille flow conditions. In vivo observations reveal that these vessels contract strongly, with diameter changes of the order of magnitude of the diameter itself over a cycle that lasts typically 2-3s. The radial velocity of the contracting vessel is on the order of the axial fluid velocity, suggesting that modeling flow in these vessels with a Poiseuille model is inappropriate. In this paper, we describe a model of a radially expanding and contracting lymphatic vessel and investigate the validity of assuming Poiseuille flow to estimate wall shear stress, which is presumably important for lymphatic endothelial cell mechanotransduction. Three different wall motions, periodic sinusoidal, skewed sinusoidal and physiologic wall motions, were investigated with steady and unsteady parabolic inlet velocities. Despite high radial velocities resulting from the wall motion, wall shear stress values were within 4% of quasi-static Poiseuille values. Therefore, Poiseuille flow is valid for the estimation of wall shear stress for the majority of the lymphangion contractile cycle.     

Pulmonary vascular compliance and viscoelasticity

When dog lung lobes were perfused at constant arterial inflow rate, occlusion of the venous outflow (VO) produced a rapid jump in venous pressure (Pv) followed by a slower rise in both arterial pressure (Pa) and Pv. During the slow rise Pa(t) and Pv(t) tended to converge and become concave upward as the volume of blood in the lungs increased. We compared the dynamic vascular volume vs. pressure curves obtained after VO with the static volume vs. pressure curves obtained by dye dilution. The slope of the static curve (the static compliance, Cst) was always larger than the slope of the dynamic curve (the dynamic compliance, Cdyn). In addition, the Cdyn decreased with increasing blood flow rate. When venous occlusion (VO) was followed after a short time interval by arterial occlusion (AO) such that the lobe was isovolumic, both Pa and Pv fell with time to a level that was below either pressure at the instant of AO. In an attempt to explain these observations a compartmental model was constructed in which the hemodynamic resistance and vascular compliance were volume dependent and the vessel walls were viscoelastic. These features of the model could account for the convergence and upward concavity of the Pa and Pv curves after VO and the pressure relaxation in the isovolumic state after AO, respectively. According to the model analysis, the difference between Cst and Cdyn and the flow dependence of Cdyn are due to wall viscosity and volume dependence of compliance, respectively. Model analysis also suggested ways of evaluating changes in the viscoelasticity of the lobar vascular bed. Hypoxic vasoconstriction that increased total vascular resistance also decreased Cst and Cdyn and appeared to increase the vessel wall viscosity.     

Mathematical modelling of dynamic adaptive tumour-induced angiogenesis: clinical implications and therapeutic targeting strategies

Angiogenesis, the growth of a network of blood vessels, is a crucial component of solid tumour growth, linking the relatively harmless avascular growth phase and the potentially fatal vascular growth phase. As a process, angiogenesis is a well-orchestrated sequence of events involving endothelial cell migration, proliferation; degradation of tissue; new capillary vessel (sprout) formation; loop formation (anastomosis) and, crucially, blood flow through the network. Once there is blood flow associated with the nascent network, the subsequent growth of the network evolves both temporally and spatially in response to the combined effects of angiogenic factors, migratory cues via the extracellular matrix and perfusion-related haemodynamic forces in a manner that may be described as both adaptive and dynamic. In this paper we present a mathematical model which simultaneously couples vessel growth with blood flow through the vessels--dynamic adaptive tumour-induced angiogenesis (DATIA). This new mathematical model presents a theoretical and computational investigation of the process and highlights a number of important new targets for therapeutic intervention. In contrast to earlier flow models, where the effects of perfusion (blood flow) were essentially evaluated a posteriori, i.e. after generating a hollow network, blood flow in the model described in this paper has a direct impact during capillary growth, with radial adaptations and network remodelling occurring as immediate consequences of primary anastomoses. Capillary network architectures resulting from the dynamically adaptive model are found to differ radically from those obtained using earlier models. The DATIA model is used to examine the effects of changing various physical and biological model parameters on the developing vascular architecture and the delivery of chemotherapeutic drugs to the tumour. Subsequent simulations of chemotherapeutic treatments under different parameter regimes lead to the identification of a number of new therapeutic targets for tumour management.     

Haemodynamic analysis of vessel remodelling in STA-MCA bypass for Moyamoya disease and its impact on bypass patency

The purpose of this study is to estimate the remodelling characteristics of STA-MCA bypass and its influence on patency via the use of computational fluid dynamic (CFD) technology. The reconstructed three-dimensional geometries from MRA were segmented to create computational domains for CFD simulations. Eleven patients, who underwent regular MRA both immediately following surgery and at the six months follow-up, were studied. The flow velocities at STA were measured via the use of quantitative MRA (QMRA) to validate simulation results. STA-MCA bypass patency was confirmed for each patient immediately following surgery. The simulation indicated that the remodelling of the arterial pedicle in nine patients was associated with a reduction in the resistance to flow through the bypass. For these cases, the modelling of a driving pressure of 10 mmHg through the bypass at 6 months post-surgery resulted in a 50% greater blood flow than those found immediately following surgery. However, two patients were found to exhibit contradictory patterns of remodelling, in which a highly curved bending at the bypass immediately post-surgery underwent progression, with increased resistance to flow through the bypass at 6 months follow-up, thereby resulting in a modelled flow rate reduction of 50% and 25%, respectively. This study revealed that STA-MCA bypass has a characteristic remodelling that usually reduces flow resistance. The initial morphology of the bypass may have had a significant effect on the outcome of vessel remodelling.     

Vascular metabolic dissipation in Murray's law

The metabolic dissipation in Murray's minimum energy hypothesis includes only the blood metabolism. The metabolic dissipation of the vascular tree, however, should also include the metabolism of passive and active components of the vessel wall. In this study, we extend the metabolic dissipation to include blood metabolism, as well as passive and active components of the vessel wall. The analysis is extended to the entire vascular arterial tree rather than a single vessel as in Murray's formulation. The calculations are based on experimentally measured morphological data of coronary artery network and the longitudinal distribution of blood pressure along the tree. Whereas the model includes multiple dissipation sources, the total metabolic consumption of a complex vascular tree is found to remain approximately proportional to the cumulative arterial volume of the unit. This implies that the previously described scaling relations for the various morphological features (volume, length, diameter, and flow) remain unchanged under the generalized condition of metabolic requirements of blood and blood vessel wall.     

Regulation of tyrosine phosphatases in the adventitia during vascular remodelling

Protein tyrosine phosphatases (PTPs) are regulators of growth factor signalling in vascular remodelling. The aim of this study was to evaluate PTP expression in the context of PDGF-signalling in the adventitia after angioplasty. Utilising a rat carotid artery model, the adventitial layers of injured and non-injured vessels were laser microdissected. The mRNA expression of the PDGF β-receptor, the ligands PDGF-A/B/C/D and the receptor-antagonising PTPs (DEP-1, TC-PTP, SHP-2, PTP1B) were determined and correlated to vascular morphometrics, proliferation markers and PDGF β-receptor phosphorylation. The levels of the PDGF β-receptor, PDGF-C and PDGF-D were upregulated concurrently with the antagonising PTPs DEP-1 and TC-PTP at day 8, and normalised at day 14 after vessel injury. Although the proliferation parameters were time-dependently altered in the adventitial layer, the phosphorylation of the PDGF β-receptor remained unchanged. The expression dynamics of specific PTPs indicate a regulatory role of PDGF-signalling also in the adventitia during vascular remodelling.     

An optimization principle for vascular radius including the effects of smooth muscle tone.

An optimization principle is proposed for the regulation of vascular morphology. This principle, which extends Murray's law, is based on the hypothesis that blood vessel diameter is controlled by a mechanism that minimizes the total energy required to drive the blood flow, to maintain the blood supply, and to support smooth muscle tone. A theoretical analysis reveals that the proposed principle predicts that the optimum shear stress on the vessel wall due to blood flow increases with blood pressure. This result agrees qualitatively with published findings that the fluid shear stress in veins is significantly smaller than it is in arteries.