Local maximal stress hypothesis and computational plaque vulnerability index for atherosclerotic plaque assessment

Noninvasive assessment of carotid plaque with subharmonic aided pressure estimation from a US contrast agent: A preliminary study

Stroke is closely associated with carotid plaques. The assessment of carotid plaque is still the key issue of stroke prevention in clinical practice. This prospective cross-sectional study included patients with carotid plaque evaluated by ultrasonography (US). The intima-media thickness (IMT), lumen stenosis severity, thickness, and length of carotid plaque were measured by the routine US, and the amplitudes of subharmonics in the upstream shoulder, top, and downstream shoulder of all plaques and corresponding lumens were observed by Subharmonic Aided Pressure Estimation (SHAPE) US examination from the US contrast agent perflubutane microbubbles (Sonazoid), which analyzed the clinical parameters of patients, the subharmonic amplitude characteristics of all plaques and lumens, and the parameter differences between the ischemic stroke (IS) group and control group. From May 2021 to February 2022, 46 carotid plaques of 23 patients were included. For plaques, the subharmonic amplitude in the plaque (-60.52 ± 4.46) was lower than that in the opposing level lumen (-56.82 ± 5.68 dB), the subharmonic gradient across the plaque cap was negatively correlated with plaque thickness (r = -0.51, p < 0.001), and with the lumen stenosis severity (r = -0.42, p = 0.003). The median IMT of the IS group was thicker than the control group. The subharmonic gradient of the intraplaque of the IS group was larger than the control group (p = 0.004). In this analysis, we use the receiver operating characteristic (ROC) curve to establish the cutoff value of the difference to predict a new monitoring method for plaque without invasion to predict IS. It still needs a large-scale study with long-term follow-up to validate these findings.

Combining IVUS + OCT Data, Biomechanical Models and Machine Learning Method for Accurate Coronary Plaque Morphology Quantification and Cap Thickness and Stress/Strain Index Predictions

Assessment and prediction of vulnerable plaque progression and rupture risk are of utmost importance for diagnosis, management and treatment of cardiovascular diseases and possible prevention of acute cardiovascular events such as heart attack and stroke. However, accurate assessment of plaque vulnerability assessment and prediction of its future changes require accurate plaque cap thickness, tissue component and structure quantifications and mechanical stress/strain calculations. Multi-modality intravascular ultrasound (IVUS), optical coherence tomography (OCT) and angiography image data with follow-up were acquired from ten patients to obtain accurate and reliable plaque morphology for model construction. Three-dimensional thin-slice finite element models were constructed for 228 matched IVUS + OCT slices to obtain plaque stress/strain data for analysis. Quantitative plaque cap thickness and stress/strain indices were introduced as substitute quantitative plaque vulnerability indices (PVIs) and a machine learning method (random forest) was employed to predict PVI changes with actual patient IVUS + OCT follow-up data as the gold standard. Our prediction results showed that optimal prediction accuracies for changes in cap-PVI (C-PVI), mean cap stress PVI (meanS-PVI) and mean cap strain PVI (meanSn-PVI) were 90.3% (AUC = 0.877), 85.6% (AUC = 0.867) and 83.3% (AUC = 0.809), respectively. The improvements in prediction accuracy by the best combination predictor over the best single predictor were 6.6% for C-PVI, 10.0% for mean S-PVI and 8.0% for mean Sn-PVI. Our results demonstrated the potential using multi-modality IVUS + OCT image to accurately and efficiently predict plaque cap thickness and stress/strain index changes. Combining mechanical and morphological predictors may lead to better prediction accuracies.

Non-Invasive Modalities in the Assessment of Vulnerable Coronary Atherosclerotic Plaques

Coronary atherosclerosis is a complex, multistep process that may lead to critical complications upon progression, revolving around plaque disruption through either rupture or erosion. Several high-risk features are associated with plaque vulnerability and may add incremental prognostic information. Although invasive imaging modalities such as optical coherence tomography or intravascular ultrasound are considered to be the gold standard in the assessment of vulnerable coronary atherosclerotic plaques (VCAPs), contemporary evidence suggests a potential role for non-invasive methods in this context. Biomarkers associated with deleterious pathophysiologic pathways, including inflammation and extracellular matrix degradation, have been correlated with VCAP characteristics and adverse prognosis. However, coronary computed tomography (CT) angiography has been the most extensively investigated technique, significantly correlating with invasive method-derived VCAP features. The estimation of perivascular fat attenuation as well as radiomic-based approaches represent additional concepts that may add incremental information. Cardiac magnetic resonance imaging (MRI) has also been evaluated in clinical studies, with promising results through the various image sequences that have been tested. As far as nuclear cardiology is concerned, the implementation of positron emission tomography in the VCAP assessment currently faces several limitations with the myocardial uptake of the radiotracer in cases of fluorodeoxyglucose use, as well as with motion correction. Moreover, the search for the ideal radiotracer and the most adequate combination (CT or MRI) is still ongoing. With a look to the future, the possible combination of imaging and circulating inflammatory and extracellular matrix degradation biomarkers in diagnostic and prognostic algorithms may represent the essential next step for the assessment of high-risk individuals.

Interferon regulatory factor-5-dependent CD11c+ macrophages contribute to the formation of rupture-prone atherosclerotic plaques

AIMS

Inflammation is a key factor in atherosclerosis. The transcription factor interferon regulatory factor-5 (IRF5) drives macrophages towards a pro-inflammatory state. We investigated the role of IRF5 in human atherosclerosis and plaque stability.

METHODS AND RESULTS

Bulk RNA sequencing from the Carotid Plaque Imaging Project biobank were used to mine associations between major macrophage associated genes and transcription factors and human symptomatic carotid disease. Immunohistochemistry, proximity extension assays, and Helios cytometry by time of flight (CyTOF) were used for validation. The effect of IRF5 deficiency on carotid plaque phenotype and rupture in ApoE-/- mice was studied in an inducible model of plaque rupture. Interferon regulatory factor-5 and ITGAX/CD11c were identified as the macrophage associated genes with the strongest associations with symptomatic carotid disease. Expression of IRF5 and ITGAX/CD11c correlated with the vulnerability index, pro-inflammatory plaque cytokine levels, necrotic core area, and with each other. Macrophages were the predominant CD11c-expressing immune cells in the plaque by CyTOF and immunohistochemistry. Interferon regulatory factor-5 immunopositive areas were predominantly found within CD11c+ areas with a predilection for the shoulder region, the area of the human plaque most prone to rupture. Accordingly, an inducible plaque rupture model of ApoE-/-Irf5-/- mice had significantly lower frequencies of carotid plaque ruptures, smaller necrotic cores, and less CD11c+ macrophages than their IRF5-competent counterparts.

CONCLUSION

Using complementary evidence from data from human carotid endarterectomies and a murine model of inducible rupture of carotid artery plaque in IRF5-deficient mice, we demonstrate a mechanistic link between the pro-inflammatory transcription factor IRF5, macrophage phenotype, plaque inflammation, and its vulnerability to rupture.

Plaque Vulnerability Index Predicts Cardiovascular Events: A Histological Study of an Endarterectomy Cohort

Background

The balance between stabilizing and destabilizing atherosclerotic plaque components is used in experimental studies and in imaging studies to identify rupture prone plaques. However, we lack the evidence that this balance predicts future cardiovascular events. Here we explore whether a calculated histological ratio, referred to as vulnerability index (VI), can predict patients at higher risk to suffer from future cardiovascular events.

Methods and Results

Carotid plaques and clinical information from 194 patients were studied. Tissue sections were used for histological analysis to calculate the VI (CD68 [cluster of differentiation 68], alpha‐actin, Oil red O, Movat pentachrome, and glycophorin A). Postoperative cardiovascular events were identified through the Swedish National Inpatient Health Register (2005–2013). During the follow‐up (60 months) 45 postoperative cardiovascular events were registered. Patients with a plaque VI in the fourth quartile compared with the first to third quartiles had significantly higher risk to suffer from a future cardiovascular event (P=0.0002). The VI was an independent predictor and none of the 5 histological variables analyzed separately predicted events. In the 13 patients who underwent bilateral carotid endarterectomy, the VI of the right plaque correlated with the VI of the left plaque and vice versa (r=0.7, P=0.01).

Conclusions

Our findings demonstrate that subjects with a high plaque VI have an increased risk of future cardiovascular events, independently of symptoms and other known cardiovascular risk factors . This strongly supports that techniques which image such plaques can facilitate risk stratification for subjects in need of more intense treatment.

