Influence of material property variability on the mechanical behaviour of carotid atherosclerotic plaques: a 3D fluid-structure interaction analysis

Differential sensitivities to blood pressure variations in internal carotid and intracranial arteries: a numerical approach to stroke prediction

Stroke remains a global health concern, necessitating early prediction for effective management. Atherosclerosis-induced internal carotid and intra cranial stenosis contributes significantly to stroke risk. This study explores the relationship between blood pressure and stroke prediction, focusing on internal carotid artery (ICA) branches: middle cerebral artery (MCA), anterior cerebral artery (ACA), and their role in hemodynamics. Computational fluid dynamics (CFD) informed by the Windkessel model were employed to simulate patient-specific ICA models with introduced stenosis. Central to our investigation is the impact of stenosis on blood pressure, flow velocity, and flow rate across these branches, incorporating Fractional Flow Reserve (FFR) analysis. Results highlight differential sensitivities to blood pressure variations, with M1 branch showing high sensitivity, ACA moderate, and M2 minimal. Comparing blood pressure fluctuations between ICA and MCA revealed heightened sensitivity to potential reverse flow compared to ICA and ACA comparisons, emphasizing MCA's role. Blood flow adjustments due to stenosis demonstrated intricate compensatory mechanisms. FFR emerged as a robust predictor of stenosis severity, particularly in the M2 branch. In conclusion, this study provides comprehensive insights into hemodynamic complexities within major intracranial arteries, elucidating the significance of blood pressure variations, flow attributes, and FFR in stenosis contexts. Subject-specific data integration enhances model reliability, aiding stroke risk assessment and advancing cerebrovascular disease understanding.

Pulse wave imaging of a stenotic artery model with plaque constituents of different stiffnesses: Experimental demonstration in phantoms and fluid-structure interaction simulation

Vulnerable plaques associated with softer components may rupture, releasing thrombotic emboli to smaller vessels in the brain, thus causing an ischemic stroke. Pulse Wave Imaging (PWI) is an ultrasound-based method that allows for pulse wave visualization while the regional pulse wave velocity (PWV) is mapped along the arterial wall to infer the underlying wall compliance. One potential application of PWI is the non-invasive estimation of plaque's mechanical properties for investigating its vulnerability. In this study, the accuracy of PWV estimation in stenotic vessels was investigated by computational simulation and PWI in validation phantoms to evaluate this modality for assessing future stroke risk. Polyvinyl alcohol (PVA) phantoms with plaque constituents of different stiffnesses were designed and constructed to emulate stenotic arteries in the experiment, and the novel fabrication process was described. Finite-element fluid-structure interaction simulations were performed in a stenotic phantom model that matched the geometry and parameters of the experiment in phantoms. The peak distension acceleration of the phantom wall was tracked to estimate PWV. PWVs of 2.57 ms-1, 3.41 ms-1, and 4.48 ms-1 were respectively obtained in the soft, intermediate, and stiff plaque material in phantoms during the experiment using PWI. PWVs of 2.10 ms-1, 3.33 ms-1, and 4.02 ms-1 were respectively found in the soft, intermediate, and stiff plaque material in the computational simulation. These results demonstrate that PWI can effectively distinguish the mechanical properties of plaque in phantoms as compared to computational simulation.

Consistency in Geometry Among Coronary Atherosclerotic Plaques Extracted From Computed Tomography Angiography

Background: The three-dimensional (3D) geometry of coronary atherosclerotic plaques is associated with plaque growth and the occurrence of coronary artery disease. However, there is a lack of studies on the 3D geometric properties of coronary plaques. We aim to investigate if coronary plaques of different sizes are consistent in geometric properties.

Methods: Nineteen cases with symptomatic stenosis caused by atherosclerotic plaques in the left coronary artery were included. Based on attenuation values on computed tomography angiography images, coronary atherosclerotic plaques and calcifications were identified, 3D reconstructed, and manually revised. Multidimensional geometric parameters were measured on the 3D models of plaques and calcifications. Linear and non-linear (i.e., power function) fittings were used to investigate the relationship between multidimensional geometric parameters (length, surface area, volume, etc.). Pearson correlation coefficient (r), R-squared, and p-values were used to evaluate the significance of the relationship. The analysis was performed based on cases and plaques, respectively. Significant linear relationship was defined as R-squared > 0.25 and p < 0.05.

Results: In total, 49 atherosclerotic plaques and 56 calcifications were extracted. In the case-based analysis, significant linear relationships were found between number of plaques and number of calcifications (r = 0.650, p = 0.003) as well as total volume of plaques (r = 0.538, p = 0.018), between number of calcifications and total volume of plaques (r = 0.703, p = 0.001) as well as total volume of calcification (r = 0.646, p = 0.003), and between the total volumes of plaques and calcifications (r = 0.872, p < 0.001). In plaque-based analysis, the power function showed higher R-squared values than the linear function in fitting the relationships of multidimensional geometric parameters. Two presumptions of plaque geometry in different growth stages were proposed with simplified geometric models developed. In the proposed models, the exponents in the power functions of geometric parameters were in accordance with the fitted values.

