 Research
 Open Access
 Published:
Hydraulic conductivity and lowdensity lipoprotein transport of the venous graft wall in an arterial bypass
BioMedical Engineering OnLinevolume 18, Article number: 50 (2019)
Abstract
Background
Blood flow condition may have influence upon the hydraulic conductivity of venous graft (L_{p,vein}) in an arterial bypass, then affecting the accumulation of lowdensity lipoproteins (LDLs) within the graft wall. To probe this possibility, we first measured in vitro the filtration rates of swine lateral saphenous vein segments under different flow rates, and the correlation of L_{p,vein} with wall shear stress (WSS) was then obtained.
Results
The experimental results showed that when WSS was very low, L_{p,vein} would increase drastically with WSS from 1.16 ± 0.15 × 10^{−11} m/s Pa at 0 dyn/cm^{2} to 2.17 ± 0.20 × 10^{−11} m/s Pa at 0.7 dyn/cm^{2}, then became constant of approximately 2.33 × 10^{−11} m/s Pa as the WSS increased further. Based on the experimental results, we assumed three different cases of L_{p,vein} and numerically simulated the LDLs transport in an arterial bypass model with venous graft. Case A: L_{p,vein} = 2.33 × 10^{−11} m/s Pa; Case B: L_{p,vein} = 1.16 × 10^{−11} m/s Pa (static condition with WSS of 0); Case C: L_{p,vein} was shear dependent. The simulation showed that the deposition/accumulation of LDLs within the venous graft wall in Case A was greatly enhanced when compared with that in Case B. However, the LDL accumulation in the graft wall was similar for Case A and Case C.
Conclusions
Our study, therefore, indicates that when the venous graft was implanted as a bypass graft, the L_{p,vein} might remain nearly constant along its whole length except for very few areas where the value of WSS was extremely low (less than 0.7 dyn/cm^{2}) and the effects of L_{p,vein} modulated by blood flow on LDL transport may be neglected.
Introduction
Autogenous vein (e.g., the greater saphenous vein) segments are widely used in vascular bypass surgery to relieve arterial occlusion. Nevertheless, when implanted into the arterial system as grafts, veins will develop a rapidly progressive and structurally diffusive form of atherosclerotic lesions which has been termed as “accelerated atherosclerosis” [1,2,3,4], and this has become the major cause of venous graft late failure [5, 6].
The mechanism of the accelerated atherogenesis in venous grafts has been studied extensively [3, 7]. It is well documented that the process of atherogenesis in venous grafts is very similar to the one in the arterial system, in which the deposition/accumulation of lipids such as lowdensity lipoproteins (LDLs) within the arterial wall is the initial event, and the high LDL concentration within the vessel wall could accelerate the atherogenesis [8].
From the viewpoint of mass transport, the hydraulic conductivity of blood vessel walls may play an important role in the accumulation of LDLs and atherosclerotic lesion development in the venous graft. In our previous study [9], we measured hydraulic conductivity of the swine lateral saphenous vein in vitro. Using the measured data, we numerically analyzed the transport of LDLs in the venous graft. Our result demonstrated that hydraulic conductivity of the venous graft was significantly higher than that of the host artery. High hydraulic conductivity could not only elevate concentration polarization of LDLs at the luminal surface of the venous graft, and then enhancing the transport of LDLs into the venous graft wall by the mechanism of concentration gradient, but also could directly lead to high influx of LDLs into the vessel wall by convective flow. As a result, LDLs would accumulate rapidly within the wall of the venous graft, in turn accelerating the genesis and development of atherosclerosis in the venous grafts.
The hydraulic conductivity of arteries and cultured endothelial cells has shown to be greatly affected by fluid flow [10]. For instance, using bovine aortic endothelial cell (BAEC) monolayers, Sill et al. [11] observed that fluid flow could cause an increase in hydraulic conductivity when compared to that in the absence of flow. Lever et al. found approximately 30% increase in hydraulic conductivity of the rabbit common carotid artery when the flow rate was increased from zero to a level of 10 mL/min [12]. However, there is a paucity of investigations on the effect of flow conditions on the hydraulic conductivity of the vein.
Since flow conditions may also affect the hydraulic conductivities of the venous graft (L_{p,vein}), in turn affecting LDL deposition/accumulation within the graft wall, in the present article, we measured in vitro the filtration rates of the swine lateral saphenous vein under different wall shear stress (WSS), from which the correlation of L_{p,vein} with wall shear stress was derived. We then numerically simulated LDL transport in a twodimensional bypass model with a venous graft and analyzed the deposition and accumulation of LDLs within the wall of venous graft.
Methods
Measurement of filtration rates under different WSS for venous wall
Preparation of vessel segments
The lateral saphenous veins obtained from farm swine were used for the study. Experiments in the present study were approved by the university ethics review board. The specifically designed metal supporting frame was used to fix the 6 to 7cmlong segment of blood excised from the swine leg. This procedure of the venous segment preparation for measurements has been described in detail previously [9].
Perfusion solution
For perfusion, freshly prepared albumin Krebs solution was used. Bovine serum albumin (Sigma Chemical Co., St. Louis, MO, USA) was dissolved in the Krebs solution (concentrations in mmol/L: NaCl, 118; KCl, 4.7; NaHCO_{3}, 25; KH_{2}PO_{4}, 1.2; MgSO_{4}, 1.2; CaCl_{2}, 2.5; glucose, 11) at a concentration of 1.0 mg/mL (density, 1.005 g/cm^{3}; viscous, 0.0116 g/cm s). The pH value of the solution was adjusted to 7.4.
Experimental procedure
The harvested vein segment fixed by a metal frame was encased in the stainless steel chamber (Fig. 1a) and connected to the experimental perfusion system (Fig. 1b), in which the upper reservoir could provide a steady flow through the vessel keeping the vessel pressurized at a preset constant pressure. Then the stainless steel chamber was filled with the Krebs solution and covered by a transparent Plexiglas plate. A calibrated capillary was installed on the Plexiglas cover plate with its orifice connecting to the inside of the chamber for measuring the increment of liquid volume within the chamber.
The filtration rate measurements of the vein were carried out under arterial pressure (hydrostatic pressure 100 mmHg, transmural pressure 70 mmHg) [13] with different flow rates, i.e., 0 mL/min, 20 mL/min, 40 mL/min, 60 mL/min, 80 mL/min, 100 mL/min, 200 mL/min and 500 mL/min. By adjusting the clamp on the outlet tubing connected to the distal end of the test blood vessel, the desired flow rate through the vessel could be obtained.
The filtration rate (v_{w}) across the test vessel wall was calculated by the following formula:
where ∆V is the liquid volume between two marked lines on the calibrated capillary, t is the time required for the liquid meniscus rising in the capillary to pass the two marked lines and S is the total area of the outer surface of the test vessel.
Based on the v_{w}, we could obtain the L_{p,vein} using Darcy’s law:
The wall shear stress (τ_{w}) acting on the luminal surface of the test vessel was given by
where μ and R represent, respectively, the viscosity of perfusion solution and the internal radius of the test vessel. Q is the flow rate through the test vessel.
Each experiment with a given flow rate was repeated three times and all experiments were carried out at a room temperature of 23 ± 0.5 °C. During the perfusing process, the length and external diameter of the vessel were measured for calculating the total area S of the vessel. With finishing each experiment turn, the test vessel was fixed at the test hemodynamics condition by perfusing with 10% formalin solution overnight to obtain the histological cross sections of the vessel, from which the internal radius R of the test vessel could be established.
Statistical analysis
SPSS19.0 (SPSS Inc.) software was applied for statistical analysis. All experiment results were expressed as mean ± standard deviation (SD). Oneway ANOVA was used to compare mean values of the experiment data. Differences were considered significant when P < 0.05.
Numerical simulation of LDL transport
Geometric model
Simplified bypass models (two dimensional) with 80% area occluded in the host artery were adopted in this simulation, as shown in Fig. 2. In the calculated model, both the host artery wall and the venous graft wall were treated as a singlelayer porous medium with an endothelium layer. \(\varGamma_{\text{adv}}\) was the interface between the media and adventitia. The outer diameter of the host artery (D_{a}) and the thickness of arterial wall (h_{a}) were set to be 5.2 and 0.35 mm, respectively. The venous graft’s outer diameter and thickness (D_{v} and h_{v}) were, respectively, 6.2 and 0.035 mm. The anastomosis angle was 45°. The host artery was 80% area occluded at the middle between points B and C.
Governing equations
To model bulk blood flow in the lumen and transmural flow in the endothelium layer, Navier–Stokes equations and Darcy’s law were employed, respectively. For modeling of mass balance, the convection–diffusion equation was employed in the lumen. To form a convection–diffusion–reaction equation in the endothelium layer, additional reaction term was added. Furthermore, the Kedem–Katchalsky equations were employed to couple the fluid dynamics and mass balance at membranes, i.e., endothelium.
Fluid flow dynamics
Blood flow in the vessel lumen was simulated by the steadystate incompressible Navier–Stokes and continuity equations [14]:
where u and p, respectively, represent fluid velocity vector and pressure. ρ and μ are the density (ρ =1050 kg/m^{3}) and the viscosity of blood (μ = 3.5 × 10^{−3} kg/m s), respectively [15]
The transmural fluid flow across the endothelium (J_{v}) was modeled by the Kedem–Katchalsky equation [16]:
where \(\Delta p\) and \(\Delta \pi\) are the hydraulic pressure and osmotic pressure differences across the endothelium; σ and L_{pend}, respectively, represent the osmotic reflection coefficient and the hydraulic conductivity of the endothelium. In the present simulation, ∆π was neglected because the influence of osmotic pressure on the fluid dynamics is too small, compared to hydraulic pressure.
The transmural fluid flow across the arterial wall and venous wall was modeled by Darcy’s law.
In a homogeneously permeable medium, Darcy’s law can be described as a simple proportional relationship between the instantaneous flow rate through a porous medium of permeability, the dynamic viscosity of the fluid and the pressure drop over a given distance. From the viewpoint of mass transport, Darcy’s law can be used to describe the hydraulic conductivity of blood vessel walls when the blood flow remains stable, laminar:
and continuity equation
where u_{w} and p_{w} represent, respectively, the velocity vector and pressure in the vessel wall, K_{w} is the hydraulic permeability of the wall and μ_{w} is the viscosity of plasma.
Mass transport
The transport of LDLs in the vessel lumen was modeled by the convection–diffusion equation as follows [17]:
where c represents the concentration of LDLs in the lumen, and D is the diffusion coefficient of LDLs in blood, D = 5.898 × 10^{−12} m^{2}/s [18].
The flux of LDLs across the endothelium (J_{s}) was described by the following equation [17]:
where ∆c is LDL concentration difference across the endothelium, P_{end} is the endothelial permeability to LDLs and \(\bar{c}\) is the mean concentration of LDLs in the endothelium. σ_{f} is the solvent reflection coefficient.
The transport of LDLs within the vessel wall was modeled by the convection–diffusion–reaction equation as follows [19]:
where c_{w} and D_{w}, respectively, represent the concentration of LDLs within vessel wall and the effective diffusivity of LDLs in the wall, r_{w} is the chemical reaction rate, and K_{lag} is the LDLs lag coefficient.
Boundary conditions
Boundary conditions for blood flow simulation
A fully developed (parabolic) velocity profile with a mean velocity of 0.22 m/s was applied at the inlet of the arterial lumen so that Reynolds number (Re) based on the internal diameter of the host artery was 300. At the outlet of arterial lumen, the pressure p was set at 100 mmHg. At the lumen side of the endothelial boundary, a wall to lumen transmural fluid velocity was prescribed by
At the wall side of the endothelial boundary, a lumentowall transmural velocity was prescribed by
where n_{l} and n_{w} represent the out normal vectors of the lumen subdomain and wall subdomain, respectively. J_{v} is the transmural fluid velocity across the endothelium.
A constant pressure boundary condition of 30 mmHg was prescribed at the media–adventitia interface (Γ_{adv}) [19].
At the occlusion of the host arterial lumen:
No flow was assumed for all other boundaries.
Boundary conditions for LDL transport
The inflow concentration of LDLs (c_{0}) at the lumen inlet boundary was taken as 28.6 × 10^{−3}nmol/mm^{3} [20]. At the occlusion of the arterial lumen, the flux of LDLs in the normal direction of the lumen domain was set at 0. At the lumen side of the endothelial boundary, a wall to lumen flux of LDLs was prescribed by
At the wall side of the endothelial boundary, a lumentowall flux of LDLs was prescribed by
For other boundaries, the concentration gradient of LDLs in the boundary normal direction was assumed to be zero.
Parameters
To solve the governing equations, we need to acquire the values of seven parameters, namely the hydraulic permeability of subendothelial wall (K_{w}), the hydraulic conductivity of endothelium (L_{pend}), the solvent reflection coefficient (σ_{f}), the endothelial permeability (P_{end}) to LDLs, the chemical reaction rate (r_{w}), the effective diffusivity of LDLs in the subendothelial layer (D_{w}) and the solute lag coefficient (K_{lag}).
Host artery
All the parameters were adopted in the literature, which are presented in Table 1.
Venous graft
The hydraulic conductivity of the venous wall (L_{p,vein}), hydraulic conductivity of the venous endothelium (L_{pend,vein}) and hydraulic conductivity of the subendothelial wall of the vein (L_{pw,vein}) were assumed to have the relation below [21]:
where L_{p,vein}= V_{w}/∆p, which could be obtained from the results of “Measurement of filtration rate for venous wall”. We assumed that the L_{pw,vein} was 4 folder of L_{pend, vein} when we performed the test experiments on vessels under no flow condition [19], then L_{pw,vein} and L_{pend,vein} under 0 dyn/cm^{2} of WSS could be acquired. Assuming that the L_{pw,vein} is constant with the flow rate, we could obtain L_{pend,vein} under other different WSS conditions.
Since \(L_{{\text{pw,vein}}} = \frac{{K_{{\text{w,vein}}} }}{{\mu_{{\text{w,vein}}} h_{{\text{vein}}} }}\) [22], and h_{vein} (the thickness of venous wall) could be obtained from the histological cross sections, we can get the hydraulic permeability of subendothelial wall of vein (K_{w,vein}). We assumed that other parameters of the vein were the same as those of the artery (Table 1).
Computation procedures
The numerical simulations were carried out using a validated finite elemental algorithm Comsol Multiphysics (COMSOL AB, Sweden). First, the pressure and the velocity fields were obtained by performing flow simulations, the solutions were used later on for the simulations of LDL transport. The threshold value for velocity and continuity residual detection was set to 1.0 × 10^{−5} to ensure the convergence of the calculated results. The numerical results for the flow and the transport of LDLs were determined to be mesh independent. In addition, we validated the mesh independence of results obtained. The criteria of mesh independence were set as the difference of J_{s} between the meshes used for computations and denser meshes was less than 3%. Based on these criteria, the final computational meshes consisted of 712,620 cubic elements.
Results
Hydraulic conductivity
In the present study, we first measured the filtration rates across (V_{w}) the walls of the venous segments under different WSS conditions. Then based on the measured data, we derived the correlation of L_{p,vein} with WSS using Eq. (2). The results are shown in Fig. 3a. It can be seen that when WSS is very low L_{p,vein} would increase dramatically from 1.16 ± 0.15 × 10^{−11} m/s Pa at WSS of 0 to 2.17 ± 0.20 × 10^{−11} m/s Pa at WSS of 0.7 dyn/cm^{2}. But beyond that L_{p,vein} would remain almost constant (approximately 2.33 × 10^{−11} m/s Pa).
Based on the values of L_{p,vein}, we calculated L_{pend,vein}, the hydraulic conductivity of the venous endothelium using Eq. (17). The obtained results are shown in Fig. 3b. By curve fitting, it can be shown that L_{pend,vein}= 2.522 × 10^{−11} × WSS^{0.1} + 1.532 × 10^{−11}.
Numerical comparison of LDL transport
Since the experimental measurements showed that when WSS was greater than 0.7 dyn/cm^{2}, L_{p,vein} would remain almost constant at a value of approximately 2.33 × 10^{−11} m/s Pa and the numerical simulation of blood flow indicated that there were very few areas in the bypass model where WSS would be less than 0.7 dyn/cm^{2}, in the present numerical study, we assumed three cases of L_{p,vein} to compute the transport and accumulation of LDLs in the bypass mode. Case A: L_{p,vein} = a constant value of 2.33 × 10^{−11} m/s Pa; Case B: L_{p,vein} = 1.16 × 10^{−11} m/s Pa (static condition with WSS of 0); Case C: L_{p,vein} is shear dependent as derived from the experimental measurements.
Figure 4a shows LDL concentration (c_{w}) distributions in the middle position of the venous graft wall for the three cases. As evident from the figure, the distribution of c_{w} for Case A and Case C is very similar. However, when compared with Case B, c_{w} for Cases A and C is generally higher in most regions of the venous graft wall. Figure 4b gives the comparison of c_{w} among the three cases in terms of percentage difference. It can be seen that the difference between Case A and Case B is prominent. For instance, the percentage difference is more than 10% in most regions, especially in the exit region of the graft inner wall (line 22), where the percentage difference reaches as high as 26%. But for Case A and Case C, there is no evident difference in all regions of the graft. The highest difference in c_{w} between the two cases is located in the middle section of the outer wall (line 11), which is less than 5%.
Effect of arterial stenosis severity on LDL transport
In clinical practice, the stenosis degree is an important indicator of necessity for bypass surgery. Usually, when the degree of arterial narrowing is greater than 70%, surgical intervention would be eminent. Besides, after the surgery, the stenosis of the host artery would continue to develop, resulting in the alteration of the flow field in the venous graft, which would, in turn, affect the transport of LDLs in the graft. To clarify the impact of stenosis severity on LDL deposition/accumulation within the venous graft wall, we numerically simulated LDL transport in bypass models with different arterial stenosis that has 70, 80 and 90% reduction in the crosssectional area of the host artery, using the sheardependent L_{p,vein}.
As shown in Fig. 5a, the distributions of c_{w} along the venous graft have much in common for the three bypass models with different stenoses. Nevertheless, it is still evident that the positions of the local minima of c_{w}, locating at the middle of the graft outer wall, are different. The percentage differences in the c_{w} for the three models are presented in Fig. 5b. From the figure, it can be seen that there is no significant difference in c_{w} between the models of 70 and 80% occlusion. Only when comparing the 90% model with the 70% one, the percentage difference is somewhat prominent, the maximum of which is approximately 15%.
Discussion
The late failure of venous grafts, which is usually caused by accelerated atherogenesis, has long been a common problem that limits the therapy of arterial bypass surgery. We previously demonstrated that enhanced L_{p,vein} could lead to fast infiltration/deposition of LDLs within the venous graft wall, and hence the accelerated genesis and development of atherosclerosis [9].
As hydraulic conductivity of blood vessels can be affected by blood flow, in the present study we investigated the effect of L_{p,vein} modulated by flow on LDL accumulations within the wall of the venous graft. We first measured the filtration rates of the swine lateral saphenous vein (hence the L_{p,vein}) under different flow/WSS conditions. To our best knowledge, this is the first in vitro study of its kind to clarify the shear dependence of venous wall hydraulic conductivity, though a few studies have investigated the hydraulic conductivities of endothelial cell monolayers, capillaries and arterial vessels [10, 12, 26,27,28,29]. From our experimental results, it was found that the luminal flow in the venous vessels investigated could cause a significant increase in L_{p,vein} when compared with the case of flow absence. This finding was consistent with the studies performed on arteries, capillaries and endothelial cell monolayers by others, which indicated that flow shear could induce increases in hydraulic conductivity via a NOdependent mechanism [10, 30]. But, our results showed that the increase in hydraulic conductivity of venous segments induced by a change of WSS (from zero to a steady level of more than 0.7 dyn/cm^{2}) was approximately 100%, which was much higher than the results obtained on intact arterial vessels (e.g., the excised rabbit carotid artery [12] that exhibited only 30% increase in hydraulic conductivity when WSS elevated from 0 to 1 dyn/cm^{2}), but was lower than those obtained on endothelial cell monolayers (e.g., exposure of cultured bovine aortic endothelial monolayers to a WSS of 12 dyn/cm^{2} elicited a 2.37fold increase in hydraulic conductivity [10]). This discrepancy might be attributed to the structural differences in the subendothelial layers of the test models. It had been documented that the subendothelial layers of blood vessels, including the intima and media, contribute most of the overall resistance to filtration flow across the vessel wall [31, 32]. Considering that the endothelial cell monolayers have no subendothelial layer and the subendothelial layer in the artery is thicker and compacter than that of the vein, the sensitivity of the venous wall hydraulic conductivity in response to WSS should be between those of the arterial wall and the endothelial cell monolayers.
Based on the obtained hydraulic conductivity, we simulated numerically the transport of LDLs within the venous graft of an arterial bypass model. It had been documented that LDL transport, including the transport in the lumen (from the bulk blood flow to the endothelium) and the transendothelial transport (across the endothelium), was mainly affected by hydraulic conductivity [18]. Indeed, our numerical results revealed a significant difference in c_{w} between Case A (with flow and L_{p,vein} = 2.33 × 10^{−11} m/s Pa) and Case B (with no flow of L_{p,vein}= 1.16 × 10^{−11} m/s Pa). This finding suggests that flow could affect L_{p,vein}, then affecting LDL deposition/accumulation within the graft wall.
In addition, our result showed that there was no evident difference in c_{w} between Case A and Case C (in which L_{p,vein} was assumed to be shear dependent), indicating that the influence of sheardependent L_{p,vein} on LDL transport was marginal. This seems different from some previous studies that demonstrated a significant influence of hydraulic conductivity with sheardependent nature on LDL transport. For instance, using an idealized model of a stenosed artery, Sun et al. [19] found that when compared with the independent case, the sheardependent L_{p} could lead to apparent differences in LDL concentration within the subendothelium, in which the difference could reach as high as 12%. We suppose that the discrepancy between our results and others could be explained as follows. Different from the artery, hydraulic conductivity of the vein was almost independent of the flowinduced wall shear stress except for when the WSS was extremely low. Since there were very few regions in the bypass model with very low WSS, the hydraulic conductivities of the venous graft wall for Case C (the sheardependent case) had to be similar to that in Case A (the independent one).
In the present study, we also investigated the effect of stenosis severity of the host artery on LDL infiltration/accumulation within the venous graft wall. The results showed that there were some moderate differences in c_{w} among the cases with different stenosis in the host artery. When compared the 90% occlusion model with the 70% one, the maximum percentage difference in the c_{w} reaches 15%. We believe that the effect of the stenosis severity on c_{w} is mostly caused by the luminal blood flow per se rather than the filtration rate across the graft wall because the hydraulic conductivity of the venous graft is almost constant along the venous graft.
Here, it should be mentioned that to simplify the simulation, the complex heterogeneous structure of the vessel wall was approximated by a simple homogeneous layer. In recent years, multilayer models were used to characterize LDL transport in blood vessel walls [18, 33, 34], which could better describe the distribution of LDLs in the venous graft. One thing that should be addressed is that the hemodynamic condition can be influenced by several factors (individual geometric parameters, material parameters, realistic blood flow condition). Investigators have paid much more attention to these aspects. For example, Peng et al. investigated the correlation between hemodynamic parameters and individual geometric factors in the patientspecific coronary arteries, and their results showed that the increasing severity of the stenosis was associated with the increased maximum wall shear stress at the stenosis region [35]. Xiong et al. calculated the blood flow field distribution in idealized 2D models, and results showed that high pressure drop would create a pathological WSS environment for the further growth of the plaque [36]. Liu et al. investigated the hemodynamic alteration in the cerebral circulation and revealed that the structural status could provide comprehensive information about the hemodynamic alterations in the pathological circle of Wills arterial structures [37]. Zhang et al. introduced a pressurebased carotid arterial functional assessment index generated from computational fluid dynamic simulation based on digital subtracted angiography (DSA) data. Their results showed that the realistic blood condition obtained through invasive approach would remain one efficient way to study the relationship between hemodynamic disorder caused by internal carotid artery (ICA) stenosis and subsequent perfusion variations in brain [14]. Based on their investigations, it can be addressed that several significant parameters in terms of realistic geometry, artery material and flow condition would certainly influence the hemodynamic behavior, and then influence the LDL transport within the blood artery from the viewpoint of mass transport. Though this can be significant limitations for our study, our study remains a preliminary study as we focus on the correlation between the wall shear stress and hydraulic conductivity of venous endothelium, and its influence on the transport phenomenon of LDL. Further studies based on more comprehensive models and realistic blood conditions still need to be carried out to strengthen our study.
Conclusion
To probe the possibility that the blood flow condition would affect the accumulation of LDLs within the graft by influencing the hydraulic conductivity of venous graft in arterial bypass, LDL transport phenomenon in a twodimensional bypass model with a venous graft was numerically simulated, while the deposition and accumulation of LDLs within the wall of venous graft were analyzed based on the filtration rates of the swine lateral saphenous vein under different WSS measured from the in vitro experiments. Our results indicate that the blood flow would affect the transport and accumulation of LDLs within the venous graft wall when the WSS is extremely low, as the L_{p,vein} can be significantly influenced by the WSS with the range of 0–0.7 dyn/cm^{2}. While when WSS is above 0.7 dyn/cm^{2}, the influence of sheardependent L_{p,vein} on LDL transport is marginal. Since the low WSS (< 0.7 dyn/cm^{2}) areas remain few, it can be addressed that the effect of L_{p,vein} modulated by blood flow on LDL transport remains weak.
Abbreviations
 LDLs:

