Background: Stress calculations in atherosclerotic coronary vulnerable plaques can aid in predicting coronary cap rupture. In vivo plaque geometry and composition of coronary arteries can merely be obtained via intravascular imaging. Only optical driven imaging techniques have sufficient resolution to visualize the fibrous cap, but due to limited penetration depth deeper components such as the backside of the necrotic core (NC) are generally not visible. The goal of this study was to investigate whether peak cap stresses can be approximated by reconstructing the backside of the NC. Methods: Manual segmentations of coronary histological cross-sections served as a geometrical ground truth and were obtained from seven patients resulting in 73 NCs. Next, the backside was removed and reconstructed according to an estimation of the relative necrotic core thickness (rNCt). The rNCt was estimated at three locations along the NC angle and based on either group averaged parameters or plaque specific parameters. Stress calculations were performed in both the ground truth geometry and the reconstructed geometries and compared. Results: Good geometrical agreement was found between the ground truth NC and the reconstructed NCs, based on group averaged rNCt estimation and plaque specific rNCt estimation, measuring the NC area difference (25.1 % IQR 14.0-41.3 % and 17.9 % IQR 9.81-32.7 %) and similarity index (0.85 IQR 0.77-0.90 and 0.88 IQR 0.79-0.91). The peak cap stresses obtained with both reconstruction methods showed a high correlation with respect to the ground truth, r² = 0.91 and r² = 0.95, respectively. For high stress plaques, the peak cap stress difference with respect to the ground truth significantly improved for the NC reconstruction based plaque specific features (6 %) compared to the reconstruction group averaged based (16 %). Conclusions: In conclusion, good geometry and stress agreement was observed between the ground truth NC geometry and the reconstructed geometries. Although group averaged rNCt estimation seemed to be sufficient for the NC reconstruction and stress calculations, including plaque specific data further improved stress predictions, especially for higher stresses.

Biomechanical models have the potential to predict plaque rupture. For reliable models, correct material properties of plaque components are a prerequisite. This study presents a new technique, where high resolution ultrasound displacement imaging and inverse finite element (FE) modeling is combined, to estimate material properties of plaque components. Iliac arteries with plaques were excised from 6 atherosclerotic pigs and subjected to an inflation test with pressures ranging from 10 to 120 mmHg. The arteries were imaged with high frequency 40 MHz ultrasound. Deformation maps of the plaques were reconstructed by cross correlation of the ultrasound radiofrequency data. Subsequently, the arteries were perfusion fixed for histology and structural components were identified. The histological data were registered to the ultrasound data to construct FE model of the plaques. Material properties of the arterial wall and the intima of the atherosclerotic plaques were estimated using a grid search method. The computed displacement fields showed good agreement with the measured displacement fields, implying that the FE models were able to capture local inhomogeneities within the plaque. On average, nonlinear stiffening of both the wall and the intima was observed, and the wall of the atheroslcerotic porcine iliac arteries was markedly stiffer than the intima (877 ± 459 vs. 100 ± 68 kPa at 100 mmHg). The large spread in the data further illustrates the wide variation of the material properties. We demonstrated the feasibility of a mixed experimental–numerical framework to determine the material properties of arterial wall and intima of atherosclerotic plaques from intact arteries, and concluded that, due to the observed variation, plaque specific properties are required for accurate stress simulations.