Reciprocal role of the inward currents ib, Na and i(f) in controlling and stabilizing pacemaker frequency of rabbit sino-atrial node cells

Genetic influences on motor learning and performance and superperforming mutants revealed by random mutational survey of the mouse genome

Evolution depends upon genetic variations that influence physiology. As defined in a genetic screen, phenotypic performance may be enhanced or degraded by such mutations. We set out to detect mutations that influence motor function, including motor learning. Thus, we tested the motor effects of 36,444 non-synonymous coding/splicing mutations induced in the germline of C57BL/6J mice with N-ethyl-N-nitrosourea by measuring changes in the performance of repetitive rotarod trials while blinded to genotype. Automated meiotic mapping was used to implicate individual mutations in causation. In total, 32,726 mice bearing all the variant alleles were screened. This was complemented with the simultaneous testing of 1408 normal mice for reference. In total, 16.3% of autosomal genes were thus rendered detectably hypomorphic or nullified by mutations in homozygosity and motor tested in at least three mice. This approach allowed us to identify superperformance mutations in Rif1, Tk1, Fan1 and Mn1. These genes are primarily related, among other less well-characterized functions, to nucleic acid biology. We also associated distinct motor learning patterns with groups of functionally related genes. These functional sets included, preferentially, histone H3 methyltransferase activity for mice that learnt at an accelerated rate relative to the remaining mutant mice. The results allow for an estimation of the fraction of mutations that can modify a behaviour influential for evolution such as locomotion. They may also enable, once the loci are further validated and the mechanisms elucidated, the harnessing of the activity of the newly identified genes to enhance motor ability or to counterbalance disability or disease. KEY POINTS: We studied the effect of chemically induced random mutations on mouse motor performance. An array of mutations influenced the rate of motor learning. DNA regulation genes predominated among these mutant loci. Several mutations in unsuspected genes led to superperformance. Assuming little-biased mutagenicity, the results allow for an estimation of the probability for any spontaneous mutation to influence a behaviour such as motor learning and ultimate performance.

Impoverished Conceptions of Gene Causation and Therapy in Developmental Neurology

We offer a primer to the modifiability of genetic neurological disease, particularly during development. One goal is to harness several unexpected observations made in the course of experimental gene modification or therapy into an explanatory conceptual context based on biological first principles. To this end, we anchor growing, disparate reports of unusual or untoward effects to a plausible framework wherein genes exhibit different degrees of modifiability and may result, when mutated or therapeutically modified, in unsuspected consequences. We propose that genetic pathogenic variant effects and modifiability depend on the number and complexity of associated protein-protein or higher-order interactions. Thus, gene malleability may range from that characteristic of the favorably modifiable primarily structural genes that subserve relatively invariant or circumscribed phenomena such as cell shape to that typical of some transcription factors, which are less functionally predictable when altered. The latter may be expressed developmentally, in compartmentalized manner, or only intermittently and yet exert vastly ramified influences sometimes circumscribed only to select species. We also argue that genetic diseases may steer the organism toward often poorly understood biological end points and co-opt multiple processes into hardly modifiable biology. Addition or modification of genes to approximate a normal state not previously experienced by the organism may lead to further aberration due to extraneous interference with the native biology of the disease state. Therefore, an understanding as perspicuous as possible of gene function, regulation, modifiability, and biological directionality down to seemingly minute but disease-relevant consequences is a prerequisite to intervention. Although we provide some groundwork steps to such an understanding, this may occasionally prove unattainable.

Genetic influences on motor learning and superperformance mutants revealed by random mutational survey of mouse locomotion

Evolution depends upon genetic variations that influence physiology. As defined in a genetic screen, phenotypic performance may be enhanced or degraded by such mutations. We set out to detect mutations that influence motor function, including motor learning. Thus, we tested the motor effects of 36,444 non-synonymous coding/splicing mutations induced in the germline of C57BL/6J mice with N-ethyl-N-nitrosourea by measuring changes in the performance of repetitive rotarod trials while blinded to genotype. Automated meiotic mapping was used to implicate individual mutations in causation. 32,726 mice bearing all the variant alleles were screened. This was complemented with the simultaneous testing of 1,408 normal mice for reference. 16.3% of autosomal genes were thus rendered detectably hypomorphic or nullified by mutations in homozygosity and motor tested in at least 3 mice. This approach allowed us to identify superperformance mutations in Rif1, Tk1, Fan1 and Mn1. These genes are primarily related, among other less well characterized functions, to nucleic acid biology. We also associated distinct motor learning patterns with groups of functionally related genes. These functional sets included preferentially histone H3 methyltransferase activity for mice that learnt at an accelerated rate relative to the rest of mutant mice. The results allow for an estimation of the fraction of mutations that can modify a behavior influential for evolution such as locomotion. They may also enable, once the loci are further validated and the mechanisms elucidated, the harnessing of the activity of the newly identified genes to enhance motor ability or to counterbalance disability or disease.

Constrained disorder principle-based variability is fundamental for biological processes: Beyond biological relativity and physiological regulatory networks

The constrained disorder principle (CDP) defines systems based on their degree of disorder bounded by dynamic boundaries. The principle explains stochasticity in living and non-living systems. Denis Noble described the importance of stochasticity in biology, emphasizing stochastic processes at molecular, cellular, and higher levels in organisms as having a role beyond simple noise. The CDP and Noble's theories (NT) claim that biological systems use stochasticity. This paper presents the CDP and NT, discussing common notions and differences between the two theories. The paper presents the CDP-based concept of taking the disorder beyond its role in nature to correct malfunctions of systems and improve the efficiency of biological systems. The use of CDP-based algorithms embedded in second-generation artificial intelligence platforms is described. In summary, noise is inherent to complex systems and has a functional role. The CDP provides the option of using noise to improve functionality.

Sodium background currents in endocrine/neuroendocrine cells: Towards unraveling channel identity and contribution in hormone secretion

In endocrine/neuroendocrine tissues, excitability of secretory cells is patterned by the repertoire of ion channels and there is clear evidence that extracellular sodium (Na⁺) ions contribute to hormone secretion. While voltage-gated channels involved in action potential generation are well-described, the background 'leak' channels operating near the resting membrane potential are much less known, and in particular the channels supporting a background entry of Na⁺ ions. These background Na⁺ currents (called here 'I_Nab') have the ability to modulate the resting membrane potential and subsequently affect action potential firing. Here we compile and analyze the data collected from three endocrine/neuroendocrine tissues: the anterior pituitary gland, the adrenal medulla and the endocrine pancreas. We also model how I_Nab can be functionally involved in cellular excitability. Finally, towards deciphering the physiological role of I_Nab in endocrine/neuroendocrine cells, its implication in hormone release is also discussed.

