The emergence of subcellular pacemaker sites for calcium waves and oscillations

Modeling Calcium Cycling in the Heart: Progress, Pitfalls, and Challenges

Intracellular calcium (Ca) cycling in the heart plays key roles in excitation-contraction coupling and arrhythmogenesis. In cardiac myocytes, the Ca release channels, i.e., the ryanodine receptors (RyRs), are clustered in the sarcoplasmic reticulum membrane, forming Ca release units (CRUs). The RyRs in a CRU act collectively to give rise to discrete Ca release events, called Ca sparks. A cell contains hundreds to thousands of CRUs, diffusively coupled via Ca to form a CRU network. A rich spectrum of spatiotemporal Ca dynamics is observed in cardiac myocytes, including Ca sparks, spark clusters, mini-waves, persistent whole-cell waves, and oscillations. Models of different temporal and spatial scales have been developed to investigate these dynamics. Due to the complexities of the CRU network and the spatiotemporal Ca dynamics, it is challenging to model the Ca cycling dynamics in the cardiac system, particularly at the tissue sales. In this article, we review the progress of modeling of Ca cycling in cardiac systems from single RyRs to the tissue scale, the pros and cons of the current models and different modeling approaches, and the challenges to be tackled in the future.

Insect nephrocyte function is regulated by a store operated calcium entry mechanism controlling endocytosis and Amnionless turnover

Insect nephrocytes are ultrafiltration cells that remove circulating proteins and exogenous toxins from the haemolymph. Experimental disruption of nephrocyte development or function leads to systemic impairment of insect physiology as evidenced by cardiomyopathy, chronic activation of immune signalling and shortening of lifespan. The genetic and structural basis of the nephrocyte's ultrafiltration mechanism is conserved between arthropods and mammals, making them an attractive model for studying human renal function and systemic clearance mechanisms in general. Although dynamic changes to intracellular calcium are fundamental to the function of many cell types, there are currently no studies of intracellular calcium signalling in nephrocytes. In this work we aimed to characterise calcium signalling in the pericardial nephrocytes of Drosophila melanogaster. To achieve this, a genetically encoded calcium reporter (GCaMP6) was expressed in nephrocytes to monitor intracellular calcium both in vivo within larvae and in vitro within dissected adults. Larval nephrocytes exhibited stochastically timed calcium waves. A calcium signal could be initiated in preparations of adult nephrocytes and abolished by EGTA, or the store operated calcium entry (SOCE) blocker 2-APB, as well as RNAi mediated knockdown of the SOCE genes Stim and Orai. Neither the presence of calcium-free buffer nor EGTA affected the binding of the endocytic cargo albumin to nephrocytes but they did impair the subsequent accumulation of albumin within nephrocytes. Pre-treatment with EGTA, calcium-free buffer or 2-APB led to significantly reduced albumin binding. Knock-down of Stim and Orai was non-lethal, caused an increase to nephrocyte size and reduced albumin binding, reduced the abundance of the endocytic cargo receptor Amnionless and disrupted the localisation of Dumbfounded at the filtration slit diaphragm. These data indicate that pericardial nephrocytes exhibit stochastically timed calcium waves in vivo and that SOCE mediates the localisation of the endocytic co-receptor Amnionless. Identifying the signals both up and downstream of SOCE may highlight mechanisms relevant to the renal and excretory functions of a broad range of species, including humans.

Kramers Rate Theory of Pacemaker Dynamics in Noisy Excitable Media

Rhythmic activities, which are usually driven by pacemakers, are common in biological systems. In noisy excitable media, pacemakers are self-organized firing clusters, but the underlying dynamics remains to be elucidated. Here we develop a Kramers rate theory of coupled cells to describe the firing properties of pacemakers and their dependence on coupling strength and system size and dimension. The theory captures accurately the simulation results of tissue models with stochastic Hodgkin-Huxley equations except when transitions from pacemakers to spiral waves occur under weak coupling.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

Size Matters: Ryanodine Receptor Cluster Size Heterogeneity Potentiates Calcium Waves

Ryanodine receptors (RyRs) mediate calcium (Ca)-induced Ca release and intracellular Ca homeostasis. In a cardiac myocyte, RyRs group into clusters of variable size from a few to several hundred RyRs, creating a spatially nonuniform intracellular distribution. It is unclear how heterogeneity of RyR cluster size alters spontaneous sarcoplasmic reticulum (SR) Ca releases (Ca sparks) and arrhythmogenic Ca waves. Here, we tested the impact of heterogeneous RyR cluster size on the initiation of Ca waves. Experimentally, we measured RyR cluster sizes at Ca spark sites in rat ventricular myocytes and further tested functional impacts using a physiologically detailed computational model with spatial and stochastic intracellular Ca dynamics. We found that the spark frequency and amplitude increase nonlinearly with the size of RyR clusters. Larger RyR clusters have lower SR Ca release threshold for local Ca spark initiation and exhibit steeper SR Ca release versus SR Ca load relationship. However, larger RyR clusters tend to lower SR Ca load because of the higher Ca leak rate. Conversely, smaller clusters have a higher threshold and a lower leak, which tends to increase SR Ca load. At the myocyte level, homogeneously large or small RyR clusters limit Ca waves (because of low load for large clusters but low excitability for small clusters). Mixtures of large and small RyR clusters potentiates Ca waves because the enhanced SR Ca load driven by smaller clusters enables Ca wave initiation and propagation from larger RyR clusters. Our study suggests that a spatially heterogeneous distribution of RyR cluster size under pathological conditions may potentiate Ca waves and thus afterdepolarizations and triggered arrhythmias.