lowdensity lipoproteins
 WSS:

wall shear stress
 BAEC:

bovine aortic endothelial cell
 SD:

standard deviation
 Re :

Reynolds number
 DSA:

digital subtracted angiography
 ICA:

internal carotid artery
References
 1.
Shukla N, Jeremy JY. Pathophysiology of saphenous vein graft failure: a brief overview of interventions. Curr Opin Pharmacol. 2012;12(2):114–20.
 2.
Raja SG, Haider Z, Ahmad M, Zaman H. Saphenous vein grafts: to use or not to use? Heart Lung Circ. 2004;13(4):403–9.
 3.
Ward AO, Caputo M, Angelini GD, George SJ, Zakkar M. Activation and inflammation of the venous endothelium in vein graft disease. Atherosclerosis. 2017;265:266–74.
 4.
Newby AC, Zaltsman AB. Molecular mechanisms in intimal hyperplasia. J Pathol. 2000;190(3):300–9.
 5.
Kim FY, Marhefka G, Ruggiero NJ, Adams S, Whellan DJ. Saphenous vein graft disease: review of pathophysiology, prevention, and treatment. Cardiol Rev. 2013;21(2):101–9.
 6.
de Vries MR, Simons KH, Jukema JW, Braun J, Quax PH. Vein graft failure: from pathophysiology to clinical outcomes. Nat Rev Cardiol. 2016;13(8):451–70.
 7.
