- Research
- Open Access

# Low-density lipoprotein concentration in the normal left coronary artery tree

- Johannes V Soulis†
^{1}Email author, - George D Giannoglou†
^{2}, - Vassilios Papaioannou†
^{1}, - George E Parcharidis†
^{2}and - George E Louridas†
^{2}

**7**:26

https://doi.org/10.1186/1475-925X-7-26

© Soulis et al; licensee BioMed Central Ltd. 2008

**Received:**17 December 2007**Accepted:**17 October 2008**Published:**17 October 2008

## Abstract

### Background

The blood flow and transportation of molecules in the cardiovascular system plays a crucial role in the genesis and progression of atherosclerosis. This computational study elucidates the Low Density Lipoprotein (LDL) site concentration in the entire normal human 3D tree of the LCA.

### Methods

A 3D geometry model of the normal human LCA tree is constructed. Angiographic data used for geometry construction correspond to end-diastole. The resulted model includes the LMCA, LAD, LCxA and their main branches. The numerical simulation couples the flow equations with the transport equation applying realistic boundary conditions at the wall.

### Results

High concentration of LDL values appears at bifurcation opposite to the flow dividers in the proximal regions of the Left Coronary Artery (LCA) tree, where atherosclerosis frequently occurs. The area-averaged normalized luminal surface LDL concentrations over the entire LCA tree are, 1.0348, 1.054 and 1.23, for the low, median and high water infiltration velocities, respectively. For the high, median and low molecular diffusivities, the peak values of the normalized LDL luminal surface concentration at the LMCA bifurcation reach 1.065, 1.080 and 1.205, respectively. LCA tree walls are exposed to a cholesterolemic environment although the applied mass and flow conditions refer to normal human geometry and normal mass-flow conditions.

### Conclusion

The relationship between WSS and luminal surface concentration of LDL indicates that LDL is elevated at locations where WSS is low. Concave sides of the LCA tree exhibit higher concentration of LDL than the convex sides. Decreased molecular diffusivity increases the LDL concentration. Increased water infiltration velocity increases the LDL concentration. The regional area of high luminal surface concentration is increased with increasing water infiltration velocity. Regions of high LDL luminal surface concentration do not necessarily co-locate to the sites of lowest WSS. The degree of elevation in luminal surface LDL concentration is mostly affected from the water infiltration velocity at the vessel wall. The paths of the velocities in proximity to the endothelium might be the most important factor for the elevated LDL concentration.

## Keywords

- Wall Shear Stress
- Molecular Diffusivity
- Left Coronary Artery
- Luminal Surface
- High Wall Shear Stress

## Background

Elucidating the blood flow and the transport of macromolecules in the cardiovascular system is essential in understanding the genesis and progression of atherosclerosis [1, 2]. Wall Shear Stress (WSS) may affect the endothelial permeability [3, 4]. Regional variations in the permeability of arterial endothelium may contribute to the localization of atherosclerosis [5]. The transportation of Low-Density Lipoproteins (LDL) across the artery wall is considered to be a step of paramount importance in atherosclerosis [6, 7]. Atherosclerosis shows a predilection in regions of the arterial tree with hemodynamic particularities, such as local disturbances of WSS in space, and locally high concentrations of lipoprotein [8, 9] and [10]. In proximal Left Coronary Artery (LCA) tree regions, where atherosclerosis frequently occurs, low WSS appears [9]. The local velocity, the molecular viscosity disturbances and the morphological (geometrical) particularities may also predispose to the formation of coronary atheromatic plaques [11]. Although the WSS has been widely proved to affect the arterial segmental biology, the near-wall localization of critical macromolecular blood particles may significantly contribute to the development of atherosclerotic plaques [12]. LDL is one of these particles, placed there by transport and diffusion.

The strategy and some of the pros and cons of computational modeling approach of the coupled fluid and mass flow (focusing on results from studies made by others) on a variety of arterial geometries have been reviewed by Ethier [1]. He concluded that elevated LDL co-localize with known sites of atherosclerotic plaque development. Besides atherosclerotic regions are low WSS regions. Consequently, it was difficult to determine whether it was low WSS, abnormal mass transfer or both that were contributing to astherosclerosis. A theoretical study of a straight artery by Wada et al. suggested the possibility that all vascular phenomena were governed by the flow-dependent concentration polarization of LDL, which carries cholesterol [13]. A multiple bend was studied to elucidate the mechanisms of localization of atherosclerotic lesions [14]. They drove into the conclusion that regions of elevated LDL luminal surface concentration did not necessarily co-located to the sites of lowest WSS. They used constant value for the blood molecular viscosity, instead of a non-linear one [15].

