1
|
Torbati S, Daneshmehr A, Pouraliakbar H, Asgharian M, Ahmadi Tafti SH, Shum-Tim D, Heidari A. Personalized evaluation of the passive myocardium in ischemic cardiomyopathy via computational modeling using Bayesian optimization. Biomech Model Mechanobiol 2024:10.1007/s10237-024-01856-0. [PMID: 38954283 DOI: 10.1007/s10237-024-01856-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2022] [Accepted: 04/28/2024] [Indexed: 07/04/2024]
Abstract
Biomechanics-based patient-specific modeling is a promising approach that has proved invaluable for its clinical potential to assess the adversities caused by ischemic heart disease (IHD). In the present study, we propose a framework to find the passive material properties of the myocardium and the unloaded shape of cardiac ventricles simultaneously in patients diagnosed with ischemic cardiomyopathy (ICM). This was achieved by minimizing the difference between the simulated and the target end-diastolic pressure-volume relationships (EDPVRs) using black-box Bayesian optimization, based on the finite element analysis (FEA). End-diastolic (ED) biventricular geometry and the location of the ischemia were determined from cardiac magnetic resonance (CMR) imaging. We employed our pipeline to model the cardiac ventricles of three patients aged between 57 and 66 years, with and without the inclusion of valves. An excellent agreement between the simulated and the target EDPVRs has been reached. Our results revealed that the incorporation of valvular springs typically leads to lower hyperelastic parameters for both healthy and ischemic myocardium, as well as a higher fiber Green strain in the viable regions compared to models without valvular stiffness. Furthermore, the addition of valve-related effects did not result in significant changes in myofiber stress after optimization. We concluded that more accurate results could be obtained when cardiac valves were considered in modeling ventricles. The present novel and practical methodology paves the way for developing digital twins of ischemic cardiac ventricles, providing a non-invasive assessment for designing optimal personalized therapies in precision medicine.
Collapse
Affiliation(s)
- Saeed Torbati
- Research Center for Advanced Technologies in Cardiovascular Medicine, Cardiovascular Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran
| | - Alireza Daneshmehr
- School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran
| | - Hamidreza Pouraliakbar
- Rajaie Cardiovascular, Medical, and Research Center, Iran University of Medical Sciences, Tehran, Iran
| | - Masoud Asgharian
- Department of Mathematics and Statistics, McGill University, Montreal, QC, Canada
| | - Seyed Hossein Ahmadi Tafti
- Research Center for Advanced Technologies in Cardiovascular Medicine, Cardiovascular Diseases Research Institute, Tehran University of Medical Sciences, Tehran, Iran.
- Department of Surgery, Tehran Heart Center, Tehran University of Medical Sciences, Tehran, Iran.
| | - Dominique Shum-Tim
- Division of Cardiac Surgery, Department of Surgery, McGill University, Montreal, QC, Canada
| | - Alireza Heidari
- Department of Mathematics and Statistics, McGill University, Montreal, QC, Canada.
- Department of Mechanical Engineering, McGill University, Montreal, QC, Canada.
- Department of Biomedical Engineering, McGill University, Montreal, QC, Canada.
| |
Collapse
|
2
|
Assila N, Begon M, Duprey S. Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion. Ann Biomed Eng 2024; 52:1240-1254. [PMID: 38376768 DOI: 10.1007/s10439-024-03449-5] [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/06/2023] [Accepted: 01/09/2024] [Indexed: 02/21/2024]
Abstract
The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.
Collapse
Affiliation(s)
- Najoua Assila
- School of Kinesiology and Exercise Sciences, Faculty of Medicine, University of Montréal, Montréal, QC, Canada.
- Research Center of the Sainte-Justine University Hospital Center, Montréal, QC, Canada.
- Univ Lyon, Univ Gustave Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR T_9406, 69622, Lyon, France.
| | - Mickaël Begon
- School of Kinesiology and Exercise Sciences, Faculty of Medicine, University of Montréal, Montréal, QC, Canada
- Research Center of the Sainte-Justine University Hospital Center, Montréal, QC, Canada
| | - Sonia Duprey
- Univ Lyon, Univ Gustave Eiffel, Univ Claude Bernard Lyon 1, LBMC UMR T_9406, 69622, Lyon, France
| |
Collapse
|
3
|
Ren M, Ong C, Buist ML, Yap CH. Biventricular biaxial mechanical testing and constitutive modelling of fetal porcine myocardium passive stiffness. J Mech Behav Biomed Mater 2022; 134:105383. [DOI: 10.1016/j.jmbbm.2022.105383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2022] [Revised: 07/02/2022] [Accepted: 07/15/2022] [Indexed: 11/15/2022]
|
4
|
Koopsen T, Van Osta N, Van Loon T, Van Nieuwenhoven FA, Prinzen FW, Van Klarenbosch BR, Kirkels FP, Teske AJ, Vernooy K, Delhaas T, Lumens J. A Lumped Two-Compartment Model for Simulation of Ventricular Pump and Tissue Mechanics in Ischemic Heart Disease. Front Physiol 2022; 13:782592. [PMID: 35634163 PMCID: PMC9130776 DOI: 10.3389/fphys.2022.782592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2021] [Accepted: 03/10/2022] [Indexed: 11/16/2022] Open
Abstract
Introduction: Computational modeling of cardiac mechanics and hemodynamics in ischemic heart disease (IHD) is important for a better understanding of the complex relations between ischemia-induced heterogeneity of myocardial tissue properties, regional tissue mechanics, and hemodynamic pump function. We validated and applied a lumped two-compartment modeling approach for IHD integrated into the CircAdapt model of the human heart and circulation. Methods: Ischemic contractile dysfunction was simulated by subdividing a left ventricular (LV) wall segment into a hypothetical contractile and noncontractile compartment, and dysfunction severity was determined by the noncontractile volume fraction (NCVF). Myocardial stiffness was determined by the zero-passive stress length (Ls0,pas) and nonlinearity (kECM) of the passive stress-sarcomere length relation of the noncontractile compartment. Simulated end-systolic pressure volume relations (ESPVRs) for 20% acute ischemia were qualitatively compared between a two- and one-compartment simulation, and parameters of the two-compartment model were tuned to previously published canine data of regional myocardial deformation during acute and prolonged ischemia and reperfusion. In six patients with myocardial infarction (MI), the NCVF was automatically estimated using the echocardiographic LV strain and volume measurements obtained acutely and 6 months after MI. Estimated segmental NCVF values at the baseline and 6-month follow-up were compared with percentage late gadolinium enhancement (LGE) at 6-month follow-up. Results: Simulation of 20% of NCVF shifted the ESPVR rightward while moderately reducing the slope, while a one-compartment simulation caused a leftward shift with severe reduction in the slope. Through tuning of the NCVF, Ls0,pas, and kECM, it was found that manipulation of the NCVF alone reproduced the deformation during acute ischemia and reperfusion, while additional manipulations of Ls0,pas and kECM were required to reproduce deformation during prolonged ischemia and reperfusion. Out of all segments with LGE>25% at the follow-up, the majority (68%) had higher estimated NCVF at the baseline than at the follow-up. Furthermore, the baseline NCVF correlated better with percentage LGE than NCVF did at the follow-up. Conclusion: We successfully used a two-compartment model for simulation of the ventricular pump and tissue mechanics in IHD. Patient-specific optimizations using regional myocardial deformation estimated the NCVF in a small cohort of MI patients in the acute and chronic phase after MI, while estimated NCVF values closely approximated the extent of the myocardial scar at the follow-up. In future studies, this approach can facilitate deformation imaging–based estimation of myocardial tissue properties in patients with cardiovascular diseases.
Collapse
Affiliation(s)
- Tijmen Koopsen
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
- *Correspondence: Tijmen Koopsen,
| | - Nick Van Osta
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| | - Tim Van Loon
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| | - Frans A. Van Nieuwenhoven
- Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| | - Frits W. Prinzen
- Department of Physiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| | - Bas R. Van Klarenbosch
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht, Utrecht, Netherlands
| | - Feddo P. Kirkels
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht, Utrecht, Netherlands
| | - Arco J. Teske
- Division of Heart and Lungs, Department of Cardiology, University Medical Center Utrecht, Utrecht, Netherlands
| | - Kevin Vernooy
- Department of Cardiology, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University Medical Center, Maastricht, Netherlands
- Department of Cardiology, Radboud University Medical Center, Nijmegen, Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht (CARIM), Maastricht University, Maastricht, Netherlands
| |
Collapse
|
5
|
Bracamonte JH, Saunders SK, Wilson JS, Truong UT, Soares JS. Patient-Specific Inverse Modeling of In Vivo Cardiovascular Mechanics with Medical Image-Derived Kinematics as Input Data: Concepts, Methods, and Applications. APPLIED SCIENCES-BASEL 2022; 12:3954. [PMID: 36911244 PMCID: PMC10004130 DOI: 10.3390/app12083954] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Inverse modeling approaches in cardiovascular medicine are a collection of methodologies that can provide non-invasive patient-specific estimations of tissue properties, mechanical loads, and other mechanics-based risk factors using medical imaging as inputs. Its incorporation into clinical practice has the potential to improve diagnosis and treatment planning with low associated risks and costs. These methods have become available for medical applications mainly due to the continuing development of image-based kinematic techniques, the maturity of the associated theories describing cardiovascular function, and recent progress in computer science, modeling, and simulation engineering. Inverse method applications are multidisciplinary, requiring tailored solutions to the available clinical data, pathology of interest, and available computational resources. Herein, we review biomechanical modeling and simulation principles, methods of solving inverse problems, and techniques for image-based kinematic analysis. In the final section, the major advances in inverse modeling of human cardiovascular mechanics since its early development in the early 2000s are reviewed with emphasis on method-specific descriptions, results, and conclusions. We draw selected studies on healthy and diseased hearts, aortas, and pulmonary arteries achieved through the incorporation of tissue mechanics, hemodynamics, and fluid-structure interaction methods paired with patient-specific data acquired with medical imaging in inverse modeling approaches.
Collapse
Affiliation(s)
- Johane H. Bracamonte
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
| | - Sarah K. Saunders
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
| | - John S. Wilson
- Department of Biomedical Engineering and Pauley Heart Center, Virginia Commonwealth University, Richmond, VA 23219, USA
| | - Uyen T. Truong
- Department of Pediatrics, School of Medicine, Children’s Hospital of Richmond at Virginia Commonwealth University, Richmond, VA 23219, USA
| | - Joao S. Soares
- Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, USA
- Correspondence:
| |
Collapse
|
6
|
Zhang Y, Adams J, Wang VY, Horwitz L, Tartibi M, Morgan AE, Kim J, Wallace AW, Weinsaft JW, Ge L, Ratcliffe MB. A finite element model of the cardiac ventricles with coupled circulation: Biventricular mesh generation with hexahedral elements, airbags and a functional mockup interface to the circulation. Comput Biol Med 2021; 137:104840. [PMID: 34508972 DOI: 10.1016/j.compbiomed.2021.104840] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 08/11/2021] [Accepted: 08/31/2021] [Indexed: 10/20/2022]
Abstract
INTRODUCTION Finite element (FE) mechanics models of the heart are becoming more sophisticated. However, there is lack of consensus about optimal element type and coupling of FE models to the circulation. We describe biventricular (left (LV) and right (RV) ventricles) FE mechanics model creation using hexahedral elements, airbags and a functional mockup interface (FMI) to lumped-parameter models of the circulation. METHODS Cardiac MRI (CMR) was performed in two healthy volunteers and a single patient with ischemic heart disease (IHD). CMR images were segmented and surfaced, meshing with hexahedral elements was performed with a "thin butterfly with septum" topology. LV and RV inflow and outflow airbags were coupled to lumped-parameter circulation models with an FMI interface. Pulmonary constriction (PAC) and vena cava occlusion (VCO) were simulated and end-systolic pressure-volume relations (ESPVR) were calculated. RESULTS Mesh construction was prompt with representative contouring and mesh adjustment requiring 32 and 26 min Respectively. The numbers of elements ranged from 4104 to 5184 with a representative Jacobian of 1.0026 ± 0.4531. Agreement between CMR-based surfaces and mesh was excellent with root-mean-squared error of 0.589 ± 0.321 mm. The LV ESPVR slope was 3.37 ± 0.09 in volunteers but 2.74 in the IHD patient. The effect of PAC and VCO on LV ESPVR was consistent with ventricular interaction (p = 0.0286). CONCLUSION Successful co-simulation using a biventricular FE mechanics model with hexahedral elements, airbags and an FMI interface to lumped-parameter model of the circulation was demonstrated. Future studies will include comparison of element type and study of cardiovascular pathologies and device therapies.
Collapse
Affiliation(s)
- Yue Zhang
- Department of Surgery, University of California, San Francisco, CA, USA; Department of Bioengineering, University of California, San Francisco, CA, USA; San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA
| | - Jennifer Adams
- School of Medicine, University of Texas Rio Grande Valley, Edinburg, TX, USA
| | - Vicky Y Wang
- Department of Surgery, University of California, San Francisco, CA, USA; Department of Bioengineering, University of California, San Francisco, CA, USA; San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA
| | - Lucas Horwitz
- San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA
| | | | - Ashley E Morgan
- Department of Surgery, University of Utah, Salt Lake City, UT, USA
| | - Jiwon Kim
- Department of Medicine, Weill Cornell Medicine, New York, NY, USA
| | - Arthur W Wallace
- Department of Anesthesia, University of California, San Francisco, CA, USA; San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA
| | | | - Liang Ge
- Department of Surgery, University of California, San Francisco, CA, USA; Department of Bioengineering, University of California, San Francisco, CA, USA; San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA
| | - Mark B Ratcliffe
- Department of Surgery, University of California, San Francisco, CA, USA; Department of Bioengineering, University of California, San Francisco, CA, USA; San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA.
| |
Collapse
|
7
|
Leong CO, Leong CN, Liew YM, Al Abed A, Aziz YFA, Chee KH, Sridhar GS, Dokos S, Lim E. The role of regional myocardial topography post-myocardial infarction on infarct extension. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3501. [PMID: 34057819 DOI: 10.1002/cnm.3501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Revised: 04/26/2021] [Accepted: 05/28/2021] [Indexed: 06/12/2023]
Abstract
Infarct extension involves necrosis of healthy myocardium in the border zone (BZ), progressively enlarging the infarct zone (IZ) and recruiting the remote zone (RZ) into the BZ, eventually leading to heart failure. The mechanisms underlying infarct extension remain unclear, but myocyte stretching has been suggested as the most likely cause. Using human patient-specific left-ventricular (LV) numerical simulations established from cardiac magnetic resonance imaging (MRI) of myocardial infarction (MI) patients, the correlation between infarct extension and regional mechanics abnormality was investigated by analysing the fibre stress-strain loops (FSSLs). FSSL abnormality was characterised using the directional regional external work (DREW) index, which measures FSSL area and loop direction. Sensitivity studies were also performed to investigate the effect of infarct stiffness on regional myocardial mechanics and potential for infarct extension. We found that infarct extension was correlated to severely abnormal FSSL in the form of counter-clockwise loop at the RZ close to the infarct, as indicated by negative DREW values. In regions demonstrating negative DREW values, we observed substantial fibre stretching in the isovolumic relaxation (IVR) phase accompanied by a reduced rate of systolic shortening. Such stretching in IVR phase in part of the RZ was due to its inability to withstand the high LV pressure that was still present and possibly caused by regional myocardial stiffness inhomogeneity. Further analysis revealed that the occurrence of severely abnormal FSSL due to IVR fibre stretching near the RZ-BZ boundary was due to a large amount of surrounding infarcted tissue, or an excessively stiff IZ.
Collapse
Affiliation(s)
- Chen Onn Leong
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur, Malaysia
| | - Chin Neng Leong
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, New South Wales, Australia
| | - Yih Miin Liew
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur, Malaysia
| | - Amr Al Abed
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, New South Wales, Australia
| | - Yang Faridah Abdul Aziz
- Department of Biomedical Imaging, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia
- University Malaya Research Imaging Centre, University of Malaya, Kuala Lumpur, Malaysia
| | - Kok Han Chee
- Department of Medicine, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia
| | | | - Socrates Dokos
- Graduate School of Biomedical Engineering, Faculty of Engineering, University of New South Wales, Sydney, New South Wales, Australia
| | - Einly Lim
- Department of Biomedical Engineering, Faculty of Engineering, University of Malaya, Kuala Lumpur, Malaysia
| |
Collapse
|
8
|
Wang VY, Tartibi M, Zhang Y, Selvaganesan K, Haraldsson H, Auger DA, Faraji F, Spaulding K, Takaba K, Collins A, Aguayo E, Saloner D, Wallace AW, Weinsaft JW, Epstein FH, Guccione J, Ge L, Ratcliffe MB. A kinematic model-based analysis framework for 3D Cine-DENSE-validation with an axially compressed gel phantom and application in sheep before and after antero-apical myocardial infarction. Magn Reson Med 2021; 86:2105-2121. [PMID: 34096083 DOI: 10.1002/mrm.28775] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Revised: 02/24/2021] [Accepted: 02/25/2021] [Indexed: 11/06/2022]
Abstract
PURPOSE Myocardial strain is increasingly used to assess left ventricular (LV) function. Incorporation of LV deformation into finite element (FE) modeling environment with subsequent strain calculation will allow analysis to reach its full potential. We describe a new kinematic model-based analysis framework (KMAF) to calculate strain from 3D cine-DENSE (displacement encoding with stimulated echoes) MRI. METHODS Cine-DENSE allows measurement of 3D myocardial displacement with high spatial accuracy. The KMAF framework uses cine cardiovascular magnetic resonance (CMR) to facilitate cine-DENSE segmentation, interpolates cine-DENSE displacement, and kinematically deforms an FE model to calculate strain. This framework was validated in an axially compressed gel phantom and applied in 10 healthy sheep and 5 sheep after myocardial infarction (MI). RESULTS Excellent Bland-Altman agreement of peak circumferential (Ecc ) and longitudinal (Ell ) strain (mean difference = 0.021 ± 0.04 and -0.006 ± 0.03, respectively), was found between KMAF estimates and idealized FE simulation. Err had a mean difference of -0.014 but larger variation (±0.12). Cine-DENSE estimated end-systolic (ES) Ecc , Ell and Err exhibited significant spatial variation for healthy sheep. Displacement magnitude was reduced on average by 27%, 42%, and 56% after MI in the remote, adjacent and MI regions, respectively. CONCLUSIONS The KMAF framework allows accurate calculation of 3D LV Ecc and Ell from cine-DENSE.
Collapse
Affiliation(s)
- Vicky Y Wang
- Veterans Affairs Medical Center, San Francisco, California, USA
| | - Mehrzad Tartibi
- Veterans Affairs Medical Center, San Francisco, California, USA
| | - Yue Zhang
- Veterans Affairs Medical Center, San Francisco, California, USA
| | - Kartiga Selvaganesan
- Department of Biomedical Engineering, University of Berkeley, Berkeley, California, USA
| | - Henrik Haraldsson
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Radiology, University of California, San Francisco, California, USA
| | - Daniel A Auger
- Department of Radiology and Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA.,Medical Metrics, Inc., Houston, Texas, USA
| | - Farshid Faraji
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Radiology, University of California, San Francisco, California, USA
| | | | - Kiyoaki Takaba
- Veterans Affairs Medical Center, San Francisco, California, USA
| | | | - Esteban Aguayo
- Veterans Affairs Medical Center, San Francisco, California, USA
| | - David Saloner
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Radiology, University of California, San Francisco, California, USA
| | - Arthur W Wallace
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Bioengineering, University of California, San Francisco, California, USA.,Department of Anesthesia, University of California, San Francisco, California, USA
| | | | - Frederick H Epstein
- Department of Radiology and Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Julius Guccione
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Bioengineering, University of California, San Francisco, California, USA.,Department of Surgery, University of California, San Francisco, California, USA
| | - Liang Ge
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Bioengineering, University of California, San Francisco, California, USA.,Department of Surgery, University of California, San Francisco, California, USA
| | - Mark B Ratcliffe
- Veterans Affairs Medical Center, San Francisco, California, USA.,Department of Bioengineering, University of California, San Francisco, California, USA.,Department of Surgery, University of California, San Francisco, California, USA.,Department of Medicine, University of California, San Francisco, California, USA
| |
Collapse
|
9
|
Finite-element based optimization of left ventricular passive stiffness in normal volunteers and patients after myocardial infarction: Utility of an inverse deformation gradient calculation of regional diastolic strain. J Mech Behav Biomed Mater 2021; 119:104431. [PMID: 33930653 DOI: 10.1016/j.jmbbm.2021.104431] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 02/23/2021] [Accepted: 02/26/2021] [Indexed: 12/15/2022]
Abstract
INTRODUCTION Left ventricular (LV) diastolic dysfunction (DD) is common after myocardial infarction (MI). Whereas current clinical assessment of DD relies on indirect markers including LV filling, finite element (FE) -based computational modeling directly measures regional diastolic stiffness. We hypothesized that an inverse deformation gradient (DG) method calculation of diastolic strain (IDGDS) allows the FE model-based calculation of regional diastolic stiffness (material parameters; MP) in post-MI patients with DD. METHODS Cardiac magnetic resonance (CMR) with tags (CSPAMM) and late gadolinium enhancement (LGE) was performed in 10 patients with post-MI DD and 10 healthy volunteers. The 3-dimensional (3D) LV DG from end-diastole (ED) to early diastolic filling (EDF; DGED→EDF) was calculated from CSPAMM. Diastolic strain was calculated from DGEDF→ED by inverting the DGED→EDF. FE models were created with MI and non-MI (remote; RM) regions determined by LGE. Guccione MPs C, and exponential fiber, bf, and transverse, bt , terms were optimized with IDGDS strain. RESULTS 3D circumferential and longitudinal diastolic strain (Ecc;Ell) calculated using IDGDS in CSPAMM obtained in volunteers and MI patients were [Formula: see text] = 0.27 ± 0.01, [Formula: see text] = 0.24 ± 0.03 and [Formula: see text] = 0.21 ± 0.02, and [Formula: see text] = 0.15 ± 0.02, respectively. MPs in the volunteer group were CH = 0.013 [0.001, 0.235] kPa, [Formula: see text] = 20.280 ± 4.994, and [Formula: see text] = 7.460 ± 2.171 and CRM = 0.0105 [0.010, 0.011] kPa, [Formula: see text] = 50.786 ± 13.511 (p = 0.0846), and [Formula: see text] = 17.355 ± 2.743 (p = 0.0208) in the remote myocardium of post-MI patients. CONCLUSION Diastolic strain, calculated from CSPAMM with IDGDS, enables calculation of FE model-based regional diastolic material parameters. Transverse stiffness of the remote myocardium, , may be a valuable new metric for determination of DD in patients after MI.
Collapse
|