Raggi P, Genest J, Giles JT, Rayner KJ, Dwivedi G, Beanlands RS, Gupta M. Role of inflammation in the pathogenesis of atherosclerosis and therapeutic interventions. Atherosclerosis. 2018;276:98–108.
 8.
Hu J, Xi D, Zhao J, Luo T, Liu J, Lu H, Li M, Xiong H, Guo Z. Highdensity lipoprotein and inflammation and its significance to atherosclerosis. Am J Med Sci. 2016;352(4):408–15.
 9.
Wang Z, Liu X, Kang H, Sun A, Fan Y, Deng X. Enhanced accumulation of LDLs within the venous graft wall induced by elevated filtration rate may account for its accelerated atherogenesis. Atherosclerosis. 2014;236(1):198–206.
 10.
Kang H, Cancel LM, Tarbell JM. Effect of shear stress on water and LDL transport through cultured endothelial cell monolayers. Atherosclerosis. 2014;233(2):682–90.
 11.
Sill HW, Chang YS, Artman JR, Frangos JA, Hollis TM, Tarbell JM. Shear stress increases hydraulic conductivity of cultured endothelial monolayers. Am J Physiol. 1995;268(2 Pt 2):H535–43.
 12.
Lever MJ, Tarbell JM, Caro CG. The effect of luminal flow in rabbit carotid artery on transmural fluid transport. Exp Physiol. 1992;77(4):553–63.
 13.
