Variability in timing of spontaneous calcium release in the intact rat heart is determined by the time course of sarcoplasmic reticulum calcium load

Spine apparatus modulates Ca 2+ in spines through spatial localization of sources and sinks

Dendritic spines are small protrusions on dendrites in neurons and serve as sites of postsynaptic activity. Some of these spines contain smooth endoplasmic reticulum (SER), and sometimes an even further specialized SER known as the spine apparatus (SA). In this work, we developed a stochastic spatial model to investigate the role of the SER and the SA in modulating Ca 2+ dynamics. Using this model, we investigated how ryanodine receptor (RyR) localization, spine membrane geometry, and SER geometry can impact Ca 2+ transients in the spine and in the dendrite. Our simulations found that RyR opening is dependent on where it is localized in the SER and on the SER geometry. In order to maximize Ca 2+ in the dendrites (for activating clusters of spines and spine-spine communication), a laminar SA was favorable with RyRs localized in the neck region, closer to the dendrite. We also found that the presence of the SER without the laminar structure, coupled with RyR localization at the head, leads to higher Ca 2+ presence in the spine. These predictions serve as design principles for understanding how spines with an ER can regulate Ca 2+ dynamics differently from spines without ER through a combination of geometry and receptor localization.

Highlights

1RyR opening in dendritic spine ER is location dependent and spine geometry dependent. Ca 2+ buffers and SERCA can buffer against runaway potentiation of spines even when CICR is activated. RyRs located towards the ER neck allow for more Ca 2+ to reach the dendrites. RyRs located towards the spine head are favorable for increased Ca 2+ in spines.

Interaction of background Ca2+ influx, sarcoplasmic reticulum threshold and heart failure in determining propensity for Ca2+ waves in sheep heart

Ventricular arrhythmias can cause death in heart failure (HF). A trigger is the occurrence of Ca²⁺ waves which activate a Na⁺ -Ca²⁺ exchange (NCX) current, leading to delayed after-depolarisations and triggered action potentials. Waves arise when sarcoplasmic reticulum (SR) Ca²⁺ content reaches a threshold and are commonly induced experimentally by raising external Ca²⁺ , although the mechanism by which this causes waves is unclear and was the focus of this study. Intracellular Ca²⁺ was measured in voltage-clamped ventricular myocytes from both control sheep and those subjected to rapid pacing to produce HF. Threshold SR Ca²⁺ content was determined by applying caffeine (10  mM) following a wave and integrating wave and caffeine-induced NCX currents. Raising external Ca²⁺ induced waves in a greater proportion of HF cells than control. The associated increase of SR Ca²⁺ content was smaller in HF due to a lower threshold. Raising external Ca²⁺ had no effect on total influx via the L-type Ca²⁺ current, I_Ca-L , and increased efflux on NCX. Analysis of sarcolemmal fluxes revealed substantial background Ca²⁺ entry which sustains Ca²⁺ efflux during waves in the steady state. Wave frequency and background Ca²⁺ entry were decreased by Gd³⁺ or the TRPC6 inhibitor BI 749327. These agents also blocked Mn²⁺ entry. Inhibiting connexin hemi-channels, TRPC1/4/5, L-type channels or NCX had no effect on background entry. In conclusion, raising external Ca²⁺ induces waves via a background Ca²⁺ influx through TRPC6 channels. The greater propensity to waves in HF results from increased background entry and decreased threshold SR content. KEY POINTS: Heart failure is a pro-arrhythmic state and arrhythmias are a major cause of death. At the cellular level, Ca²⁺ waves resulting in delayed after-depolarisations are a key trigger of arrhythmias. Ca²⁺ waves arise when the sarcoplasmic reticulum (SR) becomes overloaded with Ca²⁺ . We investigate the mechanism by which raising external Ca²⁺ causes waves, and how this is modified in heart failure. We demonstrate that a novel sarcolemmal background Ca²⁺ influx via the TRPC6 channel is responsible for SR Ca²⁺ overload and Ca²⁺ waves. The increased propensity for Ca²⁺ waves in heart failure results from an increase of background influx, and a lower threshold SR content. The results of the present study highlight a novel mechanism by which Ca²⁺ waves may arise in heart failure, providing a basis for future work and novel therapeutic targets.

High spatial and temporal resolution Ca2+ imaging of myocardial strips from human, pig and rat

Ca²⁺ handling within cardiac myocytes underpins coordinated contractile function within the beating heart. This protocol enables high spatial and temporal Ca²⁺ imaging of ex vivo multicellular myocardial strips. The endocardial surface is retained, and strips of 150-300-µm thickness are dissected, loaded with Ca²⁺ indicators and mounted within 1.5 h. A list of the equipment and reagents used and the key methodological aspects allowing the use of this technique on strips from any chamber of the mammalian heart are described. We have successfully used this protocol on human, pig and rat biopsy samples. On use of this protocol with intact endocardial endothelium, we demonstrated that the myocytes develop asynchronous spontaneous Ca²⁺ events, which can be ablated by electrically evoked Ca²⁺ transients, and subsequently redevelop spontaneously after cessation of stimulation. This protocol thus offers a rapid and reliable method for studying the Ca²⁺ signaling underpinning cardiomyocyte contraction, in both healthy and diseased tissue.

CaMKII inhibition reduces arrhythmogenic Ca2+ events in subendocardial cryoinjured rat living myocardial slices

Spontaneous Ca²⁺ release (SCR) can cause triggered activity and initiate arrhythmias. Intrinsic transmural heterogeneities in Ca²⁺ handling and their propensity to disease remodeling may differentially modulate SCR throughout the left ventricular (LV) wall and cause transmural differences in arrhythmia susceptibility. Here, we aimed to dissect the effect of cardiac injury on SCR in different regions in the intact LV myocardium using cryoinjury on rat living myocardial slices (LMS). We studied SCR under proarrhythmic conditions using a fluorescent Ca²⁺ indicator and high-resolution imaging in LMS from the subendocardium (ENDO) and subepicardium (EPI). Cryoinjury caused structural remodeling, with loss in T-tubule density and an increased time of Ca²⁺ transients to peak after injury. In ENDO LMS, the Ca²⁺ transient amplitude and decay phase were reduced, while these were not affected in EPI LMS after cryoinjury. The frequency of spontaneous whole-slice contractions increased in ENDO LMS without affecting EPI LMS after injury. Cryoinjury caused an increase in foci that generates SCR in both ENDO and EPI LMS. In ENDO LMS, SCRs were more closely distributed and had reduced latencies after cryoinjury, whereas this was not affected in EPI LMS. Inhibition of CaMKII reduced the number, distribution, and latencies of SCR, as well as whole-slice contractions in ENDO LMS, but not in EPI LMS after cryoinjury. Furthermore, CaMKII inhibition did not affect the excitation–contraction coupling in cryoinjured ENDO or EPI LMS. In conclusion, we demonstrate increased arrhythmogenic susceptibility in the injured ENDO. Our findings show involvement of CaMKII and highlight the need for region-specific targeting in cardiac therapies.

Activation of RAGE-dependent endoplasmic reticulum stress associates with exacerbated postmyocardial infarction ventricular arrhythmias in diabetes

Association between receptor for advanced glycation end products (RAGE) and postmyocardial infarction (MI) ventricular arrhythmias (VAs) in diabetes was investigated. Correlation between premature ventricular contractions (PVCs) and serum advanced glycation end products (AGEs) content was analyzed in a cohort consisting of 101 patients with ST-segment elevated MI (STEMI). MI diabetic rats were treated with anti-receptor for AGE (RAGE) antibody. Electrocardiography was used to record VAs. Myocytes were isolated from adjacent area around infracted region. Immunofluorescent stains were used to evaluate the association between FKBP12.6 (FK506-bindingprotein 12.6) and ryanodine receptor 2 (RyR2). Calcium sparks were evaluated by confocal microscope. Protein expression and phosphorylation were assessed by Western blotting. Calcineurin (CaN) enzymatic activity and RyR2 channel activity were also determined. In the cohort study, significantly increased amount of PVC was found in STEMI patients with diabetes (P < 0.05). Serum AGE concentration was significantly positively correlated with PVC amount in patients with STEMI (r = 0.416, P < 0.001). Multivariate analysis showed that serum AGE concentration was independently and positively related to frequent PVCs (adjusted hazard ratio, 1.86; 95% CI, 1.09-3.18, P = 0.022). In the animal study, increased glucose-regulated protein 78 (GRP78) expression, protein kinase RNA-like ER kinase (PERK) phosphorylation, CaN enzymatic activity, FKBP12.6-RyR2 disassociation, RyR2 channel opening, and endoplasmic reticulum (ER) calcium releasing were found in diabetic MI animals, which were attenuated by anti-RAGE antibody treatment. This RAGE blocking also significantly lowered the VA amount in diabetic MI animals. Activation of RAGE-dependent ER stress-mediated PERK/CaN/RyR2 signaling participated in post-MI VAs in diabetes.NEW & NOTEWORTHY In this study, we proposed a possible mechanism interpreting the clinical scenario that after myocardial infarction (MI) patients were more vulnerable to ventricular arrhythmias (VAs) when complicated with diabetes. A cohort study revealed that advanced glycation end products (AGEs) accumulated in patients with diabetes and closely associated post-MI VAs. In vivo and in vitro studies indicated that receptor for AGEs (RAGE)-dependent endoplasmic reticulum (ER) stress protein kinase RNA-like ER kinase (PERK) pathway triggered VAs, via ER calcium releasing, through calcineurin/RyR2 mechanism.