Conclusion: In patients with coronary artery disease, coronary plaques and calcifications are positively related in number and volume. Different coronary plaques are consistent in the relationship between geometry parameters in different dimensions.

Vessel structural stress mediates aortic media degeneration in bicuspid aortopathy: New insights based on patient-specific fluid-structure interaction analysis

This study aimed to assess the relationship between local mechanical stimuli and regional aortic tissue degeneration using fluid-structure interaction (FSI) analysis in patients with bicuspid aortic valve (BAV) disease. Nine patients underwent ascending aortic replacement were recruited. Tissues were collected to evaluate the pathology features in four regions, greater curvature (GC-region), posterior (P-region), anterior (A-region), and lesser curvature (LC-region). FSI analysis was performed to quantify vessel structural stress (VSS) and flow-induced parameters, including wall shear stress (WSS), oscillatory shear index (OSI), and particle relative residence time (RRT). The correlation between these biomechanical metrics and tissue degeneration was analyzed. Elastin in the medial layer and media thickness were thinnest and the gap between fibers was biggest in the GC-region, followed by the P-region and A-region, while the elastin and media thickness were thickest and the gap smallest in the LC-region. The collagen deposition followed a pattern with the biggest in the GC-region and least in the LC-region. There is a strong negative correlation between mean or peak VSS and elastin thickness in the arterial wall in the GC-region (r = -0.917; p = 0.001 and r = -0.899; p = 0.001), A-region (r = -0.748; p = 0.020 and r = -0.700; p = 0.036) and P-region (r = -0.773; p = 0.014 and r = -0.769; p = 0.015), and between mean VSS and fiber distance in the A-region (r = -0.702, p = 0.035). Moreover, strong negative correlation between mean or peak VSS and media thickness was also observed. No correlation was found between WSS, OSI, and RRT and aortic tissue degeneration in these four regions. These findings indicate that increased VSS correlated with local elastin degradation and aortic media degeneration, implying that it could be a potential biomechanical parameter for a refined risk stratification for patients with BAV.

Stress distribution analysis in healthy and stenosed carotid artery models reconstructed from in vivo ultrasonography

Purpose

This study investigated the accuracy of models reconstructed from ultrasound image processing by comparing the radial displacement waveforms of a subject-specific artery model and evaluated stress changes in the proximal shoulder, throat, and distal shoulder of the plaques depending on the degree of carotid artery stenosis.

Methods

Three groups of subjects (healthy and with less than 50% or more carotid stenosis) were evaluated with ultrasonography. Two-dimensional transverse imaging of the common carotid artery was performed to reconstruct the geometry. A longitudinal view of the same region was recorded to extract the Kelvin viscoelastic model parameters. The pulse pressure waveform and the effective pressure of perivascular tissue were loaded onto the internal and external walls of the model. Effective, circumferential, and principal stresses applied to the plaque throat, proximal shoulder, and distal shoulder in the transverse planes were extracted.

Results

The radial displacement waveforms of the model were closely correlated with those of image processing in all three groups. The mean of the effective, circumferential, and principal stresses of the healthy arteries were 15.01±4.93, 12.97±5.07, and 12.39±2.86 kPa, respectively. As stenosis increased from mild to significant, the mean values of the effective, circumferential, and first principal stresses increased significantly (97%, 74%, and 103% at the plaque throat, respectively) (P<0.05). The minimum effective stress was at the lipid pool. The effective stress in calcified areas was higher than in other parts of the artery wall.

Conclusion

This model can discriminate differences in stresses applied to mildly and severely stenotic plaques.

Consideration of stiffness of wall layers is decisive for patient-specific analysis of carotid artery with atheroma

The paper deals with the impact of chosen geometric and material factors on maximal stresses in carotid atherosclerotic plaque calculated using patient-specific finite element models. These stresses are believed to be decisive for the plaque vulnerability but all applied models suffer from inaccuracy of input data, especially when obtained in vivo only. One hundred computational models based on ex vivo MRI are used to investigate the impact of wall thickness, MRI slice thickness, lipid core and fibrous tissue stiffness, and media anisotropy on the calculated peak plaque and peak cap stresses. The investigated factors are taken as continuous in the range based on published experimental results, only the impact of anisotropy is evaluated by comparison with a corresponding isotropic model. Design of Experiment concept is applied to assess the statistical significance of these investigated factors representing uncertainties in the input data of the model. The results show that consideration of realistic properties of arterial wall in the model is decisive for the stress evaluation; assignment of properties of fibrous tissue even to media and adventitia layers as done in some studies may induce up to eightfold overestimation of peak stress. The impact of MRI slice thickness may play a key role when local thin fibrous cap is present. Anisotropy of media layer is insignificant, and the stiffness of fibrous tissue and lipid core may become significant in some combinations.