The funny current in genetically modified mice

Since its first description in 1979, the hyperpolarization-activated funny current (I_f) has been the object of intensive research aimed at understanding its role in cardiac pacemaker activity and its modulation by the sympathetic and parasympathetic branches of the autonomic nervous system. I_f was described in isolated tissue strips of the rabbit sinoatrial node using the double-electrode voltage-clamp technique. Since then, the rabbit has been the principal animal model for studying pacemaker activity and I_f for more than 20 years. In 2001, the first study describing the electrophysiological properties of mouse sinoatrial pacemaker myocytes and those of I_f was published. It was soon followed by the description of murine myocytes of the atrioventricular node and the Purkinje fibres. The sinoatrial node of genetically modified mice has become a very popular model for studying the mechanisms of cardiac pacemaker activity. This field of research benefits from the impressive advancement of in-vivo exploration techniques of physiological parameters, imaging, genetics, and large-scale genomic approaches. The present review discusses the influence of mouse genetic on the most recent knowledge of the funny current's role in the physiology and pathophysiology of cardiac pacemaker activity. Genetically modified mice have provided important insights into the role of I_f in determining intrinsic automaticity in vivo and in myocytes of the conduction system. In addition, gene targeting of f-(HCN) channel isoforms have contributed to elucidating the current's role in the regulation of heart rate by the parasympathetic nervous system. This review is dedicated to Dario DiFrancesco on his retirement.

A sodium background conductance controls the spiking pattern of mouse adrenal chromaffin cells in situ

KEY POINTS

Mouse chromaffin cells in acute adrenal slices exhibit two distinct spiking patterns, a repetitive mode and a bursting mode. A sodium background conductance operates at rest as demonstrated by the membrane hyperpolarization evoked by a low Na⁺ -containing extracellular saline. This sodium background current is insensitive to TTX, is not blocked by Cs⁺ ions and displays a linear I-V relationship at potentials close to chromaffin cell resting potential. Its properties are reminiscent of those of the sodium leak channel NALCN. In the adrenal gland, Nalcn mRNA is selectively expressed in chromaffin cells. The study fosters our understanding of how the spiking pattern of chromaffin cells is regulated and adds a sodium background conductance to the list of players involved in the stimulus-secretion coupling of the adrenomedullary tissue.

ABSTRACT

Chromaffin cells (CCs) are the master neuroendocrine units for the secretory function of the adrenal medulla and a finely-tuned regulation of their electrical activity is required for appropriate catecholamine secretion in response to the organismal demand. Here, we aim at deciphering how the spiking pattern of mouse CCs is regulated by the ion conductances operating near the resting membrane potential (RMP). At RMP, mouse CCs display a composite firing pattern, alternating between active periods composed of action potentials spiking with a regular or a bursting mode, and silent periods. RMP is sensitive to changes in extracellular sodium concentration, and a low Na⁺ -containing saline hyperpolarizes the membrane, regardless of the discharge pattern. This RMP drive reflects the contribution of a depolarizing conductance, which is (i) not blocked by tetrodotoxin or caesium, (ii) displays a linear I-V relationship between -110 and -40 mV, and (iii) is carried by cations with a conductance sequence g_Na  > g_K  > g_Cs . These biophysical attributes, together with the expression of the sodium-leak channel Nalcn transcript in CCs, state credible the contribution of NALCN. This inaugural report opens new research routes in the field of CC stimulus-secretion coupling, and extends the inventory of tissues in which NALCN is expressed to neuroendocrine glands.

Sinoatrial Node Structure, Mechanics, Electrophysiology and the Chronotropic Response to Stretch in Rabbit and Mouse

The rhythmic electrical activity of the heart's natural pacemaker, the sinoatrial node (SAN), determines cardiac beating rate (BR). SAN electrical activity is tightly controlled by multiple factors, including tissue stretch, which may contribute to adaptation of BR to changes in venous return. In most animals, including human, there is a robust increase in BR when the SAN is stretched. However, the chronotropic response to sustained stretch differs in mouse SAN, where it causes variable responses, including decreased BR. The reasons for this species difference are unclear. They are thought to relate to dissimilarities in SAN electrophysiology (particularly action potential morphology) between mouse and other species and to how these interact with subcellular stretch-activated mechanisms. Furthermore, species-related differences in structural and mechanical properties of the SAN may influence the chronotropic response to SAN stretch. Here we assess (i) how the BR response to sustained stretch of rabbit and mouse isolated SAN relates to tissue stiffness, (ii) whether structural differences could account for observed differences in BR responsiveness to stretch, and (iii) whether pharmacological modification of mouse SAN electrophysiology alters stretch-induced chronotropy. We found disparities in the relationship between SAN stiffness and the magnitude of the chronotropic response to stretch between rabbit and mouse along with differences in SAN collagen structure, alignment, and changes with stretch. We further observed that pharmacological modification to prolong mouse SAN action potential plateau duration rectified the direction of BR changes during sustained stretch, resulting in a positive chronotropic response akin to that of other species. Overall, our results suggest that structural, mechanical, and background electrophysiological properties of the SAN influence the chronotropic response to stretch. Improved insight into the biophysical determinants of stretch effects on SAN pacemaking is essential for a comprehensive understanding of SAN regulation with important implications for studies of SAN physiology and its dysfunction, such as in the aging and fibrotic heart.

Considering discrepancy when calibrating a mechanistic electrophysiology model

Uncertainty quantification (UQ) is a vital step in using mathematical models and simulations to take decisions. The field of cardiac simulation has begun to explore and adopt UQ methods to characterize uncertainty in model inputs and how that propagates through to outputs or predictions; examples of this can be seen in the papers of this issue. In this review and perspective piece, we draw attention to an important and under-addressed source of uncertainty in our predictions-that of uncertainty in the model structure or the equations themselves. The difference between imperfect models and reality is termed model discrepancy, and we are often uncertain as to the size and consequences of this discrepancy. Here, we provide two examples of the consequences of discrepancy when calibrating models at the ion channel and action potential scales. Furthermore, we attempt to account for this discrepancy when calibrating and validating an ion channel model using different methods, based on modelling the discrepancy using Gaussian processes and autoregressive-moving-average models, then highlight the advantages and shortcomings of each approach. Finally, suggestions and lines of enquiry for future work are provided. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Cardiac Mechano-Electric Coupling: Acute Effects of Mechanical Stimulation on Heart Rate and Rhythm

