Ion channels on the surface membrane of sinoatrial nodal pacemaker cells (SANCs) are the proximal cause of an action potential. Each individual channel type has been thoroughly characterized under ...voltage clamp, and the ensemble of the ion channel currents reconstructed in silico generates rhythmic action potentials. Thus, this ensemble can be envisioned as a surface "membrane clock" (M clock). Localized subsarcolemmal Ca(2+) releases are generated by the sarcoplasmic reticulum via ryanodine receptors during late diastolic depolarization and are referred to as an intracellular "Ca(2+) clock," because their spontaneous occurrence is periodic during voltage clamp or in detergent-permeabilized SANCs, and in silico as well. In spontaneously firing SANCs, the M and Ca(2+) clocks do not operate in isolation but work together via numerous interactions modulated by membrane voltage, subsarcolemmal Ca(2+), and protein kinase A and CaMKII-dependent protein phosphorylation. Through these interactions, the 2 subsystem clocks become mutually entrained to form a robust, stable, coupled-clock system that drives normal cardiac pacemaker cell automaticity. G protein-coupled receptors signaling creates pacemaker flexibility, ie, effects changes in the rhythmic action potential firing rate, by impacting on these very same factors that regulate robust basal coupled-clock system function. This review examines evidence that forms the basis of this coupled-clock system concept in cardiac SANCs.
The present view on heartbeat initiation is that a primary pacemaker cell or a group of cells in the sinoatrial node (SAN) center paces the rest of the SAN and the atria. However, recent ...high-resolution imaging studies show a more complex paradigm of SAN function that emerges from heterogeneous signaling, mimicking brain cytoarchitecture and function. Here, we developed and tested a new conceptual numerical model of SAN organized similarly to brain networks featuring a modular structure with small-world topology. In our model, a lower rate module leads action potential (AP) firing in the basal state and during parasympathetic stimulation, whereas a higher rate module leads during β-adrenergic stimulation. Such a system reproduces the respective shift of the leading pacemaker site observed experimentally and a wide range of rate modulation and robust function while conserving energy. Since experimental studies found functional modules at different scales, from a few cells up to the highest scale of the superior and inferior SAN, the SAN appears to feature hierarchical modularity, i.e., within each module, there is a set of sub-modules, like in the brain, exhibiting greater robustness, adaptivity, and evolvability of network function. In this perspective, our model offers a new mainframe for interpreting new data on heterogeneous signaling in the SAN at different scales, providing new insights into cardiac pacemaker function and SAN-related cardiac arrhythmias in aging and disease.
Cardiac pacemaking is a complex phenomenon that is still not completely understood. Together with experimental studies, numerical modeling has been traditionally used to acquire mechanistic insights ...in this research area. This review summarizes the present state of numerical modeling of the cardiac pacemaker, including approaches to resolve present paradoxes and controversies. Specifically we discuss the requirement for realistic modeling to consider symmetrical importance of both intracellular and cell membrane processes (within a recent “coupled-clock” theory). Promising future developments of the complex pacemaker system models include the introduction of local calcium control, mitochondria function, and biochemical regulation of protein phosphorylation and cAMP production. Modern numerical and theoretical methods such as multi-parameter sensitivity analyses within extended populations of models and bifurcation analyses are also important for the definition of the most realistic parameters that describe a robust, yet simultaneously flexible operation of the coupled-clock pacemaker cell system. The systems approach to exploring cardiac pacemaker function will guide development of new therapies such as biological pacemakers for treating insufficient cardiac pacemaker function that becomes especially prevalent with advancing age.
Each heartbeat is initiated by specialized pacemaker cells operating within the sinoatrial node (SAN). While individual cells within SAN tissue exhibit substantial heterogeneity of their ...electrophysiological parameters and Ca cycling, the role of this heterogeneity for cardiac pacemaker function remains mainly unknown. Here we investigated the problem numerically in a 25 × 25 square grid of connected coupled-clock Maltsev-Lakatta cell models. The tissue models were populated by cells with different degree of heterogeneity of the two key model parameters, maximum L-type Ca current conductance (
) and sarcoplasmic reticulum Ca pumping rate (
). Our simulations showed that in the areas of
-
parametric space at the edge of the system stability, where action potential (AP) firing is absent or dysrhythmic in SAN tissue models populated with identical cells, rhythmic AP firing can be rescued by populating the tissues with heterogeneous cells. This robust SAN function is synergistic with respect to heterogeneity in
and
and can be further strengthened by clustering of cells with similar properties. The effect of cell heterogeneity is not due to a simple summation of activity of intrinsically firing cells naturally present in heterogeneous SAN; rather AP firing cells locally and critically interact with non-firing/dormant cells. When firing cells prevail, they recruit many dormant cells to fire, strongly enhancing overall SAN function; and vice versa, prevailing dormant cells suppress AP firing in cells with intrinsic automaticity and halt SAN function. The transitions between firing and non-firing states of the system are sharp, resembling phase transitions in statistical physics. Furthermore, robust function of heterogeneous SAN tissue requires weak cell coupling, a known property of the central area of SAN where cardiac impulse emerges; stronger cell coupling reduces AP firing rate and ultimately halts SAN automaticity at the edge of stability.
Intracellular Local Ca releases (LCRs) from sarcoplasmic reticulum (SR) regulate cardiac pacemaker cell function by activation of electrogenic Na/Ca exchanger (NCX) during diastole. Prior studies ...demonstrated the existence of powerful compensatory mechanisms of LCR regulation via a complex local cross-talk of Ca pump, release and NCX. One major obstacle to study these mechanisms is that LCR exhibit complex Ca release propagation patterns (including merges and separations) that have not been characterized. Here we developed new terminology, classification, and computer algorithms for automatic detection of numerically simulated LCRs and examined LCR regulation by SR Ca pumping rate (Pup) that provides a major contribution to fight-or-flight response. In our simulations the faster SR Ca pumping accelerates action potential-induced Ca transient decay and quickly clears Ca under the cell membrane in diastole, preventing premature releases. Then the SR generates an earlier, more synchronized, and stronger diastolic LCR signal activating an earlier and larger inward NCX current. LCRs at higher Pup exhibit larger amplitudes and faster propagation with more collisions to each other. The LCRs overlap with Ca transient decay, causing an elevation of the average diastolic Ca nadir to ~200 nM (at Pup = 24 mM/s). Background Ca (in locations lacking LCRs) quickly decays to resting Ca levels (<100 nM) at high Pup, but remained elevated during slower decay at low Pup. Release propagation is facilitated at higher Pup by a larger LCR amplitude, whereas at low Pup by higher background Ca. While at low Pup LCRs show smaller amplitudes, their larger durations and sizes combined with longer transient decay stabilize integrals of diastolic Ca and NCX current signals. Thus, the local interplay of SR Ca pump and release channels regulates LCRs and Ca transient decay to insure fail-safe pacemaker cell operation within a wide range of rates.
Excitation-contraction coupling kinetics is dictated by the action potential rate of sinoatrial-nodal cells. These cells generate local Ca releases (LCRs) that activate Na/Ca exchanger current, which ...accelerates diastolic depolarization and determines the pace. LCRs are generated by clusters of ryanodine receptors, Ca release units (CRUs), residing in the sarcoplasmic reticulum. While CRU distribution exhibits substantial heterogeneity, its functional importance remains unknown. Using numerical modeling, here we show that with a square lattice distribution of CRUs, Ca-induced-Ca-release propagation during diastolic depolarization is insufficient for pacemaking within a broad range of realistic ICaL densities. Allowing each CRU to deviate randomly from its lattice position allows sparks to propagate, as observed experimentally. As disorder increases, the CRU distribution exhibits larger empty spaces and simultaneously CRU clusters, as in Poisson clumping. Propagating within the clusters, Ca release becomes synchronized, increasing action potential rate and reviving pacemaker function of dormant/nonfiring cells. However, cells with fully disordered CRU positions could not reach low firing rates and their β-adrenergic-receptor stimulation effect was substantially decreased. Inclusion of Cav1.3, a low-voltage activation L-type Ca channel isoform into ICaL, strongly increases recruitment of CRUs to fire during diastolic depolarization, increasing robustness of pacemaking and complementing effects of CRU distribution. Thus, order/disorder in CRU locations along with Cav1.3 expression regulates pacemaker function via synchronization of CRU firing. Excessive CRU disorder and/or overexpression of Cav1.3 boosts pacemaker function in the basal state, but limits the rate range, which may contribute to heart rate range decline with age and disease.
Display omitted
•Cells isolated from SA node include dormant, dysrhythmic, and rhythmic cells.•All cells generate Local Ca Releases.•Local Ca releases can be increased and synchronized by ...β-adrenergic stimulation.•Dormant & dysrhythmic cells can fire rhythmically on β-adrenergic stimulation.
Current understanding of how cardiac pacemaker cells operate is based mainly on studies in isolated single sinoatrial node cells (SANC), specifically those that rhythmically fire action potentials similar to the in vivo behavior of the intact sinoatrial node. However, only a small fraction of SANC exhibit rhythmic firing after isolation. Other SANC behaviors have not been studied.
Here, for the first time, we studied all single cells isolated from the sinoatrial node of the guinea pig, including traditionally studied rhythmically firing cells (‘rhythmic SANC’), dysrhythmically firing cells (‘dysrhythmic SANC’) and cells without any apparent spontaneous firing activity (‘dormant SANC’). Action potential-induced cytosolic Ca2+ transients and spontaneous local Ca2+ releases (LCRs) were measured with a 2D camera.
LCRs were present not only in rhythmically firing SANC, but also in dormant and dysrhythmic SANC. While rhythmic SANC were characterized by large LCRs synchronized in space and time towards late diastole, dysrhythmic and dormant SANC exhibited smaller LCRs that appeared stochastically and were widely distributed in time. β-adrenergic receptor (βAR) stimulation increased LCR size and synchronized LCR occurrences in all dysrhythmic and a third of dormant cells (25 of 75 cells tested). In response to βAR stimulation, these dormant SANC developed automaticity, and LCRs became coupled to spontaneous action potential-induced cytosolic Ca2+ transients. Conversely, dormant SANC that did not develop automaticity showed no significant change in average LCR characteristics. The majority of dysrhythmic cells became rhythmic in response to βAR stimulation, with the rate of action potential-induced cytosolic Ca2+ transients substantially increasing. In summary, isolated SANC can be broadly categorized into three major populations: dormant, dysrhythmic, and rhythmic. We interpret our results based on simulations of a numerical model of SANC operating as a coupled-clock system. On this basis, the two previously unstudied dysrhythmic and dormant cell populations have intrinsically partially or completely uncoupled clocks. Such cells can be recruited to fire rhythmically in response to βAR stimulation via increased rhythmic LCR activity and ameliorated coupling between the Ca2+ and membrane clocks.
Heart muscle contraction is normally activated by a synchronized Ca release from sarcoplasmic reticulum (SR), a major intracellular Ca store. However, under abnormal conditions, Ca leaks from the SR, ...decreasing heart contraction amplitude and increasing risk of life-threatening arrhythmia. The mechanisms and regimes of SR operation generating the abnormal Ca leak remain unclear. Here, we employed both numerical and analytical modeling to get mechanistic insights into the emergent Ca leak phenomenon. Our numerical simulations using a detailed realistic model of the Ca release unit reveal sharp transitions resulting in Ca leak. The emergence of leak is closely mapped mathematically to the Ising model from statistical mechanics. The system steady-state behavior is determined by two aggregate parameters: the analogs of magnetic field (h) and the inverse temperature (β) in the Ising model, for which we have explicit formulas in terms of SR Ca and release channel opening and closing rates. The classification of leak regimes takes the shape of a phase β-h diagram, with the regime boundaries occurring at h = 0 and a critical value of β (β∗) that we estimate using a classical Ising model and mean field theory. Our theory predicts that a synchronized Ca leak will occur when h > 0 and β >β∗, and a disordered leak occurs when β <β∗ and h is not too negative. The disorder leak is distinguished from synchronized leak (in long-lasting sparks) by larger Peierls contour lengths, an output parameter reflecting degree of disorder. Thus, in addition to our detailed numerical model approach, we also offer an instantaneous computational tool using analytical formulas of the Ising model for respective ryanodine receptor parameters and SR Ca load that describe and classify phase transitions and leak emergence.