Search

Atherosclerotic plaque rupture leading to stroke is the major cause of long-term disability as well as the third most common cause of mortality. Image-based computational models have been introduced seeking critical mechanical indicators, which may be used for plaque vulnerability assessment. This study extends the previous 2D critical stress concept to 3D by using in vivo magnetic resonance image (MRI) data of human atherosclerotic carotid plaques and 3D fluid-structure interaction (FSI) models to: identify 3D critical plaque wall stress (CPWS) and critical flow shear stress (CFSS) and to investigate their associations with plaque rupture. In vivo MRI data of carotid plaques from 18 patients scheduled for endarterectomy were acquired using histologically validated multicontrast protocols. Of the 18 plaques, histology-confirmed that six had prior rupture (group 1) as evidenced by presence of ulceration. The remaining 12 plaques (group 2) contained no rupture. The 3D multicomponent FSI models were constructed for each plaque to obtain 3D plaque wall stress (PWS) and flow shear stress (FSS) distributions. Three-dimensional CPWS and CFSS, defined as maxima of PWS and FSS from all vulnerable sites, were determined for each plaque to investigate their association with plaque rupture. Slice-based critical PWS and FSS were also calculated for all slices for more detailed analysis and comparison. The mean 3D CPWS of group 1 was 263.44 kPa, which was 100% higher than that from group 2 (132.77, p =0.03984). Five of the six ruptured plaques had 3D CPWS sites, matching the histology-confirmed rupture sites with an 83% agreement. Although the mean 3D CFSS (92.94 dyn/[cm.sup.2]) for group 1 was 76% higher than that for group 2 (52.70 dyn/[cm.sup.2]), slice-based CFSS showed no significant difference between the two groups. Only two of the six ruptured plaques had 3D CFSS sites matching the histology-confirmed rupture sites with a 33% agreement. CFSS had a good correlation with plaque stenosis severity ([R.sup.2] =0.40 with an exponential function fitting 3D CFSS data). This in vivo MRI pilot study using plaques with and without rupture demonstrates that 3D critical plaque wall stress values are more closely associated with atherosclerotic plaque rupture then critical flow shear stresses. Critical wall stress values may become indicators of high risk sites of rupture. Future work with a larger population will establish a possible CPWS-based plaque vulnerability classification. [DOI: 10.1115/1.4001028] Keywords: carotid artery, atherosclerosis, vulnerable plaque, stroke, MRI, plaque rupture, fluid-structure interaction

Heart attack and stroke are often caused by atherosclerotic plaque rupture, which happens without warning most of the time. Magnetic resonance imaging (MRl)-based atherosclerotic plaque models with fluid-structure interactions (FSIs) have been introduced to perform flow and stress/strain analysis and identify possible mechanical and morphological indices for accurate plaque vulnerability assessment. For coronary arteries, cyclic bending associated with heart motion and anisotropy of the vessel walls may have significant influence on flow and stress/strain distributions in the plaque. FSI models with cyclic bending and anisotropic vessel properties for coronary plaques are lacking in the current literature. In this paper, cyclic bending and anisotropic vessel properties were added to 3D FSI coronary plaque models so that the models would be more realistic for more accurate computational flow and stress/strain predictions. Six computational models using one ex vivo MRI human coronary plaque specimen data were constructed to assess the effects of cyclic bending, anisotropic vessel properties, pulsating pressure, plaque structure, and axial stretch on plaque stress/strain distributions. Our results indicate that cyclic bending and anisotropic properties may cause 50-800% increase in maximum principal stress (Stress-[P.sub.1]) values at selected locations. The stress increase varies with location and is higher when bending is coupled with axial stretch, nonsmooth plaque structure, and resonant pressure conditions (zero phase angle shift). Effects of cyclic bending on flow behaviors are more modest (9.8% decrease in maximum velocity, 2.5% decrease in flow rate, 15% increase in maximum flow shear stress). Inclusion of cyclic bending, anisotropic vessel material properties, accurate plaque structure, and axial stretch in computational FSI models should lead to a considerable improvement of accuracy of computational stress/strain predictions for coronary plaque vulnerability assessment. Further studies incorporating additional mechanical property data and in vivo MRI data are needed to obtain more complete and accurate knowledge about flow and stress/strain behaviors in coronary plaques and to identify critical indicators for better plaque assessment and possible rupture predictions. [DOI: 10.1115/1.3127253] Keywords: coronary artery, cardiovascular, cyclic bending, fluid-structure interaction, blood flow, atherosclerotic plaque rupture

Fiber structures and pathological features, e.g., inflammation and glycosaminoglycan (GAG) deposition, are the primary determinants of aortic mechanical properties which are associated with the development of an aneurysm. This study is designed to quantify the association of tissue ultimate strength and extensibility with the structural percentage of different components, in particular, GAG, and local fiber orientation. Thoracic aortic aneurysm (TAA) tissues from eight patients were collected. Ninety-six tissue strips of thickened intima, media, and adventitia were prepared for uni-extension tests and histopathological examination. Area ratios of collagen, elastin, macrophage and GAG, and collagen fiber dispersion were quantified. Collagen, elastin, and GAG were layer-dependent and the inflammatory burden in all layers was low. The local GAG ratio was negatively associated with the collagen ratio (r2 = 0.173, p < 0.05), but positively with elastin (r2 = 0.037, p < 0.05). Higher GAG deposition resulted in larger local collagen fiber dispersion in the media and adventitia, but not in the intima. The ultimate stretch in both axial and circumferential directions was exclusively associated with elastin ratio (axial: r2 = 0.186, p = 0.04; circumferential: r2 = 0.175, p = 0.04). Multivariate analysis showed that collagen and GAG contents were both associated with ultimate strength in the circumferential direction, but not with the axial direction (collagen: slope = 27.3, GAG: slope = -18.4, r2 = 0.438, p = 0.002). GAG may play important roles in TAA material strength. Their deposition was found to be associated positively with the local collagen fiber dispersion and negatively with ultimate strength in the circumferential direction., (Copyright © 2022 by ASME; reuse license CC-BY 4.0.)

