Does the principle of minimum work apply at the carotid bifurcation: a retrospective cohort study
© Beare et al; licensee BioMed Central Ltd. 2011
Received: 3 May 2011
Accepted: 24 August 2011
Published: 24 August 2011
There is recent interest in the role of carotid bifurcation anatomy, geometry and hemodynamic factors in the pathogenesis of carotid artery atherosclerosis. Certain anatomical and geometric configurations at the carotid bifurcation have been linked to disturbed flow. It has been proposed that vascular dimensions are selected to minimize energy required to maintain blood flow, and that this occurs when an exponent of 3 relates the radii of parent and daughter arteries. We evaluate whether the dimensions of bifurcation of the extracranial carotid artery follow this principle of minimum work.
This study involved subjects who had computed tomographic angiography (CTA) at our institution between 2006 and 2007. Radii of the common, internal and external carotid arteries were determined. The exponent was determined for individual bifurcations using numerical methods and for the sample using nonlinear regression.
Mean age for 45 participants was 56.9 ± 16.5 years with 26 males. Prevalence of vascular risk factors was: hypertension-48%, smoking-23%, diabetes-16.7%, hyperlipidemia-51%, ischemic heart disease-18.7%.
The value of the exponent ranged from 1.3 to 1.6, depending on estimation methodology.
The principle of minimum work (defined by an exponent of 3) may not apply at the carotid bifurcation. Additional factors may play a role in the relationship between the radii of the parent and daughter vessels.
There has been recent interest in the role of carotid artery anatomy, geometry and hemodynamic factors in the pathogenesis of carotid artery atherosclerosis [1–6]. The anatomy and geometry at the carotid bifurcation within the same individuals [1, 7] and between the sexes  vary greatly. The anatomy [7, 8] and geometry  of the carotid bifurcation have large influence on vortex flow at the carotid sinus. These studies support observations that plaques form preferentially at such sites as carotid artery bifurcation (extracranial site) and the carotid artery siphon (intracranial site), rather than randomly in the carotid artery. This may provide an explanation for the asymmetrical nature of carotid artery stenosis  within the same individual despite exposure to the same factors . In addition to this, there may be a role for interplay between anatomy, hemodynamic shear forces and traditional risk factors.
This work on carotid artery geometry has been inspired by the writings of Murray in 1926  and Rossitti and Lofgren in 1993 [12, 13]. A law on harmonisation of the vascular dimensions to minimize total energy cost by balancing energy lost to shear stress (frictional force) and the energy required to sustain the total blood volume (metabolic cost) throughout the vascular tree originated with the work of Murray in the 1920s and is known as the principle of minimum work  [see additional file 1]. The influence of this can be seen in the scaling law in biology describing how nutrients are transported through branching fractal network of vessels  and the scaling law of vascular trees such as the coronary tree . In simple terms, this law predicts a cubic relationship (n = 3) between the radii of the parent and daughter vessels. Murray made several assumptions: 1- blood flow is non-pulsatile, has parabolic laminar flow profile and exerts constant shear stress on the arterial wall throughout the vascular tree, 2- the vessel is straight and rigid, and 3-the blood flow is described by Poiseuille's equation for flow in tubes . The implication is that almost 94% of the frictional losses can be saved by doubling the radius of the tube while keeping length and other variables constant . Rossitti and Lofgren examine the principle of minimum work in the intracranial internal carotid artery (ICA) to explore the importance of vascular dimensions on aneurysm formation. These investigators proposed that the vascular dimensions of the cerebral arteries follow the principle of minimum work  and that the bifurcations of the cerebral arteries appear to be optimized to avoid increased hemodynamic stresses . These postulates have influenced recent works on hemodynamic factors in the development of carotid artery atherosclerotic disease  and intracranial aneurysms [17, 18].
Recent studies provide evidence suggesting that the principle of minimum work applies at branching points several orders away from the aorta [13, 17]. While the principle of minimum work may apply at the intracranial portion of the carotid artery, its applicability at the extracranial carotid artery bifurcation has not been tested . In this study, we evaluate whether the bifurcation dimensions of the extracranial carotid artery follow this principle of minimum work. An understanding of this issue may shed light on the potential role of anatomy and geometrical factors in the development of carotid artery atherosclerosis.
This study involved patients who had CT angiography at our institution from 2006 to 2007 who had attended the Stroke and Vascular clinics. The Southern Health Ethics committee approved the study
Axial CT angiography scans were performed on a fourth generation CT Light speed scanner (General Electric Medical Systems, Milwaukee, WI) using a test bolus to determine the time of maximal contrast arrival in the artery. The field of view covered the aortic arch to the circle of Willis. The images have voxel resolution of 2 mm × 0.5 mm × 0.5 mm. The CT imaging uses a tube voltage of 120 kV, an effective current of 200 mA and 75 ml of intravenous contrast delivered through a power injector at 3 ml/sec.
The method for segmentation of carotid arteries has been described previously . In brief, a set of marker points was required to define the object of interest. Two marker points were used for each of the CCA, ICA and ECA and a sketch of the arterial tree was computed automatically using minimal cost paths. This sketch was used to initialize a watershed transform which produced a three dimensional mask of the arterial tree. In a validation study, the segmentation showed an intra-class correlation of 0.96 compared to manual segmentation by an expert neuroradiologist.
Measurements of arterial radii and area
Determining the exponent of the power law (Analysis)
where r CCA is the radius of the common carotid artery (CCA), r ICA is the radius of the internal carotid artery (ICA) and r ECA is the radius of the external carotid artery (ECA). The value for the exponent n derived by Murray using the energy minimization approach is 3. A brief derivation of this power law is provided in the additional file 1. To test the applicability of the energy minimization design principle to the carotid bifurcation, we solved Equation 1 for each bifurcation using numerical root finding methods. We also used nonlinear regression to estimate the optimum exponent for our study subjects. In addition, we evaluated the average change in ICA radius necessary to satisfy the power law. Nonlinear regression was performed using a Gauss-Newton algorithm with the nls function in the R statistics package . Root estimation was performed using the zeroin Netlib algorithm  provided by the uniroot function in R.
Murray's model assumed straight arteries; hence we performed the analysis in a subset of patients with non-tortuous arteries, defined as tortuosity value less than 0.1. The tortuosity is defined as (L/D)-1, where L is the length along centreline and D is the Euclidean distance between two end points of the artery concerned .
Artery radius ratios and standard deviations for arteries and bifurcation area ratio
Artery size measurement method
Bifurcation area ratio
Estimated exponents and 95% confidence interval for all patients and arteries with low tortuosity using artery size measurements A and B
Artery size measurement method
Tortuosity < 0.1
In this analysis of the extracranial carotid bifurcation, we have found that the value of the exponent, n, in the power law describing the sizes of the parent and daughter arteries is between 1.2 and 1.6, depending on estimation method, rather than the value of 3 predicted using Murray's energy minimization argument. The estimates obtained for individual bifurcations, using numerical root finding methods, and for the subject populations, using nonlinear regression, were consistent, with very few subjects having exponents close to 3. The differences in estimate between nonlinear regression and numerical root-finding approaches remained despite removal of possible outliers and tests with alternative solvers and are most likely attributable to the weighting used during nonlinear fitting. Figures 2 and 4 suggest that the estimate based on mean of individual estimates is more appropriate in this dataset, with the line of best fit appearing better centred on the data. Only small differences in estimated exponent values were observed when different methods for estimating artery radius were used, or when the analysis was restricted to a subset of the study population with less tortuous arteries. The different methods for estimating artery diameter are likely to have a bigger impact in diseased arteries where artery cross-sections are less circular. The values obtained for n for individual arteries varied considerably and were not centred on either 2 (corresponding to preservation of cross sectional area) or 3. Changes in ICA radius necessary to satisfy the power law while holding CCA and ECA radius constant were of the order of 20% for n = 2 and 40% for n = 3. These finding suggest that additional factors such as variation in flow velocity during cardiac cycle, shape of the carotid sinus and vessel elasticity contribute to the relative sizes of arteries at the carotid bifurcation.
In this study, the measurements of arterial radii were taken from the arterial lumen due to the use of CT angiography images. The radiological contrast agent in CT angiography depicts the lumen rather than the exterior of the arterial wall. However, this use of luminal radii is similar to the measurement of radii in Rossitti and Lofgren's studies and is unlikely to have affected the analysis [12, 13]. Our use of three-dimensional angiography data made it possible to measure the cross sectional area of arteries perpendicular to the central axis. This measurement is likely to be more difficult to perform consistently using conventional two-dimensional angiography data, potentially leading to errors in arterial radii measures. Such errors may explain dramatically different exponents for the power law obtained for the intracranial artery - 2.9 (two-dimensional angiogram data)  and 1.7 (three-dimensional angiogram data) .
Our finding that the exponent of the power law is much less than 2 occurs despite efforts to perform the analysis in relatively straight arteries that are appropriate for Murray's model and using radius derived from cross sectional area. Despite removal of subjects with carotid artery stenosis, it is possible that the use of CT angiography images from patients attending vascular and stroke clinics rather than from the community may have affected our results. Previous findings of exponents less than 2 in studies using three-dimensional techniques  provide some re-assurance that this choice of subject did not greatly bias our results.
Implications for Murray's Law
The exponent of the power law at the carotid bifurcation in this study is different from intracranial component [12, 17]. This finding supports the suggestion that the exponent of the power law is not constant throughout the vascular tree . In line with this possibility, investigators suggested that the exponent of the power law is 2 in the major branches close to the aortic arch, 2.5-3 in coronary, 2.9 in MCA bifurcation, to 3 in arterioles [17, 25–27]. In-vivo measurements of wall sheer stress (WSS)  show that the assumption of constant WSS throughout the body made by Murray is unrealistic and hence the variation in power law exponent is not surprising. An exponent of 2 corresponds to preservation of area and therefore constant flow velocity into and out of the bifurcation. Exponents greater than 2 imply a decrease in flow velocity leaving the bifurcation while exponents less than 2 imply an increasing flow velocity leaving the bifurcation.
Murray's law does appear to be useful in describing size of small vessels (canine arteries with radius of 0.159 cm and capillaries, radius 3.5 × 10-4 cm) as illustrated in the original calculations. The law also successfully predicts the sizing of human coronary arteries (mean radius 1.44 mm)  and the MCA (mean radius 1.2 mm) . However the law does not apply in larger arteries.
These variations suggest that any energy minimization argument used to explain vascular size throughout the body and across species should include more terms and be more complicated than the shear stress and volume argument proposed by Murray. Possibilities include energy terms relating to energy loss due to reflections at bifurcations that may become more significant in areas where the blood flow is more pulsatile or where the flexibility of vessel walls interacts significantly with the nature of flow.
This study examined the relationship between radiuses of arteries at the extracranial carotid bifurcation using CTA imaging in 45 subjects. The power law of Murray  did not apply to this data. More complex terms are therefore needed if minimization of energy arguments can be used to describe carotid bifurcation geometry.
This work was supported by the Monash-CSIRO collaborative research support scheme (CF07/4242) and a Pfizer Cardiovascular Lipid Grant. RB was supported by Monash University, Murdoch Childrens Research Institute and the Victorian Government Operational Infrastructure Support Program. GD was supported by and MR was a student at Monash University. VS was supported by Monash University and Southern Health. Southern Health supported WC and TP. MS and JH were supported by CSIRO Mathematics, Informatics and Statistics.
- Schulz UGR, Rothwell PM: Major Variation in Carotid Bifurcation Anatomy A Possible Risk Factor for Plaque Development?. Stroke. 2001, 32: 2522-2529. 10.1161/hs1101.097391.View ArticlePubMedGoogle Scholar
- Hademenos GJ, Massoud TF: Biophysical mechanisms of stroke. Stroke. 1997, 28 (10): 2067-2077. 10.1161/01.STR.28.10.2067.View ArticlePubMedGoogle Scholar
- Schulz UGR, Rothwell PM: Sex differences in carotid bifurcation anatomy and the distribution of atherosclerotic plaque. Stroke. 2001, 32: 1525-1531. 10.1161/01.STR.32.7.1525.View ArticlePubMedGoogle Scholar
- Thomas JB, Antiga L, Che SL, Milner JS, Steinman DA, Spence JD, Rutt BK, Steinman DA: Variation in the carotid bifurcation geometry of young versus older adults: implications for geometric risk of atherosclerosis. Stroke. 2005, 36: 2450-2456. 10.1161/01.STR.0000185679.62634.0a.View ArticlePubMedGoogle Scholar
- Lee SW, Antiga L, Spence JD, Steinman DA: Geometry of the carotid bifurcation predicts its exposure to disturbed flow. Stroke. 2008, 39: 2341-2347. 10.1161/STROKEAHA.107.510644.View ArticlePubMedGoogle Scholar
- Sitzer M, Puac D, Buehler A, Steckel DA, von Kegler S, Markus HS, Steinmetz H: Internal carotid artery angle of origin: a novel risk factor for early carotid atherosclerosis. Stroke. 2003, 34 (4): 950-955. 10.1161/01.STR.0000060895.38298.C4.View ArticlePubMedGoogle Scholar
- Fisher M, Fieman S: Geometrical factors of the bifurcation in carotid atherogenesis. Stroke. 1990, 21: 227-271.Google Scholar
- Karino T, Goldsmith HL: Particle flow behavior in models of branching vessel. II. Effects of branching angle and diameter ratio on flow patterns. Biorrheology. 1985, 22: 87-104.Google Scholar
- Gnasso A, Concetta Irace C, Carallo C, De Franceschi MS, Motti C, Mattioli PL, Pujia A: In Vivo Association Between Low Wall Shear Stress and Plaque in Subjects With Asymmetrical Carotid Atherosclerosis. Stroke. 1997, 28: 993-998. 10.1161/01.STR.28.5.993.View ArticlePubMedGoogle Scholar
- Kannel WB: Contributions of the Framingham Study to the conquest of coronary artery disease. Am J Cardiol. 1988, 62: 1109-1112. 10.1016/0002-9149(88)90558-9.View ArticlePubMedGoogle Scholar
- Murray CD: The Physiological principle of minimum work. I. The vascular system abd the cost of blood voume. PROC N A S. 1926, 12: 207-214. 10.1073/pnas.12.3.207.View ArticleGoogle Scholar
- Rossitti S, Lofgren J: Vascular dimensions of the cerebral arteries follow the principle of minimum work. Stroke. 1993, 24 (3): 371-377. 10.1161/01.STR.24.3.371.View ArticlePubMedGoogle Scholar
- Rossitti S, Lofgren J: Optimality principles and flow orderliness at the branching points of cerebral arteries. Stroke. 1993, 24: 1029-1032. 10.1161/01.STR.24.7.1029.View ArticlePubMedGoogle Scholar
- West GB, Brown JH, Enquist BJ: A general model for the origin of allometric scaling laws in biology. Science. 1997, 276: 122-126. 10.1126/science.276.5309.122.View ArticlePubMedGoogle Scholar
- Kassab GS: Scaling laws of vascular trees: of form and function. Am J Physiol Heart Circ Physiol. 2006, 290: H894-H903.View ArticlePubMedGoogle Scholar
- Zamir M: Shear forces and blood vessel radii in the cardiovascular system. J Gen Physiol. 1977, 69: 449-461. 10.1085/jgp.69.4.449.View ArticlePubMedGoogle Scholar
- Ingebrigtsen T, Morgan MK, Faulder K, Ingebrigtsen L, Sparr T, Schirmer H: Bifurcation geometry and the presence of cerebral artery aneurysms. J Neurosurg. 2004, 101: 108-113. 10.3171/jns.2004.101.1.0108.View ArticlePubMedGoogle Scholar
- Bor AS, Velthuis BK, Majoie CB, Rinkel GJ: Configuration of intracranial arteries and development of aneurysms: a follow-up study. Neurology. 2008, 70 (9): 700-705. 10.1212/01.wnl.0000302176.03551.35.View ArticlePubMedGoogle Scholar
- Beare R, Das G, Chong W, Ren M, Srikanth V, Phan TG: Segmentation of carotid arteries in CTA images. International Conference on Digital Image Computing: Techniques and Applications (DICTA). 2010, Sydney, AustraliaGoogle Scholar
- Bertrand G, Malandain G: A new characterization of three-dimensional simple points. Pattern Recognition Letters. 1994, 15: 169-175. 10.1016/0167-8655(94)90046-9.View ArticleGoogle Scholar
- R_Development_Core_Team: R: A language and environment for statistical computing. 2011, R Foundation for Statistical Computing, Vienna, AustriaGoogle Scholar
- R B: Algorithms for Minimization without Derivatives. 1973, Englewood Cliffs, NJ: Prentice-HallGoogle Scholar
- Schulz UG, Rothwell PM: Major variation in carotid bifurcation anatomy: a possible risk factor for plaque development?. Stroke. 2001, 32 (11): 2522-2529. 10.1161/hs1101.097391.View ArticlePubMedGoogle Scholar
- Thomas JB, Antiga L, Che SL, Milner JS, Steinman DA, Spence JD, Rutt BK: Variation in the carotid bifurcation geometry of young versus older adults: implications for geometric risk of atherosclerosis. Stroke. 2005, 36 (11): 2450-2456. 10.1161/01.STR.0000185679.62634.0a.View ArticlePubMedGoogle Scholar
- Reneman RS, Vink H, Hoeks APG: Wall shear stress revisited. Artery Research. 2009, 3: 73-78. 10.1016/j.artres.2009.02.005.View ArticleGoogle Scholar
- Hutchins GM, Miner MM, Boitnott JK: Vessel caliber and branch-angle of human coronary artery branch-points. Circ Res. 1976, 38: 572-576.View ArticlePubMedGoogle Scholar
- Hahn JY, Gwon HC, Kwon SU, Choi SH, Choi JH, Lee SH, Hong KP, Park JE, Kim DK: Comparison of vessel geometry in bifurcation between normal and diseased segments: intravascular ultrasound analysis. Atherosclerosis. 2008, 201: 326-331. 10.1016/j.atherosclerosis.2008.02.021.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/11/17/prepub
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.