Human atrial fibrillation and genetic defects in transient outward currents: mechanistic insights from multi-scale computational models

Previous studies have linked dysfunctional I_to 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 I_to 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'.

Computational cardiac physiology for new modelers: Origins, foundations, and future

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.

Cardiac electro-mechanical activity in a deforming human cardiac tissue: modeling, existence-uniqueness, finite element computation and application to multiple ischemic disease

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.

Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time

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 AT_ell and AT_ecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (AT_v ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by AT_ta and AT_eff , 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 AT_v , 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 AT_ell and AT_ecc 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.

ECG Imaging to Detect the Site of Ventricular Ischemia Using Torso Electrodes: A Computational Study

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.

Organ-level validation of a cross-bridge cycling descriptor in a left ventricular finite element model: effects of ventricular loading on myocardial strains

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.

Effects of mechanical feedback on the stability of cardiac scroll waves: A bidomain electro-mechanical simulation study

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.

3D Finite Element Electrical Model of Larval Zebrafish ECG Signals

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.

Joint influence of transmural heterogeneities and wall deformation on cardiac bioelectrical activity: A simulation study

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.

Effects of Persistent Atrial Fibrillation-Induced Electrical Remodeling on Atrial Electro-Mechanics - Insights from a 3D Model of the Human Atria

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 Ca²⁺ content in the sarcoplasmic reticulum, resulting in a reduced amplitude of the intracellular Ca²⁺ 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 Ca²⁺ transients.

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

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.