Human atrial fibrillation: insights from computational electrophysiological models

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Populations of in silico myocytes and tissues reveal synergy of multiatrial-predominant K+ -current block in atrial fibrillation

BACKGROUND AND PURPOSE

Pharmacotherapy of atrial fibrillation (AF), the most common cardiac arrhythmia, remains unsatisfactory due to low efficacy and safety concerns. New therapeutic strategies target atrial-predominant ion-channels and involve multichannel block (poly)therapy. As AF is characterized by rapid and irregular atrial activations, compounds displaying potent antiarrhythmic effects at fast and minimal effects at slow rates are desirable. We present a novel systems pharmacology framework to quantitatively evaluate synergistic anti-AF effects of combined block of multiple atrial-predominant K⁺ currents (ultra-rapid delayed rectifier K⁺ current, I_Kur , small conductance Ca²⁺ -activated K⁺ current, I_KCa , K_2P 3.1 2-pore-domain K⁺ current, I_K2P ) in AF.

EXPERIMENTAL APPROACH

We constructed experimentally calibrated populations of virtual atrial myocyte models in normal sinus rhythm and AF-remodelled conditions using two distinct, well-established atrial models. Sensitivity analyses on our atrial populations was used to investigate the rate dependence of action potential duration (APD) changes due to blocking I_Kur , I_K2P or I_KCa and interactions caused by blocking of these currents in modulating APD. Block was simulated in both single myocytes and one-dimensional tissue strands to confirm insights from the sensitivity analyses and examine anti-arrhythmic effects of multi-atrial-predominant K⁺ current block in single cells and coupled tissue.

KEY RESULTS

In both virtual atrial myocytes and tissues, multiple atrial-predominant K⁺ -current block promoted favourable positive rate-dependent APD prolongation and displayed positive rate-dependent synergy, that is, increasing synergistic antiarrhythmic effects at fast pacing versus slow rates.

CONCLUSION AND IMPLICATIONS

Simultaneous block of multiple atrial-predominant K⁺ currents may be a valuable antiarrhythmic pharmacotherapeutic strategy for AF.

Understanding the Beat-to-Beat Variations of P-Waves Morphologies in AF Patients During Sinus Rhythm: A Scoping Review of the Atrial Simulation Studies

The remarkable advances in high-performance computing and the resulting increase of the computational power have the potential to leverage computational cardiology toward improving our understanding of the pathophysiological mechanisms of arrhythmias, such as Atrial Fibrillation (AF). In AF, a complex interaction between various triggers and the atrial substrate is considered to be the leading cause of AF initiation and perpetuation. In electrocardiography (ECG), P-wave is supposed to reflect atrial depolarization. It has been found that even during sinus rhythm (SR), multiple P-wave morphologies are present in AF patients with a history of AF, suggesting a higher dispersion of the conduction route in this population. In this scoping review, we focused on the mechanisms which modify the electrical substrate of the atria in AF patients, while investigating the existence of computational models that simulate the propagation of the electrical signal through different routes. The adopted review methodology is based on a structured analytical framework which includes the extraction of the keywords based on an initial limited bibliographic search, the extensive literature search and finally the identification of relevant articles based on the reference list of the studies. The leading mechanisms identified were classified according to their scale, spanning from mechanisms in the cell, tissue or organ level, and the produced outputs. The computational modeling approaches for each of the factors that influence the initiation and the perpetuation of AF are presented here to provide a clear overview of the existing literature. Several levels of categorization were adopted while the studies which aim to translate their findings to ECG phenotyping are highlighted. The results denote the availability of multiple models, which are appropriate under specific conditions. However, the consideration of complex scenarios taking into account multiple spatiotemporal scales, personalization of electrophysiological and anatomical models and the reproducibility in terms of ECG phenotyping has only partially been tackled so far.

A calcium transport mechanism for atrial fibrillation in Tbx5-mutant mice