The heart is vital for biological function in almost all chordates, including humans. It beats continually throughout our life, supplying the body with oxygen and nutrients while removing waste products. If it stops, so does life. The heartbeat involves precise coordination of the activity of billions of individual cells, as well as their swift and well-coordinated adaption to changes in physiological demand. Much of the vital control of cardiac function occurs at the level of individual cardiac muscle cells, including acute beat-by-beat feedback from the local mechanical environment to electrical activity (as opposed to longer term changes in gene expression and functional or structural remodeling). This process is known as mechano-electric coupling (MEC). In the current review, we present evidence for, and implications of, MEC in health and disease in human; summarize our understanding of MEC effects gained from whole animal, organ, tissue, and cell studies; identify potential molecular mediators of MEC responses; and demonstrate the power of computational modeling in developing a more comprehensive understanding of ‟what makes the heart tick.ˮ.

Patterns of cardiovascular variability after long-term sino-aortic denervation in unanesthetized adult rats

Baroreflex dysfunction is a diffuse chronic condition that is expected to be followed by a profound loss of organization of BP and HR variability. Nevertheless, long-term effects of baroreflex withdrawal are still debated. Aim of our work was to study BP and HR changes long term after sino-aortic denervation (SAD). Inter-beat-interval (IBI) and intra-arterial BP were recorded beat-by-beat in 43 Wistar-Kyoto rats (Controls, n = 33; SAD rats, n = 10). Power spectra were calculated in controls and in SAD rats within three days and at seven months from denervation. Compared to controls, chronic SAD rats showed 1) similar mean BP (control vs SAD: 95 ± 16 vs 87 ± 22 mmHg) and IBI (171 ± 22 vs 181 ± 15 ms) values, 2) dramatically higher values of BP variance (12 ± 2 vs 64 ± 2 mmHg2, p < 0.01) and of ultra- (ULF) and very-low-frequency (VLF) BP oscillations, 3) dramatically higher values of IBI variability (24 ± 2 vs 71 ± 4 ms2, p < 0.01) and of ULF-IBI oscillations that were synchronized with BP oscillations. Chronic SAD rats reveal a marked change in the pattern of cardiovascular variability characterized by the appearance of synchronized slower oscillations of BP and HR. The cardiovascular system, therefore, retains a high level of organization despite the absence of a reflex control mechanism.

Electrophysiological heterogeneity of pacemaker cells in the rabbit intercaval region, including the SA node: insights from recording multiple ion currents in each cell

Cardiac pacemaker cells, including cells of the sinoatrial node, are heterogeneous in size, morphology, and electrophysiological characteristics. The exact extent to which these cells differ electrophysiologically is unclear yet is critical to understanding their functioning. We examined major ionic currents in individual intercaval pacemaker cells (IPCs) sampled from the paracristal, intercaval region (including the sinoatrial node) that were spontaneously beating after enzymatic isolation from rabbit hearts. The beating rate was measured at baseline and after inhibition of the Ca2+ pump with cyclopiazonic acid. Thereafter, in each cell, we consecutively measured the density of funny current ( If), delayed rectifier K+ current ( IK) (a surrogate of repolarization capacity), and L-type Ca2+ current ( ICa,L) using whole cell patch clamp. The ionic current densities varied to a greater extent than previously appreciated, with some IPCs demonstrating very small or zero If . The density of none of the currents was correlated with cell size, while ICa,L and If densities were related to baseline beating rates. If density was correlated with IK density but not with that of ICa,L. Inhibition of Ca2+ cycling had a greater beating rate slowing effect in IPCs with lower If densities. Our numerical model simulation indicated that 1) IPCs with small (or zero) If or small ICa,L can operate via a major contribution of Ca2+ clock, 2) If-Ca2+-clock interplay could be important for robust pacemaking function, and 3) coupled If- IK function could regulate maximum diastolic potential. Thus, we have demonstrated marked electrophysiological heterogeneity of IPCs. This heterogeneity is manifested in basal beating rate and response to interference of Ca2+ cycling, which is linked to If. NEW & NOTEWORTHY In the present study, a hitherto unrecognized range of heterogeneity of ion currents in pacemaker cells from the intercaval region is demonstrated. Relationships between basal beating rate and L-type Ca2+ current and funny current ( If) density are uncovered, along with a positive relationship between If and delayed rectifier K+ current. Links are shown between the response to Ca2+ cycling blockade and If density.

Evolution viewed from physics, physiology and medicine

Stochasticity is harnessed by organisms to generate functionality. Randomness does not, therefore, necessarily imply lack of function or 'blind chance' at higher levels. In this respect, biology must resemble physics in generating order from disorder. This fact is contrary to Schrödinger's idea of biology generating phenotypic order from molecular-level order, which inspired the central dogma of molecular biology. The order originates at higher levels, which constrain the components at lower levels. We now know that this includes the genome, which is controlled by patterns of transcription factors and various epigenetic and reorganization mechanisms. These processes can occur in response to environmental stress, so that the genome becomes 'a highly sensitive organ of the cell' (McClintock). Organisms have evolved to be able to cope with many variations at the molecular level. Organisms also make use of physical processes in evolution and development when it is possible to arrive at functional development without the necessity to store all information in DNA sequences. This view of development and evolution differs radically from that of neo-Darwinism with its emphasis on blind chance as the origin of variation. Blind chance is necessary, but the origin of functional variation is not at the molecular level. These observations derive from and reinforce the principle of biological relativity, which holds that there is no privileged level of causation. They also have important implications for medical science.

Uncertainty and variability in computational and mathematical models of cardiac physiology

KEY POINTS

Mathematical and computational models of cardiac physiology have been an integral component of cardiac electrophysiology since its inception, and are collectively known as the Cardiac Physiome. We identify and classify the numerous sources of variability and uncertainty in model formulation, parameters and other inputs that arise from both natural variation in experimental data and lack of knowledge. The impact of uncertainty on the outputs of Cardiac Physiome models is not well understood, and this limits their utility as clinical tools. We argue that incorporating variability and uncertainty should be a high priority for the future of the Cardiac Physiome. We suggest investigating the adoption of approaches developed in other areas of science and engineering while recognising unique challenges for the Cardiac Physiome; it is likely that novel methods will be necessary that require engagement with the mathematics and statistics community.

ABSTRACT

The Cardiac Physiome effort is one of the most mature and successful applications of mathematical and computational modelling for describing and advancing the understanding of physiology. After five decades of development, physiological cardiac models are poised to realise the promise of translational research via clinical applications such as drug development and patient-specific approaches as well as ablation, cardiac resynchronisation and contractility modulation therapies. For models to be included as a vital component of the decision process in safety-critical applications, rigorous assessment of model credibility will be required. This White Paper describes one aspect of this process by identifying and classifying sources of variability and uncertainty in models as well as their implications for the application and development of cardiac models. We stress the need to understand and quantify the sources of variability and uncertainty in model inputs, and the impact of model structure and complexity and their consequences for predictive model outputs. We propose that the future of the Cardiac Physiome should include a probabilistic approach to quantify the relationship of variability and uncertainty of model inputs and outputs.

