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.

Introduction

Heart failure (HF) is often categorized into two major types (Borlaug and Paulus, 2010; MacIver, 2010a; Phan et al., 2012; Liu et al., 2013; Zouein et al., 2013): HF with reduced ejection fraction (HFrEF) and HF with preserved ejection fraction (HFpEF). The main distinguishing criterion is an arbitrary cut-off value for the left ventricular ejection fraction of >50% for HFpEF (Vasan and Levy, 2000; Zile et al., 2001; Yturralde and Gaasch, 2005; Asrar Ul Haq et al., 2014) and ≤ 50% for HFrEF (Borlaug and Paulus, 2010; Phan et al., 2012; Liu et al., 2013; Zouein et al., 2013). HFpEF and HFrEF are also commonly referred to as diastolic and systolic HF respectively and also share a variety of abnormalities (Borlaug and Paulus, 2010; Soma, 2011; MacIver and Dayer, 2012; Asrar Ul Haq et al., 2014). Currently, the prevalence of HFpEF is ~50% (Konstantinou et al., 2013; Liu et al., 2013; Asrar Ul Haq et al., 2014), but it is predicted to become the dominant form of HF within the next decade (Liu et al., 2013; Asrar Ul Haq et al., 2014).

The leading cause of HFpEF is hypertension and the dominant pathophysiological mechanism is thought to involve impaired relaxation of the left ventricle (LV) (Liu et al., 2013; Asrar Ul Haq et al., 2014). Other distinguishing features cited are increased LV stiffness and elevated LV end-diastolic pressures (Liu et al., 2013; Asrar Ul Haq et al., 2014). In addition, abnormalities of systolic shortening are common with reduced global longitudinal strain, strain rate, reduced midwall fractional shortening, reduced systolic annular motion, increased isovolumetric contraction time and reduced systolic longitudinal shortening velocities in HFpEF (Sanderson, 2007; Wang et al., 2008; Kono et al., 2009; MacIver, 2010a). HFpEF is also characterized by dyspnoea, fluid retention, exercise intolerance, coronary artery disease and atrial fibrillation (Redfield et al., 2003; Bhatia et al., 2006; Owan et al., 2006; Konstantinou et al., 2013; Liu et al., 2013). It has a greater prevalence in older people and females (Liu et al., 2013). Patients also exhibit other comorbidities such as diabetes, obesity, peptic ulcer disease, cancer and chronic obstructive pulmonary disease (Liu et al., 2013).

As relaxation is an active energy-consuming process and dependent on intracellular Ca²⁺ homeostasis (Robertson et al., 1982; Ebashi, 1984; Barry and Bridge, 1993; Bers, 2001; Konstantinou et al., 2013), impaired relaxation in HFpEF may reflect abnormal intracellular Ca²⁺ homeostasis. Indeed, experimental data from clinical electrophysiology and animal model studies have shown that HFpEF is associated with abnormal Ca²⁺ handling, including an increased SR Ca²⁺ leak current (I_leak) and a decreased Ca²⁺ release from ryanodine receptors (Borbély et al., 2005; Selby et al., 2011; Zile and Gaasch, 2011; Trenor et al., 2012; Gomez et al., 2014). In addition, some ionic currents responsible for generating cardiac action potentials are also remodeled, including those carried by the late-sodium channel, the transient outward K⁺ channel, the inward rectifier K⁺ channel, the Na⁺/K⁺ pump (I_NaK), the background Ca²⁺ channel and the Na⁺/Ca²⁺ exchanger (Borbély et al., 2005; Selby et al., 2011; Zile and Gaasch, 2011; Gomez et al., 2014; Trenor et al., 2012). Structurally, HFpEF is also associated with concentric LV hypertrophy with high LV mass/volume ratio, cardiomyocyte hypertrophy and interstitial fibrosis. This is different to HFrEF, where structural remodeling is characterized by progressive ventricular dilation, eccentric LV remodeling, low LV mass/volume ratio, cardiomyocyte loss along with replacement fibrosis (Borlaug and Paulus, 2010; Konstantinou et al., 2013; Asrar Ul Haq et al., 2014).

The functional impacts of HFpEF-related abnormal Ca²⁺ homeostasis, ion channel and structural remodeling on cardiac electro-mechanics are unclear. Thus, the present study was conducted in order: (1) to develop a novel HFpEF electromechanical single cell model based on extant experimental data; (2) to investigate the cellular mechanisms influencing myocardial calcium homeostasis in HFpEF using an electromechanical single cell model; and (3) to evaluate the functional impacts of impaired Ca²⁺ handling, ion channel remodeling and degrees of concentric LV hypertrophy on the electro-mechanical activity of the heart.

Materials and Methods

Electromechanical Single Cell Model