The accumulation of LDL in vascular districts featuring highly disturbed flow was examined [16]. Geometrical parameters such as curvature and variations of the luminal section strongly influence the LDL within the wall. Realizing that the research in macromolecular transport of curved 3D blood vessels was rare, steady and unsteady flow and mass simulation was analyzed [17]. However, the blood was treated as a Newtonian fluid, thus limiting the quantitative results. They concluded that the concentration of LDL along the aortic axis yields higher values at the outer wall (concave side) compared to inner wall (convex side).

The current study is a three-dimensional, numerical simulation that couples the flow equations with the transport equation, applying realistic boundary conditions, in terms of blood-side concentration at the wall. The blood is considered to be non-Newtonian fluid obeying to the power law. Emphasis has been put on: a) LDL (concentration) distribution, b) factors affecting mass transport from flowing blood to arterial wall using various water infiltration velocity, and molecular diffusivity, c) WSS distribution. We demonstrate that due to the semi-permeable nature of the arterial walls, high concentrations of LDL occur at the endothelium. The LDL variation is most noticeable at bifurcations opposite to the flow dividers and at the concave parts of the bent arterial segments.

## Methods

### Geometry and computational grid

### Flow equations and assumptions

*ρ*(kg/m

^{3}) is the density,

*t*(sec) is the time, $\overrightarrow{u}$ (m/s) is the velocity vector. The conservation of momentum is written,

*p*(N/m

^{2}) is the static pressure, $\overline{\tau}$ (N/m

^{2}) is the shear stress tensor and

*ρ*$\overrightarrow{g}$ (N/m

^{3}) is the gravitational body force. The shear stress tensor $\overline{\tau}$ is given by,

*μ* is the molecular viscosity, *I* is the unit tensor, and the second term in the right hand side is the effect of the volume dilation.

*η*($\dot{S}$), is given by,

*k*is 0.01691 (kg-s^n-2/m), the power-law index

*n*is 0.7,

*T*(K) and

*T*

_{ o }(K) are local and reference temperatures, respectively. The components of the WSS possibly have different effects upon endothelial cells. Some components, mainly those being diagonal, generate intercellular tension while the off-diagonal components possibly contribute to intercellular shearing forces [23]. Thus, the actual shear stress is given by,

*C*(mg/ml) is the LDL concentration, $\overrightarrow{J}$ is the diffusion flux of LDL, which arises due to concentration gradients. The diffusion flux is written as,

*D* (m^{2}/s) is the diffusion coefficient of LDL in the mixture.

### Blood diffusivity, water infiltration velocity and endothelial permeability

In the absence of systematic and reliable experimental data for human arteries, we assume the molecular diffusivity D to be isotropic and constant [13, 24]. Various molecular diffusivities and water infiltration velocities V_{w} are applied to calculate the LDL luminal surface concentration over the LCA tree. There is a limited amount of experimental data about the endothelial permeability K (mass transfer coefficient). Thus, specifications of K values are sometimes unreliable. In the current analysis the K value is set equal to 2.0 × 10^{-10} m/s [14].

### Flow boundary conditions

The blood velocity is assumed to be uniform at the orifice of LMCA. The applied inflow conditions mimic typical coronary blood averaged flow velocities of 0.05 m/s under resting conditions, corresponding to Reynolds number of 95. The blood density is set equal to 1058 kg/m^{3}. Flow discharges are set analogous to the third power of the branching vessel inlet diameter according to Murray's law [25].

### Mass transport boundary conditions

*C*

_{ o }of LDL (= 1.3 mg/ml) is applied at the coronary tree inlet. At the coronary artery vessel outlets, the gradient of LDL concentration along the vessel is set equal to zero (zero flux, Newmann condition) $\frac{\partial C}{\partial s}=0$,

*s*is the unit vector normal to the outlet surfaces of the LCA. Suitable mass transport condition must be specified at the wall. The boundary conditions at the wall can be described as,