Tissue Regeneration without Stem Cell Transplantation: Self-Healing Potential from Ancestral Chemistry and Physical Energies

The human body constantly regenerates after damage due to the self-renewing and differentiating properties of its resident stem cells. To recover the damaged tissues and regenerate functional organs, scientific research in the field of regenerative medicine is firmly trying to understand the molecular mechanisms through which the regenerative potential of stem cells may be unfolded into a clinical application. The finding that some organisms are capable of regenerative processes and the study of conserved evolutionary patterns in tissue regeneration may lead to the identification of natural molecules of ancestral species capable to extend their regenerative potential to human tissues. Such a possibility has also been strongly suggested as a result of the use of physical energies, such as electromagnetic fields and mechanical vibrations in human adult stem cells. Results from scientific studies on stem cell modulation confirm the possibility to afford a chemical manipulation of stem cell fate in vitro and pave the way to the use of natural molecules, as well as electromagnetic fields and mechanical vibrations to target human stem cells in their niche inside the body, enhancing human natural ability for self-healing.

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

Aims

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

Methods and results

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

Conclusion

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

Transverse tubular network structures in the genesis of intracellular calcium alternans and triggered activity in cardiac cells

RATIONALE

The major role of a transverse-tubular (TT) network in a cardiac cell is to facilitate effective excitation-contraction coupling and signaling. The TT network structures are heterogeneous within a single cell, and vary between different types of cells and species. They are also remodeled in cardiac diseases. However, how different TT network structures predispose cardiac cells to arrhythmogenesis remains to be revealed.

OBJECTIVE

To systematically investigate the roles of TT network structure and the underlying mechanisms in the genesis of intracellular calcium (Ca2+) alternans and triggered activity (TA).

METHODS AND RESULTS

Based on recent experimental observations, different TT network structures, including uniformly and non-uniformly random TT distributions, were modeled in a cardiac cell model consisting of a three-dimensional network of Ca2+ release units (CRUs). Our simulations showed that both Ca2+ alternans and Ca2+ wave-mediated TA were promoted when the fraction of orphaned CRUs was in an intermediate range, but suppressed in cells exhibiting either well-organized TT networks or low TT densities. Ca2+ alternans and TA could be promoted by low TT densities when the cells were small or the CRU coupling was strong. Both alternans and TA occurred more easily in uniformly random TT networks than in non-uniformly random TT networks. Subcellular spatially discordant Ca2+ alternans was promoted by non-uniformly random TT networks but suppressed by increasing CRU coupling strength. These mechanistic insights provide a holistic understanding of the effects of TT network structure on the susceptibility to arrhythmogenesis.

CONCLUSIONS

The TT network plays important roles in promoting Ca2+ alternans and TA, and different TT network structures may predispose cardiac cells differently to arrhythmogenesis.

Noise-induced effects on multicellular biopacemaker spontaneous activity: Differences between weak and strong pacemaker cells

Self-organization of spontaneous activity of a network of active elements is important to the general theory of reaction-diffusion systems as well as for pacemaking activity to initiate beating of the heart. Monolayer cultures of neonatal rat ventricular myocytes, consisting of resting and pacemaker cells, exhibit spontaneous activation of their electrical activity. Similarly, one proposed approach to the development of biopacemakers as an alternative to electronic pacemakers for cardiac therapy is based on heterogeneous cardiac cells with resting and spontaneously beating phenotypes. However, the combined effect of pacemaker characteristics, density, and spatial distribution of the pacemaker cells on spontaneous activity is unknown. Using a simple stochastic pattern formation algorithm, we previously showed a clear nonlinear dependency of spontaneous activity (occurrence and amplitude of spontaneous period) on the spatial patterns of pacemaker cells. In this study, we show that this behavior is dependent on the pacemaker cell characteristics, with weaker pacemaker cells requiring higher density and larger clusters to sustain multicellular activity. These multicellular structures also demonstrated an increased sensitivity to voltage noise that favored spontaneous activity at lower density while increasing temporal variation in the period of activity. This information will help researchers overcome the current limitations of biopacemakers.

Ambiguous interactions between diastolic and SR Ca2+ in the regulation of cardiac Ca2+ release

Sobie et al. highlight unresolved issues concerning the regulation of sarcoplasmic reticulum calcium release in cardiac myocytes.

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

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

Calcium Signaling in Interstitial Cells: Focus on Telocytes