For electrophysiology (EP), we utilized the O'Hara-Rudy (ORd) human ventricular single cell model (O'Hara et al., 2011), which was developed from undiseased human ventricular data and recapitulates human ventricular cell electrical and membrane channel properties, as well as the transmural heterogeneity of ventricular action potential (AP) across the ventricular wall (O'Hara et al., 2011). The ORd model also reproduces Ca²⁺ vs. voltage-dependent inactivation of L-type Ca²⁺ current and Ca²⁺/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca²⁺ cycling. For simulating cellular mechanics properties, we used the Tran et al. myofilament (MM) model (Tran et al., 2010). This model was chosen as it is an extension of the well-established cross-bridge cycling model of cardiac muscle contraction model of Rice et al. (2008). In addition to its ability to replicate a wide range of experimental data including steady-state force-sarcomere length (F-SL), force-calcium and sarcomere length-calcium relationships (Rice et al., 2008; Tran et al., 2010), it also reproduces many of the observed effects of MgATP, MgADP, Pi, and H⁺ on force development.

The intracellular calcium concentration from the EP model was used as the coupling link to the MM model. [Ca²⁺]_i produced as dynamic output from the EP model during the AP served as input to the MM model from which the amount of Ca²⁺ bound to troponin is calculated. The formulation of the myoplasmic Ca²⁺ concentration in the human ventricular myocyte electromechanical cell model is:

$\begin{array}{cc}\begin{array}{l}\frac{d{[C{a}^{2+}]}_i}{dt}={\beta }_{Cai}·(-(I_{pCa}+I_{Cab}-2·I_{NaCa,i})·\frac{A_{cap}}{2·F·v_{myo}}\\ -J_{up}·\frac{v_{nsr}}{v_{myo}}+J_{diff,Ca}·\frac{v_{ss}}{v_{myo}}-\frac{J_{Trop}}{1000})\end{array}& \left(1\right)\end{array}$

where β_Cai is the buffer factor for [Ca²⁺]_i, I_pCa (μA/μF) is the sarcolemmal Ca²⁺ pump current, I_Cab (μA/μF) is the Ca²⁺ background current, I_NaCa,i (μA/μF) is the myoplasmic component of Na⁺/Ca²⁺ exchange current, A_cap (cm²) is capacitive area, F (coul/mol) is the Faraday constant, v_myo (μL) is the volume of the myoplasmic compartment, v_nsr (μL)is the volume of the network sarcoplasmic reticulum compartment, v_ss (μL) is the volume of the subspace compartment, J_up (mM/ms) is the total Ca²⁺ uptake flux, via SERCA pump from myoplasm to the network sarcoplasmic reticulum, J_diff,Ca (mM/ms) is the flux of the diffusion of Ca²⁺ from the subspace to the myoplasm and J_Trop (μM/ms) is the flux of Ca²⁺ binding to troponin calculated via the MM model. β_Cai is formulated as:

$\begin{array}{cc}{\beta }_{Cai}=\frac{1}{1+\frac{[CMDN]·K_{m,CMDN}}{{\left(K_{m,CMDN}+{[C{a}^{2+}]}_i\right)}^2}}& \left(2\right)\end{array}$

where [CMDN] is the calmodulin Ca²⁺ buffer in the myoplasm and K_m,CMDN is the half-saturation concentration of calmodulin.

HFpEF Electromechanical Single Cell Model

To develop the human electromechanical single cell model of a HFpEF cardiomyocyte, we first modified parameters of the ORd model based on HF experimental data following the work of Gomez et al. (2014), Trenor et al. (2012) on human HF. This resulted in a generic human HF model. We then made alterations based on the cellular and molecular properties of HFpEF (Zile and Gaasch, 2011) to obtain a biophysically detailed model of HFpEF. For the myofilament changes relative to control, we made modifications based on experimental data on HFpEF patients following the work of Borbély et al. (2005), Selby et al. (2011). Table 1 summarizes the extensive modifications to the control ORd cardiomyocyte to produce the HFpEF cardiomyocyte model.

TABLE 1

Table 1. Changes in original ORd model to simulate HFpEF and HFrEF.

Parameter Sensitivity Analysis for Influential Cellular Processes on Relaxation in HFpEF

To our knowledge, only the work by Zile and Gaasch (2011) provides any data on the differences in the cellular and molecular processes that influence Ca²⁺ homeostasis in patients with systolic and diastolic HF. We based our HFpEF model on their work (Zile and Gaasch, 2011) as given in section HFpEF Electromechanical Single Cell Model. In simulations, the only difference between the HFrEF and HFpEF models is a reduction in NCX in the HFpEF condition (see Table 1) as compared to the HFrEF condition. In order to determine whether a reduction in NCX is responsible for the poor end-diastolic relaxation in HFpEF, we performed a parameter-dependent sensitivity analysis of the HFpEF model to NCX, by varying NCX from 70% (HFpEF) to 175% (HFrEF).

Protocols

The pacing protocols used to evaluate Ca²⁺ homeostasis were as follows:

Post-Rest Contractions

We used a post-rest contraction (PRC) protocol to evaluate sarcoplasmic reticulum (SR) Ca²⁺ content, retention, release, reuptake, and leak. After pacing the single cell for 10 min at 1 Hz to allow steady-conditions to be attained, resting intervals of 1, 2, 3, 5, and 10 s were introduced. The resting periods were then followed by a single stimulus. The varying developed indices such as active tension are a reflection of SR Ca²⁺ release.

Tissue Mechanics Model

We modeled cardiac tissue mechanics within the theoretical framework of nonlinear elasticity (Marsden and Hughes, 1994; Holzapfel, 2000) as an inhomogeneous, anisotropic, nearly incompressible nonlinear material similar to previous studies (Costa et al., 2001; Whiteley et al., 2007; Niederer and Smith, 2008; Pathmanathan and Whiteley, 2009; Adeniran et al., 2013a,b). We used a two-field variational principle with the deformation u and the hydrostatic pressure p as the two fields (Bonet and Wood, 2008; Adeniran et al., 2013a,b; Le Tallec). p is utilized as the Lagrange multiplier to enforce the near incompressibility constraint. Thus, the total potential energy function Π for the mechanics problem is formulated as:

$\begin{array}{cc}\Pi (u,p)={\Pi }_{\text{int}}(u,p)+{\Pi }_{ext}(u)& \left(3\right)\end{array}$

where Π_int (u, p) is the internal potential energy or total strain energy of the body and Π_ext (u) is the external potential energy or potential energy of the external loading of the body. As in previous studies (Niederer and Smith, 2008; Keldermann et al., 2009; Pathmanathan and Whiteley, 2009; Adeniran et al., 2013a,b), in the absence of body forces, and assuming that the body is always in instantaneous equilibrium and no inertia effects, the coordinates of the deformed body satisfies the steady-state equilibrium equation with near incompressibility enforced.

The values that minimize the total potential energy in Equation (3) are obtained by searching for its critical points in suitable admissible displacement and pressure spaces $\widehat{U}$ and $\widehat{P}$. The corresponding Euler-Lagrange equations resulting from Equation (3) lead to solving the problem (Braess and Ming, 2005; Auricchio et al., 2010; Boffi et al., 2013; Le Tallec):

Find (u, p) in $\widehat{U}$ × $\widehat{P}$ such that:

$\begin{array}{cc}\begin{array}{l}{\displaystyle {\int }_{\Omega }\frac{\partial \widehat{W}}{dF}}(x,Id+\nabla u):\nabla vdx-{\displaystyle {\int }_{\Omega }p}\frac{\partial \text{det}}{dF}(x,Id+\nabla u):\nabla vdx\\ ={\displaystyle {\int }_{\partial \Omega }g}.vds\forall v,u\in \widehat{U}\\ {\displaystyle {\int }_{\Omega }q}\left[\text{det}(Id+\nabla u)-1\right]dx=0\forall v,u\in \widehat{P}\end{array}& \left(4\right)\end{array}$

where $\widehat{U}$ and $\widehat{P}$ are the admissible variation spaces for the displacements and the pressures, respectively. F = Id + ∇u is the deformation gradient, v is a test function and $\widehat{W}$ is the material stored energy function and corresponds to the density of elastic energy locally stored in the body during the deformation.

The deformation gradient F is a tensor that maps elements from the undeformed configuration to the deformed configuration (Marsden and Hughes, 1994; Holzapfel, 2000). Following Cherubini et al. (2008), Ambrosi et al. (2011), we multiplicatively decompose F into a microscopic (active) component and a macroscopic elastic (passive) component:

$\begin{array}{cc}F=F_e{F}_0& \left(5\right)\end{array}$

The active component F_o measures the length change of the tissue due to muscle contraction while the passive component F_e accounts for the passive mechanical response of the tissue and possible tension due to external loads.

With the vector fields f, s, and n denoting the unique direction of the fibers, sheets and sheet-normals in the undeformed state of the LV, the microscopic active component of the deformation tensor F takes the form:

$\begin{array}{cc}{F}_o=I+{\gamma }_{f}f\otimes f+{\gamma }_{s}s\otimes s+\gamma n\otimes n& \left(6\right)\end{array}$

where I is the identity tensor and γ_x is a scalar field that represents the intensity of the contraction, i.e., the active strain in the appropriate direction. γ_f is defined as:

$\begin{array}{cc}{\gamma }_f=\frac{SL-S{L}_0}{S{L}_0}& \left(7\right)\end{array}$

where SL is the new sarcomeric length and SL₀ is the resting sarcomere length of the electromechanical single cell. Following the work of (Rossi et al., 2014)

${\gamma }_n=K{\gamma }_f,{\gamma }_s=\frac{1}{(1+{\gamma }_f)(1+{\gamma }_n)}-1$

Thus, γ > 0 denotes elongation, and γ < 0 denotes contraction. The parameter K is the link between the microscopic and the macroscopic active deformations (Bogaert and Rademakers, 2001; Rossi et al., 2014) and in our simulations, we used K = 4, according to experimental observations (Rademakers et al., 1994; Rossi et al., 2014).

The elastic component F_e is formulated as:

$\begin{array}{cc}{F}_e=F{F}_0^{-1}& \left(8\right)\end{array}$

and the corresponding Right Cauchy-Green strain tensor is:

$\begin{array}{cc}{C}_e=F_e^T{F}_e& \left(9\right)\end{array}$

The associated Green-Lagrange strain tensor is:

$\begin{array}{cc}{E}_e=\frac{1}{2}(C_e-I)& \left(10\right)\end{array}$

To characterize the constitutive behavior of cardiac tissue, we used a mixed formulation of the compressible neo-Hookean strain energy function W (Auricchio et al., 2010):

$\begin{array}{cc}\widehat{W}=\frac{\mu }{2}[I:\widehat{C}-d]-\mu \text{ln}\widehat{J}+p\text{ln}\widehat{J}-\frac{p^2}{2\lambda }& \left(11\right)\end{array}$

where μ and λ are positive constants, “:” represents the usual inner product for second-order tensors, $\widehat{C}$ is the right Cauchy-Green deformation tensor, d is the spatial dimension and $\widehat{J}$ = det$\widehat{F}$ is the Jacobian of the deformation gradient.

Tissue Electrophysiology Model

The monodomain representation (Colli Franzone et al., 2005; Keener and Sneyd, 2008; Adeniran et al., 2013a,b) of cardiac tissue was used for the electrophysiology model:

$\begin{array}{cc}{C}_m\frac{dv}{dt}=-(I_{ion}+I_{stim})+\nabla ·(D\nabla V)& \left(12\right)\end{array}$

where C_m is the cell capacitance per unit surface area, V is the membrane potential, I_ion is the sum of all transmembrane ionic currents from the electromechanics single cell model, I_stim is an externally applied stimulus and D is the anisotropic diffusion tensor given by:

$\begin{array}{cc}D={\sigma }_f(f\otimes f)+{\sigma }_s(s\otimes s)+{\sigma }_n(n\otimes n)& \left(13\right)\end{array}$

σ_f, σ_f, and σ_f are the intracellular conductivities in the fiber, sheet and cross-sheet directions respectively. These were set to 3.0, 0.1, and 0.31525 ms·mm⁻¹ respectively. These gave a conduction velocity of 65 cm·s⁻¹ in the fiber direction along multiple cells, which is close to the value 70 cm·s⁻¹ observed in the fiber direction in human myocardium (Taggart et al., 2000). To guard against any drift in the steady state values of the ion concentrations in the model, the electromechanical single cell models described in sections Electromechanical Single Cell Model to HFpEF Electromechanical Single Cell Model were pre-paced for a 1000 beats before being incorporated into the tissue model.

Geometry and Meshes

For unbiased comparison between increasing severity of concentric LV hypertrophy in HFpEF, we used truncated ellipsoids for the 3D simulations (Figure 1A). Each LV geometry was segmented into distinct endocardial (60%), mid-myocardial (30%), and epicardial (10%) regions. The chosen cell proportion in each region reflects experimental data for cells spanning the left ventricular wall of the human heart (Drouin et al., 1995). There exists some controversy on the existence and functional role of MCELLs in the human heart (Wilson et al., 2011). However, we included MCELLs in our model on the basis of studies that give evidence of their existence in cells isolated from the right ventricle of patients with heart failure (Li et al., 1998) and in a perfused piece of left ventricular wall (“wedge preparation”) (Drouin et al., 1995). The conditional activation sites were projected from those of a healthy 34-year old male determined empirically across the ventricle wall and were validated by reproducing the activation sequence and QRS complex in the measured 64-channel ECG (Keller et al., 2009) of that person. Figures 1B–D shows the anisotropic fiber fields reconstructed by the algorithm in Rossi et al. (2014) for the fiber (Figure 1B), sheet (Figure 1C), and cross-sheet (Figure 1D) directions.

FIGURE 1

Figure 1. Geometry and fiber directions in 3D left ventricle with varied wall thickness. (A) Control and HFpEF concentric hypertrophic geometries—NORMAL (green, 9 mm), MILD (red, 12 mm), MODERATE (blue, 15 mm) and SEVERE (yellow, 18 mm). (B) Fiber direction. (C) Sheet direction. (D) Sheet-normal or cross-sheet direction.

Solving the Electromechanics Problem

The electromechanics problem consists of two sub-problems: the electrophysiology problem and the mechanics problem. The electrophysiology problem (Equation 14) was solved with a Strang splitting method (Sundnes et al., 2005) ensuring that the solution is second-order accurate. It was discretized in time using the Crank-Nicholson method (Burnett, 1987), which is also second-order accurate and discretized in space with Finite Elements (Ciarlet, 2002; Braess, 2007; Brenner and Scott, 2010). I_ion in Equation (14) represents the single cell electromechanics model from which the active strain input to the 3D mechanics model for contraction is obtained. The system of ordinary differential equations (ODE) composing I_ion was solved with a combination of the Rush-Larsen scheme (Rush and Larsen, 1978) and the CVODE solver (Cohen and Hindmarsh, 1996; Alan and Hindmarsh, 2005).

The mechanics problem (Equation 2) was also solved using the Finite element Method using the automated scientific computing library, FEniCS (Logg et al., 2012). The resulting nonlinear system of equations was solved iteratively using Newton's method to determine the equilibrium configuration of the system. Over a typical finite element domain, P₂ elements (Braess, 2007; Brenner and Scott, 2010; Ern and Guermond, 2010; Boffi et al., 2013) were used to discretize the displacement variable u, while the pressure variable p was discretized with P₁ elements (Braess, 2007; Brenner and Scott, 2010; Ern and Guermond, 2010; Boffi et al., 2013). This P₂–P₁ mixed finite element has been proven to ensure stability (Chamberland et al., 2010; Haga et al., 2012; Logg et al., 2012) and an optimal convergence rate (Hughes, 2000; Chamberland et al., 2010; Ern and Guermond, 2010).

The algorithm for solving the full electromechanics problem is as follows:

(1) While time < tend:

(a) Solve the electrophysiology problem for Δt_mechanics = 1 ms with active strain γ as output (Δt_{electrophysiology} = 0.01 ms).

(b) Project γ from the electrophysiology mesh onto the mechanics mesh.

(2) Solve the mechanics problem with γ as input.

Results

Model Validation

To validate the HFpEF model, we computed a steady-state force-calcium (F-pCa) relation for a sarcomere length (SL) of 2.2 μm for comparison with prior experimental data (Borbély et al., 2005). Results are shown in Figure 2. The model reproduced the differences in total and passive forces between control and HFpEF, which matched experimental data (inset).

FIGURE 2

Figure 2. Force-pCa relationship. Simulated Force-pCa relation of control and HFpEF. Relative force is normalized to maximum value in control. Inset: Experimental force-pCa relation from human patients adopted from Borbély et al. (2005).

Functional Consequences of HFpEF Transmurally in the LV

Figure 3 shows the electromechanical consequences of HFpEF in a free-running cell from each of the epicardium (EPI) (Figures 3Ai–Aiv), mid-myocardium (MCELL) (Figures 3Bi–Biv) and endocardium (ENDO) (Figures 3Ci–Civ). HFpEF increased the action potential duration (APD₉₀) in each of the EPI (Figure 3Bi), MCELL (Figure 3Bii) and ENDO (Figure 3Bi) cell types. APD₉₀ was increased from 233 to 262 ms in EPI, 357 to 439 ms in MCELL and 274 to 348 ms in ENDO. These differences in APD₉₀ across the ventricular wall are suggestive of possible T-wave changes and altered dispersion of repolarization in HFpEF. The cytosolic Ca²⁺ concentration [Ca²⁺]_i was reduced in amplitude in all the cell types (Figures 3Aii–Cii), which led to a decrease in the SL shortening (Figures 3Aiii–Ciii) and consequently, a reduction in the active force (Figures 3Aiv–Civ). Of note is that the diastolic [Ca²⁺]_i level was increased in HFpEF compared to control (Figures 3Aii–Cii) leading to incomplete relaxation (Figures 3Aiv–Civ). This electromechanical model predicted a combination of reduced contractile force and stress generation and incomplete relaxation as the pathophysiological mechanism of HFpEF (Zile et al., 2004; Borlaug and Paulus, 2010; Phan et al., 2012; Konstantinou et al., 2013; Liu et al., 2013; Zouein et al., 2013; Asrar Ul Haq et al., 2014).

FIGURE 3

Figure 3. Single cell simulations of HFpEF. (Ai–Ci) Control (black) and HFpEF (green) action potentials in the EPI (Ai), MCELL (Bi), and ENDO (Ci) cell models. (Aii–Cii) Control (black) and HFpEF (green) cytosolic Ca²⁺ concentration in the EPI (Aii), MCELL (Bii), and ENDO (Cii) cell models. (Aiii–Ciii) Control (black) and HFpEF (green) sarcomere length (SL) in the EPI (Aiii), MCELL (Biii), and ENDO (Ciii) cell models. (Aiv–Civ) Control (black) and HFpEF (green) active force in the EPI (Aiv), MCELL (Biv), and ENDO (Civ) cell models. Values are normalized to Control maximum active force for each cell type.

Functional Consequences of HFpEF on Ionic Currents Compared with HFrEF

Figure 4 shows action potentials in control, HFrEF and HFpEF conditions (Figure 4A) computed from the EPI cell model, together with the corresponding ionic currents of I_CaL (Figure 4B), I_NaCa (Figure 4C), the SR content (Figure 4D), I_pCa (Figure 4E), intracellular Na⁺ ([Na⁺]_i; Figure 4F) and K⁺ ([K⁺]_i; Figure 4G) concentrations and I_Na,K. In comparison to the control condition, HFpEF and HFrEF increased the APD, which can be attributable to the increased inward Na⁺ current (I_NaL) and decreased outward K⁺ current (I_K1 and I_NaK; see Table 1) (Figure 4H) (Glitsch, 2001; Workman et al., 2003; Bueno-Orovio et al., 2014). Accompanying these changes were increased [Na⁺]_i (Figure 4F) and reduced [K⁺]_i (Figure 4G). The difference in the amplitude of I_CaL among control, HFpEF and HFrEF was small, but its duration was longer in HFpEF and longest in HFrEF (Figure 4B), possibly due to a secondary effect of a prolonged APD. In the HFpEF condition, a decreased I_NaCa was observed, partially due to reduced Na⁺-Ca²⁺ exchanger (NCX) (see Table 1) and partially due to the reduction in the Na⁺-K⁺ pump activity as it indirectly regulated Ca²⁺ extrusion by the NCX (Barry et al., 1985; Bueno-Orovio et al., 2014). This resulted in a reduced sarcoplasmic reticulum (SR) content in HFpEF (Figure 4D), and consequently leading to a smaller SR Ca²⁺ release as demonstrated by the reduced amplitude of the SR content (Figure 4D). The activity of the sarcolemmal Ca²⁺ pump was increased during the diastolic period, but reduced in the systolic period (Figure 4E). These simulation results provided a cellular basis for the abnormal Ca²⁺ handling in HFpEF. Results from the MCELL and ENDO cell models showed similar behavior (these are shown in the Supplementary Material).

FIGURE 4

Figure 4. Effects of HFpEF and HFrEF on underlying ion channel currents, concentrations and force generation. (A) Control (black), HFpEF (green), and HFrEF (red) action potentials. (B) I_CaL current profile in control (black), HFpEF (green), and HFrEF (red). (C) I_NaCa current profile in control (black), HFpEF (green), and HFrEF (red). (D) SR Ca²⁺ content profile in control (black), HFpEF (green), and HFrEF (red). (E) I_pCa current profile in control (black), HFpEF (green), and HFrEF (red). (F) [Na]_i time course in control (black), HFpEF (green), and HFrEF (red). (G) [K]_i time course in control (black), HFpEF (green), and HFrEF (red). (H) I_NaK current profile in control (black) HFpEF (green), and HFrEF (red). (I) J_up (Ca²⁺ uptake via SERCA pump) profile in control (black), HFpEF (green), and HFrEF (red). (J) Ca²⁺ concentration in Control (black) HFpEF (green), and HFrEF (red) cytosolic. (K) Sarcomere length (SL) in control (black), HFpEF (green), and HFrEF (red). (L) Active force in control (black), HFpEF (green), and HFrEF (red). Values are normalized to Control maximum active force.

Unlike in HFpEF, there was an increased I_NaCa (Figure 4C) due to the HFrEF model formulation involved a 175% NCX increase (Table 1). Figure 4C shows that NCX removed an excessive amount of Ca²⁺ from the cell in its forward mode compared to control, leading to a deficit in the SR Ca²⁺ content (Figure 4D), a greater reduction in [Ca²⁺]_i compared to HFpEF (Figure 4J) and consequently a lower active force (Figure 4L).

Post-Rest Contraction Properties of HFpEF

Further simulations were performed to investigate the rate-dependent impact of HEpEF on Ca²⁺ handling. Results are shown in Figure 5, which shows the outcome of the PRC protocol at a pacing rate of 1 Hz (Figures 5Ai–Av) and 2 Hz (Figures 5Bi–Bv). Diastolic [Ca²⁺]_i level in HFpEF increased by ~75% relative to control at all resting intervals (Figures 5Ai,Bi). This led to an increase in resting tension (Figures 5Aiii,Biii). Though there was a significant increase in diastolic Ca²⁺]_i level, the SR content was lower in HFpEF than in control (Figures 5Aii,Bii) as was peak systolic tension (Figures 5Aiv,Biv). This is counter-intuitive as with an increased diastolic concentration of Ca²⁺ in the cytosol, one would expect greater Ca²⁺ sequestration into the SR, and therefore an increased Ca²⁺ content leading to an increased Ca²⁺ release from the SR, resulting in a greater Ca²⁺ transient amplitude and a greater systolic tension. The reduction in the peak systolic Ca²⁺ level and the corresponding tension and inefficient SR Ca²⁺ activity was due to excessive leak of Ca²⁺ from the SR (Figures 5Av,Bv). These results were more pronounced at a pacing rate of 2 Hz (Figures 5Bi–Bv) because of the shorter duration between beats allowing less recovery time for Ca²⁺ cycling processes.