*C*
_{
w
}(mg/ml) is the concentration at the endothelial surface, *V*
_{
w
}(m/s) is the transmural component of the fluid velocity at the wall, *n* is the direction normal to the wall. The boundary conditions described in Eq. (9) state that the net amount of LDL (= *K C*
_{
w
}) passing from the endothelium to the vessel wall is determined from the difference of the mass flow carried to the vessel wall by infiltration flow (= *C*
_{
w
}
*V*
_{
w
}) and the amount of flow which diffuses back to the main vessel flow (= $D\frac{\partial C}{\partial n}$). It is known that the concentration of LDL at an arterial wall is higher than that in the bulk flow and it increases with increasing infiltration velocity [13].

### Solution

The user defined function capabilities of the numerical code are also incorporated to simulate the mass transport boundary condition across the endothelium, Eq. (9). For a typical satisfactory convergence solution, a total of 1000 pseudo-time steps are required. Convergence is achieved when all velocity components, fluid flow, and mass flow changes dropped below 10^{-8}.

## Results

### WSS and luminal surface LDL concentration at the LMCA bifurcation. LDL versus WSS over the entire LCA tree

All LDL luminal surface concentration values Cw are normalized with the inlet value of LDL concentration Co (= 1.3 mg/ml). The applied flow conditions refer to inlet based Reynolds number Re_{o} = 95 using V_{w} = 0.6 × 10^{-8} m/s, uniform constant LDL concentration, endothelial permeability K = 2.0 × 10^{-10} m/s and diffusivity D = 15.0 × 10^{-12} m^{2}/s.

^{2}to 2.50 N/m

^{2}, Fig. 3. High curvature affects the velocity distribution at the flow divider, giving rise to high WSS values. The normalized LDL luminal surface concentration values (= Cw/Co) are shown in Fig. 4. Elevated LDL values appear at regions opposite to the flow divider at either LMCA, or proximal LAD or LCxA branches. The spatial gradients of the LDL are considerable. The results indicate that at the main tree bifurcation the normalized LDL luminal surface concentration values range between 1.0369–1.0577.

_{o}= 640. The area averaged normalized LDL concentration over the entire LCA tree is 1.0347.

Area averaged LDL concentration Cw/Co over LCA tree segments

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0345 | 1.0347 | 1.0338 | 1.0346 | 1.0343 | 1.0350 | 1.0350 | 1.0353 | 1.0361 | 1.0349 |

### LDL at distal bend surface of the LCA

### LMCA bifurcation. Water infiltration effect in luminal surface LDL concentration

Various water infiltration velocities V_{w} are applied to calculate the LDL luminal surface concentration and WSS over the LCA tree. An inlet flow velocity of 0.05 m/s (Re_{o} = 95) is applied to simulate the fluid flow-mass transport problem using endothelial permeability K = 2.0 × 10^{-10} m/s, and molecular diffusivity D = 15.0 × 10^{-12} m^{2}/s. The applied infiltration velocities are 0.6 × 10^{-8} m/s (low), 1.0 × 10^{-8} m/s (median), and 4.0 × 10^{-8}m/s (high).

^{-11}g/s, 4.32 × 10

^{-11}g/s and 20.17 × 10

^{-11}g/s for the low, median and high infiltration velocities, respectively.

Water velocity filtration effects at V_{w} = 0.6 × 10^{-8} m/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0345 | 1.0346 | 1.0339 | 1.0347 | 1.0344 | 1.0347 | 1.0350 | 1.0352 | 1.0361 | 1.0350 |

Water velocity filtration effects at V_{w} = 1.0 × 10^{-8} m/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0548 | 1.0542 | 1.0528 | 1.0540 | 1.0535 | 1.0546 | 1.0547 | 1.0551 | 1.0567 | 1.0548 |

Water velocity filtration effects at V_{w} = 4.0 × 10^{-8} m/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.2352 | 1.2311 | 1.2219 | 1.2280 | 1.2277 | 1.2314 | 1.2340 | 1.2344 | 1.2460 | 1.2352 |

### LMCA bifurcation. Molecular diffusivity effect in luminal surface LDL concentration

^{-12}m

^{2}/s (low), 10.0 × 10

^{-12}m

^{2}/s (median) and 15.0 × 10

^{-12}m