Armstrong J, Narracott AJ, Milton R, Galea J, Cooper GJ, Lawford PV, Hose DR, Cumberland DC, Holt CM. Development of an ex vivo model to investigate the effects of altered haemodynamics on human bypass grafts. J Med Eng Technol. 2000;24(5):183–91.
 14.
Zhang D, Xu P, Qiao H, Liu X, Luo L, Huang W, Zhang H, Shi C. Carotid DSA based CFD simulation in assessing the patient with asymptomatic carotid stenosis: a preliminary study. Biomed Eng Online. 2018;17(1):31.
 15.
Liu X, Wang L, Wang Z, Li Z, Kang H, Fan Y, Sun A, Deng X. Bioinspired helical graft with taper to enhance helical flow. J Biomech. 2016;49(15):3643–50.
 16.
Kedem O, Katchalsky A. Thermodynamic analysis of the permeability of biological membranes to nonelectrolytes. Biochem Biophys Acta. 1958;27(2):229–46.
 17.
Sun N, Wood NB, Hughes AD, Thom SA, Xu XY. Effects of transmural pressure and wall shear stress on LDL accumulation in the arterial wall: a numerical study using a multilayered model. Am J Physiol Heart Circ Physiol. 2007;292(6):H3148–57.
 18.
Liu X, Fan Y, Deng X. Effect of the endothelial glycocalyx layer on arterial LDL transport under normal and high pressure. J Theor Biol. 2011;283(1):71–81.
 19.