Current concepts of anatomy and electrophysiology of the sinus node

The sinoatrial node, or sinus node, of humans is the principal pacemaker of the heart. Over the last century, studies have unraveled the complex molecular architecture of the sinus node and the expression of unique ion channels within its specialized myocytes. Aim of this review is to describe the embriology, the anatomy, the histology and the electrophisiology of the sinus node.

Characterization and influence of cardiac background sodium current in the atrioventricular node

Background inward sodium current (IB,Na) that influences cardiac pacemaking has been comparatively under-investigated. The aim of this study was to determine for the first time the properties and role of IB,Na in cells from the heart's secondary pacemaker, the atrioventricular node (AVN). Myocytes were isolated from the AVN of adult male rabbits and mice using mechanical and enzymatic dispersion. Background current was measured using whole-cell patch clamp and monovalent ion substitution with major voltage- and time-dependent conductances inhibited. In the absence of a selective pharmacological inhibitor of IB,Na, computer modelling was used to assess the physiological contribution of IB,Na. Net background current during voltage ramps was linear, reversing close to 0mV. Switching between Tris- and Na(+)-containing extracellular solution in rabbit and mouse AVN cells revealed an inward IB,Na, with an increase in slope conductance in rabbit cells at -50mV from 0.54±0.03 to 0.91±0.05nS (mean±SEM; n=61 cells). IB,Na magnitude varied in proportion to [Na(+)]o. Other monovalent cations could substitute for Na(+) (Rb(+)>K(+)>Cs(+)>Na(+)>Li(+)). The single-channel conductance with Na(+) as charge carrier estimated from noise-analysis was 3.2±1.2pS (n=6). Ni(2+) (10mM), Gd(3+) (100μM), ruthenium red (100μM), or amiloride (1mM) produced modest reductions in IB,Na. Flufenamic acid was without significant effect, whilst La(3+) (100μM) or extracellular acidosis (pH6.3) inhibited the current by >60%. Under the conditions of our AVN cell simulations, removal of IB,Na arrested spontaneous activity and, in a simulated 1D-strand, reduced conduction velocity by ~20%. IB,Na is carried by distinct low conductance monovalent non-selective cation channels and can influence AVN spontaneous activity and conduction.

The importance of Ca(2+)-dependent mechanisms for the initiation of the heartbeat

Mechanisms underlying pacemaker activity in the sinus node remain controversial, with some ascribing a dominant role to timing events in the surface membrane (“membrane clock”) and others to uptake and release of calcium from the sarcoplasmic reticulum (SR) (“calcium clock”). Here we discuss recent evidence on mechanisms underlying pacemaker activity with a particular emphasis on the many roles of calcium. There are particular areas of controversy concerning the contribution of calcium spark-like events and the importance of I(f) to spontaneous diastolic depolarisation, though it will be suggested that neither of these is essential for pacemaking. Sodium-calcium exchange (NCX) is most often considered in the context of mediating membrane depolarisation after spark-like events. We present evidence for a broader role of this electrogenic exchanger which need not always depend upon these spark-like events. Short (milliseconds or seconds) and long (minutes) term influences of calcium are discussed including direct regulation of ion channels and NCX, and control of the activity of calcium-dependent enzymes (including CaMKII, AC1, and AC8). The balance between the many contributory factors to pacemaker activity may well alter with experimental and clinical conditions, and potentially redundant mechanisms are desirable to ensure the regular spontaneous heart rate that is essential for life. This review presents evidence that calcium is central to the normal control of pacemaking across a range of temporal scales and seeks to broaden the accepted description of the “calcium clock” to cover these important influences.

Beat-to-beat cycle length variability of spontaneously beating guinea pig sinoatrial cells: relative contributions of the membrane and calcium clocks

The heartbeat arises rhythmically in the sino-atrial node (SAN) and then spreads regularly throughout the heart. The molecular mechanism underlying SAN rhythm has been attributed by recent studies to the interplay between two clocks, one involving the hyperpolarization activated cation current I_f (the membrane clock), and the second attributable to activation of the electrogenic NaCa exchanger by spontaneous sarcoplasmic releases of calcium (the calcium clock). Both mechanisms contain, in principle, sources of beat-to-beat cycle length variability, which can determine the intrinsic variability of SAN firing and, in turn, contribute to the heart rate variability. In this work we have recorded long sequences of action potentials from patch clamped guinea pig SAN cells (SANCs) perfused, in turn, with normal Tyrode solution, with the I_f inhibitor ivabradine (3 µM), then back to normal Tyrode, and again with the ryanodine channels inhibitor ryanodine (3 µM). We have found that, together with the expected increase in beating cycle length (+25%), the application of ivabradine brought about a significant and dramatic increase in beat-to-beat cycle length variability (+50%). Despite the similar effect on firing rate, ryanodine did not modify significantly beat-to-beat cycle length variability. Acetylcholine was also applied and led to a 131% increase of beating cycle length, with only a 70% increase in beat-to-beat cycle length variability. We conclude that the main source of inter-beat variability of SANCs firing rate is related to the mechanism of the calcium clock, whereas the membrane clock seems to act in stabilizing rate. Accordingly, when the membrane clock is silenced by application of ivabradine, stochastic variations of the calcium clock are free to make SANCs beating rhythm more variable.

Numerical models based on a minimal set of sarcolemmal electrogenic proteins and an intracellular Ca(2+) clock generate robust, flexible, and energy-efficient cardiac pacemaking

Recent evidence supports the idea that robust and, importantly, FLEXIBLE automaticity of cardiac pacemaker cells is conferred by a coupled system of membrane ion currents (an "M-clock") and a sarcoplasmic reticulum (SR)-based Ca(2+) oscillator ("Ca(2+)clock") that generates spontaneous diastolic Ca(2+) releases. This study identified numerical models of a human biological pacemaker that features robust and flexible automaticity generated by a minimal set of electrogenic proteins and a Ca(2+)clock. Following the Occam's razor principle (principle of parsimony), M-clock components of unknown molecular origin were excluded from Maltsev-Lakatta pacemaker cell model and thirteen different model types of only 4 or 5 components were derived and explored by a parametric sensitivity analysis. The extended ranges of SR Ca(2+) pumping (i.e. Ca(2+)clock performance) and conductance of ion currents were sampled, yielding a large variety of parameter combination, i.e. specific model sets. We tested each set's ability to simulate autonomic modulation of human heart rate (minimum rate of 50 to 70bpm; maximum rate of 140 to 210bpm) in response to stimulation of cholinergic and β-adrenergic receptors. We found that only those models that include a Ca(2+)clock (including the minimal 4-parameter model "ICaL+IKr+INCX+Ca(2+)clock") were able to reproduce the full range of autonomic modulation. Inclusion of If or ICaT decreased the flexibility, but increased the robustness of the models (a relatively larger number of sets did not fail during testing). The new models comprised of components with clear molecular identity (i.e. lacking IbNa & Ist) portray a more realistic pacemaking: A smaller Na(+) influx is expected to demand less energy for Na(+) extrusion. The new large database of the reduced coupled-clock numerical models may serve as a useful tool for the design of biological pacemakers. It will also provide a conceptual basis for a general theory of robust, flexible, and energy-efficient pacemaking based on realistic components.