FIGURE 5

Figure 5. Ca²⁺ handling and post-rest properties in HFpEF. (Ai,Bi) Diastolic Ca²⁺ level in HFpEF relative to control at 1 Hz (Ai) and 2 Hz (Bi) pacing rates. (Aii,Bii) SR Ca²⁺ content level in HFpEF relative to control at 1 Hz (Aii) and 2 Hz (Bii) pacing rates. (Aiii,Biii) Resting tension in HFpEF relative to control at 1 Hz (Aiii) and 2 Hz (Biii) pacing rates. (Aiv,Biv) Peak systolic tension in HFpEF relative to control at 1 Hz (Aiv) and 2 Hz (Biv) pacing rates. (Av,Bv) SR Ca²⁺ Ca leak in HFpEF relative to control at 1 Hz (Av) and 2 Hz (Bv) pacing rates.

Sensitivity of Diastolic Relaxation to NCX

In Figure 6, effects of a systematic change of NCX in HFpEF model from 70% (HEpEF) to 175% (HFrEF) on the action potentials (Figure 6A), I_NaCa (Figure 6B), [Ca²⁺]_i (Figure 6C) and active force (Figure 6D) are shown. Diastolic relaxation was impaired when NCX was at 70% of the control value (e.g., for HFpEF), but gradually improved with increasing NCX activity to 100 and 150% of the control value. It became normal at 175% (e.g., for HFrEF; Figure 6D) of the control value. Peak relative force during systole occurred later and with lower amplitude with increasing NCX activity.