Nonlinear multiscale analysis of coronary atherosclerotic vulnerable plaque artery: fluid-structural modeling with micromechanics

A unique three-dimensional (3D) computational multiscale modeling approach is proposed to investigate the influence of presence of microcalcification particles on the stress field distribution in the thin cap layer of a coronary atherosclerotic vulnerable plaque system. A nested 3D modeling analysis framework spanning the multiscale nature of a coronary atherosclerotic vulnerable plaque is presented. At the microscale level, a micromechanical modeling approach, which is based on computational finite-element (FE) representative unit cell, is applied to obtain the homogenized nonlinear response of the calcified tissue. This equivalent response effectively allows the integration of extremely small microcalcification inclusions in a global biomechanical FE model. Next, at the macroscale level, a 3D patient-based fluid-structure interaction FE model, reconstructing a refined coronary artery geometry with calcified plaque lesion, is generated to study the mechanical behavior of such multi-component biomechanical system. It is shown that the proposed multiscale modeling approach can generate a higher resolution of stress and strain field distributions within the coronary atherosclerotic vulnerable plaque system and allow the assessment of the local concentration stress around the microcalcifications in plaque cap layers. A comparison of stress field distributions within cap layers with and without inclusion of microcalcifications is also presented.

Predicting plaque vulnerability change using intravascular ultrasound + optical coherence tomography image-based fluid-structure interaction models and machine learning methods with patient follow-up data: a feasibility study

BACKGROUND

Coronary plaque vulnerability prediction is difficult because plaque vulnerability is non-trivial to quantify, clinically available medical image modality is not enough to quantify thin cap thickness, prediction methods with high accuracies still need to be developed, and gold-standard data to validate vulnerability prediction are often not available. Patient follow-up intravascular ultrasound (IVUS), optical coherence tomography (OCT) and angiography data were acquired to construct 3D fluid-structure interaction (FSI) coronary models and four machine-learning methods were compared to identify optimal method to predict future plaque vulnerability.

METHODS

Baseline and 10-month follow-up in vivo IVUS and OCT coronary plaque data were acquired from two arteries of one patient using IRB approved protocols with informed consent obtained. IVUS and OCT-based FSI models were constructed to obtain plaque wall stress/strain and wall shear stress. Forty-five slices were selected as machine learning sample database for vulnerability prediction study. Thirteen key morphological factors from IVUS and OCT images and biomechanical factors from FSI model were extracted from 45 slices at baseline for analysis. Lipid percentage index (LPI), cap thickness index (CTI) and morphological plaque vulnerability index (MPVI) were quantified to measure plaque vulnerability. Four machine learning methods (least square support vector machine, discriminant analysis, random forest and ensemble learning) were employed to predict the changes of three indices using all combinations of 13 factors. A standard fivefold cross-validation procedure was used to evaluate prediction results.

RESULTS

For LPI change prediction using support vector machine, wall thickness was the optimal single-factor predictor with area under curve (AUC) 0.883 and the AUC of optimal combinational-factor predictor achieved 0.963. For CTI change prediction using discriminant analysis, minimum cap thickness was the optimal single-factor predictor with AUC 0.818 while optimal combinational-factor predictor achieved an AUC 0.836. Using random forest for predicting MPVI change, minimum cap thickness was the optimal single-factor predictor with AUC 0.785 and the AUC of optimal combinational-factor predictor achieved 0.847.

CONCLUSION

This feasibility study demonstrated that machine learning methods could be used to accurately predict plaque vulnerability change based on morphological and biomechanical factors from multi-modality image-based FSI models. Large-scale studies are needed to verify our findings.

The Relation between Atherosclerosis Plaque Composition and Plaque Rupture

Background:

Intima, media, and adventitia are three layers of arteries. They have different structures and different mechanical properties. Damage to intima layer of arteries leads to an inflammatory response, which is usually the reason for atherosclerosis plaque formation. Atherosclerosis plaques mainly consist of smooth muscle cells and calcium. However, plaque geometry and mechanical properties change during time. Blood flow is the source of biomechanical stress to the plaques. Maximum stress that atherosclerosis plaque can burden before its rupture depends on fibrous cap thickness, lipid core, calcification, and artery stenosis. When atherosclerotic plaque ruptures, the blood would be in contact with coagulation factors. That is why plaque rupture is one of the main causes of fatality.

Method:

In this article, the coronary artery was modeled by ANSYS. First, fibrous cap thickness was increased from 40 μm to 250 μm by keeping other parameters constant. Then, the lipid pool percentage was incremented from 10% to 90% by keeping other parameters unchanged. Furthermore, for investigating the influence of calcium in plaque vulnerability, calcium was modeled in both agglomerated and microcalcium form.

Results:

It is proved that atherosclerosis plaque stress decreases exponentially as cap thickness increases. Larger lipid pool leads to more vulnerable plaques. In addition, the analysis showed maximum plaque stress usually increases in calcified plaque as compared with noncalcified plaque.

Conclusion:

The plaque stress is dependent on whether calcium is agglomerated near the lumen or far from it. However, in both cases, the deposition of more calcium in calcified plaque reduces maximum plaque stress.

Coronary arteries hemodynamics: effect of arterial geometry on hemodynamic parameters causing atherosclerosis

Coronary arteries have high curvatures, and hence, flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing realistic curvatures and stenosis and (ii) patient-specific coronary arteries with stenoses. Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Secondly, we have similarly determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), which were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis; our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves. This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient's coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is in our showing how in actual coronary arteries (based on their CT imaging) curvilinearity and narrowing complications affect the computed WSS and WPG, associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery. We want to go one step further and provide clinical application of our research work. For that, we are offering to cardiologists worldwide to carry out hemodynamic analysis of the medically imaged coronary arteries of their patients and compute the values of the hemodynamic parameters of WSS and WPG, so as to provide them an assessment of the risk of atherosclerosis for their patients. Graphical abstract Theme and aims: Coronary arteries have high curvatures, and hence flow through them causes disturbed flow patterns, resulting in stenosis and atherosclerosis. This in turn decreases the myocardial flow perfusion, causing myocardial ischemia and infarction. Therefore, in order to understand the mechanisms of these phenomena caused by high curvatures and branching of coronary arteries, we have conducted elaborate hemodynamic analysis for both (i) idealized coronary arteries with geometrical parameters representing curvatures and stenosis, and (ii) patient-specific coronary arteries with stenoses. Methods and results: Firstly, in idealized coronary arteries with approximated realistic arterial geometry representative of their curvedness and stenosis, we have computed the hemodynamic parameters of pressure drop, wall shear stress (WSS) and wall pressure gradient (WPG), and their association with the geometrical parameters of curvedness and stenosis. Then, we have determined the wall shear stress and wall pressure gradient distributions in four patient-specific curved stenotic right coronary arteries (RCAs), that were reconstructed from medical images of patients diagnosed with atherosclerosis and stenosis, as illustrated in Figure 1, in which the locations of the stenoses are highlighted by arrows. Figure 1: Three-dimensional CT visualization of arteries in patients with suspected coronary disease. The arteries can be seen as a combination of various curved segments with stenoses at unspecific locations highlighted by arrows. Our results show high WSS and WPG regions at the stenoses and inner wall of the arterial curves, as depicted in Figure 2. Therein, the encapsulations show (i) high WSS, and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Figure 2: WSS and WPG surface plot of realistic arteries (a), (b), (c) and (d), wherein the small squared parts are enlarged to show the detailed localized contour plots at the stenotic regions. Therein, the circular encapsulations show (i) high WSS and (ii) high WPG regions at the stenosis and inner wall of the arterial curves. Conclusion and novelty: This paper provides useful insights into the causative mechanisms of the high incidence of atherosclerosis in coronary arteries. It also provides guidelines for how simulation of blood flow in patient coronary arteries and determination of the hemodynamic parameters of WSS and WPG can provide a medical assessment of the risk of development of atherosclerosis and plaque formation, leading to myocardial ischemia and infarction. The novelty of our paper is our showing how in actual coronary arteries (based on their CT imaging), curvilinearity and narrowing complications affect the computed WSS and WPG associated with risk of atherosclerosis. This is very important for cardiologists to be able to properly take care of their patients and provide remedial measures before coronary complications lead to myocardial infarctions and necessitate stenting or coronary bypass surgery.

Detecting vulnerable plaque with vulnerability index based on convolutional neural networks