Effect of hyperpolarization-activated current I(f) on robustness of sinoatrial node pacemaking: theoretical study on influence of intracellular Na(+) concentration

To elucidate the effects of hyperpolarization-activated current I(f) on robustness of sinoatrial node (SAN) pacemaking in connection with intracellular Na(+) concentration (Na(i)) changes, we theoretically investigated 1) the impacts of I(f) on dynamical properties of SAN model cells during inhibition of L-type Ca(2+) channel currents (I(CaL)) or hyperpolarizing loads and 2) I(f)-dependent changes in Na(i) and their effects on dynamical properties of model cells. Bifurcation analyses were performed for Na(i)-variable and Na(i)-fixed versions of mathematical models for rabbit SAN cells; equilibrium points (EPs), limit cycles (LCs), and their stability were determined as functions of model parameters. Increasing I(f) conductance (g(f)) shrank I(CaL) conductance (g(CaL)) regions of unstable EPs and stable LCs (rhythmic firings) in the Na(i)-variable system but slightly broadened that of rhythmic firings at lower g(f) in the Na(i)-fixed system. In the Na(i)-variable system, increased g(f) yielded elevations in Na(i) at EPs and during spontaneous oscillations, which caused EP stabilization and shrinkage in the parameter regions of unstable EPs and rhythmic firings. As g(f) increased, parameter regions of unstable EPs and stable LCs determined for hyperpolarizing loads shrank in the Na(i)-variable system but were enlarged in the Na(i)-fixed system. These findings suggest that 1) I(f) does not enhance but rather attenuates robustness of rabbit SAN cells via facilitating EP stabilization and LC destabilization even in physiological g(f) ranges; and 2) the enhancing effect of I(f) on robustness of pacemaker activity, which could be observed at lower g(f) when Na(i) was fixed, is actually reversed by I(f)-dependent changes in Na(i).

Intracardiac origin of heart rate variability, pacemaker funny current and their possible association with critical illness

Heart rate variability (HRV) is an indirect estimator of autonomic modulation of heart rate and is considered a risk marker in critical illness, particularly in heart failure and severe sepsis. A reduced HRV has been found in critically ill patients and has been associated with neuro-autonomic uncoupling or decreased baroreflex sensitivity. However, results from human and animal experimental studies indicate that intracardiac mechanisms might also be responsible for interbeat fluctuations. These studies have demonstrated that different membrane channel proteins and especially the so-called 'funny' current (If), an hyperpolarization-activated, inward current that drives diastolic depolarization resulting in spontaneous activity in cardiac pacemaker cells, are altered during critical illness. Furthermore, membrane channels kinetics seem to have significant impact upon HRV, whose early decrease might reflect a cellular metabolic stress. In this review article we present research findings regarding intracardiac origin of HRV, at the cellular level and in both isolated sinoatrial node and whole ex vivo heart preparations. In addition, we will review results from various experimental studies that support the interrelation between If and HRV during endotoxemia. We suggest that reduced HRV during sepsis could also be associated with altered pacemaker cell membrane properties, due to ionic current remodeling.

Application of cardiac electrophysiology simulations to pro-arrhythmic safety testing

Concerns over cardiac side effects are the largest single cause of compound attrition during pharmaceutical drug development. For a number of years, biophysically detailed mathematical models of cardiac electrical activity have been used to explore how a compound, interfering with specific ion-channel function, may explain effects at the cell-, tissue- and organ-scales. With the advent of high-throughput screening of multiple ion channels in the wet-lab, and improvements in computational modelling of their effects on cardiac cell activity, more reliable prediction of pro-arrhythmic risk is becoming possible at the earliest stages of drug development. In this paper, we review the current use of biophysically detailed mathematical models of cardiac myocyte electrical activity in drug safety testing, and suggest future directions to employ the full potential of this approach.

A biological relativity view of the relationships between genomes and phenotypes

This article explores the relativistic principle that there is no privileged scale of causality in biology to clarify the relationships between genomes and phenotypes. The idea that genetic causes are primary views the genome as a program. Initially, that view was vindicated by the discovery of mutations and knockouts that have large and specific effects on the phenotype. But we now know that these form the minority of cases. Many changes at the genome level are buffered by robust networks of interactions in cells, tissues and organs. The 'differential' view of genetics therefore fails because it is too restrictive. An 'integral' view, using reverse engineering from systems biological models to quantify contributions to function, can solve this problem. The article concludes by showing that far from breaking the supervenience principle, downward causation requires that it should be obeyed.

Mechano-sensitivity of cardiac pacemaker function: pathophysiological relevance, experimental implications, and conceptual integration with other mechanisms of rhythmicity

Cardiac pacemaker cells exhibit spontaneous, rhythmic electrical excitation, termed automaticity. This automatic initiation of action potentials requires spontaneous diastolic depolarisation, whose rate determines normal rhythm generation in the heart. Pacemaker mechanisms have been split recently into: (i) cyclic changes in trans-sarcolemmal ion flows (termed the ‘membrane-clock’), and (ii) rhythmic intracellular calcium cycling (the ‘calcium-clock’). These two ‘clocks’ undoubtedly interact, as trans-sarcolemmal currents involved in pacemaking include calcium-carrying mechanisms, while intracellular calcium cycling requires trans-sarcolemmal ion flux as the mechanism by which it affects membrane potential. The split into separate ‘clocks’ is, therefore, somewhat arbitrary. Nonetheless, the ‘clock’ metaphor has been conceptually stimulating, in particular since there is evidence to support the view that either ‘clock’ could be sufficient in principle to set the rate of pacemaker activation.

Of course, the same has also been shown for sub-sets of ‘membrane-clock’ ion currents, illustrating the redundancy of mechanisms involved in maintaining such basic functionality as the heartbeat, a theme that is common for vital physiological systems.