Sensitivity of Arterial Hyperelastic Models to Uncertainties in Stress-Free Measurements

Despite major advances made in modeling vascular tissue biomechanics, the predictive power of constitutive models is still limited by uncertainty of the input data. Specifically, key measurements, like the geometry of the stress-free (SF) state, involve a definite, sometimes non-negligible, degree of uncertainty. Here, we introduce a new approach for sensitivity analysis of vascular hyperelastic constitutive models to uncertainty in SF measurements. We have considered two vascular hyperelastic models: the phenomenological Fung model and the structure-motivated Holzapfel-Gasser-Ogden (HGO) model. Our results indicate up to 160% errors in the identified constitutive parameters for a 5% measurement uncertainty in the SF data. Relative margins of errors of up to 30% in the luminal pressure, 36% in the axial force, and over 200% in the stress predictions were recorded for 10% uncertainties. These findings are relevant to the large body of studies involving experimentally based modeling and analysis of vascular tissues. The impact of uncertainties on calibrated constitutive parameters is significant in context of studies that use constitutive parameters to draw conclusions about the underlying microstructure of vascular tissues, their growth and remodeling processes, and aging and disease states. The propagation of uncertainties into the predictions of biophysical parameters, e.g., force, luminal pressure, and wall stresses, is of practical importance in the design and execution of clinical devices and interventions. Furthermore, insights provided by the present findings may lead to more robust parameters identification techniques, and serve as selection criteria in the trade-off between model complexity and sensitivity.

Progressive changes of elastic moduli of arterial wall and atherosclerotic plaque components during plaque development in human coronary arteries

Stiffness of the arterial wall and atherosclerotic plaque components is a determinant of the stress field within plaques, which has been suggested to be an indicator of plaque vulnerability. The diversity and inhomogeneous structure of atherosclerotic lesions complicate the characterization of plaque components. In the present study, stiffness of the arterial wall and atherosclerotic plaque components in human coronary arteries was examined in early and developed atherosclerotic lesions. The force-spectroscopy mode of the atomic force microscope and histological examination were used for determination of elastic moduli at specified locations within samples. Fibrous cap (E = 14.1 ± 3.8 kPa) showed lower stiffness than the fibrous tissue beneath the lipid pool (E = 17.6 ± 3.2 kPa). Calcification zones (E = 96.1 ± 18.8 kPa) and lipid pools (E = 2.7 ± 1.8 kPa) were the stiffest and softest components of atherosclerotic lesions, respectively. The increase of media stiffness (%44.8) and reduction of the elastic modulus of the internal elastic lamina (%28.9) was observed in coronary arteries. Moreover, significant differences were observed between the stiffness of medial layer in diseased parts and free-plaque segments in incomplete plaques of coronary arteries. Our results can be used for better understanding of remodeling mechanisms of the arterial wall with plaque development. Graphical abstract Stiffness alteration of the arterial wall and atherosclerotic plaque components with plaque development in coronary arteries.

Effect of distal thickening and stiffening of plaque cap on arterial wall mechanics

To investigate the effect of longitudinal variations of cap thickness and tissue properties on wall stresses and strains along the atherosclerotic stenosis, stenotic plaque models (uniformly thick, distally thickened, homogenous, and distally stiffened) were constructed and subjected to computational stress analyses with due consideration of fluid-structure interactions (FSI). The analysis considered three different cap thicknesses-45, 65, and 200 μm-and tissue properties-soft, fibrous, and hard. The maximum peak cap stress (PCS) and strain were observed in the upstream throat section and demonstrated increases of the order of 345 and 190%, respectively, as the cap thickness was reduced from 200 to 45 μm in uniformly thick models. Distal stiffening increased PCS in the downstream region; however, the overall effect of this increase was rather small. Distal thickening did not affect maximum PCS and strain values for cap thicknesses exceeding 65 μm; however, a noticeable increase in maximum PCS and corresponding longitudinal variation (or spatial gradient) in stress was observed in the very thin (45-μm-thick) cap. It was, therefore, inferred that existence of a rather thin upstream cap demonstrating distal cap thickening indicates an increased risk of plaque progression and rupture. Graphical Abstract ᅟ.

Impact of Fiber Structure on the Material Stability and Rupture Mechanisms of Coronary Atherosclerotic Plaques

