- 1Simula Research Laboratory, Department of Computational Physiology, Fornebu, Norway
- 2Departments of Bioengineering and Medicine, University of California San Diego, La Jolla, CA, United States
A wide range of arrhythmogenic phenotypes have been associated with heterogeneous mechanical dyskinesis. Pro-arrhythmic effects are often associated with dysregulated intra-cellular calcium handling, especially via the development of intra- and inter-cellular calcium waves. Experimental evidence suggests that mechanical strain can contribute to the generation and maintenance of these calcium waves via a variety of mechano-electric coupling mechanisms. Most model studies of mechano-electric coupling mechanisms have been focused on mechano-sensitive ion channels, even though experimental studies have shown that intra- and inter-cellular calcium waves triggered by mechanical perturbations are likely to be more prevalent pro-arrhythmic mechanisms in the diseased heart. A one-dimensional strongly coupled computational model of electromechanics in rabbit ventricular cardiomyocytes showed that specific myocyte stretch sequences can modulate the susceptibility threshold for delayed after-depolarizations. In simulations of mechanically-triggered calcium waves in cardiomyocytes coupled to fibroblasts, susceptibility to calcium wave propagation was reduced as the current through the gap junction caused current drain from the myocytes. In 1D multi-cellular arrays coupled via gap junctions, mechanically-induced waves may contribute to synchronizing arrhythmogenic calcium waves and after-depolarizations.
1. Introduction
Altered myocardial mechanics in disease may contribute to the triggering or maintenance of life-threatening reentrant ventricular arrhythmias via a variety of mechano-electric coupling mechanisms. In particular, altered myocyte mechanics may cause pro-arrhythmic perturbations to cardiomyocyte electrophysiology (Franz et al., 1992; Kohl et al., 2011). These perturbations may generate currents that trigger arrhythmias or create a heterogenous myocardial electrical substrate that supports reentry (Taggart and Lab, 2008). While most model studies of mechano-electric coupling mechanisms have focused on mechano-sensitive ion channels (Peyronnet et al., 2016), experimental studies (ter Keurs et al., 2006b), and recent mathematical simulations (Timmermann et al., 2019) suggest that intra- and inter-cellular calcium waves triggered by mechanical perturbations are likely to be more prevalent pro-arrhythmic mechanisms in the diseased heart.
To understand the mechanisms of mechanically-triggered arrhythmogenic currents, we must consider both how electrical excitation activates calcium-dependent muscle contraction [excitation-contraction coupling (ECC)] and the reverse processes by which stretch alters membrane currents [mechano-electric coupling (MEC)] (Kohl et al., 1999). Under normal conditions, both pathways are important for the electrical regulation of cardiac performance underlying mechanical changes (Kohl et al., 1999). At the organ scale, electrical dyssynchrony alters mechanical contraction and may further exacerbate pump dysfunction in disease (Pfeiffer et al., 2014). At the (sub)cellular level, dysregulated calcium handling is associated with proarrhythmia, notably via the development of calcium waves (Kass and Tsien, 1982; Matsuda et al., 1982; Capogrossi and Lakatta, 1985; Capogrossi et al., 1986b; Backx et al., 1989). Thus, mechanically-triggered, calcium-mediated mechanisms are likely to contribute to arrhythmogenic phenotypes. Despite intensive investigations at different levels, from single ion channels (Kohl and Sachs, 2001) to tissue (ter Keurs et al., 1998, 2006a,b; Kohl et al., 2001), the contributions of MEC to arrhythmia generation is incompletely understood.
In addition to MEC mechanisms, fibrosis and electrical remodeling are crucial modulators of long term outcomes accompanying pathological mechanical changes in the heart. Up to 70% of the total cell number in healthy rodent myocardium is comprised of fibroblasts (MacCannell et al., 2007; Maleckar et al., 2009), which are thought to electrotonically interact with myocytes via gap-junctions (MacCannell et al., 2007; Maleckar et al., 2009). These functional, dynamic fibroblasts-to-myocyte interactions can modulate electrical excitability and action potential (AP) waveform (MacCannell et al., 2007; Maleckar et al., 2009). Thus, fibroblast-mediated electrical remodeling together with mechanical heterogeneities might exacerbate the susceptibility to pro-arrhythmic calcium wave propagation by providing arrhythmogenic substrate (MacCannell et al., 2007; Maleckar et al., 2009).
The aim of this study is to mechanistically explore the critical role of altered myocyte mechanics, mechano-electric coupling, calcium homeostasis, and myocyte coupling in the generation and maintenance of ventricular arrhythmias, to increase the understanding of the complex interaction and balance of calcium homeostasis and contractile function. To address this goal, we simulated a one-dimensional (1D) cardiac myocyte, by strongly coupling a zero-dimensional AP model to sarcomere mechanics. The coupled model was modified to represent the cardiac myocyte as a 1D string of bi-directional end-to-end interactions between sarcomeres. We introduced mechanical heterogeneities by allowing stochastic variations in sarcomere resting length and by subjecting the model to a wide range of stretch protocols. Coupling mechanically heterogeneous myocytes by gap junctions increased the order of dynamic instability to trigger potentially arrhythmogenic delayed after-depolarizations (DADs), while coupling fibroblasts to the myocyte had a stabilizing affect.
2. Materials and Methods
2.1. Electro-Mechanical Myocyte Models
2.1.1. 0-Dimensional Myocyte Model
We developed a strongly coupled zero-dimensional (0D) computational model of rabbit cardiomyocyte electromechanics by coupling the well-established rabbit myocyte ionic model of Shannon et al. (2004) (distributed by Bers, 2002) with the myofilament mechanics model of Rice et al. (2008b) (downloaded from Rice et al., 2008a). The calcium transient generated by the ionic model, which represents normal calcium-induced calcium release (CICR) from the sarcoplasmic reticulum (SR) at the dyadic cleft, was input to the rabbit cardiac muscle contraction model. The formulation for cytosolic calcium binding to troponin C (TnC) in the mechanics model was used as input for the mechano-electric coupling of the mechanics model to the electrophysiology (EP) model. The contractile kinetics model described thin filament activation by intra-cellular calcium binding to TnC and the resulting crossbridge (XB) cycling.
In addition to strong electromechanical coupling, we incorporated strain-dependent changes in calcium affinity to the myofilaments to empirically approximate the triggers of mechanically-induced calcium waves. Using a parameter-fitting algorithm (Bueno-Orovio et al., 2008), the on- and off-rates for calcium binding to TnC and the constants for the strain-dependence of the XB cycling were reformulated to incorporate an exponential dependence on strain-rate. (A detailed description of the 0D model can be found in Timmermann et al., 2019).
2.1.2. 1-Dimensional Myocyte Model
We extended the 0D electromechanical model of section 2.1.1 by arraying the 0D model into 50 uniform mechanically and resistively coupled units. Each unit represents a single sarcomere of the mechanical model and a single calcium-release unit (CRU) of the EP model. This linear array model was previously constructed to investigate the ability of sub-cellular mechanical perturbations to trigger and modulate calcium waves via strain-dependent thin filament calcium affinity (Timmermann et al., 2019). As reported by Timmermann et al. (2019), time-dependent calcium diffusion was implemented for the sub-sarcolemma (SS), SR, and cytosol. The diffusion time constants for intra-CRU diffusion in the longitudinal direction were adapted to reproduce a calcium wave velocity of ~100 at Cao = 4 mM while increasing the magnitude of the calcium transient from the originally published value by Shannon et al. (2004). In particular, the calcium diffusion between the compartments of adjacent CRUs was incorporated when calculating the calcium concentrations as follows:
and for the boundaries:
with τSS = 0.4533 ms, τSR = 150 ms, τcyto = 1.2 ms, and , k ∈ {SS, SR, cyto}, the time-dependent changes calcium concentrations as in Shannon et al. (2004). (A detailed description of the one dimensional (1D) model can be found in Timmermann et al., 2019).
For this study, the model was extended by calculating the membrane potential in each of the 50 CRUs to better represent the electrical conductance of a myocyte. The first CRU in the myocyte was stimulated with a 1,000 mV pulse for 3 ms, triggering a electrical stimulation of the entire cell. Using a parameter-fitting algorithm (Bueno-Orovio et al., 2008), the time constant for voltage transfer between adjacent CRUs and the electrical pulses was chosen to allow electrical propagation while retaining a similar magnitude of the AP upstroke as reported in the original ionic model (Shannon et al., 2004).
with Cm membrane capacitance, ms conductance propagation coefficient, Istim the stimulus current, and Ii the total membrane current in CRU i ∈ {1, ..., 50}.
Furthermore, transmembrane currents with variable conductances were adapted to represent the mean experimentally reported values at the three cycle lengths 400, 600, and 1,000 ms as reported in Gemmell et al. (2016). In particular, the conductance of the transient outward potassium current (Ito) was increased by 30%, the conductance of slow delayed rectifier current (IKs) was decreased by 15%, the inward rectifier current (IK1) was decreased by 30%, and the sodium-potassium pump current was increased by 15%.
In addition to interactions between adjacent ionic CRUs in this spatially explicit model, a single sarcomere was coupled to each EP model. We used the strong coupling approach of Timmermann et al. (2019) to represent mechanical interaction between neighboring sarcomeres. Therefore, fractional changes in length between adjacent sarcomeres were included:
with the isotonic length change of sarcomere i of myocyte j.
The sarcomeres at the boundaries of each myocyte, interacted only with one other sarcomere of the same myocyte. To incorporate mechanical interactions between adjacent myocytes, adjacent sarcomeres of neighboring myocytes affected each other by their scaled fractional change in length
with the isotonic length change of the first (1) or the last (50) sarcomere, respectively, of the jth myocyte, and n = 3 or 6 for the number of coupled myocytes.
2.1.3. 1-Dimensional Multi-Cellular Model Coupled via Gap Junctions
To adapt the 1D electromechanical myocyte model for simulations of multi-cellular interactions, 3 myocytes were coupled via gap junctions. Each of the three myocytes in the coupled simulation was modeled as 1 of the 10 different single myocyte models described in section 2.1.2. Ten separate multi-cellular simulations with varying combinations of 3 of these 10 different single myocyte conditions were analyzed. The modeled combinations are given in Supplementary Tables 6, 7. To reproduce a conductance velocity of approximately 0.5 , the value of gap junction conductance was ggap = 400 nS. The corresponding equations for the linear gap junction current from the jth to the (j + 1)th cell were as follows:
with , , the change in membrane potential of the last CRU of myocyte j and the first CRU of myocyte (j + 1) in absence of the gap-junctional coupling, n, the number of myocytes, and , , the change in the potential of the membrane at last sarcomere of myocyte j and at the first sarcomere of myocyte (j + 1) in the presence of gap-junctional coupling.
For these multi-cellular simulations, the first CRU of the first myocyte (left boundary of the string of myocytes) was electrically stimulated with a 1,000 mV pulse for 3 ms. The adjacent myocytes were stimulated by the coupling via gap junctions.
2.2. Incorporation of Inter-cellular Mechanical Heterogeneities
We introduce mechanical heterogeneities by allowing sarcomeres within a single myocyte to have different resting lengths. To induce sarcomere heterogeneity, we deviated from the average resting length of 1.89 μm (Rice et al., 2008b), by 10% which are in the physiological range measured in myofibrils by Rassier et al. (2003). Therefore, for 10 random variations, the resting lengths of the 50 sarcomeres were chosen from the interval [1.701, 2.079 μm] (as shown in Supplementary Table 1) using a uniform random number generator (the rand function of MATLAB).
2.3. Incorporation of Fibroblasts
To simulate myocyte-fibroblast interactions, the active version 1 model of Maleckar et al. (2009) was incorporated to the electromechanics model. The current conductance of the fibroblasts is described by five time- and voltage-dependent equations in addition to the gap junction which allows for Na+ and K+-ions to move between the myocyte and the fibroblast (for a detailed description see Maleckar et al., 2009). Originating from the fibroblast model of MacCannell et al. (2007), the active model of Maleckar et al. (2009) includes the fibroblast membrane capacitance and four fibroblast ionic currents: time- and voltage-dependent fibroblast K+ current (IKV), inward-rectifying current (IK1) which is expressed in cardiac myofibroblasts, Na+-K+ pump current (INaK) which has expressed K+ channels resulting in K+ fluxes, and background Na+ current (Ib,Na) to balance the Na+ efflux from the Na+-K+ pump activity (for a detailed description see MacCannell et al., 2007).
Each myocyte was coupled to 1, 2, or 3 fibroblasts. Each fibroblast was randomly assigned to one CRU at 1 of the last 5 CRUs at one of the ends of the myocyte due to its relatively small size compared with the myocyte since fibroblast locate close to the gap junctions between myocytes (Goldsmith et al., 2004; Ongstad and Kohl, 2016). For each coupled fibroblast the value of the linear, fixed resistance (Gf,gap) varied randomly in the interval between the low and high ends of the range of Gf,gap, 0.5 and 8.0 nS as reported in Maleckar et al. (2009). Each of the 10 myocytes described in section 2.2, we assigned a random variation of the fibroblast locations (as shown in Supplementary Table 2).
2.4. Pacing and Stretch Protocols
In the electromechanics model we mimic various heterogeneities by subjecting the model to a wide range of stretch protocols. Myocyte stretch was modeled as an increase in sarcomere length. All stretch protocols of single cell and multi-cellular simulations were run for 200 beats to reach limit cycle. The limit cycle was ascertained by comparing the upstroke, plateau, and repolarization phase of the last two APs by calculating the mean squared error over time between both curves until mean squared error <0.01.
Under calcium overload conditions, calcium sparks may be more likely to trigger propagated waves (Backx et al., 1989; Cheng et al., 1996). Therefore, to study calcium sparks and propagating calcium waves, all simulations were performed by increasing extracellular calcium concentration in 0.1 mM steps from Cao = 2.0 mM to Cao = 4.5 mM, which results in increased calcium transient in all our simulations from the reported values of Shannon et al. (2004). Spontaneous calcium waves with wave velocities of ~ 100 have been reported at low extracellular calcium concentrations as 1.8–2.0 mM (Takamatsu and Wier, 1990; Wier and Blatter, 1991) while others reported wave velocities of ~72 at Cao = 2 mM, ~80 at Cao = 5 mM, and ~90 at Cao = 15 mM (Capogrossi et al., 1986a; Cheng et al., 1996).
We refer to normal healthy contraction as protocol H, a sinusoidal function decreasing to 20% compression, and isometric contraction as protocol I. Pathophysiological stretches were represented by a sinusoidal function decreasing to 10% compression in diastole and increasing to 10% stretch in systole (protocol S) and an exponentially increasing function, increasing from 0 to 15% at the start of contraction, and then kept constant until it decreased exponentially at the time of relaxation (protocol P) (as shown in Figure 1). These protocols were used both for simulations of single and coupled myocytes both with and without coupled fibroblasts. Additionally, all four stretch protocols were applied to the 10 different single myocyte models generated in section 2.2 as well as the 10 different single myocyte models with 1, 2, or 3 coupled fibroblasts as described in section 2.3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g001.jpg)
Figure 1. Myocyte stretch protocols. Stretch pattern H (light blue) approximates normal isotonic myocyte shortening with a sinusoidal function. Protocol I (red) is sarcomere isometric contraction. Protocol S (dark blue) simulates late systolic stretch by shifting the sine wave upwards. Protocol P (orange) simulates systolic lengthening such as occurs in acutely ischemic regions.
All single cell simulations were run in MATLAB R2018a using ode15s to solve the ordinary differential equation system, while multi-myocyte simulations were run in C++ using odeint.
2.5. Data Analysis
The statistical tests are presented as mean ± standard deviation calculated in MATLAB R2018a using the functions mean and std. The data was analyzed in MATLAB R2018a with the functions anova1 and multicompare, a one-way analysis of variance with pairwise comparison of the group means preconditioned by the Bonferroni method. A data point was assumed to be an outlier if it did not lie within three standard deviations of the mean (three-sigma rule). Significant differences between were illustrated by asterisks with three asterisks indicating a p < 0.01, two a p < 0.05, and one a p < 0.1. Statistical tests were assessed for all simulations, even though some significant differences are far below the accuracy of experimental measurement techniques. To measure the inter-cellular dyssynchrony of intra-cellular calcium dynamics, we evaluated the standard deviation between cells of the mean time-to-peak of the intra-cellular calcium transient within each cell.
2.6. Data Availability
The datasets for this study can be found in the repository: https://github.com/cmrglab/1D_MEC.git.
3. Results
3.1. The Occurrence of Calcium Waves and Wave Velocity Is Dependent on the Stretch Pattern and Cao
For all single myocyte simulations, we examined susceptibility to the generation of calcium waves after both pacing and stretching were stopped.
Even though some studies (Takamatsu and Wier, 1990; Wier and Blatter, 1991) have reported propagating calcium waves at low calcium concentrations as Cao = 1.8–2.0 mM, we did not observe spontaneous calcium waves in our model at Cao = 2.0 mM. Therefore, we increased the extracellular calcium concentration in 0.1 mM steps until in 3 of 10 control experiments (no stretch, stretch pattern I) spontaneously propagating calcium waves were observed at Cao = 4.1 mM. These findings are largely consistent with many experimental reports of calcium waves, in which the extracellular calcium concentrations varied between ~2.0 mM (Takamatsu and Wier, 1990; Wier and Blatter, 1991) and 15 mM (Capogrossi et al., 1986a; Cheng et al., 1996).
Figure 2B shows that no calcium waves were observed for the control experiments at Cao = 2.0 mM (indicated by red daggers). For isotonic stretch simulations, the first calcium waves occurred at Cao = 4.3 mM, while all for all other stretch patterns lower calcium concentrations were sufficient to trigger spontaneously propagating calcium waves (see Figure 2B). At Cao = 4.2 mM, stretch pattern P generated the fastest calcium waves which always triggered an additional beat (indicated by the red star in Figure 2B) and propagated with a mean velocity of 4596.1 ± 408.22 . This additional beat was followed by a slow calcium wave which propagated with a 175.58 ± 37.98 (shown in Supplementary Table 3). Even though the extracellular calcium concentration at which a calcium wave occurred for the different stretch patterns varied between Cao = 3.7 and 4.3 mM, the wave velocity was ~100 for all stretch patterns. While the velocity of the first calcium wave was similar between stretch patterns, the maximal calcium wave amplitude (shown in Figure 2A) varied depending on the stretch pattern due to the extracellular calcium concentration. The lowest mean of the maximal calcium wave amplitude was observed for stretch pattern H with a magnitude of 0.15 ± 0.0 μM and the highest maximal calcium wave amplitude for stretch pattern P with a magnitude of 4.46 ± 0.13 μM.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g002.jpg)
Figure 2. Maximal calcium wave amplitude magnitude and calcium wave velocity in cardiomyocyte simulations at different extracellular calcium concentrations for control (no stretch protocol I), and stretch patterns H, S, and P. (A) The magnitude of the measured maximal calcium wave amplitude increased significantly for all stretch patterns from Cao = 2.0 mM (red data points) to the extracellular calcium concentration at which a calcium wave occurred (blue data points), and 4.2 mM (green data points). Red daggers indicate that no calcium waves occurred in the data set, while red asterisks indicate that the calcium waves triggered additional beats. The related numbers specifies the amount of occurrences within the data set. The data points in light colors specify the first spontaneous event and dark colors are used for the following spontaneous event. (B) The calcium wave velocity increased for all simulations from Cao = 2.0 mM (no waves) to 4.2 mM and from control, to stretch patterns H, to S, and to P. **, *** indicate the statistical significance with p < 0.05 and p < 0.1, respectively.
At Cao = 4.2 mM, the greatest variation in calcium wave velocity was observed for stretch pattern P with a velocity of 4005.5 ± 408.8 . For all simulations at Cao = 4.2 mM, the resulting calcium waves induced suprathreshold DADs. For stretch pattern P in all 10 myocytes, a fast calcium wave triggered an extra beat followed by a slow propagating calcium wave, while for stretch pattern S in 6 of 10 myocytes, an additional beat and slow wave was observed in only two simulations. An example for the different calcium-mediated effects on AP are shown in Figure 3A and the associated calcium waves responsible for the changes in AP are illustrated in Figure 3B. At Cao = 4.2 mM a fast propagating wave triggered an additional beat, which was followed by a slower propagating wave. As described in Timmermann et al. (2019), calcium waves arose from spontaneous calcium release from the SR as a result of SR calcium overload. SR calcium overload resulted from stretch-dependent myofilament calcium dissociation that was sufficient to alter cytosolic calcium dynamics.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g003.jpg)
Figure 3. (A) Membrane potential for stretch pattern S of myocyte 3 at Cao = 4.0 mM and Cao = 4.2 mM (from left to right). After one stimulated beat, at Cao = 4.0 mM a calcium-mediated suprathreshold DAD occurred and at Cao = 4.2 mM a calcium-mediated additional beat occurred followed by a calcium wave. (B) Associating calcium wave for the supratheshold DAD, and the additional beat followed by a calcium wave (from left to right) at Cao = 4.0 mM and Cao = 4.2 mM, respectively.
Figure 2B shows that an increase in calcium wave velocity is associated with an increase in maximal calcium wave amplitude. For control (stretch protocol I), the calcium wave amplitude significantly (p < 0.01) increased for experiments at Cao = 4.1 mM to Cao = 4.2 mM from 1246.0 ± 819.0 to 1896.3 ± 52.0 (see Figure 2A). Similarly, the maximal calcium wave amplitude for simulations for stretch patterns S and P increased significantly (p < 0.01) for experiments at which calcium waves first occurred to Cao = 4.2 mM by 637.25, 217.33%, respectively. (Mean and standard deviation of the maximal calcium wave amplitude magnitude and calcium wave velocity can be found in Supplementary Table 3).
3.2. Fibroblasts Function as Current Drains in Simulations of Spontaneous Calcium Waves
Owing to the large number of fibroblasts in the myocardium, we investigated their role in modulating electrophysiological instabilities. For all stretch patterns of myocytes with coupled fibroblasts, the current through the gap junction provided a drain of electric charge from the cardiomyocyte. Therefore, we were able to replicate a significant (p < 0.01) shortening of APD90 (see Figure 4 and in Supplementary Table 4) as reported by MacCannell et al. (2007).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g004.jpg)
Figure 4. APD90 for simulations of myocytes with no (red data points), 1 (blue data points), and 3 coupled fibroblast (green data points) compared with control (no coupled fibroblast) for no contraction and stretch patterns H, S, and P. APD90 was significantly reduced for all simulations of myocytes with coupled myocytes for the different stretch patterns.
For simulations of calcium-mediated DADs the current through the gap junction between the fibroblast and the myocyte supplied an electric discharge from the cardiomyocyte as illustrated by the significantly (p < 0.01) decreasing maximal calcium wave amplitude in Figure 5A. Associated with decreased maximal calcium wave amplitude magnitude, the velocity of the induced calcium wave decreased significantly (p < 0.01) as more fibroblasts were coupled to the myocyte (see Figure 5B). For no stretch, in simulations with 3 coupled fibroblasts in 1 of 10 myocytes, no wave occurred, while for simulations of stretch patterns P and S, calcium waves were present. However, for simulations of stretch pattern P with three coupled fibroblasts no slow calcium waves occurred after the triggered beat. calcium waves in myocytes coupled to three fibroblasts did not elicit extra beats in simulations of stretch pattern S.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g005.jpg)
Figure 5. Maximal calcium wave amplitude and calcium wave velocity for simulations for control (red data points), and myocytes with one (blue data points) and two (green data points) coupled fibroblasts. The data points in light colors specify the first spontaneous event and dark colors are used for the following spontaneous event. (A) For stretch patterns I (control), S, and P the maximal calcium wave amplitude significantly (p < 0.01) decreased the more fibroblasts were coupled to the myocyte. (B) The wave velocity significantly (p < 0.01) decreased for stretch patterns I (control), S, and P for no fibroblast coupled to three coupled fibroblasts. *, **, *** indicate the statistical significance with p < 0.01, p < 0.05, and p < 0.1, respectively.
The calcium waves in our computational model are triggered through spontaneous calcium release from the SR, which in turn triggers an increase in membrane potential and potentially, a DAD. Therefore, the coupled fibroblasts not only drain electrical charge and reduce the APD90 during normal pacing, but also slow calcium waves as charge is drained from the associated calcium-mediated suprathreshold DAD.
3.3. Calcium Waves Are More Synchronized in Multi-Cellular Simulations
In simulations of 3 and 6 cardiomyocytes coupled via gap junctions, calcium wave velocity and calcium wave amplitude of the multi-cellular simulation are given as the mean of the individually measured observed calcium wave velocities and calcium wave amplitudes of each single, coupled cell. Calcium wave velocity slightly decreased for all stretch patterns and for fibroblast coupling compared to single myocytes (see Figures 5, 6, respectively). The results shown here were computed for three coupled myocytes but no differences were observed when we repeated these calculations with six coupled cells (shown in Supplementary Figure 1). Coupled myocytes behaved similarly to single myocytes. Maximal calcium wave amplitude varied between simulations of no stretch, and stretch patterns P and H with the highest maximal calcium wave amplitude occurring for stretch pattern P, which coincides with the fastest observed wave velocity (as shown in Figure 6). With stretch pattern S, fewer additional beats were triggered in the coupled myocyte simulations, showing how multi-cellular tissue can be less susceptible to extra beats than isolated cells due to the slowing and synchronization of calcium waves by electrotonic coupling effects.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g006.jpg)
Figure 6. Maximal calcium wave amplitude and calcium wave velocity for simulations of three coupled cardiomyocytes. (A) For all stretch patterns but H (from left to right: control, H, S, and P) the maximal calcium wave amplitude significantly (p < 0.01) decreased the more fibroblasts were coupled to the multi-cellular string (control red data points, 1 coupled fibroblast blue data points, and 3 coupled fibroblasts green data points; the data points in light colors specify the first spontaneous event and dark colors are used for the following spontaneous event). (B) The wave velocity significantly (p < 0.01) decreased for all stretch patterns from no coupled fibroblast to three coupled fibroblasts. *, **, *** indicate the statistical significance with p < 0.01, p < 0.05, and p < 0.1, respectively.
Electrotonic coupling of cardiomyocytes reduced inter-cellular dyssynchrony of intra-cellular calcium dynamics for all 10 coupled cell replicates. As shown in Figure 7B, we measured more synchronized DAD rise times in multi-cellular compared with isolated cell simulations (see Figure 7A). Current conduction of triggered suprathreshold DADs via gap junctions triggered calcium waves with synchronized time-to-peak calcium and wave propagation velocities in all three cells. For coupled cells, the inter-cellular dyssynchrony of intra-cellular calcium dynamics was 0 within the numerical range. By contrast, the dyssynchrony of calcium wave velocity and the DAD rise times of isolated cells varied up to 378.56 ± 170.71 and 36.46 ± 23.45 ms, respectively. However, DAD magnitude decreased in coupled cell simulations compared with isolated myocytes due to the greater synchrony.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/519951/fphys-11-519951-HTML/image_m/fphys-11-519951-g007.jpg)
Figure 7. (A) Examples of three isolated myocyte APs of the last stimulated beat and the calcium mediated DAD (top), a detailed view of the DAD (bottom left), and the calcium waves in three single cardiomyocytes coupled to three fibroblasts each (bottom right). (B) Example of the APs of three coupled myocytes showing the last stimulated beat and the DAD (top), a detailed view of the DAD (bottom left), and the calcium waves in three coupled cardiomyocytes coupled to three fibroblasts each (bottom right).
4. Discussion
This computational study examined the importance of MEC, and particular, myofilament-triggered calcium release in the context of chronic stretch. We have assessed the arrhythmogenic potential of four different stretch patterns at cellular and multi-cellular scales with and without coupled fibroblasts. Our results support the hypothesis that MEC might be pro-arrhythmic but only under specific conditions of stretch and calcium overload.
Our model suggests that healthy, isotonically contracting myocytes may be less susceptible to generating propagating calcium waves than isometric contracting myocytes, while stretched myocytes are the most susceptible to calcium waves. For myocytes that do not contract but only stretch over the duration of the twitch (stretch pattern P), the extracellular calcium concentration necessary to trigger propagating waves was 16% lower than in isotonically contracting myocytes. Nevertheless, an increase in extracellular calcium concentration (Cao ≥ 3.7 mM) was necessary for all stretch patterns to trigger calcium waves. When fibroblasts were coupled to myocytes, electrical charge was drained from the myocytes reducing APD90 during normal pacing and slowing calcium waves.
In simulations of myocytes coupled via gap junctions, the calcium wave was synchronized between cells but at a lower amplitude compared with those in the single cell simulations. Cell-cell coupling reduced DAD magnitudes, but synchronized the rise times of the triggered suprathreshold DADs due to current conduction via gap junctions.
4.1. Sub-cellular Structure and Calcium Handling
Spontaneous calcium release events have not only been associated with a calcium overloaded SR (Takamatsu and Wier, 1990; Williams et al., 1992; Cheng et al., 1993, 1996; Satoh et al., 1997; Cheng and Lederer, 2008) but also as the site of calcium wave initiation (Cheng et al., 1996). An increased frequency of calcium sparks is thought to contribute to the initiation of DADs and cellular triggered activity (Voigt et al., 2013). The wavefront of slowly propagating calcium waves are the result of discrete calcium releases from the SR, which recruit other sparks during propagation (Cheng et al., 1996). The incidence of calcium sparks, and thus propagating waves, has also been associated with mechanical heterogeneities (Miura et al., 1993, 1999; Wakayama et al., 2001; ter Keurs et al., 2006a, 2008). Experiments in damaged cardiac muscle discovered that mechanical strain can modify calcium wave velocity resulting in an increase in the propagation velocity (ter Keurs et al., 1998). Consistent with experimental and computational observations, calcium wave velocity affects the amplitude of DADs, which can act as the initiator of arrhythmia at the tissue level (Daniels and ter Keurs, 1990; Daniels et al., 1991; Miura et al., 1993; ter Keurs et al., 2008). Indeed, experiments have shown that premature beats can result from propagating calcium waves that originate close to the border zones of mechanically-heterogeneous tissues (Boyden et al., 2003; Boyden and ter Keurs, 2005; ter Keurs and Boyden, 2007). Therefore, MEC may be an important contributor to ventricular arrhythmia in cardiac pathologies associated with mechanical perturbations.
The patterns of systolic stretch used in our simulations are similar to those described in regions of myocardial mechanical heterogeneity. Experimental studies of acute ischemia (Baumeister et al., 2018) revealed that ischemic tissue shows a similar stretches as stretch pattern P used in for this computational study. Late diastolic stretches, such as stretch pattern S, have been observed in various pathological conditions such as ischemia or heart failure (Neves et al., 2016) and isometric conditions are often used in experimental set-ups to immobilize the cells for optical mapping. In pathological conditions, such as ischemia, stretch is assumed to enhance calcium spark probability (Cameron et al., 2020), which may trigger calcium waves and DADs originating from these sites (Boyden et al., 2003; Boyden and ter Keurs, 2005; ter Keurs and Boyden, 2007). In addition to calcium waves, spontaneous contractile waves have been related to pro-arrhythmic effects in arrangements of weaker damaged regions in-series with normal muscle segments. Even though the underlying arrhythmogenic mechanism is thought to be caused by calcium release, the initiation of the calcium release has been linked to stretch and release of the myofilaments. The increase in intra-cellular calcium may be conducted to adjacent myocytes triggered by tissue-level heterogeneities in mechanical loading the resulting spontaneous contractile waves (ter Keurs et al., 1998).
4.2. Comparison With Previous Models
Cardiomyocyte function is tightly regulated by calcium signaling mechanisms that are well studied. While fibroblasts are also regulated by calcium fluxes (Feng et al., 2019), the mechanisms of calcium homeostasis in these non-excitable cells, and the role of voltage-gated calcium channels are less well studied and modeled (Feng et al., 2019). Mathematical models have not yet incorporated calcium fluxes between coupled myocytes and fibroblasts (MacCannell et al., 2007; Maleckar et al., 2009), and detailed models of calcium handling in fibroblasts are independent of membrane potential (Kotwani et al., 2014; Kotwani, 2015). Thus, there is a need for better models of fibroblast calcium fluxes and their regulation by transmembrane potential.
Spontaneous calcium release has been studied extensively in experiments (Capogrossi and Lakatta, 1985; Sasse et al., 2007; Rapila et al., 2008) as well as computational studies (Means et al., 2006; Dupont et al., 2007; Solovey et al., 2008). To study myocyte calcium handling, several computational models have been developed ranging from 0D to three-dimensional (3D) implementations and from highly detailed molecular models to common pool models (Means et al., 2006; Dupont et al., 2007; Solovey et al., 2008). Since calcium handling abnormalities are often related to cardiac diseases as heart failure, models have also focused on ryanodine receptor (RyR) and L-type channel dysfunction. Both, dysfunction or location heterogeneities of RyR clusters and L-Type channels have been associated with an impact on spontaneous calcium release events, and thus, impairment of calcium waves or pro-arrhythmic calcium waves (Zahradńıková and Zahradńık, 2012; Walker et al., 2014). However, to our knowledge none of these models include MEC to study the effects of mechanical heterogeneities in cells or tissues.
In larger tissue scale studies gap junctions play an important role. While our study only considered a linear gap junction model, various models of more detailed static and dynamic gap junction models have previously been developed. Vogel and Weingart (1998) developed a static model of gap junctions in which the voltage- and time-dependent conductance changes based on the voltage in myocyte pairs. A more recent dynamic model by Henriquez et al. (2001) describes each gap junction as two hemichannels in-series, which exist in either a high or low conductance state. Hence, the model assumes that the conductive gap junction channels have four conformational states regulated by nonlinear functions of the junctional voltage (Oka et al., 2006). Experimental data (Oka et al., 2006) also suggest that gap junction gating is regulated by intra-cellular calcium and protons. While a more detailed gap junction model would produce more physiological results than our linear model, they are also more computationally expensive.
4.3. Limitations
Our model was unable to elicit calcium waves at extracellular calcium concentrations below 3.7 mM, similar to some experimental studies that have used unphysiologically high extracellular calcium concentrations (Capogrossi et al., 1986a; Cheng et al., 1996). However, this is a limitation of the model because other experimental studies have observed calcium waves at concentrations as low as 1.8–2.0 mM (Takamatsu and Wier, 1990; Wier and Blatter, 1991). There are several aspects of the calcium handling mechanisms in our model that might account for the limitation that high extracellular calcium concentrations were needed to elicit propagating calcium waves. In particular, the present model does not incorporate calcium fluxes between myocytes and fibroblasts, which have been reported in experiments (Feng et al., 2019) and have been modeled mathematically by one group (Kotwani et al., 2014; Kotwani, 2015).
Our model suggests that specific chronic stretch patterns are likely to have a pro-arrhythmic effect under calcium overload conditions. However, the myocyte model only considered a 1D representation in longitudinal direction. Therefore, the model cannot give insights into cellular anisotropy and radial wave propagation (Galice et al., 2018). Additionally, we only consider sarcomeric heterogeneity but neglect possibly important details about sub-cellular distribution of RyR clusters and L-type channels, which might have a strong impact on the calcium handling (Sutanto et al., 2018). Furthermore, our conclusions are based on a deterministic model while stochastic approaches are necessary to reproduce heterogeneous systems as RyR clusters which can open spontaneously. Therefore, our model could not generate propagating waves that do not activate the whole myocyte.
The myofilament contraction model (Rice et al., 2008b) was chosen in part because of its computational efficiency compared with spatial explicit sarcomere models. It was built largely with steady-state data and a simplified representation of calcium binding to the myofilaments, while more recent models allow for a more realistic representation of this process. The parameterization of calcium diffusion in the EP model (Shannon et al., 2004) was based on ECC data and can only partially reproduce MEC.
Finally, in this study we only considered one MEC mechanism, even though mechano-sensitive ion channels and XROS are also known to contribute to mechanics-induced arrhythmia (Kohl et al., 2011). Nevertheless, mechanically-induced calcium waves have also been reported when mechano-sensitive ion channels were blocked (Wakayama et al., 2001). XROS on the other hand might be important in some context, especially during diastolic stretch (Prosser et al., 2011, 2013). But we did not explicitly simulate XROS since we used whole muscle data for the parameterization of the myofilament-triggered calcium release events.
5. Conclusions
In this study, we investigated the ability of sub-cellular mechanical transients and heterogeneity and the effects of coupled fibroblasts to modify calcium wave velocity and the resulting calcium mediated DADs. The goal was to investigate whether MEC, transient stretches, and fibroblasts contribute to pro-arrhythmic effects in calcium overloaded myocytes and multi-cellular tissues.
The present study indicates that myofilament-triggered calcium release mechanisms may modulate the susceptibility threshold for DADs. Chronic stretch may increase cytosolic calcium concentration resulting in an increase of the magnitude of the calcium-mediated DAD. In contrast, myocyte coupled to fibroblasts are less susceptible to calcium-mediated DADs as the current through the gap junction between the myocyte and the fibroblasts may provide a drain of the electric charge. Finally, at the multi-cellular scale, mechanically-induced calcium waves may trigger synchronized calcium waves and after-depolarizations. Thus, MEC may increase the susceptibility to pro-arrhythmic intra- and inter-cellular calcium wave propagation.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
VT carried out all simulations, analyzed the data, and wrote the article. Together with VT, AM designed the research. AM commented and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This study has been supported by Simula Research Laboratory funded by the Norwegian Ministry of Education and Research. This work was also supported in part by National Institutes of Health awards (HL105242, HL137100, HL122199, HL126273) and the National Biomedical Computation Resource (P41 GM103426 to AM).
Conflict of Interest
AM is a co-founder of and has an equity interest in Insilicomed and Vektor Medical. He serves on the scientific advisory board of Insilicomed and as scientific advisor to both companies. Some of his research grants, including those acknowledged here, have been identified for conflict of interest management based on the overall scope of the project and its potential benefit to these companies.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.519951/full#supplementary-material
References
Backx, P. H., De Tombe, P. P., Van Deen, J., Mulder, B., and Ter Keurs, H. (1989). A model of propagating calcium-induced calcium release mediated by calcium diffusion. J. Gen. Physiol. 93, 963–977. doi: 10.1085/jgp.93.5.963
Baumeister, P. A., Lawen, T., Rafferty, S. A., Taeb, B., Uzelac, I., Fenton, F. H., et al. (2018). Mechanically-induced ventricular arrhythmias during acute regional ischemia. J. Mol. Cell. Cardiol. 124, 87–88. doi: 10.1016/j.yjmcc.2018.07.021
Bers, D. M. (2002). Available online at: https://somapp.ucdmc.ucdavis.edu/Pharmacology/bers/.
Boyden, P. A., Barbhaiya, C., Lee, T., and ter Keurs, H. E. (2003). Nonuniform Ca2+ transients in arrhythmogenic Purkinje cells that survive in the infarcted canine heart. Cardiovasc. Res. 57, 681–693. doi: 10.1016/S0008-6363(02)00725-3
Boyden, P. A., and ter Keurs, H. (2005). Would modulation of intracellular Ca2+ be antiarrhythmic? Pharmacol. Therap. 108, 149–179. doi: 10.1016/j.pharmthera.2005.03.011
Bueno-Orovio, A., Cherry, E. M., and Fenton, F. H. (2008). Minimal model for human ventricular action potentials in tissue. J. Theor. Biol. 253, 544–560. doi: 10.1016/j.jtbi.2008.03.029
Cameron, B. A., Kai, H., Kaihara, K., Iribe, G., and Quinn, T. A. (2020). Ischemia enhances the acute stretch-induced increase in calcium spark rate in ventricular myocytes. Front. Physiol. 11:289. doi: 10.3389/fphys.2020.00289
Capogrossi, M. C., Kort, A. A., Spurgeon, H. A., and Lakatta, E. G. (1986a). Single adult rabbit and rat cardiac myocytes retain the Ca2+-and species-dependent systolic and diastolic contractile properties of intact muscle. J. Gen. Physiol. 88, 589–613. doi: 10.1085/jgp.88.5.589
Capogrossi, M. C., and Lakatta, E. G. (1985). Frequency modulation and synchronization of spontaneous oscillations in cardiac cells. Am. J. Physiol. 248, H412–H418. doi: 10.1152/ajpheart.1985.248.3.H412
Capogrossi, M. C., Suarez-Isla, B. A., and Lakatta, E. G. (1986b). The interaction of electrically stimulated twitches and spontaneous contractile waves in single cardiac myocytes. J. Gen. Physiol. 88, 615–633. doi: 10.1085/jgp.88.5.615
Cheng, H., Lederer, M. R., Lederer, W., and Cannell, M. (1996). Calcium sparks and [Ca2+]i waves in cardiac myocytes. Am. J. Physiol. 270, C148–C159. doi: 10.1152/ajpcell.1996.270.1.C148
Cheng, H., and Lederer, W. (2008). Calcium sparks. Physiol. Rev. 88, 1491–1545. doi: 10.1152/physrev.00030.2007
Cheng, H., Lederer, W. J., and Cannell, M. B. (1993). Calcium sparks: elementary events underlying excitation-contraction coupling in heart muscle. Science 262, 740–744. doi: 10.1126/science.8235594
Daniels, M. C. G., Fedida, D., Lamont, C., and ter Keurs, H. E. D. J. (1991). Role of the sarcolemma in triggered propagated contractions in rat cardiac trabeculae. Circ. Res. 68, 1408–1421. doi: 10.1161/01.RES.68.5.1408
Daniels, M. C. G., and ter Keurs, H. E. D. J. (1990). Spontaneous contractions in rat cardiac trabeculae; trigger mechanism and propagation velocity. J. Gen. Physiol. 95, 1123–1137. doi: 10.1085/jgp.95.6.1123
Dupont, G., Combettes, L., and Leybaert, L. (2007). Calcium dynamics: spatio-temporal organization from the subcellular to the organ level. Int. Rev. Cytol. 261, 193–245. doi: 10.1016/S0074-7696(07)61005-5
Feng, J., Armillei, M. K., Yu, A. S., Liang, B. T., Runnels, L. W., and Yue, L. (2019). Ca2+ signaling in cardiac fibroblasts and fibrosis-associated heart diseases. J. Cardiovasc. Dev. Dis. 6:34. doi: 10.3390/jcdd6040034
Franz, M. R., Cima, R., Wang, D., Profitt, D., and Kurz, R. (1992). Electrophysiological effects of myocardial stretch and mechanical determinants of stretch-activated arrhythmias. Circulation 86, 968–978. doi: 10.1161/01.CIR.86.3.968
Galice, S., Xie, Y., Yang, Y., Sato, D., and Bers, D. M. (2018). Size matters: ryanodine receptor cluster size affects arrhythmogenic sarcoplasmic reticulum calcium release. J. Am. Heart Assoc. 7:e008724. doi: 10.1161/JAHA.118.008724
Gemmell, P., Burrage, K., Rodríguez, B., and Quinn, T. A. (2016). Rabbit-specific computational modelling of ventricular cell electrophysiology: using populations of models to explore variability in the response to ischemia. Prog. Biophys. Mol. Biol. 121, 169–184. doi: 10.1016/j.pbiomolbio.2016.06.003
Goldsmith, E. C., Hoffman, A., Morales, M. O., Potts, J. D., Price, R. L., McFadden, A., et al. (2004). Organization of fibroblasts in the heart. Dev. Dyn. 230, 787–794. doi: 10.1002/dvdy.20095
Henriquez, A. P., Vogel, R., Muller-Borer, B. J., Henriquez, C. S., Weingart, R., and Cascio, W. E. (2001). Influence of dynamic gap junction resistance on impulse propagation in ventricular myocardium: a computer simulation study. Biophys. J. 81, 2112–2121. doi: 10.1016/S0006-3495(01)75859-6
Kass, R. S., and Tsien, R. W. (1982). Fluctuations in membrane current driven by intracellular calcium in cardiac Purkinje fibers. Biophys. J. 38, 259–269. doi: 10.1016/S0006-3495(82)84557-8
Kohl, P., Hunter, P., and Noble, D. (1999). Stretch-induced changes in heart rate and rhythm: clinical observations, experiments and mathematical models. Prog. Biophys. Mol. Biol. 71, 91–138. doi: 10.1016/S0079-6107(98)00038-8
Kohl, P., Nesbitt, A. D., Cooper, P. J., and Lei, M. (2001). Sudden cardiac death by commotio cordis: role of mechano-electric feedback. Cardiovasc. Res. 50, 280–289. doi: 10.1016/S0008-6363(01)00194-8
Kohl, P., and Sachs, F. (2001). Mechanoelectric feedback in cardiac cells. Philos. Trans. R. Soc. Lond. Ser. A 359, 1173–1185. doi: 10.1098/rsta.2001.0824
Kohl, P., Sachs, F., and Franz, M. R. (2011). Cardiac Mechano-Electric Coupling and Arrhythmias. Oxford: Oxford University Press. doi: 10.1093/med/9780199570164.001.0001
Kotwani, M. (2015). Modeling and simulation of calcium dynamics in fibroblast cell involving excess buffer approximation (EBA), ER flux and SERCA pump. Proc. Comput. Sci. 49, 347–355. doi: 10.1016/j.procs.2015.04.263
Kotwani, M., Adlakha, N., and Mehta, M. (2014). Finite element model to study the effect of buffers, source amplitude and source geometry on spatio-temporal calcium distribution in fibroblast cell. J. Med. Imag. Health Inform. 4, 840–847. doi: 10.1166/jmihi.2014.1328
MacCannell, K. A., Bazzazi, H., Chilton, L., Shibukawa, Y., Clark, R. B., and Giles, W. R. (2007). A mathematical model of electrotonic interactions between ventricular myocytes and fibroblasts. Biophys. J. 92, 4121–4132. doi: 10.1529/biophysj.106.101410
Maleckar, M. M., Greenstein, J. L., Giles, W. R., and Trayanova, N. A. (2009). Electrotonic coupling between human atrial myocytes and fibroblasts alters myocyte excitability and repolarization. Biophys. J. 97, 2179–2190. doi: 10.1016/j.bpj.2009.07.054
Matsuda, H., Noma, A., Kurachi, Y., and Irisawa, H. (1982). Transient depolarization and spontaneous voltage fluctuations in isolated single cells from guinea pig ventricles. Calcium-mediated membrane potential fluctuations. Circul. Res. 51, 142–151. doi: 10.1161/01.RES.51.2.142
Means, S., Smith, A. J., Shepherd, J., Shadid, J., Fowler, J., Wojcikiewicz, R. J., et al. (2006). Reaction diffusion modeling of calcium dynamics with realistic ER geometry. Biophys. J. 91, 537–557. doi: 10.1529/biophysj.105.075036
Miura, M., Boyden, P. A., and Keurs, H. E. t. (1999). Ca2+ waves during triggered propagated contractions in intact trabeculae: determinants of the velocity of propagation. Circul. Res. 84, 1459–1468. doi: 10.1161/01.RES.84.12.1459
Miura, M., Ishide, N., Oda, H., Sakurai, M., Shinozaki, T., and Takishima, T. (1993). Spatial features of calcium transients during early and delayed after depolarizations. Am. J. Physiol. Heart Circul. Physiol. 265, H439–H444. doi: 10.1152/ajpheart.1993.265.2.H439
Neves, J. S., Leite-Moreira, A. M., Neiva-Sousa, M., Almeida-Coelho, J., Castro-Ferreira, R., and Leite-Moreira, A. F. (2016). Acute myocardial response to stretch: what we (don't) know. Front. Physiol. 6:408. doi: 10.3389/fphys.2015.00408
Oka, C., Matsuda, H., Sarai, N., and Noma, A. (2006). Modeling the calcium gate of cardiac gap junction channel. J. Physiol. Sci. 56, 79–85. doi: 10.2170/physiolsci.R2139
Ongstad, E., and Kohl, P. (2016). Fibroblast-myocyte coupling in the heart: potential relevance for therapeutic interventions. J. Mol. Cell. Cardiol. 91, 238–246. doi: 10.1016/j.yjmcc.2016.01.010
Peyronnet, R., Nerbonne, J. M., and Kohl, P. (2016). Cardiac mechano-gated ion channels and arrhythmias. Circul. Res. 118, 311–329. doi: 10.1161/CIRCRESAHA.115.305043
Pfeiffer, E. R., Tangney, J. R., Omens, J. H., and McCulloch, A. D. (2014). Biomechanics of cardiac electromechanical coupling and mechanoelectric feedback. J. Biomech. Eng. 136:021007. doi: 10.1115/1.4026221
Prosser, B. L., Khairallah, R. J., Ziman, A. P., Ward, C. W., and Lederer, W. (2013). X-ROS signaling in the heart and skeletal muscle: stretch-dependent local ROS regulates [Ca2+]i. J. Mol. Cell. Cardiol. 58, 172–181. doi: 10.1016/j.yjmcc.2012.11.011
Prosser, B. L., Ward, C. W., and Lederer, W. (2011). X-ROS signaling: rapid mechano-chemo transduction in heart. Science 333, 1440–1445. doi: 10.1126/science.1202768
Rapila, R., Korhonen, T., and Tavi, P. (2008). Excitation-contraction coupling of the mouse embryonic cardiomyocyte. J. Gen. Physiol. 132, 397–405. doi: 10.1085/jgp.200809960
Rassier, D. E., Herzog, W., and Pollack, G. H. (2003). Dynamics of individual sarcomeres during and after stretch in activated single myofibrils. Proc. R. Soc. Lond. Ser. B Biol. Sci. 270, 1735–1740. doi: 10.1098/rspb.2003.2418
Rice, J. J., Wang, F., Bers, D. M., and de Tombe, P. P. (2008a). Available online at: https://models.cellml.org/exposure/a31519a27f4c2fe6158e04fd40eeda98/rice_wang_bers_detombe_2008.cellml/view
Rice, J. J., Wang, F., Bers, D. M., and de Tombe, P. P. (2008b). 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
Sasse, P., Zhang, J., Cleemann, L., Morad, M., Hescheler, J., and Fleischmann, B. K. (2007). Intracellular Ca2+ oscillations, a potential pacemaking mechanism in early embryonic heart cells. J. Gen. Physiol. 130, 133–144. doi: 10.1085/jgp.200609575
Satoh, H., Blatter, L. A., and Bers, D. M. (1997). Effects of [Ca2+]i, SR Ca2+ load, and rest on Ca2+ spark frequency in ventricular myocytes. Am. J. Physiol. Heart Circul. Physiol. 272, H657–H668. doi: 10.1152/ajpheart.1997.272.2.H657
Shannon, T. R., Wang, F., Puglisi, J., Weber, C., and Bers, D. M. (2004). A mathematical treatment of integrated ca dynamics within the ventricular myocyte. Biophys. J. 87, 3351–3371. doi: 10.1529/biophysj.104.047449
Solovey, G., Fraiman, D., Pando, B., and Dawson, S. P. (2008). Simplified model of cytosolic Ca2+ dynamics in the presence of one or several clusters of Ca2+-release channels. Phys. Rev. E 78:041915. doi: 10.1103/PhysRevE.78.041915
Sutanto, H., van Sloun, B., Schönleitner, P., van Zandvoort, M. A., Antoons, G., and Heijman, J. (2018). The subcellular distribution of ryanodine receptors and l-type Ca2+ channels modulates Ca2+-transient properties and spontaneous Ca2+-release events in atrial cardiomyocytes. Front. Physiol. 9:1108. doi: 10.3389/fphys.2018.01108
Taggart, P., and Lab, M. (2008). Cardiac mechano-electric feedback and electrical restitution in humans. Prog. Biophys. Mol. Biol. 97, 452–460. doi: 10.1016/j.pbiomolbio.2008.02.021
Takamatsu, T., and Wier, W. (1990). Calcium waves in mammalian heart: quantification of origin, magnitude, waveform, and velocity. FASEB J. 4, 1519–1525. doi: 10.1096/fasebj.4.5.2307330
ter Keurs, H. E., and Boyden, P. A. (2007). Calcium and arrhythmogenesis. Physiol. Rev. 87, 457–506. doi: 10.1152/physrev.00011.2006
ter Keurs, H. E., Shinozaki, T., Zhang, Y. M., Zhang, M. L., Wakayama, Y., Sugai, Y., et al. (2008). Sarcomere mechanics in uniform and non-uniform cardiac muscle: a link between pump function and arrhythmias. Prog. Biophys. Mol. Biol. 97, 312–331. doi: 10.1016/j.pbiomolbio.2008.02.013
ter Keurs, H. E., Wakayama, Y., Miura, M., Shinozaki, T., Stuyvers, B. D., Boyden, P. A., et al. (2006a). Arrhythmogenic Ca2+ release from cardiac myofilaments. Prog. Biophys. Mol. Biol. 90, 151–171. doi: 10.1016/j.pbiomolbio.2005.07.002
ter Keurs, H. E., Wakayama, Y., Sugai, Y., Price, G., Kagaya, Y., Boyden, P. A., et al. (2006b). Role of sarcomere mechanics and Ca2+ overload in Ca2+ waves and arrhythmias in rat cardiac muscle. Ann. N. Y. Acad. Sci. 1080, 248–267. doi: 10.1196/annals.1380.020
ter Keurs, H. E., Zhang, Y. M., and Miura, M. (1998). Damage-induced arrhythmias: reversal of excitation-contraction coupling. Cardiovasc. Res. 40, 444–455. doi: 10.1016/S0008-6363(98)00263-6
Timmermann, V., Edwards, A. G., Wall, S. T., Sundnes, J., and McCulloch, A. D. (2019). Arrhythmogenic current generation by myofilament-triggered Ca2+ release and sarcomere heterogeneity. Biophys. J. 117, 2471–2485. doi: 10.1016/j.bpj.2019.11.009
Vogel, R., and Weingart, R. (1998). Mathematical model of vertebrate gap junctions derived from electrical measurements on homotypic and heterotypic channels. J. Physiol. 510, 177–189. doi: 10.1111/j.1469-7793.1998.177bz.x
Voigt, N., Heijman, J., Wang, Q., Chiang, D. Y., Li, N., Karck, M., et al. (2013). Cellular and molecular mechanisms of atrial arrhythmogenesis in patients with paroxysmal atrial fibrillation. Circulation 129, 145–156. doi: 10.1161/CIRCULATIONAHA.113.006641
Wakayama, Y., Miura, M., Sugai, Y., Kagaya, Y., Watanabe, J., ter Keurs, H. E. D. J., et al. (2001). Stretch and quick release of rat cardiac trabeculae accelerates Ca2+ waves and triggered propagated contractions. Am. J. Phys. 281, H2133–H2142. doi: 10.1152/ajpheart.2001.281.5.H2133
Walker, M. A., Williams, G. S., Kohl, T., Lehnart, S. E., Jafri, M. S., Greenstein, J. L., et al. (2014). Superresolution modeling of calcium release in the heart. Biophys. J. 107, 3018–3029. doi: 10.1016/j.bpj.2014.11.003
Wier, W., and Blatter, L. (1991). Ca2+-oscillations and Ca2+-waves in mammalian cardiac and vascular smooth muscle cells. Cell Calc. 12, 241–254. doi: 10.1016/0143-4160(91)90024-9
Williams, D. A., Delbridge, L. M., Cody, S. H., Harris, P., and Morgan, T. (1992). Spontaneous and propagated calcium release in isolated cardiac myocytes viewed by confocal microscopy. Am. J. Physiol. Cell Physiol. 262, C731–C742. doi: 10.1152/ajpcell.1992.262.3.C731
Keywords: mechno-electric feedback, intracellular calcium release, arrhythmia, calcium waves, computational model, rabbit, fibroblasts
Citation: Timmermann V and McCulloch AD (2020) Mechano-Electric Coupling and Arrhythmogenic Current Generation in a Computational Model of Coupled Myocytes. Front. Physiol. 11:519951. doi: 10.3389/fphys.2020.519951
Received: 13 December 2019; Accepted: 10 November 2020;
Published: 10 December 2020.
Edited by:
Leonid Katsnelson, Institute of Immunology and Physiology (RAS), RussiaReviewed by:
Jieyun Bai, Jinan University, ChinaKenneth Tran, The University of Auckland, New Zealand
Copyright © 2020 Timmermann and McCulloch. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Viviane Timmermann, viviane.timmermann@universitaets-herzzentrum.de