Sun N, Wood NB, Hughes AD, Thom SA, Xu XY. Fluidwall modelling of mass transfer in an axisymmetric stenosis: effects of sheardependent transport properties. Ann Biomed Eng. 2006;34(7):1119–28.
 20.
Katz MA. New formulation of water and macromolecular flux which corrects for nonideality: theory and derivation, predictions, and experimental results. J Theor Biol. 1985;112(2):369–401.
 21.
Tarbell JM, Lever MJ, Caro CG. The effect of varying albumin concentration of the hydraulic conductivity of the rabbit common carotid artery. Microvasc Res. 1988;35(2):204–20.
 22.
Wang S, Vafai K. Analysis of low density lipoprotein (LDL) transport within a curved artery. Ann Biomed Eng. 2015;43(7):1571–84.
 23.
Karner G, Perktold K, Zehentner HP. Computational modeling of macromolecule transport in the arterial wall. Comput Methods Biomech Biomed. 2001;4(6):491–504.
 24.
Meyer G, Merval R, Tedgui A. Effects of pressureinduced stretch and convection on lowdensity lipoprotein and albumin uptake in the rabbit aortic wall. Circ Res. 1996;79(3):532–40.
 25.
Sun N, Wood NB, Hughes AD, Thom SA, Xu XY. Influence of pulsatile flow on LDL transport in the arterial wall. Ann Biomed Eng. 2007;35(10):1782–90.
 26.