Plaque rupture and subsequent thrombosis are major processes of acute cardiovascular events. The Vulnerability Index is a very important indicator of whether a plaque is ruptured, and these easily ruptured or fragile plaques can be detected early. The higher the general vulnerability index, the higher the instability of the plaque. Therefore, determining a clear vulnerability index classification point can effectively reduce unnecessary interventional therapy. However, the current critical value of the vulnerability index has not been well defined. In this study, we proposed a neural network-based method to determine the critical point of vulnerability index that distinguishes vulnerable plaques from stable ones. Firstly, based on MatConvNet, the intravascular ultrasound images under different vulnerability index labels are classified. Different vulnerability indexes can obtain different accuracy rates for the demarcation points. The corresponding data points are fitted to find the existing relationship to judge the highest classification. In this way, the vulnerability index corresponding to the highest classification accuracy rate is judged. Then the article is based on the same experiment of different components of the aortic artery in the artificial neural network, and finally the vulnerability index corresponding to the highest classification accuracy can be obtained. The results show that the best vulnerability index point is 1.716 when the experiment is based on the intravascular ultrasound image, and the best vulnerability index point is 1.607 when the experiment is based on the aortic artery component data. Moreover, the vulnerability index and classification accuracy rate has a periodic relationship within a certain range, and finally the highest AUC is 0.7143 based on the obtained vulnerability index point on the verification set. In this paper, the convolution neural network is used to find the best vulnerability index classification points. The experimental results show that this method has the guiding significance for the classification and diagnosis of vulnerable plaques, further reduce interventional treatment of cardiovascular disease.

A patient-specific numerical modeling of the spontaneous coronary artery dissection in relation to atherosclerosis

The spontaneous coronary artery dissection (SCAD) is a clinical complication of angioplasty leading to an initiation of a tear/crack in the intima layer of the artery. The crack can propagate to the interface of the intima-media layer following by intramural hematoma. The relation between the SCAD and atherosclerosis is a controversial issue, as some studies stated no correlation between them while others showed that a crack can initiate in the intima but cannot propagate into the atrophied media layer. To investigate the relation between the intraluminal crack propagation in the atherosclerotic artery and SCAD, this study numerically investigated the initiation and propagation of a crack in the intraluminal and radial locations of the healthy and atherosclerotic human coronary arterial walls. The energy release rate, namely J-integral, is computed as a numerical derivative of the strain energy with respect to a crack extension using a user-defined virtual crack method (VCE) of extended finite element method (XFEM). Experimental measurements were carried out to calculate the elasto-plastic mechanical properties of the healthy and atherosclerotic human coronary arteries. The experimental data were then assigned to our own established patient-specific FE model of the coronary artery. Cracks were sketched in the intraluminal and radial locations of the arterial wall and allowed to propagate to the virtual interface of the intima-media to form a false lumen. The results revealed a higher stress at the crack tip of the healthy arterial wall compared to the atherosclerotic one. Lower crack tip opening displacement (CTOD) and crack tip opening angle (CTOA) were observed in the intraluminal crack of the atherosclerotic artery. J-integral of the atherosclerotic arterial wall was also found to be higher than the healthy one at the intraluminal crack. The results revealed that although a crack can initiate in the intraluminal of an atherosclerotic artery, it cannot propagate into the media layer due to a relatively higher rate of the strain energy release in the atherosclerotic arterial wall compared to the healthy one.

A magnetic resonance imaging-compatible small animal model under extracorporeal circulation

The impact of different extracorporeal circulation (ECC) scenarios on arterial blood flow profiles has not yet been revealed. To allow for exact measurements, magnetic resonance imaging (MRI) during ECC is required. Therefore, the present study addressed the feasibility of a high-resolution MRI-compatible animal model of ECC. For usage in New Zealand White rabbits, we developed an ECC device, the tubes of which were long enough to eliminate impacts of the magnetic field on the blood pump and heart-lung control machine. The miniaturized ECC system via thoracic access comprised an infant oxygenator, a pulsatile centrifugal pump, 1/8″ tubes, a 10-Fr aortic cannula and a 12-Fr venous cannula for vacuum-assisted drainage. This miniaturized ECC system has very low priming volume (230-255 ml) to reduce the system-inherent haemodilution to 50%. Consequently, haemoglobin rates remained high enough to guarantee adequate oxygenation (arterial pressure of oxygen >200 mmHg). Optimized venous drainage by an additionally inserted pulmonary artery vent catheter resulted in sufficient blood flow (31.6-65.8 ml/min/kg) that was maintained for 60 min with pulsatility. The current study demonstrates the feasibility of MRI-compatible ECC in rabbits, and this model allows for real-time blood flow profile measurements during different ECC scenarios in future projects.

Combining morphological and biomechanical factors for optimal carotid plaque progression prediction: An MRI-based follow-up study using 3D thin-layer models

Plaque progression prediction is of fundamental significance to cardiovascular research and disease diagnosis, prevention, and treatment. Magnetic resonance image (MRI) data of carotid atherosclerotic plaques were acquired from 20 patients with consent obtained. 3D thin-layer models were constructed to calculate plaque stress and strain. Data for ten morphological and biomechanical risk factors were extracted for analysis. Wall thickness increase (WTI), plaque burden increase (PBI) and plaque area increase (PAI) were chosen as three measures for plaque progression. Generalized linear mixed models (GLMM) with 5-fold cross-validation strategy were used to calculate prediction accuracy and identify optimal predictor. The optimal predictor for PBI was the combination of lumen area (LA), plaque area (PA), lipid percent (LP), wall thickness (WT), maximum plaque wall stress (MPWS) and maximum plaque wall strain (MPWSn) with prediction accuracy = 1.4146 (area under the receiver operating characteristic curve (AUC) value is 0.7158), while PA, plaque burden (PB), WT, LP, minimum cap thickness, MPWS and MPWSn was the best for WTI (accuracy = 1.3140, AUC = 0.6552), and a combination of PA, PB, WT, MPWS, MPWSn and average plaque wall strain (APWSn) was the best for PAI with prediction accuracy = 1.3025 (AUC = 0.6657). The combinational predictors improved prediction accuracy by 9.95%, 4.01% and 1.96% over the best single predictors for PAI, PBI and WTI (AUC values improved by 9.78%, 9.45%, and 2.14%), respectively. This suggests that combining both morphological and biomechanical risk factors could lead to better patient screening strategies.

Microcalcifications, Their Genesis, Growth, and Biomechanical Stability in Fibrous Cap Rupture

For many decades, cardiovascular calcification has been considered as a passive process, accompanying atheroma progression, correlated with plaque burden, and apparently without a major role on plaque vulnerability. Clinical and pathological analyses have previously focused on the total amount of calcification (calcified area in a whole atheroma cross section) and whether more calcification means higher risk of plaque rupture or not. However, this paradigm has been changing in the last decade or so. Recent research has focused on the presence of microcalcifications (μCalcs) in the atheroma and more importantly on whether clusters of μCalcs are located in the cap of the atheroma. While the vast majority of μCalcs are found in the lipid pool or necrotic core, they are inconsequential to vulnerable plaque. Nevertheless, it has been shown that μCalcs located within the fibrous cap could be numerous and that they behave as an intensifier of the background circumferential stress in the cap. It is now known that such intensifying effect depends on the size and shape of the μCalc as well as the proximity between two or more μCalcs. If μCalcs are located in caps with very low background stress, the increase in stress concentration may not be sufficient to reach the rupture threshold. However, the presence of μCalc(s) in the cap with a background stress of about one fifth to one half the rupture threshold (a stable plaque) will produce a significant increase in local stress, which may exceed the cap rupture threshold and thus transform a non-vulnerable plaque into a vulnerable one. Also, the classic view that treats cardiovascular calcification as a passive process has been challenged, and emerging data suggest that cardiovascular calcification may encompass both passive and active processes. The passive calcification process comprises biochemical factors, specifically circulating nucleating complexes, which would lead to calcification of the atheroma. The active mechanism of atherosclerotic calcification is a cell-mediated process via cell death of macrophages and smooth muscle cells (SMCs) and/or the release of matrix vesicles by SMCs.

Model-based cap thickness and peak cap stress prediction for carotid MRI