FIGURE 6

Figure 6. HFpEF Model sensitivity to increasing NCX activity and its influence on incomplete relaxation in the cellular model. Simulation results were compared between control (in absence of ionic current remodeling) and HFpEF condition. HFpEF simulations were performed with parameters as listed in Table 1, but with NCX activity changing from 70% (HFpEF condition) to 100, 150, and 175% (HFrEF condition) of the control value. (A) Action potential (B) I_NaCa current profile. (C) Ca²⁺ concentration. (D) Active force. Values are normalized to maximum active force in control. (Inset: magnified diastolic phase).

3D Electromechanical Consequences of HFpEF in Hypertrophic Geometries

Figure 7 shows the results of incorporating the cellular HFpEF electromechanical models into three-dimensional truncated ellipsoid representations of the LV. Effects of varying degrees of left ventricular hypertrophy on the ejection fraction are also shown. The LV was at rest before activation (0 ms). At about 200 ms, the LV was completely activated and contraction was underway in the NORMAL, MILD, MODERATE, and SEVERE conditions of hypertrophy. By 300 ms, repolarization had commenced and the LV in all conditions was undergoing relaxation. At 700 ms, repolarization was completed in all the conditions; however, relaxation was still on going in the hypertrophic cases but was complete in the NORMAL condition. The LVEF in NORMAL, MILD, MODERATE, and SEVERE was 58, 57, 63, and 71% respectively showing that LVEF was increased with increasing the end-diastolic wall thickness. These simulation results can account for incomplete relaxation and preserved LVEF at the 3D organ level, whilst at the cellular level, the activation force and cell sarcolemmal shortening are dramatically impaired in the HFpEF condition. This is attributable to the hypertrophied ventricle wall.

FIGURE 7

Figure 7. Effects of HFpEF on 3D electro-mechanics. Electrical wave propagation and mechanical contraction at 0, 200, 300, and 700 ms in NORMAL, MILD, MODERATE, and SEVERE HFpEF hypertrophic cases. (Far right) Ejection fraction in NORMAL, MILD, MODERATE, and SEVERE HFpEF hypertrophic cases.

Stress Distribution in HFpEF

Effects of LV wall thickness on the spatial distribution of stress was investigated. Figure 8 shows a 4-chamber view of the magnitude of the stress distribution across the LV in control and varied hypertrophic conditions (Figures 8A–G). At 0 ms, the stress magnitude across the ventricle was low in all the conditions as excitation had yet to commence (Figure 8A). By 100 ms, when the LV was electrically activated, there was developed stress in all cases with the greatest stress in apex of the MODERATE LV (Figure 8B). At 200 ms (Figure 8C), there was considerable stress in the LV in all conditions with the greatest stress intensity in SEVERE. The hypertrophic LV also had greater stress intensity in the LV apex and epicardium compared to NORMAL. The situation was similar at 300 ms (Figure 8D) and 400 ms (Figure 8E) except that there was progressive relaxation and the stress intensity was less than at 200 ms. By 600 ms (Figure 8F), stress in most of the NORMAL LV had reduced considerably while the hypertrophic LVs still showed ~25% stress. At 700 ms, stress in the NORMAL LV was negligible but still about ~20% in the hypertrophic cases. These simulation results showed that increased wall thickness led to increased tissue stress though, at the cellular level, the active force was reduced during the time course of action potentials in HEpEF condition.

FIGURE 8

Figure 8. Effects of HFpEF on stress magnitude distribution. Stress magnitude distribution at 0 ms (A), 100 ms (B), 200 ms (C), 300 ms (D), 400 ms (E) 600 ms (F), and 700 ms (G) in NORMAL, MILD, MODERATE, and SEVERE hypertrophic HFpEF cases.

Strain Distribution in HFpEF