Risk for Atrial Fibrillation (AF), the most common human arrhythmia, has a major genetic component. The T-box transcription factor TBX5 influences human AF risk, and adult-specific Tbx5-mutant mice demonstrate spontaneous AF. We report that TBX5 is critical for cellular Ca²⁺ homeostasis, providing a molecular mechanism underlying the genetic implication of TBX5 in AF. We show that cardiomyocyte action potential (AP) abnormalities in Tbx5-deficient atrial cardiomyocytes are caused by a decreased sarcoplasmic reticulum (SR) Ca²⁺ ATPase (SERCA2)-mediated SR calcium uptake which was balanced by enhanced trans-sarcolemmal calcium fluxes (calcium current and sodium/calcium exchanger), providing mechanisms for triggered activity. The AP defects, cardiomyocyte ectopy, and AF caused by TBX5 deficiency were rescued by phospholamban removal, which normalized SERCA function. These results directly link transcriptional control of SERCA2 activity, depressed SR Ca²⁺ sequestration, enhanced trans-sarcolemmal calcium fluxes, and AF, establishing a mechanism underlying the genetic basis for a Ca²⁺-dependent pathway for AF risk.

The human heart contains four distinct chambers that work together to pump blood around the body. In individuals with a condition called atrial fibrillation, two of the chambers (known as the atria) beat irregularly and are unable to push all the blood they hold into the other two chambers of the heart. This can cause heart failure and increases the likelihood of blood clots, which may lead to stroke and heart attacks.

Small molecules called calcium ions play a crucial role in regulating how and when the atria contract by driving electrical activity in heart cells. To contract the atria, a storage compartment within heart cells known as the sarcoplasmic reticulum releases calcium ions into the main compartment of the cells. Calcium ions also enter the cell from the surrounding tissue. As the atria relax, calcium ions are pumped back into the sarcoplasmic reticulum or out of the cell by specific transport proteins.

Individuals with mutations in a gene called Tbx5 are more likely to develop atrial fibrillation than other people, but it was not clear how such gene mutations contribute to the disease. Here, Dai, Laforest et al. used mice with a mutation in the Tbx5 gene to study how defects in Tbx5 affect electrical activity in heart cells.

The experiments found that the Tbx5 gene was critical for calcium ions to drive normal electrical activity in mouse heart cells. Compared with heart cells from normal mice, the heart cells from the mutant mice had decreased flow of calcium ions into the sarcoplasmic reticulum and increased flow of calcium ions out of the cell.

These findings provide a direct link between atrial fibrillation and the flow of calcium ions in heart cells. Together with previous work, these findings indicate that multiple different mechanisms could lead to atrial fibrillation, but that many of these involve changes in the flow of calcium ions. Therefore, personalized medicine, where clinicians uncover the specific mechanisms responsible for atrial fibrillation in individual patients, may play an important role in treating this condition in the future.

Computational modeling: What does it tell us about atrial fibrillation therapy?

Atrial fibrillation (AF) is a complex cardiac arrhythmia with diverse etiology that negatively affects morbidity and mortality of millions of patients. Technological and experimental advances have provided a wealth of information on the pathogenesis of AF, highlighting a multitude of mechanisms involved in arrhythmia initiation and maintenance, and disease progression. However, it remains challenging to identify the predominant mechanisms for specific subgroups of AF patients, which, together with an incomplete understanding of the pleiotropic effects of antiarrhythmic therapies, likely contributes to the suboptimal efficacy of current antiarrhythmic approaches. Computer modeling of cardiac electrophysiology has advanced in parallel to experimental research and provides an integrative framework to attempt to overcome some of these challenges. Multi-scale cardiac modeling and simulation integrate structural and functional data from experimental and clinical work with knowledge of atrial electrophysiological mechanisms and dynamics, thereby improving our understanding of AF mechanisms and therapy. In this review, we describe recent advances in our quantitative understanding of AF through mathematical models. We discuss computational modeling of AF mechanisms and therapy using detailed, mechanistic cell/tissue-level models, including approaches to incorporate variability in patient populations. We also highlight efforts using whole-atria models to improve catheter ablation therapies. Finally, we describe recent efforts and suggest future extensions to model clinical concepts of AF using patient-level models.

Distinct Occurrence of Proarrhythmic Afterdepolarizations in Atrial Versus Ventricular Cardiomyocytes: Implications for Translational Research on Atrial Arrhythmia

Background: Principal mechanisms of arrhythmia have been derived from ventricular but not atrial cardiomyocytes of animal models despite higher prevalence of atrial arrhythmia (e.g., atrial fibrillation). Due to significant ultrastructural and functional differences, a simple transfer of ventricular proneness toward arrhythmia to atrial arrhythmia is critical. The use of murine models in arrhythmia research is widespread, despite known translational limitations. We here directly compare atrial and ventricular mechanisms of arrhythmia to identify critical differences that should be considered in murine models for development of antiarrhythmic strategies for atrial arrhythmia.