Arteries can be considered as layered composite material. Experimental data on the stiffness of human atherosclerotic carotid arteries and their media and adventitia layers are very limited. This study used uniaxial tests to determine the stiffness (tangent modulus) of human carotid artery sections containing American Heart Association type II and III lesions. Axial and circumferential oriented adventitia, media, and full thickness specimens were prepared from six human carotid arteries (total tissue strips: 71). Each artery yielded 12 specimens with two specimens in each of the following six categories; axial full thickness, axial adventitia (AA), axial media (AM), circumferential full thickness, circumferential adventitia (CA), and circumferential media (CM). Uniaxial testing was performed using Inspec 2200 controlled by software developed using labview. The mean stiffness of the adventitia was 3570 ± 667 and 2960 ± 331 kPa in the axial and circumferential directions, respectively, while the corresponding values for the media were 1070 ± 186 and 1800 ± 384 kPa. The adventitia was significantly stiffer than the media in both the axial (p = 0.003) and circumferential (p = 0.010) directions. The stiffness of the full thickness specimens was nearly identical in the axial (1540 ± 186) and circumferential (1530 ± 389 kPa) directions. The differences in axial and circumferential stiffness of media and adventitia were not statistically significant.

In vivo magnetic resonance image (MRI)-based computational models have been introduced to calculate atherosclerotic plaque stress and strain conditions for possible rupture predictions. However, patient-specific vessel material properties are lacking in those models, which affects the accuracy of their stress/strain predictions. A noninvasive approach of combining in vivo Cine MRI, multicontrast 3D MRI, and computational modeling was introduced to quantify patient-specific carotid artery material properties and the circumferential shrinkage rate between vessel in vivo and zero-pressure geometries. In vivo Cine and 3D multicontrast MRI carotid plaque data were acquired from 12 patients after informed consent. For each patient, one nearly-circular slice and an iterative procedure were used to quantify parameter values in the modified Mooney-Rivlin model for the vessel and the vessel circumferential shrinkage rate. A sample artery slice with and without a lipid core and three material parameter sets representing stiff, median, and soft materials from our patient data were used to demonstrate the effect of material stiffness and circumferential shrinkage process on stress/strain predictions. Parameter values of the Mooney-Rivlin models for the 12 patients were quantified. The effective Young's modulus (YM, unit: kPa) values varied from 137 (soft), 431 (median), to 1435 (stiff), and corresponding circumferential shrinkages were 32%, 12.6%, and 6%, respectively. Using the sample slice without the lipid core, the maximum plaque stress values (unit: kPa) from the soft and median materials were 153.3 and 96.2, which are 67.7% and 5% higher than that (91.4) from the stiff material, while the maximum plaque strain values from the soft and median materials were 0.71 and 0.293, which are about 700% and 230% higher than that (0.089) from the stiff material, respectively. Without circumferential shrinkages, the maximum plaque stress values (unit: kPa) from the soft, median, and stiff models were inflated to 330.7, 159.2, and 103.6, which were 116%, 65%, and 13% higher than those from models with proper shrinkage. The effective Young's modulus from the 12 human carotid arteries studied varied from 137 kPa to 1435 kPa. The vessel circumferential shrinkage to the zero-pressure condition varied from 6% to 32%. The inclusion of proper shrinkage in models based on in vivo geometry is necessary to avoid over-estimating the stresses and strains by up 100%. Material stiffness had a greater impact on strain (up to 700%) than on stress (up to 70%) predictions. Accurate patient-specific material properties and circumferential shrinkage could considerably improve the accuracy of in vivo MRI-based computational stress/strain predictions.

High mechanical stress condition over the fibrous cap (FC) has been widely accepted as a contributor to plaque rupture. The relationships between the stress, lumen curvature, and FC thickness have not been explored in detail. In this study, we investigate lumen irregularity-dependent relationships between mechanical stress conditions, local FC thickness (LT(FC)), and lumen curvature (LC(lumen)). Magnetic resonance imaging slices of carotid plaque from 100 patients with delineated atherosclerotic components were used. Two-dimensional structure-only finite element simulations were performed for the mechanical analysis, and maximum principal stress (stress-P₁) at all integral nodes along the lumen was obtained. LT(FC) and LC(lumen) were computed using the segmented contour. The lumen irregularity (L-δir) was defined as the difference between the largest and the smallest lumen curvature. The results indicated that the relationship between stress-P₁, LT(FC), and LC(lumen) is largely dependent on L-δir. When L-δir ≥ .31 (irregular lumen), stress-P₁ strongly correlated with lumen curvature and had a weak/no correlation with local FC thickness, and in 73.4% of magnetic resonance (MR) slices, the critical stress (maximum of stress-P₁ over the diseased region) was found at the site where the lumen curvature was large. When L-δir ≤ 0.28 (relatively round lumen), stress-P₁ showed a strong correlation with local FC thickness but weak/no correlation with lumen curvature, and in 71.7% of MR slices, the critical stress was located at the site of minimum FC thickness. Using lumen irregularity as a method of identifying vulnerable plaque sites by referring to the lumen shape is a novel and simple method, which can be used for mechanics-based plaque vulnerability assessment.