Effects of wall thickness on the spatial distribution of strain was also analyzed. Figure 9 also shows a 4-chamber view of the magnitude of the strain distribution across the LV in control and hypertrophic conditions during the time course cardiac excitation (Figures 9A–G). Strain developed from the rest state (Figure 9A) to ~50% at 100 ms (Figure 9B) in all cases with the smallest strain magnitude in SEVERE. By 200 ms (Figure 9C), strain magnitude had reached a distribution of between 50 and 100% in all cases with SEVERE having the highest strain intensity. This was concentrated in the epicardial to mid-wall regions and the middle segment of the LV. The situation was similar at 300 ms (Figure 9D) and 400 ms (Figure 9E) except that the strain intensity was reduced as relaxation was underway. The strain magnitude reduced further at 600 ms (Figure 9F) to ~10–15% in NORMAL, ~25% in MILD, ~25–30% in MODERATE and ~30–35% in SEVERE. The strain was mainly concentrated in the LV apex in all cases. By 700 ms (Figure 9G), strain magnitude was ~5% in NORMAL, 10–25% in all the hypertrophic cases. These simulation results showed an increased residual of strain in cardiac tissue due to elevated diastolic Ca²⁺ level, which was consistent with observed incomplete relaxation of the LV in HEpEF condition.

FIGURE 9

Figure 9. Effects of HFpEF on strain magnitude distribution. Strain magnitude distribution at 0 ms (A), 100 ms (B), 200 ms (C), 300 ms (D), 400 ms (E) 600 ms (F), and 700 ms (G) in NORMAL, MILD, MODERATE, and SEVERE hypertrophic HFpEF cases.

Discussion

Summary of Major Findings

In this study we have developed, for the first time, a family of multilevel models for the electro-mechanics of the left ventricle in the setting of HFpEF, at the cellular and 3D organ levels. These models incorporated detailed HFpEF-related ion channel remodeling and impaired Ca²⁺ homeostasis at the cellular level, and concentric hypertrophy of the left ventricle wall at the organ level. Our major findings are: (i) with impaired Ca²⁺ handling and ion channel remodeling in HFpEF, the action potential duration of ventricular cells are prolonged, together with an elevated diastolic Ca²⁺ concentration, but a decreased systolic Ca²⁺ level. Such an elevated diastolic Ca²⁺ concentration provides a cellular basis for incomplete ventricular relaxation at the organ level; (ii) at the cellular level the active force and sarcomere length shortening is reduced during the time course of action potentials in HFpEF. However, at the organ level, tissue stress and strain is increased due to the increased wall thickness of concentric LV hypertrophy; (iii) the impaired Ca²⁺ homeostasis becomes more pronounced at high stimulation rates; and (iv) reduction of I_NaCa in the HEpEF model is the most influential factor on impaired relaxation dynamics. Collectively, these simulations predict the key features of HFpEF observed clinically, and also provide insights for understanding the cellular and tissue bases of impaired electro-mechanics of the heart in HFpEF.

Cellular Basis of Impaired Cardiac Electro-Mechanics

The results of our study suggest that at the cellular level the observed impaired cardiac electro-mechanics (such as reduced cell length shortening and active force) are attributable to alterations in cellular Ca²⁺ homeostasis and action potentials. In simulations, we observed an elevated diastolic intracellular Ca²⁺ concentration, but a reduced systolic Ca²⁺ concentration. The elevated diastolic Ca²⁺ concentration is due to an increased Ca²⁺ leak from the SR (see Table 1), and the reduced systolic Ca²⁺ concentration can be explained by a reduced Ca²⁺ release from the SR as a consequence of a decreased SR content (Figure 4D). The observed APD prolongation is attributable to the augmented late-Na⁺ current, reduced potassium channel currents and I_NaK as seen in the heart failure condition (see Table 1). These results were more pronounced at a pacing rate of 2 Hz (Figures 5Bi–Bv) because of the shorter duration between beats allowing less recovery time for Ca²⁺ cycling processes. This is related to changes in mechanical and relaxation restitution, which correlate physiologically to the recovery kinetics of Ca²⁺ release mechanisms and sequestration capacity of the SR (Franz et al., 1983; Burkhoff et al., 1984; Prabhu and Freeman, 1995; Zaugg et al., 1995; Kjørstad et al., 2007).

In simulations, we also observed an increased intracellular Na⁺ concentration and a reduced intracellular K⁺ concentration (Figures 4F,G), which may be attributable to the augmented I_NaL and reduced I_NaK (see Table 1). The altered Na⁺ and K⁺ homeostasis may also impair Ca²⁺ homeostasis. During the time course of an action potential, there was a reduced I_NaK, which also contributed partially to an increase of [Na⁺]_i and a decrease in [K⁺]_i within the cytosol. NCX extrudes Ca²⁺ from the cytoplasm and imports Na⁺ into it. However, with the build-up of Na⁺ in the cytoplasm as a result of the reduced activity of the Na⁺-K⁺ pump, the forward mode activity of NCX was reduced, but its reverse mode activity was enhanced. This led to an increase in Na⁺ extrusion coupled to Ca²⁺ import (see Figure 4C) leading to an increased Ca²⁺ concentration at diastolic phase but a reduced Ca²⁺ concentration at systolic phase. Consequently, this led to a decreased cell length shortening and active force in the systolic phase.

In addition, our simulation results showed that the SR Ca²⁺ content was reduced compared to control condition, though the Ca²⁺ release from the SR was compromised. This is attributable to an excessive Ca²⁺ leak from the SR. It was shown that SR Ca²⁺ leak was ~100% greater in HFpEF as compared to control (Figure 5). The excessive Ca²⁺ leak from the SR and the reduced activity of NCX produced an increased diastolic Ca²⁺ level, leading to incomplete relaxation and an increase in the active resting tone. This implies that complete relaxation will never be achieved regardless of the duration of diastole.

The simulation results discussed above were based on our HFpEF model developed from work by Zile and Gaasch (2011) and Selby et al. (2011). To our knowledge these are the only two studies that provide any notion of electrophysiological changes in HFpEF relative to control conditions. As the only difference between the HFpEF and the HFrEF models is a 30% NCX reduction in the former compared to control and a 175% increase in the latter (Table 1), we investigated whether changes in NCX activity could convert HFpEF to HFrEF and vice-versa, and therefore affect myocyte relaxation dynamics (Figure 6). This indeed proved to be the case. The results imply that a change in NCX activity is the dominant factor leading to impaired diastolic relaxation in HFpEF.

Effects of Concentric Hypertrophy of the Ventricle Wall on Cardiac Electro-Mechanics

Increased ventricular wall thickness has dramatic impacts on cardiac electro-mechanics. This is illustrated by the effects of varied degrees of hypertrophy on the ejection fraction in HFpEF condition (Figure 9). In simulations, the wall thickness of the LV varied from 9 mm (normal) to 12, 15, and 18 mm, mimicking varying degrees of hypertrophy. Though at the cellular level, cell length shortening and active force was decreased in HFpEF condition, the LVEF was preserved and increased with the increase of LV wall thickness. This agrees with what is seen clinically (Borlaug and Paulus, 2010; Phan et al., 2012; Liu et al., 2013; Zouein et al., 2013) and with previous modeling data (MacIver and Townsend, 2008; MacIver, 2010b, 2011). Our simulations confirm that the preserved LVEF in HFpEF is due to the thicker wall of the ventricle arising from concentric hypertrophy.

Taken together, our simulations suggest a possible pathway and mechanism underlying cardiac dysfunction in HFpEF. Any co-morbidities such as diabetes, hypertension, inflammation and/or hypertrophy may cause the ion channel and myofilament remodeling in myocardial cells. Cellular remodeling results in abnormal Ca²⁺ homeostasis, which in turn lead to abnormalities of contraction and incomplete relaxation of the LV with a persistent active resting tone. These abnormalities combine to increase the length of time to which the myocardium is subjected to stress prolonging systole and resulting in abnormal energy utilization and less efficient ejection of the stroke volume.

Relevance to Previous Studies

Lacombe et al. (2007) investigated the underlying mechanisms of diastolic dysfunction in type 1 diabetic rats. They observed no significant change in I_CaL, a reduction in Ca²⁺ transient amplitude and prolongation in its decay, a reduction in SR Ca²⁺ load and a decrease in the expression of sarco(endo)plasmic reticulum Ca²⁺-ATPase-2a (SERCa) protein levels. They concluded that impairment of Ca²⁺ reuptake during myocyte relaxation contributed to diastolic dysfunction, with preserved global systolic function.