Following the conceptual path of identifying individual groups of sub-mechanisms, it is important to remember that the heart is able to adapt pacemaker rate to changes in haemodynamic load, even after isolation or transplantation, and on a beat-by-beat basis. Neither the ‘membrane-’ nor the ‘calcium-clock’ do, as such, inherently account for this rapid adaptation to circulatory demand (cellular Ca²⁺ balance changes over multiple beats, while variation of sarcolemmal ion channel presence takes even longer). This suggests that a third set of mechanisms must be involved in setting the pace. These mechanisms are characterised by their sensitivity to the cyclically changing mechanical environment, and – in analogy to the above terminology – this might be considered a ‘mechanics-clock’.

In this review, we discuss possible roles of mechano-sensitive mechanisms for the entrainment of membrane current dynamics and calcium-handling. This can occur directly via stretch-activation of mechano-sensitive ion channels in the sarcolemma and/or in intracellular membrane compartments, as well as by modulation of ‘standard’ components of the ‘membrane-’ or ‘calcium-clock’. Together, these mechanisms allow rapid adaptation to changes in haemodynamic load, on a beat-by-beat basis.

Additional relevance arises from the fact that mechano-sensitivity of pacemaking may help to explain pacemaker dysfunction in mechanically over- or under-loaded tissue. As the combined contributions of the various underlying oscillatory mechanisms are integrated at the pacemaker cell level into a single output – a train of pacemaker action potentials – we will not adhere to a metaphor that implies separate time-keeping units (‘clocks’), and rather focus on cardiac pacemaking as the result of interactions of a set of coupled oscillators, whose individual contributions vary depending on the pathophysiological context. We conclude by considering the utility and limitations of viewing the pacemaker as a coupled system of voltage-, calcium-, and mechanics-modulated oscillators that, by integrating a multitude of inputs, offers the high level of functional redundancy that is vitally important for cardiac automaticity.

The "funny" current (I(f)) inhibition by ivabradine at membrane potentials encompassing spontaneous depolarization in pacemaker cells

Recent clinical trials have shown that ivabradine (IVA), a drug that inhibits the funny current (I(f)) in isolated sinoatrial nodal cells (SANC), decreases heart rate and reduces morbidity and mortality in patients with cardiovascular diseases. While IVA inhibits I(f), this effect has been reported at essentially unphysiological voltages, i.e., those more negative than the spontaneous diastolic depolarization (DD) between action potentials (APs). We tested the relative potency of IVA to block I(f) over a wide range of membrane potentials, including those that encompass DD governing to the SANC spontaneous firing rate. A clinically relevant IVA concentration of 3 μM to single, isolated rabbit SANC slowed the spontaneous AP firing rate by 15%. During voltage clamp the maximal I(f) was 18 ± 3 pA/pF (at -120 mV) and the maximal I(f) reduction by IVA was 60 ± 8% observed at -92 ± 4 mV. At the maximal diastolic depolarization (~-60 mV) I(f) amplitude was only -2.9 ± 0.4 pA/pF, and was reduced by only 41 ± 6% by IVA. Thus, I(f) amplitude and its inhibition by IVA at physiologically relevant membrane potentials are substantially less than that at unphysiological (hyperpolarized) membrane potentials. This novel finding more accurately describes how IVA affects SANC function and is of direct relevance to numerical modeling of SANC automaticity.

HCN channels and heart rate

Hyperpolarization and Cyclic Nucleotide (HCN) -gated channels represent the molecular correlates of the "funny" pacemaker current (I(f)), a current activated by hyperpolarization and considered able to influence the sinus node function in generating cardiac impulses. HCN channels are a family of six transmembrane domain, single pore-loop, hyperpolarization activated, non-selective cation channels. This channel family comprises four members: HCN1-4, but there is a general agreement to consider HCN4 as the main isoform able to control heart rate. This review aims to summarize advanced insights into the structure, function and cellular regulation of HCN channels in order to better understand the role of such channels in regulating heart rate and heart function in normal and pathological conditions. Therefore, we evaluated the possible therapeutic application of the selective HCN channels blockers in heart rate control.

A theory of biological relativity: no privileged level of causation

Must higher level biological processes always be derivable from lower level data and mechanisms, as assumed by the idea that an organism is completely defined by its genome? Or are higher level properties necessarily also causes of lower level behaviour, involving actions and interactions both ways? This article uses modelling of the heart, and its experimental basis, to show that downward causation is necessary and that this form of causation can be represented as the influences of initial and boundary conditions on the solutions of the differential equations used to represent the lower level processes. These insights are then generalized. A priori, there is no privileged level of causation. The relations between this form of 'biological relativity' and forms of relativity in physics are discussed. Biological relativity can be seen as an extension of the relativity principle by avoiding the assumption that there is a privileged scale at which biological functions are determined.

A model of cellular cardiac-neural coupling that captures the sympathetic control of sinoatrial node excitability in normotensive and hypertensive rats

Hypertension is associated with sympathetic hyperactivity. To represent this neural-myocyte coupling, and to elucidate the mechanisms underlying sympathetic control of the cardiac pacemaker, we developed a new (to our knowledge) cellular mathematical model that incorporates signaling information from cell-to-cell communications between the sympathetic varicosity and sinoatrial node (SAN) in both normotensive (WKY) and hypertensive (SHR) rats. Features of the model include 1), a description of pacemaker activity with specific ion-channel functions and Ca(2+) handling elements; 2), dynamic β-adrenergic modulation of the excitation of the SAN; 3), representation of ionic activity of sympathetic varicosity with NE release dynamics; and 4), coupling of the varicosity model to the SAN model to simulate presynaptic transmitter release driving postsynaptic excitability. This framework captures neural-myocyte coupling and the modulation of pacemaking by nitric oxide and cyclic GMP. It also reproduces the chronotropic response to brief sympathetic stimulations. Finally, the SHR model quantitatively suggests that the impairment of cyclic GMP regulation at both sides of the sympathetic cleft is crucial for development of the autonomic phenotype observed in hypertension.

Differential and integral views of genetics in computational systems biology

This article uses an integrative systems biological view of the relationship between genotypes and phenotypes to clarify some conceptual problems in biological debates about causality. The differential (gene-centric) view is incomplete in a sense analogous to using differentiation without integration in mathematics. Differences in genotype are frequently not reflected in significant differences in phenotype as they are buffered by networks of molecular interactions capable of substituting an alternative pathway to achieve a given phenotype characteristic when one pathway is removed. Those networks integrate the influences of many genes on each phenotype so that the effect of a modification in DNA depends on the context in which it occurs. Mathematical modelling of these interactions can help to understand the mechanisms of buffering and the contextual-dependence of phenotypic outcome, and so to represent correctly and quantitatively the relations between genomes and phenotypes. By incorporating all the causal factors in generating a phenotype, this approach also highlights the role of non-DNA forms of inheritance, and of the interactions at multiple levels.

The anatomy and physiology of the sinoatrial node--a contemporary review