Mechanistic link between CaM-RyR2 interactions and the genesis of cardiac arrhythmia

In this study, we develop a computational model of the interaction between ryanodine receptor type 2 (RyR2) and calmodulin (CaM) to explore the mechanistic link between CaM-RyR2 interactions and cardiac arrhythmia. Our starting point is a biophysically based computational model of CaM binding to a single RyR2 subunit, which reproduces single-channel RyR2 measurements in lipid bilayers. We then integrate this CaM-RyR2 model into a spatially distributed whole-cell model of Ca cycling, which is used to investigate the relationship between CaM and Ca cycling homeostasis. We show that a reduction in CaM concentration leads to a substantial increase in the rate of spontaneous Ca sparks, and this induces a marked reduction in sarcoplasmic reticulum Ca load during steady-state pacing. Also, we show that a reduction in CaM modifies the RyR2 open probability, which makes the cell more prone to Ca wave propagation. These results indicate that aberrant Ca cycling activity during pacing is determined by the interplay between sarcoplasmic reticulum load reduction and the threshold for Ca wave propagation. Based on these results, we show that when CaM is reduced, Ca waves can occur in a cell and induce action potential perturbations that are arrhythmogenic. Thus, this study outlines a novel, to our knowledge, mechanistic link between CaM-RyR2 binding kinetics and the induction of arrhythmias in the heart.

Cobalt oxide nanoparticles induce oxidative stress and alter electromechanical function in rat ventricular myocytes

Background

Nanotoxicology is an increasingly relevant field and sound paradigms on how inhaled nanoparticles (NPs) interact with organs at the cellular level, causing harmful conditions, have yet to be established. This is particularly true in the case of the cardiovascular system, where experimental and clinical evidence shows morphological and functional damage associated with NP exposure. Giving the increasing interest on cobalt oxide (Co₃O₄) NPs applications in industrial and bio-medical fields, a detailed knowledge of the involved toxicological effects is required, in view of assessing health risk for subjects/workers daily exposed to nanomaterials. Specifically, it is of interest to evaluate whether NPs enter cardiac cells and interact with cell function. We addressed this issue by investigating the effect of acute exposure to Co₃O₄-NPs on excitation-contraction coupling in freshly isolated rat ventricular myocytes.

Results

Patch clamp analysis showed instability of resting membrane potential, decrease in membrane electrical capacitance, and dose-dependent decrease in action potential duration in cardiomyocytes acutely exposed to Co₃O₄-NPs. Motion detection and intracellular calcium fluorescence highlighted a parallel impairment of cell contractility in comparison with controls. Specifically, NP-treated cardiomyocytes exhibited a dose-dependent decrease in the fraction of shortening and in the maximal rate of shortening and re-lengthening, as well as a less efficient cytosolic calcium clearing and an increased tendency to develop spontaneous twitches. In addition, treatment with Co₃O₄-NPs strongly increased ROS accumulation and induced nuclear DNA damage in a dose dependent manner. Finally, transmission electron microscopy analysis demonstrated that acute exposure did lead to cellular internalization of NPs.

Conclusions

Taken together, our observations indicate that Co₃O₄-NPs alter cardiomyocyte electromechanical efficiency and intracellular calcium handling, and induce ROS production resulting in oxidative stress that can be related to DNA damage and adverse effects on cardiomyocyte functionality.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12989-020-00396-6.

Stabilizer Cell Gene Therapy: A Less-Is-More Strategy to Prevent Cardiac Arrhythmias

BACKGROUND

In cardiac gene therapy to improve contractile function, achieving gene expression in the majority of cardiac myocytes is essential. In preventing cardiac arrhythmias, however, this goal may not be as important since transduction efficiencies as low as 40% suppressed ventricular arrhythmias in genetically modified mice with catecholaminergic polymorphic ventricular tachycardia.

METHODS

Using computational modeling, we simulated 1-, 2-, and 3-dimensional tissue under a variety of conditions to test the ability of genetically engineered nonarrhythmogenic stabilizer cells to suppress triggered activity due to delayed or early afterdepolarizations.

RESULTS

Due to source-sink relationships in cardiac tissue, a minority (20%-50%) of randomly distributed stabilizer cells engineered to be nonarrhythmogenic can suppress the ability of arrhythmogenic cells to generate delayed and early afterdepolarizations-related arrhythmias. Stabilizer cell gene therapy strategy can be designed to correct a specific arrhythmogenic mutation, as in the catecholaminergic polymorphic ventricular tachycardia mice studies, or more generally to suppress delayed or early afterdepolarizations from any cause by overexpressing the inward rectifier K channel Kir2.1 in stabilizer cells.

CONCLUSIONS

This promising antiarrhythmic strategy warrants further testing in experimental models to evaluate its clinical potential.

Arrhythmia mechanisms and spontaneous calcium release: Bi-directional coupling between re-entrant and focal excitation

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation.

Ventricular Endocardial Tissue Geometry Affects Stimulus Threshold and Effective Refractory Period

BACKGROUND

Understanding the biophysical processes by which electrical stimuli applied to cardiac tissue may result in local activation is important in both the experimental and clinical electrophysiology laboratory environments, as well as for gaining a more in-depth knowledge of the mechanisms of focal-trigger-induced arrhythmias. Previous computational models have predicted that local myocardial tissue architecture alone may significantly modulate tissue excitability, affecting both the local stimulus current required to excite the tissue and the local effective refractory period (ERP). In this work, we present experimental validation of this structural modulation of local tissue excitability on the endocardial tissue surface, use computational models to provide mechanistic understanding of this phenomena in relation to localized changes in electrotonic loading, and demonstrate its implications for the capture of afterdepolarizations.

METHODS AND RESULTS

Experiments on rabbit ventricular wedge preparations showed that endocardial ridges (surfaces of negative mean curvature) had a stimulus capture threshold that was 0.21 ± 0.03 V less than endocardial grooves (surfaces of positive mean curvature) for pairwise comparison (24% reduction, corresponding to 56.2 ± 6.4% of the energy). When stimulated at the minimal stimulus strength for capture, ridge locations showed a shorter ERP than grooves (n = 6, mean pairwise difference 7.4 ± 4.2 ms). When each site was stimulated with identical-strength stimuli, the difference in ERP was further increased (mean pairwise difference 15.8 ± 5.3 ms). Computational bidomain models of highly idealized cylindrical endocardial structures qualitatively agreed with these findings, showing that such changes in excitability are driven by structural modulation in electrotonic loading, quantifying this relationship as a function of surface curvature. Simulations further showed that capture of delayed afterdepolarizations was more likely in trabecular ridges than grooves, driven by this difference in loading.

CONCLUSIONS

We have demonstrated experimentally and explained mechanistically in computer simulations that the ability to capture tissue on the endocardial surface depends upon the local tissue architecture. These findings have important implications for deepening our understanding of excitability differences related to anatomical structure during stimulus application that may have important applications in the translation of novel experimental optogenetics pacing strategies. The uncovered preferential vulnerability to capture of afterdepolarizations of endocardial ridges, compared to grooves, provides important insight for understanding the mechanisms of focal-trigger-induced arrhythmias.

Spontaneous Ca2+ transients in rat pulmonary vein cardiomyocytes are increased in frequency and become more synchronous following electrical stimulation