^{2}/s (high) at Re

_{o}= 95, K = 2.0 × 10

^{-10}m/s and V

_{w}= 0.6 × 10

^{-8}m/s are applied. The normalized luminal surface LDL concentration for low, median and high molecular diffusivity values are shown in Figs. 14, 15 and 16, respectively. High molecular diffusivity results into low LDL concentration for all LCA tree segments. Furthermore, low molecular diffusivity results into increased area of elevated LDL throughout the LCA tree. High concentration values occur at regions opposite to the flow dividers. For the high, median and low molecular diffusivities, the peak normalized LDL luminal surface concentrations at the LMCA bifurcation reach 1.065, 1.080 and 1.205, respectively. The effect of molecular diffusivity upon the values and the patterns of WSS distribution are marginal (not shown).

^{-11}g/s, 2.57 × 10

^{-11}g/s and 2.54 × 10

^{-11}g/s corresponding to the low, median and high molecular diffusivities, respectively.

Molecular diffusivity effects at D = 5.0 × 10^{-12} m^{2}/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0918 | 1.0932 | 1.0909 | 1.0932 | 1.0925 | 1.0940 | 1.0947 | 1.0950 | 1.0985 | 1.0947 |

Molecular diffusivity effects at D = 10.0 × 10^{-12} m^{2}/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0479 | 1.0486 | 1.0475 | 1.0487 | 1.0482 | 1.0493 | 1.0493 | 1.0495 | 1.0510 | 1.0493 |

Molecular diffusivity effects at D = 15.0 × 10^{-12} m^{2}/s

LCxA | OM1 | OM2 | S1 | D1 | S2 | D2 | S3 | D3 | LAD |
---|---|---|---|---|---|---|---|---|---|

1.0342 | 1.0346 | 1.0339 | 1.0347 | 1.0344 | 1.0347 | 1.0350 | 1.0352 | 1.0361 | 1.0350 |

## Discussion

### Mass transfer importance

Main purpose of this computational study is to elucidate the LDL site concentration under various mass conditions. The arising are the degree of elevation of the luminal surface LDL concentration and the area of the LCA tree that it occupies in relation to water infiltration and size of macromolecules. The amount of LDL transported to the arterial wall is a function of the surface LDL elevation. However, the transported LDL is related to the permeability of the arterial wall, which in turn is a function of the WSS. The blood and the arterial wall constantly exchange substances. In normal flow conditions the mass exchange between blood flow and the near to arterial wall material is kept constant. Wherever and whenever the flow becomes abnormal, the balance between the mass flow of substances from the main blood stream to the arterial wall and vice versa is disrupted [26]. In this case the unbalanced mass transfer probably initiates the arterial wall disease. The mass transfer quantification is crucial in the genesis and the progression of atherosclerosis. Understanding the flow behaviour at the near arterial wall region is of paramount importance. However, mass transport is a highly complex physical phenomenon depending upon biological, chemical and mechanical factors. The action of mass transport occurs within a very thin layer, which is located very close to the endothelium. Thus, any perturbation of the flow in this sensitive region initiates mass flow disturbance to and from endothelium.

### Wall Shear Stress and endothelium

The endothelium permeability is particularly sensitive to local flow perturbations caused from various mechanical factors. It has been proved that endothelial cells are subject to morphological alterations which are activated via changes of WSS magnitude as well as WSS orientation [3]. Elongation of the endothelial cells occurs in regions of high WSS, with the longest cell axis oriented parallel to the flow direction [27]. Conversely, in low WSS regions, the endothelial cells have polygonal shape with no specific orientation. It is possible that these local alterations in endothelial cell morphology have different degrees of permeability for the various blood constituents. Therefore, endothelial cells from regions of various degrees of WSS exhibit various biological and biochemical functions [3]. Since the WSS interacts with the endothelial cells and subsequently with the arterial wall permeability, it is difficult to determine whether it is low WSS or abnormal mass transfer coefficient (permeability coefficient) or both that contribute to atherosclerosis.

### Wall Shear Stress and Low-Density Lipoprotein concentration