The sinoatrial node is the primary pacemaker of the heart. Nodal dysfunction with aging, heart failure, atrial fibrillation, and even endurance athletic training can lead to a wide variety of pathological clinical syndromes. Recent work utilizing molecular markers to map the extent of the node, along with the delineation of a novel paranodal area intermediate in characteristics between the node and the surrounding atrial muscle, has shown that pacemaker tissue is more widely spread in the right atrium than previously appreciated. This can explain the phenomenon of a "wandering pacemaker" and concomitant changes in the P-wave morphology. Extensive knowledge now exists regarding the molecular architecture of the node (in particular, the expression of ion channels) and how this relates to pacemaking. This review is an up-to-date summary of the current state of our appreciation of the above topics.

Integrative modeling of the cardiac ventricular myocyte

Cardiac electrophysiology is a discipline with a rich 50-year history of experimental research coupled with integrative modeling which has enabled us to achieve a quantitative understanding of the relationships between molecular function and the integrated behavior of the cardiac myocyte in health and disease. In this paper, we review the development of integrative computational models of the cardiac myocyte. We begin with a historical overview of key cardiac cell models that helped shape the field. We then narrow our focus to models of the cardiac ventricular myocyte and describe these models in the context of their subcellular functional systems including dynamic models of voltage-gated ion channels, mitochondrial energy production, ATP-dependent and electrogenic membrane transporters, intracellular Ca dynamics, mechanical contraction, and regulatory signal transduction pathways. We describe key advances and limitations of the models as well as point to new directions for future modeling research. WIREs Syst Biol Med 2011 3 392-413 DOI: 10.1002/wsbm.122

Tachybradycardia in the isolated canine right atrium induced by chronic sympathetic stimulation and pacemaker current inhibition

The mechanisms of sinoatrial node (SAN) dysfunction in patients with chronically elevated sympathetic tone and reduced pacemaker current (I(f); such as heart failure) are poorly understood. We simultaneously mapped membrane potential and intracellular Ca(2+) in the Langendorff-perfused canine right atrium (RA). Blockade of either I(f) (ZD-7288) or sarcoplasmic reticulum Ca(2+) release (ryanodine) alone decreased heart rate by 8% (n = 3) and 16% (n = 3), respectively. Combined treatment of ZD-7288 and ryanodine consistently resulted in prolonged (> or =3 s) sinus pauses (PSPs) (n = 4). However, the middle SAN remained as the leading pacemaking site after these treatments. Prolonged exposure with isoproterenol (0.01 micromol/l) followed by ZD-7288 completely suppressed SAN but triggered recurrent ectopic atrial tachycardia. Cessation of tachycardia was followed by PSPs in five of eight RAs. Isoproterenol initially increased heart rate by 75% from baseline with late diastolic intracellular Ca(2+) elevation (LDCAE) from the superior SAN. However, after a prolonged isoproterenol infusion, LDCAE disappeared in the superior SAN, the leading pacemaker shifted to the inferior SAN, and the rate reduced to 52% above baseline. Caffeine (2 ml, 20 mmol/l) injection after a prolonged isoproterenol infusion produced LDCAE in the SAN and accelerated the SAN rate, ruling out sarcoplasmic reticulum Ca(2+) depletion as a cause of Ca(2+) clock malfunction. We conclude that in an isolated canine RA preparation, chronically elevated sympathetic tone results in abnormal pacemaking hierarchy in the RA, including suppression of the superior SAN and enhanced pacemaking from ectopic sites. Combined malfunction of both membrane and Ca(2+) clocks underlies the mechanisms of PSPs.

The systems biology approach to drug development: application to toxicity assessment of cardiac drugs

Side effects account for most of the instances of failure of candidate drugs at late stages of development. These development failures contribute to the exorbitant cost of bringing new compounds to market: a single withdrawal can represent a loss of more than $1 billion. Many unwanted actions of drugs affect the heart, resulting in potentially proarrhythmic alteration of ion channel function. Because these can be fatal, potential electrophysiological cardiotoxicity is among the most stringent exclusion criteria in the licensing process.

Roles of hyperpolarization-activated current If in sinoatrial node pacemaking: insights from bifurcation analysis of mathematical models

To elucidate the roles of hyperpolarization-activated current (I(f)) in sinoatrial node (SAN) pacemaking, we theoretically investigated 1) the effects of I(f) on stability and bifurcation during hyperpolarization of SAN cells; 2) combined effects of I(f) and the sustained inward current (I(st)) or Na(+) channel current (I(Na)) on robustness of pacemaking against hyperpolarization; and 3) whether blocking I(f) abolishes pacemaker activity under certain conditions. Bifurcation analyses were performed for mathematical models of rabbit SAN cells; equilibrium points (EPs), periodic orbits, and their stability were determined as functions of parameters. Unstable steady-state potential region determined with applications of constant bias currents shrunk as I(f) density increased. In the central SAN cell, the critical acetylcholine concentration at which bifurcations, to yield a stable EP and quiescence, occur was increased by smaller I(f), but decreased by larger I(f). In contrast, the critical acetylcholine concentration and conductance of gap junctions between SAN and atrial cells at bifurcations progressively increased with enhancing I(f) in the peripheral SAN cell. These effects of I(f) were significantly attenuated by eliminating I(st) or I(Na), or by accelerating their inactivation. Under hyperpolarized conditions, blocking I(f) abolished SAN pacemaking via bifurcations. These results suggest that 1) I(f) itself cannot destabilize EPs; 2) I(f) improves SAN cell robustness against parasympathetic stimulation via preventing bifurcations in the presence of I(st) or I(Na); 3) I(f) dramatically enhances peripheral cell robustness against electrotonic loads of the atrium in combination with I(Na); and 4) pacemaker activity of hyperpolarized SAN cells could be abolished by blocking I(f).

A novel quantitative explanation for the autonomic modulation of cardiac pacemaker cell automaticity via a dynamic system of sarcolemmal and intracellular proteins