A rupture-prone carotid plaque can potentially be identified by calculating the peak cap stress (PCS). For these calculations, plaque geometry from MRI is often used. Unfortunately, MRI is hampered by a low resolution, leading to an overestimation of cap thickness and an underestimation of PCS. We developed a model to reconstruct the cap based on plaque geometry to better predict cap thickness and PCS. We used histological stained plaques from 34 patients. These plaques were segmented and served as the ground truth. Sections of these plaques contained 93 necrotic cores with a cap thickness <0.62mm which were used to generate a geometry-based model. The histological data was used to simulate in vivo MRI images, which were manually delineated by three experienced MRI readers. Caps below the MRI resolution (n=31) were (digitally removed and) reconstructed according to the geometry-based model. Cap thickness and PCS were determined for the ground truth, readers, and reconstructed geometries. Cap thickness was 0.07mm for the ground truth, 0.23mm for the readers, and 0.12mm for the reconstructed geometries. The model predicts cap thickness significantly better than the readers. PCS was 464kPa for the ground truth, 262kPa for the readers and 384kPa for the reconstructed geometries. The model did not predict the PCS significantly better than the readers. The geometry-based model provided a significant improvement for cap thickness estimation and can potentially help in rupture-risk prediction, solely based on cap thickness. Estimation of PCS estimation did not improve, probably due to the complex shape of the plaques.

NUMERICAL SIMULATION OF TRANSIENT BLOOD FLOW THROUGH THE LEFT CORONARY ARTERY WITH VARYING DEGREES OF BIFURCATION ANGLES

Atherosclerosis is a cardiovascular condition that can occur in any part of the vascular system. Especially, it can exist in bifurcated arteries such as the left and right coronary arteries, abdominal aortic bifurcation or carotid artery bifurcation. In our study, we examine the left coronary artery as an exemplification using wall shear stress (WSS) and wall pressure gradient (WPG). Then, we attempt to find the relationship between bifurcated arterial geometry and hemodynamics. Computational fluid dynamics (CFD) is a common technique applied to characterize blood flow accurately and assist us to gain an insight of atherosclerosis. In this paper, we used CFD as the computational hemodynamics analysis technique to examine flow through the left coronary artery that has variable angular bifurcation. Our results demonstrated that the region of low WSS area and magnitudes of maximum WPG increases with the angles of bifurcation. Such hemodynamic condition resulting from the large bifurcation angles has an effect on atherogenesis and is worthy of investigation for better understanding of atherosclerosis.

Cap inflammation leads to higher plaque cap strain and lower cap stress: An MRI-PET/CT-based FSI modeling approach

Plaque rupture may be triggered by extreme stress/strain conditions. Inflammation is also implicated and can be imaged using novel imaging techniques. The impact of cap inflammation on plaque stress/strain and flow shear stress were investigated. A patient-specific MRI-PET/CT-based modeling approach was used to develop 3D fluid-structure interaction models and investigate the impact of inflammation on plaque stress/strain conditions for better plaque assessment. 18FDG-PET/CT and MRI data were acquired from 4 male patients (average age: 66) to assess plaque characteristics and inflammation. Material stiffness for the fibrous cap was adjusted lower to reflect cap weakening causing by inflammation. Setting stiffness ratio (SR) to be 1.0 (fibrous tissue) for baseline, results for SR=0.5, 0.25, and 0.1 were obtained. Thin cap and hypertension were also considered. Combining results from the 4 patients, mean cap stress from 729 cap nodes was lowered by 25.2% as SR went from 1.0 to 0.1. Mean cap strain value for SR=0.1 was 0.313, 114% higher than that from SR=1.0 model. The thin cap SR=0.1 model had 40% mean cap stress decrease and 81% cap strain increase compared with SR=1.0 model. The hypertension SR=0.1 model had 19.5% cap stress decrease and 98.6% cap strain increase compared with SR=1.0 model. Differences of flow shear stress with 4 different SR values were limited (<10%). Cap inflammation may lead to large cap strain conditions when combined with thin cap and hypertension. Inflammation also led to lower cap stress. This shows the influence of inflammation on stress/strain calculations which are closely related to plaque assessment.

EFFECT OF PROGRESSIVE DEGREES OF SIDE-BRANCH BIFURCATION ANGLES ON FLOW IN ARTERY USING NUMERICAL SIMULATION

Malfunction of the cardiovascular system is a serious disease affecting human life around the world that is caused by several factors. One of the major factors is atherosclerosis that is a disease of the artery. Atherosclerosis is a serious vascular condition, which always occurs in branch vessels such as the abdominal aortic bifurcation and the carotid artery bifurcation. Wall shear stress (WSS) and wall pressure gradient (WPG) pertaining to these vessels will dramatically change when the geometry of these vessels is varied. Computational hemodynamics, as an invasive method, can be employed to understand the blood flow characteristics. In this paper, blood flow through arteries with variable side branches is presented using a computational framework. Numerical models pertaining to the different side-branch bifurcation angles are simulated to verify this. The side-branch bifurcation angle correlates positively to the maximum WSS experienced by the artery and this has an effect on atherogenesis. The low WPG regions are found to decrease with increased values of the angles, while the high WPG regions concentrate in the same region with larger values. Such hemodynamics information can be used to understand the effect of arterial geometrical variation on hemodynamics and the causes of atherosclerosis.

Morphological and Stress Vulnerability Indices for Human Coronary Plaques and Their Correlations with Cap Thickness and Lipid Percent: An IVUS-Based Fluid-Structure Interaction Multi-patient Study

Plaque vulnerability, defined as the likelihood that a plaque would rupture, is difficult to quantify due to lack of in vivo plaque rupture data. Morphological and stress-based plaque vulnerability indices were introduced as alternatives to obtain quantitative vulnerability assessment. Correlations between these indices and key plaque features were investigated. In vivo intravascular ultrasound (IVUS) data were acquired from 14 patients and IVUS-based 3D fluid-structure interaction (FSI) coronary plaque models with cyclic bending were constructed to obtain plaque wall stress/strain and flow shear stress for analysis. For the 617 slices from the 14 patients, lipid percentage, min cap thickness, critical plaque wall stress (CPWS), strain (CPWSn) and flow shear stress (CFSS) were recorded, and cap index, lipid index and morphological index were assigned to each slice using methods consistent with American Heart Association (AHA) plaque classification schemes. A stress index was introduced based on CPWS. Linear Mixed-Effects (LME) models were used to analyze the correlations between the mechanical and morphological indices and key morphological factors associated with plaque rupture. Our results indicated that for all 617 slices, CPWS correlated with min cap thickness, cap index, morphological index with r = -0.6414, 0.7852, and 0.7411 respectively (p<0.0001). The correlation between CPWS and lipid percentage, lipid index were weaker (r = 0.2445, r = 0.2338, p<0.0001). Stress index correlated with cap index, lipid index, morphological index positively with r = 0.8185, 0.3067, and 0.7715, respectively, all with p<0.0001. For all 617 slices, the stress index has 66.77% agreement with morphological index. Morphological and stress indices may serve as quantitative plaque vulnerability assessment supported by their strong correlations with morphological features associated with plaque rupture. Differences between the two indices may lead to better plaque assessment schemes when both indices were jointly used with further validations from clinical studies.

Cardiovascular diseases are closely related to atherosclerotic plaque progression and rupture. Early detection of vulnerable plaques and prediction of potential plaque rupture and related clinical events are of vital importance. Plaque vulnerability, defined as the likelihood that a plaque would rupture, is difficult to measure due to lack of in vivo plaque rupture data. Morphological and stress-based plaque vulnerability indices were introduced in this paper as alternatives to obtain quantitative vulnerability assessment with potential improvement of patient screening tools. In vivo intravascular ultrasound data were acquired from patients and computational coronary plaque models were constructed to obtain data for analysis and index assignments. For the 617 slices from the 14 patients, morphological and stress indices were assigned to each slice using methods consistent with American Heart Association plaque classification schemes. Correlation analyses were performed for all the morphological and mechanical factors considered. The stress index has 66.77% agreement with morphological index. Morphological and stress indices may serve as quantitative plaque vulnerability assessment supported by their strong correlations with morphological features associated with plaque rupture. Differences between the two indices may lead to better plaque assessment schemes when the complementary indices were jointly used with further validations from clinical studies.

Visualizing the stress distribution within vascular tissues using intravascular ultrasound elastography: a preliminary investigation

A methodology for computing the stress distribution of vascular tissue using finite element-based, intravascular ultrasound (IVUS) reconstruction elastography is described. This information could help cardiologists detect life-threatening atherosclerotic plaques and predict their propensity to rupture. The calculation of vessel stresses requires the measurement of strain from the ultrasound images, a calibrating pressure measurement and additional model assumptions. In this work, we conducted simulation studies to investigate the effect of varying the model assumptions, specifically Poisson's ratio and the outer boundary conditions, on the resulting stress fields. In both simulation and phantom studies, we created vessel geometries with two fibrous cap thicknesses to determine if we could detect a difference in peak stress (spatially) between the two. The results revealed that (i) Poisson's ratios had negligible impact on the accuracy of stress elastograms, (ii) the outer boundary condition assumption had the greatest effect on the resulting modulus and stress distributions and (iii) in simulation and in phantom experiments, our stress imaging technique was able to detect an increased peak stress for the vessel geometry with the smaller cap thickness. This work is a first step toward understanding and creating a robust stress measurement technique for evaluating atherosclerotic plaques using IVUS elastography.