In the LMCA bifurcation at regions opposite to the flow divider, predominantly low WSS values occur. High WSS values occur at the "cervix" of the origin and at the flow divider. At the "hips" of the LMCA bifurcation, i.e. at regions opposite to the flow divider, located at the outer walls of LAD and LCxA, low WSS values occur, Fig. 3. High curvatures affect the velocity distribution at vessel bends giving rise to high WSS values. Furthermore, the WSS exhibits high values in the distal regions of the LCA, where the magnitude of the mean flow velocity is relatively higher, due to vessel tapering.

The relationship between WSS and luminal surface concentration of LDL, shown in Fig. 5, indicates that LDL is elevated at locations where WSS is low. As WSS values approach zero, the LDL concentration rapidly increases. A lot of points in the graph may indicate that the LDL concentration is not only WSS dependent. Other flow parameters may seriously affect the flow pattern. Most research workers try to elucidate atherosclerosis by studying the WSS distribution. According to the prevailing theory, low WSS is responsible for atherosclerosis process. Current research analysis results that the combined fluid and mass transport are the keys to understand atherosclerosis. Under the action of the flow the LDL particles move to reach a particular site. Thereafter, it is the contact time and the interaction between LDL and endothelial surface which really matters. The net amount of LDL passing through the arterial wall depends upon the particular physical properties of the wall material (permeability).

The results indicate that LCA tree walls are exposed to a cholesterolemic environment although the applied mass and flow conditions refer to normal human geometry and mass-flow conditions. High luminal surface concentration does not necessarily denote that all the LDL molecules will be transported through the vessel wall. It is the permeability of the endothelium, which mainly determines the final amount passing through the wall.

The results also demonstrate that regions of low WSS values do not necessarily coincide with elevated LDL concentration regions. Low WSS regions exhibit elevated LDL concentrations. At the LMCA bifurcation, the peak LDL concentration is identified just upstream of the low WSS area (regions opposite the flow divider, located at the outer walls of LAD and LCxA) or just downstream of it, Figs. 3 and 4. Concave parts exhibit elevated LDL concentration. On the other hand, convex parts exhibit low LDL concentrations. For the distal LAD vessel the mechanism of elevated LDL concentration might be also attributed to high curvature effects upon the flow pattern and the subsequent spatial accumulation of the LDL. This is clearly shown in strain rate contours, Fig. 9. The paths of the velocities in proximity to the endothelium might be the most important factor for the elevated LDL concentration in either of these two high curvature areas, i.e. areas located either at the vicinity of bifurcations regions or at high curvature regions. At near wall areas, low strain rates occur at regions where elevated LDL concentration is present, Figs 8, 9 and 10. Low strain rates located near to endothelium cover larger areas than the corresponding convex part. Further downstream, increased strain rates are present at the near to wall region, while the distribution of the LDL concentration exhibits low values. The importance of fluid flow and mass transport at the near wall regions becomes apparent. Similar conclusions are reported testing a high curvature bend of RCA with flow separation [14].

As the water infiltration velocity increases, the luminal surface concentration of LDL increases as well. Henceforth, the net amount of LDL mass per second taken up by LCA tree luminal surfaces is also increased. At the low water velocity filtration value, the area averaged LDL concentration over the entire LCA tree is only 3.48% higher than at the inlet. However, at the maximum water infiltration velocity the corresponding value is 23.25%. This indicates that the water infiltration velocity parameter is an important factor for LDL concentration.

High molecular diffusivity results in low LDL concentration on all LCA tree segments. At low molecular diffusivity, the area averaged LDL concentration over the entire LCA tree is 9.39% higher than at the entrance. The concentration attains 4.89% and 3.48% values for the median and maximum molecular diffusivity, respectively. As the blood diffusivity increases, more LDL masses pass from endothelium to the main flow, reducing the LDL concentration at vessel walls.

The assumption that the permeability of the wall to LDL is constant throughout the LCA tree needs to be reconsidered. A new approach will consider the permeability of the wall to be a function of either the WSS or spatial gradient of WSS. Coupled luminal blood flow and transmural fluid flow will be a natural extension of the current research work. Furthermore, arterial wall deformation needs to be taken into consideration, since the flow pattern-mass transport is largely determined by the geometrical configuration of the artery. Arterial movement has to be taken into account since the applied acceleration force will yield a different flow field. It is expected that the flow pattern will be seriously altered, particularly at distal LCA tree regions.

## Conclusion