The rupture of an atherosclerotic plaque in the coronary circulation remains the main cause of heart attack. As a fiber-oriented structure, the fiber structure, in particular in the fibrous cap (FC), may affect both loading and material strength in the plaque. However, the role of fiber orientation and dispersion in plaque rupture is unclear. Local orientation and dispersion of fibers were calculated for the shoulder regions, mid FC, and regions with intimal thickening (IT) from histological images of 16 human coronary atherosclerotic lesions. Finite element analysis was performed to assess the effect of these properties on mechanical conditions. Fibers in shoulder regions had markedly reduced alignment (Median [interquartile range] 12.9° [6.6, 18.0], p < 0.05) compared with those in mid FC (6.1° [5.5, 9.0]) and IT regions (6.7° [5.1, 8.6]). Fiber dispersion was highest in shoulders (0.150 [0.121, 0.192]), intermediate in IT (0.119 [0.103, 0.144]), and lowest in mid FC regions (0.093 [0.081, 0.105], p < 0.05). When anisotropic properties were considered, stresses were significantly higher for the mid FC (p = 0.030) and IT regions (p = 0.002) and no difference was found for the shoulder or global regions. Shear (sliding) stress between fibers in each region and their proportion of maximum principal stress were: shoulder (25.8 kPa [17.1, 41.2], 12.4%), mid FC (13.9 kPa [5.8, 29.6], 13.8%), and IT (36.5 kPa [25.9, 47.3], 15.5%). Fiber structure within the FC has a marked effect on principal stresses, resulting in considerable shear stress between fibers. Fiber structure including orientation and dispersion may determine mechanical strength and thus rupture of atherosclerotic plaques.

Numerical simulations of a 3D fluid-structure interaction model for blood flow in an atherosclerotic artery

The inflammatory process of atherosclerosis leads to the formation of an atheromatous plaque in the intima of the blood vessel. The plaque rupture may result from the interaction between the blood and the plaque. In each cardiac cycle, blood interacts with the vessel, considered as a compliant nonlinear hyperelastic. A three dimensional idealized fluid-structure interaction (FSI) model is constructed to perform the blood-plaque and blood-vessel wall interaction studies. An absorbing boundary condition (BC) is imposed directly on the outflow in order to cope with the spurious reflexions due to the truncation of the computational domain. The difference between the Newtonian and non-Newtonian effects is highlighted. It is shown that the von Mises and wall shear stresses are significantly affected according to the rigidity of the wall. The numerical results have shown that the risk of plaque rupture is higher in the case of a moving wall, while in the case of a fixed wall the risk of progression of the atheromatous plaque is higher.

Role of biomechanical forces in the natural history of coronary atherosclerosis

Atherosclerosis remains a major cause of morbidity and mortality worldwide, and a thorough understanding of the underlying pathophysiological mechanisms is crucial for the development of new therapeutic strategies. Although atherosclerosis is a systemic inflammatory disease, coronary atherosclerotic plaques are not uniformly distributed in the vascular tree. Experimental and clinical data highlight that biomechanical forces, including wall shear stress (WSS) and plaque structural stress (PSS), have an important role in the natural history of coronary atherosclerosis. Endothelial cell function is heavily influenced by changes in WSS, and longitudinal animal and human studies have shown that coronary regions with low WSS undergo increased plaque growth compared with high WSS regions. Local alterations in WSS might also promote transformation of stable to unstable plaque subtypes. Plaque rupture is determined by the balance between PSS and material strength, with plaque composition having a profound effect on PSS. Prospective clinical studies are required to ascertain whether integrating mechanical parameters with medical imaging can improve our ability to identify patients at highest risk of rapid disease progression or sudden cardiac events.

MRI-based biomechanical parameters for carotid artery plaque vulnerability assessment

Carotid atherosclerotic plaques are a major cause of ischaemic stroke. The biomechanical environment to which the arterial wall and plaque is subjected to plays an important role in the initiation, progression and rupture of carotid plaques. MRI is frequently used to characterize the morphology of a carotid plaque, but new developments in MRI enable more functional assessment of carotid plaques. In this review, MRI based biomechanical parameters are evaluated on their current status, clinical applicability, and future developments. Blood flow related biomechanical parameters, including endothelial wall shear stress and oscillatory shear index, have been shown to be related to plaque formation. Deriving these parameters directly from MRI flow measurements is feasible and has great potential for future carotid plaque development prediction. Blood pressure induced stresses in a plaque may exceed the tissue strength, potentially leading to plaque rupture. Multi-contrast MRI based stress calculations in combination with tissue strength assessment based on MRI inflammation imaging may provide a plaque stress-strength balance that can be used to assess the plaque rupture risk potential. Direct plaque strain analysis based on dynamic MRI is already able to identify local plaque displacement during the cardiac cycle. However, clinical evidence linking MRI strain to plaque vulnerability is still lacking. MRI based biomechanical parameters may lead to improved assessment of carotid plaque development and rupture risk. However, better MRI systems and faster sequences are required to improve the spatial and temporal resolution, as well as increase the image contrast and signal-to-noise ratio.