A numerical study of the effect of varied blood pressure on the stability of carotid atherosclerotic plaque

Background

Blood pressure (BP) is associated with early atherosclerosis and plaque rupture because the BP variability can significantly affect the blood flow velocity and shear stress over the plaque. However, the mechanical response of BP variability to the plaque remains unclear. Therefore, we investigated the correlation between different maximum systolic blood pressure (SBP) and the stress distribution on plaque, as well as the stress over the plaque and blood velocity around the plaque using different BP variations, which are the BP variability in different phases during one cardiac cycle and beat-to-beat BP variability.

Method

We established a two-dimensional artery model with stenosis at the degree of 62.5%. Eight combinations of pulsatile pressure gradients between the inflow and outflow were implemented at the model. Three levels of fibrous cap thickness were taken into consideration to investigate the additional effect on the BP variability. Wall shear stress and stress/strain distribution over the plaque were derived as well as the oscillation shear index (OSI) to analyze the impact of the changing rate of BP.

Result

The stresses at diastole were 2.5% ± 1.8% lower than that at systole under the same pressure drop during one cycle. It was also found that elevated SBP might cause the immediate increment of stress in the present cycle (292% ± 72.3%), but slight reduction in the successive cycle (0.48% ± 0.4%).

Conclusion

The stress/strain distribution over the plaque is sensitive to the BP variability during one cardiac cycle, and the beat-to-beat BP variability could cause considerable impact on the progression of atherosclerosis in long-term.

Transient blood flow in elastic coronary arteries with varying degrees of stenosis and dilatations: CFD modelling and parametric study

In this paper, we have analysed pulsatile flow through partially occluded elastic arteries, to determine the haemodynamic parameters of wall shear stress (WSS), wall pressure gradient and pressure drops (ΔP), contributing to enhanced flow resistance and myocardial ischaemic regions which impair cardiac contractility and cause increased work load on the heart. In summary, it can be observed that stenoses in an artery significantly influence the haemodynamic parameters of wall shear stress and pressure drop in contrast to dilatations case. This deduces that stenosis plays a more critical role in plaque growth and vulnerability in contrast to dilatation, and should be the key element in cardiovascular pathology and diagnosis. Through quantitative analysis of WSS and ΔP, we have provided a clearer insight into the haemodynamics of atherosclerotic arteries. Determination of these parameters can be helpful to cardiologists, because it is directly implicated in the genesis and development of atherosclerosis.

Human coronary plaque wall thickness correlated positively with flow shear stress and negatively with plaque wall stress: an IVUS-based fluid-structure interaction multi-patient study

BACKGROUND

Atherosclerotic plaque progression and rupture are believed to be associated with mechanical stress conditions. In this paper, patient-specific in vivo intravascular ultrasound (IVUS) coronary plaque image data were used to construct computational models with fluid-structure interaction (FSI) and cyclic bending to investigate correlations between plaque wall thickness and both flow shear stress and plaque wall stress conditions.

METHODS

IVUS data were acquired from 10 patients after voluntary informed consent. The X-ray angiogram was obtained prior to the pullback of the IVUS catheter to determine the location of the coronary artery stenosis, vessel curvature and cardiac motion. Cyclic bending was specified in the model representing the effect by heart contraction. 3D anisotropic FSI models were constructed and solved to obtain flow shear stress (FSS) and plaque wall stress (PWS) values. FSS and PWS values were obtained for statistical analysis. Correlations with p < 0.05 were deemed significant.

RESULTS

Nine out of the 10 patients showed positive correlation between wall thickness and flow shear stress. The mean Pearson correlation r-value was 0.278 ± 0.181. Similarly, 9 out of the 10 patients showed negative correlation between wall thickness and plaque wall stress. The mean Pearson correlation r-value was -0.530 ± 0.210.

CONCLUSION

Our results showed that plaque vessel wall thickness correlated positively with FSS and negatively with PWS. The patient-specific IVUS-based modeling approach has the potential to be used to investigate and identify possible mechanisms governing plaque progression and rupture and assist in diagnosis and intervention procedures. This represents a new direction of research. Further investigations using more patient follow-up data are warranted.

The Influence of Inaccuracies in Carotid MRI Segmentation on Atherosclerotic Plaque Stress Computations

Biomechanical finite element analysis (FEA) based on in vivo carotid magnetic resonance imaging (MRI) can be used to assess carotid plaque vulnerability noninvasively by computing peak cap stress. However, the accuracy of MRI plaque segmentation and the influence this has on FEA has remained unreported due to the lack of a reliable submillimeter ground truth. In this study, we quantify this influence using novel numerical simulations of carotid MRI. Histological sections from carotid plaques from 12 patients were used to create 33 ground truth plaque models. These models were subjected to numerical computer simulations of a currently used clinically applied 3.0 T T1-weighted black-blood carotid MRI protocol (in-plane acquisition voxel size of 0.62 × 0.62 mm2) to generate simulated in vivo MR images from a known underlying ground truth. The simulated images were manually segmented by three MRI readers. FEA models based on the MRI segmentations were compared with the FEA models based on the ground truth. MRI-based FEA model peak cap stress was consistently underestimated, but still correlated (R) moderately with the ground truth stress: R = 0.71, R = 0.47, and R = 0.76 for the three MRI readers respectively (p < 0.01). Peak plaque stretch was underestimated as well. The peak cap stress in thick-cap, low stress plaques was substantially more accurately and precisely predicted (error of −12 ± 44 kPa) than the peak cap stress in plaques with caps thinner than the acquisition voxel size (error of −177 ± 168 kPa). For reliable MRI-based FEA to compute the peak cap stress of carotid plaques with thin caps, the current clinically used in-plane acquisition voxel size (∼0.6 mm) is inadequate. FEA plaque stress computations would be considerably more reliable if they would be used to identify thick-cap carotid plaques with low stresses instead.

Image-based modeling for better understanding and assessment of atherosclerotic plaque progression and vulnerability: data, modeling, validation, uncertainty and predictions

Medical imaging and image-based modeling have made considerable progress in recent years in identifying atherosclerotic plaque morphological and mechanical risk factors which may be used in developing improved patient screening strategies. However, a clear understanding is needed about what we have achieved and what is really needed to translate research to actual clinical practices and bring benefits to public health. Lack of in vivo data and clinical events to serve as gold standard to validate model predictions is a severe limitation. While this perspective paper provides a review of the key steps and findings of our group in image-based models for human carotid and coronary plaques and a limited review of related work by other groups, we also focus on grand challenges and uncertainties facing the researchers in the field to develop more accurate and predictive patient screening tools.

Effect of tissue properties, shape and orientation of microcalcifications on vulnerable cap stability using different hyperelastic constitutive models