The present study predicts increased permeation of LDL concentration at specific regions. WSS plays an important role in the wall concentration of the LDL. At low WSS regions, which appear at bifurcations opposite to the flow dividers, the LDL is elevated. Concave sides of the LCA tree exhibit elevated concentration of the LDL in comparison to the convex sides. The LCA tree walls are exposed to cholesterolemic environment although the applied mass and flow conditions refer to normal human geometry and normal mass-flow conditions. With increasing water infiltration velocity the regional area of high luminal surface concentration is increased. The normalized area-averaged concentrations of LDL over the entire LCA tree are, 3.48%, 5.4% and 23.2% higher than that at the entrance, for the low, median and high water infiltration velocities, respectively. Decreased molecular diffusivity increases the LDL concentration. The normalized area-averaged concentrations of LDL over the entire LCA tree are, 9.39%, 4.89% and 3.48% higher than that at the entrance, for the low, median and high water molecular diffusivities, respectively. The degree of increase in luminal surface concentration of LDL is mostly affected from the water infiltration velocity at the vessel wall. Regions of high LDL luminal surface concentration do not always co-locate to the sites of lowest WSS. The paths of the velocity in proximity to the endothelium might be the most important factor for the elevated LDL concentration at areas located either at the vicinity of bifurcations regions or at high curvature regions.

## Notes

## Declarations

### Acknowledgements

The authors would like to thank Vasileios A. Kamperidis for his help in editing this paper.

## Authors’ Affiliations

## References