Methods and Results: Isolated murine atrial and ventricular myocytes were analyzed by wide field microscopy and subjected to a proarrhythmic protocol during patch-clamp experiments. As expected, the spindle shaped atrial myocytes showed decreased cell area and membrane capacitance compared to the rectangular shaped ventricular myocytes. Though delayed afterdepolarizations (DADs) could be evoked in a similar fraction of both cell types (80% of cells each), these led significantly more often to the occurrence of spontaneous action potentials (sAPs) in ventricular myocytes. Interestingly, numerous early afterdepolarizations (EADs) were observed in the majority of ventricular myocytes, but there was no EAD in any atrial myocyte (EADs per cell; atrial myocytes: 0 ± 0; n = 25/12 animals; ventricular myocytes: 1.5 [0–43]; n = 20/12 animals; p < 0.05). At the same time, the action potential duration to 90% decay (APD₉₀) was unaltered and the APD₅₀ even increased in atrial versus ventricular myocytes. However, the depolarizing L-type Ca²⁺ current (I_Ca) and Na⁺/Ca²⁺-exchanger inward current (I_NCX) were significantly smaller in atrial versus ventricular myocytes.

Conclusion: In mice, atrial myocytes exhibit a substantially distinct occurrence of proarrhythmic afterdepolarizations compared to ventricular myocytes, since they are in a similar manner susceptible to DADs but interestingly seem to be protected against EADs and show less sAPs. Key factors in the generation of EADs like I_Ca and I_NCX were significantly reduced in atrial versus ventricular myocytes, which may offer a mechanistic explanation for the observed protection against EADs. These findings may be of relevance for current studies on atrial level in murine models to develop targeted strategies for the treatment of atrial arrhythmia.

In Silico Assessment of Efficacy and Safety of IKur Inhibitors in Chronic Atrial Fibrillation: Role of Kinetics and State-Dependence of Drug Binding

Current pharmacological therapy against atrial fibrillation (AF), the most common cardiac arrhythmia, is limited by moderate efficacy and adverse side effects including ventricular proarrhythmia and organ toxicity. One way to circumvent the former is to target ion channels that are predominantly expressed in atria vs. ventricles, such as K_V1.5, carrying the ultra-rapid delayed-rectifier K⁺ current (I_Kur). Recently, we used an in silico strategy to define optimal K_V1.5-targeting drug characteristics, including kinetics and state-dependent binding, that maximize AF-selectivity in human atrial cardiomyocytes in normal sinus rhythm (nSR). However, because of evidence for I_Kur being strongly diminished in long-standing persistent (chronic) AF (cAF), the therapeutic potential of drugs targeting I_Kur may be limited in cAF patients. Here, we sought to simulate the efficacy (and safety) of I_Kur inhibitors in cAF conditions. To this end, we utilized sensitivity analysis of our human atrial cardiomyocyte model to assess the importance of I_Kur for atrial cardiomyocyte electrophysiological properties, simulated hundreds of theoretical drugs to reveal those exhibiting anti-AF selectivity, and compared the results obtained in cAF with those in nSR. We found that despite being downregulated, I_Kur contributes more prominently to action potential (AP) and effective refractory period (ERP) duration in cAF vs. nSR, with ideal drugs improving atrial electrophysiology (e.g., ERP prolongation) more in cAF than in nSR. Notably, the trajectory of the AP during cAF is such that more I_Kur is available during the more depolarized plateau potential. Furthermore, I_Kur block in cAF has less cardiotoxic effects (e.g., AP duration not exceeding nSR values) and can increase Ca²⁺ transient amplitude thereby enhancing atrial contractility. We propose that in silico strategies such as that presented here should be combined with in vitro and in vivo assays to validate model predictions and facilitate the ongoing search for novel agents against AF.

Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity

The K_V1.5 potassium channel, which underlies the ultra-rapid delayed-rectifier current (I_Kur) and is predominantly expressed in atria vs. ventricles, has emerged as a promising target to treat atrial fibrillation (AF). However, while numerous K_V1.5-selective compounds have been screened, characterized, and tested in various animal models of AF, evidence of antiarrhythmic efficacy in humans is still lacking. Moreover, current guidelines for pre-clinical assessment of candidate drugs heavily rely on steady-state concentration-response curves or IC₅₀ values, which can overlook adverse cardiotoxic effects. We sought to investigate the effects of kinetics and state-dependent binding of I_Kur-targeting drugs on atrial electrophysiology in silico and reveal the ideal properties of I_Kur blockers that maximize anti-AF efficacy and minimize pro-arrhythmic risk. To this aim, we developed a new Markov model of I_Kur that describes K_V1.5 gating based on experimental voltage-clamp data in atrial myocytes from patient right-atrial samples in normal sinus rhythm. We extended the I_Kur formulation to account for state-specificity and kinetics of K_V1.5-drug interactions and incorporated it into our human atrial cell model. We simulated 1- and 3-Hz pacing protocols in drug-free conditions and with a [drug] equal to the IC₅₀ value. The effects of binding and unbinding kinetics were determined by examining permutations of the forward (k_on) and reverse (k_off) binding rates to the closed, open, and inactivated states of the K_V1.5 channel. We identified a subset of ideal drugs exhibiting anti-AF electrophysiological parameter changes at fast pacing rates (effective refractory period prolongation), while having little effect on normal sinus rhythm (limited action potential prolongation). Our results highlight that accurately accounting for channel interactions with drugs, including kinetics and state-dependent binding, is critical for developing safer and more effective pharmacological anti-AF options.

Potassium currents in the heart: functional roles in repolarization, arrhythmia and therapeutics

This is the second of the two White Papers from the fourth UC Davis Cardiovascular Symposium Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias (3-4 March 2016), a biennial event that brings together leading experts in different fields of cardiovascular research. The theme of the 2016 symposium was 'K⁺ channels and regulation', and the objectives of the conference were severalfold: (1) to identify current knowledge gaps; (2) to understand what may go wrong in the diseased heart and why; (3) to identify possible novel therapeutic targets; and (4) to further the development of systems biology approaches to decipher the molecular mechanisms and treatment of cardiac arrhythmias. The sessions of the Symposium focusing on the functional roles of the cardiac K⁺ channel in health and disease, as well as K⁺ channels as therapeutic targets, were contributed by Ye Chen-Izu, Gideon Koren, James Weiss, David Paterson, David Christini, Dobromir Dobrev, Jordi Heijman, Thomas O'Hara, Crystal Ripplinger, Zhilin Qu, Jamie Vandenberg, Colleen Clancy, Isabelle Deschenes, Leighton Izu, Tamas Banyasz, Andras Varro, Heike Wulff, Eleonora Grandi, Michael Sanguinetti, Donald Bers, Jeanne Nerbonne and Nipavan Chiamvimonvat as speakers and panel discussants. This article summarizes state-of-the-art knowledge and controversies on the functional roles of cardiac K⁺ channels in normal and diseased heart. We endeavour to integrate current knowledge at multiple scales, from the single cell to the whole organ levels, and from both experimental and computational studies.

Rabbit-specific computational modelling of ventricular cell electrophysiology: Using populations of models to explore variability in the response to ischemia

Computational modelling, combined with experimental investigations, is a powerful method for investigating complex cardiac electrophysiological behaviour. The use of rabbit-specific models, due to the similarities of cardiac electrophysiology in this species with human, is especially prevalent. In this paper, we first briefly review rabbit-specific computational modelling of ventricular cell electrophysiology, multi-cellular simulations including cellular heterogeneity, and acute ischemia. This mini-review is followed by an original computational investigation of variability in the electrophysiological response of two experimentally-calibrated populations of rabbit-specific ventricular myocyte action potential models to acute ischemia. We performed a systematic exploration of the response of the model populations to varying degrees of ischemia and individual ischemic parameters, to investigate their individual and combined effects on action potential duration and refractoriness. This revealed complex interactions between model population variability and ischemic factors, which combined to enhance variability during ischemia. This represents an important step towards an improved understanding of the role that physiological variability may play in electrophysiological alterations during acute ischemia.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Computational models of atrial cellular electrophysiology and calcium handling, and their role in atrial fibrillation

