Quantifying Patient-Specific in vivo Coronary Plaque Material Properties for Accurate Stress/Strain Calculations: An IVUS-Based Multi-Patient Study

Comparison of multilayer and single-layer coronary plaque models on stress/strain calculations based on optical coherence tomography images

Mechanical stress and strain conditions are closely related to atherosclerotic plaque progression and rupture and have been under intensive investigations in recent years. It is well known that arteries have a three-layer structure: intima, media and adventitia. However, in vivo image-based multilayer plaque models are not available in the current literature due to lack of multilayer image segmentation data. A multilayer segmentation and repairing technique was introduced to segment coronary plaque optical coherence tomography (OCT) image to obtain its three-layer vessel structure. A total of 200 OCT slices from 20 patients (13 male; 7 female) were used to construct multilayer and single-layer 3D thin-slice models to calculate plaque stress and strain and compare model differences. Our results indicated that the average maximum plaque stress values of 20 patients from multilayer and single-layer models were 385.13 ± 110.09 kPa and 270.91 ± 95.86 kPa, respectively. The relative difference was 42.2%, with single-layer stress serving as the base value. The average mean plaque stress values from multilayer and single-layer models were 129.59 ± 32.77 kPa and 93.27 ± 18.20 kPa, respectively, with a relative difference of 38.9%. The maximum and mean plaque strain values obtained from the multilayer models were 11.6% and 19.0% higher than those from the single-layer models. Similarly, the maximum and mean cap strains showed increases of 9.6% and 12.9% over those from the single-layer models. These findings suggest that use of multilayer models could improve plaque stress and strain calculation accuracy and may have large impact on plaque progression and vulnerability investigation and potential clinical applications. Further large-scale studies are needed to validate our findings.

Automated finite element approach to generate anatomical patient-specific biomechanical models of atherosclerotic arteries from virtual histology-intravascular ultrasound

Despite advancements in early detection and treatment, atherosclerosis remains the leading cause of death across all cardiovascular diseases (CVD). Biomechanical analysis of atherosclerotic lesions has the potential to reveal biomechanically instable or rupture-prone regions. Treatment decisions rarely consider the biomechanics of the stenosed lesion due in-part to difficulties in obtaining this information in a clinical setting. Previous 3D FEA approaches have incompletely incorporated the complex curvature of arterial geometry, material heterogeneity, and use of patient-specific data. To address these limitations and clinical need, herein we present a user-friendly fully automated program to reconstruct and simulate the wall mechanics of patient-specific atherosclerotic coronary arteries. The program enables 3D reconstruction from patient-specific data with heterogenous tissue assignment and complex arterial curvature. Eleven arteries with coronary artery disease (CAD) underwent baseline and 6-month follow-up angiographic and virtual histology-intravascular ultrasound (VH-IVUS) imaging. VH-IVUS images were processed to remove background noise, extract VH plaque material data, and luminal and outer contours. Angiography data was used to orient the artery profiles along the 3D centerlines. The resulting surface mesh is then resampled for uniformity and tetrahedralized to generate the volumetric mesh using TetGen. A mesh convergence study revealed edge lengths between 0.04 mm and 0.2 mm produced constituent volumes that were largely unchanged, hence, to save computational resources, a value of 0.2 mm was used throughout. Materials are assigned and finite element analysis (FEA) is then performed to determine stresses and strains across the artery wall. In a representative artery, the highest average effective stress was in calcium elements with 235 kPa while necrotic elements had the lowest average stress, reaching as low as 0.79 kPa. After applying nodal smoothening, the maximum effective stress across 11 arteries remained below 288 kPa, implying biomechanically stable plaques. Indeed, all atherosclerotic plaques remained unruptured at the 6-month longitudinal follow up diagnosis. These results suggest our automated analysis may facilitate assessment of atherosclerotic plaque stability.

Image-Based Finite Element Modeling Approach for Characterizing In Vivo Mechanical Properties of Human Arteries

Mechanical properties of the arterial walls could provide meaningful information for the diagnosis, management and treatment of cardiovascular diseases. Classically, various experimental approaches were conducted on dissected arterial tissues to obtain their stress-stretch relationship, which has limited value clinically. Therefore, there is a pressing need to obtain biomechanical behaviors of these vascular tissues in vivo for personalized treatment. This paper reviews the methods to quantify arterial mechanical properties in vivo. Among these methods, we emphasize a novel approach using image-based finite element models to iteratively determine the material properties of the arterial tissues. This approach has been successfully applied to arterial walls in various vascular beds. The mechanical properties obtained from the in vivo approach were compared to those from ex vivo experimental studies to investigate whether any discrepancy in material properties exists for both approaches. Arterial tissue stiffness values from in vivo studies generally were in the same magnitude as those from ex vivo studies, but with lower average values. Some methodological issues, including solution uniqueness and robustness; method validation; and model assumptions and limitations were discussed. Clinical applications of this approach were also addressed to highlight their potential in translation from research tools to cardiovascular disease management.

Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases

Hemodynamic factors, induced by pulsatile blood flow, play a crucial role in vascular health and diseases, such as the initiation and progression of atherosclerosis. Computational fluid dynamics, finite element analysis, and fluid-structure interaction simulations have been widely used to quantify detailed hemodynamic forces based on vascular images commonly obtained from computed tomography angiography, magnetic resonance imaging, ultrasound, and optical coherence tomography. In this review, we focus on methods for obtaining accurate hemodynamic factors that regulate the structure and function of vascular endothelial and smooth muscle cells. We describe the multiple steps and recent advances in a typical patient-specific simulation pipeline, including medical imaging, image processing, spatial discretization to generate computational mesh, setting up boundary conditions and solver parameters, visualization and extraction of hemodynamic factors, and statistical analysis. These steps have not been standardized and thus have unavoidable uncertainties that should be thoroughly evaluated. We also discuss the recent development of combining patient-specific models with machine-learning methods to obtain hemodynamic factors faster and cheaper than conventional methods. These critical advances widen the use of biomechanical simulation tools in the research and potential personalized care of vascular diseases.