- Ethier CR:
**Computational modeling of mass transfer and links to atherosclerosis.***Ann Biomed Eng*2002,**30**(4):461–471. Review. 10.1114/1.1468890View ArticleGoogle Scholar - Friedman MH, Deters OJ, Mark FF, Bargeron CB, Hutchins GM:
**Arterial geometry affects hemodynamics. A potential risk factor for atherosclerosis.***Atherosclerosis*1983,**46**(2):225–231. 10.1016/0021-9150(83)90113-2View ArticleGoogle Scholar - Reidy M, Bowyer D:
**The morphology of aortic endothelium in haemodynamically stressed areas associated with branches.***Atherosclerosis*1977,**26**(2):181–194. 10.1016/0021-9150(77)90101-0View ArticleGoogle Scholar - Sun N, Wood NB, Hughes AD, Thom SA, Xu XY:
**Fluid-wall modelling of mass transfer in an axisymmetric stenosis: effects of shear-dependent transport properties.***Ann Biomed Eng*2006,**34**(7):1119–1128. 10.1007/s10439-006-9144-2View ArticleGoogle Scholar - Ogunrinade O, Kameya GT, Trusky GA:
**Effect of fluid shear stress on the permeability of the arterial endothelium.***Ann Biomed Eng*2002,**30**(4):430–460. Review. 10.1114/1.1467924View ArticleGoogle Scholar - Nielsen LB:
**Transfer of low density lipoprotein into the arterial wall and risk of atherosclerosis.***Atherosclerosis*1996,**123**(1–2):1–15. 10.1016/0021-9150(96)05802-9View ArticleGoogle Scholar - Stangeby DK, Ethier CR:
**Computational analysis of coupled blood-wall arterial LDL transport.***J Biomech Eng*2002,**124**(1):1–8. 10.1115/1.1427041View ArticleGoogle Scholar - Glagov S, Zarins C, Giddens DP, Ku DN:
**Hemodynamics and atherosclerosis. Insights and perspectives gained from studies of human arteries.***Arch Pathol Lab Med*1988,**112**(10):1018–1031.Google Scholar - Soulis JV, Farmakis TM, Giannoglou GD, Louridas GE:
**Wall shear stress in normal left coronary artery tree.***J of Biomech*2006,**39**(4):742–749. 10.1016/j.jbiomech.2004.12.026View ArticleGoogle Scholar - Fatouraee N, Deng X, De Champlain A, Guidoin R:
**Concentration polarization of low density lipoproteins (LDL) in the arterial system.***Ann NY Acad Sci*1998,**11**(5):137–146. 10.1111/j.1749-6632.1998.tb10148.xView ArticleGoogle Scholar - Giannoglou GD, Soulis JV, Farmakis TM, Farmakis DM, Louridas GE:
**Hemodynamic factors and the important role of local low static pressure in coronary wall thickening.***Int J Cardiol*2002,**86**(1):27–40. 10.1016/S0167-5273(02)00188-2View ArticleGoogle Scholar - Wada S, Koujiya M, Karino T:
**Theoretical study of the effect of local flow disturbances on the concentration of low-density lipoproteins at the luminal surface of end-to-end anastomosed vessels.***Med Biol Eng Comput*2002,**40**(5):576–587. 10.1007/BF02345458View ArticleGoogle Scholar - Wada S, Karino T:
**Theoretical study on flow-dependent concentration polarization of low-density lipoproteins at the luminal surface of a straight artery.***Biorheology*1999,**36**(3):207–223.Google Scholar - Wada S, Karino T:
**Theoretical prediction of low-density lipoprotein concentration at the luminal surface of an artery with a multiple bend.***Ann Biomed Eng*2002,**30:**778–791. 10.1114/1.1495868View ArticleGoogle Scholar - Soulis JV, Farmakis TM, Giannoglou GD, Hatzizisis IS, Giannakoulas GA, Parcharidis GE, Louridas GE:
**Molecular viscosity in the normal left coronary arterial tree. Is it related to atherosclerosis?***Angiology*2006,**57**(1):33–40. 10.1177/000331970605700105View ArticleGoogle Scholar - Prosi M, Zunino P, Perktold K, Quarteroni A:
**Mathematical and numerical models for transfer of low-density lipoproteins through the arterial walls: a new methodology for the model set up with applications to the study of disturbed luminal flow.***J Biomech*2005,**8**(4):903–917. 10.1016/j.jbiomech.2004.04.024View ArticleGoogle Scholar - Wei L, Wen GB, Tan WC:
**Blood flow and macromolecular transport in curved blood vessels.***Applied Mathematics and Mechanics (English Edition)*2006,**27**(9):1223–1231. 10.1007/s10483-006-0909-yView ArticleGoogle Scholar - Dodge JT Jr, Brown G, Bolson E, Dodge T:
**Intrathoracic spatial location of specified coronary segments on the normal human heart. Applications in quantitative arteriography, assessment of regional risk and contraction, and anatomic display.***Circulation*1988,**78**(5 Pt 1):1167–1180.View ArticleGoogle Scholar - Dodge JT Jr, Brown BG, Bolson E, Dodge HT:
**Lumen diameter of normal human coronary arteries. Influence of age, sex, anatomic variation, and left ventricular hypertrophy or dilation.***Circulation*1992,**86**(1):232–246.View ArticleGoogle Scholar - Kelkar KM, Patankar SV:
**Development of generalized block correction procedures for the solution of discretized Navier-Stokes equations.**In*Creare Inc*. TM-459, Lebanon NH: Creare Inc; 1988.Google Scholar - Kelkar KM:
*Derivation of pressure and continuity equations for Fluent/bfc.*Lebanon NH: Creare Inc; 1989.Google Scholar - Sharma K, Bhat SV:
**Non-Newtonian rheology of leukemic blood and plasma: are n and k parameters of power Law model diagnostic?***Physiol Chem Phys Med NMR*1992,**24:**307–312.Google Scholar - Lei M, Giddens DP, Jones SA, Loth F, Bassiouny H:
**Pulsatile flow in an end-to-side vascular graft model: comparison of computations with experimental data.***J Biomech Eng*2001,**123:**80–87. 10.1115/1.1336145View ArticleGoogle Scholar - Lei M, Kleinstreuer C, Truskey G: A focal stress gradient-dependent mass transfer mechanism for atherogenesis in branching arteries. Med Eng Phys 1996, (4):326–332. 10.1016/1350-4533(95)00045-3Google Scholar
- Murray CD:
**The physiological principle of minimum work. I. The vascular system and the cost of blood volume.***Proc Natl Acad Sci*1926,**12:**207–214. 10.1073/pnas.12.3.207View ArticleGoogle Scholar - Vorp D, Steinman D, Ethier C:
**Computational modeling of arterial biomechanics.***Computing in Science and Engineering (Bioengineeing and Biophysics, IEEE)*2001, 51–63.Google Scholar - Davies PF, Tripathi SC:
**Mechanical stress mechanisms and the cell. An endothelial paradigm.***Circ Res*1993,**72**(2):239–245.View ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.