The pulmonary veins have an external sleeve of cardiomyocytes that are a widely recognised source of ectopic electrical activity that can lead to atrial fibrillation. Although the mechanisms behind this activity are currently unknown, changes in intracellular calcium (Ca²⁺) signalling are purported to play a role. Therefore, the intracellular Ca²⁺ concentration was monitored in the pulmonary vein using fluo-4 and epifluorescence microscopy. Electrical field stimulation evoked a synchronous rise in Ca²⁺ in neighbouring cardiomyocytes; asynchronous spontaneous Ca²⁺ transients between electrical stimuli were also present. Immediately following termination of electrical field stimulation at 3 Hz or greater, the frequency of the spontaneous Ca²⁺ transients was increased from 0.45 ± 0.06 Hz under basal conditions to between 0.59 ± 0.05 and 0.65 ± 0.06 Hz (P < 0.001). Increasing the extracellular Ca²⁺ concentration enhanced this effect, with the frequency of spontaneous Ca²⁺ transients increasing from 0.45 ± 0.05 Hz to between 0.75 ± 0.06 and 0.94 ± 0.09 Hz after electrical stimulation at 3 to 9 Hz (P < 0.001), and this was accompanied by a significant increase in the velocity of Ca²⁺ transients that manifested as waves. Moreover, in the presence of high extracellular Ca²⁺, the spontaneous Ca²⁺ transients occurred more synchronously in the initial few seconds following electrical stimulation. The ryanodine receptors, which are the source of spontaneous Ca²⁺ transients in pulmonary vein cardiomyocytes, were found to be arranged in a striated pattern in the cell interior, as well as along the periphery of cell. Furthermore, labelling the sarcolemma with di-4-ANEPPS showed that over 90% of pulmonary vein cardiomyocytes possessed T-tubules. These findings demonstrate that the frequency of spontaneous Ca²⁺ transients in the rat pulmonary vein are increased following higher rates of electrical stimulation and increasing the extracellular Ca²⁺ concentration.

Synchronization of Intracellular Ca2+ Release in Multicellular Cardiac Preparations

In myocardial tissue, Ca²⁺ release from the sarcoplasmic reticulum (SR) that occurs via the ryanodine receptor (RyR2) channel complex. Ca²⁺ release through RyR2 can be either stimulated by an action potential (AP) or spontaneous. The latter is often associated with triggered afterdepolarizations, which in turn may lead to sustained arrhythmias. It is believed that some synchronization mechanism exists for afterdepolarizations and APs in neighboring myocytes, possibly a similarly timed recovery of RyR2 from refractoriness, which enables RyR2s to reach the threshold for spontaneous Ca²⁺ release simultaneously. To investigate this synchronization mechanism in absence of genetic factors that predispose arrhythmia, we examined the generation of triggered activity in multicellular cardiac preparations. In myocardial trabeculae from the rat, we demonstrated that in the presence of both isoproterenol and caffeine, neighboring myocytes within the cardiac trabeculae were able to synchronize their diastolic spontaneous SR Ca²⁺ release. Using confocal Ca²⁺ imaging, we could visualize Ca²⁺ waves in the multicellular preparation, while these waves were not always present in every myocyte within the trabeculae, we observed that, over time, the Ca²⁺ waves can synchronize in multiple myocytes. This synchronized activity was sufficiently strong that it could trigger a synchronized, propagated contraction in the whole trabecula encompassing even previously quiescent myocytes. The detection of Ca²⁺ dynamics in individual myocytes in their in situ setting at the multicellular level exposed a synchronization mechanism that could induce local triggered activity in the heart in the absence of global Ca²⁺ dysregulation.

In Vivo Calcium Imaging of Cardiomyocytes in the Beating Mouse Heart With Multiphoton Microscopy

Background: Understanding the microscopic dynamics of the beating heart has been challenging due to the technical nature of imaging with micrometer resolution while the heart moves. The development of multiphoton microscopy has made in vivo, cell-resolved measurements of calcium dynamics and vascular function possible in motionless organs such as the brain. In heart, however, studies of in vivo interactions between cells and the native microenvironment are behind other organ systems. Our goal was to develop methods for intravital imaging of cardiac structural and calcium dynamics with microscopic resolution.

Methods: Ventilated mice expressing GCaMP6f, a genetically encoded calcium indicator, received a thoracotomy to provide optical access to the heart. Vasculature was labeled with an injection of dextran-labeled dye. The heart was partially stabilized by a titanium probe with a glass window. Images were acquired at 30 frames per second with spontaneous heartbeat and continuously running, ventilated breathing. The data were reconstructed into three-dimensional volumes showing tissue structure, vasculature, and GCaMP6f signal in cardiomyocytes as a function of both the cardiac and respiratory cycle.

Results: We demonstrated the capability to simultaneously measure calcium transients, vessel size, and tissue displacement in three dimensions with micrometer resolution. Reconstruction at various combinations of cardiac and respiratory phase enabled measurement of regional and single-cell cardiomyocyte calcium transients (GCaMP6f fluorescence). GCaMP6f fluorescence transients in individual, aberrantly firing cardiomyocytes were also quantified. Comparisons of calcium dynamics (rise-time and tau) at varying positions within the ventricle wall showed no significant depth dependence.

Conclusion: This method enables studies of coupling between contraction and excitation during physiological blood perfusion and breathing at high spatiotemporal resolution. These capabilities could lead to a new understanding of normal and disease function of cardiac cells.

Regional acidosis locally inhibits but remotely stimulates Ca2+ waves in ventricular myocytes

Aims

Spontaneous Ca²⁺ waves in cardiomyocytes are potentially arrhythmogenic. A powerful controller of Ca²⁺ waves is the cytoplasmic H⁺ concentration ([H⁺]_i), which fluctuates spatially and temporally in conditions such as myocardial ischaemia/reperfusion. H⁺-control of Ca²⁺ waves is poorly understood. We have therefore investigated how [H⁺]_i co-ordinates their initiation and frequency.

Methods and results

Spontaneous Ca²⁺ waves were imaged (fluo-3) in rat isolated ventricular myocytes, subjected to modest Ca²⁺-overload. Whole-cell intracellular acidosis (induced by acetate-superfusion) stimulated wave frequency. Pharmacologically blocking sarcolemmal Na⁺/H⁺ exchange (NHE1) prevented this stimulation, unveiling inhibition by H⁺. Acidosis also increased Ca²⁺ wave velocity. Restricting acidosis to one end of a myocyte, using a microfluidic device, inhibited Ca²⁺ waves in the acidic zone (consistent with ryanodine receptor inhibition), but stimulated wave emergence elsewhere in the cell. This remote stimulation was absent when NHE1 was selectively inhibited in the acidic zone. Remote stimulation depended on a locally evoked, NHE1-driven rise of [Na⁺]_i that spread rapidly downstream.

Conclusion

Acidosis influences Ca²⁺ waves via inhibitory Hi+ and stimulatory Nai+ signals (the latter facilitating intracellular Ca²⁺-loading through modulation of sarcolemmal Na⁺/Ca²⁺ exchange activity). During spatial [H⁺]_i-heterogeneity, Hi+-inhibition dominates in acidic regions, while rapid Nai+ diffusion stimulates waves in downstream, non-acidic regions. Local acidosis thus simultaneously inhibits and stimulates arrhythmogenic Ca²⁺-signalling in the same myocyte. If the principle of remote H⁺-stimulation of Ca²⁺ waves also applies in multicellular myocardium, it raises the possibility of electrical disturbances being driven remotely by adjacent ischaemic areas, which are known to be intensely acidic.

Estimating the probabilities of rare arrhythmic events in multiscale computational models of cardiac cells and tissue

Ectopic heartbeats can trigger reentrant arrhythmias, leading to ventricular fibrillation and sudden cardiac death. Such events have been attributed to perturbed Ca2+ handling in cardiac myocytes leading to spontaneous Ca2+ release and delayed afterdepolarizations (DADs). However, the ways in which perturbation of specific molecular mechanisms alters the probability of ectopic beats is not understood. We present a multiscale model of cardiac tissue incorporating a biophysically detailed three-dimensional model of the ventricular myocyte. This model reproduces realistic Ca2+ waves and DADs driven by stochastic Ca2+ release channel (RyR) gating and is used to study mechanisms of DAD variability. In agreement with previous experimental and modeling studies, key factors influencing the distribution of DAD amplitude and timing include cytosolic and sarcoplasmic reticulum Ca2+ concentrations, inwardly rectifying potassium current (IK1) density, and gap junction conductance. The cardiac tissue model is used to investigate how random RyR gating gives rise to probabilistic triggered activity in a one-dimensional myocyte tissue model. A novel spatial-average filtering method for estimating the probability of extreme (i.e. rare, high-amplitude) stochastic events from a limited set of spontaneous Ca2+ release profiles is presented. These events occur when randomly organized clusters of cells exhibit synchronized, high amplitude Ca2+ release flux. It is shown how reduced IK1 density and gap junction coupling, as observed in heart failure, increase the probability of extreme DADs by multiple orders of magnitude. This method enables prediction of arrhythmia likelihood and its modulation by alterations of other cellular mechanisms.

Calcium-Dependent Arrhythmogenic Foci Created by Weakly Coupled Myocytes in the Failing Heart