Approximately half of all cardiovascular deaths associated with acute coronary syndrome occur when the thin fibrous cap tissue overlying the necrotic core in a coronary vessel is torn, ripped or fissured under the action of high blood pressure. From a biomechanics point of view, the rupture of an atheroma is due to increased mechanical stresses in the lesion, in which the ultimate stress (i.e. peak circumferential stress (PCS) at failure) of the tissue is exceeded. Several factors including the cap thickness, morphology, residual stresses and tissue composition of the atheroma have been shown to affect the PCS. Also important, we recently demonstrated that microcalcifications (μCalcs>5 µm are a common feature in human atheroma caps, which behave as local stress concentrators, increasing the local tissue stress by at least a factor of two surpassing the ultimate stress threshold for cap tissue rupture. In the present study, we used both idealized µCalcs with spherical shape and actual µCalcs from human coronary atherosclerotic caps, to determine their effect on increasing the circumferential stress in the fibroatheroma cap using different hyperelastic constitutive models. We have found that the stress concentration factor (SCF) produced by μCalcs in the fibroatheroma cap is affected by the material tissue properties, μCalcs spacing, aspect ratio and their alignment relative to the tensile axis of the cap.

Higher critical plaque wall stress in patients who died of coronary artery disease compared with those who died of other causes: a 3D FSI study based on ex vivo MRI of coronary plaques

Mechanical forces play an important role in the rupture of vulnerable plaques. This process is often associated with cardiovascular syndromes, such as heart attack and stroke. In this study, magnetic resonance imaging (MRI)-based models were used to investigate the association between plaque wall stress (PWS) and coronary artery disease (CAD). Ex vivo MRI data of coronary plaques from 12 patients were used to construct 12 three-dimensional (3D) fluid-structure interaction (FSI) computational models. Six of the patients had died from CAD and six had died from non-CAD causes. PWS was assessed using all nodal points on the lumen surface of each plaque. The maximum PWS from all possible vulnerable sites of each plaque was defined as the 3D critical plaque wall stress (CPWS). Mean 3D CPWS in the CAD group was 94.3% higher than that in the non-CAD group (265.6 vs. 136.7 kPa, P=0.0029). There was no statistically significant difference in global maximum plaque wall stress (GMPWS) between the two groups (P=0.347). There was also no statistically significant difference in plaque burden between the CAD group (84.4±5%) and the non-CAD group (82.0±8%, P=0.552). The results indicate that plaques from patients who died from CAD were associated with higher CPWS compared with those from patients who died from non-CAD causes. With further validation, analysis of CPWS may prove to be an important component in assessment of plaque vulnerability.

Changing views of the biomechanics of vulnerable plaque rupture: a review

This review examines changing perspectives on the biomechanics of vulnerable plaque rupture over the past 25 years from the first finite element analyses (FEA) showing that the presence of a lipid pool significantly increases the local tissue stress in the atheroma cap to the latest imaging and 3D FEA studies revealing numerous microcalcifications in the cap proper and a new paradigm for cap rupture. The first part of the review summarizes studies describing the role of the fibrous cap thickness, tissue properties, and lesion geometry as main determinants of the risk of rupture. Advantages and limitations of current imaging technologies for assessment of vulnerable plaques are also discussed. However, the basic paradoxes as to why ruptures frequently did not coincide with location of PCS and why caps >65 μm thickness could rupture at tissue stresses significantly below the 300 kPa critical threshold still remained unresolved. The second part of the review describes recent studies in the role of microcalcifications, their origin, shape, and clustering in explaining these unresolved issues including the actual mechanism of rupture due to the explosive growth of tiny voids (cavitation) in local regions of high stress concentration between closely spaced microinclusions oriented along their tensile axis.

Image-based modeling and precision medicine: patient-specific carotid and coronary plaque assessment and predictions

Atherosclerotic plaques may rupture without warning and cause acute cardiovascular events such as heart attack and stroke. Current clinical screening tools are insufficient to identify those patients with risks early and prevent the adverse events from happening. Medical imaging and image-based modeling have made considerable progress in recent years in identifying plaque morphological and mechanical risk factors which may be used in developing improved patient screening strategies. The key steps and factors in image-based models for human carotid and coronary plaques were illustrated, as well as grand challenges facing the researchers in the field to develop more accurate screening tools.

A mechanistic analysis of the role of microcalcifications in atherosclerotic plaque stability: potential implications for plaque rupture

The role of microcalcifications (μCalcs) in the biomechanics of vulnerable plaque rupture is examined. Our laboratory previously proposed (Ref. 44), using a very limited tissue sample, that μCalcs embedded in the fibrous cap proper could significantly increase cap instability. This study has been greatly expanded. Ninety-two human coronary arteries containing 62 fibroatheroma were examined using high-resolution microcomputed tomography at 6.7-μm resolution and undecalcified histology with special emphasis on calcified particles <50 μm in diameter. Our results reveal the presence of thousands of μCalcs, the vast majority in lipid pools where they are not dangerous. However, 81 μCalcs were also observed in the fibrous caps of nine of the fibroatheroma. All 81 of these μCalcs were analyzed using three-dimensional finite-element analysis, and the results were used to develop important new clinical criteria for cap stability. These criteria include variation of the Young's modulus of the μCalc and surrounding tissue, μCalc size, and clustering. We found that local tissue stress could be increased fivefold when μCalcs were closely spaced, and the peak circumferential stress in the thinnest nonruptured cap (66 μm) if no μCalcs were present was only 107 kPa, far less than the proposed minimum rupture threshold of 300 kPa. These results and histology suggest that there are numerous μCalcs < 15 μm in the caps, not visible at 6.7-μm resolution, and that our failure to find any nonruptured caps between 30 and 66 μm is a strong indication that many of these caps contained μCalcs.

A numerical parametric study of the mechanical action of pulsatile blood flow onto axisymmetric stenosed arteries

In the present paper, a fluid-structure interaction model is developed, questioning how the mechanical action of the blood onto an atheromatous plaque is affected by the length and the severity of the stenosis. An axisymmetric model is considered. The fluid is assumed Newtonian. The plaque is modelled as a heterogeneous hyperelastic anisotropic solid composed of the arterial wall, the lipid core and the fibrous cap. Transient velocity and pressure conditions of actual pulsatile blood flow are prescribed. The simulation is achieved using the Arbitrary Lagrangian Eulerian scheme in the COMSOL commercial Finite Element package. The results reveal different types of behavior in function of the length (denoted L) and severity (denoted S) of the stenosis. Whereas large plaques (L>10 mm) are mostly deformed under the action of the blood pressure, it appears that shorter plaques (L<10 mm) are significantly affected by the shear stresses. The shear stresses tend to deform the plaque by pinching it. This effect is called: "the pinching effect". It has an essential influence on the mechanical response of the plaque. For two plaques having the same radius severity S=45%, the maximum stress in the fibrous cap is 50% larger for the short plaque (L=5 mm) than for a larger plaque (L=10 mm), and the maximum wall shear stress is increased by 100%. Provided that they are confirmed by experimental investigations, these results may offer some new perspectives for understanding the vulnerability of short plaques.

IVUS-based computational modeling and planar biaxial artery material properties for human coronary plaque vulnerability assessment

Image-based computational modeling has been introduced for vulnerable atherosclerotic plaques to identify critical mechanical conditions which may be used for better plaque assessment and rupture predictions. In vivo patient-specific coronary plaque models are lagging due to limitations on non-invasive image resolution, flow data, and vessel material properties. A framework is proposed to combine intravascular ultrasound (IVUS) imaging, biaxial mechanical testing and computational modeling with fluid-structure interactions and anisotropic material properties to acquire better and more complete plaque data and make more accurate plaque vulnerability assessment and predictions. Impact of pre-shrink-stretch process, vessel curvature and high blood pressure on stress, strain, flow velocity and flow maximum principal shear stress was investigated.

Effect of calcification on the mechanical stability of plaque based on a three-dimensional carotid bifurcation model

Background

This study characterizes the distribution and components of plaque structure by presenting a three-dimensional blood-vessel modelling with the aim of determining mechanical properties due to the effect of lipid core and calcification within a plaque. Numerical simulation has been used to answer how cap thickness and calcium distribution in lipids influence the biomechanical stress on the plaque.

Method

Modelling atherosclerotic plaque based on structural analysis confirms the rationale for plaque mechanical examination and the feasibility of our simulation model. Meaningful validation of predictions from modelled atherosclerotic plaque model typically requires examination of bona fide atherosclerotic lesions. To analyze a more accurate plaque rupture, fluid-structure interaction is applied to three-dimensional blood-vessel carotid bifurcation modelling. A patient-specific pressure variation is applied onto the plaque to influence its vulnerability.

Results

Modelling of the human atherosclerotic artery with varying degrees of lipid core elasticity, fibrous cap thickness and calcification gap, which is defined as the distance between the fibrous cap and calcification agglomerate, form the basis of our rupture analysis. Finite element analysis shows that the calcification gap should be conservatively smaller than its threshold to maintain plaque stability. The results add new mechanistic insights and methodologically sound data to investigate plaque rupture mechanics.

Conclusion

Structural analysis using a three-dimensional calcified model represents a more realistic simulation of late-stage atherosclerotic plaque. We also demonstrate that increases of calcium content that is coupled with a decrease in lipid core volume can stabilize plaque structurally.

Microcalcifications increase coronary vulnerable plaque rupture potential: a patient-based micro-CT fluid-structure interaction study

Asymptomatic vulnerable plaques (VP) in coronary arteries accounts for significant level of morbidity. Their main risk is associated with their rupture which may prompt fatal heart attacks and strokes. The role of microcalcifications (micro-Ca), embedded in the VP fibrous cap, in the plaque rupture mechanics has been recently established. However, their diminutive size offers a major challenge for studying the VP rupture biomechanics on a patient specific basis. In this study, a highly detailed model was reconstructed from a post-mortem coronary specimen of a patient with observed VP, using high resolution micro-CT which captured the microcalcifications embedded in the fibrous cap. Fluid-structure interaction (FSI) simulations were conducted in the reconstructed model to examine the combined effects of micro-Ca, flow phase lag and plaque material properties on plaque burden and vulnerability. This dynamic fibrous cap stress mapping elucidates the contribution of micro-Ca and flow phase lag VP vulnerability independently. Micro-Ca embedded in the fibrous cap produced increased stresses predicted by previously published analytical model, and corroborated our previous studies. The 'micro-CT to FSI' methodology may offer better diagnostic tools for clinicians, while reducing morbidity and mortality rates for patients with vulnerable plaques and ameliorating the ensuing healthcare costs.

In vivo serial MRI-based models and statistical methods to quantify sensitivity and specificity of mechanical predictors for carotid plaque rupture: location and beyond

It has been hypothesized that mechanical risk factors may be used to predict future atherosclerotic plaque rupture. Truly predictive methods for plaque rupture and methods to identify the best predictor(s) from all the candidates are lacking in the literature. A novel combination of computational and statistical models based on serial magnetic resonance imaging (MRI) was introduced to quantify sensitivity and specificity of mechanical predictors to identify the best candidate for plaque rupture site prediction. Serial in vivo MRI data of carotid plaque from one patient was acquired with follow-up scan showing ulceration. 3D computational fluid-structure interaction (FSI) models using both baseline and follow-up data were constructed and plaque wall stress (PWS) and strain (PWSn) and flow maximum shear stress (FSS) were extracted from all 600 matched nodal points (100 points per matched slice, baseline matching follow-up) on the lumen surface for analysis. Each of the 600 points was marked "ulcer" or "nonulcer" using follow-up scan. Predictive statistical models for each of the seven combinations of PWS, PWSn, and FSS were trained using the follow-up data and applied to the baseline data to assess their sensitivity and specificity using the 600 data points for ulcer predictions. Sensitivity of prediction is defined as the proportion of the true positive outcomes that are predicted to be positive. Specificity of prediction is defined as the proportion of the true negative outcomes that are correctly predicted to be negative. Using probability 0.3 as a threshold to infer ulcer occurrence at the prediction stage, the combination of PWS and PWSn provided the best predictive accuracy with (sensitivity, specificity) = (0.97, 0.958). Sensitivity and specificity given by PWS, PWSn, and FSS individually were (0.788, 0.968), (0.515, 0.968), and (0.758, 0.928), respectively. The proposed computational-statistical process provides a novel method and a framework to assess the sensitivity and specificity of various risk indicators and offers the potential to identify the optimized predictor for plaque rupture using serial MRI with follow-up scan showing ulceration as the gold standard for method validation. While serial MRI data with actual rupture are hard to acquire, this single-case study suggests that combination of multiple predictors may provide potential improvement to existing plaque assessment schemes. With large-scale patient studies, this predictive modeling process may provide more solid ground for rupture predictor selection strategies and methods for image-based plaque vulnerability assessment.

Association of lipoprotein-associated phospholipase A₂ with characteristics of vulnerable coronary atherosclerotic plaques

PURPOSE

Lipoprotein-associated phospholipase A₂ (Lp-PLA₂) is an inflammatory enzyme expressed in atherosclerotic plaques. We investigated the association of circulating Lp-PLA₂ with characteristics of vulnerable coronary atherosclerotic plaques.

MATERIALS AND METHODS

We recruited 113 patients with either unstable angina (UA, n=59) and stable angina (SA, n=54) by coronary angiography. Thirty-six healthy subjects served as controls. Intravascular ultrasound (IVUS) was used to evaluate the characteristics of coronary atherosclerotic plaque, and serum Lp-PLA₂ concentration was measured as well.

RESULTS

Lp-PLA₂ concentration was significantly higher in both UA and SA patients [(396±36) μg/L and (321±39) μg/L, respectively] compared with the controls [(127 ± 49) μg/L, p<0.01], and higher in UA than SA group. IVUS findings showed that remodeling index (RI) (0.91 ± 0.15 vs. 0.85 ± 0.11, p=0.005) and eccentricity index (EI) (0.73 ± 0.16 vs. 0.65 ± 0.22, p=0.039) were larger in UA than in SA group, and fibrous caps were thicker in SA than UA group [(0.91 ± 0.23) mm vs. (0.63 ± 0.21) mm, p=0.032]. Moreover, Lp-PLA₂ correlated positively with EI (r=0.439, p<0.01) and RI (r=0.592, p<0.05) in UA group. There was an inverse relationship between Lp-PLA₂ and fibrous cap thickness in both UA (r=-0.587, p<0.001) and SA (r=-0.318, p<0.05) groups. The independent risk factors in UA group were Lp-PLA₂ (OR=1.055, 95% CI: 1.03-1.08, p=0.013), LDL-cholesterol (OR=0.032, 95% CI: 0.00-0.05, p=0.041) and fibrous cap thickness (OR=0.008, 95% CI: 0.00-0.45, p=0.019). Lp-PLA₂ was strongly associated with both EI and fibrous cap thickness in both groups.

CONCLUSION

Serum level of Lp-PLA₂ is associated with both eccentricity index and fibrous cap thickness in both UA and SA groups. Elevated levels of circulating Lp-PLA₂ might to be a strong risk factor and more serious for unstable angina than stable angina.

Study of carotid arterial plaque stress for symptomatic and asymptomatic patients

Stroke is one of the leading causes of death in the world, resulting mostly from the sudden ruptures of atherosclerosis carotid plaques. Until now, the exact plaque rupture mechanism has not been fully understood, and also the plaque rupture risk stratification. The advanced multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components to be visualized in-vivo and reconstructed by computational modeling. In the study, plaque stress analysis using fully coupled fluid structure interaction was applied to 20 patients (12 symptomatic and 8 asymptomatic) reconstructed from in-vivo MRI, followed by a detailed biomechanics analysis, and morphological feature study. The locally extreme stress conditions can be found in the fibrous cap region, 85% at the plaque shoulder based on the present study cases. Local maximum stress values predicted in the plaque region were found to be significantly higher in symptomatic patients than that in asymptomatic patients (200 ± 43 kPa vs. 127 ± 37 kPa, p=0.001). Plaque stress level, defined by excluding 5% highest stress nodes in the fibrous cap region based on the accumulative histogram of stress experienced on the computational nodes in the fibrous cap, was also significantly higher in symptomatic patients than that in asymptomatic patients (154 ± 32 kPa vs. 111 ± 23 kPa, p<0.05). Although there was no significant difference in lipid core size between the two patient groups, symptomatic group normally had a larger lipid core and a significantly thinner fibrous cap based on the reconstructed plaques using 3D interpolation from stacks of 2D contours. Plaques with a higher stenosis were more likely to have extreme stress conditions upstream of plaque throat. The combined analyses of plaque MR image and plaque stress will advance our understanding of plaque rupture, and provide a useful tool on assessing plaque rupture risk.

Effects of intima stiffness and plaque morphology on peak cap stress

BACKGROUND

Rupture of the cap of a vulnerable plaque present in a coronary vessel may cause myocardial infarction and death. Cap rupture occurs when the peak cap stress exceeds the cap strength. The mechanical stress within a cap depends on the plaque morphology and the material characteristics of the plaque components. A parametric study was conducted to assess the effect of intima stiffness and plaque morphology on peak cap stress.

METHODS

Models with idealized geometries based on histology images of human coronary arteries were generated by varying geometric plaque features. The constructed multi-layer models contained adventitia, media, intima, and necrotic core sections. For adventitia and media layers, anisotropic hyperelastic material models were used. For necrotic core and intima sections, isotropic hyperelastic material models were employed. Three different intima stiffness values were used to cover the wide range reported in literature. According to the intima stiffness, the models were classified as stiff, intermediate and soft intima models. Finite element method was used to compute peak cap stress.

RESULTS

The intima stiffness was an essential determinant of cap stresses. The computed peak cap stresses for the soft intima models were much lower than for stiff and intermediate intima models. Intima stiffness also affected the influence of morphological parameters on cap stresses. For the stiff and intermediate intima models, the cap thickness and necrotic core thickness were the most important determinants of cap stresses. The peak cap stress increased three-fold when the cap thickness was reduced from 0.25 mm to 0.05 mm for both stiff and intermediate intima models. Doubling the thickness of the necrotic core elevated the peak cap stress by 60% for the stiff intima models and by 90% for the intermediate intima models. Two-fold increase in the intima thickness behind the necrotic core reduced the peak cap stress by approximately 25% for both intima models. For the soft intima models, cap thickness was less critical and changed the peak cap stress by 55%. However, the necrotic core thickness was more influential and changed the peak cap stress by 100%. The necrotic core angle emerged as a critical determinant of cap stresses where a larger angle lowered the cap stresses. Contrary to the stiff and intermediate intima models, a thicker intima behind the necrotic core increased the peak cap stress by approximately 25% for the soft intima models. Adventitia thickness and local media regression had limited effects for all three intima models.

CONCLUSIONS

For the stiff and intermediate intima models, the cap thickness was the most important morphological risk factor. However for soft intima models, the necrotic core thickness and necrotic core angle had a bigger impact on the peak cap stress. We therefore need to enhance our knowledge of intima material properties if we want to derive critical morphological plaque features for risk evaluation.

Numerical modeling of stress in stenotic arteries with microcalcifications: a micromechanical approximation

Most finite element models of atherosclerotic arteries do not account for the heterogeneity of the plaque constituents at the microscale. Failure of plaque lesions has been shown to be a local event, linked to stress concentrations caused by cap thinning, inflammation, macroscopic heterogeneity, and recently, the presence of microcalcifications. There is growing evidence that microcalcifications exist in the fibrous cap of plaque lesions. However, their role is not yet fully understood. The goal of the present work is to investigate the effects of localized regions of microcalcifications on the stress field of atherosclerotic plaque caps in a section of carotid artery. This is achieved by performing finite element simulations of three-dimensional fluid-structure interaction models. The material response in the region of microcalcification is modeled using a combination of finite elements, homogenization theory, and a stress concentration function that approximates the average local stresses in the fibrous tissue and microcalcification phases. The results indicate that the circumferential stress in the fibrous tissue phase increases as the volume fraction of microcalcifications is increased, and that the stress exceeds a critical threshold when the fibrous cap thickness is decreased. Furthermore, the presence of the microcalcifications significantly influences the distribution of stress by shifting the maximum circumferential stress away from the cap shoulders, where failure is most common when the effective region of microcalcification is located at the center of the cap. This is a possible explanation of why 40% of plaque ruptures occur away from the shoulder region of the cap.

Patient-Specific Carotid Plaque Progression Simulation Using 3D Meshless Generalized Finite Difference Models with Fluid-Structure Interactions Based on Serial In Vivo MRI Data

Previously, we introduced a computational procedure based on three-dimensional meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Structure-only models were used in our previous report. In this paper, fluid-stricture interaction (FSI) was added to improve on prediction accuracy. One participating patient was scanned three times (T1, T2, and T3, at intervals of about 18 months) to obtain plaque progression data. Blood flow was assumed to laminar, Newtonian, viscous and incompressible. The Navier-Stokes equations with arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations. Plaque material was assumed to be uniform, homogeneous, isotropic, linear, and nearly incompressible. The linear elastic model was used. The 3D FSI plaque model was discretized and solved using a meshless generalized finite difference (GFD) method. Growth functions with a) morphology alone; b) morphology and plaque wall stress (PWS); morphology and flow shear stress (FSS), and d) morphology, PWS and FSS were introduced to predict future plaque growth based on previous time point data. Starting from the T2 plaque geometry, plaque progression was simulated by solving the FSI model and adjusting plaque geometry using plaque growth functions iteratively until T3 is reached. Numerically simulated plaque progression agreed very well with the target T3 plaque geometry with errors ranging from 8.62%, 7.22%, 5.77% and 4.39%, with the growth function including morphology, plaque wall stress and flow shear stress terms giving the best predictions. Adding flow shear stress term to the growth function improved the prediction error from 7.22% to 4.39%, a 40% improvement. We believe this is the first time 3D plaque progression FSI simulation based on multi-year patient-tracking data was reported. Serial MRI-based progression simulation adds time dimension to plaque vulnerability assessment and will improve prediction accuracy for potential plaque rupture risk.

