1
|
Bertram CD, Macaskill C. Fluid-Dynamic Modeling of Flow in Embryonic Tissue Indicates That Lymphatic Valve Location Is Not Consistently Determined by the Local Fluid Shear or Its Gradient. Microcirculation 2024; 31:e12873. [PMID: 38953384 PMCID: PMC11303113 DOI: 10.1111/micc.12873] [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: 03/04/2024] [Revised: 05/02/2024] [Accepted: 06/08/2024] [Indexed: 07/04/2024]
Abstract
OBJECTIVE Intravascular lymphatic valves often occur in proximity to vessel junctions. It is commonly held that disturbed flow at junctions is responsible for accumulation of valve-forming cells (VFCs) at these locations as the initial step in valve creation, and the one which explains the association with these sites. However, evidence in favor is largely limited to cell culture experiments. METHODS We acquired images of embryonic lymphatic vascular networks from day E16.5, when VFC accumulation has started but the developing valve has not yet altered the local vessel geometry, stained for Prox1, which co-localizes with Foxc2. Using finite-element computational fluid mechanics, we simulated the flow through the networks, under conditions appropriate to this early development stage. Then we correlated the Prox1 distributions with the distributions of simulated fluid shear and shear stress gradient. RESULTS Across a total of 16 image sets, no consistent correlation was found between Prox1 distribution and the local magnitude of fluid shear, or its positive or negative gradient. CONCLUSIONS This, the first direct semi-empirical test of the localization hypothesis to interrogate the tissue from in vivo at the critical moment of development, does not support the idea that a feature of the local flow determines valve localization.
Collapse
Affiliation(s)
- Christopher D Bertram
- School of Mathematics and Statistics, University of Sydney, New South Wales, Australia
| | - Charlie Macaskill
- School of Mathematics and Statistics, University of Sydney, New South Wales, Australia
| |
Collapse
|
2
|
Goetz A, Jeken-Rico P, Pelissier U, Chau Y, Sédat J, Hachem E. AnXplore: a comprehensive fluid-structure interaction study of 101 intracranial aneurysms. Front Bioeng Biotechnol 2024; 12:1433811. [PMID: 39007055 PMCID: PMC11243300 DOI: 10.3389/fbioe.2024.1433811] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2024] [Accepted: 06/03/2024] [Indexed: 07/16/2024] Open
Abstract
Advances in computational fluid dynamics continuously extend the comprehension of aneurysm growth and rupture, intending to assist physicians in devising effective treatment strategies. While most studies have first modelled intracranial aneurysm walls as fully rigid with a focus on understanding blood flow characteristics, some researchers further introduced Fluid-Structure Interaction (FSI) and reported notable haemodynamic alterations for a few aneurysm cases when considering wall compliance. In this work, we explore further this research direction by studying 101 intracranial sidewall aneurysms, emphasizing the differences between rigid and deformable-wall simulations. The proposed dataset along with simulation parameters are shared for the sake of reproducibility. A wide range of haemodynamic patterns has been statistically analyzed with a particular focus on the impact of the wall modelling choice. Notable deviations in flow characteristics and commonly employed risk indicators are reported, particularly with near-dome blood recirculations being significantly impacted by the pulsating dynamics of the walls. This leads to substantial fluctuations in the sac-averaged oscillatory shear index, ranging from -36% to +674% of the standard rigid-wall value. Going a step further, haemodynamics obtained when simulating a flow-diverter stent modelled in conjunction with FSI are showcased for the first time, revealing a 73% increase in systolic sac-average velocity for the compliant-wall setting compared to its rigid counterpart. This last finding demonstrates the decisive impact that FSI modelling can have in predicting treatment outcomes.
Collapse
Affiliation(s)
- Aurèle Goetz
- Computing and Fluids Research Group, CEMEF, Mines Paris PSL, Sophia Antipolis, France
| | - Pablo Jeken-Rico
- Computing and Fluids Research Group, CEMEF, Mines Paris PSL, Sophia Antipolis, France
| | - Ugo Pelissier
- Computing and Fluids Research Group, CEMEF, Mines Paris PSL, Sophia Antipolis, France
| | - Yves Chau
- Department of Neuro-Interventional and Vascular Interventional, University Hospital of Nice, Nice, France
| | - Jacques Sédat
- Department of Neuro-Interventional and Vascular Interventional, University Hospital of Nice, Nice, France
| | - Elie Hachem
- Computing and Fluids Research Group, CEMEF, Mines Paris PSL, Sophia Antipolis, France
| |
Collapse
|
3
|
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
|
4
|
Balasubramanya A, Maes L, Rega F, Mazzi V, Morbiducci U, Famaey N, Degroote J, Segers P. Hemodynamics and wall shear metrics in a pulmonary autograft: Comparing a fluid-structure interaction and computational fluid dynamics approach. Comput Biol Med 2024; 176:108604. [PMID: 38761502 DOI: 10.1016/j.compbiomed.2024.108604] [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: 01/23/2024] [Revised: 05/02/2024] [Accepted: 05/11/2024] [Indexed: 05/20/2024]
Abstract
OBJECTIVE In young patients, aortic valve disease is often treated by placement of a pulmonary autograft (PA) which adapts to its new environment through growth and remodeling. To better understand the hemodynamic forces acting on the highly distensible PA in the acute phase after surgery, we developed a fluid-structure interaction (FSI) framework and comprehensively compared hemodynamics and wall shear-stress (WSS) metrics with a computational fluid dynamic (CFD) simulation. METHODS The FSI framework couples a prestressed non-linear hyperelastic arterial tissue model with a fluid model using the in-house coupling code CoCoNuT. Geometry, material parameters and boundary conditions are based on in-vivo measurements. Hemodynamics, time-averaged WSS (TAWSS), oscillatory shear index (OSI) and topological shear variation index (TSVI) are evaluated qualitatively and quantitatively for 3 different sheeps. RESULTS Despite systolic-to-diastolic volumetric changes of the PA in the order of 20 %, the point-by-point correlation of TAWSS and OSI obtained through CFD and FSI remains high (r > 0.9, p < 0.01) for TAWSS and (r > 0.8, p < 0.01) for OSI). Instantaneous WSS divergence patterns qualitatively preserve similarities, but large deformations of the PA leads to a decrease of the correlation between FSI and CFD resolved TSVI (r < 0.7, p < 0.01). Moderate co-localization between FSI and CFD is observed for low thresholds of TAWSS and high thresholds of OSI and TSVI. CONCLUSION FSI might be warranted if we were to use the TSVI as a mechano-biological driver for growth and remodeling of PA due to varying intra-vascular flow structures and near wall hemodynamics because of the large expansion of the PA.
Collapse
Affiliation(s)
| | - Lauranne Maes
- Department of Mechanical Engineering, KU Leuven, Leuven, Belgium
| | - Filip Rega
- Cardiac Surgery, Department of Cardiovascular Sciences, KU Leuven, Belgium
| | - Valentina Mazzi
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Umberto Morbiducci
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Nele Famaey
- Department of Mechanical Engineering, KU Leuven, Leuven, Belgium
| | - Joris Degroote
- Department of Electromechanical Systems and Metal Engineering, Ghent University, Ghent, Belgium
| | | |
Collapse
|
5
|
Candreva A, Gallo D, Munhoz D, Rizzini ML, Mizukami T, Seki R, Sakai K, Sonck J, Mazzi V, Ko B, Nørgaard BL, Jensen JM, Maeng M, Otake H, Koo BK, Shinke T, Aben JP, Andreini D, Gallinoro E, Stähli BE, Templin C, Chiastra C, De Bruyne B, Morbiducci U, Collet C. Influence of intracoronary hemodynamic forces on atherosclerotic plaque phenotypes. Int J Cardiol 2024; 399:131668. [PMID: 38141723 DOI: 10.1016/j.ijcard.2023.131668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 10/21/2023] [Accepted: 12/18/2023] [Indexed: 12/25/2023]
Abstract
BACKGROUND AND AIMS Coronary hemodynamics impact coronary plaque progression and destabilization. The aim of the present study was to establish the association between focal vs. diffuse intracoronary pressure gradients and wall shear stress (WSS) patterns with atherosclerotic plaque composition. METHODS Prospective, international, single-arm study of patients with chronic coronary syndromes and hemodynamic significant lesions (fractional flow reserve [FFR] ≤ 0.80). Motorized FFR pullback pressure gradient (PPG), optical coherence tomography (OCT), and time-average WSS (TAWSS) and topological shear variation index (TSVI) derived from three-dimensional angiography were obtained. RESULTS One hundred five vessels (median FFR 0.70 [Interquartile range (IQR) 0.56-0.77]) had combined PPG and WSS analyses. TSVI was correlated with PPG (r = 0.47, [95% Confidence Interval (95% CI) 0.30-0.65], p < 0.001). Vessels with a focal CAD (PPG above the median value of 0.67) had significantly higher TAWSS (14.8 [IQR 8.6-24.3] vs. 7.03 [4.8-11.7] Pa, p < 0.001) and TSVI (163.9 [117.6-249.2] vs. 76.8 [23.1-140.9] m-1, p < 0.001). In the 51 vessels with baseline OCT, TSVI was associated with plaque rupture (OR 1.01 [1.00-1.02], p = 0.024), PPG with the extension of lipids (OR 7.78 [6.19-9.77], p = 0.003), with the presence of thin-cap fibroatheroma (OR 2.85 [1.11-7.83], p = 0.024) and plaque rupture (OR 4.94 [1.82 to 13.47], p = 0.002). CONCLUSIONS Focal and diffuse coronary artery disease, defined using coronary physiology, are associated with differential WSS profiles. Pullback pressure gradients and WSS profiles are associated with atherosclerotic plaque phenotypes. Focal disease (as identified by high PPG) and high TSVI are associated with high-risk plaque features. CLINICAL TRIAL REGISTRATION https://clinicaltrials,gov/ct2/show/NCT03782688.
Collapse
Affiliation(s)
- Alessandro Candreva
- Department of Cardiology, Zurich University Hospital, Zurich, Switzerland; Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Diego Gallo
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Daniel Munhoz
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy; Department of internal medicine, University of Campinas (Unicamp), Campinas, Brazil
| | - Maurizio Lodi Rizzini
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Takuya Mizukami
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Cardiology, Showa University School of Medicine, Tokyo, Japan
| | - Ruiko Seki
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium
| | - Koshiro Sakai
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium
| | - Jeroen Sonck
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Valentina Mazzi
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Brian Ko
- Monash Cardiovascular Research Centre, Monash University and Monash Heart, Monash Health, Clayton, Victoria, Australia
| | | | | | - Michael Maeng
- Department of Cardiology, Aarhus University Hospital, Aarhus, Denmark
| | - Hiromasa Otake
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Cardiology, Aichi Medical University, Aichi, Japan
| | - Bon-Kwon Koo
- Department of Internal Medicine and Cardiovascular Center, Seoul National University Hospital, Seoul, South Korea
| | - Toshiro Shinke
- Department of Cardiology, Showa University School of Medicine, Tokyo, Japan
| | | | - Daniele Andreini
- Department of Cardiology, IRCCS Ospedale Galeazzi-Sant'Ambrogio, Milan, Italy and Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy
| | - Emanuele Gallinoro
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Cardiology, IRCCS Ospedale Galeazzi-Sant'Ambrogio, Milan, Italy and Department of Biomedical and Clinical Sciences, University of Milan, Milan, Italy
| | - Barbara E Stähli
- Department of Cardiology, Zurich University Hospital, Zurich, Switzerland; University of Zurich, Zurich, Switzerland
| | - Christian Templin
- Department of Cardiology, Zurich University Hospital, Zurich, Switzerland
| | - Claudio Chiastra
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Bernard De Bruyne
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Department of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Umberto Morbiducci
- PoliTo(BIO) Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Carlos Collet
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium.
| |
Collapse
|
6
|
Ledda PG, Rossi T, Badas MG, Querzoli G. Can wall shear-stress topology predict proliferative vitreoretinopathy localization following pars plana vitrectomy? J Biomech 2024; 162:111914. [PMID: 38157782 DOI: 10.1016/j.jbiomech.2023.111914] [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: 09/21/2023] [Revised: 12/08/2023] [Accepted: 12/20/2023] [Indexed: 01/03/2024]
Abstract
We numerically study the fluid dynamics of oil tamponade in models of vitrectomized eyes prompted by a subset of daily activities corresponding to movements on the horizontal plane with the patient in a standing position. Bulk flow features are related to near-wall flow topology and transport at the retinal surface through a wall shear-stress-based analysis. Proliferative VitreoRetinopathy (PVR) is the leading cause of retinal re-detachment occurring in about 20% of all cases due to the accumulation of inflammatory cells in discrete retinal regions. Signalling soluble mediators stimulate inflammatory cells' chemotaxis and studying their distribution across the retinal surface may acquire clinical relevance. In all the investigated cases, persistent and elongated regions along the retina, potentially prone to accumulate chemo-attractants and cells are observed. Gradients of soluble inflammation mediators present in the aqueous are known responsible for the so-called epithelial-mesenchymal transition that initiates PVR and favours recurrent retinal detachment prompting the proliferation of inflammatory cells with collagen matrix deposition and its contraction. The surgical apposition of encircling scleral buckling elements, known for over a century to influence PVR formation and localization, modifies the attracting regions, possibly causing an accumulation of molecules and cells along approximately vertical lines that follow the rising menisci due to the cerclage indentation. The resulting spatial pattern is compatible with clinical observations. This study may open toward rational analyses of near-wall transport to predict PVR pathogenesis by relating biochemical accumulation in certain areas of the retina to clinical conditions.
Collapse
Affiliation(s)
- P G Ledda
- DICAAR, Università degli Studi di Cagliari, Cagliari, Italy
| | - T Rossi
- IRCCS Fondazione G.B. Bietti ONLUS, Rome, Italy
| | - M G Badas
- DICAAR, Università degli Studi di Cagliari, Cagliari, Italy
| | - G Querzoli
- DICAAR, Università degli Studi di Cagliari, Cagliari, Italy.
| |
Collapse
|
7
|
Tsolaki E, Corso P, Zboray R, Avaro J, Appel C, Liebi M, Bertazzo S, Heinisch PP, Carrel T, Obrist D, Herrmann IK. Multiscale multimodal characterization and simulation of structural alterations in failed bioprosthetic heart valves. Acta Biomater 2023; 169:138-154. [PMID: 37517619 DOI: 10.1016/j.actbio.2023.07.044] [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: 02/27/2023] [Revised: 06/30/2023] [Accepted: 07/24/2023] [Indexed: 08/01/2023]
Abstract
Calcific degeneration is the most frequent type of heart valve failure, with rising incidence due to the ageing population. The gold standard treatment to date is valve replacement. Unfortunately, calcification oftentimes re-occurs in bioprosthetic substitutes, with the governing processes remaining poorly understood. Here, we present a multiscale, multimodal analysis of disturbances and extensive mineralisation of the collagen network in failed bioprosthetic bovine pericardium valve explants with full histoanatomical context. In addition to highly abundant mineralized collagen fibres and fibrils, calcified micron-sized particles previously discovered in native valves were also prevalent on the aortic as well as the ventricular surface of bioprosthetic valves. The two mineral types (fibres and particles) were detectable even in early-stage mineralisation, prior to any macroscopic calcification. Based on multiscale multimodal characterisation and high-fidelity simulations, we demonstrate that mineral occurrence coincides with regions exposed to high haemodynamic and biomechanical indicators. These insights obtained by multiscale analysis of failed bioprosthetic valves serve as groundwork for the evidence-based development of more durable alternatives. STATEMENT OF SIGNIFICANCE: Bioprosthetic valve calcification is a well-known clinically significant phenomenon, leading to valve failure. The nanoanalytical characterisation of bioprosthetic valves gives insights into the highly abundant, extensive calcification and disorganization of the collagen network and the presence of calcium phosphate particles previously reported in native cardiovascular tissues. While the collagen matrix mineralisation can be primarily attributed to a combination of chemical and mechanical alterations, the calcified particles are likely of host cellular origin. This work presents a straightforward route to mineral identification and characterization at high resolution and sensitivity, and with full histoanatomical context and correlation to hemodynamic and biomechanical indicators, hence providing design cues for improved bioprosthetic valve alternatives.
Collapse
Affiliation(s)
- Elena Tsolaki
- Laboratory for Particles-Biology Interactions, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, St. Gallen 9014, Switzerland; Nanoparticle Systems Engineering Laboratory, Department of Mechanical and Process Engineering, Institute of Energy and Process Engineering, ETH Zurich, Sonneggstrasse 3, Zurich 8092, Switzerland
| | - Pascal Corso
- ARTORG Center for Biomedical Engineering Research, University of Bern, Freiburgstrasse 3, Bern 3010, Switzerland
| | - Robert Zboray
- Center for X-Ray Analytics, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Ueberlandstrasse 129, Duebendorf 8600, Switzerland
| | - Jonathan Avaro
- Center for X-Ray Analytics, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Ueberlandstrasse 129, Duebendorf 8600, Switzerland
| | | | - Marianne Liebi
- Center for X-Ray Analytics, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Ueberlandstrasse 129, Duebendorf 8600, Switzerland; Paul Scherrer Institute, PSI, Villigen 5232, Switzerland; Department of Physics, Chalmers University of Technology, Gothenburg 41296, Sweden
| | - Sergio Bertazzo
- Department of Medical Physics and Biomedical Engineering, University College London, WC1E 6BT, UK; London Centre for Nanotechnology, University College London, WC1E 6BT, UK
| | - Paul Philipp Heinisch
- Department of Cardiovascular Surgery, Inselspital, University of Bern, Freiburgstrasse 18, Bern 3010, Switzerland; Department of Congenital and Pediatric Heart Surgery, German Heart Center Munich, Technische Universität München, Germany
| | - Thierry Carrel
- Department of Cardiovascular Surgery, Inselspital, University of Bern, Freiburgstrasse 18, Bern 3010, Switzerland; Department of Cardiac Surgery, University Hospital Zurich (USZ), Rämistrasse 101, Zürich 8091, Switzerland.
| | - Dominik Obrist
- ARTORG Center for Biomedical Engineering Research, University of Bern, Freiburgstrasse 3, Bern 3010, Switzerland.
| | - Inge K Herrmann
- Laboratory for Particles-Biology Interactions, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, St. Gallen 9014, Switzerland; Nanoparticle Systems Engineering Laboratory, Department of Mechanical and Process Engineering, Institute of Energy and Process Engineering, ETH Zurich, Sonneggstrasse 3, Zurich 8092, Switzerland.
| |
Collapse
|
8
|
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
|
9
|
Zalud NC, Bulusu KV, Plesniak MW. Shear stress metrics associated with pro-atherogenic high-risk anatomical features in a carotid artery bifurcation model. Clin Biomech (Bristol, Avon) 2023; 105:105956. [PMID: 37098301 DOI: 10.1016/j.clinbiomech.2023.105956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/13/2022] [Revised: 04/06/2023] [Accepted: 04/13/2023] [Indexed: 04/27/2023]
Abstract
BACKGROUND Diseases associated with atherosclerotic plaques in the carotid artery are a major cause of deaths in the United States. Blood-flow-induced shear-stresses are known to trigger plaque formation. Prior literature suggests that the internal carotid artery sinus is prone to atherosclerosis, but there is limited understanding of why only certain patients are predisposed towards plaque formation. METHODS We computationally investigate the effect of vessel geometry on wall-shear-stress distribution by comparing flowfields and wall-shear-stress-metrics between a low-risk and a novel predisposed high-risk carotid artery bifurcation anatomy. Both models were developed based on clinical risk estimations and patient-averaged anatomical features. The high-risk geometry has a larger internal carotid artery branching angle and a lower internal-to-carotid-artery-diameter-ratio. A patient-averaged physiological carotid artery inflow waveform is used. FINDINGS The high-risk geometry experiences stronger flow separation in the sinus. Furthermore, it experiences a more equal flow split at the bifurcation, thereby reducing internal carotid artery flowrate and increasing atherosclerosis-prone low-velocity areas. Lowest time-averaged-wall-shear-stresses are present at the sinus outer wall, where plaques are often found, for both geometries. The high-risk geometry has significantly high, unfavorable oscillatory-shear-index values not found in the low-risk geometry. High oscillatory-shear-index areas are located at the vessels outside walls distal to the bifurcation and on the sinus wall. INTERPRETATION These results highlight the effectiveness of oscillatory-shear-index, to augment classical time-averaged-wall-shear-stress, in evaluating pro-atherogenic geometry features. Furthermore, the flow split at the bifurcation is a promising clinical indicator for atherosclerosis risk as it can be directly accessed using clinical imaging, whereas shear-stress-metrics cannot.
Collapse
Affiliation(s)
- Nora C Zalud
- Department of Mechanical and Aerospace Engineering, The George Washington University, 800 22nd Street NW, Science & Engineering Hall, Suite 3000, Washington, DC 20052, United States
| | - Kartik V Bulusu
- Department of Mechanical and Aerospace Engineering, The George Washington University, 800 22nd Street NW, Science & Engineering Hall, Suite 3000, Washington, DC 20052, United States
| | - Michael W Plesniak
- Department of Mechanical and Aerospace Engineering, The George Washington University, 800 22nd Street NW, Science & Engineering Hall, Suite 3000, Washington, DC 20052, United States; Department of Biomedical Engineering, The George Washington University, 800 22nd Street NW, Science & Engineering Hall, Suite 5000, Washington, DC 20052, United States.
| |
Collapse
|
10
|
Sache A, Reymond P, Brina O, Jung B, Farhat M, Vargas MI. Near-wall hemodynamic parameters quantification in in vitro intracranial aneurysms with 7 T PC-MRI. MAGMA (NEW YORK, N.Y.) 2023; 36:295-308. [PMID: 37072539 PMCID: PMC10140017 DOI: 10.1007/s10334-023-01082-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2022] [Revised: 03/20/2023] [Accepted: 03/21/2023] [Indexed: 04/20/2023]
Abstract
OBJECTIVE Wall shear stress (WSS) and its derived spatiotemporal parameters have proven to play a major role on intracranial aneurysms (IAs) growth and rupture. This study aims to demonstrate how ultra-high field (UHF) 7 T phase contrast magnetic resonance imaging (PC-MRI) coupled with advanced image acceleration techniques allows a highly resolved visualization of near-wall hemodynamic parameters patterns in in vitro IAs, paving the way for more robust risk assessment of their growth and rupture. MATERIALS AND METHODS We performed pulsatile flow measurements inside three in vitro models of patient-specific IAs using 7 T PC-MRI. To this end, we built an MRI-compatible test bench, which faithfully reproduced a typical physiological intracranial flow rate in the models. RESULTS The ultra-high field 7 T images revealed WSS patterns with high spatiotemporal resolution. Interestingly, the high oscillatory shear index values were found in the core of low WSS vortical structures and in flow stream intersecting regions. In contrast, maxima of WSS occurred around the impinging jet sites. CONCLUSIONS We showed that the elevated signal-to-noise ratio arising from 7 T PC-MRI enabled to resolve high and low WSS patterns with a high degree of detail.
Collapse
Affiliation(s)
- Antoine Sache
- Department of Mechanical Engineering, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.
| | - Philippe Reymond
- Division of Neuroradiology, Geneva University Hospital, University of Geneva, Geneva, Switzerland
| | - Olivier Brina
- Division of Neuroradiology, Geneva University Hospital, University of Geneva, Geneva, Switzerland
| | - Bernd Jung
- Department of Diagnostic, Interventional and Paediatric Radiology, Inselspital, Bern University Hospital, University of Bern, Bern, Switzerland
| | - Mohamed Farhat
- Department of Mechanical Engineering, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Maria Isabel Vargas
- Division of Neuroradiology, Geneva University Hospital, University of Geneva, Geneva, Switzerland
| |
Collapse
|
11
|
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
|
12
|
Wen J, Huang H, Su Z, Jiang L, Gao Q, Chen X, Yan T, Peng L. Predicting the Risk of Type B Aortic Dissection Using Hemodynamic Parameters in Aortic Arches: A Comparative Study between Healthy and Repaired Aortas. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 230:107326. [PMID: 36608431 DOI: 10.1016/j.cmpb.2022.107326] [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: 09/03/2022] [Revised: 11/19/2022] [Accepted: 12/24/2022] [Indexed: 06/17/2023]
Abstract
BACKGROUND AND OBJECTIVE The development of acute aortic dissection (AD) remains unpredictable due to the intricate nature of the AD mechanism and the varied patient-specific aortic anatomy. The aim of this study was to simulate the hemodynamic parameters in the aortas before the onset of TBAD with healthy controls. METHODS This study numerically assessed the effectiveness of hemodynamic indicators in predicting the risk of type B AD (TBAD) by investigating the differences in hemodynamic parameters between healthy and repaired aortas (aortas before TBAD development). Four wall shear stress (WSS)-based indicators and three helicity-based indicators were adopted and analyzed. RESULTS The results showed that more pathological anatomical feathers can be observed in the repaired aortas. For WSS-based indicators, only averaged cross flow index (CFI) and oscillatory shear index OSI (CFI, 1.03 ± 0.07 vs. 0.83 ± 0.10 and OSI, 0.12 ± 0.03 vs. 0.04 ± 0.02) (all p<0.001) were significantly higher in the repaired aortas than those in the healthy aortas. On the other hand, average helicity in the repaired aortas also showed a significant difference compared with that in healthy aortas (h1, 3.88 ± 5.55 vs. -8.03 ± 14.16) (p<0.05). Furthermore, the skewed helical structure and flow disturbance was found in the repaired aortas. CONCLUSION 1) There are marked differences in pathological anatomical features, such as aortic dilation, elongation and tortuosity between the healthy aortas and repaired aortas, and the corresponding hemodynamic indicators also have also been significantly changed. 2) Compared with anatomical characteristics, hemodynamic indicators may be more accurate for predicting the risk and location of TBAD, such as the OSI and CFI index were significantly enhanced in the region where the entry tears have occurred. 3) In clinical practice, anatomical features remain important factors for assessing the risk for development of TBAD; however, hemodynamic analyses with quantitative data and more visualizing characteristics have showed promising potential in this aspect.
Collapse
Affiliation(s)
- Jun Wen
- Department of Computer Science and technology, Southwest University of Science and Technology, Mianyang 621010, China
| | - Haodi Huang
- Institute of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
| | - Zhiqiao Su
- Institute of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
| | - Linke Jiang
- Institute of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
| | - Qi Gao
- Institute of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China
| | - Xiaoyi Chen
- Department of Radiology, West China Hospital, Sichuan University, Chengdu 610041, China
| | - Tingli Yan
- Department of Radiology, West China Hospital, Sichuan University, Chengdu 610041, China
| | - Liqing Peng
- Department of Radiology, West China Hospital, Sichuan University, Chengdu 610041, China.
| |
Collapse
|
13
|
Mazzi V, De Nisco G, Calò K, Chiastra C, Daemen J, Steinman DA, Wentzel JJ, Morbiducci U, Gallo D. Divergence of the normalized wall shear stress as an effective computational template of low-density lipoprotein polarization at the arterial blood-vessel wall interface. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2022; 226:107174. [PMID: 36223707 DOI: 10.1016/j.cmpb.2022.107174] [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: 08/17/2022] [Revised: 09/20/2022] [Accepted: 10/03/2022] [Indexed: 06/16/2023]
Abstract
BACKGROUND AND OBJECTIVE Near-wall transport of low-density lipoproteins (LDL) in arteries plays a relevant role in the initiation of atherosclerosis. Although it can be modelled in silico by coupling the Navier-Stokes equations with the 3D advection-diffusion (AD) equation, the associated computational cost is high. As wall shear stress (WSS) represents a first-order approximation of the near-wall velocity in arteries, we aimed at identifying computationally convenient WSS-based quantities to infer LDL near-wall transport based on the underlying near-wall hemodynamics in five models of three human arterial districts (aorta, carotid bifurcations, coronary arteries). The simulated LDL transport and its WSS-based surrogates were qualitatively compared with in vivo longitudinal measurements of wall thickness growth on the coronary artery models. METHODS Numerical simulations of blood flow coupled with AD equations for LDL transport and blood-wall transfer were performed. The co-localization of the simulated LDL concentration polarization patterns with luminal surface areas characterized by low cycle-average WSS, near-wall flow stagnation and WSS attracting patterns was quantitatively assessed by the similarity index (SI). In detail, the latter two represent features of the WSS topological skeleton, obtained respectively through the Lagrangian tracking of surface-born particles, and the Eulerian analysis of the divergence of the normalized cycle-average WSS vector field. RESULTS Convergence of the solution of the AD problem required the simulation of 3 (coronary artery) to 10 (aorta) additional cardiac cycles with respect to the Navier-Stokes problem. Co-localization results underlined that WSS topological skeleton features indicating near-wall flow stagnation and WSS attracting patterns identified LDL concentration polarization profiles more effectively than low WSS, as indicated by higher SI values (SI range: 0.17-0.50 for low WSS; 0.24-0.57 for WSS topological skeleton features). Moreover, the correspondence between the simulated LDL uptake and WSS-based quantities profiles with the in vivo measured wall thickness growth in coronary arteries appears promising. CONCLUSIONS The recently introduced Eulerian approach for identifying WSS attracting patterns from the divergence of normalized WSS provides a computationally affordable template of the LDL polarization at the arterial blood-wall interface without simulating the AD problem. It thus candidates as an effective biomechanical tool for elucidating the mechanistic link amongst LDL transfer at the arterial blood-wall interface, WSS and atherogenesis.
Collapse
Affiliation(s)
- Valentina Mazzi
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Giuseppe De Nisco
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Karol Calò
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Claudio Chiastra
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Joost Daemen
- Department of Cardiology, Biomedical Engineering, Erasmus MC, 3000 CA Rotterdam, the Netherlands
| | - David A Steinman
- Biomedical Simulation Lab, Department of Mechanical and Industrial Engineering, University of Toronto, Toronto, Ontario, Canada
| | - Jolanda J Wentzel
- Department of Cardiology, Biomedical Engineering, Erasmus MC, 3000 CA Rotterdam, the Netherlands
| | - Umberto Morbiducci
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Diego Gallo
- PoliTo(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy.
| |
Collapse
|
14
|
Atherogenic potential of microgravity hemodynamics in the carotid bifurcation: a numerical investigation. NPJ Microgravity 2022; 8:39. [PMID: 36085153 PMCID: PMC9463447 DOI: 10.1038/s41526-022-00223-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Accepted: 08/22/2022] [Indexed: 11/08/2022] Open
Abstract
Long-duration spaceflight poses multiple hazards to human health, including physiological changes associated with microgravity. The hemodynamic adaptations occurring upon entry into weightlessness have been associated with retrograde stagnant flow conditions and thromboembolic events in the venous vasculature but the impact of microgravity on cerebral arterial hemodynamics and function remains poorly understood. The objective of this study was to quantify the effects of microgravity on hemodynamics and wall shear stress (WSS) characteristics in 16 carotid bifurcation geometries reconstructed from ultrasonography images using computational fluid dynamics modeling. Microgravity resulted in a significant 21% increase in flow stasis index, a 22-23% decrease in WSS magnitude and a 16-26% increase in relative residence time in all bifurcation branches, while preserving WSS unidirectionality. In two anatomies, however, microgravity not only promoted flow stasis but also subjected the convex region of the external carotid arterial wall to a moderate increase in WSS bidirectionality, which contrasted with the population average trend. This study suggests that long-term exposure to microgravity has the potential to subject the vasculature to atheroprone hemodynamics and this effect is modulated by subject-specific anatomical features. The exploration of the biological impact of those microgravity-induced WSS aberrations is needed to better define the risk posed by long spaceflights on cardiovascular health.
Collapse
|
15
|
Shar JA, Keswani SG, Grande-Allen KJ, Sucosky P. Significance of aortoseptal angle anomalies to left ventricular hemodynamics and subaortic stenosis: A numerical study. Comput Biol Med 2022; 146:105613. [PMID: 35751200 PMCID: PMC10570849 DOI: 10.1016/j.compbiomed.2022.105613] [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: 04/05/2022] [Revised: 05/02/2022] [Accepted: 05/10/2022] [Indexed: 11/03/2022]
Abstract
PURPOSE Discrete subaortic stenosis (DSS) is an obstructive cardiac disease caused by a membranous lesion in the left ventricular (LV) outflow tract (LVOT). Although its etiology is unknown, the higher prevalence of DSS in LVOT anatomies featuring a steep aortoseptal angle (AoSA) suggests a potential role for hemodynamics. Therefore, the objective of this study was to quantify the impact of AoSA steepening on the LV three-dimensional (3D) hemodynamic stress environment. METHODS A 3D LV model reconstructed from cardiac cine-magnetic resonance imaging was connected to four LVOT geometrical variations spanning the clinical AoSA range (115°-160°). LV hemodynamic stresses were characterized in terms of cycle-averaged pressure, temporal shear magnitude (TSM), and oscillatory shear index. The wall shear stress (WSS) topological skeleton was further analyzed by computing the scaled divergence of the WSS vector field. RESULTS AoSA steepening caused an increasingly perturbed subaortic flow marked by LVOT flow skewness and complex 3D secondary flow patterns. These disturbances generated WSS overloads (>45% increase in TSM vs. 160° model) on the inferior LVOT wall, and increased WSS contraction (>66% decrease in WSS divergence vs. 160° model) in regions prone to DSS membrane formation. CONCLUSIONS AoSA steepening generated substantial hemodynamic stress abnormalities in LVOT regions prone to DSS formation. Further studies are needed to assess the possible impact of such mechanical abnormalities on the tissue and cellular responses.
Collapse
Affiliation(s)
- Jason A Shar
- Department of Mechanical Engineering, Kennesaw State University, 840 Polytechnic Lane, Marietta, GA, 30060, USA.
| | - Sundeep G Keswani
- Division of Pediatric Surgery, Texas Children's Hospital, Department of Surgery, Baylor College of Medicine, USA.
| | | | - Philippe Sucosky
- Department of Mechanical Engineering, Kennesaw State University, 840 Polytechnic Lane, Marietta, GA, 30060, USA.
| |
Collapse
|
16
|
He Y, Northrup H, Le H, Cheung AK, Berceli SA, Shiu YT. Medical Image-Based Computational Fluid Dynamics and Fluid-Structure Interaction Analysis in Vascular Diseases. Front Bioeng Biotechnol 2022; 10:855791. [PMID: 35573253 PMCID: PMC9091352 DOI: 10.3389/fbioe.2022.855791] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2022] [Accepted: 04/08/2022] [Indexed: 01/17/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Yong He
- Division of Vascular Surgery and Endovascular Therapy, University of Florida, Gainesville, FL, United States
| | - Hannah Northrup
- Department of Biomedical Engineering, University of Utah, Salt Lake City, UT, United States
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
| | - Ha Le
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
| | - Alfred K. Cheung
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
- Veterans Affairs Salt Lake City Healthcare System, Salt Lake City, UT, United States
| | - Scott A. Berceli
- Division of Vascular Surgery and Endovascular Therapy, University of Florida, Gainesville, FL, United States
- Vascular Surgery Section, Malcom Randall Veterans Affairs Medical Center, Gainesville, FL, United States
| | - Yan Tin Shiu
- Division of Nephrology and Hypertension, Department of Internal Medicine, University of Utah, Salt Lake City, UT, United States
- Veterans Affairs Salt Lake City Healthcare System, Salt Lake City, UT, United States
- *Correspondence: Yan Tin Shiu,
| |
Collapse
|
17
|
Franco P, Sotelo J, Guala A, Dux-Santoy L, Evangelista A, Rodríguez-Palomares J, Mery D, Salas R, Uribe S. Identification of hemodynamic biomarkers for bicuspid aortic valve induced aortic dilation using machine learning. Comput Biol Med 2021; 141:105147. [PMID: 34929463 DOI: 10.1016/j.compbiomed.2021.105147] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 12/13/2021] [Accepted: 12/13/2021] [Indexed: 01/06/2023]
Abstract
Recent advances in medical imaging have confirmed the presence of altered hemodynamics in bicuspid aortic valve (BAV) patients. Therefore, there is a need for new hemodynamic biomarkers to refine disease monitoring and improve patient risk stratification. This research aims to analyze and extract multiple correlation patterns of hemodynamic parameters from 4D Flow MRI data and find which parameters allow an accurate classification between healthy volunteers (HV) and BAV patients with dilated and non-dilated ascending aorta using machine learning. Sixteen hemodynamic parameters were calculated in the ascending aorta (AAo) and aortic arch (AArch) at peak systole from 4D Flow MRI. We used sequential forward selection (SFS) and principal component analysis (PCA) as feature selection algorithms. Then, eleven machine-learning classifiers were implemented to separate HV and BAV patients (non- and dilated ascending aorta). Multiple correlation patterns from hemodynamic parameters were extracted using hierarchical clustering. The linear discriminant analysis and random forest are the best performing classifiers, using five hemodynamic parameters selected with SFS (velocity angle, forward velocity, vorticity, and backward velocity in AAo; and helicity density in AArch) a 96.31 ± 1.76% and 96.00 ± 0.83% accuracy, respectively. Hierarchical clustering revealed three groups of correlated features. According to this analysis, we observed that features selected by SFS have a better performance than those selected by PCA because the five selected parameters were distributed according to 3 different clusters. Based on the proposed method, we concluded that the feature selection method found five potentially hemodynamic biomarkers related to this disease.
Collapse
Affiliation(s)
- Pamela Franco
- Biomedical Imaging Center, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile; Electrical Engineering Department, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile; Millennium Nucleus for Cardiovascular Magnetic Resonance, Cardio, MR, Chile
| | - Julio Sotelo
- Biomedical Imaging Center, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile; Millennium Nucleus for Cardiovascular Magnetic Resonance, Cardio, MR, Chile; School of Biomedical Engineering, Universidad de Valparaíso, Valparaíso, Chile
| | - Andrea Guala
- Department of Cardiology, Hospital Universitari Vall d'Hebron, Vall d'Hebron Institut de Recerca (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Lydia Dux-Santoy
- Department of Cardiology, Hospital Universitari Vall d'Hebron, Vall d'Hebron Institut de Recerca (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Arturo Evangelista
- Department of Cardiology, Hospital Universitari Vall d'Hebron, Vall d'Hebron Institut de Recerca (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - José Rodríguez-Palomares
- Department of Cardiology, Hospital Universitari Vall d'Hebron, Vall d'Hebron Institut de Recerca (VHIR), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Domingo Mery
- Department of Computer Science, Pontificia Universidad Católica de Chile, Santiago, Chile; Instituto Milenio Intelligent Healthcare Engineering, Chile
| | - Rodrigo Salas
- School of Biomedical Engineering, Universidad de Valparaíso, Valparaíso, Chile; Instituto Milenio Intelligent Healthcare Engineering, Chile
| | - Sergio Uribe
- Biomedical Imaging Center, School of Engineering, Pontificia Universidad Católica de Chile, Santiago, Chile; Millennium Nucleus for Cardiovascular Magnetic Resonance, Cardio, MR, Chile; Radiology Department, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile; Instituto Milenio Intelligent Healthcare Engineering, Chile.
| |
Collapse
|
18
|
Pillalamarri NR, Piskin S, Patnaik SS, Murali S, Finol EA. Patient-Specific Computational Analysis of Hemodynamics in Adult Pulmonary Hypertension. Ann Biomed Eng 2021; 49:3465-3480. [PMID: 34799807 DOI: 10.1007/s10439-021-02884-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2021] [Accepted: 10/29/2021] [Indexed: 11/25/2022]
Abstract
Pulmonary hypertension (PH) is a progressive disease characterized by elevated pressure and vascular resistance in the pulmonary arteries. Nearly 250,000 hospitalizations occur annually in the US with PH as the primary or secondary condition. A definitive diagnosis of PH requires right heart catheterization (RHC) in addition to a chest computed tomography, a walking test, and others. While RHC is the gold standard for diagnosing PH, it is invasive and posseses inherent risks and contraindications. In this work, we characterized the patient-specific pulmonary hemodynamics in silico for diverse PH WHO groups. We grouped patients on the basis of mean pulmonary arterial pressure (mPAP) into three disease severity groups: at-risk ([Formula: see text], denoted with A), mild ([Formula: see text], denoted with M), and severe ([Formula: see text], denoted with S). The pulsatile flow hemodynamics was simulated by evaluating the three-dimensional Navier-Stokes system of equations using a flow solver developed by customizing OpenFOAM libraries (v5.0, The OpenFOAM Foundation). Quasi patient-specific boundary conditions were implemented using a Womersley inlet velocity profile and transient resistance outflow conditions. Hemodynamic indices such as spatially averaged wall shear stress ([Formula: see text]), wall shear stress gradient ([Formula: see text]), time-averaged wall shear stress ([Formula: see text]), oscillatory shear index ([Formula: see text]), and relative residence time ([Formula: see text]), were evaluated along with the clinical metrics pulmonary vascular resistance ([Formula: see text]), stroke volume ([Formula: see text]) and compliance ([Formula: see text]), to assess possible spatiotemporal correlations. We observed statistically significant decreases in [Formula: see text], [Formula: see text], and [Formula: see text], and increases in [Formula: see text] and [Formula: see text] with disease severity. [Formula: see text] was moderately correlated with [Formula: see text] and [Formula: see text] at the mid-notch stage of the cardiac cycle when these indices were computed using the global pulmonary arterial geometry. These results are promising in the context of a long-term goal of identifying computational biomarkers that can serve as surrogates for invasive diagnostic protocols of PH.
Collapse
Affiliation(s)
- Narasimha R Pillalamarri
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
| | - Senol Piskin
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
- Department of Mechanical Engineering, Istinye University, Istanbul, Turkey
| | - Sourav S Patnaik
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA
- Department of Bioengineering, University of Texas at Dallas, Richardson, TX, USA
| | - Srinivas Murali
- Division of Cardiovascular Medicine, Allegheny Health Network, Pittsburgh, PA, USA
| | - Ender A Finol
- Department of Mechanical Engineering, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX, 78249, USA.
- UTSA/UTHSA Joint Graduate Program in Biomedical Engineering, University of Texas at San Antonio, San Antonio, TX, USA.
| |
Collapse
|
19
|
Candreva A, Pagnoni M, Rizzini ML, Mizukami T, Gallinoro E, Mazzi V, Gallo D, Meier D, Shinke T, Aben JP, Nagumo S, Sonck J, Munhoz D, Fournier S, Barbato E, Heggermont W, Cook S, Chiastra C, Morbiducci U, De Bruyne B, Muller O, Collet C. Risk of myocardial infarction based on endothelial shear stress analysis using coronary angiography. Atherosclerosis 2021; 342:28-35. [PMID: 34815069 DOI: 10.1016/j.atherosclerosis.2021.11.010] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 08/30/2021] [Revised: 10/25/2021] [Accepted: 11/09/2021] [Indexed: 11/18/2022]
Abstract
BACKGROUND AND AIMS Wall shear stress (WSS) has been associated with atherogenesis and plaque progression. The present study assessed the value of WSS analysis derived from conventional coronary angiography to detect lesions culprit for future myocardial infarction (MI). METHODS AND RESULTS Three-dimensional quantitative coronary angiography (3DQCA), was used to calculate WSS and pressure drop in 80 patients. WSS descriptors were compared between 80 lesions culprit of future MI and 108 non-culprit lesions (controls). Endothelium-blood flow interaction was assessed by computational fluid dynamics (10.8 ± 1.41 min per vessel). Median time between baseline angiography and MI was 25.9 (21.9-29.8) months. Mean patient age was 70.3 ± 12.7. Clinical presentation was STEMI in 35% and NSTEMI in 65%. Culprit lesions showed higher percent area stenosis (%AS), translesional vFFR difference (ΔvFFR), time-averaged WSS (TAWSS) and topological shear variation index (TSVI) compared to non-culprit lesions (p < 0.05 for all). TSVI was superior to TAWSS in predicting MI (AUC-TSVI = 0.77, 95%CI 0.71-0.84 vs. AUC-TAWSS = 0.61, 95%CI 0.53-0.69, p < 0.001). The addition of TSVI increased predictive and reclassification abilities compared to a model based on %AS and ΔvFFR (NRI = 1.04, p < 0.001, IDI = 0.22, p < 0.001). CONCLUSIONS A 3DQCA-based WSS analysis was feasible and can identify lesions culprit for future MI. The combination of area stenoses, pressure gradients and WSS predicted the occurrence of MI. TSVI, a novel WSS descriptor, showed strong predictive capacity to detect lesions prone to cause MI.
Collapse
Affiliation(s)
- Alessandro Candreva
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Cardiology, Zurich University Hospital, Zurich, Switzerland; Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Mattia Pagnoni
- Dept. of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Maurizio Lodi Rizzini
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Takuya Mizukami
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Cardiology, Showa University School of Medicine, Tokyo, Japan
| | - Emanuele Gallinoro
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Translational Medical Sciences, University of Campania 'Luigi Vanvitelli', Naples, Italy
| | - Valentina Mazzi
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - Diego Gallo
- Polito(BIO)Med Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| | - David Meier
- Dept. of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Toshiro Shinke
- Dept. of Cardiology, Showa University School of Medicine, Tokyo, Japan
| | | | - Sakura Nagumo
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Cardiology, Showa University School of Medicine, Tokyo, Japan
| | - Jeroen Sonck
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Daniel Munhoz
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy; Department of Internal Medicine, University of Campinas (Unicamp), Campinas, Brazil
| | - Stephane Fournier
- Dept. of Cardiology, Lausanne University Hospital, Lausanne, Switzerland; Dept. of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | - Emanuele Barbato
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Advanced Biomedical Sciences, University of Naples Federico II, Naples, Italy
| | | | - Stephane Cook
- Department of Cardiology, HFR Fribourg, Fribourg, 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
| | - Bernard De Bruyne
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium; Dept. of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Oliver Muller
- Dept. of Cardiology, Lausanne University Hospital, Lausanne, Switzerland
| | - Carlos Collet
- Cardiovascular Center Aalst, OLV-Clinic, Aalst, Belgium.
| |
Collapse
|
20
|
Stokes C, Bonfanti M, Li Z, Xiong J, Chen D, Balabani S, Díaz-Zuccarini V. A novel MRI-based data fusion methodology for efficient, personalised, compliant simulations of aortic haemodynamics. J Biomech 2021; 129:110793. [PMID: 34715606 PMCID: PMC8907869 DOI: 10.1016/j.jbiomech.2021.110793] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 08/24/2021] [Accepted: 09/30/2021] [Indexed: 01/24/2023]
Abstract
We present a novel, cost-efficient methodology to simulate aortic haemodynamics in a patient-specific, compliant aorta using an MRI data fusion process. Based on a previously-developed Moving Boundary Method, this technique circumvents the high computational cost and numerous structural modelling assumptions required by traditional Fluid-Structure Interaction techniques. Without the need for Computed Tomography (CT) data, the MRI images required to construct the simulation can be obtained during a single imaging session. Black Blood MR Angiography and 2D Cine-MRI data were used to reconstruct the luminal geometry and calibrate wall movement specifically to each region of the aorta. 4D-Flow MRI and non-invasive pressure measurements informed patient-specific inlet and outlet boundary conditions. Luminal area closely matched 2D Cine-MRI measurements with a mean error of less than 4.6% across the cardiac cycle, while physiological pressure and flow distributions were simulated to within 3.3% of patient-specific targets. Moderate agreement with 4D-Flow MRI velocity data was observed. Despite lower peak velocity, an equivalent rigid-wall simulation predicted a mean Time-Averaged Wall Shear Stress (TAWSS) 13% higher than the compliant simulation. The agreement observed between compliant simulation results and MRI data is testament to the accuracy and efficiency of this MRI-based simulation technique.
Collapse
Affiliation(s)
- Catriona Stokes
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Mirko Bonfanti
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Zeyan Li
- School of Life Science, Beijing Institute of Technology, Beijing, China.
| | - Jiang Xiong
- Department of Vascular and Endovascular Surgery, Chinese PLA General Hospital, Beijing, China.
| | - Duanduan Chen
- School of Life Science, Beijing Institute of Technology, Beijing, China.
| | - Stavroula Balabani
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| | - Vanessa Díaz-Zuccarini
- Mechanical Engineering Department, Roberts Engineering Building, University College London, Torrington Place, London, WC1E 7JE, United Kingdom; Wellcome/EPSRC Centre for Interventional and Surgical Sciences (WEISS), Charles Bell House, London, W1W 7TY, United Kingdom.
| |
Collapse
|
21
|
Mazzi V, De Nisco G, Hoogendoorn A, Calò K, Chiastra C, Gallo D, Steinman DA, Wentzel JJ, Morbiducci U. Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability. Ann Biomed Eng 2021; 49:2606-2621. [PMID: 34324092 PMCID: PMC8455396 DOI: 10.1007/s10439-021-02829-5] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.
Collapse
Affiliation(s)
- Valentina Mazzi
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Giuseppe De Nisco
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Ayla Hoogendoorn
- Department of Cardiology, Biomedical Engineering, Erasmus MC, 3000 CA, Rotterdam, The Netherlands
| | - Karol Calò
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Claudio Chiastra
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Diego Gallo
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - David A Steinman
- Biomedical Simulation Laboratory, Department of Mechanical & Industrial Engineering, University of Toronto, Toronto, Canada
| | - Jolanda J Wentzel
- Department of Cardiology, Biomedical Engineering, Erasmus MC, 3000 CA, Rotterdam, The Netherlands
| | - Umberto Morbiducci
- PoliToBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy.
| |
Collapse
|
22
|
Mahmoudi M, Farghadan A, McConnell DR, Barker AJ, Wentzel JJ, Budoff MJ, Arzani A. The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction. J Biomech Eng 2021; 143:041002. [PMID: 33156343 DOI: 10.1115/1.4049026] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Indexed: 12/20/2022]
Abstract
Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.
Collapse
Affiliation(s)
- Mostafa Mahmoudi
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ 86011
| | - Ali Farghadan
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ 86011
| | - Daniel R McConnell
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ 86011
| | - Alex J Barker
- Department of Pediatric Radiology, University of Colorado Denver Anschutz Medical Campus, Aurora, CO 80045
| | - Jolanda J Wentzel
- Department of Cardiology, Biomedical Engineering, Erasmus MC, Rotterdam, The Netherlands
| | | | - Amirhossein Arzani
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ 86011
| |
Collapse
|
23
|
Wall Shear Stress Topological Skeleton Analysis in Cardiovascular Flows: Methods and Applications. MATHEMATICS 2021. [DOI: 10.3390/math9070720] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A marked interest has recently emerged regarding the analysis of the wall shear stress (WSS) vector field topological skeleton in cardiovascular flows. Based on dynamical system theory, the WSS topological skeleton is composed of fixed points, i.e., focal points where WSS locally vanishes, and unstable/stable manifolds, consisting of contraction/expansion regions linking fixed points. Such an interest arises from its ability to reflect the presence of near-wall hemodynamic features associated with the onset and progression of vascular diseases. Over the years, Lagrangian-based and Eulerian-based post-processing techniques have been proposed aiming at identifying the topological skeleton features of the WSS. Here, the theoretical and methodological bases supporting the Lagrangian- and Eulerian-based methods currently used in the literature are reported and discussed, highlighting their application to cardiovascular flows. The final aim is to promote the use of WSS topological skeleton analysis in hemodynamic applications and to encourage its application in future mechanobiology studies in order to increase the chance of elucidating the mechanistic links between blood flow disturbances, vascular disease, and clinical observations.
Collapse
|
24
|
Ahmadpour-B M, Nooraeen A, Tafazzoli-Shadpour M, Taghizadeh H. Contribution of atherosclerotic plaque location and severity to the near-wall hemodynamics of the carotid bifurcation: an experimental study and FSI modeling. Biomech Model Mechanobiol 2021; 20:1069-1085. [PMID: 33609192 DOI: 10.1007/s10237-021-01431-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 02/03/2021] [Indexed: 01/22/2023]
Abstract
Atherosclerosis is initiated by endothelial injury that is related to abnormal values of hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI) and stress phase angle (SPA), which are more common in arterial bifurcations due to the complex structure. An experimental model of human carotid bifurcation with accurate geometrical and mechanical features was set up, and using realistic pulsatile flow rates, the inlet and outlet pressure pulses were measured for normal and stenosed models with 40% and 80% severities at common carotid (CCA), internal carotid (ICA) and external carotid (ECA) arteries. Based on the obtained experimental data, fluid-structure models were developed to obtain WSS, OSI, and SPA and evaluate pathological consequences at different locations. Mild severity had minor impact, however, inducing severe 80% stenosis in each branch led to considerable localized changes of hemodynamic parameters both in the stenosis site and other locations. This included sharp increases in WSS values accompanied by very low values close to zero before and after the peaks. Severe stenosis not only caused significant changes in the local artery, but also in other branches. OSI and SPA were less sensitive to stenosis, although high peaks were observed on bifurcation site for the stenosis at ECA. The interconnection of arteries at carotid bifurcation results in altered pressure/flow patterns in all branches when a stenosis is applied in any site. Such effect confirms pathological findings that atherosclerotic plaques are observed simultaneously in different carotid branches, although with different degrees of plaque growth and severity.
Collapse
Affiliation(s)
- Mahyar Ahmadpour-B
- Cardiovascular Engineering Lab, Faculty of Biomedical Engineering, Amirkabir University of Technology, Tehran, Iran
| | - Ahmad Nooraeen
- Tissue Mechanics Lab, Faculty of Biomedical Engineering, Sahand University of Technology, Tabriz, Iran
| | - Mohammad Tafazzoli-Shadpour
- Cardiovascular Engineering Lab, Faculty of Biomedical Engineering, Amirkabir University of Technology, Tehran, Iran.
| | - Hadi Taghizadeh
- Tissue Mechanics Lab, Faculty of Biomedical Engineering, Sahand University of Technology, Tabriz, Iran.
| |
Collapse
|
25
|
Morbiducci U, Mazzi V, Domanin M, De Nisco G, Vergara C, Steinman DA, Gallo D. Wall Shear Stress Topological Skeleton Independently Predicts Long-Term Restenosis After Carotid Bifurcation Endarterectomy. Ann Biomed Eng 2020; 48:2936-2949. [PMID: 32929560 PMCID: PMC7723943 DOI: 10.1007/s10439-020-02607-9] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Accepted: 09/02/2020] [Indexed: 12/26/2022]
Abstract
Wall Shear Stress (WSS) topological skeleton, composed by fixed points and the manifolds linking them, reflects the presence of blood flow features associated to adverse vascular response. However, the influence of WSS topological skeleton on vascular pathophysiology is still underexplored. This study aimed to identify direct associations between the WSS topological skeleton and markers of vascular disease from real-world clinical longitudinal data of long-term restenosis after carotid endarterectomy (CEA). Personalized computational hemodynamic simulations were performed on a cohort of 13 carotid models pre-CEA and at 1 month after CEA. At 60 months after CEA, intima-media thickness (IMT) was measured to detect long-term restenosis. The analysis of the WSS topological skeleton was carried out by applying a Eulerian method based on the WSS vector field divergence. To provide objective thresholds for WSS topological skeleton quantitative analysis, a computational hemodynamic dataset of 46 ostensibly healthy carotid bifurcation models was considered. CEA interventions did not completely restore physiological WSS topological skeleton features. Significant associations emerged between IMT at 60 months follow-up and the exposure to (1) high temporal variation of WSS contraction/expansion (R2 = 0.51, p < 0.05), and (2) high fixed point residence times, weighted by WSS contraction/expansion strength (R2 = 0.53, p < 0.05). These WSS topological skeleton features were statistically independent from the exposure to low WSS, a previously reported predictor of long-term restenosis, therefore representing different hemodynamic stimuli and potentially impacting differently the vascular response. This study confirms the direct association between WSS topological skeleton and markers of vascular disease, contributing to elucidate the mechanistic link between flow disturbances and clinical observations of vascular lesions.
Collapse
Affiliation(s)
- Umberto Morbiducci
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Valentina Mazzi
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Maurizio Domanin
- Department of Clinical Sciences and Community Health, Università di Milano, Milan, Italy.,Unità Operativa di Chirurgia Vascolare, Fondazione I.R.C.C.S. Cà Granda Ospedale Maggiore Policlinico, Milan, Italy
| | - Giuseppe De Nisco
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy
| | - Christian Vergara
- Laboratory of Biological Structure Mechanics (LaBS), Dipartimento di Chimica, Materiali e Ingegneria Chimica ''Giulio Natta'', Politecnico di Milano, Milan, Italy
| | - David A Steinman
- Biomedical Simulation Laboratory, Department of Mechanical & Industrial Engineering, University of Toronto, Toronto, ON, Canada
| | - Diego Gallo
- PolitoBIOMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi, 24, 10129, Turin, Italy.
| |
Collapse
|
26
|
Deciphering ascending thoracic aortic aneurysm hemodynamics in relation to biomechanical properties. Med Eng Phys 2020; 82:119-129. [DOI: 10.1016/j.medengphy.2020.07.003] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 06/19/2020] [Accepted: 07/09/2020] [Indexed: 12/20/2022]
|