1
|
Alrabghi G, Liu Y, Hu W, Hancox JC, Zhang H. Human atrial fibrillation and genetic defects in transient outward currents: mechanistic insights from multi-scale computational models. Philos Trans R Soc Lond B Biol Sci 2023; 378:20220166. [PMID: 37122220 PMCID: PMC10150223 DOI: 10.1098/rstb.2022.0166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023] Open
Abstract
Previous studies have linked dysfunctional Ito arising from mutations to KCND3-encoded Kv4.3 and KCND2-encoded Kv4.2 to atrial fibrillation. Using computational models, this study aimed to investigate the mechanisms underlying pro-arrhythmic effects of the gain-of-function Kv4.3 (T361S, A545P) and Kv4.2 (S447R) mutations. Wild-type and mutant Ito formulations were developed from and validated against experimental data and incorporated into the Colman et al. model of human atrial cells. Single-cell models were incorporated into one- (1D) and two-dimensional (2D) models of atrial tissue, and a three-dimensional (3D) realistic model of the human atria. The three gain-of-function mutations had similar, albeit quantitatively different, effects: shortening of the action potential duration; lowering the plateau membrane potential, abbreviating the effective refractory period (ERP) and the wavelength (WL) of atrial excitation at the tissue level. Restitution curves for the WL, the ERP and the conduction velocity were leftward shifted, facilitating the conduction of atrial excitation waves at high excitation rates. The mutations also increased lifespan and stationarity of re-entry in both 2D and 3D simulations, which further highlighted a mutation-induced increase in spatial dispersion of repolarization. Collectively, these changes account for pro-arrhythmic effects of these Kv4.3 and Kv4.2 mutations in facilitating AF. This article is part of the theme issue 'The heartbeat: its molecular basis and physiological mechanisms'.
Collapse
Affiliation(s)
- Ghadah Alrabghi
- Biological Physics Group, Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
- Department of Physics, Faculty of Science, University of Jeddah, 21959 Jeddah, Saudi Arabia
| | - Yizhou Liu
- Biological Physics Group, Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
| | - Wei Hu
- Biological Physics Group, Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
| | - Jules C Hancox
- Biological Physics Group, Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
- School of Physiology, Pharmacology and Neuroscience, Medical Sciences Building, University Walk, Bristol BS8 1TD, UK
| | - Henggui Zhang
- Biological Physics Group, Department of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK
- Key Laboratory of Medical Electrophysiology of Ministry of Education and Medical Electrophysiological Key Laboratory of Sichuan Province, Institute of Cardiovascular Research, Southwest Medical University, 646099 Luzhou, People's Republic of China
| |
Collapse
|
2
|
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
|
3
|
Pargaei M, Kumar BVR, Pavarino LF, Scacchi S. Cardiac electro-mechanical activity in a deforming human cardiac tissue: modeling, existence-uniqueness, finite element computation and application to multiple ischemic disease. J Math Biol 2022; 84:17. [PMID: 35142929 DOI: 10.1007/s00285-022-01717-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Revised: 12/11/2021] [Accepted: 01/10/2022] [Indexed: 11/24/2022]
Abstract
In this study, the cardiac electro-mechanical model in a deforming domain is taken with the addition of mechanical feedback and stretch-activated channel current coupled with the ten Tusscher human ventricular cell level model that results in a coupled PDE-ODE system. The existence and uniqueness of such a coupled system in a deforming domain is proved. At first, the existence of a solution is proved in the deformed domain. The local existence of the solution is proved using the regularization and the Faedo-Galerkin technique. Then, the global existence is proved using the energy estimates in appropriate Banach spaces, Gronwall lemma, and the compactness procedure. The existence of the solution in an undeformed domain is proved using the lower semi-continuity of the norms. Uniqueness is proved using Young's inequality, Gronwall lemma, and the Cauchy-Schwartz inequality. For the application purpose, this model is applied to understand the electro-mechanical activity in ischemic cardiac tissue. It also takes care of the development of active tension, conductive, convective, and ionic feedback. The Second Piola-Kirchoff stress tensor arising in Lagrangian mapping between reference and moving frames is taken as a combination of active, passive, and volumetric components. We investigated the effect of varying strength of hyperkalemia and hypoxia, in the ischemic subregions of human cardiac tissue with local multiple ischemic subregions, on the electro-mechanical activity of healthy and ischemic zones. This system is solved numerically using the [Formula: see text] finite element method in space and the implicit-explicit Euler method in time. Discontinuities arising with the modeled multiple ischemic regions are treated to the desired order of accuracy by a simple regularization technique using the interpolating polynomials. We examined the cardiac electro-mechanical activity for several cases in multiple hyperkalemic and hypoxic human cardiac tissue. We concluded that local multiple ischemic subregions severely affect the cardiac electro-mechanical activity more, in terms of action potential (v) and mechanical parameters, intracellular calcium ion concentration [Formula: see text], active tension ([Formula: see text]), stretch ([Formula: see text]) and stretch rate ([Formula: see text]), of a healthy cell in its vicinity, compared to a single Hyperkalemic or Hypoxic subregion. The four moderate hypoxically generated ischemic subregions affect the waveform of the stretch along the fiber and the stretch rate more than a single severe ischemic subregion.
Collapse
Affiliation(s)
- Meena Pargaei
- Department of Mathematics and Statistics, Indian Institute of Technology, Kanpur, India.,Govt. Post Graduate College, Champawat, Uttarakhand, India
| | - B V Rathish Kumar
- Department of Mathematics and Statistics, Indian Institute of Technology, Kanpur, India
| | - Luca F Pavarino
- Department of Mathematics, University of Pavia, Pavia, Italy
| | - Simone Scacchi
- Department of Mathematics, University of Milan, Milan, Italy
| |
Collapse
|
4
|
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
|
5
|
Physicochemical Characterization of PHBV Nanoparticles Functionalized with Multiple Bioactives Designed to be Theranostics for Lung Cancer. J CLUST SCI 2020. [DOI: 10.1007/s10876-020-01912-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
6
|
Kara V, Ni H, Perez Alday EA, Zhang H. ECG Imaging to Detect the Site of Ventricular Ischemia Using Torso Electrodes: A Computational Study. Front Physiol 2019; 10:50. [PMID: 30804799 PMCID: PMC6378918 DOI: 10.3389/fphys.2019.00050] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2018] [Accepted: 01/17/2019] [Indexed: 12/02/2022] Open
Abstract
Electrocardiography provides some information useful for ischemic diagnosis. However, more recently there has been substantial growth in the area of ECG imaging, which by solving the inverse problem of electrocardiography aims to produce high-resolution mapping of the electrical and magnetic dynamics of the heart. Most inverse studies use the full resolution of the body surface potential (BSP) to reconstruct the epicardial potentials, however using a limited number of torso electrodes to interpolate the BSP is more clinically relevant and has an important effect on the reconstruction which must be quantified. A circular ischemic lesion on the right ventricle lateral wall 27 mm in radius is reconstructed using three Tikhonov methods along with 6 different electrode configurations ranging from 32 leads to 1,024 leads. The 2nd order Tikhonov solution performed the most accurately (~80% lesion identified) followed by the 1st (~50% lesion identified) and then the 0 order Tikhonov solution performed the worst with a maximum of ~30% lesion identified regardless of how many leads were used. With an increasing number of leads the solution produces less error, and the error becomes more localised around the lesion for all three regularisation methods. In noisy conditions, the relative performance gap of the 1st and 2nd order Tikhonov solutions was reduced, and determining an accurate regularisation parameter became relatively more difficult. Lesions located on the left ventricle walls were also able to be identified but comparatively to the right ventricle lateral wall performed marginally worse with lesions located on the interventricular septum being able to be indicated by the reconstructions but not successfully identified against the error. The quality of reconstruction was found to decrease as the lesion radius decreased, with a lesion radius of <20 mm becoming difficult to correctly identify against the error even when using >512 torso electrodes.
Collapse
Affiliation(s)
- Vinay Kara
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester, Manchester, United Kingdom
| | - Haibo Ni
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester, Manchester, United Kingdom.,Department of Pharmacology, The University of California, Davis, Davis, CA, United States
| | - Erick Andres Perez Alday
- Division of Cardiovascular Medicine, Oregon Health and Science University, Portland, OR, United States
| | - Henggui Zhang
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester, Manchester, United Kingdom.,School of Computer Science and Technology, Harbin Institute of Technology, Harbin, China.,China Space Institute of Southern China, Shenzhen, China
| |
Collapse
|
7
|
Shavik SM, Wall ST, Sundnes J, Burkhoff D, Lee LC. Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains. Physiol Rep 2018; 5:5/21/e13392. [PMID: 29122952 PMCID: PMC5688770 DOI: 10.14814/phy2.13392] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2017] [Accepted: 07/23/2017] [Indexed: 12/22/2022] Open
Abstract
Although detailed cell‐based descriptors of cross‐bridge cycling have been applied in finite element (FE) heart models to describe ventricular mechanics, these multiscale models have never been tested rigorously to determine if these descriptors, when scaled up to the organ‐level, are able to reproduce well‐established organ‐level physiological behaviors. To address this void, we here validate a left ventricular (LV) FE model that is driven by a cell‐based cross‐bridge cycling descriptor against key organ‐level heart physiology. The LV FE model was coupled to a closed‐loop lumped parameter circulatory model to simulate different ventricular loading conditions (preload and afterload) and contractilities. We show that our model is able to reproduce a linear end‐systolic pressure volume relationship, a curvilinear end‐diastolic pressure volume relationship and a linear relationship between myocardial oxygen consumption and pressure–volume area. We also show that the validated model can predict realistic LV strain‐time profiles in the longitudinal, circumferential, and radial directions. The predicted strain‐time profiles display key features that are consistent with those measured in humans, such as having similar peak strains, time‐to‐peak‐strain, and a rapid change in strain during atrial contraction at late‐diastole. Our model shows that the myocardial strains are sensitive to not only LV contractility, but also to the LV loading conditions, especially to a change in afterload. This result suggests that caution must be exercised when associating changes in myocardial strain with changes in LV contractility. The methodically validated multiscale model will be used in future studies to understand human heart diseases.
Collapse
Affiliation(s)
| | | | | | - Daniel Burkhoff
- Cardiovascular Research Foundation and Department of Medicine, Columbia University, New York, New York
| | - Lik Chuan Lee
- Department of Mechanical Engineering, Michigan State University, East Lansing, Michigan
| |
Collapse
|
8
|
In-silico investigations of the functional impact of KCNA5 mutations on atrial mechanical dynamics. J Mol Cell Cardiol 2017; 111:86-95. [DOI: 10.1016/j.yjmcc.2017.08.005] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 07/19/2017] [Accepted: 08/04/2017] [Indexed: 12/20/2022]
|
9
|
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
|
10
|
Crowcombe J, Dhillon SS, Hurst RM, Egginton S, Müller F, Sík A, Tarte E. 3D Finite Element Electrical Model of Larval Zebrafish ECG Signals. PLoS One 2016; 11:e0165655. [PMID: 27824910 PMCID: PMC5100939 DOI: 10.1371/journal.pone.0165655] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2016] [Accepted: 10/14/2016] [Indexed: 01/08/2023] Open
Abstract
Assessment of heart function in zebrafish larvae using electrocardiography (ECG) is a potentially useful tool in developing cardiac treatments and the assessment of drug therapies. In order to better understand how a measured ECG waveform is related to the structure of the heart, its position within the larva and the position of the electrodes, a 3D model of a 3 days post fertilisation (dpf) larval zebrafish was developed to simulate cardiac electrical activity and investigate the voltage distribution throughout the body. The geometry consisted of two main components; the zebrafish body was modelled as a homogeneous volume, while the heart was split into five distinct regions (sinoatrial region, atrial wall, atrioventricular band, ventricular wall and heart chambers). Similarly, the electrical model consisted of two parts with the body described by Laplace's equation and the heart using a bidomain ionic model based upon the Fitzhugh-Nagumo equations. Each region of the heart was differentiated by action potential (AP) parameters and activation wave conduction velocities, which were fitted and scaled based on previously published experimental results. ECG measurements in vivo at different electrode recording positions were then compared to the model results. The model was able to simulate action potentials, wave propagation and all the major features (P wave, R wave, T wave) of the ECG, as well as polarity of the peaks observed at each position. This model was based upon our current understanding of the structure of the normal zebrafish larval heart. Further development would enable us to incorporate features associated with the diseased heart and hence assist in the interpretation of larval zebrafish ECGs in these conditions.
Collapse
Affiliation(s)
- James Crowcombe
- School of Engineering, University of Birmingham, Birmingham, United Kingdom
| | - Sundeep Singh Dhillon
- Institute of Cancer and Genomic Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Institute of Clinical Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Rhiannon Mary Hurst
- Institute of Cancer and Genomic Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- Institute of Clinical Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Stuart Egginton
- School of Biomedical Sciences, Faculty of Biological Sciences, University of Leeds, Leeds, United Kingdom
| | - Ferenc Müller
- Institute of Cancer and Genomic Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Attila Sík
- Institute of Clinical Sciences, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
| | - Edward Tarte
- School of Engineering, University of Birmingham, Birmingham, United Kingdom
| |
Collapse
|
11
|
Colli Franzone P, Pavarino LF, Scacchi S. Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study. Math Biosci 2016; 280:71-86. [PMID: 27545966 DOI: 10.1016/j.mbs.2016.08.003] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2016] [Revised: 05/25/2016] [Accepted: 08/04/2016] [Indexed: 10/21/2022]
Abstract
The aim of this work is to investigate, by means of numerical simulations, the influence of myocardial deformation due to muscle contraction and relaxation on the cardiac repolarization process in presence of transmural intrinsic action potential duration (APD) heterogeneities. The three-dimensional electromechanical model considered consists of the following four coupled components: the quasi-static transversely isotropic finite elasticity equations for the deformation of the cardiac tissue; the active tension model for the intracellular calcium dynamics and cross-bridge binding; the anisotropic Bidomain model for the electrical current flow through the deforming cardiac tissue; the membrane model of ventricular myocytes, including stretch-activated channels. The numerical simulations are based on our finite element parallel solver, which employs Multilevel Additive Schwarz preconditioners for the solution of the discretized Bidomain equations and Newton-Krylov methods for the solution of the discretized non-linear finite elasticity equations. Our findings show that: (i) the presence of intrinsic transmural cellular APD heterogeneities is not fully masked by electrotonic current flow or by the presence of the mechanical deformation; (ii) despite the presence of transmural APD heterogeneities, the recovery process follows the activation sequence and there is no significant transmural repolarization gradient; (iii) with or without transmural APD heterogeneities, epicardial electrograms always display the same wave shape and discordance between the polarity of QRS complex and T-wave; (iv) the main effects of the mechanical deformation are an increase of the dispersion of repolarization time and APD, when computed over the total cardiac domain and over the endo- and epicardial surfaces, while there is a slight decrease along the transmural direction.
Collapse
Affiliation(s)
- P Colli Franzone
- Dipartimento di Matematica, Università di Pavia, Via Ferrata 1, Pavia 27100, Italy.
| | - L F Pavarino
- Dipartimento di Matematica, Università di Milano, Via Saldini 50, Milano 20133, Italy.
| | - S Scacchi
- Dipartimento di Matematica, Università di Milano, Via Saldini 50, Milano 20133, Italy.
| |
Collapse
|
12
|
Adeniran I, MacIver DH, Garratt CJ, Ye J, Hancox JC, Zhang H. Effects of Persistent Atrial Fibrillation-Induced Electrical Remodeling on Atrial Electro-Mechanics - Insights from a 3D Model of the Human Atria. PLoS One 2015; 10:e0142397. [PMID: 26606047 PMCID: PMC4659575 DOI: 10.1371/journal.pone.0142397] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2015] [Accepted: 10/21/2015] [Indexed: 11/28/2022] Open
Abstract
Aims Atrial stunning, a loss of atrial mechanical contraction, can occur following a successful cardioversion. It is hypothesized that persistent atrial fibrillation-induced electrical remodeling (AFER) on atrial electrophysiology may be responsible for such impaired atrial mechanics. This simulation study aimed to investigate the effects of AFER on atrial electro-mechanics. Methods and Results A 3D electromechanical model of the human atria was developed to investigate the effects of AFER on atrial electro-mechanics. Simulations were carried out in 3 conditions for 4 states: (i) the control condition, representing the normal tissue (state 1) and the tissue 2–3 months after cardioversion (state 2) when the atrial tissue recovers its electrophysiological properties after completion of reverse electrophysiological remodelling; (ii) AFER-SR condition for AF-remodeled tissue with normal sinus rhythm (SR) (state 3); and (iii) AFER-AF condition for AF-remodeled tissue with re-entrant excitation waves (state 4). Our results indicate that at the cellular level, AFER (states 3 & 4) abbreviated action potentials and reduced the Ca2+ content in the sarcoplasmic reticulum, resulting in a reduced amplitude of the intracellular Ca2+ transient leading to decreased cell active force and cell shortening as compared to the control condition (states 1 & 2). Consequently at the whole organ level, atrial contraction in AFER-SR condition (state 3) was dramatically reduced. In the AFER-AF condition (state 4) atrial contraction was almost abolished. Conclusions This study provides novel insights into understanding atrial electro-mechanics illustrating that AFER impairs atrial contraction due to reduced intracellular Ca2+ transients.
Collapse
Affiliation(s)
- Ismail Adeniran
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
| | - David H. MacIver
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
- Taunton & Somerset Hospital, Somerset, United Kingdom
| | - Clifford J. Garratt
- Manchester Heart Centre, Manchester Royal Infirmary, Manchester, United Kingdom
| | - Jianqiao Ye
- Department of Engineering, Lancaster University, Lancaster, United Kingdom
| | - Jules C. Hancox
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
- School of Physiology and Pharmacology, and Cardiovascular Research Laboratories, University of Bristol, Bristol, United Kingdom
| | - Henggui Zhang
- Biological Physics Group, School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
- * E-mail:
| |
Collapse
|
13
|
Adeniran I, MacIver DH, Hancox JC, Zhang H. Abnormal calcium homeostasis in heart failure with preserved ejection fraction is related to both reduced contractile function and incomplete relaxation: an electromechanically detailed biophysical modeling study. Front Physiol 2015; 6:78. [PMID: 25852567 PMCID: PMC4367530 DOI: 10.3389/fphys.2015.00078] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2014] [Accepted: 02/26/2015] [Indexed: 01/08/2023] Open
Abstract
Heart failure with preserved ejection fraction (HFpEF) accounts for about 50% of heart failure cases. It has features of incomplete relaxation and increased stiffness of the left ventricle. Studies from clinical electrophysiology and animal experiments have found that HFpEF is associated with impaired calcium homeostasis, ion channel remodeling and concentric left ventricle hypertrophy (LVH). However, it is still unclear how the abnormal calcium homeostasis, ion channel and structural remodeling affect the electro-mechanical dynamics of the ventricles. In this study we have developed multiscale models of the human left ventricle from single cells to the 3D organ, which take into consideration HFpEF-induced changes in calcium handling, ion channel remodeling and concentric LVH. Our simulation results suggest that at the cellular level, HFpEF reduces the systolic calcium level resulting in a reduced systolic contractile force, but elevates the diastolic calcium level resulting in an abnormal residual diastolic force. In our simulations, these abnormal electro-mechanical features of the ventricular cells became more pronounced with the increase of the heart rate. However, at the 3D organ level, the ejection fraction of the left ventricle was maintained due to the concentric LVH. The simulation results of this study mirror clinically observed features of HFpEF and provide new insights toward the understanding of the cellular bases of impaired cardiac electromechanical functions in heart failure.
Collapse
Affiliation(s)
- Ismail Adeniran
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester Manchester, UK
| | - David H MacIver
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester Manchester, UK ; Department of Cardiology, Taunton and Somerset Hospital Musgrove Park, Taunton, UK
| | - Jules C Hancox
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester Manchester, UK ; School of Physiology and Pharmacology and Cardiovascular Research Laboratories, School of Medical Sciences Bristol, UK
| | - Henggui Zhang
- Biological Physics Group, School of Physics and Astronomy, The University of Manchester Manchester, UK ; School of Computer Science and Technology, Harbin Institute of Technology Harbin, China
| |
Collapse
|