Kim MH, Harris NR, Tarbell JM. Regulation of capillary hydraulic conductivity in response to an acute change in shear. Am J Physiol Heart Circ Physiol. 2005;289(5):H2126–35.
 27.
Lakshminarayanan S, Gardner TW, Tarbell JM. Effect of shear stress on the hydraulic conductivity of cultured bovine retinal microvascular endothelial cell monolayers. Curr Eye Res. 2000;21(6):944–51.
 28.
LopezQuintero SV, Amaya R, Pahakis M, Tarbell JM. The endothelial glycocalyx mediates shearinduced changes in hydraulic conductivity. Am J Physiol Heart Circ Physiol. 2009;296(5):H1451–6.
 29.
LopezQuintero SV, Ji XY, Antonetti DA, Tarbell JM. A threepore model describes transport properties of bovine retinal endothelial cells in normal and elevated glucose. Investig Ophthalmol Vis Sci. 2011;52(2):1171–80.
 30.
Chang YS, Yaccino JA, Lakshminarayanan S, Frangos JA, Tarbell JM. Shearinduced increase in hydraulic conductivity in endothelial cells is mediated by a nitric oxidedependent mechanism. Arterioscler Thromb Vasc Biol. 2000;20(1):35–42.
 31.
Dabagh M, Jalali P, Tarbell JM. The transport of LDL across the deformable arterial wall: the effect of endothelial cell turnover and intimal deformation under hypertension. Am J Physiol Heart Circ Physiol. 2009;297(3):H983–96.
 32.