Selby et al. (2011) carried out a study to evaluate tachycardia-induced relaxation abnormalities in myocardium from patients with a normal ejection fraction. They observed incomplete relaxation with increased diastolic tension development at increasing pacing rates, significant resting tone and disproportionately elevated Ca²⁺ loads due to reduced sarcolemmal Ca²⁺ extrusion reserve. However, their patients did not carry a clinical diagnosis of heart failure.

MacIver and Townsend (2008) performed a mathematical study on HFpEF to determine the effect of changes in LV hypertrophy on stroke volume and LVEF. They concluded from their model that the preserved LVEF in HFpEF patients was a result of LV hypertrophy, which amplified absolute radial wall thickening in the setting of reduced long-axis shortening. MacIver also showed, using a simple abstract model, that remodeling was necessary to normalize stroke volume and suggested that regulation of end-diastolic volume was a primary compensatory mechanism in heart failure (MacIver, 2010b). However, their models did not consider the contribution of cardiomyocytes, coupled electrical wave propagation or nonlinear anisotropic cardiac mechanics.

Our simulation results are in agreement with and extend the findings of these previous studies (Lacombe et al., 2007; MacIver and Townsend, 2008; Selby et al., 2011), adding new evidence that impaired Ca²⁺ homeostasis (such as reduced systolic Ca²⁺ concentration and elevated diastolic Ca²⁺ concentration) together with hypertrophied wall underlie the key features of HFpEF with preserved ejection fraction and incomplete end-diastolic relaxation.

Limitations