The complexity of the heart makes an intuitive understanding of the relative contribution of ion channels, transporters and signalling pathways to cardiac electrophysiology challenging. Computational modelling of cardiac cellular electrophysiology has proven useful to integrate experimental findings, extrapolate results obtained in expression systems or animal models to other systems, test quantitatively ideas based on experimental data and provide novel hypotheses that are experimentally testable. While the bulk of computational modelling has traditionally been directed towards ventricular bioelectricity, increasing recognition of the clinical importance of atrial arrhythmias, particularly atrial fibrillation, has led to widespread efforts to apply computational approaches to understanding atrial electrical function. The increasing availability of detailed, atrial-specific experimental data has stimulated the development of novel computational models of atrial-cellular electrophysiology and Ca(2+) handling. To date, more than 300 studies have employed mathematical simulations to enhance our understanding of atrial electrophysiology, arrhythmogenesis and therapeutic responses. Future modelling studies are likely to move beyond current whole-cell models by incorporating new data on subcellular architecture, macromolecular protein complexes, and localized ion-channel regulation by signalling pathways. At the same time, more integrative multicellular models that take into account regional electrophysiological and Ca(2+) handling properties, mechano-electrical feedback and/or autonomic regulation will be needed to investigate the mechanisms governing atrial arrhythmias. A combined experimental and computational approach is expected to provide the more comprehensive understanding of atrial arrhythmogenesis that is required to develop improved diagnostic and therapeutic options. Here, we review this rapidly expanding area, with a particular focus on Ca(2+) handling, and provide ideas about potential future directions.

Human atrial cell models to analyse haemodialysis-related effects on cardiac electrophysiology: work in progress

During haemodialysis (HD) sessions, patients undergo alterations in the extracellular environment, mostly concerning plasma electrolyte concentrations, pH, and volume, together with a modification of sympathovagal balance. All these changes affect cardiac electrophysiology, possibly leading to an increased arrhythmic risk. Computational modeling may help to investigate the impact of HD-related changes on atrial electrophysiology. However, many different human atrial action potential (AP) models are currently available, all validated only with the standard electrolyte concentrations used in experiments. Therefore, they may respond in different ways to the same environmental changes. After an overview on how the computational approach has been used in the past to investigate the effect of HD therapy on cardiac electrophysiology, the aim of this work has been to assess the current state of the art in human atrial AP models, with respect to the HD context. All the published human atrial AP models have been considered and tested for electrolytes, volume changes, and different acetylcholine concentrations. Most of them proved to be reliable for single modifications, but all of them showed some drawbacks. Therefore, there is room for a new human atrial AP model, hopefully able to physiologically reproduce all the HD-related effects. At the moment, work is still in progress in this specific field.

Mathematical approaches to understanding and imaging atrial fibrillation: significance for mechanisms and management

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans. The mechanisms that govern AF initiation and persistence are highly complex, of dynamic nature, and involve interactions across multiple temporal and spatial scales in the atria. This article aims to review the mathematical modeling and computer simulation approaches to understanding AF mechanisms and aiding in its management. Various atrial modeling approaches are presented, with descriptions of the methodological basis and advancements in both lower-dimensional and realistic geometry models. A review of the most significant mechanistic insights made by atrial simulations is provided. The article showcases the contributions that atrial modeling and simulation have made not only to our understanding of the pathophysiology of atrial arrhythmias, but also to the development of AF management approaches. A summary of the future developments envisioned for the field of atrial simulation and modeling is also presented. The review contends that computational models of the atria assembled with data from clinical imaging modalities that incorporate electrophysiological and structural remodeling could become a first line of screening for new AF therapies and approaches, new diagnostic developments, and new methods for arrhythmia prevention.

Comparative assessment of nonlinear metrics to quantify organization-related events in surface electrocardiograms of atrial fibrillation

Atrial fibrillation (AF) is today the most common sustained arrhythmia, its treatment being not completely satisfactory. Electrical activity organization analysis within the atria could play a key role in the improvement of current AF therapies. The application of a nonlinear regularity index, such as sample entropy (SampEn), to the atrial activity (AA) fundamental waveform has proven to be a successful noninvasive AF organization estimator. However, the use of alternative nonlinear metrics within this context is a pending issue. The present work analyzes the ability of several nonlinear indices to assess regularity of patterns and, thus, organization, in the AA signal and its fundamental waveform, defined as the main atrial wave (MAW). Precisely, Fuzzy Entropy, Spectral Entropy, Lempel-Ziv Complexity and Hurst Exponents were studied, achieving more robust and accurate AF organization estimates than SampEn. Results also provided better AF organization estimates from the MAW than from the AA signal for all the tested nonlinear metrics, which agrees with previous works only focused on SampEn. Furthermore, some of these indices reported a discriminant ability close to 95% in the classification of AF organization-dependent events, thus outperforming the diagnostic accuracy of SampEn and other widely used noninvasive estimators, such as the dominant atrial frequency (DAF). As a conclusion, these nonlinear metrics could be considered as promising estimators of noninvasive AF organization and could be helpful in making appropriate decisions on the patients' management.