Classical numerical models have attributed the regulation of normal cardiac automaticity in sinoatrial node cells (SANCs) largely to G protein-coupled receptor (GPCR) modulation of sarcolemmal ion currents. More recent experimental evidence, however, has indicated that GPCR modulation of SANCs automaticity involves spontaneous, rhythmic, local Ca(2+) releases (LCRs) from the sarcoplasmic reticulum (SR). We explored the GPCR rate modulation of SANCs using a unique and novel numerical model of SANCs in which Ca(2+)-release characteristics are graded by variations in the SR Ca(2+) pumping capability, mimicking the modulation by phospholamban regulated by cAMP-mediated, PKA-activated signaling. The model faithfully predicted the entire range of physiological chronotropic modulation of SANCs by the activation of beta-adrenergic receptors or cholinergic receptors only when experimentally documented changes of sarcolemmal ion channels are combined with a simultaneous increase/decrease in SR Ca(2+) pumping capability. The novel numerical mechanism of GPCR rate modulation is based on numerous complex synergistic interactions between sarcolemmal and intracellular processes via membrane voltage and Ca(2+). Major interactions include changes of diastolic Na(+)/Ca(2+) exchanger current that couple earlier/later diastolic Ca(2+) releases (predicting the experimentally defined LCR period shift) of increased/decreased amplitude (predicting changes in LCR signal mass, i.e., the product of LCR spatial size, amplitude, and number per cycle) to the diastolic depolarization and ultimately to the spontaneous action potential firing rate. Concomitantly, larger/smaller and more/less frequent activation of L-type Ca(2+) current shifts the cellular Ca(2+) balance to support the respective Ca(2+) cycling changes. In conclusion, our model simulations corroborate recent experimental results in rabbit SANCs pointing to a new paradigm for GPCR heart rate modulation by a complex system of dynamically coupled sarcolemmal and intracellular proteins.

Induction of atrial ectopic beats with calcium release inhibition: Local hierarchy of automaticity in the right atrium

BACKGROUND

Recent evidence indicates that spontaneous sarcoplasmic reticulum (SR) calcium (Ca) release underlies the mechanism of sinoatrial node (SAN) acceleration during beta-stimulation, indicating the importance of the Ca clock in SAN automaticity. Whether or not the same mechanism applies to atrial ectopic pacemakers (AEPs) remains unclear.

OBJECTIVE

The purpose of this study was to assess the mechanism of AEP.

METHODS

We simultaneously mapped intracellular calcium (Ca(i)) and membrane potential in 12 isolated canine right atria. The late diastolic Ca(i) elevation (LDCAE) was used to detect the Ca clock activity. Pharmacological interventions with isoproterenol (ISO), ryanodine, and ZD7288, a blocker of the I(f) membrane current, were performed.

RESULTS

Ryanodine, which inhibits SR Ca release, reduced LDCAE in SAN, resulting in an inferior shift of the pacemaking site. Cycle length increased significantly in a dose-dependent way. In the presence of 3 to 10 mumol/l of ryanodine, ISO infusion consistently induces AEPs from the lower crista terminalis. All ectopic beats continuing over 30 seconds were located at the lower crista terminalis. These AEPs were resistant to ryanodine treatment even at high doses. Subsequent blockade of I(f) inhibited the AEP and resulted in profound bradycardia.

CONCLUSION

Spontaneous SR Ca release underlies ISO-induced increase of superior SAN activity. As compared with SAN, the AEP is less dependent on the Ca clock and more dependent on the membrane clock for its automaticity. AEPs outside the SAN can effectively serve as backup pacemakers when the Ca clock functionality is reduced.

Pacemaker activity of the human sinoatrial node: role of the hyperpolarization-activated current, I(f)

The mechanism of primary, spontaneous cardiac pacemaker activity of the sinoatrial node (SAN) has extensively been studied in several animal species, but is virtually unexplored in man. Understanding the mechanisms of human SAN pacemaker activity is important for developing new therapeutic approaches for controlling the heart rate in the sick sinus syndrome and in diseased myocardium. Here we review the functional role of the hyperpolarization-activated 'funny' current, I(f), in human SAN pacemaker activity. Despite the many animal studies performed over the years, the contribution of I(f) to pacemaker activity is still controversial and not fully established. However, recent clinical data on mutations in the I(f) encoding HCN4 gene, which is thought to be the most abundant isoform of the HCN gene family in SAN, suggest a functional role of I(f) in human pacemaker activity. These clinical findings are supported by recent experimental data from single isolated human SAN cells that provide direct evidence that I(f) contributes to human SAN pacemaker activity. Therefore, controlling heart rate in clinical practice via I(f) blockers offers a valuable approach to lowering heart rate and provides an attractive alternative to conventional treatment for a wide range of patients with confirmed stable angina, while upregulation or artificial expression of I(f) may relieve disease-causing bradycardias.

Synergism of coupled subsarcolemmal Ca2+ clocks and sarcolemmal voltage clocks confers robust and flexible pacemaker function in a novel pacemaker cell model

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.

Genesis and regulation of the heart automaticity

The heart automaticity is a fundamental physiological function in higher organisms. The spontaneous activity is initiated by specialized populations of cardiac cells generating periodical electrical oscillations. The exact cascade of steps initiating the pacemaker cycle in automatic cells has not yet been entirely elucidated. Nevertheless, ion channels and intracellular Ca(2+) signaling are necessary for the proper setting of the pacemaker mechanism. Here, we review the current knowledge on the cellular mechanisms underlying the generation and regulation of cardiac automaticity. We discuss evidence on the functional role of different families of ion channels in cardiac pacemaking and review recent results obtained on genetically engineered mouse strains displaying dysfunction in heart automaticity. Beside ion channels, intracellular Ca(2+) release has been indicated as an important mechanism for promoting automaticity at rest as well as for acceleration of the heart rate under sympathetic nerve input. The potential links between the activity of ion channels and Ca(2+) release will be discussed with the aim to propose an integrated framework of the mechanism of automaticity.

Hyperpolarization-activated cyclic nucleotide-modulated 'HCN' channels confer regular and faster rhythmicity to beating mouse embryonic stem cells

The hyperpolarization-activated cation current (I(f)), and the hyperpolarization-activated cyclic nucleotide-modulated 'HCN' subunits that underlie it, are important components of spontaneous activity in the embryonic mouse heart, but whether they contribute to this activity in mouse embryonic stem cell-derived cardiomyocytes has not been investigated. We address this issue in spontaneously beating cells derived from mouse embryonic stem cells (mESCs) over the course of development in culture. I(f) and action potentials were recorded from single beating cells at early, intermediate and late development stages using perforated whole-cell voltage- and current-clamp techniques. Our data show that the proportion of cells expressing I(f), and the density of I(f) in these cells, increased during development and correlated with action potential frequency and the rate of diastolic depolarization. The I(f) blocker ZD7288 (0.3 microm) reduced I(f) and the beating rate of embryoid bodies. Taken together, the activation kinetics of I(f) and results from Western blots are consistent with the presence of the HCN2 and HCN3 isoforms. At all stages of development, isoproterenol (isoprenaline) and acetylcholine shifted the voltage dependence of I(f) to more positive and negative voltages, respectively, and they also increased and decreased the beating rate of embryonic cell bodies, respectively. Together, the data suggest that current through HCN2 and HCN3 channels confers regular and faster rhythmicity to mESCs, which mirrors the developing embryonic mouse heart, and contributes to modulation of rhythmicity by autonomic stimulation.