Tedgui A, Lever MJ. Filtration through damaged and undamaged rabbit thoracic aorta. Am J Physiol. 1984;247(5 Pt 2):H784–91.
 33.
Yang N, Vafai K. Modeling of lowdensity lipoprotein (LDL) transport in the arteryeffects of hypertension. Int J Heat Mass Transf. 2006;49(5–6):850–67.
 34.
Kenjeres S, de Loor A. Modelling and simulation of lowdensity lipoprotein transport through multilayered wall of an anatomically realistic carotid artery bifurcation. J R Soc Interface. 2014;11(91):20130941.
 35.
Peng C, Wang X, Xian Z, Liu X, Huang W, Xu P, Wang J. The impact of the geometric characteristics on the hemodynamics in the stenotic coronary artery. PLoS ONE. 2016;11(6):e0157490.
 36.
Xiong H, Liu X, Tian X, Pu L, Zhang H, Lu M, Huang W, Zhang YT. A numerical study of the effect of varied blood pressure on the stability of carotid atherosclerotic plaque. Biomed Eng online. 2014;13:152.
 37.
Liu X, Gao Z, Xiong H, Ghista D, Ren L, Zhang H, Wu W, Huang W, Hau WK. Threedimensional hemodynamics analysis of the circle of Willis in the patientspecific nonintegral arterial structures. Biomech Model Mechanobiol. 2016;15(6):1439–56.
Authors’ contributions
ZW and ML participated in the design of the study and drafted the manuscript; XL and AS revised the manuscript; YF and XD coordinated the study and helped draft the manuscript. All authors read and approved the final manuscript.
Acknowledgements
We would like to acknowledge the funding agencies for the support of the work.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Consent for publication
We consent for the publication of this work.
Ethics approval and consent to participate
Not applicable.
Funding
This work was supported by GrantsinAid from the National Key Research and Development Plan of China (2017YFB0702501), National Natural Science Foundation of China (Nos. 11332003, 31170904, 31570947,11602064,11421202,11472031), and Natural Science Foundation of Jiangsu Province (BK20161366).
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Venous graft
 Lowdensity lipoprotein
 Hydraulic conductivity
 Wall shear stress