In this review, we describe the current knowledge on calcium signaling pathways in interstitial cells with a special focus on interstitial cells of Cajal (ICCs), interstitial Cajal-like cells (ICLCs), and telocytes. In detail, we present the generation of Ca²⁺ oscillations, the inositol triphosphate (IP₃)/Ca²⁺ signaling pathway and modulation exerted by cytokines and vasoactive agents on calcium signaling in interstitial cells. We discuss the physiology and alterations of calcium signaling in interstitial cells, and in particular in telocytes. We describe the physiological contribution of calcium signaling in interstitial cells to the pacemaking activity (e.g., intestinal, urinary, uterine or vascular pacemaking activity) and to the reproductive function. We also present the pathological contribution of calcium signaling in interstitial cells to the aortic valve calcification or intestinal inflammation. Moreover, we summarize the current knowledge of the role played by calcium signaling in telocytes in the uterine, cardiac and urinary physiology, and also in various pathologies, including immune response, uterine and cardiac pathologies.

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

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

Nonlinear onset of calcium wave propagation in cardiac cells

Spontaneous calcium (Ca) waves in cardiac myocytes are known to underlie a wide range of cardiac arrhythmias. However, it is not understood which physiological parameters determine the onset of waves. In this study, we explore the relationship between Ca signaling between ion channels and the nucleation of Ca waves. In particular, we apply a master equation approach to analyze the stochastic interaction between neighboring clusters of ryanodine receptor (RyR) channels. Using this analysis, we show that signaling between clusters can be described as a barrier hopping process with exponential sensitivity to system parameters. A consequence of this feature is that the probability that Ca release at a cluster induces release at a neighboring cluster exhibits a sigmoid dependence on the Ca content in the cell. This nonlinearity originates from the regulation of RyR opening due to more than one Ca ion binding site, in conjunction with Ca mediated cooperativity between RyR channels in clusters. We apply a spatially distributed stochastic model of Ca cycling to analyze the physiological consequences of this nonlinearity, and show that it explains the sharp onset of Ca wave nucleation in cardiac cells. Furthermore, we show that this sharp onset can serve as a mechanism for Ca alternans under physiologically relevant conditions. Thus our findings identify the nonlinear features of Ca signaling which potentially underlie the onset of Ca waves and Ca alternans in cardiac cells.

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

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

Role of connectivity and fluctuations in the nucleation of calcium waves in cardiac cells

Spontaneous calcium release (SCR) occurs when ion channel fluctuations lead to the nucleation of calcium waves in cardiac cells. This phenomenon is important since it has been implicated as a cause of various cardiac arrhythmias. However, to date, it is not understood what determines the timing and location of spontaneous calcium waves within cells. Here, we analyze a simplified model of SCR in which calcium release is modeled as a stochastic processes on a two-dimensional network of randomly distributed sites. Using this model we identify the essential parameters describing the system and compute the phase diagram. In particular, we identify a critical line which separates pinned and propagating fronts, and show that above this line wave nucleation is governed by fluctuations and the spatial connectivity of calcium release units. Using a mean-field analysis we show that the sites of wave nucleation are predicted by localized eigenvectors of a matrix representing the network connectivity of release sites. This result provides insight on the interplay between connectivity and fluctuations in the genesis of SCR in cardiac myocytes.

Acute reversal of phospholamban inhibition facilitates the rhythmic whole-cell propagating calcium waves in isolated ventricular myocytes

Phospholamban (PLB) inhibits the activity of cardiac sarcoplasmic reticulum (SR) Ca(2+)-ATPase (SERCA2a). Phosphorylation of PLB during sympathetic activation reverses SERCA2a inhibition, increasing SR Ca(2+) uptake. However, sympathetic activation also modulates multiple other intracellular targets in ventricular myocytes (VMs), making it impossible to determine the specific effects of the reversal of PLB inhibition on the spontaneous SR Ca(2+) release. Therefore, it remains unclear how PLB regulates rhythmic activity in VMs. Here, we used the Fab fragment of 2D12, a monoclonal anti-PLB antibody, to test how acute reversal of PLB inhibition affects the spontaneous SR Ca(2+) release in normal VMs. Ca(2+) sparks and spontaneous Ca(2+) waves (SCWs) were recorded in the line-scan mode of confocal microscopy using the Ca(2+) fluorescent dye Fluo-4 in isolated permeabilized mouse VMs. Fab, which reverses PLB inhibition, significantly increased the frequency, amplitude, and spatial/temporal spread of Ca(2+) sparks in VMs exposed to 50 nM free [Ca(2+)]. At physiological diastolic free [Ca(2+)] (100-200 nM), Fab facilitated the formation of whole-cell propagating SCWs. At higher free [Ca(2+)], Fab increased the frequency and velocity, but decreased the decay time of the SCWs. cAMP had little additional effect on the frequency or morphology of Ca(2+) sparks or SCWs after Fab addition. These findings were complemented by computer simulations. In conclusion, acute reversal of PLB inhibition alone significantly increased the spontaneous SR Ca(2+) release, leading to the facilitation and organization of whole-cell propagating SCWs in normal VMs. PLB thus plays a key role in subcellular Ca(2+) dynamics and rhythmic activity of VMs.

Mechanisms of ventricular arrhythmias: from molecular fluctuations to electrical turbulence

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

Nonlinear and Stochastic Dynamics in the Heart

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