Bifurcations, sustained oscillations and torus bursting involving ionic concentrations dynamics in a canine atrial cell model

Atrial fibrillation is a disorganization of the electrical propagation in the atria often initiated by ectopic beats. This spontaneous activity might be associated with the appearance of sustained oscillations in some portion of the tissue. Adrenergic stress and specific gene polymorphisms known to promote atrial fibrillation are notably related to calcium and potassium channel conductances. We performed codimension-one and two bifurcation analysis along these conductances in an ionic canine atrial myocyte model. Two Hopf bifurcations were found, related to two distinct mechanisms: (1) a fast calcium gating-driven oscillator, and (2) a slow concentration-driven oscillator. These two mechanisms interact through a double Hopf bifurcation (HH) in a neighborhood of which a torus (Neimark-Sacker) bifurcation leads to bursting. A complex codimension-two theoretical scenario was identified around HH, through systematic comparison with the attractors found numerically. The concentration oscillator was further decomposed to reveal the minimal oscillating subnetwork, in which the Na(+)/Ca(2+) exchanger plays a prominent role.

Fibroblast-myocyte electrotonic coupling: does it occur in native cardiac tissue?

Heterocellular electrotonic coupling between cardiac myocytes and non-excitable connective tissue cells has been a long-established and well-researched fact in vitro. Whether or not such coupling exists in vivo has been a matter of considerable debate. This paper reviews the development of experimental insight and conceptual views on this topic, describes evidence in favour of and against the presence of such coupling in native myocardium, and identifies directions for further study needed to resolve the riddle, perhaps less so in terms of principal presence which has been demonstrated, but undoubtedly in terms of extent, regulation, patho-physiological context, and actual relevance of cardiac myocyte–non-myocyte coupling in vivo. This article is part of a Special Issue entitled "Myocyte-Fibroblast Signalling in Myocardium."

•

Electrical coupling of cardiomyocytes and fibroblasts is well-established in vitro

•

Whether such hetero-cellular coupling exists in vivo has been a matter of debate

•

We review the development of experimental and conceptual insight into the topic

•

Conclusion 1: hetero-cellular coupling in heart tissue has been shown in principle

•

Conclusion 2: extent, regulation, context, and relevance remain to be established

Sodium-calcium exchangers (NCX): molecular hallmarks underlying the tissue-specific and systemic functions

NCX proteins explore the electrochemical gradient of Na(+) to mediate Ca(2+)-fluxes in exchange with Na(+) either in the Ca(2+)-efflux (forward) or Ca(2+)-influx (reverse) mode, whereas the directionality depends on ionic concentrations and membrane potential. Mammalian NCX variants (NCX1-3) and their splice variants are expressed in a tissue-specific manner to modulate the heartbeat rate and contractile force, the brain's long-term potentiation and learning, blood pressure, renal Ca(2+) reabsorption, the immune response, neurotransmitter and insulin secretion, apoptosis and proliferation, mitochondrial bioenergetics, etc. Although the forward mode of NCX represents a major physiological module, a transient reversal of NCX may contribute to EC-coupling, vascular constriction, and synaptic transmission. Notably, the reverse mode of NCX becomes predominant in pathological settings. Since the expression levels of NCX variants are disease-related, the selective pharmacological targeting of tissue-specific NCX variants could be beneficial, thereby representing a challenge. Recent structural and biophysical studies revealed a common module for decoding the Ca(2+)-induced allosteric signal in eukaryotic NCX variants, although the phenotype variances in response to regulatory Ca(2+) remain unclear. The breakthrough discovery of the archaebacterial NCX structure may serve as a template for eukaryotic NCX, although the turnover rates of the transport cycle may differ ~10(3)-fold among NCX variants to fulfill the physiological demands for the Ca(2+) flux rates. Further elucidation of ion-transport and regulatory mechanisms may lead to selective pharmacological targeting of NCX variants under disease conditions.