RATIONALE

Intercellular uncoupling and Ca²⁺ (Ca) mishandling can initiate triggered ventricular arrhythmias. Spontaneous Ca release activates inward current which depolarizes membrane potential (V_m) and can trigger action potentials in isolated myocytes. However, cell-cell coupling in intact hearts limits local depolarization and may protect hearts from this arrhythmogenic mechanism. Traditional optical mapping lacks the spatial resolution to assess coupling of individual myocytes.

OBJECTIVE

We investigate local intercellular coupling in Ca-induced depolarization in intact hearts, using confocal microscopy to measure local V_m and intracellular [Ca] simultaneously.

METHODS AND RESULTS

We used isolated Langendorff-perfused hearts from control (CTL) and heart failure (HF) mice (HF induced by transaortic constriction). In CTL hearts, 1.4% of myocytes were poorly synchronized with neighboring cells and exhibited asynchronous (AS) Ca transients. These AS myocytes were much more frequent in HF (10.8% of myocytes, P<0.05 versus CTL). Local Ca waves depolarized V_m in HF but not CTL hearts, suggesting weaker gap junction coupling in HF-AS versus CTL-AS myocytes. Cell-cell coupling was assessed by calcein fluorescence recovery after photobleach during intracellular [Ca] recording. All regions in CTL hearts exhibited faster calcein diffusion than in HF, with HF-AS myocyte being slowest. In HF-AS, enhancing gap junction conductance (with rotigaptide) increased coupling and suppressed V_m depolarization during Ca waves. Conversely, in CTL hearts, gap junction inhibition (carbenoxolone) decreased coupling and allowed Ca wave-induced depolarizations. Synchronization of Ca wave initiation and triggered action potentials were observed in HF hearts and computational models.

CONCLUSIONS

Well-coupled CTL myocytes are effectively voltage-clamped during Ca waves, protecting the heart from triggered arrhythmias. Spontaneous Ca waves are much more common in HF myocytes and these AS myocytes are also poorly coupled, enabling local Ca-induced inward current of sufficient source strength to overcome a weakened current sink to depolarize V_m and trigger action potentials.

A model of cardiac ryanodine receptor gating predicts experimental Ca2+-dynamics and Ca2+-triggered arrhythmia in the long QT syndrome

Abnormal Ca²⁺ handling is well-established as the trigger of cardiac arrhythmia in catecholaminergic polymorphic ventricular tachycardia and digoxin toxicity, but its role remains controversial in Torsade de Pointes (TdP), the arrhythmia associated with the long QT syndrome (LQTS). Recent experimental results show that early afterdepolarizations (EADs) that initiate TdP are caused by spontaneous (non-voltage-triggered) Ca²⁺ release from Ca²⁺-overloaded sarcoplasmic reticulum (SR) rather than the activation of the L-type Ca²⁺-channel window current. In bradycardia and long QT type 2 (LQT2), a second, non-voltage triggered cytosolic Ca²⁺ elevation increases gradually in amplitude, occurs before overt voltage instability, and then precedes the rise of EADs. Here, we used a modified Shannon-Puglisi-Bers model of rabbit ventricular myocytes to reproduce experimental Ca²⁺ dynamics in bradycardia and LQT2. Abnormal systolic Ca²⁺-oscillations and EADs caused by SR Ca²⁺-release are reproduced in a modified 0-dimensional model, where 3 gates in series control the ryanodine receptor (RyR2) conductance. Two gates control RyR2 activation and inactivation and sense cytosolic Ca²⁺ while a third gate senses luminal junctional SR Ca²⁺. The model predicts EADs in bradycardia and low extracellular [K⁺] and cessation of SR Ca²⁺-release terminate salvos of EADs. Ca²⁺-waves, systolic cell-synchronous Ca²⁺-release, and multifocal diastolic Ca²⁺ release seen in subcellular Ca²⁺-mapping experiments are observed in the 2-dimensional version of the model. These results support the role of SR Ca²⁺-overload, abnormal SR Ca²⁺-release, and the subsequent activation of the electrogenic Na⁺/Ca²⁺-exchanger as the mechanism of TdP. The model offers new insights into the genesis of cardiac arrhythmia and new therapeutic strategies.

Stochastic spontaneous calcium release events and sodium channelopathies promote ventricular arrhythmias

Premature ventricular complexes (PVCs), the first initiating beats of a variety of cardiac arrhythmias, have been associated with spontaneous calcium release (SCR) events at the cell level. However, the mechanisms underlying the degeneration of such PVCs into arrhythmias are not fully understood. The objective of this study was to investigate the conditions under which SCR-mediated PVCs can lead to ventricular arrhythmias. In particular, we sought to determine whether sodium (Na⁺) current loss-of-function in the structurally normal ventricles provides a substrate for unidirectional conduction block and reentry initiated by SCR-mediated PVCs. To achieve this goal, a stochastic model of SCR was incorporated into an anatomically accurate compute model of the rabbit ventricles with the His-Purkinje system (HPS). Simulations with reduced Na⁺ current due to a negative-shift in the steady-state channel inactivation showed that SCR-mediated delayed afterdepolarizations led to PVC formation in the HPS, where the electrotonic load was lower, conduction block, and reentry in the 3D myocardium. Moreover, arrhythmia initiation was only possible when intrinsic electrophysiological heterogeneity in action potential within the ventricles was present. In conclusion, while benign in healthy individuals SCR-mediated PVCs can lead to life-threatening ventricular arrhythmias when combined with Na⁺ channelopathies.

The role of luminal Ca regulation in Ca signaling refractoriness and cardiac arrhythmogenesis

Györke et al. discuss the role of sarcoplasmic reticulum Ca2+ in cardiac refractoriness and pathological implications.

Multiscale Determinants of Delayed Afterdepolarization Amplitude in Cardiac Tissue

Spontaneous calcium (Ca) waves in cardiac myocytes underlie delayed afterdepolarizations (DADs) that trigger cardiac arrhythmias. How these subcellular/cellular events overcome source-sink factors in cardiac tissue to generate DADs of sufficient amplitude to trigger action potentials is not fully understood. Here, we evaluate quantitatively how factors at the subcellular scale (number of Ca wave initiation sites), cellular scale (sarcoplasmic reticulum (SR) Ca load), and tissue scale (synchrony of Ca release in populations of myocytes) determine DAD features in cardiac tissue using a combined experimental and computational modeling approach. Isolated patch-clamped rabbit ventricular myocytes loaded with Fluo-4 to image intracellular Ca were rapidly paced during exposure to elevated extracellular Ca (2.7 mmol/L) and isoproterenol (0.25 μmol/L) to induce diastolic Ca waves and subthreshold DADs. As the number of paced beats increased from 1 to 5, SR Ca content (assessed with caffeine pulses) increased, the number of Ca wave initiation sites increased, integrated Ca transients and DADs became larger and shorter in duration, and the latency period to the onset of Ca waves shortened with reduced variance. In silico analysis using a computer model of ventricular tissue incorporating these experimental measurements revealed that whereas all of these factors promoted larger DADs with higher probability of generating triggered activity, the latency period variance and SR Ca load had the greatest influences. Therefore, incorporating quantitative experimental data into tissue level simulations reveals that increased intracellular Ca promotes DAD-mediated triggered activity in tissue predominantly by increasing both the synchrony (decreasing latency variance) of Ca waves in nearby myocytes and SR Ca load, whereas the number of Ca wave initiation sites per myocyte is less important.

Stochastic initiation and termination of calcium-mediated triggered activity in cardiac myocytes

Cardiac myocytes normally initiate action potentials in response to a current stimulus that depolarizes the membrane above an excitation threshold. Aberrant excitation can also occur due to spontaneous calcium (Ca²⁺) release (SCR) from intracellular stores after the end of a preceding action potential. SCR drives the Na⁺/Ca²⁺ exchange current inducing a "delayed afterdepolarization" that can in turn trigger an action potential if the excitation threshold is reached. This "triggered activity" is known to cause arrhythmias, but how it is initiated and terminated is not understood. Using computer simulations of a ventricular myocyte model, we show that initiation and termination are inherently random events. We determine the probability of those events from statistical measurements of the number of beats before initiation and before termination, respectively, which follow geometric distributions. Moreover, we elucidate the origin of randomness by a statistical analysis of SCR events, which do not follow a Poisson process observed in other eukaryotic cells. Due to synchronization of Ca²⁺ releases during the action potential upstroke, waiting times of SCR events after the upstroke are narrowly distributed, whereas SCR amplitudes follow a broad normal distribution with a width determined by fluctuations in the number of independent Ca²⁺ wave foci. This distribution enables us to compute the probabilities of initiation and termination of bursts of triggered activity that are maintained by a positive feedback between the action potential upstroke and SCR. Our results establish a theoretical framework for interpreting complex and varied manifestations of triggered activity relevant to cardiac arrhythmias.

