Fully Automated Construction of Three-dimensional Finite Element Simulations from Optical Coherence Tomography.

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

Fully automated construction of three-dimensional finite element simulations from Optical Coherence Tomography.

Despite recent advances in diagnosis and treatment, atherosclerotic coronary artery diseases remain a leading cause of death worldwide. Various imaging modalities and metrics can detect lesions and predict patients at risk; however, identifying unstable lesions is still difficult. Current techniques cannot fully capture the complex morphology-modulated mechanical responses that affect plaque stability, leading to catastrophic failure and mute the benefit of device and drug interventions. Finite Element (FE) simulations utilizing intravascular imaging OCT (Optical Coherence Tomography) are effective in defining physiological stress distributions. However, creating 3D FE simulations of coronary arteries from OCT images is challenging to fully automate given OCT frame sparsity, limited material contrast, and restricted penetration depth. To address such limitations, we developed an algorithmic approach to automatically produce 3D FE-ready digital twins from labeled OCT images. The 3D models are anatomically faithful and recapitulate mechanically relevant tissue lesion components, automatically producing morphologies structurally similar to manually constructed models whilst including more minute details. A mesh convergence study highlighted the ability to reach stress and strain convergence with average errors of just 5.9% and 1.6% respectively in comparison to FE models with approximately twice the number of elements in areas of refinement. Such an automated procedure will enable analysis of large clinical cohorts at a previously unattainable scale and opens the possibility for in-silico methods for patient specific diagnoses and treatment planning for coronary artery disease.

A platform for high-fidelity patient-specific structural modelling of atherosclerotic arteries: from intravascular imaging to three-dimensional stress distributions.

The pathophysiology of atherosclerotic lesions, including plaque rupture triggered by mechanical failure of the vessel wall, depends directly on the plaque morphology-modulated mechanical response. The complex interplay between lesion morphology and structural behaviour can be studied with high-fidelity computational modelling. However, construction of three-dimensional (3D) and heterogeneous models is challenging, with most previous work focusing on two-dimensional geometries or on single-material lesion compositions. Addressing these limitations, we here present a semi-automatic computational platform, leveraging clinical optical coherence tomography images to effectively reconstruct a 3D patient-specific multi-material model of atherosclerotic plaques, for which the mechanical response is obtained by structural finite-element simulations. To demonstrate the importance of including multi-material plaque components when recovering the mechanical response, a computational case study was conducted in which systematic variation of the intraplaque lipid and calcium was performed. The study demonstrated that the inclusion of various tissue components greatly affected the lesion mechanical response, illustrating the importance of multi-material formulations. This platform accordingly provides a viable foundation for studying how plaque micro-morphology affects plaque mechanical response, allowing for patient-specific assessments and extension into clinically relevant patient cohorts.

Biomechanics of diastolic dysfunction: a one-dimensional computational modeling approach.

Diastolic dysfunction (DD) is a major component of heart failure with preserved ejection fraction (HFpEF). Accordingly, a profound understanding of the underlying biomechanical mechanisms involved in DD is needed to elucidate all aspects of HFpEF. In this study, we have developed a computational model of DD by leveraging the power of an advanced one-dimensional arterial network coupled to a four-chambered zero-dimensional cardiac model. The two main pathologies investigated were linked to the active relaxation of the myocardium and the passive stiffness of the left ventricular wall. These pathologies were quantified through two parameters for the biphasic delay of active relaxation, which simulate the early and late-phase relaxation delay, and one parameter for passive stiffness, which simulates the increased nonlinear stiffness of the ventricular wall. A parameter sensitivity analysis was conducted on each of the three parameters to investigate their effect in isolation. The three parameters were then concurrently adjusted to produce the three main phenotypes of DD. It was found that the impaired relaxation phenotype can be replicated by mainly manipulating the active relaxation, the pseudo-normal phenotype was replicated by manipulating both the active relaxation and passive stiffness, and, finally, the restricted phenotype was replicated by mainly changing the passive stiffness. This article presents a simple model producing a holistic and comprehensive replication of the main DD phenotypes and presents novel biomechanical insights on how key parameters defining the relaxation and stiffness properties of the myocardium affect the development and manifestation of DD.NEW & NOTEWORTHY This study uses a complete and validated computational model of the cardiovascular system to simulate the two main pathologies involved in diastolic dysfunction (DD), i.e., abnormal active relaxation and increased ventricular diastolic stiffness. The three phenotypes of DD were successfully replicated according to literature data. We elucidate the biomechanical effect of the relaxation pathologies involved and how these pathologies interact to create the various phenotypes of DD.