1
|
Vuong TNAM, Bartolf‐Kopp M, Andelovic K, Jungst T, Farbehi N, Wise SG, Hayward C, Stevens MC, Rnjak‐Kovacina J. Integrating Computational and Biological Hemodynamic Approaches to Improve Modeling of Atherosclerotic Arteries. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2024; 11:e2307627. [PMID: 38704690 PMCID: PMC11234431 DOI: 10.1002/advs.202307627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Revised: 03/12/2024] [Indexed: 05/07/2024]
Abstract
Atherosclerosis is the primary cause of cardiovascular disease, resulting in mortality, elevated healthcare costs, diminished productivity, and reduced quality of life for individuals and their communities. This is exacerbated by the limited understanding of its underlying causes and limitations in current therapeutic interventions, highlighting the need for sophisticated models of atherosclerosis. This review critically evaluates the computational and biological models of atherosclerosis, focusing on the study of hemodynamics in atherosclerotic coronary arteries. Computational models account for the geometrical complexities and hemodynamics of the blood vessels and stenoses, but they fail to capture the complex biological processes involved in atherosclerosis. Different in vitro and in vivo biological models can capture aspects of the biological complexity of healthy and stenosed vessels, but rarely mimic the human anatomy and physiological hemodynamics, and require significantly more time, cost, and resources. Therefore, emerging strategies are examined that integrate computational and biological models, and the potential of advances in imaging, biofabrication, and machine learning is explored in developing more effective models of atherosclerosis.
Collapse
Affiliation(s)
| | - Michael Bartolf‐Kopp
- Department of Functional Materials in Medicine and DentistryInstitute of Functional Materials and Biofabrication (IFB)KeyLab Polymers for Medicine of the Bavarian Polymer Institute (BPI)University of WürzburgPleicherwall 297070WürzburgGermany
| | - Kristina Andelovic
- Department of Functional Materials in Medicine and DentistryInstitute of Functional Materials and Biofabrication (IFB)KeyLab Polymers for Medicine of the Bavarian Polymer Institute (BPI)University of WürzburgPleicherwall 297070WürzburgGermany
| | - Tomasz Jungst
- Department of Functional Materials in Medicine and DentistryInstitute of Functional Materials and Biofabrication (IFB)KeyLab Polymers for Medicine of the Bavarian Polymer Institute (BPI)University of WürzburgPleicherwall 297070WürzburgGermany
- Department of Orthopedics, Regenerative Medicine Center UtrechtUniversity Medical Center UtrechtUtrecht3584Netherlands
| | - Nona Farbehi
- Graduate School of Biomedical EngineeringUniversity of New South WalesSydney2052Australia
- Tyree Institute of Health EngineeringUniversity of New South WalesSydneyNSW2052Australia
- Garvan Weizmann Center for Cellular GenomicsGarvan Institute of Medical ResearchSydneyNSW2010Australia
| | - Steven G. Wise
- School of Medical SciencesUniversity of SydneySydneyNSW2006Australia
| | - Christopher Hayward
- St Vincent's HospitalSydneyVictor Chang Cardiac Research InstituteSydney2010Australia
| | | | - Jelena Rnjak‐Kovacina
- Graduate School of Biomedical EngineeringUniversity of New South WalesSydney2052Australia
- Tyree Institute of Health EngineeringUniversity of New South WalesSydneyNSW2052Australia
- Australian Centre for NanoMedicine (ACN)University of New South WalesSydneyNSW2052Australia
| |
Collapse
|
2
|
Ninno F, Chiastra C, Colombo M, Dardik A, Strosberg D, Aboian E, Tsui J, Bartlett M, Balabani S, Díaz-Zuccarini V. Modelling lower-limb peripheral arterial disease using clinically available datasets: impact of inflow boundary conditions on hemodynamic indices for restenosis prediction. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2024; 251:108214. [PMID: 38759252 DOI: 10.1016/j.cmpb.2024.108214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Revised: 04/22/2024] [Accepted: 05/01/2024] [Indexed: 05/19/2024]
Abstract
BACKGROUND AND OBJECTIVES The integration of hemodynamic markers as risk factors in restenosis prediction models for lower-limb peripheral arteries is hindered by fragmented clinical datasets. Computed tomography (CT) scans enable vessel geometry reconstruction and can be obtained at different times than the Doppler ultrasound (DUS) images, which provide information on blood flow velocity. Computational fluid dynamics (CFD) simulations allow the computation of near-wall hemodynamic indices, whose accuracy depends on the prescribed inlet boundary condition (BC), derived from the DUS images. This study aims to: (i) investigate the impact of different DUS-derived velocity waveforms on CFD results; (ii) test whether the same vessel areas, subjected to altered hemodynamics, can be detected independently of the applied inlet BC; (iii) suggest suitable DUS images to obtain reliable CFD results. METHODS CFD simulations were conducted on three patients treated with bypass surgery, using patient-specific DUS-derived inlet BCs recorded at either the same or different time points than the CT scan. The impact of the chosen inflow condition on bypass hemodynamics was assessed in terms of wall shear stress (WSS)-derived quantities. Patient-specific critical thresholds for the hemodynamic indices were applied to identify critical luminal areas and compare the results with a reference obtained with a DUS image acquired in close temporal proximity to the CT scan. RESULTS The main findings indicate that: (i) DUS-derived inlet velocity waveforms acquired at different time points than the CT scan led to statistically significantly different CFD results (p<0.001); (ii) the same luminal surface areas, exposed to low time-averaged WSS, could be identified independently of the applied inlet BCs; (iii) similar outcomes were observed for the other hemodynamic indices if the prescribed inlet velocity waveform had the same shape and comparable systolic acceleration time to the one recorded in close temporal proximity to the CT scan. CONCLUSIONS Despite a lack of standardised data collection for diseased lower-limb peripheral arteries, an accurate estimation of luminal areas subjected to altered near-wall hemodynamics is possible independently of the applied inlet BC. This holds if the applied inlet waveform shares some characteristics - derivable from the DUS report - as one matching the acquisition time of the CT scan.
Collapse
Affiliation(s)
- Federica Ninno
- Department of Medical Physics and Biomedical Engineering, University College London, London, UK; Wellcome-EPSRC Centre for Interventional Surgical Sciences, London, UK
| | - Claudio Chiastra
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Monika Colombo
- Department of Mechanical and Production Engineering, Aarhus University, Aarhus, Denmark
| | - Alan Dardik
- Vascular Biology and Therapeutics, Yale University School of Medicine, New Haven, Connecticut, USA; Division of Vascular Surgery and Endovascular Therapy, Department of Surgery, Yale University School of Medicine, New Haven, Connecticut, USA
| | - David Strosberg
- Department of Surgery, VA Connecticut Healthcare Systems, West Haven, Connecticut, USA; Department of Vascular Surgery, Royal Free Hospital NHS Foundation Trust, London, UK
| | - Edouard Aboian
- Department of Surgery, VA Connecticut Healthcare Systems, West Haven, Connecticut, USA
| | - Janice Tsui
- Department of Vascular Surgery, Royal Free Hospital NHS Foundation Trust, London, UK; Division of Surgery & Interventional Science, Department of Surgical Biotechnology, Faculty of Medical Sciences, University College London, London, UK
| | - Matthew Bartlett
- Division of Surgery & Interventional Science, Department of Surgical Biotechnology, Faculty of Medical Sciences, University College London, London, UK; Department of Mechanical Engineering, University College London, London, UK
| | - Stavroula Balabani
- Wellcome-EPSRC Centre for Interventional Surgical Sciences, London, UK; Department of Mechanical Engineering, University College London, London, UK
| | - Vanessa Díaz-Zuccarini
- Wellcome-EPSRC Centre for Interventional Surgical Sciences, London, UK; Department of Mechanical Engineering, University College London, London, UK.
| |
Collapse
|
3
|
Montino Pelagi G, Baggiano A, Regazzoni F, Fusini L, Alì M, Pontone G, Valbusa G, Vergara C. Personalized Pressure Conditions and Calibration for a Predictive Computational Model of Coronary and Myocardial Blood Flow. Ann Biomed Eng 2024; 52:1297-1312. [PMID: 38334838 PMCID: PMC10995040 DOI: 10.1007/s10439-024-03453-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 01/14/2024] [Indexed: 02/10/2024]
Abstract
Predictive modeling of hyperemic coronary and myocardial blood flow (MBF) greatly supports diagnosis and prognostic stratification of patients suffering from coronary artery disease (CAD). In this work, we propose a novel strategy, using only readily available clinical data, to build personalized inlet conditions for coronary and MBF models and to achieve an effective calibration for their predictive application to real clinical cases. Experimental data are used to build personalized pressure waveforms at the aortic root, representative of the hyperemic state and adapted to surrogate the systolic contraction, to be used in computational fluid-dynamics analyses. Model calibration to simulate hyperemic flow is performed in a "blinded" way, not requiring any additional exam. Coronary and myocardial flow simulations are performed in eight patients with different clinical conditions to predict FFR and MBF. Realistic pressure waveforms are recovered for all the patients. Consistent pressure distribution, blood velocities in the large arteries, and distribution of MBF in the healthy myocardium are obtained. FFR results show great accuracy with a per-vessel sensitivity and specificity of 100% according to clinical threshold values. Mean MBF shows good agreement with values from stress-CTP, with lower values in patients with diagnosed perfusion defects. The proposed methodology allows us to quantitatively predict FFR and MBF, by the exclusive use of standard measures easily obtainable in a clinical context. This represents a fundamental step to avoid catheter-based exams and stress tests in CAD diagnosis.
Collapse
Affiliation(s)
- Giovanni Montino Pelagi
- LABS, Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico di Milano, 20133, Milan, Italy.
| | - Andrea Baggiano
- Perioperative Cardiology and Cardiovascular Imaging Department, Centro Cardiologico Monzino IRCCS, Via Carlo Parea 4, 20138, Milan, Italy
- Department of Clinical Sciences and Community Health, University of Milan, Milan, Italy
| | - Francesco Regazzoni
- MOX, Dipartimento di Matematica, Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milan, Italy
| | - Laura Fusini
- Perioperative Cardiology and Cardiovascular Imaging Department, Centro Cardiologico Monzino IRCCS, Via Carlo Parea 4, 20138, Milan, Italy
- Department of Electronics, Information and Biomedical Engineering, Politecnico di Milano, 20133, Milan, Italy
| | - Marco Alì
- Bracco Imaging S.p.A., Via Caduti di Marcinelle 13, 20134, Milan, Italy
- Department of Diagnostic Imaging and Stereotactic Radiosurgery, Centro Diagnostico Italiano S.p.A., Via Saint Bon 20, 20147, Milan, Italy
| | - Gianluca Pontone
- Perioperative Cardiology and Cardiovascular Imaging Department, Centro Cardiologico Monzino IRCCS, Via Carlo Parea 4, 20138, Milan, Italy
- Department of Biomedical, Surgical and Dental Sciences, University of Milan, 20134, Milan, Italy
| | - Giovanni Valbusa
- Bracco Imaging S.p.A., Via Caduti di Marcinelle 13, 20134, Milan, Italy
| | - Christian Vergara
- LABS, Dipartimento di Chimica, Materiali e Ingegneria Chimica, Politecnico di Milano, 20133, Milan, Italy
| |
Collapse
|
4
|
Menon K, Khan MO, Sexton ZA, Richter J, Nguyen PK, Malik SB, Boyd J, Nieman K, Marsden AL. Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees. NPJ IMAGING 2024; 2:9. [PMID: 38706558 PMCID: PMC11062925 DOI: 10.1038/s44303-024-00014-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 02/25/2024] [Indexed: 05/07/2024]
Abstract
Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary anatomies combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. We propose an optimization-based method to personalize multiscale coronary flow simulations by assimilating clinical CT myocardial perfusion imaging and cardiac function measurements to yield patient-specific flow distributions and model parameters. Using this proof-of-concept study on a cohort of six patients, we reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based purely on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.
Collapse
Affiliation(s)
- Karthik Menon
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA USA
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA USA
| | - Muhammed Owais Khan
- Department of Electrical, Computer, and Biomedical Engineering, Toronto Metropolitan University, Toronto, ON Canada
| | | | - Jakob Richter
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA USA
| | - Patricia K. Nguyen
- VA Palo Alto Healthcare System, Palo Alto, CA USA
- Division of Cardiovascular Medicine, Stanford School of Medicine, Stanford, CA USA
| | | | - Jack Boyd
- Department of Cardiothoracic Surgery, Stanford School of Medicine, Stanford, CA USA
| | - Koen Nieman
- Division of Cardiovascular Medicine, Stanford School of Medicine, Stanford, CA USA
- Department of Radiology, Stanford School of Medicine, Stanford, CA USA
| | - Alison L. Marsden
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA USA
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA USA
- Department of Bioengineering, Stanford University, Stanford, CA USA
| |
Collapse
|
5
|
Lodi Rizzini M, Candreva A, Mazzi V, Pagnoni M, Chiastra C, Aben JP, Fournier S, Cook S, Muller O, De Bruyne B, Mizukami T, Collet C, Gallo D, Morbiducci U. Blood Flow Energy Identifies Coronary Lesions Culprit of Future Myocardial Infarction. Ann Biomed Eng 2024; 52:226-238. [PMID: 37733110 DOI: 10.1007/s10439-023-03362-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2023] [Accepted: 09/02/2023] [Indexed: 09/22/2023]
Abstract
The present study establishes a link between blood flow energy transformations in coronary atherosclerotic lesions and clinical outcomes. The predictive capacity for future myocardial infarction (MI) was compared with that of established quantitative coronary angiography (QCA)-derived predictors. Angiography-based computational fluid dynamics (CFD) simulations were performed on 80 human coronary lesions culprit of MI within 5 years and 108 non-culprit lesions for future MI. Blood flow energy transformations were assessed in the converging flow segment of the lesion as ratios of kinetic and rotational energy values (KER and RER, respectively) at the QCA-identified minimum lumen area and proximal lesion sections. The anatomical and functional lesion severity were evaluated with QCA to derive percentage area stenosis (%AS), vessel fractional flow reserve (vFFR), and translesional vFFR (ΔvFFR). Wall shear stress profiles were investigated in terms of topological shear variation index (TSVI). KER and RER predicted MI at 5 years (AUC = 0.73, 95% CI 0.65-0.80, and AUC = 0.76, 95% CI 0.70-0.83, respectively; p < 0.0001 for both). The predictive capacity for future MI of KER and RER was significantly stronger than vFFR (p = 0.0391 and p = 0.0045, respectively). RER predictive capacity was significantly stronger than %AS and ΔvFFR (p = 0.0041 and p = 0.0059, respectively). The predictive capacity for future MI of KER and RER did not differ significantly from TSVI. Blood flow kinetic and rotational energy transformations were significant predictors for MI at 5 years (p < 0.0001). The findings of this study support the hypothesis of a biomechanical contribution to the process of plaque destabilization/rupture leading to MI.
Collapse
Affiliation(s)
- Maurizio Lodi Rizzini
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
| | - Alessandro Candreva
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
- Department of Cardiology, Zurich University Hospital, Zurich, Switzerland
| | - Valentina Mazzi
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
| | - Mattia Pagnoni
- Department of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Claudio Chiastra
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
| | | | - Stephane Fournier
- Department of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
- Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Stephane Cook
- Department of Cardiology, HFR Fribourg, Fribourg, Switzerland
| | - Olivier Muller
- Department of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | | | | | - Carlos Collet
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium
| | - Diego Gallo
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
| | - Umberto Morbiducci
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy.
| |
Collapse
|
6
|
Shi Y, Zheng J, Zhang Y, Sun Q, Shen J, Gao Y, Sun J, Yang N, Zhou X, Li S, Weir-McCall JR, Xia P, Teng Z. The influence of flow distribution strategy for the quantification of pressure- and wall shear stress-derived parameters in the coronary artery: A CTA-based computational fluid dynamics analysis. J Biomech 2023; 161:111857. [PMID: 37939424 DOI: 10.1016/j.jbiomech.2023.111857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2023] [Revised: 10/15/2023] [Accepted: 11/01/2023] [Indexed: 11/10/2023]
Abstract
For image-based computational fluid dynamics (CFD) analysis to characterize the local coronary hemodynamic environment, the accuracy depends on the flow rate which is in turn associated with outlet branches' morphology. A good flow distribution strategy is important to mitigate the effect when certain branches cannot be considered. In this study, stenotic coronary arteries from 13 patients were used to analyze the effect of missing branches and different flow distribution strategies. Pressure- and wall shear stress (WSS)-derived parameters around the stenotic region (ROI) were compared, including fractional flow reserve (CT-FFR), instantaneous wave-free ratio (CT-iFR), resting distal to aortic coronary pressure (CT-Pd/Pa), time-averaged WSS, oscillatory shear index (OSI) and relative residence time (RRT). Three flow distribution strategies were the Huo-Kassab model at distal outlets (Type I), flow distribution based on outlet resistances (Type II), and a developed algorithm distributing flow at each bifurcation until the final outlets (Type III). Results showed that Type III strategy for models with truncated branch(es) had a good agreement in both pressure- and WSS-related results (interquatile range less than 0.12% and 4.02%, respectively) with the baseline model around the ROI. The relative difference of pressure- and WSS-related results were correlated with the flow differences in the ROI to the baseline mode. Type III strategy had the best performance in maintaining the flow in intermediate branches. It is recommended for CFD analysis. Removal of branches distal to a stenosis can be undertaken with an improved performance and maintained accuracy, while those proximal to the ROI should be kept.
Collapse
Affiliation(s)
- Yibing Shi
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Jin Zheng
- Department of Radiology, University of Cambridge, UK
| | - Ying Zhang
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Quanlin Sun
- Department of Radiology, University of Cambridge, UK; Nanjing Jingsan Medical Science and Technology, Ltd., Jiangsu, China
| | - Jinhua Shen
- Nanjing Jingsan Medical Science and Technology, Ltd., Jiangsu, China
| | - Yongguang Gao
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Jingxi Sun
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Ning Yang
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Xuanxuan Zhou
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Suqing Li
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China
| | - Jonathan R Weir-McCall
- Department of Radiology, University of Cambridge, UK; Department of Radiology, Royal Papworth Hospital, Cambridge, UK
| | - Ping Xia
- Department of Radiology, Xuzhou Central Hospital, Xuzhou Clinical School of Xuzhou Medical University, Jiangsu, China.
| | - Zhongzhao Teng
- Department of Radiology, University of Cambridge, UK; Nanjing Jingsan Medical Science and Technology, Ltd., Jiangsu, China.
| |
Collapse
|
7
|
De Nisco G, Lodi Rizzini M, Verardi R, Chiastra C, Candreva A, De Ferrari G, D'Ascenzo F, Gallo D, Morbiducci U. Modelling blood flow in coronary arteries: Newtonian or shear-thinning non-Newtonian rheology? COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 242:107823. [PMID: 37757568 DOI: 10.1016/j.cmpb.2023.107823] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 09/14/2023] [Accepted: 09/18/2023] [Indexed: 09/29/2023]
Abstract
BACKGROUND The combination of medical imaging and computational hemodynamics is a promising technology to diagnose/prognose coronary artery disease (CAD). However, the clinical translation of in silico hemodynamic models is still hampered by assumptions/idealizations that must be introduced in model-based strategies and that necessarily imply uncertainty. This study aims to provide a definite answer to the open question of how to properly model blood rheological properties in computational fluid dynamics (CFD) simulations of coronary hemodynamics. METHODS The geometry of the right coronary artery (RCA) of 144 hemodynamically stable patients with different stenosis degree were reconstructed from angiography. On them, unsteady-state CFD simulations were carried out. On each reconstructed RCA two different simulation strategies were applied to account for blood rheological properties, implementing (i) a Newtonian (N) and (ii) a shear-thinning non-Newtonian (non-N) rheological model. Their impact was evaluated in terms of wall shear stress (WSS magnitude, multidirectionality, topological skeleton) and helical flow (strength, topology) profiles. Additionally, luminal surface areas (SAs) exposed to shear disturbances were identified and the co-localization of paired N and non-N SAs was quantified in terms of similarity index (SI). RESULTS The comparison between paired N vs. shear-thinning non-N simulations revealed remarkably similar profiles of WSS-based and helicity-based quantities, independent of the adopted blood rheology model and of the degree of stenosis of the vessel. Statistically, for each paired N and non-N hemodynamic quantity emerged negligible bias from Bland-Altman plots, and strong positive linear correlation (r > 0.94 for almost all the WSS-based quantities, r > 0.99 for helicity-based quantities). Moreover, a remarkable co-localization of N vs. non-N luminal SAs exposed to disturbed shear clearly emerged (SI distribution 0.95 [0.93, 0.97]). Helical flow topology resulted to be unaffected by blood rheological properties. CONCLUSIONS This study, performed on 288 angio-based CFD simulations on 144 RCA models presenting with different degrees of stenosis, suggests that the assumptions on blood rheology have negligible impact both on WSS and helical flow profiles associated with CAD, thus definitively answering to the question "is Newtonian assumption for blood rheology adequate in coronary hemodynamics simulations?".
Collapse
Affiliation(s)
- Giuseppe De Nisco
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Maurizio Lodi Rizzini
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Roberto Verardi
- Hemodynamic Laboratory, Department of Medical Sciences, University of Turin, Turin, Italy
| | - Claudio Chiastra
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Alessandro Candreva
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Gaetano De Ferrari
- Hemodynamic Laboratory, Department of Medical Sciences, University of Turin, Turin, Italy
| | - Fabrizio D'Ascenzo
- Hemodynamic Laboratory, Department of Medical Sciences, University of Turin, Turin, Italy
| | - Diego Gallo
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy.
| | - Umberto Morbiducci
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| |
Collapse
|
8
|
Menon K, Khan MO, Sexton ZA, Richter J, Nieman K, Marsden AL. Personalized coronary and myocardial blood flow models incorporating CT perfusion imaging and synthetic vascular trees. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.08.17.23294242. [PMID: 37645850 PMCID: PMC10462196 DOI: 10.1101/2023.08.17.23294242] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/31/2023]
Abstract
Computational simulations of coronary artery blood flow, using anatomical models based on clinical imaging, are an emerging non-invasive tool for personalized treatment planning. However, current simulations contend with two related challenges - incomplete anatomies in image-based models due to the exclusion of arteries smaller than the imaging resolution, and the lack of personalized flow distributions informed by patient-specific imaging. We introduce a data-enabled, personalized and multi-scale flow simulation framework spanning large coronary arteries to myocardial microvasculature. It includes image-based coronary models combined with synthetic vasculature for arteries below the imaging resolution, myocardial blood flow simulated using Darcy models, and systemic circulation represented as lumped-parameter networks. Personalized flow distributions and model parameters are informed by clinical CT myocardial perfusion imaging and cardiac function using surrogate-based optimization. We reveal substantial differences in flow distributions and clinical diagnosis metrics between the proposed personalized framework and empirical methods based on anatomy; these errors cannot be predicted a priori. This suggests virtual treatment planning tools would benefit from increased personalization informed by emerging imaging methods.
Collapse
Affiliation(s)
- Karthik Menon
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA, USA
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA, USA
| | - Muhammed Owais Khan
- Department of Electrical, Computer, and Biomedical Engineering, Toronto Metropolitan University, Toronto, Ontario, Canada
| | - Zachary A Sexton
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| | - Jakob Richter
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA, USA
| | - Koen Nieman
- Departments of Radiology and Medicine (Cardiovascular Medicine), Stanford School of Medicine, Stanford, CA, USA
| | - Alison L Marsden
- Department of Pediatrics (Cardiology), Stanford School of Medicine, Stanford, CA, USA
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA, USA
- Department of Bioengineering, Stanford University, Stanford, CA, USA
| |
Collapse
|
9
|
Candreva A, Rizzini ML, Schweiger V, Gallo D, Montone RA, Würdinger M, Stehli J, Gilhofer T, Gotschy A, Frank R, Stähli BE, Chiastra C, Morbiducci U, Templin C. Is spontaneous coronary artery dissection (SCAD) related to local anatomy and hemodynamics? An exploratory study. Int J Cardiol 2023:S0167-5273(23)00657-5. [PMID: 37201616 DOI: 10.1016/j.ijcard.2023.05.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Revised: 04/16/2023] [Accepted: 05/05/2023] [Indexed: 05/20/2023]
Abstract
AIMS Spontaneous coronary artery dissection (SCAD) is an increasingly diagnosed cause of myocardial infarction with unclear pathophysiology. The aim of the study was to test if vascular segments site of SCAD present distinctive local anatomy and hemodynamic profiles. METHODS Coronary arteries with spontaneously healed SCAD (confirmed by follow-up angiography) underwent three-dimensional reconstruction, morphometric analysis with definition of vessel local curvature and torsion, and computational fluid dynamics (CFD) simulations with derivation of time-averaged wall shear stress (TAWSS) and topological shear variation index (TSVI). The (reconstructed) healed proximal SCAD segment was visually inspected for co-localization with curvature, torsion, and CFD-derived quantities hot spots. RESULTS Thirteen vessels with healed SCAD underwent the morpho-functional analysis. Median time between baseline and follow-up coronary angiograms was 57 (interquartile range [IQR] 45-95) days. In seven cases (53.9%), SCAD was classified as type 2b and occurred in the left anterior descending artery or near a bifurcation. In all cases (100%), at least one hot spot co-localized within the healed proximal SCAD segment, in 9 cases (69.2%) ≥3 hot spots were identified. Healed SCAD in proximity of a coronary bifurcation presented lower TAWSS peak values (6.65 [IQR 6.20-13.2] vs. 3.81 [2.53-5.17] Pa, p = 0.008) and hosted less frequently TSVI hot spots (100% vs. 57.1%, p = 0.034). CONCLUSION Vascular segments of healed SCAD were characterized by high curvature/torsion and WSS profiles reflecting increased local flow disturbances. Hence, a pathophysiological role of the interaction between vessel anatomy and shear forces in SCAD is hypothesized.
Collapse
Affiliation(s)
- Alessandro Candreva
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland; PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Maurizio Lodi Rizzini
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Victor Schweiger
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland
| | - Diego Gallo
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Rocco A Montone
- Department of Cardiovascular Medicine, Fondazione Policlinico Universitario A. Gemelli IRCCS, Rome, Italy
| | - Michael Würdinger
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland
| | - Julia Stehli
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland
| | - Thomas Gilhofer
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland
| | - Alexander Gotschy
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland
| | - Ruschitzka Frank
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland; University of Zurich, Zurich, Switzerland
| | - Barbara E Stähli
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland; University of Zurich, Zurich, Switzerland
| | - Claudio Chiastra
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Umberto Morbiducci
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Christian Templin
- Department of Cardiology, University Heart Center, Zurich University Hospital, Zurich, Switzerland; University of Zurich, Zurich, Switzerland.
| |
Collapse
|
10
|
Carpenter HJ, Ghayesh MH, Zander AC, Psaltis PJ. On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 231:107418. [PMID: 36842347 DOI: 10.1016/j.cmpb.2023.107418] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 02/01/2023] [Accepted: 02/09/2023] [Indexed: 06/18/2023]
Abstract
BACKGROUND AND OBJECTIVE In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied. METHODS We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology. RESULTS 22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately. CONCLUSIONS These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.
Collapse
Affiliation(s)
- Harry J Carpenter
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia.
| | - Mergen H Ghayesh
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia.
| | - Anthony C Zander
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Peter J Psaltis
- Vascular Research Centre, Lifelong Health Theme, South Australian Health and Medical Research Institute (SAHMRI), Adelaide, South Australia 5000, Australia; Adelaide Medical School, University of Adelaide, Adelaide, South Australia 5005, Australia; Department of Cardiology, Central Adelaide Local Health Network, Adelaide, South Australia 5000, Australia
| |
Collapse
|