Regional distribution of T-tubule density in left and right atria in dogs

BACKGROUND

The peculiarities of transverse tubule (T-tubule) morphology and distribution in the atrium-and how they contribute to excitation-contraction coupling-are just beginning to be understood.

OBJECTIVES

The objectives of this study were to determine T-tubule density in the intact, live right and left atria in a large animal and to determine intraregional differences in T-tubule organization within each atrium.

METHODS

Using confocal microscopy, T-tubules were imaged in both atria in intact, Langendorf-perfused normal dog hearts loaded with di-4-ANEPPS. T-tubules were imaged in large populations of myocytes from the endocardial surface of each atrium. Computerized data analysis was performed using a new MatLab (Mathworks, Natick, MA) routine, AutoTT.

RESULTS

There was a large percentage of myocytes that had no T-tubules in both atria with a higher percentage in the right atrium (25.1%) than in the left atrium (12.5%) (P < .02). The density of transverse and longitudinal T-tubule elements was low in cells that did contain T-tubules, but there were no significant differences in density between the left atrial appendage, the pulmonary vein-posterior left atrium, the right atrial appendage, and the right atrial free wall. In contrast, there were significant differences in sarcomere spacing and cell width between different regions of the atria.

CONCLUSION

There is a sparse T-tubule network in atrial myocytes throughout both dog atria, with significant numbers of myocytes in both atria-the right atrium more so than the left atrium-having no T-tubules at all. These regional differences in T-tubule distribution, along with differences in cell width and sarcomere spacing, may have implications for the emergence of substrate for atrial fibrillation.

Calcium-voltage coupling in the genesis of early and delayed afterdepolarizations in cardiac myocytes

Early afterdepolarizations (EADs) and delayed afterdepolarizations (DADs) are voltage oscillations known to cause cardiac arrhythmias. EADs are mainly driven by voltage oscillations in the repolarizing phase of the action potential (AP), while DADs are driven by spontaneous calcium (Ca) release during diastole. Because voltage and Ca are bidirectionally coupled, they modulate each other's behaviors, and new AP and Ca cycling dynamics can emerge from this coupling. In this study, we performed computer simulations using an AP model with detailed spatiotemporal Ca cycling incorporating stochastic openings of Ca channels and ryanodine receptors to investigate the effects of Ca-voltage coupling on EAD and DAD dynamics. Simulations were complemented by experiments in mouse ventricular myocytes. We show that: 1) alteration of the Ca transient due to increased ryanodine receptor leakiness and/or sarco/endoplasmic reticulum Ca ATPase activity can either promote or suppress EADs due to the complex effects of Ca on ionic current properties; 2) spontaneous Ca waves also exhibit complex effects on EADs, but cannot induce EADs of significant amplitude without the participation of ICa,L; 3) lengthening AP duration and the occurrence of EADs promote DADs by increasing intracellular Ca loading, and two mechanisms of DADs are identified, i.e., Ca-wave-dependent and Ca-wave-independent; and 4) Ca-voltage coupling promotes complex EAD patterns such as EAD alternans that are not observed for solely voltage-driven EADs. In conclusion, Ca-voltage coupling combined with the nonlinear dynamical behaviors of voltage and Ca cycling play a key role in generating complex EAD and DAD dynamics observed experimentally in cardiac myocytes, whose mechanisms are complex but analyzable.

The link between abnormal calcium handling and electrical instability in acquired long QT syndrome--Does calcium precipitate arrhythmic storms?

Release of Ca(2+) ions from sarcoplasmic reticulum (SR) into myocyte cytoplasm and their binding to troponin C is the final signal form myocardial contraction. Synchronous contraction of ventricular myocytes is necessary for efficient cardiac pumping function. This requires both shuttling of Ca(2+) between SR and cytoplasm in individual myocytes, and organ-level synchronization of this process by means of electrical coupling among ventricular myocytes. Abnormal Ca(2+) release from SR causes arrhythmias in the setting of CPVT (catecholaminergic polymorphic ventricular tachycardia) and digoxin toxicity. Recent optical mapping data indicate that abnormal Ca(2+) handling causes arrhythmias in models of both repolarization impairment and profound bradycardia. The mechanisms involve dynamic spatial heterogeneity of myocardial Ca(2+) handling preceding arrhythmia onset, cell-synchronous systolic secondary Ca(2+) elevation (SSCE), as well as more complex abnormalities of intracellular Ca(2+) handling detected by subcellular optical mapping in Langendorff-perfused hearts. The regional heterogeneities in Ca(2+) handling cause action potential (AP) heterogeneities through sodium-calcium exchange (NCX) activation and eventually overwhelm electrical coupling of the tissue. Divergent Ca(2+) dynamics among different myocardial regions leads to temporal instability of AP duration and - on the patient level - in T wave lability. Although T-wave alternans has been linked to cardiac arrhythmias, non-alternans lability is observed in pre-clinical models of the long QT syndrome (LQTS) and CPVT, and in LQTS patients. Analysis of T wave lability may provide a real-time window on the abnormal Ca(2+) dynamics causing specific arrhythmias such as Torsade de Pointes (TdP).

An integrated finite element simulation of cardiomyocyte function based on triphasic theory

In numerical simulations of cardiac excitation-contraction coupling, the intracellular potential distribution and mobility of cytosol and ions have been mostly ignored. Although the intracellular potential gradient is small, during depolarization it can be a significant driving force for ion movement, and is comparable to diffusion in terms of net flux. Furthermore, fluid in the t-tubules is thought to advect ions to facilitate their exchange with the extracellular space. We extend our previous finite element model that was based on triphasic theory to examine the significance of these factors in cardiac physiology. Triphasic theory allows us to study the behavior of solids (proteins), fluids (cytosol) and ions governed by mechanics and electrochemistry in detailed subcellular structures, including myofibrils, mitochondria, the sarcoplasmic reticulum, membranes, and t-tubules. Our simulation results predicted an electrical potential gradient inside the t-tubules at the onset of depolarization, which corresponded to the Na(+) channel distribution therein. Ejection and suction of fluid between the t-tubules and the extracellular compartment during isometric contraction were observed. We also examined the influence of t-tubule morphology and mitochondrial location on the electrophysiology and mechanics of the cardiomyocyte. Our results confirm that the t-tubule structure is important for synchrony of Ca(2+) release, and suggest that mitochondria in the sub-sarcolemmal region might serve to cancel Ca(2+) inflow through surface sarcolemma, thereby maintaining the intracellular Ca(2+) environment in equilibrium.

Delayed afterdepolarizations generate both triggers and a vulnerable substrate promoting reentry in cardiac tissue

BACKGROUND

Delayed afterdepolarizations (DADs) have been well characterized as arrhythmia triggers, but their role in generating a tissue substrate vulnerable to reentry is not well understood.

OBJECTIVE

The purpose of this study was to test the hypothesis that random DADs can self-organize to generate both an arrhythmia trigger and a vulnerable substrate simultaneously in cardiac tissue as a result of gap junction coupling.

METHODS

Computer simulations in 1-dimensional cable and 2-dimensional tissue models were performed. The cellular DAD amplitude was varied by changing the strength of sarcoplasmic reticulum calcium release. Random DAD latency and amplitude in different cells were simulated using gaussian distributions.

RESULTS

Depending on the strength of spontaneous sarcoplasmic reticulum calcium release and other conditions, random DADs in cardiac tissue resulted in the following behaviors: (1) triggered activity (TA); (2) a vulnerable tissue substrate causing unidirectional conduction block and reentry by inactivating sodium channels; (3) both triggers and a vulnerable substrate simultaneously by generating TA in regions next to regions with subthreshold DADs susceptible to unidirectional conduction block and reentry. The probability of the latter 2 behaviors was enhanced by reduced sodium channel availability, reduced gap junction coupling, increased tissue heterogeneity, and less synchronous DAD latency.

CONCLUSION

DADs can self-organize in tissue to generate arrhythmia triggers, a vulnerable tissue substrate, and both simultaneously. Reduced sodium channel availability and gap junction coupling potentiate this mechanism of arrhythmias, which are relevant to a variety of heart disease conditions.

Synchronous systolic subcellular Ca2+-elevations underlie ventricular arrhythmia in drug-induced long QT type 2

BACKGROUND

Repolarization delay is a common clinical problem, which can promote ventricular arrhythmias. In myocytes, abnormal sarcoplasmic reticulum Ca(2+)-release is proposed as the mechanism that causes early afterdepolarizations, the cellular equivalent of ectopic-activity in drug-induced long-QT syndrome. A crucial missing link is how such a stochastic process can overcome the source-sink mismatch to depolarize sufficient ventricular tissue to initiate arrhythmias.

