Background and Objective:Coronary plaque rupture is a precipitating event responsible for two thirds of myocardial infarctions. Currently, the risk of plaque rupture is computed based on demographic, clinical, and image-based adverse features. However, using these features the absolute event rate per single higher-risk lesion remains low. This work studies the power of a novel framework based on biomechanical markers accounting for material uncertainty to stratify vulnerable and non-vulnerable coronary plaques. Methods:Virtual histology intravascular ultrasounds from 55 patients, 29 affected by acute coronary syndrome and 26 affected by stable angina pectoris, were included in this study. Two-dimensional vessel cross-sections for finite element modeling (10 sections per plaque) incorporating plaque structure (medial tissue, loose matrix, lipid core and calcification) were reconstructed. A Montecarlo finite element analysis was performed on each section to account for material variability on three biomechanical markers: peak plaque structural stress at diastolic and systolic pressure, and peak plaque stress difference between systolic and diastolic pressures, together with the luminal pressure. Machine learning decision tree classifiers were trained on 75% of the dataset and tested on the remaining 25% with a combination of feature selection techniques. Performance against classification trees based on geometric markers (i.e., luminal, external elastic membrane and plaque areas) was also performed. Results:Our results indicate that the plaque structural stress outperforms the classification capacity of the combined geometric markers only (0.82 vs 0.51 area under curve) when accounting for uncertainty in material parameters. Furthermore, the results suggest that the combination of the peak plaque structural stress at diastolic and systolic pressures with the maximum plaque structural stress difference between systolic and diastolic pressures together with the systolic pressure and the diastolic to systolic pressure gradient is a robust classifier for coronary plaques when the intrinsic variability in material parameters is considered (area under curve equal to [0.91–0.93]). Conclusion:In summary, our results emphasize that peak plaque structural stress in combination with the patient’s luminal pressure is a potential classifier of plaque vulnerability as it independently considers stress in all directions and incorporates total geometric and compositional features of atherosclerotic plaques.