Hemodynamics and wall mechanics in human carotid bifurcation and its consequences for atherogenesis: investigation of inter-individual variation.

Finite element simulations of fluid-solid interactions were used to investigate inter-individual variations in flow dynamics and wall mechanics at the carotid artery bifurcation, and its effects on atherogenesis, in three healthy humans (normal volunteers: NV1, NV2, NV4). Subject-specific calculations were based on MR images of structural anatomy and ultrasound measurements of flow at domain boundaries. For all subjects, the largest contiguous region of low wall shear stress (WSS) occurred at the carotid bulb, WSS was high (6-10 Pa) at the apex, and a small localized region of WSS > 10 Pa occurred close to the inner wall of the external carotid artery (ECA). NV2 and NV4 had a "spot" of low WSS distal to the bifurcation at the inner wall of the ECA. Low WSS patches in the common carotid artery (CCA) were contiguous with the carotid bulb low WSS region in NV1 and NV2, but not in NV4. In all three subjects, areas of high oscillatory shear index (OSI) were confined to regions of low WSS. Only NV4 exhibited high levels of OSI on the external adjoining wall of the ECA and CCA. For all subjects, the maximum wall shear stress temporal gradient (WSSTG) was highest at the flow divider (reaching 1,000 Pa/s), exceeding 300 Pa/s at the walls connecting the ECA and CCA, but remaining below 250 Pa/s outside of the ECA. In all subjects, (maximum principle) cyclic strain (CS) was greatest at the apex (NV1: 14%; NV2: 11%; NV4: 6%), and a second high CS region occurred at the ECA-CCA adjoining wall (NV1: 11%, NV2: 9%, NV4: 5%). Wall deformability was included in one simulation (NV2) to verify that it had little influence on the parameters studied. Location and magnitude of low WSS were similar, except for the apex (differences of up to 25%). Wall distensibility also influenced OSI, doubling it in most of the CCA, separating the single high OSI region of the carotid bulb into two smaller regions, and shrinking the ECA internal and external walls' high OSI regions. These observations provide further evidence that significant intra-subject variability exists in those factors thought to impact atherosclerosis.

Characterization of the atherosclerotic carotid bifurcation using MRI, finite element modeling, and histology.

Atherogenesis is known to be associated with the stresses that act on or within the arterial wall. Still, the uneven distribution of atherosclerotic lesions and the impact of vessel remodeling on disease progression are poorly understood. A methodology is proposed to study the correlations between fluid dynamic parameters and histological markers of atherosclerosis. Trends suggested by preliminary data from four patients with advanced carotid bifurcation arterial disease are examined and compared to hypotheses in the literature. Four patients were scanned using MRI and ultrasound, and subsequently underwent carotid endarterectomy. For each patient. a geometric model and a numerical mesh were constructed from MR data, and velocity boundary conditions established. Computations yield values for average wall shear stress (WSS), maximum wall shear stress temporal gradient (WSSTG), and Oscillatory Shear Index (OSI). Following surgery, the excised plaques were sectioned, stained for smooth muscle cells (SMC), macrophages (M phi), lipid (LIP), and collagen (COL), and analyzed quantitatively. Correlations attempted between the various fluid dynamic variables and the biological markers were interesting but inconclusive. Tendencies of WSSTG and WSS to correlate negatively with M phi and LIP, and positively with COL and SMC, as well as tendencies of OSI to correlate positively with Mphi and LIP and negatively with COL and SMC, were observed. These trends agree with hypotheses in the literature, which are based on ex vivo and in vitro experimental studies.