METHODS AND RESULTS

Optical maps of action potentials and Ca(2+)-transients from Langendorff rabbit hearts were measured at low (150×150 μm(2)/pixel) and high (1.5×1.5 μm(2)/pixel) resolution before and during arrhythmias. Drug-induced long QT type 2, elicited with dofetilide inhibition of IKr (the rapid component of rectifying K+ current), produced spontaneous Ca(2+)-elevations during diastole and systole, before the onset of arrhythmias. Diastolic Ca(2+-)waves appeared randomly, propagated within individual myocytes, were out-of-phase with adjacent myocytes, and often died-out. Systolic secondary Ca(2+-)elevations were synchronous within individual myocytes, appeared 188±30 ms after the action potential-upstroke, occurred during high cytosolic Ca(2+) (40%-60% of peak-Ca(2+)-transients), appeared first in small islands (0.5×0.5 mm(2)) that enlarged and spread throughout the epicardium. Synchronous systolic Ca(2+-)elevations preceded voltage-depolarizations (9.2±5 ms; n=5) and produced pronounced Spatial Heterogeneities of Ca(2+)-transient-durations and action potential-durations. Early afterdepolarizations originating from sites with the steepest gradients of membrane-potential propagated and initiated arrhythmias. Interestingly, more complex subcellular Ca(2+)-dynamics (multiple chaotic Ca(2+)-waves) occurred during arrhythmias. K201, a ryanodine receptor stabilizer, eliminated Ca(2+)-elevations and arrhythmias.

CONCLUSIONS

The results indicate that systolic and diastolic Ca(2+)-elevations emanate from sarcoplasmic reticulum Ca(2+)-release and systolic Ca(2+)-elevations are synchronous because of high cytosolic and luminal-sarcoplasmic reticulum Ca(2+), which overcomes source-sink mismatch to trigger arrhythmias in intact hearts.

Stochastic spontaneous calcium release events trigger premature ventricular complexes by overcoming electrotonic load

AIMS

Premature ventricular complexes (PVCs) due to spontaneous calcium (Ca) release (SCR) events at the cell level can precipitate ventricular arrhythmias. However, the mechanistic link between SCRs and PVC formation remains incompletely understood. The aim of this study was to investigate the conditions under which delayed afterdepolarizations resulting from stochastic subcellular SCR events can overcome electrotonic source-sink mismatch, leading to PVC initiation.

METHODS AND RESULTS

A stochastic subcellular-scale mathematical model of SCR was incorporated in a realistic model of the rabbit ventricles and Purkinje system (PS). Elevated levels of diastolic sarcoplasmic reticulum Ca(2+) (CaSR) were imposed until triggered activity was observed, allowing us to compile statistics on probability, timing, and location of PVCs. At CaSR≥ 1500 µmol/L PVCs originated in the PS. When SCR was incapacitated in the PS, PVCs also emerged in the ventricles, but at a higher CaSR (≥1550 µmol/L) and with longer waiting times. For each model configuration tested, the probability of PVC occurrence increased from 0 to 100% within a well-defined critical CaSR range; this transition was much more abrupt in organ-scale models (∼50 µmol/L CaSR range) than in the tissue strand (∼100 µmol/L) or single-cell (∼450 µmol/L) models. Among PVCs originating in the PS, ∼68% were located near Purkinje-ventricular junctions (<1 mm).

CONCLUSION

SCR events overcome source-sink mismatch to trigger PVCs at a critical CaSR threshold. Above this threshold, PVCs emerge due to increased probability and reduced variability in timing of SCR events, leading to significant diastolic depolarization. Sites of lower electronic load, such as the PS, are preferential locations for triggering.

Approximate analytical solutions for excitation and propagation in cardiac tissue

It is well known that a variety of cardiac arrhythmias are initiated by a focal excitation in heart tissue. At the single cell level these currents are typically induced by intracellular processes such as spontaneous calcium release (SCR). However, it is not understood how the size and morphology of these focal excitations are related to the electrophysiological properties of cardiac cells. In this paper a detailed physiologically based ionic model is analyzed by projecting the excitation dynamics to a reduced one-dimensional parameter space. Based on this analysis we show that the inward current required for an excitation to occur is largely dictated by the voltage dependence of the inward rectifier potassium current (I(K1)), and is insensitive to the detailed properties of the sodium current. We derive an analytical expression relating the size of a stimulus and the critical current required to induce a propagating action potential (AP), and argue that this relationship determines the necessary number of cells that must undergo SCR in order to induce ectopic activity in cardiac tissue. Finally, we show that, once a focal excitation begins to propagate, its propagation characteristics, such as the conduction velocity and the critical radius for propagation, are largely determined by the sodium and gap junction currents with a substantially lesser effect due to repolarizing potassium currents. These results reveal the relationship between ion channel properties and important tissue scale processes such as excitation and propagation.

Perspective: a dynamics-based classification of ventricular arrhythmias

Despite key advances in the clinical management of life-threatening ventricular arrhythmias, culminating with the development of implantable cardioverter-defibrillators and catheter ablation techniques, pharmacologic/biologic therapeutics have lagged behind. The fundamental issue is that biological targets are molecular factors. Diseases, however, represent emergent properties at the scale of the organism that result from dynamic interactions between multiple constantly changing molecular factors. For a pharmacologic/biologic therapy to be effective, it must target the dynamic processes that underlie the disease. Here we propose a classification of ventricular arrhythmias that is based on our current understanding of the dynamics occurring at the subcellular, cellular, tissue and organism scales, which cause arrhythmias by simultaneously generating arrhythmia triggers and exacerbating tissue vulnerability. The goal is to create a framework that systematically links these key dynamic factors together with fixed factors (structural and electrophysiological heterogeneity) synergistically promoting electrical dispersion and increased arrhythmia risk to molecular factors that can serve as biological targets. We classify ventricular arrhythmias into three primary dynamic categories related generally to unstable Ca cycling, reduced repolarization, and excess repolarization, respectively. The clinical syndromes, arrhythmia mechanisms, dynamic factors and what is known about their molecular counterparts are discussed. Based on this framework, we propose a computational-experimental strategy for exploring the links between molecular factors, fixed factors and dynamic factors that underlie life-threatening ventricular arrhythmias. The ultimate objective is to facilitate drug development by creating an in silico platform to evaluate and predict comprehensively how molecular interventions affect not only a single targeted arrhythmia, but all primary arrhythmia dynamics categories as well as normal cardiac excitation-contraction coupling.

Distinguishing between overdrive excited and suppressed ventricular beats in guinea pig ventricular myocardium

Rapid ventricular pacing rates induces two types of beats following pacing cessation: recovery cycle length (RCL) prolongation (overdrive suppression) and RCL shortening (overdrive excitation). The goals of this study were to compare common experimental protocols for studying triggered activity in whole-heart preparations and differentiate between recovery beats using a new methodology. Post-pacing recovery beat cycle length (RCL) and QRS were normalized to pre-paced R-R and QRS intervals and analyzed using a K-means clustering algorithm. Control hearts only produced suppressed beats: RCL ratio increased with rapid pacing (25 ± 4.0%, n = 10) without changing QRS duration. Rapid pacing during hypercalcemia + hypothermia (5.5 mM and 34°C) produced significantly earlier excited beats (53 ± 14%, n = 5) with wider QRS durations (58 ± 6.3%, n = 5) than suppressed beats. Digoxin + hypothermia (0.75 μM) produced the most excited beats with significantly earlier RCL (44 ± 3.2%, n = 6) and wider QRS (60 ± 3.1%, n = 6) ratios relative to suppressed beats. Increasing pacing further shortened RCL (30 ± 7.8%, n = 6). In a prospective study, TTX (100 nM) increased RCL ratio (15 ± 6.0%, n = 10) without changing the QRS duration of excited beats. The algorithm was compared to a cross-correlation analysis with 93% sensitivity and 94% specificity. This ECG based algorithm distinguishes between triggered and automatic activity.

Depolarization of cardiac membrane potential synchronizes calcium sparks and waves in tissue