In addition to acknowledged limitations of both the ORd electrophysiology model (O'Hara et al., 2011) and the (Tran et al., 2010) myofilament model, as experimental data show a reduction in NCX activity in diastolic dysfunction compared to control (Zile and Gaasch, 2011), we made a 30% reduction in NCX activity in our HFpEF model as there was no quantitative data available at the time of this study. In HFpEF patients, collagen production results in interstitial fibrosis (Heymans et al., 2005; Konstantinou et al., 2013), which we considered by reducing the intracellular electrical conductivities by 20% due to a lack of quantitative data. The HFpEF model also relaxes somewhat faster than the experimental data, which would imply that the effects of incomplete relaxation would be even greater than we have shown. As there are no data available, we assumed the same degree of ion channel remodeling for the hypertrophied ventricles and also assumed the same distribution of ENDO, MCELL, and EPI cell types across the ventricular wall. We also did not consider electromechanical feedback between the electrical wave propagation and mechanically contracting ventricles. Finally, the use of a fluid-structure interaction model with the interaction of blood and the myocardial wall to determine pressure boundary conditions would allow a more realistic pressure profile. Whilst it is important that these potential limitations are stated, they do not fundamentally alter the principal conclusions of this study.

Conclusion

We have developed a novel, biophysically detailed model of HFpEF and used it to investigate the cellular mechanisms underlying myocardial Ca²⁺ homeostasis in HFpEF. We observed an elevated diastolic [Ca²⁺]_i level, a reduction in SR Ca²⁺ content and reduced SR Ca²⁺ release, a reduction in SR Ca²⁺ sequestration, an increase in resting tension, incomplete relaxation, reduced systolic stress and prolonged stress and strain durations. These mechanisms suggest that in HFpEF patients, impaired Ca²⁺ handling principally caused by reduction in NCX activity is a dominant abnormality in the condition that explains the mechanisms for impaired cardiac electro-mechanics in HFpEF.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by project grants from Engineering and Physical Science Research Council UK (EP/J00958X/1; EP/I029826/1).

Supplementary Material

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fphys.2015.00078/abstract

References

Adeniran, I., Hancox, J. C., and Zhang, H. (2013a). Effect of cardiac ventricular mechanical contraction on the characteristics of the ECG: a simulation study. J. Biomed. Sci. Eng. 6, 47–60. doi: 10.4236/jbise.2013.612A007

CrossRef Full Text | Google Scholar

Adeniran, I., Hancox, J., and Zhang, H. (2013b). In silico investigation of the short QT syndrome, using human ventricle models incorporating electromechanical coupling. Front. Physiol. 4:166. doi: 10.3389/fphys.2013.00166

CrossRef Full Text | Google Scholar

Alan, C., and Hindmarsh, P. N. B. (2005). SUNDIALS: suite of nonlinear and differential/algebraic equation solvers. ACM Trans. Math. Softw. 31, 363–396. doi: 10.1145/1089014.1089020

CrossRef Full Text | Google Scholar

Ambrosi, D., Arioli, G., Nobile, F., and Quarteroni, A. (2011). Electromechanical coupling in cardiac dynamics: the active strain approach. SIAM J. Appl. Math. 71, 605–621. doi: 10.1137/100788379

CrossRef Full Text | Google Scholar

Asrar Ul Haq, M., Mutha, V., Rudd, N., Hare, D. L., and Wong, C. (2014). Heart failure with preserved ejection fraction—unwinding the diagnosis mystique. Am. J. Cardiovasc. Dis. 4, 100–113.

PubMed Abstract | Full Text | Google Scholar

Auricchio, F., Beirão da Veiga, L., Lovadina, C., and Reali, A. (2010). The importance of the exact satisfaction of the incompressibility constraint in nonlinear elasticity: mixed FEMs versus NURBS-based approximations. Comput. Methods Appl. Mech. Eng. 199, 314–323. doi: 10.1016/j.cma.2008.06.004

CrossRef Full Text | Google Scholar

Barry, W. H., and Bridge, J. H. (1993). Intracellular calcium homeostasis in cardiac myocytes. Circulation 87, 1806–1815. doi: 10.1161/01.CIR.87.6.1806

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Barry, W. H., Hasin, Y., and Smith, T. W. (1985). Sodium pump inhibition, enhanced calcium influx via sodium-calcium exchange, and positive inotropic response in cultured heart cells. Circ. Res. 56, 231–241. doi: 10.1161/01.RES.56.2.231

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Bers, D. (2001). Excitation-Contraction Coupling and Cardiac Contractile Force, 2nd Edn. Dordrecht: Springer.

Google Scholar

Bhatia, R. S., Tu, J. V., Lee, D. S., Austin, P. C., Fang, J., Haouzi, A., et al. (2006). Outcome of heart failure with preserved ejection fraction in a population-based study. N. Engl. J. Med. 355, 260–269. doi: 10.1056/NEJMoa051530

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Boffi, D., Brezzi, F., and Fortin, M. (2013). Mixed Finite Element Methods and Applications, 2013 Edn. New York, NY: Springer.

Bogaert, J., and Rademakers, F. E. (2001). Regional nonuniformity of normal adult human left ventricle. Am. J. Physiol. Heart Circ. Physiol. 280, H610–H620.

PubMed Abstract | Full Text | Google Scholar

Bonet, J., and Wood, R. D. (2008). Nonlinear Continuum Mechanics for Finite Element Analysis, 2nd Edn. New York, NY: Cambridge University Press.

Google Scholar

Borbély, A., van der Velden, J., Papp, Z., Bronzwaer, J. G., Edes, I., Stienen, G. J., et al. (2005). Cardiomyocyte stiffness in diastolic heart failure. Circulation 111, 774–781. doi: 10.1161/01.CIR.0000155257.33485.6D

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Borlaug, B. A., and Paulus, W. J. (2010). Heart failure with preserved ejection fraction: pathophysiology, diagnosis, and treatment. Eur. Heart J. 32, 670–679. doi: 10.1093/eurheartj/ehq426

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Braess, D. (2007). Finite Elements: Theory, Fast Solvers, and Applications in Solid Mechanics, 3rd Edn. New York, NY: Cambridge University Press.

Google Scholar

Braess, D., and Ming, P. (2005). A finite element method for nearly incompressible elasticity problems. Math. Comput. 74, 25–52. doi: 10.1090/S0025-5718-04-01662-X

CrossRef Full Text | Google Scholar

Brenner, S. C., and Scott, R. (2010). The Mathematical Theory of Finite Element Methods, 3rd Edn. 2008. New York, NY: Springer.

Google Scholar

Bueno-Orovio, A., Sánchez, C., Pueyo, E., and Rodriguez, B. (2014). Na/K pump regulation of cardiac repolarization: insights from a systems biology approach. Pflügers. Arch. Eur. J. Physiol. 466, 183–193. doi: 10.1007/s00424-013-1293-1

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Burkhoff, D., Yue, D. T., Franz, M. R., Hunter, W. C., and Sagawa, K. (1984). Mechanical restitution of isolated perfused canine left ventricles. Am. J. Physiol. 246, H8–H16.

PubMed Abstract | Full Text | Google Scholar

Burnett, D. S. (1987). Finite Element Analysis: From Concepts to Applications, 1st Edn. Reading, MA: Addison Wesley.

Chamberland, E., Fortin, A., and Fortin, M. (2010). Comparison of the performance of some finite element discretizations for large deformation elasticity problems. Comput. Struct. 88, 664–673. doi: 10.1016/j.compstruc.2010.02.007

CrossRef Full Text | Google Scholar

Cherubini, C., Filippi, S., Nardinocchi, P., and Teresi, L. (2008). An electromechanical model of cardiac tissue: constitutive issues and electrophysiological effects. Prog. Biophys. Mol. Biol. 97, 562–573. doi: 10.1016/j.pbiomolbio.2008.02.001

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Ciarlet, P. G. (2002). The Finite Element Method for Elliptic Problems, 2nd Edn. New York, NY: Society for Industrial and Applied Mathematics.

Google Scholar

Cohen, S., and Hindmarsh, A. C. (1996). “Cvode, a stiff/nonstiff ode solver in C,” in C. Computers in Physics, ed L. M. Holmes (New York, NY: American Institute of Physics Inc.), 138–143.

Colli Franzone, P., Pavarino, L. F., and Taccardi, B. (2005). Simulating patterns of excitation, repolarization and action potential duration with cardiac Bidomain and Monodomain models. Math. Biosci. 197, 35–66. doi: 10.1016/j.mbs.2005.04.003

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Costa, K. D., Holmes, J. W., and McCulloch, A. D. (2001). Modelling cardiac mechanical properties in three dimensions. Philos. Trans. R. Soc. Lond. Ser. Math. Phys. Eng. Sci. 359, 1233–1250. doi: 10.1098/rsta.2001.0828

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Drouin, E., Charpentier, F., Gauthier, C., Laurent, K., and Le Marec, H. (1995). Electrophysiologic characteristics of cells spanning the left ventricular wall of human heart: evidence for presence of M cells. J. Am. Coll. Cardiol. 26, 185–192. doi: 10.1016/0735-1097(95)00167-X

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Ebashi, S. (1984). Ca²⁺ and the contractile proteins. J. Mol. Cell. Cardiol. 16, 129–136. doi: 10.1016/S0022-2828(84)80701-4

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Ern, A., and Guermond, J.-L. (2010). Theory and Practice of Finite Elements, 1st Edn. 2004. New York, NY: Springer.

Franz, M. R., Schaefer, J., Schöttler, M., Seed, W. A., and Noble, M. I. (1983). Electrical and mechanical restitution of the human heart at different rates of stimulation. Circ. Res. 53, 815–822. doi: 10.1161/01.RES.53.6.815

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Glitsch, H. G. (2001). Electrophysiology of the sodium-potassium-ATPase in cardiac cells. Physiol. Rev. 81, 1791–1826.

PubMed Abstract | Full Text | Google Scholar

Gomez, J. F., Cardona, K., Romero, L., Ferrero, J. M., and Trenor, B. (2014). Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study. PLoS ONE 9:e106602. doi: 10.1371/journal.pone.0106602

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Haga, J. B., Osnes, H., and Langtangen, H. P. (2012). On the causes of pressure oscillations in low-permeable and low-compressible porous media. Int. J. Numer. Anal. Methods Geomech. 36, 1507–1522. doi: 10.1002/nag.1062

CrossRef Full Text | Google Scholar

Heymans, S., Schroen, B., Vermeersch, P., Milting, H., Gao, F., Kassner, A., et al. (2005). Increased cardiac expression of tissue inhibitor of metalloproteinase-1 and tissue inhibitor of metalloproteinase-2 is related to cardiac fibrosis and dysfunction in the chronic pressure-overloaded human heart. Circulation 112, 1136–1144. doi: 10.1161/CIRCULATIONAHA.104.516963

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Holzapfel, G. A. (2000). Nonlinear Solid Mechanics: A Continuum Approach for Engineering, 1st Edn. Chichester: Wiley.

Hughes, T. J. R. (2000). The Finite Element Method: Linear Static and Dynamic Finite Element Analysis. New York, NY: Dover Publications.

Google Scholar

Keener, J., and Sneyd, J. (2008). Mathematical Physiology: II: Systems Physiology, 2nd Edn. New York, NY: Springer.

Google Scholar

Keldermann, R. H., Nash, M. P., and Panfilov, A. V. (2009). Modeling cardiac mechano-electrical feedback using reaction-diffusion-mechanics systems. Phys. Nonlinear Phenom. 238, 1000–1007. doi: 10.1016/j.physd.2008.08.017

CrossRef Full Text | Google Scholar

Keller, D. U. J., Kalayciyan, R., Dössel, O., and Seemann, G. (2009). “Fast creation of endocardial stimulation profiles for the realistic simulation of body surface ECGs,” in IFMBE Proceedings, 145–148. Available online at: http://www.scopus.com/inward/record.url?eid=2-s2.0-77950134835&partnerID=40&md5=2db2b89cc19de5e28ec01ebbf0dffff0 (Accessed July 20, 2011).

Kjørstad, K. E., Nordhaug, D. O., Korvald, C., Müller, S., Steensrud, T., and Myrmel, T. (2007). Mechanical restitution curves: a possible load independent assessment of contractile function. Eur. J. Cardio-Thorac. Surg. 31, 677–684. doi: 10.1016/j.ejcts.2007.01.013

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Kono, M., Kisanuki, A., Takasaki, K., Nakashiki, K., Yuasa, T., Kuwahara, E., et al. (2009). Left ventricular systolic function is abnormal in diastolic heart failure: re-assessment of systolic function using cardiac time interval analysis. J. Cardiol. 53, 437–446. doi: 10.1016/j.jjcc.2009.02.014

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Konstantinou, D. M., Chatzizisis, Y. S., and Giannoglou, G. D. (2013). Pathophysiology-based novel pharmacotherapy for heart failure with preserved ejection fraction. Pharmacol. Ther. 140, 156–166. doi: 10.1016/j.pharmthera.2013.05.012

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Lacombe, V. A., Viatchenko-Karpinski, S., Terentyev, D., Sridhar, A., Emani, S., Bonagura, J. D., et al. (2007). Mechanisms of impaired calcium handling underlying subclinical diastolic dysfunction in diabetes. Am. J. Physiol. Regul. Integr. Comp. Physiol. 293, R1787–R1797. doi: 10.1152/ajpregu.00059.2007

CrossRef Full Text | Google Scholar

Le Tallec, P. (1994). “Numerical methods for nonlinear three-dimensional elasticity,” in Handbook of Numerical Analysis, Vol. III, eds P. G. Ciarlet and J. L. Lions (North-Holland: Amsterdam), 465–622.

Li, G. R., Feng, J., Yue, L., and Carrier, M. (1998). Transmural heterogeneity of action potentials and Ito1 in myocytes isolated from the human right ventricle. Am. J. Physiol. 275, H369–H377.

PubMed Abstract | Full Text | Google Scholar

Liu, Y., Haddad, T., and Dwivedi, G. (2013). Heart failure with preserved ejection fraction: current understanding and emerging concepts. Curr. Opin. Cardiol. 28, 187–196. doi: 10.1097/HCO.0b013e32835c5492

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Logg, A., Mardal, K.-A., and Wells, G. (eds.) (2012). Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book, 2012th Edn. Heidelberg: Springer.

Google Scholar

MacIver, D. H. (2010a). Current controversies in heart failure with a preserved ejection fraction. Future Cardiol. 6, 97–111. doi: 10.2217/fca.09.56

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

MacIver, D. H. (2010b). Is remodeling the dominant compensatory mechanism in both chronic heart failure with preserved and reduced left ventricular ejection fraction? Basic Res. Cardiol. 105, 227–234. doi: 10.1007/s00395-009-0063-x

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

MacIver, D. H. (2011). A new method for quantification of left ventricular systolic function using a corrected ejection fraction. Eur. J. Echocardiogr. 12, 228–234. doi: 10.1093/ejechocard/jeq185

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

MacIver, D. H., and Dayer, M. J. (2012). An alternative approach to understanding the pathophysiological mechanisms of chronic heart failure. Int. J. Cardiol. 154, 102–110. doi: 10.1016/j.ijcard.2011.05.075

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

MacIver, D. H., and Townsend, M. (2008). A novel mechanism of heart failure with normal ejection fraction. Heart Br. Card. Soc. 94, 446–449. doi: 10.1136/hrt.2006.114082

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Marsden, J. E., and Hughes, T. J. R. (1994). Mathematical Foundations of Elasticity. New Jersey: Dover Publications.

Google Scholar

Niederer, S. A., and Smith, N. P. (2008). An improved numerical method for strong coupling of excitation and contraction models in the heart. Prog. Biophys. Mol. Biol. 96, 90–111. doi: 10.1016/j.pbiomolbio.2007.08.001

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

O'Hara, T., Virág, L., Varró, A., and Rudy, Y. (2011). Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput. Biol. 7:e1002061. doi: 10.1371/journal.pcbi.1002061

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Owan, T. E., Hodge, D. O., Herges, R. M., Jacobsen, S. J., Roger, V. L., and Redfield, M. M. (2006). Trends in prevalence and outcome of heart failure with preserved ejection fraction. N. Engl. J. Med. 355, 251–259. doi: 10.1056/NEJMoa052256

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Pathmanathan, P., and Whiteley, J. P. (2009). A numerical method for cardiac mechanoelectric simulations. Ann. Biomed. Eng. 37, 860–873. doi: 10.1007/s10439-009-9663-8

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Phan, T. T., Shivu, G. N., Abozguia, K., Sanderson, J. E., and Frenneaux, M. (2012). The pathophysiology of heart failure with preserved ejection fraction: from molecular mechanisms to exercise haemodynamics. Int. J. Cardiol. 158, 337–343. doi: 10.1016/j.ijcard.2011.06.113

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Prabhu, S. D., and Freeman, G. L. (1995). Effect of tachycardia heart failure on the restitution of left ventricular function in closed-chest dogs. Circulation 91, 176–185. doi: 10.1161/01.CIR.91.1.176

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Rademakers, F. E., Rogers, W. J., Guier, W. H., Hutchins, G. M., Siu, C. O., Weisfeldt, M. L., et al. (1994). Relation of regional cross-fiber shortening to wall thickening in the intact heart. Three-dimensional strain analysis by NMR tagging. Circulation 89, 1174–1182. doi: 10.1161/01.CIR.89.3.1174

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Redfield, M. M., Jacobsen, S. J., Burnett, J. C., Mahoney, D. W., Bailey, K. R., and Rodeheffer, R. J. (2003). Burden of systolic and diastolic ventricular dysfunction in the community: appreciating the scope of the heart failure epidemic. JAMA 289, 194–202. doi: 10.1001/jama.289.2.194

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Rice, J. J., Wang, F., Bers, D. M., and de Tombe, P. P. (2008). Approximate model of cooperative activation and crossbridge cycling in cardiac muscle using ordinary differential equations. Biophys. J. 95, 2368–2390. doi: 10.1529/biophysj.107.119487

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Robertson, S. P., Johnson, J. D., Holroyde, M. J., Kranias, E. G., Potter, J. D., and Solaro, R. J. (1982). The effect of troponin I phosphorylation on the Ca²⁺-binding properties of the Ca²⁺-regulatory site of bovine cardiac troponin. J. Biol. Chem. 257, 260–263.

PubMed Abstract | Full Text | Google Scholar

Rossi, S., Lassila, T., Ruiz-Baier, R., Sequeira, A., and Quarteroni, A. (2014). Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechanics. Eur. J. Mech. ASolids. 48, 129–142. doi: 10.1016/j.euromechsol.2013.10.009

CrossRef Full Text | Google Scholar

Rush, S., and Larsen, H. (1978). A practical algorithm for solving dynamic membrane equations. IEEE Trans. Biomed. Eng. 25, 389–392. doi: 10.1109/TBME.1978.326270

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sanderson, J. E. (2007). Heart failure with a normal ejection fraction. Heart 93, 155–158. doi: 10.1136/hrt.2005.074187

CrossRef Full Text | Google Scholar

Selby, D. E., Palmer, B. M., LeWinter, M. M., and Meyer, M. (2011). Tachycardia-induced diastolic dysfunction and resting tone in myocardium from patients with a normal ejection fraction. J. Am. Coll. Cardiol. 58, 147–154. doi: 10.1016/j.jacc.2010.10.069

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Soma, J. (2011). Heart failure with preserved left ventricular ejection fraction: concepts, misconceptions and future directions. Blood Press. 20, 129–133. doi: 10.3109/08037051.2010.542642

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Sundnes, J., Lines, G. T., and Tveito, A. (2005). An operator splitting method for solving the bidomain equations coupled to a volume conductor model for the torso. Math. Biosci. 194, 233–248. doi: 10.1016/j.mbs.2005.01.001

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Taggart, P., Sutton, P. M., Opthof, T., Coronel, R., Trimlett, R., Pugsley, W., et al. (2000). Inhomogeneous transmural conduction during early ischaemia in patients with coronary artery disease. J. Mol. Cell. Cardiol. 32, 621–630. doi: 10.1006/jmcc.2000.1105

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Tran, K., Smith, N. P., Loiselle, D. S., and Crampin, E. J. (2010). A metabolite-sensitive, thermodynamically constrained model of cardiac cross-bridge cycling: implications for force development during ischemia. Biophys. J. 98, 267–276. doi: 10.1016/j.bpj.2009.10.011

CrossRef Full Text | Google Scholar

Trenor, B., Cardona, K., Gomez, J. F., Rajamani, S., Ferrero, J. M., Belardinelli, L., et al. (2012). Simulation and mechanistic investigation of the arrhythmogenic role of the late sodium current in human heart failure. PLoS ONE 7:e32659. doi: 10.1371/journal.pone.0032659

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Vasan, R. S., and Levy, D. (2000). Defining diastolic heart failure: a call for standardized diagnostic criteria. Circulation 101, 2118–2121. doi: 10.1161/01.CIR.101.17.2118

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Wang, J., Khoury, D. S., Yue, Y., Torre-Amione, G., and Nagueh, S. F. (2008). Preserved left ventricular twist and circumferential deformation, but depressed longitudinal and radial deformation in patients with diastolic heart failure. Eur. Heart J. 29, 1283–1289. doi: 10.1093/eurheartj/ehn141

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Whiteley, J. P., Bishop, M. J., and Gavaghan, D. J. (2007). Soft tissue modelling of cardiac fibres for use in coupled mechano-electric simulations. Bull. Math. Biol. 69, 2199–2225. doi: 10.1007/s11538-007-9213-1

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Wilson, L. D., Jennings, M. M., and Rosenbaum, D. S. (2011). Point: M cells are present in the ventricular myocardium. Heart Rhythm 8, 930–933. doi: 10.1016/j.hrthm.2011.01.026

PubMed Abstract | Full Text | CrossRef Full Text

Workman, A. J., Kane, K. A., and Rankin, A. C. (2003). Characterisation of the Na, K pump current in atrial cells from patients with and without chronic atrial fibrillation. Cardiovasc. Res. 59, 593–602. doi: 10.1016/S0008-6363(03)00466-8

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Yturralde, R. F., and Gaasch, W. H. (2005). Diagnostic criteria for diastolic heart failure. Prog. Cardiovasc. Dis. 47, 314–319. doi: 10.1016/j.pcad.2005.02.007

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zaugg, C. E., Kojima, S., Wu, S. T., Wikman Coffelt, J., Parmley, W. W., and Buser, P. T. (1995). 997-86 intracellular calcium transients underlie mechanical restitution in whole rat hearts. J. Am. Coll. Cardiol. 25, 325A. doi: 10.1016/0735-1097(95)92811-I

CrossRef Full Text | Google Scholar

Zile, M. R., Baicu, C. F., and Gaasch, W. H. (2004). Diastolic heart failure—abnormalities in active relaxation and passive stiffness of the left ventricle. N. Engl. J. Med. 350, 1953–1959. doi: 10.1056/NEJMoa032566

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zile, M. R., and Gaasch, W. H. (2011). Abnormal calcium homeostasis: one mechanism in diastolic heart failure. J. Am. Coll. Cardiol. 58, 155–157. doi: 10.1016/j.jacc.2010.10.068

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zile, M. R., Gaasch, W. H., Carroll, J. D., Feldman, M. D., Aurigemma, G. P., Schaer, G. L., et al. (2001). Heart failure with a normal ejection fraction: is measurement of diastolic function necessary to make the diagnosis of diastolic heart failure? Circulation 104, 779–782. doi: 10.1161/hc3201.094226

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar

Zouein, F. A., de Castro Brás, L. E., da Costa, D. V., Lindsey, M. L., Kurdi, M., and Booz, G. W. (2013). Heart failure with preserved ejection fraction: emerging drug strategies. J. Cardiovasc. Pharmacol. 62, 13–21. doi: 10.1097/FJC.0b013e31829a4e61

PubMed Abstract | Full Text | CrossRef Full Text | Google Scholar