1
|
Syomin FA, Galushka VA, Tsaturyan AK. Effect of strain-dependent conduction slowing on the re-entry formation and maintenance in cardiac muscle: 2D computer simulation. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3676. [PMID: 36562353 DOI: 10.1002/cnm.3676] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Revised: 12/06/2022] [Accepted: 12/19/2022] [Indexed: 06/17/2023]
Abstract
The effect of mechano-electrical feedback on re-entry formation and maintenance was studied using a model of myocardial electromechanics that accounts for two components of myocardial conductivity and delayed strain-dependent changes in membrane capacitance that causes a conduction slowing. Two scenarios were simulated in 2D numerical experiments: (i) propagation of an excitation-contraction wave beyond the edge of a nonconductive nonexcitable obstacle; (ii) circulation of a re-entry wave around a nonconductive nonexcitable obstacle. The simulations demonstrated that the delayed strain-dependent deceleration of the conduction waves promotes the detachment of the excitation-contraction waves from the sharp edge of an elongated obstacle and modulates the re-entry waves rotating around a compact obstacle. The data show that the mechano-electrical feedback, together with an increase in the stimulation frequency and an increase in the excitation threshold, is an arrhythmogenic factor that must be taken into account when analyzing the possibility of the re-entry formation.
Collapse
Affiliation(s)
- Fyodor A Syomin
- Institute of Mechanics, Lomonosov Moscow State University, Moscow, Russia
| | | | - Andrey K Tsaturyan
- Institute of Mechanics, Lomonosov Moscow State University, Moscow, Russia
| |
Collapse
|
2
|
Buonocunto M, Lyon A, Delhaas T, Heijman J, Lumens J. Electrophysiological effects of stretch-activated ion channels: a systematic computational characterization. J Physiol 2023. [PMID: 37665242 DOI: 10.1113/jp284439] [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: 01/30/2023] [Accepted: 08/07/2023] [Indexed: 09/05/2023] Open
Abstract
Cardiac electrophysiology and mechanics are strongly interconnected. Their interaction is, among others, mediated by mechano-electric feedback through stretch-activated ion channels (SACs). The electrophysiological changes induced by SACs may contribute to arrhythmogenesis, but the precise SAC-induced electrophysiological changes remain incompletely understood. Here, we provide a systematic characterization of stretch effects through three distinguished SACs on cardiac electrophysiology using computational modelling. We implemented potassium-selective, calcium-selective and non-selective SACs in the Tomek-Rodriguez-O'Hara-Rudy model of human ventricular electrophysiology. The model was calibrated to experimental data from isolated cardiomyocytes undergoing stretch, considering inter-species differences, and disease-related remodelling of SACs. SAC-mediated effects on the action potential (AP) were analysed by varying stretch amplitude, application timing and/or duration. Afterdepolarizations of different amplitudes were observed with transient 10-ms stretch stimuli of 15-18% applied during phase 4, while stretch ≥18% during phase 4 elicited triggered APs. Longer stimuli shifted the threshold of AP trigger during phase 4 to lower amplitudes, while shorter stimuli increased it. Continuous stretch provoked electrophysiological remodelling. Furthermore, stretch shortened duration or changed morphology of a subsequent electrically evoked AP, and, if applied during a vulnerable time window with sufficient amplitude, prevented its occurrence because of stretch-induced modulation of sodium and L-type calcium channel gating. These effects were more pronounced with disease-related SAC remodelling due to increased stretch sensitivity of diseased hearts. We showed that SACs may induce afterdepolarizations and triggered activities, and prevent subsequent AP generation or change its morphology. These effects depend on cardiomyocyte stretch characteristics and disease-related SACs remodelling and may contribute to cardiac arrhythmogenesis. KEY POINTS: The interplay between cardiac electrophysiology and mechanics is mediated by mechano-electric feedback through stretch-activated ion channels (SACs). These channels may be pro-arrhythmic, but their precise effect on electrophysiology remains unclear. Here we present a systematic in silico characterization of stretch effects through three SACs by implementing inter-species differences as well as disease-related remodelling of SACs in a novel computational model of human ventricular cardiomyocyte electrophysiology. Our simulations showed that, at the cellular level, SACs may provoke electrophysiological remodelling, afterdepolarizations, triggered activities, change the morphology or shorten subsequent electrically evoked action potentials. The model further suggests that a vulnerable window exists in which stretch prevents the following electrically triggered beat occurrence. The pro-arrhythmic effects of stretch strongly depend on disease-related SAC remodelling as well as on stretch characteristics, such as amplitude, time of application and duration. Our study helps in understanding the role of stretch in cardiac arrhythmogenesis and revealing the underlying cellular mechanisms.
Collapse
Affiliation(s)
- Melania Buonocunto
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, Netherlands
| | - Aurore Lyon
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, Netherlands
| | - Tammo Delhaas
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, Netherlands
| | - Jordi Heijman
- Department of Cardiology, Cardiovascular Research Institute Maastricht, Maastricht University Medical Centre, Maastricht, Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, Netherlands
| |
Collapse
|
3
|
Asiri F, Haque Siddiqui MI, Ali MA, Alam T, Dobrotă D, Chicea R, Dobrotă RD. Mathematical modeling of active contraction of the human cardiac myocyte: A review. Heliyon 2023; 9:e20065. [PMID: 37809539 PMCID: PMC10559823 DOI: 10.1016/j.heliyon.2023.e20065] [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: 06/06/2023] [Revised: 08/26/2023] [Accepted: 09/10/2023] [Indexed: 10/10/2023] Open
Abstract
Background and objective In this present research paper, a mathematical model has been developed to study myocyte contraction in the human cardiac muscle, using the Land model. Different parts of the human heart with a focus on the composition of the myocyte cells have been explored numerically to enabling us to determine the interaction of various parameters in the heart muscle. The main objective of the work is to direct the study of the Land model, which has been exploited to simulate the contraction of real human myocytes. Methods Mathematical models has been developed based on the Hill model and Huxley model. Myocyte contraction for different scenarios, such as in isometric tension and isotonic tension have been studied. Results It is found that increase in stretch, the peak active tension increases, in line with well-established length-dependent tension generation. Five parameters are selected: [Ca2+]T50, Tref, TRPN50, β0, and β1, which have been varied in between the range of -50%-100%, to examine the isometric effects of each parameter on the behavior of the tension developed in the intact myocyte cells, with the most sensitive parameter being [Ca2+]T50. Conclusion In conclusion, it is found that the Land model provides a good platform for the analysis of the active contraction of the human cardiac myocyte.
Collapse
Affiliation(s)
- Fisal Asiri
- Department of Mathematics, Taibah University, Medina, 42353, Saudi Arabia
| | | | - Masood Ashraf Ali
- Department of Industrial Engineering, College of Engineering, Prince Sattam Bin Abdulaziz University, Al-Kharj, 16273, Saudi Arabia
| | - Tabish Alam
- CSIR-Central Building Research Institute, Roorkee, 247667, India
| | - Dan Dobrotă
- Faculty of Engineering, Lucian Blaga University of Sibiu, 550024, Sibiu, Romania
| | - Radu Chicea
- Faculty of Medicine, Lucian Blaga University of Sibiu, 550024, Sibiu, Romania
| | | |
Collapse
|
4
|
Mehdi RR, Kumar M, Mendiola EA, Sadayappan S, Avazmohammadi R. Machine learning-based classification of cardiac relaxation impairment using sarcomere length and intracellular calcium transients. Comput Biol Med 2023; 163:107134. [PMID: 37379617 PMCID: PMC10525035 DOI: 10.1016/j.compbiomed.2023.107134] [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/06/2023] [Revised: 05/29/2023] [Accepted: 06/04/2023] [Indexed: 06/30/2023]
Abstract
Impaired relaxation of cardiomyocytes leads to diastolic dysfunction in the left ventricle. Relaxation velocity is regulated in part by intracellular calcium (Ca2+) cycling, and slower outflux of Ca2+ during diastole translates to reduced relaxation velocity of sarcomeres. Sarcomere length transient and intracellular calcium kinetics are integral parts of characterizing the relaxation behavior of the myocardium. However, a classifier tool that can separate normal cells from cells with impaired relaxation using sarcomere length transient and/or calcium kinetics remains to be developed. In this work, we employed nine different classifiers to classify normal and impaired cells, using ex-vivo measurements of sarcomere kinematics and intracellular calcium kinetics data. The cells were isolated from wild-type mice (referred to as normal) and transgenic mice expressing impaired left ventricular relaxation (referred to as impaired). We utilized sarcomere length transient data with a total of n = 126 cells (n = 60 normal cells and n = 66 impaired cells) and intracellular calcium cycling measurements with a total of n = 116 cells (n = 57 normal cells and n = 59 impaired cells) from normal and impaired cardiomyocytes as inputs to machine learning (ML) models for classification. We trained all ML classifiers with cross-validation method separately using both sets of input features, and compared their performance metrics. The performance of classifiers on test data showed that our soft voting classifier outperformed all other individual classifiers on both sets of input features, with 0.94 and 0.95 area under the receiver operating characteristic curves for sarcomere length transient and calcium transient, respectively, while multilayer perceptron achieved comparable scores of 0.93 and 0.95, respectively. However, the performance of decision tree, and extreme gradient boosting was found to be dependent on the set of input features used for training. Our findings highlight the importance of selecting appropriate input features and classifiers for the accurate classification of normal and impaired cells. Layer-wise relevance propagation (LRP) analysis demonstrated that the time to 50% contraction of the sarcomere had the highest relevance score for sarcomere length transient, whereas time to 50% decay of calcium had the highest relevance score for calcium transient input features. Despite the limited dataset, our study demonstrated satisfactory accuracy, suggesting that the algorithm can be used to classify relaxation behavior in cardiomyocytes when the potential relaxation impairment of the cells is unknown.
Collapse
Affiliation(s)
- Rana Raza Mehdi
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
| | - Mohit Kumar
- Heart, Lung, and Vascular Institute, Division of Cardiovascular Health and Disease, Department of Internal Medicine, Cincinnati, OH, USA; Department of Pharmacology and Systems Physiology, University of Cincinnati, Cincinnati, OH, USA
| | - Emilio A Mendiola
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA
| | - Sakthivel Sadayappan
- Heart, Lung, and Vascular Institute, Division of Cardiovascular Health and Disease, Department of Internal Medicine, Cincinnati, OH, USA; Department of Pharmacology and Systems Physiology, University of Cincinnati, Cincinnati, OH, USA
| | - Reza Avazmohammadi
- Department of Biomedical Engineering, Texas A&M University, College Station, TX 77843, USA; Department of Cardiovascular Sciences, Houston Methodist Research Institute, Houston, TX 77030, USA; J. Mike Walker '66 Department of Mechanical Engineering, Texas A&M University, College Station, TX 77843, USA.
| |
Collapse
|
5
|
Strocchi M, Longobardi S, Augustin CM, Gsell MAF, Petras A, Rinaldi CA, Vigmond EJ, Plank G, Oates CJ, Wilkinson RD, Niederer SA. Cell to whole organ global sensitivity analysis on a four-chamber heart electromechanics model using Gaussian processes emulators. PLoS Comput Biol 2023; 19:e1011257. [PMID: 37363928 DOI: 10.1371/journal.pcbi.1011257] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Accepted: 06/09/2023] [Indexed: 06/28/2023] Open
Abstract
Cardiac pump function arises from a series of highly orchestrated events across multiple scales. Computational electromechanics can encode these events in physics-constrained models. However, the large number of parameters in these models has made the systematic study of the link between cellular, tissue, and organ scale parameters to whole heart physiology challenging. A patient-specific anatomical heart model, or digital twin, was created. Cellular ionic dynamics and contraction were simulated with the Courtemanche-Land and the ToR-ORd-Land models for the atria and the ventricles, respectively. Whole heart contraction was coupled with the circulatory system, simulated with CircAdapt, while accounting for the effect of the pericardium on cardiac motion. The four-chamber electromechanics framework resulted in 117 parameters of interest. The model was broken into five hierarchical sub-models: tissue electrophysiology, ToR-ORd-Land model, Courtemanche-Land model, passive mechanics and CircAdapt. For each sub-model, we trained Gaussian processes emulators (GPEs) that were then used to perform a global sensitivity analysis (GSA) to retain parameters explaining 90% of the total sensitivity for subsequent analysis. We identified 45 out of 117 parameters that were important for whole heart function. We performed a GSA over these 45 parameters and identified the systemic and pulmonary peripheral resistance as being critical parameters for a wide range of volumetric and hemodynamic cardiac indexes across all four chambers. We have shown that GPEs provide a robust method for mapping between cellular properties and clinical measurements. This could be applied to identify parameters that can be calibrated in patient-specific models or digital twins, and to link cellular function to clinical indexes.
Collapse
Affiliation(s)
- Marina Strocchi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
| | - Stefano Longobardi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | | | | | - Argyrios Petras
- Johann Radon Institute for Computational and Applied Mathematics (RICAM), Linz, Austria
| | - Christopher A Rinaldi
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- Guy's and St Thomas' NHS Foundation Trust, London, United Kingdom
| | - Edward J Vigmond
- University of Bordeaux, CNRS, Bordeaux, Talence, France
- IHU Liryc, Bordeaux, Talence, France
| | - Gernot Plank
- Medical University of Graz, Graz, Austria
- BioTechMed-Graz, Graz, Austria
| | - Chris J Oates
- Newcastle University, Newcastle upon Tyne, United Kingdom
| | | | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
- National Heart and Lung Institute, Imperial College London, London, United Kingdom
- Alan Turing Institute, London, United Kingdom
| |
Collapse
|
6
|
Mora MT, Zaza A, Trenor B. Insights from an electro-mechanical heart failure cell model: Role of SERCA enhancement on arrhythmogenesis and myocyte contraction. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 230:107350. [PMID: 36689807 DOI: 10.1016/j.cmpb.2023.107350] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 12/27/2022] [Accepted: 01/10/2023] [Indexed: 06/17/2023]
Abstract
BACKGROUND AND OBJECTIVE Structural and electrical remodeling in heart failure predisposes the heart to ventricular arrhythmias. Computer modeling approaches, used to complement experimental results, can provide a more mechanistic knowledge of the biophysical phenomena underlying cardiac pathologies. Indeed, previous in-silico studies have improved the understanding of the electrical correlates of heart failure involved in arrhythmogenesis; however, information on the crosstalk between electrical activity, intracellular Ca2+ and contraction is still incomplete. This study aims to investigate the electro-mechanical behavior of virtual failing human ventricular myocytes to help in the development of therapies, which should ideally target pump failure and arrhythmias at the same time. METHODS We implemented characteristic remodeling of heart failure with reduced ejection fraction by including reported changes in ionic conductances, sarcomere function and cell structure (e.g. T-tubules disarray). Model parametrization was based on published experimental data and the outcome of simulations was validated against experimentally observed patterns. We focused on two aspects of myocardial dysfunction central in heart failure: altered force-frequency relationship and susceptibility to arrhythmogenic early afterdepolarizations. Because biological variability is a major problem in the generalization of in-silico findings based on a unique set of model parameters, we generated and evaluated a population of models. RESULTS The population-based approach is crucial in robust identification of parameters at the core of abnormalities and in generalizing the outcome of their correction. As compared to non-failing ones, failing myocytes had prolonged repolarization, a higher incidence of early afterdepolarizations, reduced contraction and a shallower force-frequency relationship, all features peculiar of heart failure. Component analysis applied to the model population identified reduced SERCA function as a relevant contributor to most of these derangements, which were largely reverted or diminished by restoration of SERCA function alone. CONCLUSIONS These simulated results encourage the development of strategies comprising SERCA stimulation and highlight the need to evaluate both electrical and mechanical outcomes.
Collapse
Affiliation(s)
- Maria Teresa Mora
- Centro de Investigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain
| | - Antonio Zaza
- Dipartimento di Biotecnologie e Bioscienze, Università degli Studi Milano-Bicocca, Italy; Unità di Fisiologia Cardiovascolare, IRCCs Istituto Auxologico Italiano, Italy
| | - Beatriz Trenor
- Centro de Investigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain.
| |
Collapse
|
7
|
Koivumäki JT, Hoffman J, Maleckar MM, Einevoll GT, Sundnes J. Computational cardiac physiology for new modelers: Origins, foundations, and future. Acta Physiol (Oxf) 2022; 236:e13865. [PMID: 35959512 DOI: 10.1111/apha.13865] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 08/04/2022] [Accepted: 08/05/2022] [Indexed: 01/29/2023]
Abstract
Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.
Collapse
Affiliation(s)
- Jussi T Koivumäki
- Faculty of Medicine and Health Technology, and Centre of Excellence in Body-on-Chip Research, Tampere University, Tampere, Finland
| | - Johan Hoffman
- Division of Computational Science and Technology, KTH Royal Institute of Technology, Stockholm, Sweden
| | - Mary M Maleckar
- Computational Physiology Department, Simula Research Laboratory, Oslo, Norway
| | - Gaute T Einevoll
- Centre for Integrative Neuroplasticity, University of Oslo, Oslo, Norway.,Department of Physics, University of Oslo, Oslo, Norway.,Department of Physics, Norwegian University of Life Sciences, Ås, Norway
| | - Joakim Sundnes
- Computational Physiology Department, Simula Research Laboratory, Oslo, Norway
| |
Collapse
|
8
|
Kim N, Pronto JD, Nickerson DP, Taberner AJ, Hunter PJ. A novel modular modeling approach for understanding different electromechanics between left and right heart in rat. Front Physiol 2022; 13:965054. [PMID: 36176770 PMCID: PMC9513479 DOI: 10.3389/fphys.2022.965054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 08/22/2022] [Indexed: 12/01/2022] Open
Abstract
While ion channels and transporters involved in excitation-contraction coupling have been linked and constructed as comprehensive computational models, validation of whether each individual component of a model can be reused has not been previously attempted. Here we address this issue while using a novel modular modeling approach to investigate the underlying mechanism for the differences between left ventricle (LV) and right ventricle (RV). Our model was developed from modules constructed using the module assembly principles of the CellML model markup language. The components of three existing separate models of cardiac function were disassembled as to create smaller modules, validated individually, and then the component parts were combined into a new integrative model of a rat ventricular myocyte. The model was implemented in OpenCOR using the CellML standard in order to ensure reproducibility. Simulated action potential (AP), Ca2+ transient, and tension were in close agreement with our experimental measurements: LV AP showed a prolonged duration and a more prominent plateau compared with RV AP; Ca2+ transient showed prolonged duration and slow decay in LV compared to RV; the peak value and relaxation of tension were larger and slower, respectively, in LV compared to RV. Our novel approach of module-based mathematical modeling has established that the ionic mechanisms underlying the APs and Ca2+ handling play a role in the variation in force production between ventricles. This simulation process also provides a useful way to reuse and elaborate upon existing models in order to develop a new model.
Collapse
Affiliation(s)
- Nari Kim
- NLRL for Innovative Cardiovascular Engineering, Department of Physiology, College of Medicine, Inje University, Busan, South Korea
- Cardiovascular and Metabolic Disease Center, Inje University, Busan, South Korea
- *Correspondence: Nari Kim,
| | - Julius D. Pronto
- NLRL for Innovative Cardiovascular Engineering, Department of Physiology, College of Medicine, Inje University, Busan, South Korea
- Cardiovascular and Metabolic Disease Center, Inje University, Busan, South Korea
| | - David P. Nickerson
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Andrew J. Taberner
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Peter J. Hunter
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
9
|
Tossas-Betancourt C, Li NY, Shavik SM, Afton K, Beckman B, Whiteside W, Olive MK, Lim HM, Lu JC, Phelps CM, Gajarski RJ, Lee S, Nordsletten DA, Grifka RG, Dorfman AL, Baek S, Lee LC, Figueroa CA. Data-driven computational models of ventricular-arterial hemodynamics in pediatric pulmonary arterial hypertension. Front Physiol 2022; 13:958734. [PMID: 36160862 PMCID: PMC9490558 DOI: 10.3389/fphys.2022.958734] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Accepted: 08/01/2022] [Indexed: 11/13/2022] Open
Abstract
Pulmonary arterial hypertension (PAH) is a complex disease involving increased resistance in the pulmonary arteries and subsequent right ventricular (RV) remodeling. Ventricular-arterial interactions are fundamental to PAH pathophysiology but are rarely captured in computational models. It is important to identify metrics that capture and quantify these interactions to inform our understanding of this disease as well as potentially facilitate patient stratification. Towards this end, we developed and calibrated two multi-scale high-resolution closed-loop computational models using open-source software: a high-resolution arterial model implemented using CRIMSON, and a high-resolution ventricular model implemented using FEniCS. Models were constructed with clinical data including non-invasive imaging and invasive hemodynamic measurements from a cohort of pediatric PAH patients. A contribution of this work is the discussion of inconsistencies in anatomical and hemodynamic data routinely acquired in PAH patients. We proposed and implemented strategies to mitigate these inconsistencies, and subsequently use this data to inform and calibrate computational models of the ventricles and large arteries. Computational models based on adjusted clinical data were calibrated until the simulated results for the high-resolution arterial models matched within 10% of adjusted data consisting of pressure and flow, whereas the high-resolution ventricular models were calibrated until simulation results matched adjusted data of volume and pressure waveforms within 10%. A statistical analysis was performed to correlate numerous data-derived and model-derived metrics with clinically assessed disease severity. Several model-derived metrics were strongly correlated with clinically assessed disease severity, suggesting that computational models may aid in assessing PAH severity.
Collapse
Affiliation(s)
| | - Nathan Y. Li
- Department of Pharmaceutical Sciences, University of Michigan, Ann Arbor, MI, United States
| | - Sheikh M. Shavik
- Department of Mechanical Engineering, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh
| | - Katherine Afton
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Brian Beckman
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Wendy Whiteside
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Mary K. Olive
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Heang M. Lim
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Jimmy C. Lu
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Christina M. Phelps
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Robert J. Gajarski
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - Simon Lee
- Department of Pediatrics, Nationwide Children’s Hospital, Columbus, OH, United States
| | - David A. Nordsletten
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
- Department of Surgery, University of Michigan, Ann Arbor, MI, United States
- School of Biomedical Engineering and Imaging Sciences, King’s College London, London, United Kingdom
| | - Ronald G. Grifka
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Adam L. Dorfman
- Department of Pediatrics, Division of Pediatric Cardiology, University of Michigan, Ann Arbor, MI, United States
| | - Seungik Baek
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, MI, United States
| | - C. Alberto Figueroa
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
- Department of Surgery, University of Michigan, Ann Arbor, MI, United States
| |
Collapse
|
10
|
Niarakis A, Waltemath D, Glazier J, Schreiber F, Keating SM, Nickerson D, Chaouiya C, Siegel A, Noël V, Hermjakob H, Helikar T, Soliman S, Calzone L. Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology. Brief Bioinform 2022; 23:bbac212. [PMID: 35671510 PMCID: PMC9294410 DOI: 10.1093/bib/bbac212] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Revised: 04/20/2022] [Accepted: 05/06/2022] [Indexed: 11/14/2022] Open
Abstract
Computational models are often employed in systems biology to study the dynamic behaviours of complex systems. With the rise in the number of computational models, finding ways to improve the reusability of these models and their ability to reproduce virtual experiments becomes critical. Correct and effective model annotation in community-supported and standardised formats is necessary for this improvement. Here, we present recent efforts toward a common framework for annotated, accessible, reproducible and interoperable computational models in biology, and discuss key challenges of the field.
Collapse
Affiliation(s)
- Anna Niarakis
- Université Paris-Saclay, Laboratoire Européen de Recherche pour la Polyarthrite rhumatoïde - Genhotel, Univ Evry, Evry, France
- Lifeware Group, Inria, Saclay-île de France, 91120 Palaiseau, France
| | - Dagmar Waltemath
- Department of Medical Informatics, University Medicine Greifswald, Greifswald, Germany
| | - James Glazier
- Biocomplexity Institute and Department of Intelligent Systems Engineering, Indiana University, Bloomington, IN, USA
| | - Falk Schreiber
- Department of Computer and Information Science, University of Konstanz, Konstanz, Germany
- Faculty of Information Technology, Monash University, Clayton, Australia
| | | | - David Nickerson
- Auckland Bioengineering Institute, University of Auckland, Auckland, New Zealand
| | | | - Anne Siegel
- Univ Rennes, CNRS, Inria - IRISA lab. Rennes
| | - Vincent Noël
- Institut Curie, PSL Research University, Paris, France
- INSERM, U900, Paris, France
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, Paris, France
| | - Henning Hermjakob
- EMBL-European Bioinformatics Institute, Wellcome Genome Campus, Cambridge, UK
| | - Tomáš Helikar
- Department of Biochemistry, University of Nebraska-Lincoln, Lincoln, NE, 68588, USA
| | - Sylvain Soliman
- Lifeware Group, Inria, Saclay-île de France, 91120 Palaiseau, France
| | - Laurence Calzone
- Institut Curie, PSL Research University, Paris, France
- INSERM, U900, Paris, France
- MINES ParisTech, PSL Research University, CBIO-Centre for Computational Biology, Paris, France
| |
Collapse
|
11
|
Guan D, Gao H, Cai L, Luo X. A new active contraction model for the myocardium using a modified hill model. Comput Biol Med 2022; 145:105417. [DOI: 10.1016/j.compbiomed.2022.105417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 02/11/2022] [Accepted: 02/21/2022] [Indexed: 11/16/2022]
|
12
|
The role of mechano-electric feedbacks and hemodynamic coupling in scar-related ventricular tachycardia. Comput Biol Med 2022; 142:105203. [DOI: 10.1016/j.compbiomed.2021.105203] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 12/29/2021] [Accepted: 12/29/2021] [Indexed: 12/17/2022]
|
13
|
Colli Franzone P, Pavarino LF, Scacchi S. Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3285. [PMID: 31808301 DOI: 10.1002/cnm.3285] [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: 01/17/2019] [Revised: 10/11/2019] [Accepted: 11/10/2019] [Indexed: 06/10/2023]
Abstract
Recent advances in the development of noninvasive cardiac imaging technologies have made it possible to measure longitudinal and circumferential strains at a high spatial resolution also at intramural level. Local mechanical activation times derived from these strains can be used as noninvasive estimates of electrical activation, in order to determine, eg, the origin of premature ectopic beats during focal arrhythmias or the pathway of reentrant circuits. The aim of this work is to assess the reliability of mechanical activation time markers derived from longitudinal and circumferential strains, denoted by ATell and ATecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (ATv ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by ATta and ATeff , respectively. Our numerical simulations are based on a strongly coupled electromechanical model, including bidomain representation of the cardiac tissue, mechanoelectric (ie, stretch-activated channels) and geometric feedbacks, transversely isotropic strain energy function for the description of passive mechanics and detailed membrane and excitation-contraction coupling models. The results have shown that, during endocardial and epicardial ectopic stimulations, all the mechanical markers considered are highly correlated with ATv , exhibiting correlation coefficients larger than 0.8. However, during multiple endocardial stimulations, mimicking the ventricular sinus rhythm, the mechanical markers are less correlated with the electrical activation time, because of the more complex resulting excitation sequence. Moreover, the inspection of the endocardial and epicardial isochrones has shown that the ATell and ATecc mechanical activation sequences reproduce only some qualitative features of the electrical activation sequence, such as the areas of early and late activation, but in some cases, they might yield wrong excitation sources and significantly different isochrones patterns.
Collapse
Affiliation(s)
| | - Luca F Pavarino
- Dipartimento di Matematica, Università di Milano, Milano, Italy
| | - Simone Scacchi
- Dipartimento di Matematica, Università di Milano, Milano, Italy
| |
Collapse
|
14
|
Forouzandehmehr M, Koivumäki JT, Hyttinen J, Paci M. A mathematical model of hiPSC cardiomyocytes electromechanics. Physiol Rep 2021; 9:e15124. [PMID: 34825519 PMCID: PMC8617339 DOI: 10.14814/phy2.15124] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 10/18/2021] [Accepted: 11/02/2021] [Indexed: 01/21/2023] Open
Abstract
Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are becoming instrumental in cardiac research, human-based cell level cardiotoxicity tests, and developing patient-specific care. As one of the principal functional readouts is contractility, we propose a novel electromechanical hiPSC-CM computational model named the hiPSC-CM-CE. This model comprises a reparametrized version of contractile element (CE) by Rice et al., 2008, with a new passive force formulation, integrated into a hiPSC-CM electrophysiology formalism by Paci et al. in 2020. Our simulated results were validated against in vitro data reported for hiPSC-CMs at matching conditions from different labs. Specifically, key action potential (AP) and calcium transient (CaT) biomarkers simulated by the hiPSC-CM-CE model were within the experimental ranges. On the mechanical side, simulated cell shortening, contraction-relaxation kinetic indices (RT50 and RT25 ), and the amplitude of tension fell within the experimental intervals. Markedly, as an inter-scale analysis, correct classification of the inotropic effects due to non-cardiomyocytes in hiPSC-CM tissues was predicted on account of the passive force expression introduced to the CE. Finally, the physiological inotropic effects caused by Verapamil and Bay-K 8644 and the aftercontractions due to the early afterdepolarizations (EADs) were simulated and validated against experimental data. In the future, the presented model can be readily expanded to take in pharmacological trials and genetic mutations, such as those involved in hypertrophic cardiomyopathy, and study arrhythmia trigger mechanisms.
Collapse
Affiliation(s)
| | - Jussi T. Koivumäki
- Faculty of Medicine and Health TechnologyTampere UniversityTampereFinland
| | - Jari Hyttinen
- Faculty of Medicine and Health TechnologyTampere UniversityTampereFinland
| | - Michelangelo Paci
- Faculty of Medicine and Health TechnologyTampere UniversityTampereFinland
| |
Collapse
|
15
|
Yoneda K, Okada JI, Watanabe M, Sugiura S, Hisada T, Washio T. A Multiple Step Active Stiffness Integration Scheme to Couple a Stochastic Cross-Bridge Model and Continuum Mechanics for Uses in Both Basic Research and Clinical Applications of Heart Simulation. Front Physiol 2021; 12:712816. [PMID: 34483965 PMCID: PMC8414591 DOI: 10.3389/fphys.2021.712816] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2021] [Accepted: 07/06/2021] [Indexed: 11/25/2022] Open
Abstract
In a multiscale simulation of a beating heart, the very large difference in the time scales between rapid stochastic conformational changes of contractile proteins and deterministic macroscopic outcomes, such as the ventricular pressure and volume, have hampered the implementation of an efficient coupling algorithm for the two scales. Furthermore, the consideration of dynamic changes of muscle stiffness caused by the cross-bridge activity of motor proteins have not been well established in continuum mechanics. To overcome these issues, we propose a multiple time step scheme called the multiple step active stiffness integration scheme (MusAsi) for the coupling of Monte Carlo (MC) multiple steps and an implicit finite element (FE) time integration step. The method focuses on the active tension stiffness matrix, where the active tension derivatives concerning the current displacements in the FE model are correctly integrated into the total stiffness matrix to avoid instability. A sensitivity analysis of the number of samples used in the MC model and the combination of time step sizes confirmed the accuracy and robustness of MusAsi, and we concluded that the combination of a 1.25 ms FE time step and 0.005 ms MC multiple steps using a few hundred motor proteins in each finite element was appropriate in the tradeoff between accuracy and computational time. Furthermore, for a biventricular FE model consisting of 45,000 tetrahedral elements, one heartbeat could be computed within 1.5 h using 320 cores of a conventional parallel computer system. These results support the practicality of MusAsi for uses in both the basic research of the relationship between molecular mechanisms and cardiac outputs, and clinical applications of perioperative prediction.
Collapse
Affiliation(s)
- Kazunori Yoneda
- Section Solutions Division, Healthcare Solutions Development Unit, Fujitsu Japan Ltd., Tokyo, Japan
| | - Jun-ichi Okada
- UT-Heart Inc., Kashiwa, Japan
- Future Center Initiative, University of Tokyo, Kashiwa, Japan
| | - Masahiro Watanabe
- Section Solutions Division, Healthcare Solutions Development Unit, Fujitsu Japan Ltd., Tokyo, Japan
| | | | | | - Takumi Washio
- UT-Heart Inc., Kashiwa, Japan
- Future Center Initiative, University of Tokyo, Kashiwa, Japan
| |
Collapse
|
16
|
Feng L, Gao H, Qi N, Danton M, Hill NA, Luo X. Fluid-structure interaction in a fully coupled three-dimensional mitral-atrium-pulmonary model. Biomech Model Mechanobiol 2021; 20:1267-1295. [PMID: 33770307 PMCID: PMC8298265 DOI: 10.1007/s10237-021-01444-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Accepted: 03/01/2021] [Indexed: 12/17/2022]
Abstract
This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.
Collapse
Affiliation(s)
- Liuyang Feng
- School of Mathematics and Statistics, University of Glasgow, Glasgow, G12 8SQ, UK.
| | - Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, G12 8SQ, UK
| | - Nan Qi
- Institute of Marine Science and Technology, Shandong University, Shangdong, 266237, People's Republic of China
| | - Mark Danton
- Department of Cardiac Surgery, Royal Hospital for Children, Glasgow, UK
| | - Nicholas A Hill
- School of Mathematics and Statistics, University of Glasgow, Glasgow, G12 8SQ, UK
| | - Xiaoyu Luo
- School of Mathematics and Statistics, University of Glasgow, Glasgow, G12 8SQ, UK
| |
Collapse
|
17
|
Mijailovich SM, Prodanovic M, Poggesi C, Geeves MA, Regnier M. Multiscale modeling of twitch contractions in cardiac trabeculae. J Gen Physiol 2021; 153:e202012604. [PMID: 33512405 PMCID: PMC7852458 DOI: 10.1085/jgp.202012604] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2020] [Revised: 08/31/2020] [Accepted: 11/20/2020] [Indexed: 12/23/2022] Open
Abstract
Understanding the dynamics of a cardiac muscle twitch contraction is complex because it requires a detailed understanding of the kinetic processes of the Ca2+ transient, thin-filament activation, and the myosin-actin cross-bridge chemomechanical cycle. Each of these steps has been well defined individually, but understanding how all three of the processes operate in combination is a far more complex problem. Computational modeling has the potential to provide detailed insight into each of these processes, how the dynamics of each process affect the complexity of contractile behavior, and how perturbations such as mutations in sarcomere proteins affect the complex interactions of all of these processes. The mechanisms involved in relaxation of tension during a cardiac twitch have been particularly difficult to discern due to nonhomogeneous sarcomere lengthening during relaxation. Here we use the multiscale MUSICO platform to model rat trabecular twitches. Validation of computational models is dependent on being able to simulate different experimental datasets, but there has been a paucity of data that can provide all of the required parameters in a single experiment, such as simultaneous measurements of force, intracellular Ca2+ transients, and sarcomere length dynamics. In this study, we used data from different studies collected under similar experimental conditions to provide information for all the required parameters. Our simulations established that twitches either in an isometric sarcomere or in fixed-length, multiple-sarcomere trabeculae replicate the experimental observations if models incorporate a length-tension relationship for the nonlinear series elasticity of muscle preparations and a scheme for thick-filament regulation. The thick-filament regulation assumes an off state in which myosin heads are parked onto the thick-filament backbone and are unable to interact with actin, a state analogous to the super-relaxed state. Including these two mechanisms provided simulations that accurately predict twitch contractions over a range of different conditions.
Collapse
Affiliation(s)
| | - Momcilo Prodanovic
- Bioengineering Research and Development Center, Kragujevac, Serbia
- Faculty of Engineering, University of Kragujevac, Kragujevac, Serbia
| | - Corrado Poggesi
- Department of Experimental & Clinical Medicine, University of Florence, Florence, Italy
| | | | - Michael Regnier
- Department of Bioengineering, University of Washington, Seattle, WA
| |
Collapse
|
18
|
Hazim A, Belhamadia Y, Dubljevic S. A Simulation Study of the Role of Mechanical Stretch in Arrhythmogenesis during Cardiac Alternans. Biophys J 2020; 120:109-121. [PMID: 33248131 PMCID: PMC7820729 DOI: 10.1016/j.bpj.2020.11.018] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 11/03/2020] [Accepted: 11/13/2020] [Indexed: 12/20/2022] Open
Abstract
The deformation of the heart tissue due to the contraction can modulate the excitation, a phenomenon referred to as mechanoelectrical feedback (MEF), via stretch-activated channels. The effects of MEF on the electrophysiology at high pacing rates are shown to be proarrhythmic in general. However, more studies need to be done to elucidate the underlying mechanism. In this work, we investigate the effects of MEF on cardiac alternans, which is an alternation in the width of the action potential that typically occurs when the heart is paced at high rates, using a biophysically detailed electromechanical model of cardiac tissue. We observe that the transition from spatially concordant alternans to spatially discordant alternans, which is more arrhythmogenic than concordant alternans, may occur in the presence of MEF and when its strength is sufficiently large. We show that this transition is due to the increase of the dispersion of conduction velocity. In addition, our results also show that the MEF effects, depending on the stretch-activated channels’ conductances and reversal potentials, can result in blocking action potential propagation.
Collapse
Affiliation(s)
- Azzam Hazim
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada
| | - Youssef Belhamadia
- Department of Mathematics and Statistics, American University of Sharjah, Sharjah, United Arab Emirates
| | - Stevan Dubljevic
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta, Canada.
| |
Collapse
|
19
|
Rumindo GK, Ohayon J, Croisille P, Clarysse P. In vivo estimation of normal left ventricular stiffness and contractility based on routine cine MR acquisition. Med Eng Phys 2020; 85:16-26. [PMID: 33081960 DOI: 10.1016/j.medengphy.2020.09.003] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2019] [Revised: 07/03/2020] [Accepted: 09/10/2020] [Indexed: 10/23/2022]
Abstract
Post-myocardial infarction remodeling process is known to alter the mechanical properties of the heart. Biomechanical parameters, such as tissue stiffness and contractility, would be useful for clinicians to better assess the severity of the diseased heart. However, these parameters are difficult to obtain in the current clinical practice. In this paper, we estimated subject-specific in vivo myocardial stiffness and contractility from 21 healthy volunteers, based on left ventricle models constructed from data acquired from routine cardiac MR acquisition only. The subject-specific biomechanical parameters were quantified using an inverse finite-element modelling approach. The personalized models were evaluated against relevant clinical metrics extracted from the MR data, such as circumferential strain, wall thickness and fractional thickening. We obtained the ranges of healthy biomechanical indices of 1.60 ± 0.22 kPa for left ventricular stiffness and 95.13 ± 14.56 kPa for left ventricular contractility. These reference normal values can be used for future model-based investigation on the stiffness and contractility of ischemic myocardium.
Collapse
Affiliation(s)
- Gerardo Kenny Rumindo
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1206, Lyon, France
| | - Jacques Ohayon
- University Savoie Mont-Blanc, Polytech Annecy-Chambéry and Laboratory TIMC-IMAG, UGA, CNRS UMR 5525, Grenoble, France
| | - Pierre Croisille
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1206, Lyon, France
| | - Patrick Clarysse
- Univ Lyon, INSA-Lyon, Université Claude Bernard Lyon 1, UJM-Saint Etienne, CNRS, Inserm, CREATIS UMR 5220, U1206, Lyon, France.
| |
Collapse
|
20
|
Regazzoni F, Dedè L, Quarteroni A. Biophysically detailed mathematical models of multiscale cardiac active mechanics. PLoS Comput Biol 2020; 16:e1008294. [PMID: 33027247 PMCID: PMC7571720 DOI: 10.1371/journal.pcbi.1008294] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2020] [Revised: 10/19/2020] [Accepted: 08/27/2020] [Indexed: 11/19/2022] Open
Abstract
We propose four novel mathematical models, describing the microscopic mechanisms of force generation in the cardiac muscle tissue, which are suitable for multiscale numerical simulations of cardiac electromechanics. Such models are based on a biophysically accurate representation of the regulatory and contractile proteins in the sarcomeres. Our models, unlike most of the sarcomere dynamics models that are available in the literature and that feature a comparable richness of detail, do not require the time-consuming Monte Carlo method for their numerical approximation. Conversely, the models that we propose only require the solution of a system of PDEs and/or ODEs (the most reduced of the four only involving 20 ODEs), thus entailing a significant computational efficiency. By focusing on the two models that feature the best trade-off between detail of description and identifiability of parameters, we propose a pipeline to calibrate such parameters starting from experimental measurements available in literature. Thanks to this pipeline, we calibrate these models for room-temperature rat and for body-temperature human cells. We show, by means of numerical simulations, that the proposed models correctly predict the main features of force generation, including the steady-state force-calcium and force-length relationships, the length-dependent prolongation of twitches and increase of peak force, the force-velocity relationship. Moreover, they correctly reproduce the Frank-Starling effect, when employed in multiscale 3D numerical simulation of cardiac electromechanics.
Collapse
Affiliation(s)
- Francesco Regazzoni
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Luca Dedè
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Alfio Quarteroni
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
- Mathematics Institute, École Polytechnique Fédérale de Lausanne, Av. Piccard, CH-1015 Lausanne, Switzerland (Professor Emeritus)
| |
Collapse
|
21
|
Regazzoni F, Dedè L, Quarteroni A. Active Force Generation in Cardiac Muscle Cells: Mathematical Modeling and Numerical Simulation of the Actin-Myosin Interaction. VIETNAM JOURNAL OF MATHEMATICS 2020; 49:87-118. [PMID: 34722731 PMCID: PMC8549950 DOI: 10.1007/s10013-020-00433-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/18/2019] [Accepted: 05/21/2020] [Indexed: 06/13/2023]
Abstract
Cardiac in silico numerical simulations are based on mathematical models describing the physical processes involved in the heart function. In this review paper, we critically survey biophysically-detailed mathematical models describing the subcellular mechanisms behind the generation of active force, that is the process by which the chemical energy of ATP (adenosine triphosphate) is transformed into mechanical work, thus making the muscle tissue contract. While presenting these models, that feature different levels of biophysical detail, we analyze the trade-off between the accuracy in the description of the subcellular mechanisms and the number of parameters that need to be estimated from experiments. Then, we focus on a generalized version of the classic Huxley model, built on the basis of models available in the literature, that is able to reproduce the main experimental characterizations associated to the time scales typical of a heartbeat-such as the force-velocity relationship and the tissue stiffness in response to small steps-featuring only four independent parameters. Finally, we show how those parameters can be calibrated starting from macroscopic measurements available from experiments.
Collapse
Affiliation(s)
- Francesco Regazzoni
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Luca Dedè
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
| | - Alfio Quarteroni
- MOX - Dipartimento di Matematica, Politecnico di Milano, P.zza Leonardo da Vinci 32, 20133 Milano, Italy
- Mathematics Institute, École Polytechnique Fédérale de Lausanne (EPFL), Av. Piccard, CH-1015 Lausanne, Switzerland
| |
Collapse
|
22
|
Li W. Biomechanics of infarcted left ventricle: a review of modelling. Biomed Eng Lett 2020; 10:387-417. [PMID: 32864174 DOI: 10.1007/s13534-020-00159-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2020] [Revised: 05/06/2020] [Accepted: 05/26/2020] [Indexed: 11/26/2022] Open
Abstract
Mathematical modelling in biomechanics of infarcted left ventricle (LV) serves as an indispensable tool for remodelling mechanism exploration, LV biomechanical property estimation and therapy assessment after myocardial infarction (MI). However, a review of mathematical modelling after MI has not been seen in the literature so far. In the paper, a systematic review of mathematical models in biomechanics of infarcted LV was established. The models include comprehensive cardiovascular system model, essential LV pressure-volume and stress-stretch models, constitutive laws for passive myocardium and scars, tension models for active myocardium, collagen fibre orientation optimization models, fibroblast and collagen fibre growth/degradation models and integrated growth-electro-mechanical model after MI. The primary idea, unique characteristics and key equations of each model were identified and extracted. Discussions on the models were provided and followed research issues on them were addressed. Considerable improvements in the cardiovascular system model, LV aneurysm model, coupled agent-based models and integrated electro-mechanical-growth LV model are encouraged. Substantial attention should be paid to new constitutive laws with respect to stress-stretch curve and strain energy function for infarcted passive myocardium, collagen fibre orientation optimization in scar, cardiac rupture and tissue damage and viscoelastic effect post-MI in the future.
Collapse
Affiliation(s)
- Wenguang Li
- School of Engineering, University of Glasgow, Glasgow, G12 8QQ UK
| |
Collapse
|
23
|
In-silico human electro-mechanical ventricular modelling and simulation for drug-induced pro-arrhythmia and inotropic risk assessment. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2020; 159:58-74. [PMID: 32710902 PMCID: PMC7848595 DOI: 10.1016/j.pbiomolbio.2020.06.007] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Revised: 06/08/2020] [Accepted: 06/28/2020] [Indexed: 12/28/2022]
Abstract
Human-based computational modelling and simulation are powerful tools to accelerate the mechanistic understanding of cardiac patho-physiology, and to develop and evaluate therapeutic interventions. The aim of this study is to calibrate and evaluate human ventricular electro-mechanical models for investigations on the effect of the electro-mechanical coupling and pharmacological action on human ventricular electrophysiology, calcium dynamics, and active contraction. The most recent models of human ventricular electrophysiology, excitation-contraction coupling, and active contraction were integrated, and the coupled models were calibrated using human experimental data. Simulations were then conducted using the coupled models to quantify the effects of electro-mechanical coupling and drug exposure on electrophysiology and force generation in virtual human ventricular cardiomyocytes and tissue. The resulting calibrated human electro-mechanical models yielded active tension, action potential, and calcium transient metrics that are in agreement with experiments for endocardial, epicardial, and mid-myocardial human samples. Simulation results correctly predicted the inotropic response of different multichannel action reference compounds and demonstrated that the electro-mechanical coupling improves the robustness of repolarisation under drug exposure compared to electrophysiology-only models. They also generated additional evidence to explain the partial mismatch between in-silico and in-vitro experiments on drug-induced electrophysiology changes. The human calibrated and evaluated modelling and simulation framework constructed in this study opens new avenues for future investigations into the complex interplay between the electrical and mechanical cardiac substrates, its modulation by pharmacological action, and its translation to tissue and organ models of cardiac patho-physiology.
Collapse
|
24
|
Guidry ME, Nickerson DP, Crampin EJ, Nash MP, Loiselle DS, Tran K. Insights From Computational Modeling Into the Contribution of Mechano-Calcium Feedback on the Cardiac End-Systolic Force-Length Relationship. Front Physiol 2020; 11:587. [PMID: 32547426 PMCID: PMC7273927 DOI: 10.3389/fphys.2020.00587] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Accepted: 05/11/2020] [Indexed: 11/23/2022] Open
Abstract
In experimental studies on cardiac tissue, the end-systolic force-length relation (ESFLR) has been shown to depend on the mode of contraction: isometric or isotonic. The isometric ESFLR is derived from isometric contractions spanning a range of muscle lengths while the isotonic ESFLR is derived from shortening contractions across a range of afterloads. The ESFLR of isotonic contractions consistently lies below its isometric counterpart. Despite the passing of over a hundred years since the first insight by Otto Frank, the mechanism(s) underlying this protocol-dependent difference in the ESFLR remain incompletely explained. Here, we investigate the role of mechano-calcium feedback in accounting for the difference between these two ESFLRs. Previous studies have compared the dynamics of isotonic contractions to those of a single isometric contraction at a length that produces maximum force, without considering isometric contractions at shorter muscle lengths. We used a mathematical model of cardiac excitation-contraction to simulate isometric and force-length work-loop contractions (the latter being the 1D equivalent of the whole-heart pressure-volume loop), and compared Ca2+ transients produced under equivalent force conditions. We found that the duration of the simulated Ca2+ transient increases with decreasing sarcomere length for isometric contractions, and increases with decreasing afterload for work-loop contractions. At any given force, the Ca2+ transient for an isometric contraction is wider than that during a work-loop contraction. By driving simulated work-loops with wider Ca2+ transients generated from isometric contractions, we show that the duration of muscle shortening was prolonged, thereby shifting the work-loop ESFLR toward the isometric ESFLR. These observations are explained by an increase in the rate of binding of Ca2+ to troponin-C with increasing force. However, the leftward shift of the work-loop ESFLR does not superimpose on the isometric ESFLR, leading us to conclude that while mechano-calcium feedback does indeed contribute to the difference between the two ESFLRs, it does not completely account for it.
Collapse
Affiliation(s)
- Megan E Guidry
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - David P Nickerson
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| | - Edmund J Crampin
- Systems Biology Laboratory, School of Mathematics and Statistics, The University of Melbourne, Melbourne, VIC, Australia.,ARC Centre of Excellence in Convergent Bio-Nano Science and Technology, School of Chemical and Biomedical Engineering, The University of Melbourne, Melbourne, VIC, Australia
| | - Martyn P Nash
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.,Department of Engineering Science, The University of Auckland, Auckland, New Zealand
| | - Denis S Loiselle
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand.,Department of Physiology, The University of Auckland, Auckland, New Zealand
| | - Kenneth Tran
- Auckland Bioengineering Institute, The University of Auckland, Auckland, New Zealand
| |
Collapse
|
25
|
Whiteley JP. An evaluation of some assumptions underpinning the bidomain equations of electrophysiology. MATHEMATICAL MEDICINE AND BIOLOGY-A JOURNAL OF THE IMA 2020; 37:262-302. [PMID: 31680135 DOI: 10.1093/imammb/dqz014] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Revised: 04/11/2019] [Accepted: 08/02/2019] [Indexed: 11/13/2022]
Abstract
Tissue level cardiac electrophysiology is usually modelled by the bidomain equations, or the monodomain simplification of the bidomain equations. One assumption made when deriving the bidomain equations is that both the intracellular and extracellular spaces are in electrical equilibrium. This assumption neglects the disturbance of this equilibrium in thin regions close to the cell membrane known as Debye layers. We first demonstrate that the governing equations at the cell, or microscale, level may be adapted to take account of these Debye layers with little additional complexity, provided the permittivity within the Debye layers satisfies certain conditions that are believed to be satisfied for biological cells. We then homogenize the microscale equations using a technique developed for an almost periodic microstructure. Cardiac tissue is usually modelled as sheets of cardiac fibres stacked on top of one another. A common assumption is that an orthogonal coordinate system can be defined at each point of cardiac tissue, where the first axis is in the fibre direction, the second axis is orthogonal to the first axis but lies in the sheet of cardiac fibres and the third axis is orthogonal to the cardiac sheet. It is assumed further that both the intracellular and extracellular conductivity tensors are diagonal with respect to these axes and that the diagonal entries of these tensors are constant across the whole tissue. Using the homogenization technique we find that this assumption is usually valid for cardiac tissue, but highlight situations where the assumption may not be valid.
Collapse
Affiliation(s)
- Jonathan P Whiteley
- Department of Computer Science, University of Oxford, Wolfson Building, Parks Road, Oxford OX1 3QD, UK
| |
Collapse
|
26
|
Cardiomyocyte calcium handling in health and disease: Insights from in vitro and in silico studies. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2020; 157:54-75. [PMID: 32188566 DOI: 10.1016/j.pbiomolbio.2020.02.008] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/29/2019] [Revised: 12/31/2019] [Accepted: 02/29/2020] [Indexed: 02/07/2023]
Abstract
Calcium (Ca2+) plays a central role in cardiomyocyte excitation-contraction coupling. To ensure an optimal electrical impulse propagation and cardiac contraction, Ca2+ levels are regulated by a variety of Ca2+-handling proteins. In turn, Ca2+ modulates numerous electrophysiological processes. Accordingly, Ca2+-handling abnormalities can promote cardiac arrhythmias via various mechanisms, including the promotion of afterdepolarizations, ion-channel modulation and structural remodeling. In the last 30 years, significant improvements have been made in the computational modeling of cardiomyocyte Ca2+ handling under physiological and pathological conditions. However, numerous questions involving the Ca2+-dependent regulation of different macromolecular complexes, cross-talk between Ca2+-dependent regulatory pathways operating over a wide range of time scales, and bidirectional interactions between electrophysiology and mechanics remain to be addressed by in vitro and in silico studies. A better understanding of disease-specific Ca2+-dependent proarrhythmic mechanisms may facilitate the development of improved therapeutic strategies. In this review, we describe the fundamental mechanisms of cardiomyocyte Ca2+ handling in health and disease, and provide an overview of currently available computational models for cardiomyocyte Ca2+ handling. Finally, we discuss important uncertainties and open questions about cardiomyocyte Ca2+ handling and highlight how synergy between in vitro and in silico studies may help to answer several of these issues.
Collapse
|
27
|
Sampedro-Puente DA, Fernandez-Bes J, Szentandrássy N, Nánási P, Taggart P, Pueyo E. Time Course of Low-Frequency Oscillatory Behavior in Human Ventricular Repolarization Following Enhanced Sympathetic Activity and Relation to Arrhythmogenesis. Front Physiol 2020; 10:1547. [PMID: 32009971 PMCID: PMC6971219 DOI: 10.3389/fphys.2019.01547] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Accepted: 12/09/2019] [Indexed: 11/21/2022] Open
Abstract
Background and Objectives: Recent studies in humans and dogs have shown that ventricular repolarization exhibits a low-frequency (LF) oscillatory pattern following enhanced sympathetic activity, which has been related to arrhythmic risk. The appearance of LF oscillations in ventricular repolarization is, however, not immediate, but it may take up to some minutes. This study seeks to characterize the time course of the action potential (AP) duration (APD) oscillatory behavior in response to sympathetic provocations, unveil its underlying mechanisms and establish a potential link to arrhythmogenesis under disease conditions. Materials and Methods: A representative set of human ventricular computational models coupling cellular electrophysiology, calcium dynamics, β-adrenergic signaling, and mechanics was built. Sympathetic provocation was modeled via phasic changes in β-adrenergic stimulation (β-AS) and mechanical stretch at Mayer wave frequencies within the 0.03–0.15 Hz band. Results: Our results show that there are large inter-individual differences in the time lapse for the development of LF oscillations in APD following sympathetic provocation, with some cells requiring just a few seconds and other cells needing more than 3 min. Whereas, the oscillatory response to phasic mechanical stretch is almost immediate, the response to β-AS is much more prolonged, in line with experimentally reported evidences, thus being this component the one driving the slow development of APD oscillations following enhanced sympathetic activity. If β-adrenoceptors are priorly stimulated, the time for APD oscillations to become apparent is remarkably reduced, with the oscillation time lapse being an exponential function of the pre-stimulation level. The major mechanism underlying the delay in APD oscillations appearance is related to the slow IKs phosphorylation kinetics, with its relevance being modulated by the IKs conductance of each individual cell. Cells presenting short oscillation time lapses are commonly associated with large APD oscillation magnitudes, which facilitate the occurrence of pro-arrhythmic events under disease conditions involving calcium overload and reduced repolarization reserve. Conclusions: The time course of LF oscillatory behavior of APD in response to increased sympathetic activity presents high inter-individual variability, which is associated with different expression and PKA phosphorylation kinetics of the IKs current. Short time lapses in the development of APD oscillations are associated with large oscillatory magnitudes and pro-arrhythmic risk under disease conditions.
Collapse
Affiliation(s)
| | | | - Norbert Szentandrássy
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Department of Dental Physiology and Pharmacology, Faculty of Dentistry, University of Debrecen, Debrecen, Hungary
| | - Péter Nánási
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary.,Department of Dental Physiology and Pharmacology, Faculty of Dentistry, University of Debrecen, Debrecen, Hungary
| | - Peter Taggart
- Department of Cardiovascular Sciences, University College London, London, United Kingdom
| | - Esther Pueyo
- BSICOS Group, I3A, IIS Aragón, University of Zaragoza, Zaragoza, Spain.,Center for Biomedical Research in the Network in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Zaragoza, Spain
| |
Collapse
|
28
|
Feng L, Gao H, Griffith B, Niederer S, Luo X. Analysis of a coupled fluid-structure interaction model of the left atrium and mitral valve. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3254. [PMID: 31454470 PMCID: PMC7003446 DOI: 10.1002/cnm.3254] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/08/2019] [Revised: 08/20/2019] [Accepted: 08/20/2019] [Indexed: 05/17/2023]
Abstract
We present a coupled left atrium-mitral valve model based on computed tomography scans with fibre-reinforced hyperelastic materials. Fluid-structure interaction is realised by using an immersed boundary-finite element framework. Effects of pathological conditions, eg, mitral valve regurgitation and atrial fibrillation, and geometric and structural variations, namely, uniform vs non-uniform atrial wall thickness and rule-based vs atlas-based fibre architectures, on the system are investigated. We show that in the case of atrial fibrillation, pulmonary venous flow reversal at late diastole disappears, and the filling waves at the left atrial appendage orifice during systole have reduced magnitude. In the case of mitral regurgitation, a higher atrial pressure and disturbed flows are seen, especially during systole, when a large regurgitant jet can be found with the suppressed pulmonary venous flow. We also show that both the rule-based and atlas-based fibre defining methods lead to similar flow fields and atrial wall deformations. However, the changes in wall thickness from non-uniform to uniform tend to underestimate the atrial deformation. Using a uniform but thickened wall also lowers the overall strain level. The flow velocity within the left atrial appendage, which is important in terms of appendage thrombosis, increases with the thickness of the left atrial wall. Energy analysis shows that the kinetic and dissipation energies of the flow within the left atrium are altered differently by atrial fibrillation and mitral valve regurgitation, providing a useful indication of the atrial performance in pathological situations.
Collapse
Affiliation(s)
- Liuyang Feng
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| | - Hao Gao
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| | - Boyce Griffith
- Departments of Mathematics, Applied Physical Sciences, and Biomedical EngineeringUniversity of North CarolinaChapel HillNorth CarolinaUSA
| | - Steven Niederer
- Department of Biomedical EngineeringKing's College LondonLondonUK
| | - Xiaoyu Luo
- School of Mathematics and StatisticsUniversity of GlasgowGlasgowUK
| |
Collapse
|
29
|
Pewowaruk RJ, Philip JL, Tewari SG, Chen CS, Nyaeme MS, Wang Z, Tabima DM, Baker AJ, Beard DA, Chesler NC. Multiscale Computational Analysis of Right Ventricular Mechanoenergetics. J Biomech Eng 2019; 140:2679646. [PMID: 30003251 DOI: 10.1115/1.4040044] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Indexed: 11/08/2022]
Abstract
Right ventricular (RV) failure, which occurs in the setting of pressure overload, is characterized by abnormalities in mechanical and energetic function. The effects of these cell- and tissue-level changes on organ-level RV function are unknown. The primary aim of this study was to investigate the effects of myofiber mechanics and mitochondrial energetics on organ-level RV function in the context of pressure overload using a multiscale model of the cardiovascular system. The model integrates the mitochondria-generated metabolite concentrations that drive intracellular actin-myosin cross-bridging and extracellular myocardial tissue mechanics in a biventricular heart model coupled with simple lumped parameter circulations. Three types of pressure overload were simulated and compared to experimental results. The computational model was able to capture a wide range of cardiovascular physiology and pathophysiology from mild RV dysfunction to RV failure. Our results confirm that, in response to pressure overload alone, the RV is able to maintain cardiac output (CO) and predict that alterations in either RV active myofiber mechanics or RV metabolite concentrations are necessary to decrease CO.
Collapse
Affiliation(s)
- Ryan J Pewowaruk
- Mem. ASME Biomedical Engineering, University of Wisconsin-Madison, 2145 Engineering Centers Building, 1550 Engineering Drive, Madison, WI 53706 e-mail:
| | - Jennifer L Philip
- Surgery, University of Wisconsin-Madison, , 1550 Engineering Drive, Madison, WI 53706 e-mail:
| | - Shivendra G Tewari
- Molecular & Integrative Physiology, University of Michigan-Ann Arbor, , North Campus Research Center, Ann Arbor, MI 48109-5622 e-mail:
| | - Claire S Chen
- Mechanical Engineering, University of Wisconsin-Madison, , 1550 Engineering Drive, Madison, WI 53706 e-mail:
| | - Mark S Nyaeme
- Biomedical Engineering, University of Wisconsin-Madison, , 1550 Engineering Drive, Madison, WI 53706 e-mail:
| | - Zhijie Wang
- Mechanical Engineering, Colorado State University, , Fort Collins, CO 80521 e-mail:
| | - Diana M Tabima
- Biomedical Engineering, University of Wisconsin-Madison, , 1550 Engineering Drive, Madison, WI 53706 e-mail:
| | - Anthony J Baker
- Medicine, University of California-San Francisco, , San Francisco, CA 94121; VA Medical Center, 4150 Clement St., San Francisco, CA 94121 e-mail:
| | - Daniel A Beard
- Molecular & Integrative Physiology, University of Michigan-Ann Arbor, , North Campus Research Center, Ann Arbor, MI 48109-5622 e-mail:
| | - Naomi C Chesler
- Fellow ASME Biomedical Engineering, University of Wisconsin-Madison Medicine, , 1550 Engineering Drive, Madison, WI 53706 e-mail:
| |
Collapse
|
30
|
Sampedro-Puente DA, Fernandez-Bes J, Porter B, van Duijvenboden S, Taggart P, Pueyo E. Mechanisms Underlying Interactions Between Low-Frequency Oscillations and Beat-to-Beat Variability of Celullar Ventricular Repolarization in Response to Sympathetic Stimulation: Implications for Arrhythmogenesis. Front Physiol 2019; 10:916. [PMID: 31427979 PMCID: PMC6687852 DOI: 10.3389/fphys.2019.00916] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2018] [Accepted: 07/04/2019] [Indexed: 12/12/2022] Open
Abstract
Background and Objectives: Enhanced beat-to-beat variability of ventricular repolarization (BVR) has been linked to arrhythmias and sudden cardiac death. Recent experimental studies on human left ventricular epicardial electrograms have shown that BVR closely interacts with low-frequency (LF) oscillations of activation recovery interval during sympathetic provocation. In this work human ventricular computational cell models are developed to reproduce the experimentally observed interactions between BVR and its LF oscillations, to assess underlying mechanisms and to establish a relationship with arrhythmic risk. Materials and Methods: A set of human ventricular action potential (AP) models covering a range of experimental electrophysiological characteristics was constructed. These models incorporated stochasticity in major ionic currents as well as descriptions of β-adrenergic stimulation and mechanical effects to investigate the AP response to enhanced sympathetic activity. Statistical methods based on Automatic Relevance Determination and Canonical Correlation Analysis were developed to unravel individual and common factors contributing to BVR and LF patterning of APD in response to sympathetic provocation. Results: Simulated results reproduced experimental evidences on the interactions between BVR and LF oscillations of AP duration (APD), with replication of the high inter-individual variability observed in both phenomena. ICaL, IKr and IK1 currents were identified as common ionic modulators of the inter-individual differences in BVR and LF oscillatory behavior and were shown to be crucial in determining susceptibility to arrhythmogenic events. Conclusions: The calibrated family of human ventricular cell models proposed in this study allows reproducing experimentally reported interactions between BVR and LF oscillations of APD. Ionic factors involving ICaL, IKr and IK1 currents are found to underlie correlated increments in both phenomena in response to sympathetic provocation. A link to arrhythmogenesis is established for concomitantly elevated levels of BVR and its LF oscillations.
Collapse
Affiliation(s)
| | | | - Bradley Porter
- Department of Imaging Sciences and Biomedical Engineering, Kings College London, London, United Kingdom
| | | | - Peter Taggart
- Department of Cardiovascular Sciences, University College London, London, United Kingdom
| | - Esther Pueyo
- BSICOS Group, I3A, IIS Aragón, University of Zaragoza, Zaragoza, Spain.,CIBER-BBN, Madrid, Spain
| |
Collapse
|
31
|
Lewalle A, Land S, Merken JJ, Raafs A, Sepúlveda P, Heymans S, Kleinjans J, Niederer SA. Balance of Active, Passive, and Anatomical Cardiac Properties in Doxorubicin-Induced Heart Failure. Biophys J 2019; 117:2337-2348. [PMID: 31447110 PMCID: PMC6990149 DOI: 10.1016/j.bpj.2019.07.033] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Revised: 07/12/2019] [Accepted: 07/23/2019] [Indexed: 01/31/2023] Open
Abstract
Late-onset heart failure (HF) is a known side effect of doxorubicin chemotherapy. Typically, patients are diagnosed when already at an irreversible stage of HF, which allows few or no treatment options. Identifying the causes of compromised cardiac function in this patient group may improve early patient diagnosis and support treatment selection. To link doxorubicin-induced changes in cardiac cellular and tissue mechanical properties to overall cardiac function, we apply a multiscale biophysical biomechanics model of the heart to measure the plausibility of changes in model parameters representing the passive, active, or anatomical properties of the left ventricle for reproducing measured patient phenotypes. We create representative models of healthy controls (N = 10) and patients with HF induced by (N = 22) or unrelated to (N = 25) doxorubicin therapy. The model predicts that HF in the absence of doxorubicin is characterized by a 2- to 3-fold stiffness increase, decreased tension (0–20%), and ventricular dilation (of order 10–30%). HF due to doxorubicin was similar but showed stronger bias toward reduced active contraction (10–30%) and less dilation (0–20%). We find that changes in active, passive, and anatomical properties all play a role in doxorubicin-induced cardiotoxicity phenotypes. Differences in parameter changes between patient groups are consistent with doxorubicin cardiotoxicity having a greater dependence on reduced cellular contraction and less anatomical remodeling than HF not caused by doxorubicin.
Collapse
Affiliation(s)
- Alexandre Lewalle
- Department of Biomedical Engineering, St Thomas's Hospital, King's College London, London, United Kingdom
| | - Sander Land
- Department of Biomedical Engineering, St Thomas's Hospital, King's College London, London, United Kingdom
| | - Jort J Merken
- Department of Cardiology, Maastricht University, Maastricht, the Netherlands
| | - Anne Raafs
- Department of Cardiology, Maastricht University, Maastricht, the Netherlands
| | - Pilar Sepúlveda
- Regenerative Medicine and Heart Transplantation Unit, Instituto de Investigación Sanitaria La Fe, Valencia, Spain
| | - Stéphane Heymans
- Department of Cardiology, Maastricht University, Maastricht, the Netherlands
| | - Jos Kleinjans
- Department of Toxicogenomics, Maastricht University, Maastricht, the Netherlands
| | - Steven A Niederer
- Department of Biomedical Engineering, St Thomas's Hospital, King's College London, London, United Kingdom.
| |
Collapse
|
32
|
Syomin FA, Zberia MV, Tsaturyan AK. Multiscale simulation of the effects of atrioventricular block and valve diseases on heart performance. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2019; 35:e3216. [PMID: 31083764 DOI: 10.1002/cnm.3216] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Revised: 04/12/2019] [Accepted: 05/08/2019] [Indexed: 06/09/2023]
Abstract
A new mathematical model of the cardiovascular system is proposed. The left ventricle is described by an axisymmetric multiscale model where myocardium is treated as an incompressible transversely isotropic medium with a realistic distribution of fibre orientation. Active tension and its regulation by Ca2+ ions are described by our recent kinetic model. A lumped parameter model is used for the simulation of blood circulation, in which the left and right atria and the right ventricle are described by a system of ordinary differential equations for active pressure-volume relationships. The stress and strain of the left ventricle myocardium were calculated by the finite element method implemented by the authors. The changes in the haemodynamics upon changes in preload of a healthy heart, upon physical exercise, and in case of atrioventricular block with different types of arrhythmias were simulated. To simulate the effect of stenosis or regurgitation of the aortic or mitral valves, the hydraulic and inertial flow resistances of the heart valves were set as functions of their orifice areas. The model reproduced a number of phenomena observed in clinical practice, including the classification of the severity of valve disease.
Collapse
Affiliation(s)
- Fyodor A Syomin
- Department of Biomechanics, Institute of Mechanics, M.V. Lomonosov Moscow State University, 1 Mitchurinsky Prosp., Moscow, 119192, Russian Federation
- Peoples' Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow, 117198, Russian Federation
| | - Maria V Zberia
- Department of Biomechanics, Institute of Mechanics, M.V. Lomonosov Moscow State University, 1 Mitchurinsky Prosp., Moscow, 119192, Russian Federation
| | - Andrey K Tsaturyan
- Department of Biomechanics, Institute of Mechanics, M.V. Lomonosov Moscow State University, 1 Mitchurinsky Prosp., Moscow, 119192, Russian Federation
| |
Collapse
|
33
|
Hazim A, Belhamadia Y, Dubljevic S. Effects of mechano-electrical feedback on the onset of alternans: A computational study. CHAOS (WOODBURY, N.Y.) 2019; 29:063126. [PMID: 31266317 DOI: 10.1063/1.5095778] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2019] [Accepted: 06/05/2019] [Indexed: 06/09/2023]
Abstract
Cardiac alternans is a heart rhythm instability that is associated with cardiac arrhythmias and may lead to sudden cardiac death. The onset of this instability, which is linked to period-doubling bifurcation and may be a route to chaos, is of particular interest. Mechano-electric feedback depicts the effects of tissue deformation on cardiac excitation. The main effect of mechano-electric feedback is delivered via the so-called stretch-activated ion channels and is caused by stretch-activated currents. Mechano-electric feedback, which is believed to have proarrhythmic and antiarrhythmic effects on cardiac electrophysiology, affects the action potential duration in a manner dependent on cycle length, but the mechanisms by which this occurs remain to be elucidated. In this study, a biophysically detailed electromechanical model of cardiac tissue is employed to show how a stretch-activated current can affect the action potential duration at cellular and tissue levels, illustrating its effects on the onset of alternans. Also, using a two-dimensional iterated map that incorporates stretch-activated current effects, we apply linear stability analysis to study the stability of the bifurcation. We show that alternans bifurcation can be prevented depending on the strength of the stretch-activated current.
Collapse
Affiliation(s)
- Azzam Hazim
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta T6G 2V2, Canada
| | - Youssef Belhamadia
- Department of Mathematics and Statistics, American University of Sharjah, Sharjah, United Arab Emirates
| | - Stevan Dubljevic
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta T6G 2V4, Canada
| |
Collapse
|
34
|
Niederer SA, Campbell KS, Campbell SG. A short history of the development of mathematical models of cardiac mechanics. J Mol Cell Cardiol 2019; 127:11-19. [PMID: 30503754 PMCID: PMC6525149 DOI: 10.1016/j.yjmcc.2018.11.015] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/30/2018] [Revised: 11/02/2018] [Accepted: 11/21/2018] [Indexed: 11/15/2022]
Abstract
Cardiac mechanics plays a crucial role in atrial and ventricular function, in the regulation of growth and remodelling, in the progression of disease, and the response to treatment. The spatial scale of the critical mechanisms ranges from nm (molecules) to cm (hearts) with the fastest events occurring in milliseconds (molecular events) and the slowest requiring months (growth and remodelling). Due to its complexity and importance, cardiac mechanics has been studied extensively both experimentally and through mathematical models and simulation. Models of cardiac mechanics evolved from seminal studies in skeletal muscle, and developed into cardiac specific, species specific, human specific and finally patient specific calculations. These models provide a formal framework to link multiple experimental assays recorded over nearly 100 years into a single unified representation of cardiac function. This review first provides a summary of the proteins, physiology and anatomy involved in the generation of cardiac pump function. We then describe the evolution of models of cardiac mechanics starting with the early theoretical frameworks describing the link between sarcomeres and muscle contraction, transitioning through myosin-level models to calcium-driven systems, and ending with whole heart patient-specific models.
Collapse
Affiliation(s)
| | - Kenneth S Campbell
- Department of Physiology and Division of Cardiovascular Medicine, University of Kentucky, Lexington, USA
| | - Stuart G Campbell
- Departments of Biomedical Engineering and Cellular and Molecular Physiology, Yale University, New Haven, USA
| |
Collapse
|
35
|
Arens S, Dierckx H, Panfilov AV. GEMS: A Fully Integrated PETSc-Based Solver for Coupled Cardiac Electromechanics and Bidomain Simulations. Front Physiol 2018; 9:1431. [PMID: 30386252 PMCID: PMC6198176 DOI: 10.3389/fphys.2018.01431] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2018] [Accepted: 09/20/2018] [Indexed: 01/23/2023] Open
Abstract
Cardiac contraction is coordinated by a wave of electrical excitation which propagates through the heart. Combined modeling of electrical and mechanical function of the heart provides the most comprehensive description of cardiac function and is one of the latest trends in cardiac research. The effective numerical modeling of cardiac electromechanics remains a challenge, due to the stiffness of the electrical equations and the global coupling in the mechanical problem. Here we present a short review of the inherent assumptions made when deriving the electromechanical equations, including a general representation for deformation-dependent conduction tensors obeying orthotropic symmetry, and then present an implicit-explicit time-stepping approach that is tailored to solving the cardiac mono- or bidomain equations coupled to electromechanics of the cardiac wall. Our approach allows to find numerical solutions of the electromechanics equations using stable and higher order time integration. Our methods are implemented in a monolithic finite element code GEMS (Ghent Electromechanics Solver) using the PETSc library that is inherently parallelized for use on high-performance computing infrastructure. We tested GEMS on standard benchmark computations and discuss further development of our software.
Collapse
Affiliation(s)
- Sander Arens
- Department of Physics and Astronomy, Ghent University, Ghent, Belgium
| | - Hans Dierckx
- Department of Physics and Astronomy, Ghent University, Ghent, Belgium
| | - Alexander V Panfilov
- Department of Physics and Astronomy, Ghent University, Ghent, Belgium.,Laboratory of Computational Biology and Medicine, Ural Federal University, Ekaterinburg, Russia
| |
Collapse
|
36
|
Hazim A, Belhamadia Y, Dubljevic S. Mechanical perturbation control of cardiac alternans. Phys Rev E 2018; 97:052407. [PMID: 29906969 DOI: 10.1103/physreve.97.052407] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2017] [Indexed: 12/27/2022]
Abstract
Cardiac alternans is a disturbance in heart rhythm that is linked to the onset of lethal cardiac arrhythmias. Mechanical perturbation control has been recently used to suppress alternans in cardiac tissue of relevant size. In this control strategy, cardiac tissue mechanics are perturbed via active tension generated by the heart's electrical activity, which alters the tissue's electric wave profile through mechanoelectric coupling. We analyze the effects of mechanical perturbation on the dynamics of a map model that couples the membrane voltage and active tension systems at the cellular level. Therefore, a two-dimensional iterative map of the heart beat-to-beat dynamics is introduced, and a stability analysis of the system of coupled maps is performed in the presence of a mechanical perturbation algorithm. To this end, a bidirectional coupling between the membrane voltage and active tension systems in a single cardiac cell is provided, and a discrete form of the proposed control algorithm, that can be incorporated in the coupled maps, is derived. In addition, a realistic electromechanical model of cardiac tissue is employed to explore the feasibility of suppressing alternans at cellular and tissue levels. Electrical activity is represented in two detailed ionic models, the Luo-Rudy 1 and the Fox models, while two active contractile tension models, namely a smooth variant of the Nash-Panfilov model and the Niederer-Hunter-Smith model, are used to represent mechanical activity in the heart. The Mooney-Rivlin passive elasticity model is employed to describe passive mechanical behavior of the myocardium.
Collapse
Affiliation(s)
- Azzam Hazim
- Department of Biomedical Engineering, University of Alberta, Edmonton, Alberta, Canada T6G 2V2
| | - Youssef Belhamadia
- Department of Mathematics and Statistics, American University of Sharjah, Sharjah, United Arab Emirates
| | - Stevan Dubljevic
- Department of Chemical and Materials Engineering, University of Alberta, Edmonton, Alberta, Canada T6G 2V4
| |
Collapse
|
37
|
High tension in sarcomeres hinders myocardial relaxation: A computational study. PLoS One 2018; 13:e0204642. [PMID: 30286135 PMCID: PMC6171862 DOI: 10.1371/journal.pone.0204642] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 09/12/2018] [Indexed: 11/19/2022] Open
Abstract
Experiments have shown that the relaxation phase of cardiac sarcomeres during an isometric twitch is prolonged in muscles that reached a higher peak tension. However, the mechanism is not completely understood. We hypothesize that the binding of calcium to troponin is enhanced by the tension in the thin filament, thus contributing to the prolongation of contraction upon higher peak tension generation. To test this hypothesis, we developed a computational model of sarcomere mechanics that incorporates tension-dependence of calcium binding. The model was used to simulate isometric twitch experiments with time dependency in the form of a two-state cross-bridge cycle model and a transient intracellular calcium concentration. In the simulations, peak isometric twitch tension appeared to increase linearly by 51.1 KPa with sarcomere length from 1.9 μm to 2.2 μm. Experiments showed an increase of 47.3 KPa over the same range of sarcomere lengths. The duration of the twitch also increased with both sarcomere length and peak intracellular calcium concentration, likely to be induced by the inherently coupled increase of the peak tension in the thin filament. In the model simulations, the time to 50% relaxation (tR50) increased over the range of sarcomere lengths from 1.9 μm to 2.2 μm by 0.11s, comparable to the increased duration of 0.12s shown in experiments. Model simulated tR50 increased by 0.12s over the range of peak intracellular calcium concentrations from 0.87 μM to 1.45 μM. Our simulation results suggest that the prolongation of contraction at higher tension is a result of the tighter binding of Ca2+ to troponin in areas under higher tension, thus delaying the deactivation of the troponin.
Collapse
|
38
|
Electromechanical effects of concentric hypertrophy on the left ventricle: A simulation study. Comput Biol Med 2018; 99:236-256. [DOI: 10.1016/j.compbiomed.2018.06.004] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2018] [Revised: 06/06/2018] [Accepted: 06/07/2018] [Indexed: 11/19/2022]
|
39
|
Regazzoni F, Dedè L, Quarteroni A. Active contraction of cardiac cells: a reduced model for sarcomere dynamics with cooperative interactions. Biomech Model Mechanobiol 2018; 17:1663-1686. [DOI: 10.1007/s10237-018-1049-0] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Accepted: 06/28/2018] [Indexed: 11/24/2022]
|
40
|
Gao H, Aderhold A, Mangion K, Luo X, Husmeier D, Berry C. Changes and classification in myocardial contractile function in the left ventricle following acute myocardial infarction. J R Soc Interface 2018; 14:rsif.2017.0203. [PMID: 28747397 PMCID: PMC5550971 DOI: 10.1098/rsif.2017.0203] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2017] [Accepted: 07/04/2017] [Indexed: 01/05/2023] Open
Abstract
In this research, we hypothesized that novel biomechanical parameters are discriminative in patients following acute ST-segment elevation myocardial infarction (STEMI). To identify these biomechanical biomarkers and bring computational biomechanics ‘closer to the clinic’, we applied state-of-the-art multiphysics cardiac modelling combined with advanced machine learning and multivariate statistical inference to a clinical database of myocardial infarction. We obtained data from 11 STEMI patients (ClinicalTrials.gov NCT01717573) and 27 healthy volunteers, and developed personalized mathematical models for the left ventricle (LV) using an immersed boundary method. Subject-specific constitutive parameters were achieved by matching to clinical measurements. We have shown, for the first time, that compared with healthy controls, patients with STEMI exhibited increased LV wall active tension when normalized by systolic blood pressure, which suggests an increased demand on the contractile reserve of remote functional myocardium. The statistical analysis reveals that the required patient-specific contractility, normalized active tension and the systolic myofilament kinematics have the strongest explanatory power for identifying the myocardial function changes post-MI. We further observed a strong correlation between two biomarkers and the changes in LV ejection fraction at six months from baseline (the required contractility (r = − 0.79, p < 0.01) and the systolic myofilament kinematics (r = 0.70, p = 0.02)). The clinical and prognostic significance of these biomechanical parameters merits further scrutinization.
Collapse
Affiliation(s)
- Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Andrej Aderhold
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Kenneth Mangion
- British Heart Foundation Glasgow Cardiovascular Research Centre, University of Glasgow, Glasgow, UK
| | - Xiaoyu Luo
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Colin Berry
- British Heart Foundation Glasgow Cardiovascular Research Centre, University of Glasgow, Glasgow, UK
| |
Collapse
|
41
|
Satriano A, Vigmond EJ, Schwartzman DS, Di Martino ES. Mechano-electric finite element model of the left atrium. Comput Biol Med 2018. [PMID: 29529527 DOI: 10.1016/j.compbiomed.2018.02.010] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Mechanical stretch plays a major role in modulating atrial function, being responsible for beat-by-beat responses to changes in chamber preload, enabling a prompt regulation of cardiac function. Mechano-electric coupling (MEC) operates through many mechanisms and has many targets, making it experimentally difficult to isolate causes and effects especially under sinus conditions where effects are more transient and subtle. Therefore, modelling is a powerful tool to help understand the role of MEC with respect to the atrial electromechanical interaction. We propose a cellular-based computational model of the left atrium that includes a strongly coupled MEC component and mitral flow component to account for correct pressure generation in the atrial chamber as a consequence of blood volume and contraction. The method was applied to a healthy porcine left atrium. Results of the strongly coupled simulation show that strains are higher in the areas adjacent to the mitral annulus, the rim of the appendage, around the pulmonary venous trunks and at the location of the Bachmann's bundle, approximately between the mitral annulus and the region where the venous tissue transitions into atrial. These are regions where arrhythmias are likely to originate. The role of stretch-activated channels was very small for sinus rhythm for the single cardiac beat simulation, although tension development was very sensitive to stretch. The method could be applied to investigate potential therapeutic interventions acting on the mechano-electrical properties of the left atrium.
Collapse
Affiliation(s)
- Alessandro Satriano
- Stephenson Cardiac Imaging Centre, The University of Calgary, 2500 University Drive NW, Calgary, Alberta, T2N 1N4, Canada
| | - Edward J Vigmond
- Department of Electrical and Computer Engineering, University of Calgary, 2500 University Dr NW, Calgary, Alberta T2N 1N4, Canada; LIRYC, Electrophysiology and Heart Modelling Institute, PTIB-Hopital Xavier Arnozan, Avenue Haut-Lévèque, Pessac, 33600, France; IMB, University of Bordeaux, 351 Cours de la Liberation, Talence, 33405, France
| | - David S Schwartzman
- Heart and Vascular Institute, University of Pittsburgh, UPMC Presbyterian, B535, Pittsburgh, PA 15213 2582, United States
| | - Elena S Di Martino
- Department of Civil Engineering, Libin Cardiovascular Institute of Alberta and Centre for Bioengineering Research and Education, The University of Calgary, 2500 University Drive NW, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
42
|
Aboelkassem Y, Campbell SG. Acute Optogenetic Modulation of Cardiac Twitch Dynamics Explored Through Modeling. J Biomech Eng 2017; 138:2552973. [PMID: 27618140 DOI: 10.1115/1.4034655] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Indexed: 11/08/2022]
Abstract
Optogenetic approaches allow cellular membrane potentials to be perturbed by light. When applied to muscle cells, mechanical events can be controlled through a process that could be termed "optomechanics." Besides functioning as an optical on/off switch, we hypothesized that optomechanical control could include the ability to manipulate the strength and duration of contraction events. To explore this possibility, we constructed an electromechanical model of the human ventricular cardiomyocyte while adding a representation of channelrhodopsin-2 (ChR2), a light-activated channel commonly used in optogenetics. Two hybrid stimulus protocols were developed that combined light-based stimuli with traditional electrical current (all-or-none) excitation. The first protocol involved delivery of a subthreshold optical stimulus followed 50-90 ms later by an electrical stimulus. The result was a graded inhibition of peak cellular twitch force in concert with a prolongation of the intracellular Ca2+ transient. The second protocol was comprised of an electrical stimulus followed by a long light pulse (250-350 ms) that acted to prolong the cardiac action potential (AP). This created a pulse duration-dependent prolongation of the intracellular Ca2+ transient that in turn altered the rate of muscle relaxation without changing peak twitch force. These results illustrate the feasibility of acute, optomechanical manipulation of cardiomyocyte contraction and suggest that this approach could be used to probe the dynamic behavior of the cardiac sarcomere without altering its intrinsic properties. Other experimentally meaningful stimulus protocols could be designed by making use of the optomechanical cardiomyocyte model presented here.
Collapse
Affiliation(s)
- Yasser Aboelkassem
- Institute for Computational Medicine, Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD 21218 e-mail:
| | - Stuart G Campbell
- Department of Biomedical Engineering, Yale University, New Haven, CT 06511 e-mail:
| |
Collapse
|
43
|
Mechano-electrical feedback in the clinical setting: Current perspectives. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2017; 130:365-375. [DOI: 10.1016/j.pbiomolbio.2017.06.001] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Revised: 06/01/2017] [Accepted: 06/02/2017] [Indexed: 12/13/2022]
|
44
|
Gao H, Mangion K, Carrick D, Husmeier D, Luo X, Berry C. Estimating prognosis in patients with acute myocardial infarction using personalized computational heart models. Sci Rep 2017; 7:13527. [PMID: 29051544 PMCID: PMC5648923 DOI: 10.1038/s41598-017-13635-2] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2017] [Accepted: 09/12/2017] [Indexed: 02/06/2023] Open
Abstract
Biomechanical computational models have potential prognostic utility in patients after an acute ST-segment-elevation myocardial infarction (STEMI). In a proof-of-concept study, we defined two groups (1) an acute STEMI group (n = 6, 83% male, age 54 ± 12 years) complicated by left ventricular (LV) systolic dysfunction; (2) an age- and sex- matched hyper-control group (n = 6, 83% male, age 46 ± 14 years), no prior history of cardiovascular disease and normal systolic blood pressure (SBP < 130 mmHg). Cardiac MRI was performed in the patients (2 days & 6 months post-STEMI) and the volunteers, and biomechanical heart models were synthesized for each subject. The candidate parameters included normalized active tension (AT norm) and active tension at the resting sarcomere length (T req, reflecting required contractility). Myocardial contractility was inversely determined from personalized heart models by matching CMR-imaged LV dynamics. Compared with controls, patients with recent STEMI exhibited increased LV wall active tension when normalized by SBP. We observed a linear relationship between T req 2 days post-MI and global longitudinal strain 6 months later (r = 0.86; p = 0.03). T req may be associated with changes in LV function in the longer term in STEMI patients complicated by LV dysfunction. Further studies seem warranted.
Collapse
Affiliation(s)
- Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Kenneth Mangion
- British Heart Foundation, Glasgow Cardiovascular Research Centre, Institute of Cardiovascular and Medical Science, University of Glasgow, Glasgow, UK
- Golden Jubilee National Hospital, Clydebank, UK
| | - David Carrick
- British Heart Foundation, Glasgow Cardiovascular Research Centre, Institute of Cardiovascular and Medical Science, University of Glasgow, Glasgow, UK
- Golden Jubilee National Hospital, Clydebank, UK
| | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Xiaoyu Luo
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Colin Berry
- British Heart Foundation, Glasgow Cardiovascular Research Centre, Institute of Cardiovascular and Medical Science, University of Glasgow, Glasgow, UK.
- Golden Jubilee National Hospital, Clydebank, UK.
| |
Collapse
|
45
|
Weise LD, Panfilov AV. Mechanism for Mechanical Wave Break in the Heart Muscle. PHYSICAL REVIEW LETTERS 2017; 119:108101. [PMID: 28949179 DOI: 10.1103/physrevlett.119.108101] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Indexed: 06/07/2023]
Abstract
Using a reaction-diffusion-mechanics model we identify a mechanism for mechanical wave break in the heart muscle. For a wide range of strengths and durations an external mechanical load causes wave front dissipation leading to formation and breakup of spiral waves. We explain the mechanism, and discuss under which conditions it can cause or abolish cardiac arrhythmias.
Collapse
Affiliation(s)
- L D Weise
- Department of Physics and Astronomy, Ghent University, Krijgslaan 281, S9, Ghent 9000, Belgium
- Theoretical Biology, Utrecht University, Padualaan 8, Utrecht 3584 CH, Netherlands
| | - A V Panfilov
- Department of Physics and Astronomy, Ghent University, Krijgslaan 281, S9, Ghent 9000, Belgium
| |
Collapse
|
46
|
Colli Franzone P, Pavarino LF, Scacchi S. Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study. CHAOS (WOODBURY, N.Y.) 2017; 27:093905. [PMID: 28964121 DOI: 10.1063/1.4999465] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
In this work, we investigate the influence of cardiac tissue deformation on re-entrant wave dynamics. We have developed a 3D strongly coupled electro-mechanical Bidomain model posed on an ideal monoventricular geometry, including fiber direction anisotropy and stretch-activated currents (SACs). The cardiac mechanical deformation influences the bioelectrical activity with two main mechanical feedback: (a) the geometric feedback (GEF) due to the presence of the deformation gradient in the diffusion coefficients and in a convective term depending on the deformation rate and (b) the mechano-electric feedback (MEF) due to SACs. Here, we investigate the relative contribution of these two factors with respect to scroll wave stability. We extend the previous works [Keldermann et al., Am. J. Physiol. Heart Circ. Physiol. 299, H134-H143 (2010) and Hu et al., PLoS One 8(4), e60287 (2013)] that were based on the Monodomain model and a simple non-selective linear SAC, while here we consider the full Bidomain model and both selective and non-selective components of SACs. Our simulation results show that the stability of cardiac scroll waves is influenced by MEF, which in case of low reversal potential of non-selective SACs might be responsible for the onset of ventricular fibrillation; GEF increases the scroll wave meandering but does not determine the scroll wave stability.
Collapse
Affiliation(s)
- P Colli Franzone
- Dipartimento di Matematica, Università di Pavia, Via Ferrata 1, 27100 Pavia, Italy
| | - L F Pavarino
- Dipartimento di Matematica, Università di Pavia, Via Ferrata 1, 27100 Pavia, Italy
| | - S Scacchi
- Dipartimento di Matematica, Università di Milano, Via Saldini 50, 20133 Milano, Italy
| |
Collapse
|
47
|
Gao H, Feng L, Qi N, Berry C, Griffith BE, Luo X. A coupled mitral valve-left ventricle model with fluid-structure interaction. Med Eng Phys 2017; 47:128-136. [PMID: 28751011 PMCID: PMC6779302 DOI: 10.1016/j.medengphy.2017.06.042] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2017] [Revised: 06/13/2017] [Accepted: 06/24/2017] [Indexed: 12/16/2022]
Abstract
Understanding the interaction between the valves and walls of the heart is important in assessing and subsequently treating heart dysfunction. This study presents an integrated model of the mitral valve (MV) coupled to the left ventricle (LV), with the geometry derived from in vivo clinical magnetic resonance images. Numerical simulations using this coupled MV-LV model are developed using an immersed boundary/finite element method. The model incorporates detailed valvular features, left ventricular contraction, nonlinear soft tissue mechanics, and fluid-mediated interactions between the MV and LV wall. We use the model to simulate cardiac function from diastole to systole. Numerically predicted LV pump function agrees well with in vivo data of the imaged healthy volunteer, including the peak aortic flow rate, the systolic ejection duration, and the LV ejection fraction. In vivo MV dynamics are qualitatively captured. We further demonstrate that the diastolic filling pressure increases significantly with impaired myocardial active relaxation to maintain a normal cardiac output. This is consistent with clinical observations. The coupled model has the potential to advance our fundamental knowledge of mechanisms underlying MV-LV interaction, and help in risk stratification and optimisation of therapies for heart diseases.
Collapse
Affiliation(s)
- Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK.
| | - Liuyang Feng
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Nan Qi
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Colin Berry
- Institute of Cardiovascular and Medical Science, University of Glasgow, Glasgow, UK
| | - Boyce E Griffith
- Departments of Mathematics and Biomedical Engineering and McAllister Heart Institute, University of North Carolina, Chapel Hill, NC, USA
| | - Xiaoyu Luo
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| |
Collapse
|
48
|
Cai L, Wang Y, Gao H, Li Y, Luo X. A mathematical model for active contraction in healthy and failing myocytes and left ventricles. PLoS One 2017; 12:e0174834. [PMID: 28406991 PMCID: PMC5391010 DOI: 10.1371/journal.pone.0174834] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2016] [Accepted: 03/15/2017] [Indexed: 02/07/2023] Open
Abstract
Cardiovascular disease is one of the leading causes of death worldwide, in particular myocardial dysfunction, which may lead to heart failure eventually. Understanding the electro-mechanics of the heart will help in developing more effective clinical treatments. In this paper, we present a multi-scale electro-mechanics model of the left ventricle (LV). The Holzapfel-Ogden constitutive law was used to describe the passive myocardial response in tissue level, a modified Grandi-Pasqualini-Bers model was adopted to model calcium dynamics in individual myocytes, and the active tension was described using the Niederer-Hunter-Smith myofilament model. We first studied the electro-mechanics coupling in a single myocyte in the healthy and diseased left ventricle, and then the single cell model was embedded in a dynamic LV model to investigate the compensation mechanism of LV pump function due to myocardial dysfunction caused by abnormality in cellular calcium dynamics. The multi-scale LV model was solved using an in-house developed hybrid immersed boundary method with finite element extension. The predictions of the healthy LV model agreed well with the clinical measurements and other studies, and likewise, the results in the failing states were also consistent with clinical observations. In particular, we found that a low level of intracellular Ca2+ transient in myocytes can result in LV pump function failure even with increased myocardial contractility, decreased systolic blood pressure, and increased diastolic filling pressure, even though they will increase LV stroke volume. Our work suggested that treatments targeted at increased contractility and lowering the systolic blood pressure alone are not sufficient in preventing LV pump dysfunction, restoring a balanced physiological Ca2+ handling mechanism is necessary.
Collapse
Affiliation(s)
- Li Cai
- NPU-UoG International Cooperative Lab for Computation & Application in Cardiology, Northwestern Polytechnical University, Xi’an, Shanxi Province, China
| | - Yongheng Wang
- NPU-UoG International Cooperative Lab for Computation & Application in Cardiology, Northwestern Polytechnical University, Xi’an, Shanxi Province, China
| | - Hao Gao
- School of Mathematics and Statistics, University of Glasgow, Glasgow, Scotland, United Kingdom
| | - Yiqiang Li
- NPU-UoG International Cooperative Lab for Computation & Application in Cardiology, Northwestern Polytechnical University, Xi’an, Shanxi Province, China
| | - Xiaoyu Luo
- School of Mathematics and Statistics, University of Glasgow, Glasgow, Scotland, United Kingdom
| |
Collapse
|
49
|
Land S, Park-Holohan SJ, Smith NP, Dos Remedios CG, Kentish JC, Niederer SA. A model of cardiac contraction based on novel measurements of tension development in human cardiomyocytes. J Mol Cell Cardiol 2017; 106:68-83. [PMID: 28392437 DOI: 10.1016/j.yjmcc.2017.03.008] [Citation(s) in RCA: 60] [Impact Index Per Article: 8.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/14/2016] [Revised: 01/12/2017] [Accepted: 03/31/2017] [Indexed: 11/18/2022]
Abstract
Experimental data from human cardiac myocytes at body temperature is crucial for a quantitative understanding of clinically relevant cardiac function and development of whole-organ computational models. However, such experimental data is currently very limited. Specifically, important measurements to characterize changes in tension development in human cardiomyocytes that occur with perturbations in cell length are not available. To address this deficiency, in this study we present an experimental data set collected from skinned human cardiac myocytes, including the passive and viscoelastic properties of isolated myocytes, the steady-state force calcium relationship at different sarcomere lengths, and changes in tension following a rapid increase or decrease in length, and after constant velocity shortening. This data set is, to our knowledge, the first characterization of length and velocity-dependence of tension generation in human skinned cardiac myocytes at body temperature. We use this data to develop a computational model of contraction and passive viscoelasticity in human myocytes. Our model includes troponin C kinetics, tropomyosin kinetics, a three-state crossbridge model that accounts for the distortion of crossbridges, and the cellular viscoelastic response. Each component is parametrized using our experimental data collected in human cardiomyocytes at body temperature. Furthermore we are able to confirm that properties of length-dependent activation at 37°C are similar to other species, with a shift in calcium sensitivity and increase in maximum tension. We revise our model of tension generation in the skinned isolated myocyte to replicate reported tension traces generated in intact muscle during isometric tension, to provide a model of human tension generation for multi-scale simulations. This process requires changes to calcium sensitivity, cooperativity, and crossbridge transition rates. We apply this model within multi-scale simulations of biventricular cardiac function and further refine the parametrization within the whole organ context, based on obtaining a healthy ejection fraction. This process reveals that crossbridge cycling rates differ between skinned myocytes and intact myocytes.
Collapse
Affiliation(s)
- Sander Land
- Department of Biomedical Engineering, King's College London, UK.
| | - So-Jin Park-Holohan
- Cardiovascular Division, King's College London British Heart Foundation Centre of Research Excellence, UK
| | - Nicolas P Smith
- Department of Engineering Science, University of Auckland, New Zealand
| | | | - Jonathan C Kentish
- Cardiovascular Division, King's College London British Heart Foundation Centre of Research Excellence, UK
| | | |
Collapse
|
50
|
A simple model of cardiac muscle for multiscale simulation: Passive mechanics, crossbridge kinetics and calcium regulation. J Theor Biol 2017; 420:105-116. [PMID: 28223172 DOI: 10.1016/j.jtbi.2017.02.021] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2016] [Revised: 02/09/2017] [Accepted: 02/16/2017] [Indexed: 11/22/2022]
Abstract
A simple model of cardiac muscle was designed for multiscale simulation of heart mechanics. Relaxed cardiac muscle was described as a transversally isotropic hyperelastic material. Active tension caused by actin-myosin crossbridges depends on the ensemble averaged strain of myosin heads bound to actin. Calcium activation was modeled by Ca2+ binding to troponin-C. To account for the dependence of troponin affinity for Ca2+ on myosin heads strongly bound to actin, the kinetics of troponin binding to Ca2+ in the overlap zone of the thin and thick filaments and outside it were separated. The changes in the length of these zones during muscle shortening or lengthening were accounted for explicitly. Simplified version of the model contains only 5 ordinary differential equations (ODE). Model parameters were estimated from a limited set of experiments with skeletal and cardiac muscle. Simulations have shown that model reproduces qualitatively a number of experimental observations: steady-state force-velocity and stiffness-velocity relations; mechanical responses to step changes in muscle length or load; steep Ca2+-tension relationship and its dependence on sarcomere length tension (the Frank-Starling mechanism); tension, shortening and Ca2+-transients in twitch isometric and isotonic contractions, tension development and redevelopment upon fast change in Ca2+ concentration or muscle release followed by re-stretch. We believe that the model can be effectively used for modeling contraction and relaxation of the heart.
Collapse
|