The diastolic membrane potential (Vm) can be hyperpolarized or depolarized by various factors such as hyperkalemia or hypokalemia in the long term, or by delayed afterdepolarizations in the short term. In this study, we investigate how Vm affects Ca sparks and waves. We use a physiologically detailed mathematical model to investigate individual factors that affect Ca spark generation and wave propagation. We focus on the voltage range of -90 ∼ -70 mV, which is just below the Vm for sodium channel activation. We find that Vm depolarization promotes Ca wave propagation and hyperpolarization prevents it. This finding is directly validated in voltage clamp experiments with Ca waves using isolated rat ventricular myocytes. Ca transport by the sodium-calcium exchanger (NCX) is determined by Vm as well as Na and Ca concentrations. Depolarized Vm reduces NCX-mediated efflux, elevating [Ca]i, and thus promoting Ca wave propagation. Moreover, depolarized Vm promotes spontaneous Ca releases that can cause initiation of multiple Ca waves. This indicates that during delayed afterdepolarizations, Ca release units (CRUs) interact with not just the immediately adjacent CRUs via Ca diffusion, but also further CRUs via fast (∼0.1 ms) changes in Vm mediated by the voltage and Ca-sensitive NCX. This may contribute significantly to synchronization of Ca waves among multiple cells in tissue.

Variations in local calcium signaling in adjacent cardiac myocytes of the intact mouse heart detected with two-dimensional confocal microscopy

Dyssynchronous local Ca release within individual cardiac myocytes has been linked to cellular contractile dysfunction. Differences in Ca kinetics in adjacent cells may also provide a substrate for inefficient contraction and arrhythmias. In a new approach we quantify variation in local Ca transients between adjacent myocytes in the whole heart. Langendorff-perfused mouse hearts were loaded with Fluo-8 AM to detect Ca and Di-4-ANEPPS to visualize cell membranes. A spinning disc confocal microscope with a fast camera allowed us to record Ca signals within an area of 465 μm by 315 μm with an acquisition speed of 55 fps. Images from multiple transients recorded at steady state were registered to their time point in the cardiac cycle to restore averaged local Ca transients with a higher temporal resolution. Local Ca transients within and between adjacent myocytes were compared with regard to amplitude, time to peak and decay at steady state stimulation (250 ms cycle length). Image registration from multiple sequential Ca transients allowed reconstruction of high temporal resolution (2.4 ± 1.3 ms) local CaT in 2D image sets (N = 4 hearts, n = 8 regions). During steady state stimulation, spatial Ca gradients were homogeneous within cells in both directions and independent of distance between measured points. Variation in CaT amplitudes was similar across the short and the long side of neighboring cells. Variations in TAU and TTP were similar in both directions. Isoproterenol enhanced the CaT but not the overall pattern of spatial heterogeneities. Here we detected and analyzed local Ca signals in intact mouse hearts with high temporal and spatial resolution, taking into account 2D arrangement of the cells. We observed significant differences in the variation of CaT amplitude along the long and short axis of cardiac myocytes. Variations of Ca signals between neighboring cells may contribute to the substrate of cardiac remodeling.

Micron-scale voltage and [Ca(2+)]i imaging in the intact heart

Studies in isolated cardiomyocytes have provided tremendous information at the cellular and molecular level concerning regulation of transmembrane voltage (V_m) and intracellular calcium ([Ca²⁺]_i). The ability to use the information gleaned to gain insight into the function of ion channels and Ca²⁺ handling proteins in a more complex system, e.g., the intact heart, has remained a challenge. We have developed laser scanning fluorescence microscopy-based approaches to monitor, at the sub-cellular to multi-cellular level in the immobilized, Langendorff-perfused mouse heart, dynamic changes in [Ca²⁺]_i and V_m. This article will review the use of single- or dual-photon laser scanning microscopy [Ca²⁺]_i imaging in conjunction with transgenic reporter technology to (a) interrogate the extent to which transplanted, donor-derived myocytes or cardiac stem cell-derived de novo myocytes are capable of forming a functional syncytium with the pre-existing myocardium, using entrainment of [Ca²⁺]_i transients by the electrical activity of the recipient heart as a surrogate for electrical coupling, and (b) characterize the Ca²⁺ handling phenotypes of cellular implants. Further, we will review the ability of laser scanning fluorescence microscopy in conjunction with a fast-response voltage-sensitive to resolve, on a subcellular level in Langendorff-perfused mouse hearts, V_m dynamics that typically occur during the course of a cardiac action potential. Specifically, the utility of this technique to measure microscopic-scale voltage gradients in the normal and diseased heart is discussed.

Alternating membrane potential/calcium interplay underlies repetitive focal activity in a genetic model of calcium-dependent atrial arrhythmias

KEY POINTS

Atrial fibrillation is often initiated and perpetuated by abnormal electrical pulses repetitively originating from regions outside the heart's natural pacemaker. In this study we examined the causal role of abnormal calcium releases from the sarcoplasmic reticulum in producing repetitive electrical discharges in atrial cells and tissues. Calsequestrin2 is a protein that stabilizes the closed state of calcium release channels, i.e. the ryanodine receptors. In the atria from mice predisposed to abnormal calcium releases secondary to the absence of calsequestrin2, we observed abnormal repetitive electrical discharges that may lead to atrial fibrillation. Here, we report a novel pathological rhythm generator. Specifically, abnormal calcium release leads to electrical activation, which in turn results in another abnormal calcium release. This process repeats itself and thus sustains the repetitive electrical discharges. These results suggest that improving the stability of ryanodine receptors might be useful to treat atrial fibrillation.

ABSTRACT

Aberrant diastolic calcium (Ca) release due to leaky ryanodine receptors (RyR2s) has been recently associated with atrial fibrillation (AF) and catecholaminergic polymorphic ventricular tachycardia (CPVT). However, it remains unclear how diastolic Ca release contributes to the rising of rapid repetitive focal activity, which is considered as a common AF triggering mechanism. To address this question, we conducted simultaneous voltage/Ca optical mapping in atrial tissue and one-/two-dimensional confocal imaging in atrial tissue and myocytes from wild-type (WT, n = 15) and CPVT mice lacking calsequestrin 2 (Casq2(-/-), n = 45), which promotes diastolic Ca release. During β-adrenergic stimulation (100 nM isoproterenol), only Casq2(-/-) atrial myocytes showed pacing-induced self-sustained repetitive activity (31 ± 21 s vs. none in WT). Importantly, in atrial tissue, this repetitive activity could translate to Ca-dependent focal arrhythmia. Ectopic action potential (AP) firing during repetitive activity occurred only when diastolic Ca release achieved a sufficient level of synchronization. The AP, in turn, synchronized subsequent diastolic Ca release by temporally aligning multiple sources of Ca waves both within individual myocytes and throughout the atrial tissue. This alternating interplay between AP and diastolic Ca release perpetuates the self-sustaining repetitive activity. In fact, pharmacological disruption of synchronized diastolic Ca release (by ryanodine) prevented aberrant APs; and vice versa, the inhibition of AP (by TTX or 0 Na, 0 Ca solution) de-synchronized diastolic Ca release. Taken together, these results suggest that a cyclical interaction between synchronized diastolic Ca release and AP forms a pathological rhythm generator that is involved in Ca-dependent atrial arrhythmias in CPVT.

Mechanisms of ventricular arrhythmias: from molecular fluctuations to electrical turbulence

Ventricular arrhythmias have complex causes and mechanisms. Despite extensive investigation involving many clinical, experimental, and computational studies, effective biological therapeutics are still very limited. In this article, we review our current understanding of the mechanisms of ventricular arrhythmias by summarizing the state of knowledge spanning from the molecular scale to electrical wave behavior at the tissue and organ scales and how the complex nonlinear interactions integrate into the dynamics of arrhythmias in the heart. We discuss the challenges that we face in synthesizing these dynamics to develop safe and effective novel therapeutic approaches.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

Cardiac sarcoplasmic reticulum calcium leak: basis and roles in cardiac dysfunction

Synchronized SR calcium (Ca) release is critical to normal cardiac myocyte excitation-contraction coupling, and ideally this release shuts off completely between heartbeats. However, other SR Ca release events are referred to collectively as SR Ca leak (which includes Ca sparks and waves as well as smaller events not detectable as Ca sparks). Much, but not all, of the SR Ca leak occurs via ryanodine receptors and can be exacerbated in pathological states such as heart failure. The extent of SR Ca leak is important because it can (a) reduce SR Ca available for release, causing systolic dysfunction; (b) elevate diastolic [Ca]i, contributing to diastolic dysfunction; (c) cause triggered arrhythmias; and (d) be energetically costly because of extra ATP used to repump Ca. This review addresses quantitative aspects and manifestations of SR Ca leak and its measurement, and how leak is modulated by Ca, associated proteins, and posttranslational modifications in health and disease.

The emergence of subcellular pacemaker sites for calcium waves and oscillations

