Heart rate (HR) variability (HRV; beat-to-beat changes in the R-wave to R-wave interval) has attracted considerable attention during the past 30+ years (PubMed currently lists >17 000 publications). ...Clinically, a decrease in HRV is correlated to higher morbidity and mortality in diverse conditions, from heart disease to fetal distress. It is usually attributed to fluctuation in cardiac autonomic nerve activity. We calculated HRV parameters from a variety of cardiac preparations (including humans, living animals, Langendorff-perfused heart, and single sinoatrial nodal cell) in diverse species, combining this with data from previously published articles. We show that regardless of conditions, there is a universal exponential decay-like relationship between HRV and HR. Using 2 biophysical models, we develop a theory for this and confirm that HRV is primarily dependent on HR and cannot be used in any simple way to assess autonomic nerve activity to the heart. We suggest that the correlation between a change in HRV and altered morbidity and mortality is substantially attributable to the concurrent change in HR. This calls for re-evaluation of the findings from many articles that have not adjusted properly or at all for HR differences when comparing HRV in multiple circumstances.
Whether intracellular Ca2+ cycling dynamics regulate cardiac pacemaker cell function on a beat-to-beat basis remains unknown. Here we show that under physiological conditions, application of low ...concentrations of caffeine (2–4 mM) to isolated single rabbit sinoatrial node cells acutely reduces their spontaneous action potential cycle length (CL) and increases Ca2+ transient amplitude for several cycles. Numerical simulations, using a modified Maltsev-Lakatta coupled-clock model, faithfully reproduced these effects, and also the effects of CL prolongation and dysrhythmic spontaneous beating (produced by cytosolic Ca2+ buffering) and an acute CL reduction (produced by flash-induced Ca2+ release from a caged Ca2+ buffer), which we had reported previously. Three contemporary numerical models (including the original Maltsev-Lakatta model) failed to reproduce the experimental results. In our proposed new model, Ca2+ releases acutely change the CL via activation of the Na+/Ca2+ exchanger current. Time-dependent CL reductions after flash-induced Ca2+ releases (the memory effect) are linked to changes in Ca2+ available for pumping into sarcoplasmic reticulum which, in turn, changes the sarcoplasmic reticulum Ca2+ load, diastolic Ca2+ releases, and Na+/Ca2+ exchanger current. These results support the idea that Ca2+ regulates CL in cardiac pacemaker cells on a beat-to-beat basis, and suggest a more realistic numerical mechanism of this regulation.
Local Ca2+ Releases (LCRs) are crucial events involved in cardiac pacemaker cell function. However, specific algorithms for automatic LCR detection and analysis have not been developed in live, ...spontaneously beating pacemaker cells. In the present study we measured LCRs using a high-speed 2D-camera in spontaneously contracting sinoatrial (SA) node cells isolated from rabbit and guinea pig and developed a new algorithm capable of detecting and analyzing the LCRs spatially in two-dimensions, and in time. Our algorithm tracks points along the midline of the contracting cell. It uses these points as a coordinate system for affine transform, producing a transformed image series where the cell does not contract. Action potential-induced Ca2+ transients and LCRs were thereafter isolated from recording noise by applying a series of spatial filters. The LCR birth and death events were detected by a differential (frame-to-frame) sensitivity algorithm applied to each pixel (cell location). An LCR was detected when its signal changes sufficiently quickly within a sufficiently large area. The LCR is considered to have died when its amplitude decays substantially, or when it merges into the rising whole cell Ca2+ transient. Ultimately, our algorithm provides major LCR parameters such as period, signal mass, duration, and propagation path area. As the LCRs propagate within live cells, the algorithm identifies splitting and merging behaviors, indicating the importance of locally propagating Ca2+-induced-Ca2+-release for the fate of LCRs and for generating a powerful ensemble Ca2+ signal. Thus, our new computer algorithms eliminate motion artifacts and detect 2D local spatiotemporal events from recording noise and global signals. While the algorithms were developed to detect LCRs in sinoatrial nodal cells, they have the potential to be used in other applications in biophysics and cell physiology, for example, to detect Ca2+ wavelets (abortive waves), sparks and embers in muscle cells and Ca2+ puffs and syntillas in neurons.
Recent experimental studies have demonstrated that sinoatrial node cells (SANC) generate spontaneous, rhythmic, local subsarcolemmal Ca(2+) releases (Ca(2+) clock), which occur during late diastolic ...depolarization (DD) and interact with the classic sarcolemmal voltage oscillator (membrane clock) by activating Na(+)-Ca(2+) exchanger current (I(NCX)). This and other interactions between clocks, however, are not captured by existing essentially membrane-delimited cardiac pacemaker cell numerical models. Using wide-scale parametric analysis of classic formulations of membrane clock and Ca(2+) cycling, we have constructed and initially explored a prototype rabbit SANC model featuring both clocks. Our coupled oscillator system exhibits greater robustness and flexibility than membrane clock operating alone. Rhythmic spontaneous Ca(2+) releases of sarcoplasmic reticulum (SR)-based Ca(2+) clock ignite rhythmic action potentials via late DD I(NCX) over much broader ranges of membrane clock parameters e.g., L-type Ca(2+) current (I(CaL)) and/or hyperpolarization-activated ("funny") current (I(f)) conductances. The system Ca(2+) clock includes SR and sarcolemmal Ca(2+) fluxes, which optimize cell Ca(2+) balance to increase amplitudes of both SR Ca(2+) release and late DD I(NCX) as SR Ca(2+) pumping rate increases, resulting in a broad pacemaker rate modulation (1.8-4.6 Hz). In contrast, the rate modulation range via membrane clock parameters is substantially smaller when Ca(2+) clock is unchanged or lacking. When Ca(2+) clock is disabled, the system parametric space for fail-safe SANC operation considerably shrinks: without rhythmic late DD I(NCX) ignition signals membrane clock substantially slows, becomes dysrhythmic, or halts. In conclusion, the Ca(2+) clock is a new critical dimension in SANC function. A synergism of the coupled function of Ca(2+) and membrane clocks confers fail-safe SANC operation at greatly varying rates.
At the beginning of this century, debates regarding "what are the main control mechanisms that ignite the action potential (AP) in heart pacemaker cells" dominated the electrophysiology field. The ...original theory which prevailed for over 50 years had advocated that the ensemble of surface membrane ion channels (i.e., "M-clock") is sufficient to ignite rhythmic APs. However, more recent experimental evidence in a variety of mammals has shown that the sarcoplasmic reticulum (SR) acts as a "Ca(2+)-clock" rhythmically discharges diastolic local Ca(2+) releases (LCRs) beneath the cell surface membrane. LCRs activate an inward current (likely that of the Na(+)/Ca(2+) exchanger) that prompts the surface membrane "M-clock" to ignite an AP. Theoretical and experimental evidence has mounted to indicate that this clock "crosstalk" operates on a beat-to-beat basis and determines both the AP firing rate and rhythm. Our review is focused on the evolution of experimental definition and numerical modeling of the coupled-clock concept, on how mechanisms intrinsic to pacemaker cell determine both the heart rate and rhythm, and on future directions to develop further the coupled-clock pacemaker cell concept.
The current dogma about the heartbeat origin is based on "the pacemaker cell," a specialized cell residing in the sinoatrial node (SAN) that exhibits spontaneous diastolic depolarization triggering ...rhythmic action potentials (APs). Recent high-resolution imaging, however, demonstrated that Ca signals and APs in the SAN are heterogeneous, with many cells generating APs of different rates and rhythms or even remaining non-firing (dormant cells), i.e., generating only subthreshold signals. Here we numerically tested a hypothesis that a community of dormant cells can generate normal automaticity, i.e., "the pacemaker cell" is not required to initiate rhythmic cardiac impulses. Our model includes 1) non-excitable cells generating oscillatory local Ca releases and 2) an excitable cell lacking automaticity. While each cell in isolation was not "the pacemaker cell", the cell system generated rhythmic APs: The subthreshold signals of non-excitable cells were transformed into respective membrane potential oscillations
electrogenic Na/Ca exchange and further transferred and integrated (computed) by the excitable cells to reach its AP threshold, generating rhythmic pacemaking. Cardiac impulse is an emergent property of the SAN cellular network and can be initiated by cells lacking intrinsic automaticity. Cell heterogeneity, weak coupling, subthreshold signals, and their summation are critical properties of the new pacemaker mechanism, i.e., cardiac pacemaker can operate
a signaling process basically similar to that of "temporal summation" happening in a neuron with input from multiple presynaptic cells. The new mechanism, however, does not refute the classical pacemaker cell-based mechanism: both mechanisms can co-exist and interact within SAN tissue.
Recent data suggest that cardiac pacemaker cell function is determined by numerous time-, voltage-, and Ca-dependent interactions of cell membrane electrogenic proteins (M-clock) and intracellular Ca ...cycling proteins (Ca-clock), forming a coupled-clock system. Many aspects of the coupled-clock system, however, remain underexplored. The key players of the system are Ca release channels (ryanodine receptors), generating local Ca releases (LCRs) from sarcoplasmic reticulum, electrogenic Na/Ca exchanger (NCX) current, and L-type Ca current (ICaL). We combined numerical model simulations with experimental simultaneous recordings of action potentials (APs) and Ca to gain further insight into the complex interactions within the system. Our simulations revealed a positive feedback mechanism, dubbed AP ignition, which accelerates the diastolic depolarization (DD) to reach AP threshold. The ignition phase begins when LCRs begin to occur and the magnitude of inward NCX current begins to increase. The NCX current, together with funny current and T-type Ca current accelerates DD, bringing the membrane potential to ICaL activation threshold. During the ignition phase, ICaL-mediated Ca influx generates more LCRs via Ca-induced Ca release that further activates inward NCX current, creating a positive feedback. Simultaneous recordings of membrane potential and confocal Ca images support the model prediction of the positive feedback among LCRs and ICaL, as diastolic LCRs begin to occur below and continue within the voltage range of ICaL activation. The ignition phase onset (identified within the fine DD structure) begins when DD starts to notably accelerate (∼0.15 V/s) above the recording noise. Moreover, the timing of the ignition onset closely predicted the duration of each AP cycle in the basal state, in the presence of autonomic receptor stimulation, and in response to specific inhibition of either the M-clock or Ca-clock, thus indicating general importance of the new coupling mechanism for regulation of the pacemaker cell cycle duration, and ultimately the heart rate.