Intraplaque hemorrhage is associated with higher structural stresses in human atherosclerotic plaques: an in vivo MRI-based 3D fluid-structure interaction study

Background

Studies using medical images have shown that intraplaque hemorrhage may accelerate plaque progression and may produce a stimulus for atherosclerosis development by increasing lipid core and plaque volume and creating new destabilizing factors. Image-based 3D computational models with fluid-structure interactions (FSI) will be used to perform plaque mechanical analysis and investigate possible associations between intraplaque hemorrhage and both plaque wall stress (PWS) and flow shear stress (FSS).

Methods

In vivo MRI data of carotid plaques from 5 patients with intraplaque hemorrhage confirmed by histology were acquired. 3D multi-component FSI models were constructed for each plaque to obtain mechanical stresses. Plaque Wall Stress (PWS) and Flow Shear Stress (FSS) were extracted from all nodal points on the lumen surface of each plaque for analysis.

Results

The mean PWS value from all hemorrhage nodes of the 5 plaques combined was higher than that from non-hemorrhage nodes (75.6 versus 68.1 kPa, P = 0.0003). The mean PWS values from hemorrhage nodes for each of the 5 plaques were all significantly higher (5 out of 5) than those from non-hemorrhage nodes (P < 0.05). The mean FSS value from all hemorrhage nodes of the 5 plaques combined was 30.4% higher than that from all non-hemorrhage nodes (15.0 versus 11.5 dyn/cm², P = 0.0002). However, the mean flow shear stress values from individual cases showed mixed results: only one out of five plaques showed mean FSS value from hemorrhage nodes was higher than that from non-hemorrhage nodes; three out of five plaques showed that their mean FSS values from hemorrhage nodes were lower than those from non-hemorrhage nodes; and one plaque showed that the difference had no statistical significance.