Calcium (Ca(2+)) waves generating oscillatory Ca(2+) signals are widely observed in biological cells. Experimental studies have shown that under certain conditions, initiation of Ca(2+) waves is random in space and time, while under other conditions, waves occur repetitively from preferred locations (pacemaker sites) from which they entrain the whole cell. In this study, we use computer simulations to investigate the self-organization of Ca(2+) sparks into pacemaker sites generating Ca(2+) oscillations. In both ventricular myocyte experiments and computer simulations of a heterogeneous Ca(2+) release unit (CRU) network model, we show that Ca(2+) waves occur randomly in space and time when the Ca(2+) level is low, but as the Ca(2+) level increases, waves occur repetitively from the same sites. Our analysis indicates that this transition to entrainment can be attributed to the fact that random Ca(2+) sparks self-organize into Ca(2+) oscillations differently at low and high Ca(2+) levels. At low Ca(2+), the whole cell Ca(2+) oscillation frequency of the coupled CRU system is much slower than that of an isolated single CRU. Compared to a single CRU, the distribution of interspike intervals (ISIs) of the coupled CRU network exhibits a greater variation, and its ISI distribution is asymmetric with respect to the peak, exhibiting a fat tail. At high Ca(2+), however, the coupled CRU network has a faster frequency and lesser ISI variation compared to an individual CRU. The ISI distribution of the coupled network no longer exhibits a fat tail and is well-approximated by a Gaussian distribution. This same Ca(2+) oscillation behaviour can also be achieved by varying the number of ryanodine receptors per CRU or the distance between CRUs. Using these results, we develop a theory for the entrainment of random oscillators which provides a unified explanation for the experimental observations underlying the emergence of pacemaker sites and Ca(2+) oscillations.

Decreased RyR2 refractoriness determines myocardial synchronization of aberrant Ca2+ release in a genetic model of arrhythmia

Dysregulated intracellular Ca(2+) signaling is implicated in a variety of cardiac arrhythmias, including catecholaminergic polymorphic ventricular tachycardia. Spontaneous diastolic Ca(2+) release (DCR) can induce arrhythmogenic plasma membrane depolarizations, although the mechanism responsible for DCR synchronization among adjacent myocytes required for ectopic activity remains unclear. We investigated the synchronization mechanism(s) of DCR underlying untimely action potentials and diastolic contractions (DCs) in a catecholaminergic polymorphic ventricular tachycardia mouse model with a mutation in cardiac calsequestrin. We used a combination of different approaches including single ryanodine receptor channel recording, optical imaging (Ca(2+) and membrane potential), and contractile force measurements in ventricular myocytes and intact cardiac muscles. We demonstrate that DCR occurs in a temporally and spatially uniform manner in both myocytes and intact myocardial tissue isolated from cardiac calsequestrin mutation mice. Such synchronized DCR events give rise to triggered electrical activity that results in synchronous DCs in the myocardium. Importantly, we establish that synchronization of DCR is a result of a combination of abbreviated ryanodine receptor channel refractoriness and the preceding synchronous stimulated Ca(2+) release/reuptake dynamics. Our study reveals how aberrant DCR events can become synchronized in the intact myocardium, leading to triggered activity and the resultant DCs in the settings of a cardiac rhythm disorder.

The timing statistics of spontaneous calcium release in cardiac myocytes

A variety of cardiac arrhythmias are initiated by a focal excitation that disrupts the regular beating of the heart. In some cases it is known that these excitations are due to calcium (Ca) release from the sarcoplasmic reticulum (SR) via propagating subcellular Ca waves. However, it is not understood what are the physiological factors that determine the timing of these excitations at both the subcellular and tissue level. In this paper we apply analytic and numerical approaches to determine the timing statistics of spontaneous Ca release (SCR) in a simplified model of a cardiac myocyte. In particular, we compute the mean first passage time (MFPT) to SCR, in the case where SCR is initiated by spontaneous Ca sparks, and demonstrate that this quantity exhibits either an algebraic or exponential dependence on system parameters. Based on this analysis we identify the necessary requirements so that SCR occurs on a time scale comparable to the cardiac cycle. Finally, we study how SCR is synchronized across many cells in cardiac tissue, and identify a quantitative measure that determines the relative timing of SCR in an ensemble of cells. Using this approach we identify the physiological conditions so that cell-to-cell variations in the timing of SCR is small compared to the typical duration of an SCR event. We argue further that under these conditions inward currents due to SCR can summate and generate arrhythmogenic triggered excitations in cardiac tissue.

Ca²⁺ waves in the heart

Ca(2+) waves were probably first observed in the early 1940s. Since then Ca(2+) waves have captured the attention of an eclectic mixture of mathematicians, neuroscientists, muscle physiologists, developmental biologists, and clinical cardiologists. This review discusses the current state of mathematical models of Ca(2+) waves, the normal physiological functions Ca(2+) waves might serve in cardiac cells, as well as how the spatial arrangement of Ca(2+) release channels shape Ca(2+) waves, and we introduce the idea of Ca(2+) phase waves that might provide a useful framework for understanding triggered arrhythmias.

Store-dependent deactivation: cooling the chain-reaction of myocardial calcium signaling

In heart cells, Ca(2+) released from the internal storage unit, the sarcoplasmic reticulum (SR) through ryanodine receptor (RyR2) channels is the predominant determinant of cardiac contractility. Evidence obtained in recent years suggests that SR Ca(2+) release is tightly regulated not only by cytosolic Ca(2+) but also by intra-store Ca(2+) concentration. Specifically, Ca(2+)-induced Ca(2+) release (CICR) that relies on auto-catalytic action of Ca(2+) at the cytosolic side of RyR2s is precisely balanced and counteracted by RyR2 deactivation dependent on a reciprocal decrease of Ca(2+) at the luminal side of RyR2s. Dysregulation of this inherently unstable Ca(2+) signaling is considered to be an underlying cause of triggered arrhythmias, and is associated with genetic and acquired forms of sudden cardiac death. In this article, we present an overview of recent advances in our understanding of the regulatory role luminal Ca(2+) plays in Ca(2+) handling, with a particular emphasis on the role of Ca(2+)release refractoriness in aberrant Ca(2+) release.

Slow Calcium-Depolarization-Calcium waves may initiate fast local depolarization waves in ventricular tissue

Intercellular calcium waves in cardiac myocytes are a well-recognized, if incompletely understood, phenomenon. In a variety of preparations, investigators have reported multi-cellular calcium waves or triggered propagated contractions, but the mechanisms of propagation and pathological importance of these events remain unclear. Here, we review existing experimental data and present a computational approach to investigate the mechanisms of multi-cellular calcium wave propagation. Over the past 50 years, the standard modeling paradigm for excitable cardiac tissue has seen increasingly detailed models of the dynamics of individual cells coupled in tissue solely by intercellular and interstitial current flow. Although very successful, this modeling regime has been unable to capture two important phenomena: 1) the slow intercellular calcium waves observed experimentally, and 2) how intercellular calcium events resulting in delayed after depolarizations at the cellular level could overcome a source-sink mismatch to initiate depolarization waves in tissue. In this paper, we introduce a mathematical model with subcellular spatial resolution, in which we allow both inter- and intracellular current flow and calcium diffusion. In simulations of coupled cells employing this model, we observe: a) slow inter-cellular calcium waves propagating at about 0.1 mm/s, b) faster Calcium-Depolarization-Calcium (CDC) waves, traveling at about 1 mm/s, and c) CDC-waves that can set off fast depolarization-waves (50 cm/s) in tissue with varying gap-junction conductivity.

Criticality in intracellular calcium signaling in cardiac myocytes

Calcium (Ca) is a ubiquitous second messenger that regulates many biological functions. The elementary events of local Ca signaling are Ca sparks, which occur randomly in time and space, and integrate to produce global signaling events such as intra- and intercellular Ca waves and whole-cell Ca oscillations. Despite extensive experimental characterization in many systems, the transition from local random to global synchronous events is still poorly understood. Here we show that criticality, a ubiquitous dynamical phenomenon in nature, is responsible for the transition from local to global Ca signaling. We demonstrate this first in a computational model of Ca signaling in a cardiac myocyte and then experimentally in mouse ventricular myocytes, complemented by a theoretical agent-based model to delineate the underlying dynamics. We show that the interaction between the Ca release units via Ca-induced Ca release causes self-organization of Ca spark clusters. When the coupling between Ca release units is weak, the cluster-size distribution is exponential. As the interactions become strong, the cluster-size distribution changes to a power-law distribution, which is characteristic of criticality in thermodynamic and complex nonlinear systems, and facilitates the formation and propagation of Ca waves and whole-cell Ca oscillations. Our findings illustrate how criticality is harnessed by a biological cell to regulate Ca signaling via self-organization of random subcellular events into cellular-scale oscillations, and provide a general theoretical framework for the transition from local Ca signaling to global Ca signaling in biological cells.