Conclusion

The results of this study suggested that intraplaque hemorrhage nodes were associated with higher plaque wall stresses. Compared to flow shear stress, plaque wall stress has a better correlation with plaque component feature (hemorrhage) linked to plaque progression and vulnerability. With further validation, plaque stress analysis may provide additional stress indicators for image-based vulnerability assessment.

Effect of posture change on the geometric features of the healthy carotid bifurcation

Segmented cross-sectional MRI images were used to construct 3-D virtual models of the carotid bifurcation in ten healthy volunteers. Geometric features, such as bifurcation angle, internal carotid artery (ICA) angle, planarity angle, asymmetry angle, tortuosity, curvature, bifurcation area ratio, ICA/common carotid artery (CCA), external carotid artery (ECA)/CCA, and ECA/ICA diameter ratios, were calculated for both carotids in two head postures: 1) the supine neutral position; and 2) the prone sleeping position with head rotation to the right ( ∼ 80°). The results obtained have shown that head rotation causes 1) significant variations in bifurcation angle [32% mean increase for the right carotid (RC) and 21% mean decrease for the left carotid (LC)] and internal carotid artery angle (97% mean increase for the RC, 43% mean decrease for the LC); 2) a slight increase in planarity and asymmetry angles for both RC and LC; 3) minor and variable curvature changes for the CCA and for the branches; 4) slight tortuosity changes for the braches but not for the CCA; and 5) unsubstantial alterations in area and diameter ratios (percentage changes %). The significant geometric changes observed in most subjects with head posture may also cause significant changes in bifurcation hemodynamics and warrant future investigation of the hemodynamic parameters related to the development of atherosclerotic disease such as low oscillating wall shear stress and particle residence times.

Hemodynamics of ulcerated plaques: before and after

The hemodynamics and fluid mechanical forces in blood vessels have long been implicated in the deposition and growth of atherosclerotic plaque. Detailed information about the hemodynamics in vessels affected by significant plaque deposits can also provide insight into the mechanisms and likelihood of plaque weakening and rupture. In the current study, the governing equations are solved in their finite volume formulation in several patient-specific stenotic geometries. Of specific interest are the flow patterns and forces near ulcerations in the plaque. The flow patterns and forces in vessels with ulcerated plaques are compared with those in stenotic vessels without evidence of ulceration and to the hemodynamics in the same vessels as they likely appeared prior to ulceration. Hemodynamics "before" and "after" hemorrhage may suggest fluid mechanical and morphological factors of diagnostic and predictive value. Recirculation zones, vortex shedding, and secondary flows are captured by both laminar and turbulent solutions. The forces on vessel walls are shown to correlate with unstable plaque deposits. Performing before and after studies of vessels in long-term radiology studies may illuminate mechanisms of hemorrhage and other vessel remodeling.